实验基于MATLAB的FIR数字滤波器设计.docx
- 文档编号:623170
- 上传时间:2022-10-11
- 格式:DOCX
- 页数:17
- 大小:218.73KB
实验基于MATLAB的FIR数字滤波器设计.docx
《实验基于MATLAB的FIR数字滤波器设计.docx》由会员分享,可在线阅读,更多相关《实验基于MATLAB的FIR数字滤波器设计.docx(17页珍藏版)》请在冰豆网上搜索。
实验基于MATLAB的FIR数字滤波器设计
实验4基于MATLAB的FIR数字滤波器设计
实验目的:
加深对数字滤波器的常用指标和设计过程的理解。
实验原理:
低通滤波器的常用指标:
通带边缘频率
,阻带边缘频率
,通带起伏
,
通带峰值起伏
,
阻带起伏
,最小阻带衰减
。
数字滤波器有IIR和FIR两种类型,它们的特点和设计方法不同。
在MATLAB中,可以用b=fir1(N,Wn,’ftype’,taper)等函数辅助设计FIR数字滤波器。
N代表滤波器阶数;Wn代表滤波器的截止频率(归一化频率),当设计带通和带阻滤波器时,Wn为双元素相量;ftype代表滤波器类型,如’high’高通,’stop’带阻等;taper为窗函数类型,默认为海明窗,窗系数需要实现用窗函数blackman,hamming,hanningchebwin,kaiser产生。
例1用凯塞窗设计一FIR低通滤波器,通带边界频率
,阻带边界频率
,阻带衰减
不小于50dB。
解首先由过渡带宽和阻带衰减
来决定凯塞窗的N和
,
,
上图给出了以上设计的频率特性,(a)为N=30直接截取的频率特性(b)为凯塞窗设计的频率特性。
凯塞窗设计对应的MATLAB程序为:
wn=kaiser(30,4.55);
nn=[0:
1:
29];
alfa=(30-1)/2;
hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa));
h=hd.*wn;
[h1,w1]=freqz(h,1);
或者:
b=fir1(29,0.4,kaiser(30,4.55));
[h1,w1]=freqz(b,1);
plot(w1/pi,20*log10(abs(h1)));
axis([0,1,-80,10]);
grid;
xlabel('归一化频率/π');
ylabel('幅度/dB');
还可以使用[n,Wn,beta,ftype]=kaiserord(f,a,dev)函数来估计滤波器阶数等,得到凯塞窗滤波器:
fcuts=[0.30.5];%归一化频率omega/pi
mags=[10];
devs=[0.0510^(-2.5)];
[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs);%计算出凯塞窗N,beta的值
hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');
freqz(hh);
实际中,一般调用MATLAB信号处理工具箱函数remezord来计算等波纹滤波器阶数N和加权函数W(ω),调用函数remez可进行等波纹滤波器的设计,直接求出滤波器系数。
函数remezord中的数组fedge为通带和阻带边界频率,数组mval是两个边界处的幅值,而数组dev是通带和阻带的波动,fs是采样频率单位为Hz。
例2利用雷米兹交替算法设计等波纹滤波器,设计一个线性相位低通FIR数字滤波器,其指标为:
通带边界频率fc=800Hz,阻带边界fr=1000Hz,通带波动
阻带最小衰减At=40dB,采样频率fs=4000Hz。
解
在MATLAB中可以用remezord和remez两个函数设计,其结果如图2,MATLAB程序如下:
fedge=[8001000];
mval=[10];
dev=[0.05590.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)));
grid;
xlabel('频率/Hz');
ylabel('幅度/dB');
一、实验内容:
利用MATLAB编程设计一个数字带通滤波器,指标要求如下:
通带边缘频率:
,
,通带峰值起伏:
。
阻带边缘频率:
,
,最小阻带衰减:
。
分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器。
实验要求:
给出FIR数字滤波器的冲激响应,绘出它们的幅度和相位频响曲线,讨论它们各自的实现形式和特点。
1-1)用窗函数法实现:
调用函数[n,wn,bta,ftype]=kaiserord(f,a,dev,fs)
参数:
f=[0.30.450.650.8]为对应数字频率
,
,
,,
a=[010]为由f指定的各个频带上的幅值向量,一般只有0和1表示;和f长度关系为(2*a的长度)—2=(f的长度)
devs=[0.010.10870.01]用于指定各个频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或偏差,长度与a相等,计算公式:
阻带衰减误差=
,通带波动衰减误差=
fs缺省值为2HZ
>>[n,wn,bta,ftype]=kaiserord([0.30.450.650.8],[010],[0.010.10870.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=
Columns1through8
0.00410.0055-0.0091-0.0018-0.0056-0.00000.0391-0.0152
Columns9through16
-0.03810.0077-0.02930.09400.0907-0.2630-0.05170.3500
Columns17through24
-0.0517-0.26300.09070.0940-0.02930.0077-0.0381-0.0152
Columns25through31
0.0391-0.0000-0.0056-0.0018-0.00910.00550.0041
图4-1
如果直接用freqz(h1,1,256);
画图得:
1-2)用等波纹法设计:
调用函数[n,fpts,mag,wt]=remezord(f,a,dev)
f=[0.30.450.650.8]
a=[010]
dev=[0.010.10870.01]
其含义同函数[n,wn,bta,ftype]=kaiserord(f,a,dev,fs)中的参数相同。
>>[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)
>>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
h2=
Columns1through9
-0.00130.0092-0.0255-0.06420.11770.0922-0.2466-0.04660.3116
Columns10through17
-0.0466-0.24660.09220.1177-0.0642-0.02550.0092-0.0013
图4-2
用freqz(h2,1,256);直接得图:
二、对课本作业9.23画图。
2-1)用汉宁窗实现:
1)手动计算
>>n=0.001:
58.001;
>>hd=sin(0.18125*pi*(n-29))./(pi*(n-29));
>>win=0.5+0.5*cos(2*pi*(n-29)/58);
>>h1=2*cos(pi*(n-29)/2).*hd.*win;
>>[hh1,w1]=freqz(h1,1,256);
>>figure
(1)
>>subplot(2,1,1)
>>plot(w1,20*log10(abs(hh1)))
>>grid
>>xlabel('数字频率w/rad');ylabel('幅度/db');
>>subplot(2,1,2)
>>plot(w1,angle(hh1))
>>grid
>>xlabel('数字频率w/rad');ylabel('相位/rad');
图4-3
2-2)用自带的fir1函数:
>>n=59;
>>wn=[3/85/8];
>>h2=fir1(n,wn,'bandpass',hann(n+1));
>>[hh2,w2]=freqz(h2,1,256);
>>figure
(2)
>>subplot(2,1,1)
>>plot(w2,20*log10(abs(hh2)))
>>xlabel('数字频率w/rad');ylabel('幅度/db');
>>grid
>>subplot(2,1,2)
>>plot(w2,angle(hh2))
>>xlabel('数字频率w/rad');ylabel('相角/rad');
>>grid
图4-4
2-3)用等波纹法设计:
调用函数[n,fpts,mag,wt]=remezord(f,a,dev)
f=[0.26250.3750.6250.7375]
a=[010]
dev=[0.010.10870.01]
>>[n,fpts,mag,wt]=remezord([0.26250.3750.6250.7375],[010],[0.010.10870.01]);>>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('归一化频率
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 基于 MATLAB FIR 数字滤波器 设计