實際可運行的FIR濾波器系數設計的C語言程序

#include"math.h"

#include"stdio.h"

/*n:濾波器的階數,n是奇數時可用來設計各種濾波器

band:濾波器類型1、2、3、4分別對應低通，高通，帶通，帶阻

fln:通帶下邊界頻率

fhn:通帶上邊界頻率

wn:窗函數的類型1-7分別對應

h:長度為n+1,存放濾波器的系數

fs:采樣頻率*/

void firwin(n,band,fln,fhn,wn,h)

int n,band,wn;

double fln,fhn,h[];

{

int i,n2,mid;

double s,pi,wc1,wc2,beta,delay,fs;

double window();//窗函數的計算

beta=0.0;

if(wn==7)//只有凱塞窗需要輸入系數beta

{printf("input beta parameter of Kaiser window(2

scanf("%1f",&beta);

}

pi=4.0*atan(1.0);//pi=PI;

if((n%2)==0)/*如果階數n是偶數*/

{n2=n/2+1;/**/

mid=1;//

}

else

{n2=n/2;//n是奇數,則窗口長度為偶數

mid=0;//

}

delay=n/2.0;

wc1=2.0*pi*fln;//

if(band>=3) /*先判斷用戶輸入的數據，如果band參數大于3*/

{wc2=2.0*pi*fhn;}

switch(band)

{

case 1:

{for(i=0;i<=n2;i++)

{s=i-delay;//

h[i]=(sin(wc1*s/fs)/(pi*s))*window(wn,n+1,i,beta);//低通,窗口長度=階數+1，故為n+1

h[n-i]=h[i];

}

if(mid==1)h[n/2]=wc1/pi;//n為偶數時，修正中間值系數

break;

}

case 2:

{for(i=0;i<=n2;i++)

{s=i-delay;

h[i]=(sin(pi*s)-sin(wc1*s/fs))/(pi*s);//高通-//對

h[i]=h[i]*window(wn,n+1,i,beta);

h[n-i]=h[i];

}

if(mid==1)h[n/2]=1.0-wc1/pi;//對

break;

}

case 3:

{for(i=0;i

{s=i-delay;

h[i]=(sin(wc2*s/fs)-sin(wc1*s/fs))/(pi*s);//帶通-//對

h[i]=h[i]*window(wn,n+1,i,beta);

h[n-i]=h[i];

}

if(mid==1)h[n/2]=(wc2-wc1)/pi;//對

break;

}

case 4:

{for(i=0;i<=n2;i++)

{s=i-delay;

h[i]=(sin(wc1*s/fs)+sin(pi*s)-sin(wc2*s/fs))/(pi*s);//帶阻-//對

h[i]=h[i]*window(wn,n+1,i,beta);

h[n-i]=h[i];

}

if(mid==1)h[n/2]=(wc1+pi-wc2)/pi;

break;

}

}

}

/*返回窗函數的表達式

n:窗口長度

type:選擇窗函數的類型

beta:生成凱塞窗的系數*/

static double window(type,n,i,beta)

int i,n,type;

double beta;

{

int k;

double pi,w;

double kaiser();

pi=4.0*atan(1.0);//pi=PI;

w=1.0;

switch(type)

{

case 1:

{w=1.0;//矩形窗

break;

}

case 2:

{k=(n-2)/10;

if(i<=k)

w=0.5*(1.0-cos(i*pi/(k+1)));//圖基窗

break;

}

case 3:

{w=1.0-fabs(1.0-2*i/(n-1.0));//三角窗

break;

}

case 4:

{w=0.5*(1.0-cos(2*i*pi/(n-1)));//漢寧窗

break;

}

case 5:

{w=0.54-0.46*cos(2*i*pi/(n-1));//海明窗

break;

}

case 6:

{w=0.42-0.5*cos(2*i*pi/(n-1))+0.08*cos(4*i*pi/(n-1));//布萊克曼窗

break;

}

case 7:

{w=kaiser(i,n,beta);//凱塞窗

break;

}

}

return(w);

}

static double kaiser(i,n,beta)//因凱塞窗的表達式復雜，調用貝塞爾窗表達式

int i,n;

double beta;

{

double a,w,a2,b1,b2,beta1;

double bessel0();

b1=bessel0(beta);

a

總結

以上是生活随笔為你收集整理的fir滤波器算法c语言程序,FIR滤波器设计C语言程序的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。