ex19模拟域频率变换法试题.docVIP

数字信号处理实验 第 十一 次实验 实验名称： 模拟域频率变换法学生班级： 电信学生姓名：学生学号：指导教师： zgx实验目的 加深对模拟域频率变换法的了解 掌握使用模拟低通滤波器原型进行频率变换及设计低通、高通、带通、带阻滤波器的方法 了解MATLAB有关模拟域频率变换的子函数及其使用方法 实验任务 用频率变换法设计一个切比雪夫二型模拟低通滤波器，要求通带截止频率fp3.5kHz，通带最大衰减Rp≤1dB，阻带截止频率fs6kHz，阻带最小衰减As≥40dB。绘制归一化的模拟滤波器原型和实际的模拟低通滤波器的频率特性。 用频率变换法设计一个切比雪夫二型模拟高通滤波器，要求通带截止频率fp6kHz，通带最大衰减Rp≤1dB，阻带截止频率fs3.5kHz，阻带最小衰减As≥40dB。绘制归一化的模拟滤波器原型和实际的模拟高通滤波器的频率特性。 用频率变换法设计一个椭圆模拟带通滤波器，要求通带截止频率fp13.5kHz，fp25.5kHz，通带最大衰减Rp≤1dB，阻带下截止频率fs13kHz，阻带上截止频率fs6kHz，阻带最小衰减As≥40dB。绘制归一化的模拟滤波器原型和实际的模拟带通滤波器的频率特性。 用频率变换法设计一个切比雪夫一型模拟带阻滤波器，要求通带下截止频率fp13kHz，通带上截止频率fp27kHz，通带最大衰减Rp≤1dB，阻带下截止频率fs14kHz，阻带上截止频率fs6kHz，阻带最小衰减As≥35dB。绘制归一化的模拟滤波器原型和实际的模拟阻带滤波器的频率特性。 实验原理 模拟域频率变换法 MATLAB提供的lp2lp等子函数能将模拟滤波器原型转换到实际的模拟滤波器中 实验步骤及结果 实验程序清单： fp3500;Rp1; fs6000;As40; wp2*pi*fp; ws2*pi*fs; [n,wc]cheb2ordwp,ws,Rp,As,s [z0,p0,k0]cheb2apn,As b0k0*realpolyz0; a0realpolyp0; [H,w0]freqsb0,a0; dbH20*log10absH+eps/maxabsH; [b1,a1]lp2lpb0,a0,wc; [Ha,wa]freqsb1,a1; dbHa20*log10absHa+eps/maxabsHa; wp0fp/wc; wc0wc/2/pi/wc; ws0fs/wc; fcfloorwc/2/pi; subplot2,2,1,plotw0/2*pi,dbH;grid title归一化模拟原型幅度;ylabeldB; subplot2,2,2,plotw0/2*pi,angleH/pi*180;grid title归一化模拟原型相位;ylabel\phi; subplot2,2,3,plotwa/2*pi,dbHa;grid title实际模拟低通幅度;ylabeldB;xlabel频率Hz; subplot2,2,4,plotwa/2*pi,angleHa/pi*180;grid title实际模拟低通相位;ylabel\phi;xlabel频率Hz; 运行结果如图1所示： 图 1fs3500;As40; wp2*pi*fp;ws2*pi*fs; [n,wc]cheb2ordwp,ws,Rp,As,s [z0,p0,k0]cheb2apn,As b0k0*realpolyz0; a0realpolyp0; [H,w0]freqsb0,a0; dbH20*log10absH+eps/maxabsH; [b1,a1]lp2hpb0,a0,wc; [Ha,wa]freqsb1,a1; dbHa20*log10absHa+eps/maxabsHa; wp0fp/wc;wc0wc/2/pi/wc; ws0fs/wc;fcfloorwc/2/pi; subplot2,2,1,plotw0/2*pi,dbH;grid title归一化模拟原型幅度;ylabeldB; subplot2,2,2,plotw0/2*pi,angleH/pi*180;grid title归一化模拟原型相位;ylabel\phi; subplot2,2,3,plotwa/2*pi,dbHa;grid title实际模拟高通幅度;ylabeldB;xlabel频率Hz; subplot2,2,4,plotwa/2*pi,angleHa/pi*180;grid title实际模拟高通相位;ylabel\phi;xlabel频率Hz;运行结果如图2所示： 图 2 fs13000;fs26000; As40;Rp1; wp12*pi*fp1;wp22*pi*fp2; ws12*pi*fs1;ws22*pi*fs2; wp