实时fir滤波器的输出长度是否和输入长度一样长？

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
滤波函数：
#define  NUM_TAPS   48
& h! g# u6 S% U( o) E' O5 ~#define  NUM_DATA   80


6 H8 h7 X& _+ L. vvoid floatPointFir(float *x, float *h, short order, float *y, float *w)
{
    short i;
    float sum;
, t  Z( ]( c- I2 \1 ~ w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
  sum += h * w;
}
' Z. C5 d% C$ e- @ *y++ = sum;                         // Save filter output
7 Z+ p7 X  H- R for (i=order-1; i>0; i--)           // Update data delay line
+ ]* i, e7 D1 s, e, ~" y9 S {
  w = w[i-1] ;
& a& K0 t/ o9 T6 N$ B9 f- c- X+ U }
}


$ i' h; d: k/ I7 g) E测试代码：
float w[NUM_TAPS];
void main()
{
    FILE  *fpIn,*fpOut;
! U$ M3 H/ q' p( s    short i;
) J% o  Z& I  a& G5 K( m  C    char  temp[2];
    float x,  // Input data
- H3 A" C' C& }- D          y;  // Output data
  
    fpIn = fopen("..\\data\\input.pcm", "rb");
& y/ g' b3 v6 o% l    fpOut = fopen("..\\data\\output.pcm", "wb");
    if (fpIn == NULL)
    {
        printf("Can't open input file\n");
- g4 U& N- w' O2 m$ W! Z3 J        exit(0);
2 \2 p* I& N6 q# v" q$ t    }

/ K! x- J, x. }    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
$ y$ i8 D1 r) s4 N4 T    {
        w = 0.0;
    }
    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
  W. t1 n4 @: G% o( R    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
; A2 Y5 S5 e4 \6 u: ?6 S  i, D        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
        temp[0] = (char)((short)y&0xff);
        temp[1] = (char)(((short)y>>8)&0xff);
        fwrite(&temp, sizeof(char), 2, fpOut);
    }
6 l* j2 C* Q' c/ {1 W3 K    fclose(fpIn);
7 B9 T  N  A7 c& c+ X    fclose(fpOut);
" i, x) U% d8 L2 z}

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾