C語(yǔ)言代碼 - 巴特沃斯濾波器c語(yǔ)言代碼

計(jì)算出穩(wěn)定的極點(diǎn)之后，就可以進(jìn)行傳遞函數(shù)的計(jì)算了。傳遞的函數(shù)的計(jì)算，就像下式一樣

這里，為了得到模擬濾波器的系數(shù)，需要將分母乘開(kāi)。很顯然，這里的極點(diǎn)不一定是整數(shù)，或者來(lái)說(shuō)，這里的乘開(kāi)需要做復(fù)數(shù)運(yùn)算。其復(fù)數(shù)的乘法代碼如下，

int Complex_Multiple（COMPLEX a，COMPLEX b，

double *Res_Real，double *Res_Imag）

{

*（Res_Real） = （a.Real_part）*（b.Real_part） - （a.Imag_Part）*（b.Imag_Pa

rt）;

*（Res_Imag）= （a.Imag_Part）*（b.Real_part） + （a.Real_part）*（b.Imag_Par

t）;

return （int）1;

}

有了乘法代碼之后，我們現(xiàn)在簡(jiǎn)單的情況下，看看其如何計(jì)算其濾波器系數(shù)。我們做如下假設(shè)

這個(gè)時(shí)候，其傳遞函數(shù)為

將其乘開(kāi)，其大致的關(guān)系就像下圖所示一樣。

Res［0］.Real_part = poles［0］.Real_part;

Res［0］.Imag_Part= poles［0］.Imag_Part;

Res［1］.Real_part = 1;

Res［1］.Imag_Part= 0; 5.

for（count_1 = 0;count_1 < N-1;count_1++）

{

for（count = 0;count <= count_1 + 2;count++）

{

if（0 == count）

{

Complex_Multiple（Res［count］， poles［count_1+1］，

&（Res_Save［count］.Real_part），

&（Res_Save［count］.Imag_Part））;

}

else if（（count_1 + 2） == count）

{

Res_Save［count］.Real_part += Res［count - 1］.Real_part;

Res_Save［count］.Imag_Part += Res［count - 1］.Imag_Part;

}

else

{

Complex_Multiple（Res［count］， poles［count_1+1］，

&（Res_Save［count］.Real_part），

&（Res_Save［count］.Imag_Part））;

1 Res_Save［count］.Real_part += Res［count - 1］.Real_part;

Res_Save［count］.Imag_Part += Res［count - 1］.Imag_Part;

}

}

*（b+N） = *（a+N）;