MATLAB实现数字FIR高通带通等滤波器源程序31471.docVIP

MATLAB实现数字FIR高通带通等滤波器源程序31471.doc

  1. 1、本文档共12页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
利用汉宁窗设计型数字高通滤波器 clear all; Wp=0.6*pi; Ws=0.4*pi; tr_width=Wp-Ws;???????????????????????? %过渡带宽度 N=ceil(6.2*pi/tr_width)??????????????????????? %滤波器长度 n=0:1:N-1; Wc=(Ws+Wp)/2;???????????????????????? %理想低通滤波器的截止频率 hd=ideal_hp1(Wc,N);????????????????????? %理想低通滤波器的单位冲激响应 w_han=(hanning(N));???????????????????? %汉宁窗 h=hd.*w_han;??????????????????????????? %截取得到实际的单位脉冲响应 [db,mag,pha,w]=freqz_m2(h,[1]);????????????? %计算实际滤波器的幅度响应 delta_w=2*pi/1000; Ap=-(min(db(Wp/delta_w+1:1:501)))???????? %实际通带纹波 As=-round(max(db(1:1:Ws/delta_w+1)))???????? %实际阻带纹波 subplot(221) stem(n,hd) title(理想单位脉冲响应hd(n)) subplot(222) stem(n,w_han) title(汉宁窗w(n)) subplot(223) stem(n,h) title(实际单位脉冲响应hd(n)) subplot(224) plot(w/pi,db) title(幅度响应(dB)) axis([0,1,-100,10]) clear all; Wp=0.6*pi; Ws=0.4*pi; tr_width=Wp-Ws;???????????????????????? %过渡带宽度 N=ceil(6.2*pi/tr_width)??????????????????????? %滤波器长度 n=0:1:N-1; Wc=(Ws+Wp)/2;???????????????????????? %理想低通滤波器的截止频率 hd=ideal_hp1(Wc,N);????????????????????? %理想低通滤波器的单位冲激响应 w_han=(hanning(N));???????????????????? %汉宁窗 h=hd.*w_han;??????????????????????????? %截取得到实际的单位脉冲响应 [db,mag,pha,w]=freqz_m2(h,[1]);????????????? %计算实际滤波器的幅度响应 delta_w=2*pi/1000; Ap=-(min(db(Wp/delta_w+1:1:501)))???????? %实际通带纹波 As=-round(max(db(1:1:Ws/delta_w+1)))???????? %实际阻带纹波 subplot(221) stem(n,hd) title(理想单位脉冲响应hd(n)) subplot(222) stem(n,w_han) title(汉宁窗w(n)) subplot(223) stem(n,h) title(实际单位脉冲响应hd(n)) subplot(224) plot(w/pi,db) title(幅度响应(dB)) axis([0,1,-100,10])基于切比雪夫一致逼近法设计FIR数字低通滤波器 clear all;?????????????? f=[0 0.6 0.7 1];???????????????????? %给定频率轴分点 A=[1 1 0 0];???????????????????? %给定在这些频率分点上理想的幅频响应 weigh=[1 10];????????????????? %给定在这些频率分点上的加权 b=remez(32,f,A,weigh);??????????? %设计出切比雪夫最佳一致逼近滤波器 [h,w]=freqz(b,1,256,1); h=abs(h); h=20*log10(h); subplot(211) stem(b,.); grid; title(切比雪夫逼近滤波器的抽样值) subplot(212) plot(w,h); grid; title(滤波器幅频特性(dB)) 利用汉宁窗设计型数字带阻滤波器 clear all;???????????????????? Wpl=0.2*pi; Wph=0.8*pi; Wsl=0.4*pi; Wsh=0.6*pi; tr_width=min((Wsl-Wpl),(Wph-Wsh));????

文档评论(0)

ipad0c + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档