实验一 连续时间系统的时域和频域分析.docVIP

信息科学与工程学院信号与系统实验报告 专业班级 电信0802班 姓 名 殷 玮 学 号 200804135052 实验时间 指导教师 成 绩 实验名称 实验三 模拟滤波器及IIR数字滤波器的设计 实 验 题 目 、 源 程 序 、 实 验 结 果 设计一个巴特沃斯模拟低通滤波器，以满足：通带截止频率，通带最大衰减，阻带截止频率，阻带最小衰减。要求绘出滤波器的幅频特性曲线。（幅度用分贝值表示） fp=5; wp=2*pi*fp; fs=12; ws=2*pi*fs; rp=2; rs=30; [N,Wn]=buttord(wp,ws,rp,rs,s); [b,a]=butter(N,Wn,s); freqs(b,a) 设计一个巴特沃斯模拟高通滤波器，以满足：通带截止频率，通带最大衰减，阻带截止频率，阻带最小衰减。要求绘出滤波器的幅频特性曲线。（幅度用分贝值表示） fp=20; wp=2*pi*fp; fs=10; ws=2*pi*fs; rp=3; rs=15; [N,omegac]=buttord(wp,ws,rp,rs,s); [z0,p0,k0]=buttap(N); b0=k0*real(poly(z0)); a0=real(poly(p0)); [H,omega0]=freqs(b0,a0); dbH=20*log10((abs(H)+eps)/max(abs(H))); [b1,a1]=lp2hp(b0,a0,omegac); [H1,omega1]=freqs(b1,a1); dbH1=20*log10((abs(H1)+eps)/max(abs(H1))); subplot(221),plot(omega0/(2*pi),dbH);grid axis([0,1,-50,1]);title(归一化模拟低通原型幅度响应); ylabel(dB); subplot(222),plot(omega0/(2*pi),angle(H));grid axis([0,1,-4,4]);title(归一化模拟低通原型相位响应); ylabel(弧度); subplot(223),plot(omega1/(2*pi),dbH1);grid axis([0,2*fs,-50,1]);title(实际模拟高通幅度响应); ylabel(dB);xlabel(频率(HZ)); subplot(224),plot(omega1/(2*pi),angle(H1));grid axis([0,2*fs,-4,4]);title(实际模拟高通相位响应); ylabel(弧度);xlabel(频率(HZ)); 设计一个巴特沃斯模拟带通滤波器，以满足：通带范围为10Hz～25Hz，阻带截止频率分别为5Hz、30Hz，通带最大衰减为3dB，阻带最小衰减为30dB。要求绘出滤波器的幅频特性曲线。（幅度用分贝值表示） fp1=10;wp1=2*pi*fp1; fp2=25;wp2=2*pi*fp2; omegap=[wp1 wp2]; fs1=5;ws1=2*pi*fs1; fs2=30;ws2=2*pi*fs2; omegas=[ws1 ws2]; bw=wp2-wp1;w0=sqrt(ws1*ws2); rp=3;rs=30; [N,omegac]=buttord(omegap,omegas,rp,rs,s); [z0,p0,k0]=buttap(N); b0=k0*real(poly(z0)); a0=real(poly(p0)); [H,omega0]=freqs(b0,a0); dbH=20*log10((abs(H)+eps)/max(abs(H))); [b1,a1]=lp2bp(b0,a0,w0,bw); [H1,omega1]=freqs(b1,a1); abH1=20*log10((abs(H1)+eps)/max(abs(H1))); subplot(221),plot(omega0/(2*pi),dbH);grid axis([0,1,-50,1]);title(归一化模拟低通原型幅度响应); ylabel(dB); subplot(222),plot(omega0/(2*pi),angle(H));grid axis([0,1,-4,4]);title(归一化模拟低通原型相位响应); ylabel(弧度); subplot(223),plot(omega1/(2*pi),dbH1);grid axis([0,40,-50,1]);title(实际模拟带通幅度响应); yl