- 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
. . 一、设计目的 1. 掌握窗函数法设计FIR滤波器的原理和方法,观察用几种常用窗函数设计的FIR数字滤波器技术指标; 2. 掌握FIR滤波器的线性相位特性; 3. 了解各种窗函数对滤波特性的影响。 二、设计原理 如果所希望的滤波器的理想频率响应函数为Hd(ejω),则其对应的单位脉冲响应为,用窗函数wN(n)将hd(n)截断,并进行加权处理,得到实际滤波器的单位脉冲响应h(n)=hd(n)wN(n),其频率响应函数为。如果要求线性相位特性,则h(n)还必须满足。可根据具体情况选择h(n)的长度及对称性。 可以调用MATLAB工具箱函数fir1实现本实验所要求的线性相位FIR-DF的设计,调用一维快速傅立叶变换函数fft来计算滤波器的频率响应函数。 fir1是用窗函数法设计线性相位FIRDF的工具箱函数,调用格式如下: hn=fir1(N, wc, ‘ftype’, window) fir1实现线性相位FIR滤波器的标准窗函数法设计。 hn=fir1(N,wc)可得到6 dB截止频率为wc的N阶(单位脉冲响应h(n)长度为N+1)FIR低通滤波器,默认(缺省参数windows)选用hammiing窗。其单位脉冲响应h(n)满足线性相位条件:h(n)=h(N-1-n) 其中wc为对π归一化的数字频率,0≤wc≤1。 当wc=[wc1, wc2]时,得到的是带通滤波器。 hn=fir1(N,wc,’ftype’)可设计高通和带阻滤波器。 当ftype=high时,设计高通FIR滤波器; 当ftype=stop时,设计带阻FIR滤波器。 应当注意,在设计高通和带阻滤波器时,阶数N只能取偶数(h(n)长度N+1为奇数)。不过,当用户将N设置为奇数时,fir1会自动对N加1。 hn=fir1(N,wc,window)可以指定窗函数向量window。如果缺省window参数,则fir1默认为hamming窗。可用的其他窗函数有Boxcar, Hanning, Bartlett, Blackman, Kaiser和Chebwin窗。例如: hn=fir1(N,wc,bartlett(N+1))使用Bartlett窗设计; hn=fir1(N,wc,chebwin(N+1,R))使用Chebyshev窗设计。 hn=fir1(N,wc,’ftype’,window)通过选择wc、 ftype和window参数(含义同上),可以设计各种加窗滤波器。 详细设计步骤 1、用窗函数法设计一个线性相位FIR低通数字滤波器 FIR低通数字滤波器指标为: 因为衰减为50dB,所以选择海明窗。 过渡带宽为Ws-Wp=0.2π,由公式N>6.6π÷0.2π=33,所以 N=34。 所以程序如下: N=34; Wc=pi/5; %通带截止频率 wc=Wc/pi;%频率归一化 h=fir1(N,wc); [H,m]=freqz(h,[1],1024,whole); %频率响应 mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); subplot(2,2,1) n=0:N; stem(n,h,.) axis([0 N -0.1 0.3]) hold on n=0:N; x=zeros(N+1); plot(n,x,-) hold off xlabel(n) ylabel(h(n)) title(实际低通滤波器的h(n)) subplot(2,2,2) plot(m/pi,db) axis([0 1 -100 0]) xlabel(w/pi) ylabel(dB) title(副频衰减特性) grid on subplot(2,2,3) plot(m,pha) hold on n=0:7; x=zeros(8); plot(n,x,-) hold off axis([0 3.15 -4 4]) xlabel(频率(rad)) ylabel(相位(rad)) title(相频特性) subplot(2,2,4) plot(m,mag) axis([0 6.15 0 1.5]) xlabel(频率W(rad)) ylabel(幅值) title(幅频特性) 2、用窗函数法设计一个线性相位FIR高通数字滤波器。 要求: FIR高通数字滤波器指标为: 因为衰减为40dB,所以选择汉宁窗。 过渡带宽为Wp-Ws=0.2π,由公式N>6.2π÷0.2π=31,所以 N=32。 程序如下: N=32; Wc=pi/2; wc=Wc/pi;%频率归一化 h=fir1(N,wc, high, Hanning(N+1)); [H,m]=freqz(h,[1],1
文档评论(0)