非线性方程的数值解法

本文格式为Word版，下载可任意编辑——非线性方程的数值解法

《计算方法》

期末论文

论文题目非线性方程的数值解法

学院专业班级姓名学号指导教师日期

目录

摘要第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",插值点序列",iix,yi=0,1,…,n",要计算的函数点x。step2.FORi=0,1,…,n",i制拉格朗日基函数序列问题的提出和研究目的和意义非线性方程的问题在工程实践中有好多用途",研究其数值解法是当前一个研究方向。目前已有相当一部分算法在广泛使用于工程实践中。非线性方程组和无约束最优化的数值解法",一直是数值优化领域中热门的研究课题。本文对传统的方法进行改进和提出新的算法",该算法不仅有重要的论价值，而且有很高的实用价值。例如在天体力学中，有如下Kepler",开普勒方程",x-t-",sinx=0,0<",<1,其中t表示时间",x表示弧度，行星运动的轨道x是t的函数。也就是说，对每个时刻it",上述方程有唯一解ix，运动轨道位置。

国内外相关研究综述随着科学技术的高速发展和计算机的广泛应用",求解形如F(x)=0的非线性方程组问题越来越多的被提出来了",其中F是的连续可微函数。例如非线性有限元问题、非线性断裂问题、弹塑性问题、电路问题、电子系统计算以及经济与非线性规划问题等都可转化为非线性方程组的求解问题。只要包含有未知函数及其导函数的非线性项的微分方程，无论是用差分方法或有限元方法，离散化后得到的方程组都是非线性方程组。与线性方程组相比,非线性方程组的求解问题无论在理论上还是在解法上都不如线性方程组成熟和有效.例如,非线性方程组是否有解,有多少解,理论上都没有很好的

解法,而对于非线性方程组,除了形式极为特别的小型方程组以外,直接解法几乎是不可能的.因而,我们主要考虑迭代解法.一般都是采用线性化的方法去构造各种形式的迭代系列.寻常都要探讨以下几个基本问题:第一个问题是,迭代点列的适定性问题,即要求迭代点列是有意义的.例如对于牛顿法,Jacobi矩阵必需是非奇异的.其次个问题,也是最基本的问题,生成的迭代点列的收敛性以及极限点是否为方程组的解.最终一个问题是,迭代点列的收敛速度问题.

早在七十年代以前,大量学者在理论上和数值解法上都对非线性方程组做了大量的研究.OrtegaRheinboldt系统的介绍了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",1p=g(kp)，可得到序列值{kp}。(3)牛顿——雷扶生法假如f(x)f‘(x)和f\在根p附近连续则可将它作为f(x)的特性,用于开发产生收敛到根p的序列{kp}的算法。而且这种算法产生序列{kp}的速度比二分法快。牛顿——雷扶生法依靠于f’(x)和f\(x)的连续性,是这类方法中已知的最有用和最好的方法之一。

(4)迭代法的收敛阶和收敛方法、割线法只计算f(x)不计算f’(x)",而且在单根上的收敛阶R1.618033989。割线法比牛顿法收敛速度慢一些??顿法的收敛阶为2。