之前学习的内容，因为没有开发板，所以只能ccs联合MATLAB一起做数字滤波器的设计。这个对于刚刚接触的新手还是非常有难度的，像我就熬了一个通宵才做出来，当然还是有点菜。详细的程序都在安装包里。可以参考。而且ccs这个软件有点迷，有时候没有图或者一些错误，可能重启一下软件就好了。
源程序如下:
- /*
- * main.c
- */
- #include "math.h"
- #include "stdio.h"
- #include <stdint.h>
- #include "stdlib.h"
- #include "string.h"
- #include "limits.h"
- //#include"fdacoefsiir900.h"
- //#include "fdacoefsiir600.h"
- //#include "fdacoefsiir50+900.h"
- //#include "fdacoefsiir50.h"
- #include "fdacoefs50001.h"
- #define length 2048
- #define pi 3.1415926
- long fs=10000;
- int f1=50;
- int f2=200;
- int f3=600;
- int f4=900;
- #define w1 2*pi*f1/fs
- #define w2 2*pi*f2/fs
- #define w3 2*pi*f3/fs
- #define w4 2*pi*f4/fs
- double input[length];
- double output[length];
- static double xlast[2];
- static double mlast[2];
- static double IIR_DR2(double x,double *plast,const double (*A)[3],const double (*B)[3])
- {
- double tmp,last;
- tmp = x*B[0][0];
- last = tmp - (A[1][1]*plast[0] + A[1][2]*plast[1]);
- tmp = last + (B[1][1] * plast[0] + B[1][2]*plast[1]);
- plast[1] = plast[0];
- plast[0] = last;
- return tmp;
- }
- double IIR_Filter(double x)
- {
- double mid,y;
- mid = IIR_DR2(x,xlast,DEN,NUM);
- y = IIR_DR2(mid,mlast,&DEN[2],&NUM[2]);//二阶滤波器组合成更高阶数的滤波器
- //更多阶数...
- return y;
- }
- void Init_Filter(void)//初始化中间数值
- {
- xlast[0] = 0;
- xlast[1] = 0;
- mlast[0] = 0;
- mlast[1] = 0;
- }
- int main(void)
- {
- unsigned int i,n;
- Init_Filter();
- for(i=0;i<length;i++)
- {
- input[i]=2048*sin(w1*i)+2048*sin(w2*i)+2048*sin(w3*i)+2048*sin(w4*i);
- }
- for(n=0;n<length;n++)
- {
- output[n]=IIR_Filter(input[n]);
- }
- for(;;);
- }
所有资料51hei提供下载:
7.zip
(101.16 KB, 下载次数: 48)
|
