




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、设计流程图如下: 设计思想: 首先设计一个源信号和一个混合信号,通过其频谱对比得出最大和最小通 带,最大和最小阻带 ; 然后再根据得到的参数来设计切比雪夫滤波器, 最后通过 切比雪夫型滤波器和切比雪夫型滤波器的对比来得出那种效果好。 切比雪夫滤波器设计原理: 切比雪夫滤波器的振幅具有等波纹特性,它有两种形式: 1) 振幅特性在通 带内是等波纹的、在阻带内是单调的切比雪夫 I 型滤波器; 2) 振幅特性在通带 内是单调的、在阻带内是等波纹的切比雪夫 II 型滤波器,采用何种形式的切比 雪夫滤波器取决于实际用途 . 切比雪夫滤波器的设计方法就是将逼近精确度均匀分布在整个通带内,或 者均匀分布在整个
2、阻带内,或者均匀分布在两者之内,这样就可以使滤波器阶 数大大降低。 切比雪夫I 型滤波器平方幅度响应函数表示为 : G( j ) 2=1+ 2C2N( ) 1/2 其中 1. 而切比雪夫 II 型滤波器平方幅度响应函数表示为: G(j ) 2=1+ 2 C 2N( )/ 2N( / c) 2 1 其中 1(正数) ,表示波纹变化情况; c为截止频率; N为滤波器的阶次 , 也是 CN ( / N) 的阶次。 源信号编码及其图形: t=-1:1 y=(cos(2*pi*10*t)+cos(2*pi*40*t); N=length(y); fx=fft(y); df=100/N; n=0:N/2;
3、 f=n*df; subplot(2,1,1); plot(f,abs(fx(n+1)*2/N); grid; title( 源波形频谱 ) 图(一) 混合信号编码及其图形: t=-1:1; X=(cos(2*pi*10*t)+cos(2*pi*25*t)+cos(2*pi*40*t); N=length(X); fx=fft(X); df=100/N; n=0:N/2; f=n*df; subplot(2,1,2); plot(f,abs(fx(n+1)*2/N); grid; title( 混合波形频谱 ) 图(二) 从图(一)和图(二)对比可以得出:为了能达到和满足我们的要求,我们取以
4、下的参数,最大通带 wp2:,最小通带 wp1:,最大阻带 ws2: ,最小阻带 ws1: 切比雪夫型滤波器设计如下: ws1=*pi;ws2=*pi; % 滤波器的阻带截止频率 wp1=*pi;wp2=*pi; % 滤波器的通带截止频率 Rp=1;As=20; % 滤波器的通阻带衰减指标 %转换为模拟滤波器的技术指标 T=;Fs=1/T; Omgp1=(2/T)*tan(wp1/2);Omgp2=(2/T)*tan(wp2/2); Omgp=Omgp1,Omgp2; Omgs1=(2/T)*tan(ws1/2);Omgs2=(2/T)*tan(ws2/2); Omgs=Omgs1,Omgs2
5、; bw=Omgp2-Omgp1;w0=sqrt(Omgp1*Omgp2);模 %拟通带带宽和中心频率 ripple=10(-Rp/20); % 滤波器的通带衰减对应的幅度值 Attn=10(-As/20); % 滤波器的阻带衰减对应的幅度值 %模拟原型滤波器计算 n,Omgn=cheb1ord(Omgp,Omgs,Rp,As,s) % 计算阶数 n 和截止频率 z0,p0,k0=cheb1ap(n,Rp); % 设计归一化的模拟滤波器原型 ba1=k0*real(poly(z0); % 求原型滤波器的系数 b aa1=real(poly(p0); % 求原型滤波器的系数 a ba,aa=lp
6、2bs(ba1,aa1,w0,bw); %用双线性变换法计算数字滤波器系数 bd,ad=bilinear(ba,aa,Fs) %求数字系统的频率特性 H,w=freqz(bd,ad); dbH=20*log10(abs(H)+eps)/max(abs(H); subplot(2,2,1);plot(w/pi,abs(H); ylabel(|H|);xlabel(频率(pi);title(幅度响应 );axis(0,1,0,); set(gca,XTickMode,manual,XTick,0,); set(gca,YTickMode,manual,YTick,0,Attn,ripple,1);
7、grid 图(三) n =3 Omgn = bd = ad = 分析:由图(三)运行结果可知,最大通带,最小通带,最大阻带,最小阻带 ; 切比雪夫型滤波器的设计的个性技术指标精确度是均匀分布的。 而其幅度特性 在通带内是等波纹的,在阻带内是单调下降的。 虽然达到了所设计滤波器的要求, 滤除了我们不需要的波形,但是存在一些问题,如误差等。 切比雪夫型滤波器设计如下 : ws1=*pi;ws2=*pi; % 滤波器的阻带截止频率 wp1=*pi;wp2=*pi; % 滤波器的通带截止频率 Rp=1;As=20; % 滤波器的通阻带衰减指标 %转换为模拟滤波器的技术指标 T=;Fs=1/T; Omg
8、p1=(2/T)*tan(wp1/2);Omgp2=(2/T)*tan(wp2/2); Omgp=Omgp1,Omgp2; Omgs1=(2/T)*tan(ws1/2);Omgs2=(2/T)*tan(ws2/2); Omgs=Omgs1,Omgs2; bw=Omgp2-Omgp1;w0=sqrt(Omgp1*Omgp2);模 %拟通带带宽和中心频率 ripple=10(-Rp/20); % 滤波器的通带衰减对应的幅度值 Attn=10(-As/20); % 滤波器的阻带衰减对应的幅度值 %模拟原型滤波器计算 n,Omgn=cheb2ord(Omgp,Omgs,Rp,As,s) z0,p0,k
9、0=cheb2ap(n,As); % 设计归一化的 cheb2 型模拟滤波器原型 ba1=k0*real(poly(z0); %求原型滤波器的系数 b aa1=real(poly(p0); %求原型滤波器的系数 a ba,aa=lp2bs(ba1,aa1,w0,bw); % 变换为模拟带通滤波器 %用双线性变换法计算数字滤波器系数 bd,ad=bilinear(ba,aa,Fs) %求数字系统的频率特性 H,w=freqz(bd,ad); dbH=20*log10(abs(H)+eps)/max(abs(H); subplot(2,2,1);plot(w/pi,abs(H); ylabel(|
10、H|);xlabel(频率(pi);title(幅度响应 );axis(0,1,0,); set(gca,XTickMode,manual,XTick,0,); set(gca,YTickMode,manual,YTick,0,Attn,ripple,1);grid 图(四) n =3 Omgn = bd = ad = 分析:由图(四)的运行结果可知, 最大通带,最小通带,最大阻带,最小阻带 ; 所设计的切比雪夫型滤波器的各项技术指标精确度是均匀分布的。 而其幅度特 性在通带内是单调的, 在阻带内是等波纹的。虽然也达到了所设计滤波器的要求, 滤除了我们不需要的并行,但是也还是存在误差,这是可以理解的。 由图(三)和图(四)的运行
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年中国木麻黄栲胶项目商业计划书
- 哈尔滨市中医院心脏功能测定考核
- 通辽市中医院溶栓药物不良反应处理考核
- 中国四丁基四氟硼酸铵项目商业计划书
- 天津市中医院老年疾病影像考核
- 运城市人民医院肝静脉成形技术考核
- 邯郸市人民医院现金流量表编制与分析试题
- 唐山市人民医院文化建设在人事工作中的落地试题
- 大庆市中医院肿瘤遗传咨询与家族风险评估考核
- 中国紫外线涂料项目经营分析报告
- 预防医学考试题+答案
- 二年级上册数学北师大版课件第5课时 小熊开店
- 跌倒坠床原因分析预防措施
- 52206马工程组织行为学课件
- 我和我的祖国课件
- 各类食物营养与配餐(蛋类的营养)课件
- 公司内账管理系统
- 全国细菌耐药监测网信息系统-附件
- 妇产科产前诊断技术服务临床医师考核题(附答案)
- 校园欺凌工作台账(完整资料)
- DB33∕T 1146-2018 浙江省城市轨道交通规范
评论
0/150
提交评论