fir滤波器设计c语言

C语言编写FIR滤波器，代码如下：

#include

#include

#include

/**********************************宏定义**************************************/

#define PI 3.1415926535

#define Size_Datas 512

/*******************************数据类型构造***********************************/

int i=0,j=0, nNow=0;

// h[]参数使用matlab计算获得,用fopen，接着fprintf,最后fclose等命令保存为dat文件，不要直接用save命令。

float h[]={

#include "b_fir.dat"

};

const int ORDER = sizeof(h)/sizeof(h[0]);

//输入数据形式与滤波器参数形式相同

float Xin[] = {

#include "inputwave0_fir.dat"

};

float sum=0;

float output[Size_Datas],w[Size_Datas];

int n1;

void main(void)

{

for (i=0; i

{

for(j=0;j

{

if(i>=j)

sum+=h[j]*Xin[i-j]; // 时域卷积

else

;

}

output[i] = sum ;

}

// 考虑到前33个点为不完全累加和，故抛去前33个点。应该是前后各32个点//都是不完全累加和，

//都应该去掉，对于数据分段进入滤波器的情况，应该把前一段的后面数据放到//下一段的前面

for(i=0;i

{

output[i]=output[i+ORDER];

}

}