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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
. p# }0 _! w) G滤波函数：
" G; e6 Z8 h% l  S" ^: X, z$ d6 j#define  NUM_TAPS   48
#define  NUM_DATA   80


$ C9 Q1 c- }) F" b+ kvoid floatPointFir(float *x, float *h, short order, float *y, float *w)
# I2 f9 Y) i0 H{
/ {& E0 {  q$ z8 I    short i;
    float sum;
6 Q/ `' F- W; R' b w[0] = *x++;                        // Get the current data to delay line
& M0 D8 p$ ~/ V6 Y: N8 L for (sum=0, i=0; i<order; i++)      // FIR filter processing
1 D; O3 g4 P  O8 |7 ^ {
  sum += h * w;
}
*y++ = sum;                         // Save filter output
& U$ b! B6 F  ] for (i=order-1; i>0; i--)           // Update data delay line
0 }6 s  [  p7 f {
  w = w[i-1] ;
}
; q  H/ O; y5 Z4 |9 U}

( w  R2 P1 h0 L. O9 i" F7 l0 _
测试代码：
float w[NUM_TAPS];
void main()
{
    FILE  *fpIn,*fpOut;
    short i;
    char  temp[2];
    float x,  // Input data
          y;  // Output data
' M) Q# A6 j: }  
    fpIn = fopen("..\\data\\input.pcm", "rb");
. H* ~$ M# ]9 ~5 x9 Z7 [0 G7 _: ]    fpOut = fopen("..\\data\\output.pcm", "wb");
& Z# H! t& K  e( X" e  d/ s' r    if (fpIn == NULL)
, R' B+ [' Q7 N+ w    {
8 g/ {& j7 y3 q; e9 }' C6 |/ }        printf("Can't open input file\n");
9 @5 r. ^3 U. ]: e9 h8 B        exit(0);
! `  s8 q$ }, O& i; B    }

    // Initialize for filtering process
* L) z& u5 m. Z+ F8 b    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
& \2 k2 x5 T2 j  b+ _. N9 K    }
    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
, Q  [5 B' d- l/ B1 w" y4 x        // 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);
  E* n$ p$ g7 V* u6 W6 b" u    }
    fclose(fpIn);
    fclose(fpOut);
}

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾