Создайте lowpass allpole БИХ-фильтр с полюсом в 0,9. Вычислите эффективную импульсную продолжительность ответа. Получите импульсный ответ. Постройте результат.

b = 1;
a = [1 -0.9];
len = impzlength(b,a)

len = 93

[h,t] = impz(b,a);
stem(t,h)

h(len)

ans = 6.1704e-05

Разработайте 4-й порядок lowpass эллиптический фильтр с частотой среза 0.4π рад/выборка. Задайте 1 дБ пульсации полосы пропускания и 60 дБ затухания полосы задерживания. Разработайте фильтр в нулевом усилении полюса, формируют и получают матрицу раздела второго порядка использование zp2sos. Определите эффективную импульсную длину последовательности ответа из матрицы раздела второго порядка.

[z,p,k] = ellip(4,1,60,.4);
[sos,g] = zp2sos(z,p,k);
len = impzlength(sos)

len = 80

Используйте designfilt, чтобы разработать 4-й порядок lowpass эллиптический фильтр с нормированной частотой полосы пропускания 0.4π рад/выборка. Задайте 1 дБ пульсации полосы пропускания и 60 дБ затухания полосы задерживания. Определите эффективную импульсную длину последовательности ответа и визуализируйте ее.

d = designfilt('lowpassiir','FilterOrder',4,'PassbandFrequency',0.4, ...
               'PassbandRipple',1,'StopbandAttenuation',60, ...
               'DesignMethod','ellip');
len = impzlength(d)

len = 80

impz(d)

Этот пример требует программного обеспечения DSP System Toolbox™.

Разработайте 4-й порядок lowpass эллиптический фильтр с частотой среза 0.4*pi рад/выборка. Задайте 1 дБ пульсации полосы пропускания и 60 дБ затухания полосы задерживания. Разработайте фильтр в нулевом усилении полюса, формируют и получают матрицу раздела второго порядка с помощью zp2sos. Создайте biquad, фильтруют Систему object™ и вводят Системный объект к impzlength.

[z,p,k] = ellip(4,1,60,.4);
[sos,g] = zp2sos(z,p,k);
hBqdFilt = dsp.BiquadFilter('Structure','Direct form I',...
                                            'SOSMatrix', sos,...
                                            'ScaleValues',g);
len = impzlength(hBqdFilt)

len = 80

Разработайте БИХ Баттерворт и КИХ equiripple фильтры для данных, выбранных на уровне 1 кГц. Частота полосы пропускания составляет 100 Гц, и частота полосы задерживания составляет 150 Гц. Пульсация полосы пропускания составляет 0,5 дБ и существует 60 дБ затухания полосы задерживания. Получите Системные объекты для фильтров и сравните импульсные длины последовательности ответа фильтра.

d = fdesign.lowpass('Fp,Fst,Ap,Ast',100,150,0.5,60,1000);
Hd1 = design(d,'butter','SystemObject',true);
Hd2 = design(d,'equiripple','SystemObject',true);
len = [impzlength(Hd1) impzlength(Hd2)]

len = 1×2

   183    49