Спроектируйте минимальный порядок фильтр lowpass с частотой среза полосы пропускания на 500 Гц и частотой среза полосы задерживания на 600 Гц. Задайте частоту дискретизации 2 000 Гц. Потребуйте по крайней мере 40 дБ затухания в полосе задерживания и меньше чем 3 дБ пульсации в полосе пропускания.

rp = 3;           % Passband ripple in dB 
rs = 40;          % Stopband ripple in dB
fs = 2000;        % Sampling frequency
f = [500 600];    % Cutoff frequencies
a = [1 0];        % Desired amplitudes

Преобразуйте отклонения в линейные модули. Спроектируйте фильтр и визуализируйте его величину и фазовые отклики.

dev = [(10^(rp/20)-1)/(10^(rp/20)+1) 10^(-rs/20)]; 
[n,fo,ao,w] = firpmord(f,a,dev,fs);
b = firpm(n,fo,ao,w);
freqz(b,1,1024,fs)
title('Lowpass Filter Designed to Specifications')

Фильтр немного далек от встречи технические требования неравномерности в полосе пропускания и затухание в полосе задерживания. Используя n+1 вместо n в вызове firpm достигает желаемых амплитудных характеристик.

Спроектируйте фильтр lowpass с частотой среза полосы пропускания на 1 500 Гц и частотой среза полосы задерживания на 2 000 Гц. Задайте частоту дискретизации 8 000 Гц. Потребуйте максимальной амплитуды полосы задерживания 0,1 и максимальной ошибки полосы пропускания (пульсация) 0,01.

[n,fo,ao,w] = firpmord([1500 2000],[1 0],[0.01 0.1],8000);
b = firpm(n,fo,ao,w);

Получите эквивалентный результат при наличии firpmord сгенерируйте массив ячеек. Визуализируйте частотную характеристику фильтра.

c = firpmord([1500 2000],[1 0],[0.01 0.1],8000,'cell');
B = firpm(c{:});
freqz(B,1,1024,8000)