fir滤波器设计,c语言,.docVIP

C语言编写FIR滤波器，代码如下： #include stdio.h #include math.h #includestdlib.h /**********************************宏定义**************************************/ #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; iSize_Datas; i++) // FIR运算 { for(j=0;jORDER;j++) { if(i=j) sum+=h[j]*Xin[i-j]; // 时域卷积 else ; } output[i] = sum ; } // 考虑到前33个点为不完全累加和，故抛去前33个点。应该是前后各32个点//都是不完全累加和， //都应该去掉，对于数据分段进入滤波器的情况，应该把前一段的后面数据放到//下一段的前面 for(i=0;iSize_Datas-ORDER;i++) { output[i]=output[i+ORDER]; } }