非线性方程几种数值解法与MATLAB 程序.docVIP

. . 　　 本文档下载自 HYPERLINK / 文库下载网，内容可能不完整，您可以点击以下网址继续阅读或下载：  HYPERLINK /doc/1acb024c767f5acfa1c7cdda.html /doc/1acb024c767f5acfa1c7cdda.html 非线性方程几种数值解法的MATLAB程序 无 第１期２００８年１月 岳阳职业技术学院学报ＪＯＵＲＮＡＬＯＦＹＵＥＹＡＮＧＶＯＣＡＴＩＯＮＡＬＴＥＣＨＮＩＣＡＬ ＣＯＬＬＥＧＥ Ｎｏ．１Ｊａｎ．２００８ 非线性方程几种数值解法的ＭＡＴＬＡＢ程序 卢翼飞 （岳阳职业技术学院公共课部，湖南 摘 岳阳 ４１４０００） 要：研究非线性方程数值解，给出了二分法、简单迭代法、牛顿迭代法在ＭＡＴＬＡＢ的程序，并进 行了近似计算，可知牛顿迭代法收敛最快。 关键词：非线性方程；ＭＡＴＬＡＢ程序；二分???；迭代法中图分类号：Ｏ２４１．７ 文献标识码：Ａ 文章编号：１６７２－７３８Ｘ（２００８）０１－００９０－０４ １引言 在科学研究和科学计算中常常碰到非线性方程求解问题。非线性方程的解一般不能解析求出。所以数值解法显得非常重要，而数值解法在实际中的实现则更为重要，下面将介绍几种数值解法在Ｍａｔｌａｂ里的实现程序。 ２二分法 设ｆ（ｘ）在［ａ，ｂ］连续，假定ｆ（ａ）＜０，ｆ（ｂ）＞０，取中点ｘ０＝ａ＋ｂ，检查ｆ（ｘ０）符号。若ｆ（ｘ０）＝０，则ｘ０就是一个 根；若ｆ（ｘ０）＞０，记ａ为ａ１，ｘ０为ｂ１，则得有根区间［ａ１，ｂ１］；若ｆ（ｘ０）＜０，记ｘ０为ａ１，ｂ为ｂ１，则得有根区间［ａ１，ｂ１］。后两种情况都得到有根区间［ａ１，ｂ１］，它的长度为原区间的一半。对［ａ１，ｂ１］，令ｘ０＝ａ＋ｂ，再施以同样的方 法，可得新的有根区间［ａ２，ｂ２］，它的长度为［ａ１，ｂ１］的一半，如此反复进行下去，其中每一个区间是前一区间的一半。有 这就是方程的根。而ｘｎ＝ａ＋ｂ即为方程的近似根，且有估计误差 ｂ－ａ ｘｎ－ｘ≤２ ＊ （２．１）（２．２） 下面用二分法求ｆ（ｘ）＝ｘ３＋４ｘ２－１０ 在区间［１，２］上的根：因为二分法只能求单根，首先我们可以/doc/1acb024c767f5acfa1c7cdda.html搜索函数（２．２）在区间［１，２］的根的情况。在ＭＡＴＬＡＢ命令窗中输入：［２］ ｘ＝１：０．０１：２；ｙ＝ｘ．＾３＋４＊ｘ．＾２－１０；ｐｌｏｔ（ｘ，ｙ） 收稿日期：２００７－０８－１５ 作者简介：卢翼飞（１９６７－），女，湖南岳阳人，岳阳职业技术学院公共课部数学讲师。研究方向：经济数学。 总第６９期卢翼飞：非线性方程几种数值解法的ＭＡＴＬＡＢ程序 ９１ 图１ ｆ（ｘ）在区间［１，２］的图像从图１中可以看出式（２．２）在区间［１，２］间有惟一的一个大于１．３５而小于１．４的单根． ?ＭＡＴＬＡＢ实现先建立ｅｒｆｅｎ．ｍ ｆｕｎｃｔｉｏｎｙ＝ｅｒｆｅｎ（ｆｕｎ，ａ，ｂ，ｅｓｐ） ｉｆｎａｒｇｉｎ＜４ｅｓｐ＝１ｅ－４；ｅｎｄ ｉｆｆｅｖａｌ （ｆｕｎ，ａ）＊ｆｅｖａｌ （ｆｕｎ，ｂ）＜０ ｎ＝１； ｃ＝（ａ＋ｂ）／２；ｗｈｉｌｅｃ＞ｅｓｐ ｉｆｆｅｖａｌ（ｆｕｎ，ａ）＊ｆｅｖａｌ（ｆｕｎ，ｃ）＜０ｂ＝ｃ；ｃ＝（ａ＋ｂ）／２； ｅｌｓｅｉｆｆｅｖａｌ（ｆｕｎ，ｃ）＊ｆｅｖａｌ（ｆｕｎ，ｂ）＜０ａ＝ｃ；ｃ＝（ａ＋ｂ）／２； ｅｌｓｅｙ＝ｃ；ｅｓｐ＝１００００； ｅｎｄｎ＝ｎ＋１；ｅｎｄｆｏｒｍａｔｌｏｎｇ ｙ＝ｃ； ｅｌｓｅｉｆｆｅｖａｌ（ｆｕｎ，ａ）＝＝０ｙ＝ａ； ｅｌｓｅｉｆｆｅｖａｌ（ｆｕｎ，ｂ）＝＝０ ｙ＝ｂ； ｅｌｓｅｄｉｓｐ（＇ｍａｙｎｏｔｂｅａｒｏｏｔ＇）；ｅｎｄｎ 再建立ｆｃ．ｍ函数文件 ｆｕｎｃｔｉｏｎｙ＝ｆｃ（ｘ）；ｙ＝ｘ＾３＋４＊ｘ＾２－１０； 最后在ＭＡＴＬＡＢ命令窗中输入： ｅｒｆｅｎ（＇ｆｃ＇，１，２，０．０００１）ｎ＝５０ａｎｓ＝/doc/1acb024c767f5acfa1c7cdda.html １．３６５２３００１３４１４１０ 即迭代５０次计算结果为１．３６５２３００１３４１４１０。精度已经非常之高。 如果希望精度达到，只需要迭代１３次，结果是 １．３６５１７３３９０。 二分法优点是计算简单，收敛性有保证。缺点是收敛不快，不能求复根及重根。 ３简单迭代法 将函数ｆ（ｘ）＝０非线性方程变为同解方程：ｘ＝!（ｘ），其中函数!（ｘ）为连续函数。然后构造迭代格式 ｘｋ＋１＝!（ｘｋ） （３．１） 并选取初值ｘ０。根据初值ｘ０及迭代格式（３．１）可产生迭代序列｛ｘｋ｝．［３］ 下面用简单迭代法计算式（２．２）．由式（２．２）可以得到下列等价方程： ｘ＝!１（ｘ）＝ｘ－ｆ（ｘ）＝ｘ－ｘ３－４ｘ２＋１０１ｘ＝!２（ｘ）＝（１０－４ｘ） １ｘ＝!３（ｘ