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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
+ E" r- G( l# b+ r! B8 ~# r滤波函数：
; B* [6 b3 ~: B" g9 v#define  NUM_TAPS   48
#define  NUM_DATA   80

( ]2 ~. S8 V  w2 V9 z  n# b
void floatPointFir(float *x, float *h, short order, float *y, float *w)
{
    short i;
    float sum;
  ]+ o0 _  h& b5 H w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
( u) ~" {" b& @. }5 r  i& P) O  sum += h * w;
5 J, Y, B- [: X* K  f; I; u }
/ T% w8 u  H, X *y++ = sum;                         // Save filter output
$ C. h. T+ @- \0 k# L# Q for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
. ?" H2 a: L- g }
9 A4 j3 |" F# ^# J1 `}
4 J" I$ H3 h$ C7 {2 a4 Y

* @$ u5 S/ F* i* W测试代码：
float w[NUM_TAPS];
7 d2 u$ Z! P( d: [0 Lvoid main()
+ o2 r" S. t9 ^{
    FILE  *fpIn,*fpOut;
& U6 ?! a) o! R. J2 ]% w/ I2 w    short i;
    char  temp[2];
    float x,  // Input data
          y;  // Output data
3 a* D5 ]) s  D( `; Y+ G0 j. S' A2 {  
/ Q! d* T  T! `. {" c0 i+ B: l    fpIn = fopen("..\\data\\input.pcm", "rb");
5 q4 y! o/ y/ k  t0 S9 T    fpOut = fopen("..\\data\\output.pcm", "wb");
2 q3 [( F2 y/ u    if (fpIn == NULL)
    {
7 m- d% A& D$ J2 H; a. ^        printf("Can't open input file\n");
        exit(0);
    }

    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
' a0 j! ?3 s3 |6 j    {
        w = 0.0;
% D8 d1 G: ]6 _7 ?1 }    }
    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
6 f* i6 ?4 T5 k( v    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
) G; W" R" Z8 m" i& W# }/ [        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
3 R7 Z' W* G$ b  e" h. n        temp[0] = (char)((short)y&0xff);
0 m* P$ U$ P* ]( L& m8 L. x        temp[1] = (char)(((short)y>>8)&0xff);
8 D, _' ^6 w9 Q) |. I        fwrite(&temp, sizeof(char), 2, fpOut);
    }
    fclose(fpIn);
    fclose(fpOut);
# ~- n: T7 A/ M4 c}

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾