极小值搜索与平方和转化求解多元非线性方程组.pdfVIP

№ ．3 陕西科技大学学报 Jun．2o10 V0I．28 JOURNALOFSHAANXIUNIVERSITYOFSCIENCE TECHNOL0GY · 143 · 文章编号 ：1000—5811(2010)03—0043—05 极小值搜索与平方和转化求解多元非线性方程组 李建文 ，张成现 ，李 颀 ，马学宗 (1．陕西科技大学 电气与信息工程学院，陕西 西安 71002l；2．西安工程大学计算机科学学院，陕西 西安 710048) 摘 要 ：论述 了一种通过极值搜索法求解多元非线性方程组的通用算法．该算法主要通过移动 中心点与缩放域半径的方法搜索多元函数 函数极值 ，其所需要 的条件 十分宽松 ，无需求解多 元非线性方程组以及其它过程，因此推导过程简单 ，且收敛快，易于通过计算机编程 实现．此 外，利用平方和将 方程组求根 问题转化为多元非线性函数的极小值搜索，解决了多元非线性方 程组的求根 问题．经过推广 ，又解决了任意隐式微分方程 的数值计算 问题．由于这种算法通用 性极强，应该将它作为数值方法教科书 中的基本方法． 关键词 ：极小值搜 索；多元非线性方程组；平方和 中图分类号 ：TP301．6 文献标识码 ：A 0 引言 多元非线性函数的极值数值解问题是一个非常普遍的数学问题，在理论上往往把多元非线性极值问 题归结为多元非线性方程组的求根问题 [1。]，但实际上多元非线性方程组的求根问题 比多元非线性函数极 值问题更难解决．在计算数学中为了解决多元极值问题引出了许多新的方法 ，如最速下降法、遗传算法、边 缘检测算法等 4【]． 现在 ，由于微型计算机的普及和微型计算机软件 的快速发展，人们普遍使用 Matlab来解决数值计算 问题_7 ．但 Matlab的求根能力很有限，它不但速度慢 ，而且还对许多简单的方程组无法求解． 作者在解决石英晶振体热敏网络温度补偿问题时总结出了一种新的用于计算机查找极值的方法 “多 元非线性函数极值的微分搜索算法”，以后又将这种方法用在其它许多领域 ，如语音基频分析 、彩色电视机 荫罩曲面求解等 ．通过实际使用证明，这种方法比教科书里所提到的所有查找极值的方法都方便 ，而 且还可以将多元非线性函数极值 、多元非线性方程组、任意隐式微分方程等一系列的数值计算问题一并解 决．在此 ，作者将这种新方法与传统的牛顿迭代法进行比较． 1 传统迭代法及其缺点 设多元非线性 函数为 f(x ，z，…， )，根据微分原理 ，在多元非线性 函数的极值点上各个分量 的微 分值必然为 0，所以通过解如下的非线性方程组 一 0 (忌一0，1，2，… ，n) (1) * 收稿 日期：2010—03—10 作者简介：李建文(1959一)，男，陕西省西安市人，教授，研究方 向：计算机及其应用软件 的研究与开发 基金项 目：国家 自然科学基金项 目“声一电转换皮肤听声的进一步研究” · 144 · 陕西科技大学学报 第 28卷 必然获得到极值点的坐标．传统上都使用牛顿迭代法求解这个非线性方程组 ]． 求解过程可 以划分为 ：求偏导 、列方程组 、非线性方程组线性化、得出含 (+1)× 矩阵的 阶线性方 程组 、选迭代初值 、迭代口]，这种迭代的运算量级为 ，z，因此其过程十分复杂 ，尤其是 当一个非线性多元 函数的表达式很复杂时它的偏导函数的表达式更加复杂，对于一个非常复杂的非线性导函数方程组线性 化需要忽略许多相关因素，这必然造成结果的不准确性 ，之后的迭代初值如果选取不合适会造成迭代过程 发散 ，等等 ，并且其中每一个过程的实现都是极其 困难的． 实际上这种求极值的方法把问题反而复杂化 了，在实际中很难被使用，很多情况是不能使用．后面还 有实例说明． 另一种常用的算法称为最速下降法．最速下降法是根据多元函数的梯度的反方向进行搜索，仍然需要 求偏导，还需要选择一个步长．这种方法的