c语言描述fir.docVIP

FIR滤波器c语言源程序设计 #include stdio.h #ifdef WIN32 #include conio.h #endif #define SAMPLE double /* define the type used for data samples */ void clear(int ntaps, SAMPLE z[]) { int ii; for (ii = 0; ii ntaps; ii++) { z[ii] = 0; } } SAMPLE fir_basic(SAMPLE input, int ntaps, const SAMPLE h[], SAMPLE z[]) { int ii; SAMPLE accum; /* store input at the beginning of the delay line */ z[0] = input; /* calc FIR */ accum = 0; for (ii = 0; ii ntaps; ii++) { accum += h[ii] * z[ii]; } /* shift delay line */ for (ii = ntaps - 2; ii = 0; ii--) { z[ii + 1] = z[ii]; } return accum; } SAMPLE fir_circular(SAMPLE input, int ntaps, const SAMPLE h[], SAMPLE z[], int *p_state) { int ii, state; SAMPLE accum; state = *p_state; /* copy the filters state to a local */ /* store input at the beginning of the delay line */ z[state] = input; if (++state = ntaps) { /* incr state and check for wrap */ state = 0; } /* calc FIR and shift data */ accum = 0; for (ii = ntaps - 1; ii = 0; ii--) { accum += h[ii] * z[state]; if (++state = ntaps) { /* incr state and check for wrap */ state = 0; } } *p_state = state; /* return new state to caller */ return accum; } SAMPLE fir_shuffle(SAMPLE input, int ntaps, const SAMPLE h[], SAMPLE z[]) { int ii; SAMPLE accum; /* store input at the beginning of the delay line */ z[0] = input; /* calc FIR and shift data */ accum = h[ntaps - 1] * z[ntaps - 1]; for (ii = ntaps - 2; ii = 0; ii--) { accum += h[ii] * z[ii]; z[ii + 1] = z[ii]; } return accum; } SAMPLE fir_split(SAMPLE input, int ntaps, const SAMPLE h[], SAMPLE z[], int *p_state) { int ii, end_ntaps, state = *p