连续正弦函数状态方程离散化

假设连续的正弦函数状态方程为：

dx/dt = Ax + Bu
y = Cx + Du

其中，`x` 是状态向量，`u` 是输入向量，`y` 是输出向量，`A`、`B`、`C`、`D` 是系数矩阵。为了将该连续系统离散化，可以使用 Euler 离散化方法，其基本思想是将时间区间划分为若干个小区间，并在每个时间点上对连续系统进行近似。

具体而言，对于上述连续系统，Euler 离散化方法的实现如下：

T = 0.1; % 时间间隔
Ad = eye(size(A)) + A * T; % 状态转移矩阵
Bd = B * T; % 输入转移矩阵
Cd = C; % 输出转移矩阵
Dd = D; % 直通矩阵

其中，`Ad`、`Bd`、`Cd`、`Dd` 分别是离散化后的状态、输入、输出、直通矩阵。

需要注意的是，Euler 离散化方法的精度不高，当时间间隔 `T` 较大时可能会引入较大的误差。因此，如果需要更高精度的离散化方法，可以尝试使用其他方法，比如 ZOH（Zero-Order Hold）离散化方法或 Tustin（Bilinear）离散化方法等。