Курсовая работа: Реализация и анализ ЦФ с КИХ

i=1; % инициализация программного счетчика

Ni=120; % задание количества отсчётов импульса

h=[1:Ni]; % подготовка массива отсчетов ИХ к последующему заполнению

ti=120; % задание длины импульса в точках

Fd=120; % задание частоты дискретизации

F0=Fd/4; %задание несущей частоты

Wd=2*pi*Fd; % задание угловой частоты дискретизации

W0=2*pi*F0; %задание угловой несущей частоты

dF=Fd/4; % задание девиации частоты

dW=dF*2*pi; % задание угловой девиации частоты

mju=dF/ti; % задание скорости роста частоты (скорость частотной модуляции)

n=0:1/Fd:1-1/Fd; % подготовка дискретов времени с шагом, равным периоду дискретизации

Sn=chirp(n,0,1,dF); % заполнение массива отсчетов дискретизированного сигнала

figure(1)

bar(Sn,0.1); % вывод на экран графика входного сигнала

f=abs(fft(Sn,Ni)); % вычисление спектра входного сигнала с использованием БПФ

figure(2);

stem(f); % постоение графика спектра входного сигнала

w=hann(Ni); % вычисление Ni-точечного окна Ханна

figure(3);

stem(w); % построение графика вычисленного окна

Ww=abs(fft(w,Ni)); % вычисление амплитудного спектра вычисленного окна с помощью БПФ

figure(4);

plot(Ww); % построени найденного спектра

while i<=Ni, % цикл нахождения ИХ СФ (путем зеркального отражения входного сигнала)

h(i)=Sn(Ni-i+1);

i=i+1;

end

hw=h.*w'; % взвешивание фильтра окном полученным окном Ханна