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

freebigfish

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

" o' ~% i, P, p
# D8 [( a; `+ k+ [/ ?void floatPointFir(float *x, float *h, short order, float *y, float *w)
- k+ [3 H6 H) |% @! _. v* M3 w{
    short i;
    float sum;
) f$ t& ~, a: y+ }- K, @ w[0] = *x++;                        // Get the current data to delay line
for (sum=0, i=0; i<order; i++)      // FIR filter processing
{
, {0 z% E0 q* L  sum += h * w;
}
" n1 P7 J4 k2 t% a5 q *y++ = sum;                         // Save filter output
for (i=order-1; i>0; i--)           // Update data delay line
, i1 |2 j4 Y& H  l6 n1 ` {
0 ~: Y9 H4 ?8 r3 [: b  w = w[i-1] ;
! m; C5 R3 u& p) _$ y# S }
}
: n, g: P- O( E6 ]7 M
: g0 u* o4 {# H/ W" o; o: y4 ~5 e
5 U9 c9 D# b% c9 t+ p测试代码：
3 I5 ~/ L: d( d' h, L" s' L& Rfloat w[NUM_TAPS];
void main()
{
    FILE  *fpIn,*fpOut;
    short i;
  P& g% T3 s' l9 |5 m0 R3 u    char  temp[2];
    float x,  // Input data
( _) A9 `! w+ v! X" |! L' V, @          y;  // Output data
; B; Y  A' u/ c  
4 t' X" L* h' ~+ h4 H    fpIn = fopen("..\\data\\input.pcm", "rb");
, ~; o/ _5 `, Q2 B5 z    fpOut = fopen("..\\data\\output.pcm", "wb");
    if (fpIn == NULL)
+ h; K9 j8 r' r9 H, G4 h    {
3 w3 J& o, f5 ]3 ~- Y) N# |        printf("Can't open input file\n");
        exit(0);
    }

, l% V' W2 H# p( j- l! Z    // Initialize for filtering process
# B; z3 V3 l9 G- S1 p  H2 g; ?    for (i=0; i<NUM_TAPS; i++)
    {
8 |  g' L# D3 S' R: s2 |        w = 0.0;
7 `+ g! p' d+ |8 M8 h4 i    }
    // Begin filtering the data
+ n% R) r2 I5 w6 ~7 v  G) H    while (fread(&temp, sizeof(char), 2, fpIn) == 2)
( T; k& r$ g# P    {
        x = (float)((temp[1]<<8)|(temp[0]&0xff));
        // Filter the data x and save output y
: r/ q9 A, c* r, a        floatPointFir(&x, firCoefFloatingPoint, NUM_TAPS, &y, w);
* M$ `3 V  W) b7 [8 Z5 R7 W3 p        temp[0] = (char)((short)y&0xff);
6 x! S$ [8 N& g; U        temp[1] = (char)(((short)y>>8)&0xff);
        fwrite(&temp, sizeof(char), 2, fpOut);
# h% I* l/ ?0 O# K    }
6 o" J8 _( i% J4 M/ E    fclose(fpIn);
    fclose(fpOut);
}
- ~; L) M+ Q2 g0 \& T
) t+ r5 ^4 k+ S# F! e3 w

helendcany

路过。

freebigfish

怎么看不到回复啊

freebigfish

大神们都忙，不屑一顾