firfilterdesign.m

这是matlab中的一些有用并且常见的原代码

wp=0.2;ws=0.4;Rp=1;As=20;                    % 定义滤波器参数(可修改)
[n,wc]=buttord(wp,ws,Rp,As);                 % 计算滤波器阶次和3dB截止频率
[b,a]=butter(n,wc);                          % 设计巴特沃斯滤波器
w=[0:pi/200:pi];
H=freqz(b,a,w);                              % 求出[0，pi]范围上的频率响应
magH=abs(H);phaH=angle(H);unwrap(phaH);      % 计算幅度，相位并对相位进行抖动消除
magHdB=10*log10(magH);                       % 计算相对幅度
delta=[1,zeros(1,100)];
h=filter(b,a,delta);                         % 计算脉冲响应
% 以下为绘图命令，将其补充的更完善
subplot(2,2,1);plot(h);grid on;
subplot(2,2,2);plot(w/pi,magHdB);grid on
subplot(2,2,3);plot(w/pi,magH);grid on
subplot(2,2,4);plot(w/pi,phaH/pi);grid on