Курсовая работа: Порівняння характеристик аналогового та цифрового фільтрів
Wsz=[Ws1/Wd Ws2/Wd]
[n,wn]=cheb1ord (Wpz, Wsz, Rp, Rs);
В результаті розрахунків маємо:
Мінімальний порядок фільтра n=9.
Частоти зрізу (відносні) Wn = [0.1567 0.2402]
2.2 Визначення коефіцієнтів фільтра та отримання передаточної функції
Для визначення коефіцієнтів фільтра використаємо середовище MatLab, при цьому можливості MatLab дозволяють виключити етап розрахунку фільтра прототипу і одразу визначити коефіцієнти фільтра згідно з ТЗ.
Wp 1 z = Wp 1/ Wd
Wp 2 z = Wp 2/ Wd
[b,a]=cheby1(n,Rp,[Wp1z Wp2z]);
Hrz=tf(b,a,1/(Wd*2*pi))
В результаті розрахунків отримуємо:
Вектор коефіцієнтів чисельника
b = 1.0e-008 * [0.0053 0 -0.0480 0 0.1920 0 -0.4481 0 0.6721 0 -0.6721 0 0.4481 0 -0.1920 0 0.0480 0 -0.0053]
a = 1.0e+004 * [0.0001 -0.0014 0.0102 -0.0462 0.1510 -0.3780 0.7496 -1.2033 1.5843 -1.7234 1.5525 -1.1555 0.7054 -0.3486 0.1365 -0.0409 0.0088 -0.0012 0.0001]
Передаточна функція фільтра:
Побудуємо характеристики отриманого фільтра за допомогою MatLab.
[b,a]=cheby1(n,Rp,[Wp1z Wp2z]);
[H,W]=freqz(b,a,'whole');
figure;
plot (W,abs(H));grid on;
Hrz=tf(b,a,1/(Wd*2*pi))
figure;
pzmap(Hrz)