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

本文介绍如何将连续时间系统的传递函数离散化,并使用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];
}



评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值