designfir.m

来自「这是我编的FIR滤波器的dsp程序,可调成滤波器阶数,汇编语言编程,执行效率高,」· M 代码 · 共 12 行

rp = 3;          % Passband ripple
rs = 50;         % Stopband ripple
fs = 2000;       % Sampling frequency
f = [500 600];   % Cutoff frequencies
a = [1 0];       % Desired amplitudes
% Compute deviations,[0.01 0.1]
dev = [(10^(rp/20)-1)/(10^(rp/20)+1)  10^(-rs/20)]; 
[n,fo,ao,w] = firpmord(f,a,[0.01 0.1],fs);
b = firpm(n,fo,ao,w);
display(b);
freqz(b,1,1024,fs);
title('Lowpass Filter Designed to Specifications');