Спектральные преобразования сигналов


Спектральные преобразования сигналов

Выполните спектральный и корреляционный анализ импульсного сигнала в соответствии с вашим вариантом (см. таблицу 1 и рисунок 1).

Определите спектральную плотность одиночного импульса, его амплитудный, фазовый и энергетический спектр. Постройте графики спектров.

Определите энергию импульса.

Определите АКФ импульса, постройте ее график. Покажите связь АКФ с энергетическим спектром импульса.

Найдите амплитудный и фазовый спектры, разложив в ряд Фурье последовательность импульсов с периодом. Постройте графики зависимостей амплитуд и фаз гармоник ряда от частоты [Гц].

Выполните проверку расчетов пунктов 1-4 в Matlab.

Отчет должен содержать: ход и результаты математических расчетов (пункты 1 - 4), код программы в Matlab и результаты ее работы.

Сигнал:

10;

Спектральный сигнал импульс математический

1) Спектральная плотность сигнала:

Амплитудный спектр:

Фазовый спектр:

Энергетический спектр:

    2) Энергия сигнала: 3) 4) АКФ импульса:

Период сигнала T=5*tau=12.5*10-3

BN=0 , так как функция четная.

    5) Код программы Matlab: 1) clear; clc; 2) t=-2.5*10^(-3):0.001:2.5*10^(-3); 3) Сигнал: 4) Y=zeros(size(t)); 5) s=Y+10; 6) plot(t, s), grid on; 7) Спектральная плотность сигнала: 8) w=fft(s); 9) A=abs(w); 10) 11) Фазовый спектр: 12) FI=angle(w); 13) Энергитический спектр: 14) W=(abs(w)).^2; 15) 16) subplot(1,1,1), plot(t, s), grid on; 17) title('Сигнал'); 18) xlabel('t, сек'); 19) ylabel('s(t)'); 20) freq=-100:100; 21) subplot(1,1,1), plot(freq(100:102),fftshift(A(1:3))), grid on; 22) title('Cпектральная плотность сигнала(АЧС)'); 23) xlabel('f, Гц'); 24) ylabel('|S(w)|'); 25) 26) subplot(1,1,1), plot(freq(100:102),FI(1:3)), grid on; 27) title('Фазовый спектр'); 28) xlabel('f, Гц'); 29) ylabel('FI(w)'); 30) 31) subplot(1,1,1), plot(freq(100:102),fftshift(W(1:3))), grid on; 32) title('Энерегетический спектр'); 33) xlabel('f, Гц'); 34) ylabel('W(f)'); 35) c=xcorr(s); 36) Tau=0:10; 37) subplot(1,1,1), plot(Tau, c), grid on; 38) title('Корреляция'); 39) xlabel('Tau'); 40) ylabel('B(Tau)'); 41) 42) w=0.16*pi; 43) for n=1:5 44) w1(n)=w*n; 45) A(n)=3200.*(sin(0.001.*n.*pi))./(0.001.*pi.*n); 46) plot([A(n) A(n)],[0 w1(n)]), grid on; 47) hold on; 48) end 49) title('Амплитудный спектр периодического сигнала'); 50) xlabel('f, Гц'); 51) ylabel('|A(f)|');

Похожие статьи




Спектральные преобразования сигналов

Предыдущая | Следующая