《 计 算 方 法 》
期 末 论 文
论文题目 非线性方程的数值解法
学 院 专 业 班 级 姓 名 学 号 指 导 教 师 日 期
目 录
摘要 第1 章 绪论
1.1 问题的提出和研究目的和意义 1.2 国内外相关研究综述 1.3 论文的结构与研究方法 第2 章 非线性方程的数值解法 2.1 二分法 2.2 迭代法
2.3 迭代法的局部收敛性及收敛的阶 2.4 牛顿迭代法 2.5 牛顿法的改进 2.6 插值
摘要",
数值计算方法，是一种研究解决数学问题的数值近似解方法，它的计算对象是那些。
在理论上有解而又无法用手工计算的数学问题。在科学研究和工程技术中都要用到各种计算方法。例如",在地质勘探、汽车制造、桥梁设计、天气预报和汉字设计中都有计算方法的踪影。本文讨论了非线性方程的数值解法：非线性方程的二分法、迭代法原理、牛顿迭代法，迭代法的收敛性条件及适合非线性方程的插值法等等。
第1 章 绪论
可以证明插值多项式L (x) n 存在并唯一。拉格朗日插值多项式的算法",step1.输入",插值节点控制数n",插值点序列", i i x , y i=0,1,…,n",要计算的函数点x。step2. FOR i =0,1,…,n ",i 制拉格朗日基函数序列问题的提出和研究目的和意义非线性方程的问题在工程实践中有很多用途",研究其数值解法是当前一个研究方向。目前已有相当一部分算法在广泛使用于工程实践中。非线性方程组和无约束最优化的数值解法",一直是数值优化领域中热门的研究课题。本文对传统的方法进行改进和提出新的算法",该算法不仅有重要的论价值，而且有很高的实用价值。例如在天体力学中，有如下Kepler",开普勒方程",x-t- ",sin x=0,0< ",<1,其中t 表示时间",x 表示弧度，行星运动的轨道x 是t 的函数。也就是说，对每个时刻i t ",上述方程有唯一解i x ，运动轨道位置。
国内外相关研究综述随着科学技术的高速发展和计算机的广泛应用",求解形如F(x)=0 的非线性方程组问题越来越多的被提出来了",其中F 是的连续可微函数。例如非线性有限元问题、非线性断裂问题、弹塑性问题、电路问题、电子系统计算以及经济与非线性规划问题等都可转化为非线性方程组的求解问题。只要包含有未知函数及其导函数的非线性项的微分方程，无论是用差分方法或有限元方法，离散化后得到的方程组都是非线性方程组。与线性方程组相比,非线性方程组的求解问题无论在理论上还是在解法上都不如线性方程组成熟和有效.例如,非线性方程组是否有解,有多少解,理论上都没有很好的
解法,而对于非线性方程组,除了形式极为特殊的小型方程组以外,直接解法几乎是不可能的.因而,我们主要考虑迭代解法.一般都是采用线性化的方法去构造各种形式的迭代系列.通常都要讨论以下几个基本问题:第一个问题是,迭代点列的适定性问题,即要求迭代点列是有意义的.例如对于牛顿法,Jacobi 矩阵必须是非奇异的.第二个问题,也是最基本的问题,生成的迭代点列的收敛性以及极限点是否为方程组的解.最后一个问题是,迭代点列的收敛速度问题.
早在七十年代以前,许多学者在理论上和数值解法上都对非线性方程组做了大量的研究.Ortega Rheinboldt 系统的介绍了n 阶非线性方程组的基本理论成果,并对牛顿法,延拓法等几种主要迭代法作了详尽的分析.另外,也有一些学者把非线性方程组的求解问题转化为极小化问题, 得到一类称为极小化方法的迭代法, 如下降法, 共轭方向法,Gauss-Newton 法等,李,莫&祁详细介绍了一些适合在计算机上求解的有效算法,如Broyden 算法,以及近十几年来发展的新方法,如区间迭代法,单调迭代法和单纯形法等. 论文的结构与研究方法
1.欲解决的主要问题是：综合当前各类非线性方程的数值解法，通过比较分析，二分法，迭代法，牛顿——雷扶生方法，迭代法的收敛阶和加速收敛方法，解非线性方程的插值方法，这以上五种的算法应用对某个具体实际问题选择相应的数值解法。
2.比较各类数值算法",分析其优缺点",并应用到具体的实际问题中。
3．利用计算机MATLAB 语言对非线性方程的数值解法进行程序设计。
研究的基本思路是结合目标所提出的问题针对各种方法来具体分析比较",
(1) 二分法 起始区间[a,b]必须满足f(a)与f(b)符号相反的条件。二分法的第一部是选择中点c=(a+b)/2，然后分析可能存在的三种情况如果f(a)和f(c)符号相反，则在区间[a,c]内存在零点。如果f(c)和f(b)符号相反",则在区间[c,b]内存在零点。如果f(c)=0，则c是零点。 （2）迭代法 迭代是指重复执行一个计算过程，直到找到答案。首先需要有一个用于逐项计算的规划或函数g(x)，并且有一个起始po。然后通过迭代规则k ",1 p =g( k p )，可得到序列值{ k p }。 (3)牛顿——雷扶生法 如果f(x) f ‘(x)和f \在根p 附近连续则可将它作为f(x)的特性,用于开发产生收敛到根p 的序列{ k p }的算法。而且这种算法产生序列{ k p }的速度比二分法快。牛顿——雷扶生法依赖于f’(x)和f \(x)的连续性,是这类方法中已知的最有用和最好的方法之一。
(4)迭代法的收敛阶和收敛方法、割线法只计算f(x)不计算f ’(x)",而且在单根上的收敛阶R 1.618033989。割线法比牛顿法收敛速度慢一些??顿法的收敛阶为2。当p 是一个M 阶根时??要更好的求根技术以获得比线性收敛更快的速度。最终结果显示??过对牛顿法进行改进??使其在重根的情况下的收敛阶为2。加速收敛方法有 Aitken 加速法和Steffensen 加速法。Steffensen 算法是促使迭代