西南交大数值分析上机非线性方程MATLAB求解.docxVIP

数值分析 摘要 非线性方程问题是经常出现的，求解此类方程往往很难得到精确解，经常需 要求其近似解。相应的，其求近似解的方法也逐渐得到大家的重视。本文通过对 《数值分析》的学习，对二分法、不动点迭代法、埃特金加速收敛方法、史蒂芬 森迭代法、牛顿法、简化牛顿法和弦截法共7种方法进行MATLAB编程，通过 求解算例来探究各种方法的适用条件和优缺点。 关键词：非线性方程；数值方法；MATLAB编程 Ll录 TOC \o 1-5 \h \z \o Current Document 第1章引言 3 \o Current Document 第2章 非线性方程的数值解法 4 二分法 4 1二分法简介 4 1.2二分法MATLAB主程序： 4 \o Current Document 不动点迭代法 6 1.不动点迭代法 6 2.不动点迭代法MATLAB主程序： 6 \o Current Document 埃特金加速收敛方法 7 1.埃特金加速收敛方法简介 7 2.埃特金方法MATLAB主程序： 7 \o Current Document 斯特芬森迭代法 8 1.斯特芬森迭代法简介 8 2.史蒂芬森迭代法MATLAB主程序： 8 \o Current Document 牛顿法 9 1.牛顿法简介 9 2.牛顿法的MATLAB主程序： 10 \o Current Document 简化牛顿法 11 1.简化牛顿法简介 11 2.简化牛顿法的MATLAB主程序： 11 \o Current Document 弦截法 12 1.弦截法简介 12 2.弦截法的MATLAB主程序： 12 \o Current Document 第3章 非线性方程求解算例 14 例1.求方程 14 例2.求方程 15 例3.求方程 16 例4.求方程 18 例5.求方程 19 例6.求方程 20 例7.求方程 22 例求方程 23 \o Current Document 第4章结果分析 25 \o Current Document 第5章结论 26 \o Current Document 参考文献 27 第1章引言 非线性是实际问题经常出现的，并且在科学与工程计算中的地位越来越重要, 很多我们熟悉的线性模型都是在一定条件下有非线性问题简化而来的，为了得到 更符合实际的解答，往往需要研究非线性模型，从而产生非线性科学，它是21世 纪科学技术发展的重要支柱。非线性问题有无限维的如微分方程，也有有限维的。 但要用计算机进行科学计算都要转化为非线性的单个方程或方程组求解。本文采 用7种方法即二分法、不动点迭代法、牛顿法、简易牛顿法、埃特金加速收敛方 法、弦截法、史蒂芬森法，进行MATLAB编程求解。然后对八个例子进行计算 和比较分析，从而发现各种算法的优势与不足，进而对各种算法有更深刻的理解。 本文的计算平台如表1-1所示： 表1-1计算平台 处理器（CPU） Illtel(R) Core(TM) i5-10400F CPU @ 2.9OGHZ 运行内存（RAM） 8.0OGB 操作系统 WindoWS 10教育版 软件版本 MATLAB R2018b 第2章非线性方程的数值解法 L二分法 1.1二分法简介 若f(Q是区间［a,b］上的连续函数，且f (α) f (h) 0,则/XQ在［a,b］内必有 一个零点。 因为f (Q)f(b) 0,所以函数fα)在区间［a,b］上改变符号，因此它在这个区 间内至少存在一个零点。这是中值定理的结论，因此利用区间减半的方法对这一 数值计算思想进行实践运算。 给定精确度ξ,用二分法求函数f(x)零点近似值的步骤如下： 1确定区间［a,b］,验证f(α) ?∕^(b) 0,给定精确度ξ. 2求区间(a, b)的中点c. 3计算f(c). (1)若f(c) = 0,则C就是函数的零点； ⑵若f(α)?f(c) V0,则令 Zc; ⑶ 若∕?(c)?∕?(b)V0,则令 a=c. (4)判断是否达到精确度g:即若Ia-Z)I t则得到零点近似值a (或b),否则重 复步骤2-4. 1.2二分法MATLAB主程序： EUICtiOn ［xp, eιτ, numj =eff (a, b, tol, max) %输入端点值、精度、最大迭代次数，返回根近似值、误差、迭代次数。 tic; %求解计时起点。 fa=fιιn (a) ; fb-fιιn (b); if fa?fbO disp (, eιτor,) ; %判断区间是否合法。 end k二 O; tt=100; WIIiIe (tt=tol∕2) (kmax) %^ :分法循环求解直到满足停机准则。 k=k+l; C=(a÷b)/2; fc=