Проектирование фильтра нижних частот минимального порядка с частотой отсечки полосы пропускания 500 Гц и частотой отсечки полосы отсечки 600 Гц. Укажите частоту дискретизации 2000 Гц. Требуется не менее 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 достигается желаемая амплитудная характеристика.

Проектирование фильтра нижних частот с частотой отсечки полосы пропускания 1500 Гц и частотой отсечки полосы отсечки 2000 Гц. Укажите частоту дискретизации 8000 Гц. Требуется максимальная амплитуда полосы останова 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)