利用矩阵函数的导数公式求解一阶常系数微分方程组的解

考虑以下形式的线性常系数微分方程组：

其中， 是 n×1 的状态向量，A 是 n×n 的常系数矩阵。

步骤

1. 特征值和特征向量 ：求出矩阵 A 的特征值和特征向量。如果 A 是对角化的，那么可以将其对角化为 ，其中 D 是一个对角矩阵，包含 A 的特征值，P 包含对应的特征向量。

2. 一般解的形式：由于矩阵微分方程的解可以写成矩阵指数的形式，所以解为

   

   其中，表示矩阵 A 的指数。

3. 计算矩阵指数 ： 对于对角化矩阵 A ，矩阵指数 可以写成

   

   其中是对角矩阵，每个元素为 （ 是 A 的特征值）。

4. 解的求得 ：通过矩阵指数 和初始条件 x(0) 可以求得解 x(t) 。

注：的具体求法参考文章：求解矩阵函数值的方法

示例

假设有一个 2x2 的常系数矩阵微分方程：

步骤 1：求矩阵 A 的特征值和特征向量

给定矩阵为：

1. 求特征值：特征值 λ 满足 det⁡(A−λI)=0 。

   

   解此方程得到特征值：
   

2. 求特征向量：

   对于 ，我们解 ：
   

   这给出特征向量。

   对于 ，我们解 ：
   

   这给出特征向量 。

因此，矩阵 A 的特征值和特征向量分别为：

,

步骤 2：写出矩阵 P 和对角矩阵 D

将特征向量按列排成矩阵 P，则：

对角矩阵 D 包含特征值：

步骤 3：计算矩阵指数

由于，所以

首先，计算 ：

然后求 ：

因此，

进行矩阵乘法得到

步骤 4：写出通解

假设初始条件为 ，则解为

结论

通过矩阵指数的方法，我们得到了该微分方程组的通解：