电子技术论坛_中国专业的电子工程师学习交流社区-中电网技术论坛

标题: 基于Matlab的FIR滤波器设计与实现（2） [打印本页]

作者: yuyang911220 时间: 2015-9-21 20:24 标题: 基于Matlab的FIR滤波器设计与实现（2）

方法一：窗函数法

程序如下：
[n,wn,bta,ftype]=kaiserord([0.3 0.45 0.65 0.8],[0 1 0],[0.01 0.1087 0.01]);%用kaiserord函数估计出滤波器阶数n和beta参数
h1=fir1(n,wn,ftype,kaiser(n+1,bta),'noscale');
[hh1,w1]=freqz(h1,1,256);
figure(1)
subplot(2,1,1)
plot(w1/pi,20*log10(abs(hh1)))
grid
xlabel('归一化频率w');ylabel('幅度/db');
subplot(2,1,2)
plot(w1/pi,angle(hh1))
grid
xlabel('归一化频率w');ylabel('相位/rad');

波形如下：

滤波器系数为：
h1 =

  Columns 1 through 8

0.0041 0.0055 -0.0091 -0.0018 -0.0056 -0.0000 0.0391 -0.0152

  Columns 9 through 16

-0.0381 0.0077 -0.0293 0.0940 0.0907 -0.2630 -0.0517 0.3500

  Columns 17 through 24

-0.0517 -0.2630 0.0907 0.0940 -0.0293 0.0077 -0.0381 -0.0152

  Columns 25 through 31

0.0391 -0.0000 -0.0056 -0.0018 -0.0091 0.0055 0.0041

如果直接用freqz(h1,1,256)，得幅频特性和相频特性曲线：

方法二：等波纹法设计
程序如下：
[n,fpts,mag,wt]=remezord([0.3 0.45 0.65 0.8],[0 1 0],[0.01 0.1087 0.01]);%用remezord函数估算出remez函数要用到的阶n、归一化频带边缘矢量fpts、频带内幅值响应矢量mag及加权矢量w，使remez函数设计出的滤波器满足f、a及dev指定的性能要求。
h2=remez(n,fpts,mag,wt);%设计出等波纹滤波器
[hh2,w2]=freqz(h2,1,256);
figure(2)
subplot(2,1,1)
plot(w2/pi,20*log10(abs(hh2)))
grid
xlabel('归一化频率w');ylabel('幅度/db');
subplot(2,1,2)
plot(w2/pi,angle(hh2))
grid
xlabel('归一化频率w');ylabel('相位/rad');
h2

波形如下：

欢迎光临电子技术论坛_中国专业的电子工程师学习交流社区-中电网技术论坛 (http://bbs.eccn.com/)