状态转移矩阵(STMs)是线性微分方程的齐次部分的解。 线性微分方程是指: 它的齐次部分指: 解齐次部分要比解整个方程容易,齐次部分的解可以用来解非齐次部分。 Φ(t)被称为齐次方程的基本解,如果它满足下面的条件: 其中,I是nxn的单位矩阵。 对于任何可能的x(0)来说,都有如下表达式成立: 所以,如果我们想知道从某个0时刻开始的x(t),只需要知道 Φ(t)即可: x(t)= Φ(t)x(0) 如果是离散线性系统,则有: 对于离散线性系统来说, Φ是 它的状态转移矩阵(STM) ,C是输入耦合矩阵。 为什么 Φ被称为状态转移矩阵?因为已知了x(0)和 Φ(t),我们就能知道x(t)。 Φ将0时刻的状态,转移到了t时刻。 假如不是从0时刻开始,而是从 τ时刻转移到t时刻?我们有如下表达式成立: Φ有下面这些特征: 这些特征不难想象,使用实际的物理现象也很容易印证。 例:如何求解弹簧阻尼振荡器的状态转移矩阵 Φ? 解:对于一个弹簧阻尼振荡器来说,我们已有它的F矩阵表达式: F表达式由牛顿力学导出,x表示位移,x的导数表示速度。 对于一个简单的 线性 动力系统,获取F矩阵非常容易,如果不会的话,让牛顿的苹果砸一砸先^_^ 然后,求解 λI-F的行列式: 乃们还记得行列式是怎么求的么? 对于2x2和3x3矩阵分别有: 所以,对于这里的 λI-F有: 此二次多项式的根: 接下来,我们使用 δ来表示 位移,则x矩阵可以表示为: δ可以写成通用形式: 对于一个欠阻尼系统,有下面的表达式成立: 阻尼为0的振荡器,会永远振荡下去。 过阻尼的振荡器,会因为阻力过大而回不到平衡位置。 欠阻尼的振荡器,会因为阻力很小而越过平衡位置,但是振动幅度越来越小。 欠阻尼系统的δ是非实数,所以又可以写成: 这里的a和b是实数。 再获取x矩阵: 我们可以使用初始化的值求解a和b: 最后得出来状态转移矩阵 Φ: 所以,求解 弹簧阻尼振荡器的状态转移矩阵的基本步骤是: 先通过牛顿力学方程建立F矩阵表达式 ; 然后求解 λI-F的特征值,以及特征值的二次项平方根; 接着使用特征值表示x矩阵,通过初始化条件求解参数a和b; 最后代入a和b,求得状态转移矩阵。