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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
. {$ p# V6 z, D  k0 A滤波函数：
& w7 J9 F. }3 S: i6 O# k, F#define  NUM_TAPS   48
, F! J/ B; X  N7 R/ P#define  NUM_DATA   80
5 }4 V5 a% S3 r# h: n/ ]7 w% \
& x" L4 n6 P) f( A
void floatPointFir(float *x, float *h, short order, float *y, float *w)
) z3 }9 O9 l! o, n- p{
4 V- v  R; n8 h2 }. G$ F    short i;
    float sum;
# O5 \7 @5 [7 Q w[0] = *x++;                        // Get the current data to delay line
) R  u+ E* d2 A0 m8 B for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
: ?$ E# T& n: {+ B  sum += h * w;
}
*y++ = sum;                         // Save filter output
4 b2 ?0 v+ b% d& V# v for (i=order-1; i>0; i--)           // Update data delay line
{
  w = w[i-1] ;
& B0 q: ~9 u9 K- G0 B7 j$ K. _: v }
4 O" G3 O/ y: K}


5 n# ?6 t; X  i5 a% t* E' \7 v测试代码：
float w[NUM_TAPS];
& H  x& L7 i5 g. ?void main()
( ]  H9 q2 [# E{
7 p. r: V4 V0 t4 i3 [3 A# I    FILE  *fpIn,*fpOut;
1 F" |  c6 x. s7 ?0 U: x: ?4 ~    short i;
    char  temp[2];
! c2 R8 |7 O7 l! Y/ {/ c    float x,  // Input data
          y;  // Output data
* T- h/ J$ l5 Z/ R  
    fpIn = fopen("..\\data\\input.pcm", "rb");
    fpOut = fopen("..\\data\\output.pcm", "wb");
    if (fpIn == NULL)
4 I- [4 ]6 X, e6 |! ]" @) R. H  N    {
        printf("Can't open input file\n");
' K, r: n. ?- v  S        exit(0);
8 E3 L  h4 C5 N$ \# Y0 r    }

    // Initialize for filtering process
% |# l3 r/ l) b4 v1 Z1 E3 _    for (i=0; i<NUM_TAPS; i++)
1 _. J1 H) b5 z# z0 B; Z    {
        w = 0.0;
4 n6 w! a* [; e  r7 e    }
1 x2 W2 {$ Z- N    // Begin filtering the data
8 r4 P% e3 O+ r( L    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
1 L( l( `% |, T; @( i0 c+ T/ l    {
9 \; e0 h1 }" t        x = (float)((temp[1]<<8)|(temp[0]&0xff));
( _0 \" t1 v# \# }        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
+ I1 i3 I& h+ g        temp[0] = (char)((short)y&0xff);
        temp[1] = (char)(((short)y>>8)&0xff);
0 ]/ ^  R# l3 e9 f% a  W: Q( W        fwrite(&temp, sizeof(char), 2, fpOut);
  a% W. N0 ^% _' a5 V4 X; c  ^    }
    fclose(fpIn);
    fclose(fpOut);
. n  u1 Q9 V* q4 x8 K}

; q: d) i  s. B. D/ Z
- u* A  C0 \2 [* u' w3 K

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾