%Samp6_7
wp=0.2*pi;ws=0.3*pi;Rp=1;Rs=15; %數字濾波器截止頻率通帶波紋和阻帶衰減
Fs=100;Ts=1/Fs;Nn=128;  %采樣頻率
Wp=2/Ts*tan(wp/2.);Ws=2/Ts*tan(ws/2.);   %按頻率轉換公式進行轉換
[N,Wn]=ellipord(Wp,Ws,Rp,Rs,'s');  %計算模擬濾波器的最小階數
[z,p,k]=ellipap(N,Rp,Rs);  %設計模擬原型濾波器
[Bap,Aap]=zp2tf(z,p,k);  %零點極點增益形式轉換為傳遞函數形式
[b,a]=lp2lp(Bap,Aap,Wn); %低通轉換為低通濾波器的頻率轉換
[bz,az]=bilinear(b,a,Fs);  %運用雙線性變換法得到數字濾波器傳遞函數
[H,f]=freqz(bz,az,Nn,Fs);   %求出頻率特性
subplot(2,1,1),plot(f,20*log10(abs(H)));
xlabel('頻率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('頻率/Hz');ylabel('相位/^o');grid on;