一类非线性代数方程组的反问题.pdfVIP

维普资讯 大 庆 石 油 学 院 学 报 第32卷 第 4期 2008年 8月 JOURNALOFDAQINGPETROLEUM INSTITUTE Vo1．32 No．4 Aug． 2008 一 类非线性代数方程组的反 问题 王玉学，刘春兰，周少华，王银凤 (大庆石油学院 数学科学与技术学院，黑龙江 大庆 163318) 摘 要：给出了部分解已知条件下求非线性代数方程组非齐次项和系数项的方法．基于优化算法，给出了方便的计 算公式． 关 键 词：非线性代数方程组；反演；梯度公式 中圈分类号 ：O177 文献标识码 ：A 文章编号 ：i000—1891(2008)04—0101—03 用数值方法求解数学物理方程 的正问题时，均可化为求解一个线性或非线性代数方程组 ．其相应 的反问题的数值化求解就可归结为对线性或非线性代数方程组的右端项或系数矩阵的参数项进行反演． 文献[3，4]对线性代数方程组进行研究，给出了反演线性代数方程组的右端项和系数矩阵参数项的优化算 法．笔者研究了一类非线性代数方程组参数项和右端项反演的优化算法 ，并总假定正问题是适定的． 1 反演右端项 非线性代数方程组为 A () — b ， (1) 式中：系数矩阵A 中的每个元素 口 ()均为 的函数．若式 (1)中的部分解可实测 出，记为 z ，z点，…， z：，其 中{il，i2，…，is)c {1，2，…，n)．利用这些实测解可反演式(1)中的右端项．记 ()： ∑(z一z)． (2) 由(1)可知 ，是b的函数，故式(2)中J通过 成为b的复合 函数． 若 b=b。时，()取得极小值，则 b 就是反问题 的解．若采用直接法求解时 rainJ()，每迭代一步就 需求 n次至 2n次的正问题，而当n较大时，其计算量大的惊人．采用解析法求 rainJ()，这时需要知道 grad ()，通过分析可给 出只需解一次正问题就可求得任一点处 grad ()的公式． 引理 1 方程组 (1)的解 对b 的偏导数满足 C 一 (3) 式中：矩阵c=(f，c，…，c)，c= +口，口为矩阵A的第 列， 是第 个元素为 1的单位向量． 证明：方程组 (1)两端对 b求偏导得 OA+A象： (4) A OAA 易证 一 3x _ _ _ _ _ ，A象： ，代入式c4得 收稿 日期：2007—03—19；审稿人 ：孔令彬 ；编辑：关开澄 基金项 目：黑龙江省 自然科学基金项 目(AOl--O4) 作者简介：王玉学(1970一)，男，博士生 。副教授 ，主要从事微分方程数值解语最优化方法的研究 维普资讯