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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
- T( q* d9 `9 n( ^滤波函数：
#define  NUM_TAPS   48
#define  NUM_DATA   80


" q8 w/ g# `) z" Dvoid floatPointFir(float *x, float *h, short order, float *y, float *w)
{
3 ?1 F, U* ~1 e' J* j7 w5 f6 V2 O    short i;
' y9 x8 v+ h+ C2 Z. y    float sum;
" J( T8 p& e- T5 U8 X1 o. x w[0] = *x++;                        // Get the current data to delay line
1 V: n3 k, T8 ^ for (sum=0, i=0; i<order; i++)      // FIR filter processing
: j7 S; y8 ]3 G# W  b' j( u/ M0 ^ {
- |& q, i. M' Q" @  sum += h * w;
}
( |$ \+ d! t" h5 _& e9 F *y++ = sum;                         // Save filter output
for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
}
}


) a2 Z, q) ~% c* ~% Z4 \: f测试代码：
4 [4 z( x; f7 O7 J! j# h1 qfloat w[NUM_TAPS];
void main()
{
& C# g/ {6 T2 S! Z  B# `    FILE  *fpIn,*fpOut;
, i/ [& ]" @7 e    short i;
8 S* e5 P/ e& b. G  g    char  temp[2];
- V3 J) }3 _! ~) I; j& p    float x,  // Input data
6 G5 A$ p  ?' t/ ?          y;  // Output data
  
    fpIn = fopen("..\\data\\input.pcm", "rb");
7 f: k7 n8 O. \    fpOut = fopen("..\\data\\output.pcm", "wb");
: f3 p2 D/ e* t3 M5 F7 n* P    if (fpIn == NULL)
8 E6 ~* Q3 Q3 h( f0 b6 f    {
4 ^  U" O! q, i, U# d        printf("Can't open input file\n");
' G9 a3 ]' q9 h        exit(0);
    }

+ Z& c; {3 }0 }& U3 _6 R7 u    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
    }
    // Begin filtering the data
3 O0 @& b- \8 C  l# {    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
        // Filter the data x and save output y
5 _0 v: N2 L4 D/ r6 j  `( m; M        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
4 ^( ]7 X! \$ Y  R- i1 e6 |        temp[0] = (char)((short)y&0xff);
4 ~& s4 j* l: O, M/ {/ _* {        temp[1] = (char)(((short)y>>8)&0xff);
        fwrite(&temp, sizeof(char), 2, fpOut);
    }
    fclose(fpIn);
6 l0 P  U. ]' j6 `* L8 l    fclose(fpOut);
}

  C2 l2 @0 `' p3 V

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾