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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
滤波函数：
* B9 Z, f4 b' ~' \% S9 i& ~#define  NUM_TAPS   48
#define  NUM_DATA   80
/ j* a+ `4 Z8 a# y# \

void floatPointFir(float *x, float *h, short order, float *y, float *w)
{
    short i;
8 p# |- X) |) V: E    float sum;
5 b5 w3 ]# A) l% q1 `% M5 @) c4 ? w[0] = *x++;                        // Get the current data to delay line
& y( \: l& l* n for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
0 Z2 a. b1 q& Y9 J  sum += h * w;
( a* b0 v" Q) M& R. Q. t1 m! } }
/ N8 }+ l$ o: ^ *y++ = sum;                         // Save filter output
9 r+ M  D9 k0 l) k( q for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
}
}
2 f0 O1 J9 s3 H4 L- {8 E! u* P
- W% Q5 t; G8 J. W8 P" z! T
" \3 O0 g  Y" ^测试代码：
2 D1 y6 j* R% W  Cfloat w[NUM_TAPS];
void main()
{
    FILE  *fpIn,*fpOut;
    short i;
    char  temp[2];
    float x,  // Input data
& K9 o5 O7 q1 |8 r9 @          y;  // Output data
/ t. V; n8 |5 K" Z, Q6 |/ I  
; p2 s- O# D& [/ @    fpIn = fopen("..\\data\\input.pcm", "rb");
    fpOut = fopen("..\\data\\output.pcm", "wb");
1 @9 ?1 E( E- l4 S! E    if (fpIn == NULL)
    {
        printf("Can't open input file\n");
7 W) h$ S/ b% ^, G        exit(0);
    }

2 \7 o- G/ f2 D" h9 p; s; o    // Initialize for filtering process
3 a. K: ~% e+ n, D) Q( A% B    for (i=0; i<NUM_TAPS; i++)
& v3 j+ b8 n- x* C    {
        w = 0.0;
) c$ z& [* V7 u# h; P5 p    }
    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
! _9 z; ~+ _) _% C! h. X    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
. u3 E- D- v* R7 x        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
$ `/ z: h, P% B. _$ A& s        temp[0] = (char)((short)y&0xff);
; M/ m0 D% q6 a& F$ q1 c        temp[1] = (char)(((short)y>>8)&0xff);
, f$ x* W, i* ?# C+ q        fwrite(&temp, sizeof(char), 2, fpOut);
    }
. h! w3 ~7 F2 h  `( k& K    fclose(fpIn);
    fclose(fpOut);
) @; W. C! X0 ?  `}
, y4 Z$ a' m& j7 M! Z4 E7 ~4 A: E! j

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾