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

freebigfish

各位大神，我有一个疑问。fir的输出序列长度按道理是输入序列的长度+滤波器的阶数。可最近看了一个代码，它的输出序列长度冒似与输入序列长度相同，恳请指点一下。
滤波函数：
#define  NUM_TAPS   48
#define  NUM_DATA   80
( b: u$ Z3 f6 X% ^
3 C& W3 c( E& u1 H. T/ N$ {; _
& q6 {( T, Y- B3 Jvoid floatPointFir(float *x, float *h, short order, float *y, float *w)
1 c) Y) z( U8 k8 B{
8 L' \+ @8 x* F/ p    short i;
7 X2 L% s0 B; B    float sum;
w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
, T# S  R2 }+ q8 e" |  sum += h * w;
4 @# p! c1 T% [; |: k: w& D% W }
*y++ = sum;                         // Save filter output
  ^$ i# K- Z8 _7 Z  I4 v for (i=order-1; i>0; i--)           // Update data delay line
9 ~% v; G2 _" X+ v* s. ~9 m {
  w = w[i-1] ;
8 n1 e2 R: Z7 M% a1 U: U }
, Y+ f' r2 M: @/ T& i' }6 J4 h}

' V% ~, d. Z7 q8 [. F% c, ]3 P
  N* i! G: h) _8 L3 t5 h# d; P. E测试代码：
float w[NUM_TAPS];
4 b3 `1 a' r- b" Qvoid main()
{
/ d3 t3 y7 @  b' G4 W! d    FILE  *fpIn,*fpOut;
! [; W8 v1 }% x( [( I9 b4 q    short i;
    char  temp[2];
! S2 H1 S7 p2 A! w$ c1 {& v0 H    float x,  // Input data
* G1 K8 S/ |6 r  B          y;  // Output data
  
    fpIn = fopen("..\\data\\input.pcm", "rb");
    fpOut = fopen("..\\data\\output.pcm", "wb");
    if (fpIn == NULL)
. |. U8 R7 D2 M) W! x1 M    {
        printf("Can't open input file\n");
        exit(0);
! G2 ?; ]5 \8 G    }

    // Initialize for filtering process
1 d( m( i/ l$ ?3 ?9 ~    for (i=0; i<NUM_TAPS; i++)
    {
        w = 0.0;
! v. Z* o& d: `6 O    }
& C6 w- h2 R( f( B* V+ y$ J9 h1 t    // Begin filtering the data
    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
    {
9 s& p% s7 Y0 |" O- ~* d6 U0 [& X        x = (float)((temp[1]<<8)|(temp[0]&0xff));
+ |# ?" m$ A, e        // Filter the data x and save output y
        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
        temp[0] = (char)((short)y&0xff);
        temp[1] = (char)(((short)y>>8)&0xff);
9 v) }3 R+ c5 f. l' x        fwrite(&temp, sizeof(char), 2, fpOut);
    }
    fclose(fpIn);
- G  R! j1 E  B" B& ^! R    fclose(fpOut);
4 f) W" z8 ?; z}


, Z9 \. ]) c) H6 q; Q, _" F

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾