2017年西电电院数字信号处理上机实验报告五

2017年西电电院数字信号处理上机实验报告五

ID:38974162

大小:423.30 KB

页数:13页

时间:2019-06-22

2017年西电电院数字信号处理上机实验报告五_第1页
2017年西电电院数字信号处理上机实验报告五_第2页
2017年西电电院数字信号处理上机实验报告五_第3页
2017年西电电院数字信号处理上机实验报告五_第4页
2017年西电电院数字信号处理上机实验报告五_第5页
资源描述:

《2017年西电电院数字信号处理上机实验报告五》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、数字信号处理上机报告实验五、IIR数字滤波器设计及其网络结构班级:学号:姓名:成绩:1实验目的(1)熟悉数字滤波的基本概念、数字滤波器的主要技术指标及其物理意义;(2)掌握巴特沃斯和切比雪夫模拟低通滤波器的设计方法和IIR数字低通滤波器的脉冲响应不变设计法、双线性变换法设计方法。(3)了解模拟和数字滤波器的频率变换、IIR数字滤波器的直接(优化)设计方法;2实验内容(1)设计计算机程序,根据滤波器的主要技术指标设计IIR数字巴特沃斯和切比雪夫低通、高通、带通和带阻滤波器;(2)绘制滤波器的幅频特性和相频特性曲线,验证滤波器的设计结果是否

2、达到设计指标要求;(3)画出数字滤波器的直接型、级联型、并联型网络结构信号流图。。3实验步骤(1)设计相应的八种滤波器的MATLAB程序;(2)画出幅频相频特性曲线;(3)画出信号流图。4程序设计%%巴特沃斯低通wp=0.2;ws=0.35;rp=1;rs=10;[N,wc]=buttord(wp,ws,rp,rs);[B,A]=butter(N,wc);w=0:0.001:pi;[H,w]=freqz(B,A,w);H1=20*log10(abs(H))subplot(2,1,1)plot(w/pi,H1),gridon;xlabel

3、('omega/pi'),ylabel('

4、H(e^i^omega)

5、')subplot(2,1,2)数字信号处理上机报告plot(w/pi,angle(H)/pi);xlabel('omega/pi'),ylabel('phi(omega)/pi')%%巴特沃斯高通wp=0.8;ws=0.6;rp=1;rs=10;[N,wc]=buttord(wp,ws,rp,rs);[B,A]=butter(N,wc,'high');w=0:0.001:pi;[H,w]=freqz(B,A,w);H1=20*log10(abs(H)

6、);subplot(2,1,1)plot(w/pi,H1),gridon;xlabel('omega/pi'),ylabel('

7、H(e^i^omega)

8、')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('omega/pi'),ylabel('phi(omega)/pi')%%巴特沃斯带通wpl=0.4;wpu=0.6;wsl=0.2;wsu=0.8wp=[wpl,wpu];ws=[wsl,wsu];rp=1;rs=20;[N,wc]=buttord(wp,ws,rp,rs);

9、[B,A]=butter(N,wc);w=0:0.001:pi;[H,w]=freqz(B,A,w);H1=20*log10(abs(H));subplot(2,1,1)plot(w/pi,H1),gridon;xlabel('omega/pi'),ylabel('

10、H(e^i^omega)

11、')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('omega/pi'),ylabel('phi(omega)/pi')%%巴特沃斯带阻wpl=0.2;wpu=0.8;wsl=0.4;wsu

12、=0.6wp=[wpl,wpu];ws=[wsl,wsu];rp=1;rs=20;[N,wc]=buttord(wp,ws,rp,rs);[B,A]=butter(N,wc,'stop');w=0:0.001:pi;[H,w]=freqz(B,A,w);H1=20*log10(abs(H));subplot(2,1,1)plot(w/pi,H1),gridon;xlabel('omega/pi'),ylabel('

13、H(e^i^omega)

14、')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel

15、('omega/pi'),ylabel('phi(omega)/pi')%%切比雪夫低通wp=0.2;ws=0.5;rp=1;rs=40;[N,wpo]=cheb1ord(wp,ws,rp,rs);[B,A]=cheby1(N,rp,wpo);w=0:0.001:pi;数字信号处理上机报告[H,w]=freqz(B,A,w);H1=20*log10(abs(H));subplot(2,1,1)plot(w/pi,H1),gridon;xlabel('omega/pi'),ylabel('

16、H(e^i^omega)

17、')s

18、ubplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('omega/pi'),ylabel('phi(omega)/pi')%%切比雪夫高通wp=0.7;ws=0.

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。