求解线性方程组的初参数方法.pdfVIP

维普资讯 上 海 交 通 大 学 学 报 第 28卷 第 3期 JOURNALOFSHANGHAIJIAOTONGUNIVERSITY Vol 一 ／0 求解线性方程组的初参数方法 王纪林 周 钢 p If5f、 (应用数学系) 摘 要 本文提 出一种求解线性方程组 的直接解法．该方法把较大方程组 的求解化 A 为可并行 的较小方程组求解 ，可适用于大系统 问题或高维数值解 问题 的求解 ，也可处 理混合 问题的不 同区域．方法 的实施与方程组的可解性无关，这对于求解变系数不定 常 问题是有意义 的， 关键词 、 垒 查堡垒，塑查塾之墨，!苎±丝一 中圉法分类号 0241，6 0 引 言 初参数方法的思想来源于结构力学 ，是分析子结构的静力和动力性态时常用的方法．把初 参数概念在数学上加以提取，并应用到线性方程组的求解问题上 ，得到一种新的直接求解法 ， 称之为初参数方法．方法中对方程组的分解有很大的 自由度，这对于提高并行计算的效率和处 理混合问题都有着实际意义． 1 初参数方法的原理和计算步骤 考虑下列问题 ：求 ∈ ，使Az=6．其中A为 阶复数方阵 b~o为 维复数列 向量 为说明初参数方法，设 A 一 ㈦ ，一 其 中A 为r阶方 阵 (1≤r≤ 一1)，det(A。)≠0，z、b。都 是r维 列 向量 ，则 问题 化 为求 z一 ( ，Y)，使 收稿 日期：1993—4—30 维普资讯 王纪林等；求解线性方程组 的初参数方法 lol fAlz+ A2Y= b1 (1a) A【az+A^Y—b2 (1b) 称式 (1a)为 r阶选取方程．任取 Y的初值Y。，由式 (I1a)解得 X。= AI (6】一 A2Y。) (2) 把 。；(z ，yro)代入式 (1b)，并令 b = b2一 (A3zo+ A ‘。) (3) 若 bI一0，则可得z。是式 (1)的解．否则选取 Y的 —r个初值 Y， ，…，… 记y (1，Y2，… ， 一 ， )，满足 det(Y)：~0．由x：--A1AY得r×(—r)矩阵X，构造 —r阶方阵 A = A3X + A Y (4) 然后，求解以下 —r维方程组(称为剩余方程组)，即 A ： b (5) 下面将证明式 (5)与式(1)的可解是等价的． 定理 如果 det(Y)：~0，则有： (1)式 (5)有解 = 的充分必要条件是式 (1)有解 一 = xo+