




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、用双线性变换法设计切比雪夫II型的数字IIR带通滤波器用双线性变换法设计原型低通为切比雪夫II型的数字IIR带通滤波器,要求通带边界频率为400Hz,500Hz,阻带边界频率分别为350Hz,550Hz,通带最大衰减1dB,阻带最小衰减40dB,抽样频率为2000Hz,用MATLAB画出幅频特性,画出并分析滤波器传输函数的零极点;信号经过该滤波器,其中450Hz,600Hz,滤波器的输出是什么?用Matlab验证你的结论并给出的图形。Matlab 详细设计:% Design of a Cheb II Bandpass Digital Filter by using bilinear metho
2、dclc;clear all;Rp = 1; % bandpass attenuation in dBRs = 40; % bandstop attenuation in dBOmegaS1_1=350; OmegaS1_2=550;OmegaP1_1=400; OmegaP1_2=500;Fp=2000; % samling frequencyWp1=2*pi*OmegaP1_1/Fp; Wp2=2*pi*OmegaP1_2/Fp; Ws1=2*pi*OmegaS1_1/Fp; Ws2=2*pi*OmegaS1_2/Fp; OmegaP1=2*Fp*tan(Wp1/2); % nonline
3、arlizationOmegaP2=2*Fp*tan(Wp2/2); % nonlinearlizationOmegaS1=2*Fp*tan(Ws1/2); % nonlinearlizationOmegaS2=2*Fp*tan(Ws2/2); % nonlinearlizationOmegaP0=sqrt(OmegaP1*OmegaP2);% equivalent mid frequencyBw=OmegaP2-OmegaP1; % bandwithEta_P0=OmegaP0/Bw; % NormalizationEta_P1=OmegaP1/Bw; % NormalizationEta_
4、P2=OmegaP2/Bw; % NormalizationEta_S1=OmegaS1/Bw; % NormalizationEta_S2=OmegaS2/Bw; % NormalizationLemta_P_EquivalentLowPass=Eta_P2/(Eta_P22-Eta_P02); % change to the equivalent Lowpass patameterLemta_S1_EquivalentLowPass=-Eta_S1/(Eta_S12-Eta_P02); % change to the equivalent Lowpass patameterLemta_S2
5、_EquivalentLowPass=Eta_S2/(Eta_S22-Eta_P02); % change to the equivalent Lowpass patameterLemta_S_EquivalentLowPass=min(Lemta_S1_EquivalentLowPass,Lemta_S2_EquivalentLowPass); % get the smallest % Estimate the Filter OrderN, Wn=cheb2ord(Lemta_P_EquivalentLowPass, Lemta_S_EquivalentLowPass, Rp, Rs,s);
6、 % Design the Filternum1,den1=cheby2(N,Rs,Wn,s); num2,den2=lp2bp(num1,den1,OmegaP0,Bw);num,den=bilinear(num2,den2,Fp);% Compute the gain responsew = 0:pi/255:pi;h = freqz(num,den,w);g = 20*log10(abs(h);% Plot the gain responsefigure;plot(w/pi,g);grid axis(0 1 -60 5);xlabel(omega /pi); ylabel(Gain in
7、 dB);title(Gain Response of a Cheb II Bandpass Filter);%Plot the poles and zerosz,p,k=tf2zp(num,den);figure;zplane(z,p); %title()f1=450;f2=600;t=0:0.0001:1x1=sin(2*pi*f1*t);x2=sin(2*pi*f2*t);x=x1+x2;figure;subplot(2,2,1)%x1plot(x1);grid on;axis(0,50*pi,-3,3);xlabel(t);ylabel(x1(t);title(x1);subplot(2,2,2)%x1plot(x2);grid on;axis(0,50*pi,-3,3);xlabel(t);ylabel(x2(t);title(x2);subplot(2,2,3)%xplot(x);grid on;axis(0,50*pi,-3,3);xlabel(t);ylabel(x(t);title(x)%X=fft(x);y=filter(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 牛奶公司年会活动方案
- 生鲜柜促销活动方案
- 珍珠钻石活动方案
- 玩具共享活动方案
- 环保督察活动方案
- 瓷砖年底活动方案
- 生物安全调研活动方案
- 物交会志愿服务活动方案
- 爱心涌动活动方案
- 物资采购工会活动方案
- 2025届黑龙江省哈尔滨四十七中学七年级英语第二学期期末统考试题含答案
- 译林版(2024)七年级下册英语期末复习:完形填空+阅读理解 练习题(含答案)
- 第5章 相交线与平行线 复习课件
- 机械租赁投标服务方案
- 《工程勘察设计收费标准》(2002年修订本)
- 人工智能知到章节答案智慧树2023年复旦大学
- 建筑安全生产自查台账(建筑施工)
- 人教版 小学音乐下册 一至六年级全套精品教案(1-6年级全套合集)
- 承台模板及支架检验批质量验收记录表
- 完整word版ZEMAX光学设计超级学习手册第1章
- 人教版小学数学三年级下册期末复习计算题专项训练—口算与估算【含答案】
评论
0/150
提交评论