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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
0 Q+ j  ]) i' Y( a: X0 L滤波函数：
#define  NUM_TAPS   48
#define  NUM_DATA   80
0 }. f0 A8 W. P! M+ B; s2 U+ l
8 r% v: f5 U8 H
void floatPointFir(float *x, float *h, short order, float *y, float *w)
% s% W0 W( _3 R/ O$ Y- x2 r{
    short i;
! k& i8 P# m: _' v8 S    float sum;
, c! V; a3 L0 S, z7 G3 W w[0] = *x++;                        // Get the current data to delay line
6 e  I. }5 V1 V) ?2 G for (sum=0, i=0; i<order; i++)      // FIR filter processing
, q  L2 s$ k4 L9 t, ?0 {* ?, L {
  sum += h * w;
}
: U: ]3 ]0 o( a *y++ = sum;                         // Save filter output
. U: C0 ~9 _4 m5 \: M. m for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
1 r  c$ c: ]7 i! x! A. Z+ x% I# J }
3 T% u" R" a+ G}


测试代码：
float w[NUM_TAPS];
4 X, K& W) A4 u' @4 C! g; O0 evoid main()
{
    FILE  *fpIn,*fpOut;
' e; E. d+ M% m    short i;
    char  temp[2];
  ?' H5 @* u* {3 W0 W    float x,  // Input data
9 [. a# \4 d8 u' ^; K7 j5 @          y;  // Output data
  
    fpIn = fopen("..\\data\\input.pcm", "rb");
/ B7 v" l; Q: _/ r' t8 ^& }8 T& C( }    fpOut = fopen("..\\data\\output.pcm", "wb");
$ I2 B) v" C6 m9 E    if (fpIn == NULL)
    {
        printf("Can't open input file\n");
        exit(0);
! W+ G7 y; L' ?0 t9 v# F: B    }

    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
: L, u- j/ q6 _3 @3 V    {
        w = 0.0;
; S! p- X* X; ?0 z- o  D8 X) ?    }
    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
    {
# g/ T5 V' z3 l/ Z        x = (float)((temp[1]<<8)|(temp[0]&0xff));
8 X/ |, {' ], U% X' g/ y& g# i( l' C        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
        temp[0] = (char)((short)y&0xff);
2 m3 `5 ^) u1 U8 K        temp[1] = (char)(((short)y>>8)&0xff);
        fwrite(&temp, sizeof(char), 2, fpOut);
6 ~( f0 K+ n6 _8 [! B    }
( _2 |& ^+ X) _0 W  B    fclose(fpIn);
    fclose(fpOut);
}

6 z: {: @/ X; |
# P& B' G' D6 p- G8 c; D

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾