%考爾Cauer濾波器的離散濾波算法 模擬程序
%程序設(shè)計(jì):劉輝
%pangelaca@yahoo.com.cn (email/msn)
%由于我在柏林工業(yè)大學(xué)學(xué)習(xí),源程序的一些變量名由德語縮寫成.(但其實(shí)閱讀上沒有影響)
%各種參數(shù)都可以替換,產(chǎn)生不同階,不同極點(diǎn),不同零點(diǎn)的Cauer濾波器.
%程序會顯示其對數(shù)振幅特性曲線,相位特性曲線,沖激響應(yīng),階躍響應(yīng)等,涉及了信號處理的方方面面.
%將圖像顯示部分去掉,利用產(chǎn)生的分子系數(shù)和分母系數(shù),可以實(shí)現(xiàn)直接對任何輸入信號的濾波.

%各參數(shù)說名:
%N: 階數(shù)
%Pol_Cauer: N個極點(diǎn)的位置
%Null_Cauer: N個零點(diǎn)的位置
%kCauer: Cauer系數(shù)

%源程序里給出的7階的Cauer濾波器各零極點(diǎn)參數(shù),是一份比較精確的珍貴數(shù)據(jù).

N=7;

Null_Cauer=[0.903710+0.428145i, 0.903710-0.428145i, 0.965298+0.261152i, 0.965298-0.261152i, 0.975060+0.221964i, 0.975060-0.221964i, 0]
Pol_Cauer=[0.977990+0.076022i, 0.977990-0.076022i, 0.981010+0.130168i, 0.981010-0.130168i, 0.984920+0.154134i, 0.984920-0.154134i, 0.976644]
kCauer=9.11e-5;


KoefNull_Cauer=kCauer* poly (Null_Cauer);
KoefPol_Cauer=poly (Pol_Cauer);

[HzC,w]=freqz(KoefNull_Cauer,KoefPol_Bassel,2048);

AmpCauer=log(abs(HzC));
subplot(4,1,1);
plot (w,AmpCauer);
title ('對數(shù)幅頻曲線');

PhaCauer=angle(HzC);
subplot(4,1,2);
plot (PhaCauer);
title ('相位特性曲線');

ImpulsantwortCauer=filter(KoefNull_Cauer, KoefPol_Cauer,ImpulsFolge);
subplot(4,1,3);
bar (ImpulsantwortCauer);
title ('沖激響應(yīng)');


SprungantwortCauer=filter(KoefNull_Cauer,KoefPol_Cauer,SprungFolge);
subplot(4,1,4);
bar (SprungantwortCauer);
title ('階躍響應(yīng)');