DSP的IIR数字滤波器实现

sdsdwwwe22

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

, I% F/ L% X5 u  D5 g
1 E$ L* ~* N9 y: T. K
附录（C程序源码）：

• #include "DSP2833x_Device.h"     // DSP2833x HeadeRFile Include File
• #include "DSP2833x_Examples.h"   // DSP2833x Examples Include File
• #include  <math.h>
• #define IIRNUMBER 2
• #define SIGNAL1F 1000
• #define SIGNAL2F 4500
• #define SAMPLEF  10000
• #define PI       3.1415926
• float InputWave();
• float IIR();
• float fBn[IIRNUMBER]={ 0.0,0.7757 };
• float fAn[IIRNUMBER]={ 0.1122,0.1122 };
• float fXn[IIRNUMBER]={ 0.0 };
• float fYn[IIRNUMBER]={ 0.0 };
• float fInput,fOutput;
• float fSignal1,fSignal2;
• float fStepSignal1,fStepSignal2;
• float f2PI;
• int i;
• float fin[256],fout[256];
• //输入fIn毛刺正弦波，输出fOut较光滑正弦波,Graph: 256,32-bit floating integer,
• int nIn,nOut;
• main()
• {
•     nIn=0; nOut=0;
•     fInput=fOutput=0;
•     f2PI=2*PI;
•     fSignal1=0.0;
•     fSignal2=PI*0.1;
• //  fStepSignal1=2*PI/30;
• //  fStepSignal2=2*PI*1.4;
•     fStepSignal1=2*PI/50;
•     fStepSignal2=2*PI/2.5;
•     while ( 1 )
•     {
•         fInput=InputWave();
•         fin[nIn]=fInput;
•         nIn++; nIn%=256;
•         fOutput=IIR();
•         fout[nOut]=fOutput;
•         nOut++;             // break point
•         if ( nOut>=256 )
•         {
•             nOut=0;
•         }
•     }
• }
• float InputWave()
• {
•     for ( i=IIRNUMBER-1;i>0;i-- )
•     {
•         fXn=fXn[i-1];
•         fYn=fYn[i-1];
•     }
•     fXn[0]=sin((double)fSignal1)+cos((double)fSignal2)/6.0;
•     fYn[0]=0.0;
•     fSignal1+=fStepSignal1;
•     if ( fSignal1>=f2PI )   fSignal1-=f2PI;
•     fSignal2+=fStepSignal2;
•     if ( fSignal2>=f2PI )   fSignal2-=f2PI;
•     return(fXn[0]);
• }
• float IIR()
• {
•     float fSum;
•     fSum=0.0;
•     for ( i=0;i<IIRNUMBER;i++ )
•     {
•         fSum+=(fXn*fAn);
•         fSum+=(fYn*fBn);
•     }
•     return(fSum);
  5 T6 a& [( a; C6 E

% l* k& ~4 |( |+ f. A

  T% K1 ]9 G) V; ~. [! u# k2 }6 M

道法自然

这波形不错呀

kiygb

   感谢分享

hunterccc

毛刺怎么解决

Heaven_1

这个滤波器比较常见，也很重要