直接型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;
}
}