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

bk_out(i)=(2*dt/T1)*sum(out_sign.*sin(t*i*W1));

end;

a0_out=dt*sum(out_sign)/T1;

Ak_out=sqrt(ak_out.^2+bk_out.^2);

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

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

figure;

hold on

stem(Wn,[a0 Akn],'ro');

stem(Wn,[a0_out Akn_out],'b+');

xlabel('w'); ylabel('Magnitude');

axis([0 150000 0 70]);

hold off

grid on

Рис. 9. Спектри вхідного та вихідного сигналу.

2 . Розрахунок та дослідження цифрового фільтра

2 .1 Визначення порядку та частоти зрізу цифрового фільтра

Як і у випадку проектування аналогових фільтрів для визначення порядку цифрових фільтрів та їх частоти зрізу можна скористатися функціями пакета Communication середовища MatLab.

Згідно з теоремою Котельникова частота дискретизації має бути принаймні вдвічі більшою ніж максимальна частота сигналу. Для моделі фільтра використаємо частоту , де - частота дискретизації, - верхня межи смуги затримки.

% Початкові дані

Wp1=8442; % Нижня межа смуги пропускання рад / сек

Wp2=12940; % Верхня межа смуги пропускання рад / сек

Ws1=7914; % Нижня межа смуги затримки рад / сек

Ws2=13468; % Верхня межа смуги затримки рад / сек

Rp=2; % Коефіцієнт пульсацій дБ

Rs=37; %Затухан ня дБ

Wd=Ws2*4 % Частота дискретизації рад / сек