- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验五 IIR滤波器设计
一、实验目的:
1、熟悉各种模拟原型滤波器的滤波特性;
2、掌握模数滤波器变换时的脉冲响应不变法和双线性变换法;
3、低通滤波器变换成其他类型滤波器的方法;
4、能够根据指标选择合适的原型滤波器和合适的方法设计IIR滤波器。
二、实验内容:
1、自定设计指标(通带截止频率、通带波纹、阻带截止频率、阻带衰减),选择合适的模拟原型低通滤波器和合适的设计方法(脉冲响应不变法或双线性变换法),设计符合要求的IIR滤波器;并验证设计好的滤波器是否满足要求。
2、将此低通滤波器映射为高通、带通或带阻滤波器,并验证变换后的滤波器是否满足指标。
3、求输入x(n)=[cos(w1n)+cos(w2n)]u(n)经过系统后的输出y(n)。其中w1位于通带内,w2位于阻带内,要求做一个两行两列的子图,第一个做x的时域波形,第二个做x的幅度谱,第三个做y的时域波形,第四个做y的幅度谱。
三、实验平台: MATLAB集成系统
四、四、实验要求
1、学习参考程序,理解程序内容,在重要程序段后作出文字说明;
2、比较写出设计过程与步骤;
3、查看相应函数的MATLAB函数的帮助文档,学会使用该函数。
参考程序:
窗口设计法:
ws1=0.2*pi;wp1=0.35*pi;wp2=0.65*pi;ws2=0.8*pi;As=60;
tr_width=min((wp1-ws1),(ws2-wp2));M=ceil(11*pi/tr_width)+1;
n=[0:1:M-1];wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
alpha=(M-1)/2;m1=n-alpha;fc1=wc1/pi;hd1=wc1*sinc(fc1*m1);
alpha=(M-1)/2;m2=n-alpha;fc2=wc2/pi;hd2=wc2*sinc(fc2*m2);
hd=hd1-hd2;
w_bla=(blackman(M));h=hd.*w_bla;b=h;a=[1];
[H,w]=freqz(b,a,1000,whole);
mag=abs(H);db=20*log10((mag+eps)/max(mag));
pha=angle(H);grd=grpdelay(b,a,w);
delta_w=2*pi/1000;
Rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))
As=-round(max(db(ws2/delta_w+1:1:501)))
subplot(2,2,1);stem(n,-hd);xlabel(n);ylabel(hd(n)); axis([0,M-1,-1.2,2]);
subplot(2,2,2);stem(n,w_bla);xlabel(n);ylabel(w(n)); axis([0,M-1,0,1.1]);
subplot(2,2,3);stem(n,-h);xlabel(n);ylabel(h(n)); axis([0,M-1, -1.2,2]);
subplot(2,2,4);plot(w/pi,db);xlabel(w);ylabel(db); axis([0,1, -139,10]);
接上面的程序:
ws1=0.2*pi;wp1=0.35*pi;wp2=0.65*pi;ws2=0.8*pi;As=60;
tr_width=min((wp1-ws1),(ws2-wp2));M=ceil(11*pi/tr_width)+1;
n=[0:1:M-1];wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
alpha=(M-1)/2;m1=n-alpha;fc1=wc1/pi;hd1=wc1*sinc(fc1*m1);
alpha=(M-1)/2;m2=n-alpha;fc2=wc2/pi;hd2=wc2*sinc(fc2*m2);
hd=hd1-hd2;
w_bla=(blackman(M));h=hd.*w_bla;
w1=0.5*pi;w2=0.1*pi;
x=cos(w1*n)+cos(w2*n);
nyb=n(1)+n(1);nye=n(length(x))+n(length(h));
ny=[nyb:nye];y=conv(x,h);
s=length(n);
w=0:pi/s:2*pi-pi/s;
k=-s:s-1;
X=x*(exp(-j*pi/s)).^(n*k);
N=length(ny);
g=0:pi/N:2*pi-pi/N;
p=-N:N-1;
Y=y*(exp(-j*pi/N)).^(ny*p);
magX=abs(X);magY=abs(Y);
subplot(2,2,1);stem(n,x);xla
文档评论(0)