资源描述:
《柱子信号处理滤波器的设计》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验三:滤波器的设计第五题:fp=6000;fs=12000;ap=3;as=25;wp=2*pi*fp;ws=2*pi*fs;[N,WC]=buttord(wp,ws,ap,as,'s')[B,A]=butter(N,WC,'s')x=20*log10(abs(H));y=angle(H);f=0:20;w=2*pi*f;H=freqs(B,A,w);subplot(2,1,1);plot(f,x)axis([0,20,-10.0e-16,1.0e-16]);xlabel('f/HZ');ylabel(
2、'幅度/dB');gridonsubplot(2,1,2);plot(f,y)gridonaxis([0,5,-9.8e-04,0]);xlabel('f/HZ');ylabel('相位');第六题:clear,clc,closeallfp=12000;fs=6000;ap=1;as=40;wp=2*pi*fp;ws=2*pi*fs;[N,WC]=cheb1ord(wp,ws,ap,as,'s')[B,A]=cheby1(N,ap,WC,'high','s')f=2000:18000;w=2*pi*f;H
3、=freqs(B,A,w);x=20*log10(abs(H));y=angle(H);subplot(2,1,1);plot(f,x)axis([2000,18000,-50,2]);xlabel('f/HZ');ylabel('·幅度/dB');gridonsubplot(2,1,2);plot(f,y)gridonaxis([2000,18000,-5,5]);xlabel('f/HZ');ylabel('相位');第七题:clear,clc,closeallfp=[20000,40000];fs=
4、[10000,50000];ap=1;as=40;wp=2*pi*fp;ws=2*pi*fs;[N,WC]=ellipord(wp,ws,ap,as,'s')[B,A]=ellip(N,ap,as,WC,'s')f=10000:60000;w=2*pi*f;H=freqs(B,A,w);x=20*log10(abs(H));y=angle(H);subplot(2,1,1);plot(f,x)axis([10000,60000,-100,100]);xlabel('f/HZ');ylabel('·幅度/d
5、B');gridonsubplot(2,1,2);plot(f,y)gridonaxis([10000,60000,0.2:0.01:0.3]);xlabel('f/HZ');ylabel('相位');第八题:clear,clc,closeallfp=[20000,40000];fs=[10000,50000];ap=1;as=40;wp=2*pi*fp;ws=2*pi*fs;[N,WC]=buttord(wp,ws,ap,as,'s')[B,A]=butter(N,WC,'s')f=10000:6000
6、0;w=2*pi*f;H=freqs(B,A,w);x=20*log10(abs(H));y=angle(H);subplot(2,1,1);plot(f,x)axis([10000,60000,-110,110]);xlabel('f/HZ');ylabel('·幅度/dB');gridonsubplot(2,1,2);plot(f,y)gridonaxis([10000,60000,-4,4]);xlabel('f/HZ');ylabel('相位');第九题:ap=0.5;as=50;wp=0.4;w
7、s=0.6;[N,WC]=buttord(wp,ws,ap,as);[B,A]=butter(N,WC);f=10:40;w=2*pi*f;H=freqz(B,A,w);[hn,n]=impz(B,A,20);y=20*log10(abs(H));subplot(3,1,1);stem(f,y)axis([10,20,-7e-15,-4e-15]);xlabel('f/HZ');ylabel('幅度/dB');gridonsubplot(3,1,2);stem(f,angle(H))xlabel('f/H
8、Z');ylabel('相位');gridonsubplot(3,1,3);stem(n,hn)xlabel('f/HZ');ylabel('单位脉冲响应');gridon第十题:ap=1;as=40;wp=[0.25,0.75];ws=[0.35,0.65];[N,WC]=cheb1ord(wp,ws,ap,as);[B,A]=cheby1(N,ap,WC,'stop');t=10:40;w=2*pi*t;H=freqz(B