线性微分方程组的矩阵分块解法.pdfVIP

维普资讯 ， — - -． 1996年 青海师范大学学报 (自熊科学版) 1996 第 2期 JournalofQinghaiNormalUniversit(NaturalScience) No．2 (2) 线性微分方程组的矩阵分块解法 J_] 左光纪 Q!-1|＼ (-It海 民族学院) A摘要本文讨论矩阵接可约性或嵇的分块问题．对可约阵或奇异阵分别给出了线性赣分方程组的矩阵分 块解法．还作了物理解释 ． 胖的生墨柏图 解I ， 1引言 ”。， EJ阵，『” ，， 设有线性微分方程组 xI：A (f)x 其中 阶函数矩阵 A(f)在 tECa，b3上连续 周知 ，求 (1)的通解归结为求它的基解矩阵。当A(f)为常矩阵时，可用矩阵特征值法求基 解矩阵；但是，对高阶矩阵求特征值仍有困难，当A(f)为变矩阵时尚无求解方法c1]。看来 ，求 解 (1)的难点在于 A(f)的阶太高。许多由实际问题得出的方程组 (1)中，A(f)往往含有大量零 元素 (所谓稀疏矩阵)，或者不是满秩的。这时，可否将 A(f)分块使求解问题降阶?本文将结合 矩阵的组合性质和代数性质，对A(f)可约或奇异两种情形，讨论 (1)的矩阵分块解法。 2 矩阵的不可约分块与图论方法 对稀疏矩阵，零元素的位置分布是至关重要的。文C23把仅与矩阵的零元位置分布有关，而 与非零元的具体数值无关的那些性质称为矩阵的组合性质 。矩阵的可约性就是一种组合性质， 它在置换相似 (即行列同步交换)下不变。 。 定义 l 设A 为 阶方阵，若存在 阶置换矩阵P，使 P P 具有如下分块形式 r ． 0 1 PAP：Ll ^．1，AIA2为非空方阵； 2J 则称 A是可约的；否则 ，称A是不可约的。 定义 2 设A=Ca]是 ≈阶方阵，A的伴随图G(A)=(，E)是以 ={1，2，…，n}为点集 ， 以E= {(，)Iaij≠O}为弧集的一个有向图。 命题 1 方阵A为不可约的充要条件是伴随图G(A)强连通 。 命题 2 设 A是可约方阵，则存在置换阵P，使 PAP 具有分块下三角形： PAP 一 维普资讯 青海师范大学学报 (自然科学版) 1996l 其中对角块 ^ ，^ “，A，为不可约子方阵。 我们称 (2)为方阵 A的不可约块标准形。下述例子显示了化方阵为不可约块标准形的图 论方法 ，用到的术语可参见文C33C4]。 例 1 求矩阵^的不可约块标准形 0 O 0 0 0 0 5 0 O 0