Лабораторная работа: Расчет цифровых фильтров с бесконечными импульсными характеристиками

plot(w,abs(h),'g');

subplot(312);

plot(w,unwrap(angle(h)),'g');

subplot(313);

grpdelay(bt,at,64,Fs);%%%групповая задержка

Рисунок 5 – Диаграмма нулей и полюсов

Рисунок 6 – АЧХ , ФЧХ и групповая задержка фильтра

Задание №4: Рассчитайте цифровой ФВЧ с параметрами заданными для каждого варианта. Определить какой порядок будут иметь разные типы фильтров. Сравнить эффективность аппроксимаций при более жестких и более мягких требованиях к АЧХ.

close all;

clear all;

clear all;

Wp = 0.4; % нормированная граничная частота подавления

Ws = 0.7; % нормированная граничная частота пропускания

Rp = 0.5; % допустимая неравномерность в полосе пропускания

Rs = 70; % минимальное затухание в полосе подавления

%%%%%%%%%%%%%%Фильтр Баттерворта%%%%%%%%%%%%

[n1, Wn1] = buttord(Wp, Ws, Rp, Rs);

sprintf( 'battervorta - %d',n1)

[z,p,k] = butter(n1,Wn1,'high');

[b,a]=zp2tf(z,p,k);

figure;

zplane(z,p);

figure;

title('battervorta');

freqs(z,p);

%%%%%%%%%%%%%%Фильтр Чебышева 1 рода%%%%%%%%

[n2, Wn2] = cheb1ord(Wp, Ws, Rp, Rs);