基于Matlab的FIR滤波器设计与实现Word文档格式.docx
《基于Matlab的FIR滤波器设计与实现Word文档格式.docx》由会员分享，可在线阅读，更多相关《基于Matlab的FIR滤波器设计与实现Word文档格式.docx（13页珍藏版）》请在冰豆网上搜索。
手动计算滤波器阶数N和β值，之后在通过程序设计出滤波器。
第一步：
通过过渡带宽度和阻带衰减，计算滤波器的阶数B和β值。
第二步：
通过程序设计滤波器。
程序如下：
b
=
fir1（29,0.4,kaiser（30,4.55））;
[h1,w1]=freqz（b,1）;
figure
（1）
plot（w1/pi,abs（h1））;
grid;
xlabel（'
归一化频率/p'
）
;
ylabel（'
幅度/dB'
figure
（2）
plot（w1/pi,angle（h1））;
相位'
波形如下：
例2
利用雷米兹交替算法设计等波纹滤波器，设计一个线性相位低通FIR数字滤波器，其指标为：
通带边界频率fc=800Hz，阻带边界fr=1000Hz，通带波动
阻带最小衰减At=40dB，采样频率fs=4000Hz。
一般调用MATLAB信号处理工具箱函数remezord来计算等波纹滤波器阶数N和加权函数W（ω），调用函数remez可进行等波纹滤波器的设计，直接求出滤波器系数。
函数remezord中的数组fedge为通带和阻带边界频率，数组mval是两个边界处的幅值，而数组dev是通带和阻带的波动，fs是采样频率单位为Hz。
解：
在MATLAB中可以用remezord
和remez两个函数设计
fedge=[800
1000];
mval=[1
0];
dev=[0.0559
0.01];
fs=4000;
[N,fpts,mag,wt]=remezord（fedge,mval,dev,fs）;
b=remez（N,fpts,mag,wt）;
[h,w]=freqz（b,1,256）;
plot（w*2000/pi,20*log10（abs（h）））;
频率/Hz'
）;
例3
利用MATLAB编程设计一个数字带通滤波器，指标要求如下：
通带边缘频率：
Ωp1=0.45pi，Ωp2=0.65pi，通带峰值起伏：
δ1<
=1[dB]。
阻带边缘频率：
Ωs1=0.3pi，Ωs2=0.8pi，最小阻带衰减：
δ2>
=40[dB]
。
窗函数法
[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
归一化频率w'
幅度/db'
subplot（2,1,2）
plot（w1/pi,angle（hh1））
相位/rad'
滤波器系数为：
h1
=
Columns
through
8
0.0041
0.0055
-0.0091
-0.0018
-0.0056
-0.0000
0.0391
-0.0152
9
16
-0.0381
0.0077
-0.0293
0.0940
0.0907
-0.2630
-0.0517
0.3500
17
24
25
31
0.0041
如果直接用freqz（h1,1,256），得幅频特性和相频特性曲线：
方法二：
等波纹法设计
[n,fpts,mag,wt]=remezord（[0.30.450.650.8],[010],[0.010.10870.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）
plot（w2/pi,20*log10（abs（hh2）））
plot（w2/pi,angle（hh2））
h2
滤波器系数如下：
h2
-0.0013
0.0092
-0.0255
-0.0642
0.1177
0.0922
-0.2466
-0.0466
0.3116
10
-0.0013
如果直接用freqz（h2,1,256）;
得幅频特性和相频特性曲线：
方法三：
采用FDATool工具
这种方法需要事先计算出滤波器的阶数，bate值，然后设置相应参数，最后生成滤波器。
设置界面如下图所示：
将上述圈圈的区域设置好之后，生成滤波器，最后通过analysis菜单可以观察生成的滤波器的各种特性曲线和滤波器系数。
这里的滤波器系数跟方法一的一样。
五、结果分析
5.1
滤波器设计总结
FIR滤波器实现一般采用窗函数法和等纹波设计法。
窗函数法还包含两个分支，一种是用公式先手动算出N值和其他对应得窗函数参数值，再代入窗函数和fir1实现，一种是用函数*rord估算出N和相应参数再用fir1实现。
不过要注意*rord会低估或高估阶次n，可能会使滤波器达不到指定的性能，这时应稍微增加或降低阶次。
如果截止频率在0或Nyquist频率附近，或者设定的dev值较大，则得不到正确结果。
滤波器实现形式及特点：
由于一般的滤波器在利用窗函数是其通带波纹和阻带波纹不同（一般为第一个阻带波纹最大）因此，在满足第一个阻带衰减旁瓣时，比其频率高的旁瓣，它们的衰减都大大超出要求。
而根据阻带衰减与项数的近似关系N=P（δ2）*fs/TW，可得当阻带衰减越大，所需项数越多。
5.2
窗函数法和等波纹设计的不同之处
窗函数设计是通过最小平方积分办法来设计的，即该滤波器的误差为：
即要求最小方法来设计滤波器，这样的滤波器更忠实于理想滤波器（即滤波系数更接近于理想滤波器）。
证明如下：
因此，幅度频谱差值越小，实际滤波器就越接近理想滤波器。
而等波纹滤波器是通过最大加权误差最小化来实现，其误差为：
要求该误差最小来实现滤波器，得出来的滤波系数较窗函数设计相差较远。
以下通过对例3中的h1及h2作比较。
%sigsum是用来对数组各元素进行求和
function
y=sigsum（n1,n2,n,x）;
y=0;
for
i=n1+1-min（n）:
n2+1-min（n）
y=y+x（i）;
end
n=0.001:
30.001;
h=2*cos（0.55*pi*（n-15））.*sin（0.175*pi*（n-15））./（pi*（n-15））;
delta1=h-h1;
16.001;
delta2=h-h2;
y1=sigsum（0,30,[0:
30],（abs（delta1）.^2））/31;
y2=sigsum（0,16,[0:
16],（abs（delta2）.^2））/17;
结果如下：
y1
1.9099e-004
y2
0.0278
由此得到用窗函数实现的滤波系数比用等波纹滤波器系数的每一项更接近于理想滤波器（y1为用窗函数实现的与理想滤波器的差值，y2为用等波纹滤波器实现的与理想滤波器的差值）；
对比二者的幅度频谱可知，等波纹滤波器阻带边缘比用窗函数实现的更平滑（理想滤波器为垂直下降的）。
从设计的角度考虑，由于窗函数设计法都是通过已有的窗函数对理想滤波器的改造，因此，可以用手算的办法方便的设计滤波器。
而等波纹滤波器，其实现是通过大量的迭代运算来实现，这样的方法一般只能通过软件来设计。
项数的问题由于等波纹滤波器能较平均的分布误差，因此对于相同的阻带衰减，其所需的滤波系数比窗函数的要少。
5.3
几点说明
1.相频特性曲线形状不同说明
上面第一个图是用角度为单位画出来的，下面的图是用rad单位画出来的。
从图形可以观察到在0.3到0.8数字频率间两个图都是严格的线性相位，至于下面的图为什么在这个区间会有跳变是因为rad的区间只有-pi——pi，当相位由-pi继续增加时只能跳到pi而不能大于pi，而角度表示则可以连续增大。
2.调用firl或者reme函数时，用scale（缺省方式）对滤波器进行归一化，即滤波器通带中心频率处的响应幅值为0db。
用noscale不对滤波器归一化。
升级会员 