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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
- O' \0 e: ]% c滤波函数：
#define  NUM_TAPS   48
% j4 f/ U* C  p: c2 E1 S7 |7 e#define  NUM_DATA   80
. [3 f4 n8 m$ O/ ~3 f0 @! W6 D

. @$ p) q2 B) M, E. {% K4 Mvoid floatPointFir(float *x, float *h, short order, float *y, float *w)
- Z9 c' P* x% b{
    short i;
" ?" j0 {5 l0 P0 D; g    float sum;
. M! \. w& ~5 `- i+ u w[0] = *x++;                        // Get the current data to delay line
; ~/ s4 g3 I# ^8 N% o0 I for (sum=0, i=0; i<order; i++)      // FIR filter processing
0 ~# F. D* y: b* p6 ^ {
  sum += h * w;
}
*y++ = sum;                         // Save filter output
for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
}
}


测试代码：
float w[NUM_TAPS];
9 B2 e" C2 ~5 ~6 k6 G  p9 ?void main()
  }1 z% T* g  ~1 n- |{
4 }  m* d# U. O6 c    FILE  *fpIn,*fpOut;
+ Q+ d7 x9 G) A: [8 v$ x    short i;
2 ]6 y2 v; V) h- s+ r" ^8 x    char  temp[2];
    float x,  // Input data
          y;  // Output data
+ q7 g7 E" o! h$ R7 F  
    fpIn = fopen("..\\data\\input.pcm", "rb");
6 I% d, t) Y* `    fpOut = fopen("..\\data\\output.pcm", "wb");
" S1 X/ {& M0 E+ B& s2 K& d    if (fpIn == NULL)
- x& B+ ]5 t3 w2 D9 k& v' Q* j) t6 f    {
        printf("Can't open input file\n");
/ j, `9 R- s; S# y, J4 ]. z        exit(0);
9 N3 G3 {6 u7 i  f; Q. N& }. N    }

    // Initialize for filtering process
8 p4 e; g3 U& K, M* [. I$ Q" l    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
    }
2 a7 l* _: H# @# L5 a+ n    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
    {
; H. S4 }+ y* @! L8 p; _        x = (float)((temp[1]<<8)|(temp[0]&0xff));
, [/ f7 ^4 E0 l6 j* d4 r        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
$ E' P! b* G3 r) {/ Y: g6 y* a        temp[0] = (char)((short)y&0xff);
4 _) J# N" A- Y1 f+ O        temp[1] = (char)(((short)y>>8)&0xff);
! n! z  W; I6 q2 B8 j3 I. ?        fwrite(&temp, sizeof(char), 2, fpOut);
3 u) ^* R% q5 J" l7 m" a1 z% |    }
$ t% p. O" i* T& E0 B! U6 r    fclose(fpIn);
    fclose(fpOut);
' f" Z+ U* u  A8 t}

) W: m2 K( U+ z' g2 q& q
- T0 O  o4 N8 }( S( `/ e

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾