%題目：利用MATLAB結(jié)合雙線(xiàn)性變換法設(shè)計(jì)?一個(gè)數(shù)字切比雪夫帶通IIR濾波
%姓名：羅揚(yáng)靜
%班別：通信172班
%學(xué)號(hào)：201708402208
clear
%數(shù)字濾波器的技術(shù)指標(biāo)?即任務(wù)設(shè)計(jì)要求
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）;

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

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

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

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

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

%數(shù)字帶通濾波器幅頻響應(yīng)曲線(xiàn)
subplot（323）;plot（w/piabs（h））;
xlabel（‘w（rad）‘）;ylabel（‘|H（z）|‘）;
title（‘?dāng)?shù)字帶通濾波器幅頻特性曲線(xiàn)‘）;

%數(shù)字帶通濾波器幅頻響應(yīng)曲線(xiàn)（dB）
subplot（324）;plot（w/pi20*log10（abs（h）））;axis（[-12-25050]）;
xlabel（‘w（rad）‘）;ylabel（‘|H（z）|.dB‘）;
title（‘?dāng)?shù)字帶通濾波器幅頻特性曲線(xiàn)‘）;

%數(shù)字帶通濾波器相頻特性曲線(xiàn)（dB）
subplot（325）;plot（w/piphz）;grid?on;
xlabel（‘w（rad）‘）;
ylabel（‘H（z）‘）;
title（‘?dāng)?shù)字帶通濾波器相頻特性曲線(xiàn)‘）;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????1689??2019-11-15?17:10??DSP課設(shè)\IIR.m
?????文件??????460800??2019-11-15?20:46??DSP課設(shè)\MATLAB結(jié)合雙線(xiàn)性變換法設(shè)計(jì)一個(gè)數(shù)字切比雪夫帶同IIR濾波器.doc