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

freebigfish

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

% K  x9 A9 u1 \+ L4 R0 [  |' \4 q
7 v2 ?# F. m0 @4 \7 Q/ ]void floatPointFir(float *x, float *h, short order, float *y, float *w)
{
8 Q; b) r; f# P    short i;
    float sum;
w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
, Z, z8 X8 E0 x {
  sum += h * w;
5 `, o6 a5 ]! W. Q }
*y++ = sum;                         // Save filter output
7 W8 m4 F8 Y. H3 X/ \# J for (i=order-1; i>0; i--)           // Update data delay line
6 q9 c3 f+ _4 y* I0 {. H# a$ ^3 j {
# ^+ l% _: N" ^# ^  w = w[i-1] ;
9 U* U' C( A/ @, R8 U8 A) q }
}
0 w  i! C  r* r+ G. W

$ z' I) f* `3 D" |6 Q% f  h# g测试代码：
float w[NUM_TAPS];
, U! M, a  Z' {( e: a' I% K3 c) svoid main()
{
# l7 t# ^2 L8 V& v    FILE  *fpIn,*fpOut;
    short i;
    char  temp[2];
    float x,  // Input data
, ~. q, _& D' x8 |& p          y;  // Output data
  
( w" F& t) e+ t: X, I    fpIn = fopen("..\\data\\input.pcm", "rb");
    fpOut = fopen("..\\data\\output.pcm", "wb");
, y' y8 o5 }2 ~8 }    if (fpIn == NULL)
    {
        printf("Can't open input file\n");
' G8 ]4 {6 d, R: Q- W! z5 l- V        exit(0);
" s; X( G' y9 @5 o7 h    }

    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
    }
4 q* W, }9 u. {2 `& H6 |    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
* V4 X+ M6 s. @8 |8 V+ [6 {    {
  Y( w5 v# _3 J6 m+ R        x = (float)((temp[1]<<8)|(temp[0]&0xff));
% S, Y& c; o2 Y6 `        // Filter the data x and save output y
$ a- b) _( G- Z9 E5 z        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
        temp[0] = (char)((short)y&0xff);
        temp[1] = (char)(((short)y>>8)&0xff);
. N. [( b6 G/ }9 I2 M        fwrite(&temp, sizeof(char), 2, fpOut);
    }
2 O9 [8 Y* \: X    fclose(fpIn);
    fclose(fpOut);
}
% F4 G8 A  y  _; O. j

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾