资源描述：
《Matlab Simulink建模与仿真例题源代码》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、〔實例3.1〕試設計一個模擬低通濾波器，fp=2400Hz，fs=5000Hz，Rp=3dB，Rs=25dB。分別用巴特沃斯和橢圓濾波器原型，求出其3dB截止頻率和濾波器階數，傳遞函數，並作出幅頻、相頻特性曲線。巴特沃斯濾波器設計的程序代碼如下：%ch3example1A.mclear;f_p=2400;f_s=5000;R_p=3;R_s=25;%設計要求指標[n,fn]=buttord(f_p,f_s,R_p,R_s,'s');%計算階數和截止頻率Wn=2*pi*fn;%轉換為角頻率[b,a]=b
2、utter(n,Wn,'s');%計算H(s)f=0:100:10000;%計算頻率點和頻率範圍s=j*2*pi*f;%s=jw=j*2*pi*fH_s=polyval(b,s)./polyval(a,s);%計算相應頻率點處H(s)的值figure(1);subplot(2,1,1);plot(f,20*log10(abs(H_s)));%幅頻特性axis([010000-401]);xlabel('頻率Hz');ylabel('幅度dB');subplot(2,1,2);plot(f,angle(
3、H_s));%相頻特性xlabel('頻率Hz');ylabel('相角rad');figure(2);freqs(b,a);%也可用指令freqs直接畫出H(s)的頻率響應曲線。橢圓濾波器設計的程序代碼如下：〔程序代碼〕ch3example1B.m%ch3example1B.mclear;f_p=2400;f_s=5000;R_p=3;R_s=25;%設計要求指標[n,fn]=ellipord(f_p,f_s,R_p,R_s,'s');%計算階數和截止頻率Wn=2*pi*fn;%轉換為角頻率[b,a
4、]=ellip(n,R_p,R_s,Wn,'s');%計算H(s)f=0:100:10000;%計算頻率點和頻率範圍s=j*2*pi*f;%s=jw=j*2*pi*fH_s=polyval(b,s)./polyval(a,s);%計算相應頻率點處H(s)的值figure(1);subplot(2,1,1);plot(f,20*log10(abs(H_s)));%幅頻特性axis([010000-401]);xlabel('頻率Hz');ylabel('幅度dB');subplot(2,1,2);plo
5、t(f,angle(H_s));%相頻特性xlabel('頻率Hz');ylabel('相角rad');figure(2);freqs(b,a);%也可用指令freqs直接畫出H(s)的頻率響應曲線。〔實例3.2〕試設計一個巴特沃斯型數字低通濾波器，設採樣率為8000Hz，fp=2100Hz，fs=2500Hz，Rp=3dB，Rs=25dB。設計程序代碼如下：〔程序代碼〕ch3example2A.m%ch3example2A.mf_N=8000;%採樣率f_p=2100;f_s=2500;R_p=3;
6、R_s=25;%設計要求指標Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);%計算歸一化頻率[n,Wn]=buttord(Wp,Ws,R_p,R_s);%計算階數和截止頻率[b,a]=butter(n,Wn);%計算H(z)figure(1);freqz(b,a,1000,8000)%作出H(z)的幅頻相頻圖,freqz(b,a,計算點數,採樣率)subplot(2,1,1);axis([04000-303])figure(2);%第二種作圖方法f=0:40:4000;%計算頻率點和頻率範
7、圍z=exp(j*2*pi*f./(f_N));%H_z=polyval(b,z)./polyval(a,z);%計算相應頻率點處H(s)的值subplot(2,1,1);plot(f,20*log10(abs(H_z)));%幅頻特性axis([04000-401]);xlabel('頻率Hz');ylabel('幅度dB');subplot(2,1,2);plot(f,angle(H_z));%相頻特性xlabel('頻率Hz');ylabel('相角rad');〔實例3.3〕試設計一個切比雪夫1
8、型高通數字濾波器，採樣率為8000Hz，fp=1000Hz，fs=700Hz，Rp=3dB，Rs=20dB。設計程序代碼如下：〔程序代碼〕ch3example3A.m%ch3example3A.mf_N=8000;%採樣率〔程序代碼〕ch3example3A.m%ch3example3A.mf_N=8000;%採樣率〔實例3.4〕試設計一橢圓型帶通數字濾波器。設採樣率為10000Hz，fp=[1000;1500]Hz，fs=[600;1900]Hz，Rp
