基于凝聚函数法的非线性方程组求解.pdfVIP

第10卷第6期2010年2月 科学技术与工程 V01．10No．6 Science and ◎2010 1671—1815(2010)6—1494-04 TechnologyEngineering 基于凝聚函数法的非线性方程组求解 陈海霞杨铁贵 (郑州职业技术学院基础教育处，郑州450121) 摘要给出了求解非线性方程组问题的一种有效方法，称为凝聚函数法。首先把非线性方程组转化为一个不可微优化问 题，然后用一个称之为凝聚函数的光滑函数直接代替不可微的极大值函数，从而可把非线性方程组的求解转化为无约束优化 问题，因此可以直接利用现有的无约束优化算法软件求解。在此基础上，给出了相应算法，并做了数值实验，数值实验结果表 明了该算法具有收敛稳定，算法简单及计算效率高等优点。 关键词非线性方程组 凝聚函数 不可微优化 无约束优化 中图法分类号0241．7； 文献标志码A 求解非线性方程组是科学技术和工程应用中的 的求解算法更快。文献[3]探讨了用模拟退火算法 常见问题。非线性方程组的解法长期以来是数学和 求解非线性方程及方程组的方法，该方法把方程组 工程应用中的一个重要的研究内容。很多实际问题 求解问题转化为函数优化问题。文献[4]提出了一 的数学模型最终都归结为求解形如F(搿)=0的非种求解线性和非线性方程组的通用算法——蒙特卡 线性方程组，求解非线性方程组的常用方法是牛顿 罗算法。关于各种求解非线性方程组的智能算法正 法。牛顿法的收敛性和性能特征在很大程度上依赖 在不断的发展当中，随着计算机技术的发展，相信会 于初始点，在适当的条件下，该方法具有二阶收敛速 有更多更好的智能算法出现。现给出求解非线性方 度。牛顿法对迭代初始值的要求相当苛刻，如果初 程组的一种新方法，称为凝聚函数法。首先把非线 始值不在牛顿迭代的收敛域内，该方法失效。但是 性方程组问题转化为一个不可微优化问题，然后用 对于很多非线性方程组，选择好的初始点是一件非 一个称之为凝聚函数的光滑函数直接代替不可微的 常困难的事情；鉴于牛顿法的缺点，文献[1]结合遗 极大值函数，使用这一技术，可把非线性方程组的求 传算法和经典算法的优点，提出了一种用于求解非 解转化为无约束优化问题，因此可以直接利用现有 线性方程组的混合遗传算法。该混合算法充分发挥 的无约束优化算法软件求解，因此该算法特别易于 了遗传算法的群体搜索和全局收敛性，有效地克服 了经典算法的初始点敏感问题；同时在遗传算法中 数值实验结果表明了该算法的有效性。 引入经典算法(Powell法、拟牛顿迭代法)作局部搜 索，克服了遗传算法收敛速度慢和精度差的缺点。 1凝聚函数光滑化方法 文献[2]提出了基于微粒群优化算法的非线性方程 组求解，测试实例显示了该算法在求解非线性方程 组时具有简单性、高效性和普适性，比基于遗传算法 F(x)=∽(搿)以(搿)，…五(搿))1 可微的实值函数。若用0·忆表示向量空间p一范 2009年11月17日收到 数，于是非线性方程组可以转化为求解 第一作者简介：陈海霞(1970一)，女，河南省荥阳市人，讲师，研究方 万方数据 min 0 F(x)0。。 向：算法设计与研究。E—mail：y65558@163．tom。