exp6b.c

(Ebook-Pdf) Dsp - Real Time Digital Signal Processing (Usando Tms320-55Xx). 有书

/* 
   exp6b.c - Direct form II - IIR function implementation 
             in Fixed-point C and using signal generator
*/                           

#define M       128    /* Number of samples per block */
#define Ns      2      /* Number od 2nd order sections */

/* Low-pass IIR filter coefficients in Q14 format */
int C[Ns*5]={ /* i=section number */
              /* A[i][1],A[i][2],B[i][2],B[i][0],B[i][1] */
                  -14187,   3505,   1624,   1624,   3249,
                  -18430,   9454,   1624,   1624,   3249};
/* IIR filter delay line:
w[]=w[i][n-1],w[i+1][n-1],...,w[i][n-2],w[i+1][n-2],... */
int w[Ns*2];	 

int out[M];  
int in[M];  

/* IIR filter function */
extern void iir(int*, int, int *, int *, int, int *);

/* Software signal generator */
extern void signal_gen2(int *, unsigned int);

void main(void)
{
    int i;

    /* Initialize IIR filter delay line and signal generator */
    for (i=0; i<Ns*2; i++)
        w[i]=0; 

    /* IIR filter experiment start */  
    for (;;)           
    {
        signal_gen2(in,M);     /* Generate a buffer of samples */
        iir(in,M,out,C,Ns,w);  /* Filter a buffer of samples */                         
    }
}