我有一个在 C 语言中实现 FIR 滤波器的家庭作业,我想知道你是否认为我正确理解了作业。我认为可以解决问题的程序是:
#include <stdio.h>
float FIRfloats[5];
void floatFIR(float newsample)
{
int i;
float sum=0;
FIRfloats[0]=newsample*0.0299;
FIRfloats[1]=FIRfloats[2]*0.4701;
FIRfloats[2]=FIRfloats[3]*0.4701;
FIRfloats[3]=FIRfloats[4]*0.0299;
/* sum */
for(i=0;i<5;i++)
{
sum=sum+FIRfloats[i];
}
printf("Sum: %fn", sum);
}
int main ()
{
float n=0.0f;
while (scanf("%f", &n) > 0)
{
floatFIR(n);
}
return 0;
}
规格是
系数应为 c0=0.0299、c1=0.4701、c2=在新样本 xk 到达之前,旧样本被转移到 右边,然后每个样本在 计算结果 YK,即所有缩放样本的总和
0.4701、c3=0.0299。
你认为我正确解决了作业吗?我认为这似乎太容易了,因此我想知道。
恐
怕问题中提供的实现不会提供正确的结果。
在具有 4 个系数的 FIR(有限脉冲响应)滤波器中,输入序列 (x) 的输出序列 (y) 为:
y[t] = c0*x[t] + c1*x[t-1] + c2*x[t-2] + c3*x[t-3]
因此,实现应类似于:
/* add includes (stdio.h and whatever else you'll need...) */
float floatFIR(float inVal, float* x, float* coef, int len)
{
float y = 0.0;
for (int i = (len-1) ; i > 0 ; i--)
{
x[i] = x[i-1];
y = y + (coef[i] * x[i]);
}
x[0] = inVal;
y = y + (coef[0] * x[0]);
return y;
}
main(int argc, char** argv)
{
float coef[4] = {0.0299, 0.4701, 0.4701, 0.0299};
float x[4] = {0, 0, 0, 0}; /* or any other initial condition*/
float y;
float inVal;
while (scanf("%f", &inVal) > 0)
{
y = floatFIR(inVal, x, coef, 4);
}
return 0;
}
这在同一循环中完成移位和乘法(这不会影响结果 - 只是更有效。如果你想完全遵循规范,你可以像这样更改 floatFir:
float floatFIR(float inVal, float* x, float* coef, int len)
{
float y = 0.0;
for (int i = (len-1) ; i > 0 ; i--)
{
x[i] = x[i-1];
}
x[0] = inVal;
for (int i = 0 ; i < len ; i++)
{
y = y + (coef[i] * x[i]);
}
return y;
}