实验四 数字低通巴特沃斯滤波器的设计.doc

实验四 数字低通巴特沃斯滤波器的设计.doc

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

实验四 数字低通巴特沃斯滤波器的设计 1.掌握IIR数字滤波器的设计方法 IIR数字滤波器的设计,主要采用间接法,即:首先设计出低通模拟滤波器H(S);进行频率变换,将其转换为高通、带通、带阻滤波器;再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤波器。 对单极点的N阶H(S)用部分展开式: 冲激响应不变法取H(S)的单阶极点Sk的指数函数作为H(Z)的极点 双线性变换法是用代换H(S)中的S得到H(Z),双线性变换法可完全消除频率混叠失真但存在非线性频率失真,而冲激响应不变法存在混叠失真。 在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的滤波器设计函数,这些函数代表了不同类型的逼近函数的滤波器,常用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。 下图所示为IIR滤波器设计过程及所用到的MATLAB工具箱中的函数。 合成一步的设计函数 (butter,cheb1,cheb2,ellip,besself) 1.设计一个低通巴特沃斯模拟滤波器,绘制滤波器的的幅频响应及零极点图。指标如下: 通带截止频率:WP=1000HZ, 通带最大衰减:RP=3dB 阻带截止频率:Ws=2000HZ, 阻带最小衰减:Rs=40 dB 参考程序butter1.m 2. 用冲激响应不变法和双线性变换法将一模拟低通滤波器转换为数字低通滤波器 并图释H(S)和H(Z),采样频率Fs=1000Hz 参考程序butter2.m 1.列出本次实验编写的所有文件及各项实验结果的曲线,加注必要的说明 频率特性应包括幅频、相频特性 2.理论计算模拟低通原型滤波器的阶数和极点分布,说明实验中所用的冲激响应不变法和双线性变换法的原理及其编程根据 3.总结实验体会及实验中存在的问题 1.双线性变换法和冲激响应不变法比较有何优点? 巴特沃斯滤波器为例介绍设计函数的功能和用法,其它设计函数的用法类似。 (1)[N, Wn] = BUTTORD(Wp, Ws, Rp, Rs, 's') 功能:求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn 说明:Wp, Ws,为通带截止频率和阻带截止频率如式所示; , Rp, Rs为通带最大衰减 Apass(dB)和阻带最小衰减Astop(dB)如式 's'对应模拟滤波器。 (2)[Z,P,K] = BUTTAP(N) 功能:按给定的阶数N设计出巴特沃斯模拟低通滤波器 说明:Z,P,K为返回的N阶模拟滤波器的零点、极点和增益系数。如: [式3.5.4] 其中k为增益常数,zj和pi为系统函数H(Z)的N个零点和N个极点。 (3)模拟数字变换函数-双线性变换函数bilinear或脉冲响应不变法函数impinvar 调用格式:[NUMd,DENd] = BILINEAR(NUM,DEN,Fs) 功能: 把模拟滤波器系数为NUM,DEN变为近似等价的数字滤波器系数为NUMd,DENd 调用格式:[BZ,AZ] = IMPINVAR(B,A,Fs) 功能:设计出数字滤波器其单位抽样响应h(n)为模拟滤波器的单位冲激响应h(t)的采样值,采样频率为Fs 说明:B,A为模拟滤波器的系统函数H(S)的分子和分母多项式的系数,如下所示: B=[b1,b2,…,bm,bm+1],A=[a1,a2,…,an,an+1] (4) 合为一步的数字滤波器设计函数butter 调用格式1:[B,A] = BUTTER(N,Wn) 功能:设计N阶低通数字滤波器,系数矢量B,A按Z的降幂排列,当Wn为二元矢量时,Wn = [W1 W2],返回一个2N阶带通滤波器,通带范围:W1 < W < W2,如式3.5.6 调用格式2:[B,A] = BUTTER(N, Wn, 'high') 设计高通数字滤波器. 调用格式3:[B,A] = BUTTER(N, Wn, 'stop') 设计带阻滤波器 Wn = [W1,W2]. 调用格式4:[Z,P,K] = BUTTER(N,Wn) 返回低通数字滤波器的N阶零、极点矢量. FREQS 调用格式: [H,W] = FREQS(B,A,M) 功能:此函数可以求出系统频率响应的数值解,并列出此系统的幅频及相频响应曲线。 1.butter1.m %巴特沃兹滤波器的幅频响应图 subplot(1,2,1); %分两个窗口,幅频图在第一个窗口 wp=1000;ws=2000;rp=3;rs=40; %设置指标 [N,wn]=buttord(wp,ws,rp,rs,'s') %计算巴特沃斯低通滤波器的阶数和3dB截止频率 [B,A]=butter(N,wn,'s'); %代入N和Wn设计巴特沃斯模拟低通滤波器 [Z,P,K]=b

文档评论(0)

xcs88858 + 关注
实名认证
内容提供者

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

版权声明书
用户编号:8130065136000003

1亿VIP精品文档

相关文档