基于matlab用c语言实现传递函数

最新推荐文章于 2024-02-22 11:02:56 发布

原创最新推荐文章于 2024-02-22 11:02:56 发布 · 9.1k 阅读

78 ·

CC 4.0 BY-SA版权

文章标签：

#c语言 #传递函数 #离散化 #matlab

本文介绍如何将连续时间系统的传递函数离散化，并使用C语言进行实现。通过具体例子展示了离散化过程及C语言代码实现细节。

以传递函数sys = tf([1, 1], [1, 1, 1])为例：

1、生成系统

sys = tf([1, 1], [1, 1, 1])

2、设置离散时间

ts = 0.001

3、离散化

dsys = c2d(sys, ts, 'method')

4、查看系数

[num, den] = tfdata(dsys, 'v')

num = 1.0e-03 * 0    1.0000   -0.9990
den = 1.0000   -1.9990    0.9990

5、c语言实现，函数ts时间调用一次

float function(float input)
{
        /* n为传递函数分子分母最高阶次 */
	const uint8_t n = 2;
	/* 分子num */
	const float b[n + 1] = {1.0e-03 * 0, 1.0e-03 * 1.0000, -1.0e-03 * 0.9990};
        /* 分母den */
	const float a[n + 1] = {1.0000, -1.9990, 0.9990};

	static float x[n + 1] = {0};
	static float y[n + 1] = {0};
	
	for (uint8_t i = n; i > 0; i--)
	{
		y[i] = y[i - 1];
		x[i] = x[i - 1];
	}
	x[0] = input;
	y[0] = 0;
	
	for (uint8_t i = n; i > 0; i--)
	{
		y[0] = y[0] + b[i] * x[i];
		y[0] = y[0] - a[i] * y[i];
	}
	
	return y[0] + b[0] * x[0];
}