直接型FIR滤波器的C语言实现

设输入数据x[N],输出数据y[N],滤波器系数h[n]

1.直接法(由y(m)=h(0)*x(m)+h(1)*x(m-1)+...+h(N-1)*x(m-n-1));

void fir(short x[], short h[], short y[])
{
  int i,j;
  long long sum;  
for (j = 0; j < N; j++)
{
sum = 0;
for (i = 0; i < n; i++)
sum += x[j-i] * h[i];
y[j] = sum >> 15;
}
}

乘法器使用次数：N*n

2.逆推法：

void fir(short x[], short h[], short y[])
{
  int i,j;
  long sum;  
for (j = 0; j < n; j++)
{
for (i = 0; i < N; i++)
{
sum = 0;
sum = h[j] * x[i]
y[i] += sum >> 15;
}
}
}

乘法器使用次数：N*n

3.倒序法：（输入输出可以是同一量）

void fir(short x[], short h[], short y[])
{
  int i,j;
  long long sum;  
for (j = N; j > 0; j–)
{
sum = 0;
for (i = n; i > 0; i–)
sum += x[j-i] * h[i];
y[j] = sum >> 15;
}
}