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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
# ?! K  X7 U2 h' N% J- Q1 i& T滤波函数：
#define  NUM_TAPS   48
#define  NUM_DATA   80

8 n% s5 r  s' E$ ^" |
void floatPointFir(float *x, float *h, short order, float *y, float *w)
{
    short i;
    float sum;
" h9 k/ L) H# F/ z: i w[0] = *x++;                        // Get the current data to delay line
$ A3 _1 |' c3 \8 m) n1 e6 y% } for (sum=0, i=0; i<order; i++)      // FIR filter processing
+ x5 A  S4 {0 }  Q+ m, b {
0 n( O; Z! |9 n+ L" z' e  sum += h * w;
}
& t* e' B0 L3 p$ L. I3 W0 U5 f! g *y++ = sum;                         // Save filter output
- X: X2 n) q# y* n) g3 Q& X for (i=order-1; i>0; i--)           // Update data delay line
{
' s" R7 Q$ I( z9 y2 T$ J1 x  w = w[i-1] ;
, |+ |! P0 c. E6 w- H9 z' z5 Y }
$ O+ m6 L2 D7 V}


! i5 ]* ^; m/ h# D4 j测试代码：
2 ?' k- w% T( Y# M4 S1 Y6 |/ U4 ?float w[NUM_TAPS];
0 Z& X$ ]( s) i3 hvoid main()
5 E! f- J9 j  c' L9 B/ x" T8 A7 M{
    FILE  *fpIn,*fpOut;
1 Q4 }0 `; y" i( m2 Z0 t    short i;
7 d* D* n9 R% l  y1 T( ~6 K    char  temp[2];
    float x,  // Input data
& p/ ^& A7 }( ~* X9 ^4 ]5 v          y;  // Output data
  
    fpIn = fopen("..\\data\\input.pcm", "rb");
    fpOut = fopen("..\\data\\output.pcm", "wb");
    if (fpIn == NULL)
    {
8 B- M5 [9 N- n* P" f+ Y! `        printf("Can't open input file\n");
  Z8 V& \7 O+ Q2 Q  I: ?3 w2 c3 T, }) J        exit(0);
    }

    // Initialize for filtering process
7 [# x: O' O1 j% M( F- u    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
+ b4 b/ Q5 M+ Q2 Y! F    }
6 k8 h+ H" ^- {; N* F' b0 n- Z    // Begin filtering the data
5 Q7 m; M0 O* d$ [( H- Z% o    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
9 U+ Q% s* J& d% |0 S) [# L+ N- u    {
4 y1 U0 w( D# y1 F* b        x = (float)((temp[1]<<8)|(temp[0]&0xff));
        // Filter the data x and save output y
3 @2 Q/ O# i2 a* ]. J6 w        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
  \# W% o, D+ F0 u" u9 u% x        temp[0] = (char)((short)y&0xff);
9 X( A8 P+ d# ], ~' l8 [        temp[1] = (char)(((short)y>>8)&0xff);
        fwrite(&temp, sizeof(char), 2, fpOut);
6 E9 h- y& G  w+ b, T2 o% h    }
    fclose(fpIn);
" |( W. M+ _# A+ f  I- Z7 X    fclose(fpOut);
: S' W0 q% u6 q) M7 Y}
8 o  j/ m* H* i, r6 q* \/ @

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾