%Bilinear函數：模擬濾波器轉換成為數字濾波器的雙線性變換法
Wp=2*pi*1400;Wp1=2*pi*Wp;
Ws=2*pi*1000;Ws1=2*pi*Ws;
Rp=0.3;
Rs=15;
Fs=20000;
%選擇濾波器的最小階數
[N,Wn]=cheb2ord(Wp1,Ws1,Rp,Rs,'s');
%創建低通巴特沃斯濾波器
[Z,P,K]=cheb2ap(N,Rs);
[A,B,C,D]=zp2ss(Z,P,K);
%低通到高通的轉變
[At1,Bt1,Ct1,Dt1]=lp2hp(A,B,C,D,Wn);
%雙線性變換法把模擬濾波器轉換為數字濾波器
[At2,Bt2,Ct2,Dt2]=bilinear(At1,Bt1,Ct1,Dt1,Fs);
%繪制幅頻響應曲線
[num,den]=ss2tf(At2,Bt2,Ct2,Dt2);
[H,W]=freqz(num,den);
plot(W*Fs/2/pi,abs(H));grid;
xlabel('頻率/Hz'),ylabel('幅度');