%題目：利用MATLAB結合雙線性變換法設計?一個數字切比雪夫帶通IIR濾波
%姓名：羅揚靜
%班別：通信172班
%學號：201708402208
clear
%數字濾波器的技術指標?即任務設計要求
wp1=0.45*pi;???wp2=0.65*pi;????wp0=0.55*pi;???
ws1=0.3*pi;????ws2=0.75*pi;
Ap=1;??????????As=40;?T=2;

%帶通到低通的頻率變換
Wp1=（2/T）*tan（wp1/2）;???Wp2=（2/T）*tan（wp2/2）;
Ws1=（2/T）*tan（ws1/2）;???Ws2=（2/T）*tan（ws2/2）;?Wp=（2/T）*tan（wp0/2）;
BW=Wp2-Wp1;?%帶通濾波器的通帶寬度??????
W0=Wp1*Wp2;?
WP=1;?%歸一化處理????
WS=WP*（W0^2-Ws1^2）/（Ws1*BW）;

%切比雪夫模擬低通原型濾波器設計
[NWn]=cheb1ord（WPWSApAs‘s‘）;
[B1A1]=cheby1（NApWn‘s‘）;

%模擬低通原型濾波器幅頻特性曲線（dB）
[h1w1]=freqs（B1A1）;
subplot（321）;plot（w1/pi20*log（abs（h1）））;grid?on;
xlabel（‘w（rad）‘）;
ylabel（‘|H（jw）|.dB‘）;
title（‘模擬低通濾波器幅頻特性曲線‘）;

%?由模擬低通原型濾波器變換為模擬帶通濾波器
[B2A2]=lp2bp（B1A1WpBW）;
[h2w2]=?freqs（B2A2）;

%模擬帶通濾波器幅頻特性曲線（dB）
subplot（322）;plot（w220*log（abs（h2）））;axis（[03-40050]）;grid?on;
xlabel（‘w（rad）‘）;
ylabel（‘|H（jw）|.dB‘）;
title（‘模擬帶通濾波器幅頻特性曲線‘）;

%雙線性變換：由模擬濾波器向數字濾波器的變換
[B3A3]=bilinear（B2A20.5）;
[hw]=freqz（B3A364）;
phz=unwrap（angle（h））;

%數字帶通濾波器幅頻響應曲線
subplot（323）;plot（w/piabs（h））;
xlabel（‘w（rad）‘）;ylabel（‘|H（z）|‘）;
title（‘數字帶通濾波器幅頻特性曲線‘）;

%數字帶通濾波器幅頻響應曲線（dB）
subplot（324）;plot（w/pi20*log10（abs（h）））;axis（[-12-25050]）;
xlabel（‘w（rad）‘）;ylabel（‘|H（z）|.dB‘）;
title（‘數字帶通濾波器幅頻特性曲線‘）;

%數字帶通濾波器相頻特性曲線（dB）
subplot（325）;plot（w/piphz）;grid?on;
xlabel（‘w（rad）‘）;
ylabel（‘H（z）‘）;
title（‘數字帶通濾波器相頻特性曲線‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1689??2019-11-15?17:10??DSP課設\IIR.m
?????文件??????460800??2019-11-15?20:46??DSP課設\MATLAB結合雙線性變換法設計一個數字切比雪夫帶同IIR濾波器.doc