DSP的IIR数字滤波器实现程序源码

亲测有效的IIR数字滤波器C程序实现，输入与输出有明显的滤波效果，如图



附录（C程序源码）：

1. #include "DSP2833x_Device.h"     // DSP2833x Headerfile Include File
   
2. #include "DSP2833x_Examples.h"   // DSP2833x Examples Include File
   
3. #include  <math.h>
   
4. #define IIRNUMBER 2
   
5. #define SIGNAL1F 1000
   
6. #define SIGNAL2F 4500
   
7. #define SAMPLEF  10000
   
8. #define PI       3.1415926
   
9. 
   
10. float InputWave();
    
11. float IIR();
    
12. 
    
13. float fBn[IIRNUMBER]={ 0.0,0.7757 };
    
14. float fAn[IIRNUMBER]={ 0.1122,0.1122 };
    
15. float fXn[IIRNUMBER]={ 0.0 };
    
16. float fYn[IIRNUMBER]={ 0.0 };
    
17. float fInput,fOutput;
    
18. float fSignal1,fSignal2;
    
19. float fStepSignal1,fStepSignal2;
    
20. float f2PI;
    
21. int i;
    
22. float fin[256],fout[256];
    
23. //输入fIn毛刺正弦波，输出fOut较光滑正弦波,Graph: 256,32-bit floating integer,
    
24. int nIn,nOut;
    
25. 
    
26. main()
    
27. {
    
28.     nIn=0; nOut=0;
    
29.     fInput=fOutput=0;
    
30.     f2PI=2*PI;
    
31.     fSignal1=0.0;
    
32.     fSignal2=PI*0.1;
    
33. //  fStepSignal1=2*PI/30;
    
34. //  fStepSignal2=2*PI*1.4;
    
35.     fStepSignal1=2*PI/50;
    
36.     fStepSignal2=2*PI/2.5;
    
37.     while ( 1 )
    
38.     {
    
39.         fInput=InputWave();
    
40.         fin[nIn]=fInput;
    
41.         nIn++; nIn%=256;
    
42.         fOutput=IIR();
    
43.         fout[nOut]=fOutput;
    
44.         nOut++;             // break point
    
45.         if ( nOut>=256 )
    
46.         {
    
47.             nOut=0;
    
48.         }
    
49.     }
    
50. }
    
51. 
    
52. float InputWave()
    
53. {
    
54.     for ( i=IIRNUMBER-1;i>0;i-- )
    
55.     {
    
56.         fXn[i]=fXn[i-1];
    
57.         fYn[i]=fYn[i-1];
    
58.     }
    
59.     fXn[0]=sin((double)fSignal1)+cos((double)fSignal2)/6.0;
    
60.     fYn[0]=0.0;
    
61.     fSignal1+=fStepSignal1;
    
62.     if ( fSignal1>=f2PI )   fSignal1-=f2PI;
    
63.     fSignal2+=fStepSignal2;
    
64.     if ( fSignal2>=f2PI )   fSignal2-=f2PI;
    
65.     return(fXn[0]);
    
66. }
    
67. 
    
68. float IIR()
    
69. {
    
70.     float fSum;
    
71.     fSum=0.0;
    
72.     for ( i=0;i<IIRNUMBER;i++ )
    
73.     {
    
74.         fSum+=(fXn[i]*fAn[i]);
    
75.         fSum+=(fYn[i]*fBn[i]);
    
76.     }
    
77.     return(fSum);
    
78. }

复制代码

全部资料51hei下载地址：
IIR滤波C程序.docx (13.53 KB, 下载次数: 78)

2018-11-2 13:35 上传

点击文件名下载附件
下载积分: 黑币 -5

值得参考的例程