Проектируйте lowpass минимального порядка с частотой отключения полосы пропускания 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 достигает желаемых амплитудных характеристик.

Проектируйте lowpass с частотой отключения полосы пропускания 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)