Курсовая работа: Порівняння характеристик аналогового та цифрового фільтрів

figure(7);

lsim(Hrez,Analog_sign,t);

grid on

axis([0 0.02 -40 40]);

Input signal

Рис. 8. Вхідний та вихідний сигнали.

1 .5 Побудова спектрів вхідного та вихідного сигналів

Побудуємо спектри вхідного та вихідного сигналів. Для побудови спектрів використаємо алгоритм неперервного перетворення Фур’є.

Для дійсних функцій ряд Фур’є часто означається, як ряд по синусах та косинусах, аргументи яких обираються в такий спосіб, аби на відрізку означення функції вкладалась ціла кількість періодів. Легко зрозуміти, що перелік таких синусів та косинусів вичерпується множинами та , де n є довільним натуральним числом, а для може приймати також значення 0.

%------------------------- Побудова спектрів ---------------------------%

[ y , x ]= lsim ( Hrez , Analog _ sign , t );

a0=dt*sum(Analog_sign)/T1;

for i=1:k

ak(i)=(2*dt/T1)*sum(Analog_sign.*cos(t*i*W1));

bk(i)=(2*dt/T1)*sum(Analog_sign.*sin(t*i*W1));

end;

Ak=sqrt(ak.^2+bk.^2);

k1=1:i; % номера відліків

W0=k1*W1; % масштабовані відліки

Wn=zeros(1,50); % завдання розміру вектора відліків

Wn =0: W 0:50* W 0; % заповнення вектора відліків (перші 50 елементів вектора W0)

Akn=zeros(1,50); % завдання розміру вектора амплітуд з урахуванням а0

Akn(1,1:50)=Ak(1,1:50); %заполнение вектора амплитуд

out_sign=y';

for i=1:k