IIR 滤波测试

最新推荐文章于 2024-04-11 13:04:02 发布

原创最新推荐文章于 2024-04-11 13:04:02 发布 · 384 阅读

2 ·

CC 4.0 BY-SA版权

算法专栏收录该内容

7 篇文章

订阅专栏

本文详细介绍了如何使用MATLAB生成滤波器参数，并将其转换为C语言代码进行实时信号处理。通过对比静态浮点测试函数与MATLAB生成的IIR滤波器，验证了两者在信号处理效果上的高度一致性。

1.使用matlab产生产生参数。

static float test_f(float input)
{
	static float xv[2+1], yv[2+1];
	xv[0] = xv[1]; xv[1] = xv[2]; 
	xv[2] = input/ 1.949595788e+02;
	yv[0] = yv[1]; yv[1] = yv[2]; 
	yv[2] =   (xv[0] + xv[2]) + 2 * xv[1]
				 + ( -0.8079495914 * yv[0]) + (  1.7874325180 * yv[1]);
	return yv[2];
}

static float iir_filter_matab(float input)
{
	//y[n] = 1*x[n] + 2*x[n-1] + 1*x[n-2] - (-1.1429805025)*y[n-1] - (0.412801596)*y[n-2]
	// matlab 都是二阶滤波级联
	double gain =  0.005129268366107;
	float table_N[3] = {1,2,1};
	float table_D[3] = { 1,   -1.787432517956,   0.8079495914209};
	static float xv[3], yv[3];
	xv[2] = input*gain;
	xv[0] = xv[1]; xv[1] = xv[2]; 
	yv[0] = yv[1]; yv[1] = yv[2]; 
	yv[2] = table_N[0]*xv[2]+table_N[1]*xv[1]+table_N[2]*xv[0] - table_D[1]*yv[1] - table_D[2]*yv[0];
	
	return yv[2]; 
}

y[n]=b0⋅x[n]+b1⋅x[n−1]+b2⋅x[n−2]−a1⋅y[n−1]−a2⋅y[n−2]

matlab 中高阶滤波是用二阶滤波级联。

测试和网页版生成的代码滤波效果基本上一样。