




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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年企业首席质量官培训考核试题(含答案)
- 《离骚》拼音及注释
- 宝钢设备大修管理办法
- (高清版)DB46∕T 707-2025 榴莲栽培技术规程
- 成人阻塞性睡眠呼吸暂停多学科诊疗指南主要内容
- 迁地种群遗传多样性-洞察及研究
- Q-CSG1211016-2025 光伏发电站接入电网技术规范
- 社保费培训课件税务局
- 2025-2030年古董行业市场深度调研及发展趋势与投资研究报告
- 《ISO 37001-2025 反贿赂管理体系要求及使用指南》专业深度解读和应用培训指导材料之1:2范围+3术语和定义(雷泽佳编制-2025A1)
- 音乐节与音乐会策划项目可行性分析报告
评论
0/150
提交评论