实验五IIR数字滤波器设计与滤波附思考题程序教学文案Word文档下载推荐.docx
- 文档编号:16302305
- 上传时间:2022-11-22
- 格式:DOCX
- 页数:8
- 大小:184.12KB
实验五IIR数字滤波器设计与滤波附思考题程序教学文案Word文档下载推荐.docx
《实验五IIR数字滤波器设计与滤波附思考题程序教学文案Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《实验五IIR数字滤波器设计与滤波附思考题程序教学文案Word文档下载推荐.docx(8页珍藏版)》请在冰豆网上搜索。
第二种常用指标是用参数和A表示通带和阻带要求,如图所示:
出另外两个参数d,k,分别称为判别因子和选择性因子。
字频率和等效的模拟频率之间的映射关系:
tan©
),由于二者不是线性关
BUTTERWOR低通滤波器:
幅度平方函数定义为
Ha(J)1/、、2,为3DB带宽。
BUTTERWOR低通滤波器系统函数有以下形式:
由模拟滤波器设计IIR数字滤波器,必须建立好s平面和z平面的映射关
系。
使模拟系统函数Ha(s)变换成数字滤波器的系统函数H(z),通常采用冲激
相应不变法和双线性变换法。
冲激相应不变法存在频谱混叠现象,双线性变换
法消除了这一线象,在IIR数字滤波器的设计中得到了更广泛的应用。
系,所以称为预畸变。
3.实验内容及其步骤
实验的步骤:
(1)给定数字滤波器的幅度相应参数。
(2)用预畸变公式将数字滤波器参数变换为相应的等效模拟滤波器参数。
(3)采用模拟滤波器设计方法设计等效模拟滤波器Ha(s)
(4)采用双线性变换公式把等效模拟滤波器映射为所期望的数字滤波器。
其中第三步中模拟滤波器设计步骤为:
首先,根据滤波器指标求选择因子k和判别因子d
再次,设3db截止频率c最后由表查出归一化巴特沃斯滤波器系数。
设计举例:
%技术指标
动,在50rad/s处具有至少30dB的衰减。
求出级联形式的系统函数,画出滤波器的幅度响应、对数幅度响应、相位响应和脉冲响应图
MATLA参考程序:
Wp=30;
Ws=50;
Rp=1;
As=30;
Ripple=10A(-Rp/20);
Attn=10*As/20);
[b,a]=afd_butt(Wp,Ws,Rp,As)
[db,mag,pha,w]=freqs_m(b,a,50);
figure1);
clf;
subplot(2,2,1);
plot(w,mag);
title('
MagnitudeResponse'
);
xlabel('
Analogfrequencyinrad/s'
);
ylabel('
H'
axis([0,50,0,1.1]);
grid;
subplot(2,2,2);
plot(w,db);
MagnitudeindB'
ylabel('
decibels'
axis([0,50,-40,5])
grid
subplot(2,2,3);
plot(w,pha/pi);
PhaseResponse);
radians'
axis([0,50,-1.1,1.1])
gridsubplot(2,2,4);
plot(t,ha);
ImpulseResponse'
timeinseconds'
ylabel('
ha(t)'
axis([0,max(t)+0.05,min(ha),max(ha)+0.025]);
grid
function[b,a]=afd_butt(Wp,Ws,Rp,As);
ifWp<
=0
error('
PassbandedgemustbelargerthanO'
)
end
ifWs<
=Wp
StopbandedgemustbelargerthanPassededge'
if(Rp<
=0)|(As<
0)
PBrippleand/0rSBattenuationmustbelargerthan0'
N=ceil((log10((10A(Rp/10)-1)/(10A(As/10)-1)))/(2*log10(Wp/Ws)));
OmegaC=Wp/((10A(Rp/10)-1F(1/(2*N)));
[b,a]=u_buttap(N,OmegaC);
%设计非归一化巴特沃兹模拟低通滤波器原型子程序
function[b,a]=u_buttap(N,OmegaC);
[z,p,k]=buttap(N);
z、p和k分别是设计出的G(p)的极点、零点及增益
p=p*OmegaC;
%非归一化
k=k*OmegaCAN;
B=real(poly(z));
bO=k;
b=k*B;
分子向量
a=real(poly(p));
分母向量
%+算系统函数的幅度响应和相位响应子程序
function[db,mag,pha,w]=freqs_m(b,a,wmax);
w=[0:
1:
500]*wmax/500;
H=freqs(b,a,w);
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
例2设计如下指标的数字低通滤波器。
相应的MATLAB程序如下:
fp=100;
fst=300;
Fs=1000;
rp=3;
rs=20;
wp=2*pi*fp/Fs;
ws=2*pi*fst/Fs;
Fs=Fs/Fs;
%letFs=1
wap=tan(wp/2);
was=tan(ws/2);
预畸变公式
[n,wn]=buttord(wap,was,rp,rs,'
s'
设计等效模拟滤波器
[z,p,k]=buttap(n);
[bp,ap]=zp2tf(z,p,k)设计模拟低通原型滤波器
[bs,as]=lp2lp(bp,ap,wap);
将模拟低通原型滤波器转换为低通滤波器
[bz,az]=bilinear(bs,as,Fs/2)实现双线性变换,即由模拟滤波器Hs得到数
字滤波器
[h,w]=freqz(bz,az,256,Fs*1000);
plot(w,abs(h));
gridon;
调试运行该程序,运行结果如下:
bp=[1,0,0],ap=[1,1.4142,1],bs=[0.1056,0,0],as=[1,0.4595,0.1056]
bz=[0.0675,0.1349,0.06745],az=[1,-1.143,0.4128]
4.实验用MATLA函数介绍
在实验过程中,MATLA函数命令有数字滤波器函数
[N,Wn]=buttord(Wp,Ws,Rp,Rs)、模拟滤波器函数[N,Wn]=buttord(Wp,Ws,Rp,Rs,s’)。
式中Wp,Ws分别上通带和阻带的截止频
率,实际上它们是归一化频率,其值在0~1之间,Rp,Rs分别是通带和阻带的衰减,单位为dB。
N是求出的相应低通滤波器的阶次,Wn是求出的3dB频率,第二个式子的单位为rad/s,因此,它们实际上是频率。
用来设计模拟低通原型滤波器G(p),其调用格式是[z,p,k]=buttap(N)N是欲设计的低通原型滤波器的阶次,z、p和k分别是设计出的G(p)的极点、零点及增益。
以下4个文件用来将模拟低通原型滤波器G(p)分别转换为低通、高通、带通、及带阻滤波器。
其调用格式分别是
(1)[B,A]=lp2lp(b,a,Wo)或[B,A]=lp2hp(b,a,Wo)
(2)[B,A]=lp2bp(b,a,Wo,Bw)或[B,A]=lp2bs(b,a,Wo,Bw)
式中b,a分别是模拟低通原型滤波器G(p)有分子、分母多项式的系数向
量,B,A分别是转换后的Hs有分子、分母多项式的系数向量;
在格式
(1)中,Wo是低通或高通滤波器的截止频率;
在格式
(2)中Wo是带通或带阻滤波器的中心频率,Bw是其带宽。
实现双线性变换,即由模拟滤波器Hs得到数字滤波器Hz。
其调用格式是
[Bz,Az]=bilinear(B,A,Fs)
式中B、A分别是Hs的分子、分母多项式的系数向量;
Bz、Az分别是Hz的分子、分母多项式的系数向量,Fs是抽样频率。
5.思考题
(1)IIR滤波器设计步骤。
(2)对实验过程中所涉及的问题进行分析,试编写和修改相应的程序,得出
最终正确的结果和波形图,并对实验报告进行整理分析。
(3)设计低通数字滤波器,要求在通带内频带低于0.2rad时,允许幅度误差在1dB以内,在频率0.3rad~rad之间的阻带衰减大于15dB。
用双线性设计数字滤波器,T1,模拟滤波器采用巴特沃兹滤波器原型。
>
rp=1;
rs=15;
wp=0.2*pi;
ws=0.3*pi;
Fs=1;
wap=tan(wp/2);
was=tan(ws/2);
[n,wn]=buttord(wap,was,rp,rs,'
[z,p,k]=buttap(n);
[bp,ap]=zp2tf(z,p,k);
[bs,as]=lp2lp(bp,ap,wap);
[bz,az]=bilinear(bs,as,Fs/2);
[h,w]=freqz(bz,az,256);
plot(w,abs(h));
(4)设计一个巴特沃兹高通滤波器,要求通带截止频率为0.6,通带内衰减不大于1dB,阻带起始频率为0.4,阻带内衰减不小于15dB,T1。
(选做)
wp=0.6*pi;
ws=0.4*pi;
[bs,as]=lp2hp(bp,ap,wap);
6.实验报告要求
(1)明确实验目的以及实验的原理。
(2)通过实验内容掌握滤波器的设计。
(3)完成思考题的内容,对实验结果及其波形图进行分析,总结主要结论
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 IIR 数字滤波器 设计 滤波 思考题 程序 教学 文案