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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
9 h- _6 w8 r9 T3 K1 N$ ^滤波函数：
0 N% r& k6 q0 }  F#define  NUM_TAPS   48
% d( b. Q  }" N+ p6 I8 K# O9 S- N#define  NUM_DATA   80
* R* V. u6 Y( Y+ m
! _8 c8 x. l- K8 V8 x  ?1 w
void floatPointFir(float *x, float *h, short order, float *y, float *w)
{
+ Q4 k# H6 J( m1 \/ F    short i;
3 f: `6 `* W0 ^+ {; F    float sum;
w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
2 b+ }* \) z+ ~7 a  D# E {
  sum += h * w;
0 t1 \9 ?  U7 U }
0 A. H1 f7 v' W4 p+ ]0 s4 o8 j *y++ = sum;                         // Save filter output
for (i=order-1; i>0; i--)           // Update data delay line
* f7 P7 r% ?& z) n {
  w = w[i-1] ;
}
}

' B) H- k: U: k6 [( u
" l  I6 q5 r9 B. R3 S1 ]9 m+ J测试代码：
! Y2 A, P. `. z. Afloat w[NUM_TAPS];
void main()
7 ~( K1 e- b# i7 a* e! W$ p; j{
6 @, f7 S3 I- G3 N( |" N    FILE  *fpIn,*fpOut;
    short i;
    char  temp[2];
4 |! h2 ?1 Z5 W! f2 `$ _    float x,  // Input data
1 m6 t4 m! O- [3 x9 r4 f0 ^. c; T          y;  // Output data
  
/ a  I& ?. [2 L    fpIn = fopen("..\\data\\input.pcm", "rb");
+ J# D2 x7 }$ M" @$ b4 W' @# p( t    fpOut = fopen("..\\data\\output.pcm", "wb");
5 a6 x8 U& [2 H4 B; _1 a9 F    if (fpIn == NULL)
; E& O4 {, N) Y    {
        printf("Can't open input file\n");
) m1 h( p) ~/ N# j        exit(0);
; A* X5 g& c6 B. l1 x; k    }
+ {4 y3 m8 A; T; J
) s: e+ X; [' ^1 ]. d/ v4 A    // Initialize for filtering process
    for (i=0; i<NUM_TAPS; i++)
1 W  F1 h& C' q& Q8 s0 l, s    {
. |& P3 {9 s. p% _. w0 z: B        w = 0.0;
  I( V  }5 \9 Z: K9 l    }
; Q" {( ]' t7 }6 ~1 f    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
8 D8 ?( W' }. u/ z# P5 q" Y. i    {
& w0 E, w; y6 X5 O5 P        x = (float)((temp[1]<<8)|(temp[0]&0xff));
        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
        temp[0] = (char)((short)y&0xff);
" R7 c  [) _9 O* s, }        temp[1] = (char)(((short)y>>8)&0xff);
% }* L; X) j& h$ D! i  H1 a        fwrite(&temp, sizeof(char), 2, fpOut);
1 ]+ W6 b2 d6 k6 ~: P, O' n    }
    fclose(fpIn);
+ V3 c1 @+ g9 S) N8 ^4 Z6 B& ^- i& s    fclose(fpOut);
}

* x. t# C0 j; ~/ Z! n2 u7 P
# f0 v1 U; [& n& \) z$ n

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾