标题: DSP的IIR数字滤波器实现程序源码 [打印本页]
作者: Mr_SYC    时间: 2018-11-2 13:35
标题: 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

  10. float InputWave();
  11. float IIR();

  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;

  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. }

  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. }

  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, 下载次数: 84)

作者: JNypf    时间: 2024-5-5 14:45
值得参考的例程



欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1