考虑以下形式的线性常系数微分方程组:
其中, 是 n×1 的状态向量,A 是 n×n 的常系数矩阵。
步骤
-
特征值和特征向量 :求出矩阵 A 的特征值和特征向量。如果 A 是对角化的,那么可以将其对角化为 ,其中 D 是一个对角矩阵,包含 A 的特征值,P 包含对应的特征向量。
-
一般解的形式:由于矩阵微分方程的解可以写成矩阵指数的形式,所以解为
其中,表示矩阵 A 的指数。
-
计算矩阵指数 : 对于对角化矩阵 A ,矩阵指数 可以写成
其中是对角矩阵,每个元素为 ( 是 A 的特征值)。
-
解的求得 :通过矩阵指数 和初始条件 x(0) 可以求得解 x(t) 。
注:的具体求法参考文章:求解矩阵函数值的方法
示例
假设有一个 2x2 的常系数矩阵微分方程:
步骤 1:求矩阵 A 的特征值和特征向量
给定矩阵为:
-
求特征值:特征值 λ 满足 det(A−λI)=0 。
解此方程得到特征值:
-
求特征向量:
对于 ,我们解 :
这给出特征向量。
对于 ,我们解 :
这给出特征向量 。
因此,矩阵 A 的特征值和特征向量分别为:
,
步骤 2:写出矩阵 P 和对角矩阵 D
将特征向量按列排成矩阵 P,则:
对角矩阵 D 包含特征值:
步骤 3:计算矩阵指数
由于,所以
首先,计算 :
然后求 :
因此,
进行矩阵乘法得到
步骤 4:写出通解
假设初始条件为 ,则解为
结论
通过矩阵指数的方法,我们得到了该微分方程组的通解: