一种新的基于Matlab环境的同伦路径跟踪算法62266.pdfVIP

第28卷第5期 东北大学学报(自然科学版) V01．28．No．5 2007年5月 ofNortheastern 2007 Journal University(NaturalScience) May 一种新的基于Matlab环境的 同伦路径跟踪算法 惠兴杰1，李 霞2，陈一呜2 (1．东北大学秦皇岛分校，河北秦皇岛066004；2．燕山大学理学院，河北秦皇岛066004) 摘 要：由于同伦延拓算法不需要初始值就能求出非线性方程组的解，但计算量大，为此利用Matlab平 台和它里面的函数“＼”，采用割线预估和牛顿校正，提出了一种新的路径跟踪算法——基于Matlab环境的同 伦路径跟踪算法．算例证明：该算法正确、简便易行、可靠性高，为同伦算法提供了简便的实现方法． 关键词：同伦延拓算法；“＼”；割线预估；Matlab；跟踪 中图分类号：TH113 文献标识码：A 文章编号：1005．3026(2007)05—0620—03 ANew BasedonMatlab HomotopicPath-TrackingAlgorithm HUl Xia2．CHEN Xing-jiel．Ll Yi—min92 of of (1．BranchQinhuangdao，NortheasternUniversity，Qinhuangdao066004，China；2．SchoolSciences， Yanshan University，Qinhuangdao066004，China．Correspondent：HUI cn) Abstract：Thecontinuationcansolvenon—linear withoutinitial homotopy algorithm equations lotsof shallbedone．Therefore．anew values,but computation homotopicpath—tracking onMatlab based anditsfunction“＼”is andNewton— algorithm proposedusingsecant—prediction showedthatthe is and correction．Examplesalgorithmcorrect，simple，convenientreliable，thus fl to the offeringwayimplementeasilyhomotopyalgorithm． continuation Keywords：homotopy 对于非线性方程组的求解问题，传统数值迭 Matlab和它里面的函数，提出了跟踪同伦路径的 代法需要选定合适的初始值，而区间分析法虽是 一个新算法．它的基础是割线预估、牛顿校正思 一种大范围收敛的算法，理论上可以求出全部解， 想．因为在文献[7一10]中用的欧拉预估，即切线 但仍需要确定合适的初值区间，同伦延拓算 预估，每步都需要矩阵求逆，而编程求一个矩阵的 法L