IIR数字滤波器讲述.ppt

下载文档 降价啦

3
0
约 15页
2017-03-04 发布于湖北
举报
版权申诉
保障服务

IIR数字滤波器讲述.ppt

1、本文档共15页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

IIR数字滤波器讲述

MATLAB实现 * 1、matlab中的交互式图形用户界面工具—sptool 导入信号： [x,fs,bits]=wavread(lx.wav,[60000 120000]); sound(x,fs,bits); x=x(:,1); 利用sptool工具箱进行滤波去噪处理 sptool MATLAB实现 * MATLAB实现 * MATLAB实现 * MATLAB实现 * MATLAB实现 * 2、matlab中的滤波器直接设计IIR滤波器的函数 freqz函数:求解数字滤波器频率响应 [h, f]=freqz{b, a, n, fsy=abs buttord函数:选择巴特沃斯滤波器最小阶数 [n, wn]=buttord(wp, ws, rp, rs) Cheblord函数:选择切笔雪夫I型滤波器最小阶数 [n,wn]= Cheb1ord （wp, ws, rp, rs） Cheb2ord函数:选择切笔雪夫II型滤波器最小阶数 [n,wn]= Cheb2ord （wp, ws, rp, rs） Ellipard函数:选择椭圆滤波器最小阶数 [n,wn]=Ellipord(wp, ws, rp, rs) MATLAB实现 * N = length(x); t=1/fs*(0:N-1) f=fs/N*(-N/2:N/2-1); w = 2*f/fs; X=fftshift(fft(x)); magX=abs(X); angX=angle(X); figure(1) plot(t,x);title(原始信号波形); figure(2) plot(f,magX);title(原始信号幅值);axis([-5000 5000 -inf inf]);xlabel(频率); figure(3) plot(w,magX);title(原始信号幅值);axis([-1 1 -inf inf]);xlabel(数字角频率); 画出原信号波形及频谱图 MATLAB实现 * MATLAB实现 * 程序设计切笔雪夫II型滤波器 wp=0.1; ws=0.2; Rp=1; Rs=80; [N,wc]=cheb2ord(wp,ws,Rp,Rs); [B,A]=cheby2(N,Rs,wc); figure(4) freqz(B,A); MATLAB实现 * 进行滤波 y = filter(B,A,x); Y=fftshift(fft(y)); magY=abs(Y); angY=angle(Y); figure(5) plot(t,y);title(原始信号波形); figure(6) plot(f,magY);title(原始信号幅值);axis([-5000 5000 -inf inf]);xlabel(频率); figure(7) plot(w,magY);title(原始信号幅值);axis([-1 1 -inf inf]);xlabel(数字角频率); MATLAB实现 * MATLAB实现 * 滤波前后比较 MATLAB实现 * 滤波前后比较原信号滤波后 Wavwrite(x,fs,bits); Wavwrite(y,fs,bits); 总结 * 利用MATLAB基本实现了去除音频文件中的高频部分。程序仿真中还有待于进一步优化。谢谢观看！ *