第4章线微分方程组.ppt

2. 一阶线性齐次方程组解的结构 ???? ??? 显见，第一个方程等于第二与第三个方程之和.故上述方程组中仅有两个方程是独立的，即 求它的一个非零解.不妨令 则 于是 故，原方程组的通解为 A:=1 | -1 | -1 , 1 | 1 | 0 , 3 | 0 | 1: Eigenvalues(A): Eigenvectors(A): T1:=0 , -1 , 1: T2:=2/3*I , 1/3 , 1: T3:=-2/3*I , 1/3 , 1: e^t*(cos(2*t)+I*sin(2*t))*T2: simplify(%): U:=-2/3*e^t*(sin(2*t)) , 1/3*e^t*(cos(2*t)) , e^t*(cos(2*t)): V:=2/3*e^t*(cos(2*t)) , 1/3*e^t*(sin(2*t)) , e^t*(sin(2*t)): Y:=C1*e^(1*t)*T1+C2*U+C3*V; with(LinearAlgebra); 2、 矩阵A的特征根有重根的情形 例3 求解方程组 解 系数矩阵为 特征方程为 特征根为 对应的解是 下面求 所对应的两个线性无关解.由定理，其解形如 并且 满足 由于 那么由 可解出两个线性无关向量 可解出两个线性无关向量 将上述两个向量分别代入 中，均得到 为零向量.于是 对应的两个线性无关解是 最后得到通解 with(LinearAlgebra); A:=0 | 1 | 1 , 1 | 0 | 1 , 1 | 1 | 0: E:=1 | 0 | 0 , 0 | 1 | 0 , 0 | 0 | 1: Eigenvalues(A): Eigenvectors(A): Y1:=e^(2*x)*1 , 1 , 1: C:=evalm(A+E): AA:=evalm(C*C): X:=a , b , c: f:=evalm(AA*X): eqn1:=f[1,1]: eqn2:=f[2,1]: eqn3:=f[3,1]: solve({eqn1,eqn2,eqn3},{a,b,c}): R0[1]:=-1 , 1 , 0: R0[2]:=-1 , 0 , 1: R11:=evalm(C*R0[1]): R12:=evalm(C*R0[2]): Y2:=e^(-x)*R0[1]: Y3:=e^(-x)*R0[2]: Y:=C1*Y1+C2*Y2+C3*Y3; 例4 求解方程组 解 系数矩阵为 特征方程为 特征根为 对应的解是 定理3.4 定义3.3 即 n阶齐次方程的n个线性无关解称为基本解组。 即 由基本解组构成的矩阵称为基解矩阵 定理3.5 齐次方程（3.10）必存在基本解组。 定理 刘维尔公式 三、 一阶线性非齐次方程组解的结构 定理3.6 叠加性: 方程(3.8)的任何解与(3.10)的解的和仍是(3.8)的解. 方程(3.8)的任何两个解的差是(3.10)的解. § 3.3常系数线性微分方程组的解法 但是对于一般的方程组(3.10)，如何求出基本解组，至今尚无 一般方法. 求线性齐次方程组(3.10)的通解问题，归结到求其基本解组. 然而对于常系数线性齐次方程组 约当标准型的形式与矩阵A的特征方程的根的情况有关 上述方程也称为常系数齐次方程组(3.24)的特征方程式 (3.21) 1、 矩阵A的特征根均是单根的情形. 2、 矩阵A的特征根有重单根的情形. 易见方程组有n个解 1、 矩阵A的特征根均是单根的情形. 把这n个解代回变换之中，便得到方程组(3.24)的n 个解 即：方程组(3.24)的系数阵A的n个特征根 是方程组(3.24)的一个基本解组 定理 [例題1] 求以下微分方程組的通解 [解] 因此基本矩阵为 故通解为 或以分量表示为 ?例2 试求方程组的通解 解 它的系数矩阵是 特征方程是 先求 对应的特征向量 a, b, c满足方程 , 这三个特征根所对应的特征向量分别是 故方程组的通解是 可求出另两个特征根所对应的特征向量， with(LinearAlgebra); A:=3 | -1 | 1 , -1 | 5 | -1 , 1 | -1 | 3: Eigenvalues(A): Eigenvectors(A): T1:=1 , 1 , 1: T2:=-1 , 0 , 1: T3:=1 , -2 , 1: Y:=C1*e^(3*t)*T1+C2*e^(2*t)*T2+C3*e^(6*t)*T3; 现在考虑复根情形． 因为A