数值分析作业资料.docVIP

数 值 分 析 作 业 ——非线性方程的求解方法与分析 学院： 学号： 姓名： 摘 要 本文主要阐述了五种非线性方程的求解方法，分别为二分法、简易牛顿法、牛顿迭代法、牛顿下山法与弦截法。并分别对五种求解方法的计算结果进行了相应地分析。二分法运用函数有根区间中点与端点的函数值，缩小根区间，从而得到较快的收敛速度。牛顿迭代法，是一种常见的求解具有单重零点的非线性方程的数值方法，具有局部二阶收敛性。简易牛顿法便是简化的牛顿迭代法，将迭代点的导数值固定为初始值点的导数值，从而简化计算次数。牛顿下山法，为避免初值选取不当而使得迭代不收敛而在牛顿迭代法改进的方法。弦截法，克服了牛顿迭代法需求零点处函数导数的缺点，使用两次迭代点的差商替代了函数的导数值。本文非线性方程的求解方法均运用MATLAB编程及实现。 关键词：非线性方程；二分法；牛顿迭代法；牛顿下山法；弦截法 目 录 第一章 非线性方程 1 1.1 非线性方程简介 1 1.2 非线性方程求解方法简介 1 1.2.1 二分法 1 1.2.2 牛顿迭代法 2 1.2.3 牛顿下山法 4 1.2.4 简易牛顿法 4 1.2.5 弦截法 5 第二章 计算机配置 7 2.1 处理器 7 2.2 存储设备 7 2.3 显卡 8 2.4 显示屏 8 2.5 操作系统 8 第三章 算法的MATLAB实现及结果分析 9 3.1 二分法 9 3.2 牛顿迭代法 12 3.3 简易牛顿法 14 3.4 牛顿下山法 17 3.5 弦截法 20 结 论 24 第一章 非线性方程 1.1 非线性方程简介 1.2 非线性方程求解方法简介 1.2.1 二分法 1.2.2 牛顿迭代法 下述个定理分别讨论了牛顿法的收敛性质： 定理1：对于方程，设在上有二阶连续导数且满足下述条件： ； ，，对任意的； 选取，满足 则牛顿法产生的序列收敛于在内的唯一根 定理2：对于方程在上的根，且，则存在的一个邻域，使任意初值，牛顿迭代收敛于，且满足。 定理3：设是方程的根，在的某个开区间内连续且，则存在当时，由牛顿迭代法产生的序列是以不低于二阶的收敛速度收敛到 1.2.3 牛顿下山法 1.2.4 简易牛顿法 1.2.5 弦截法 的一种十分重要的插值方法。用牛顿迭代法求解非线性方程的根时，每一步迭代都要计算一次导数值。当函数较为复杂时，计算导数往往较为困难，并且，在计算机上，计算一次导数的近似值比计算函数的近似值要麻烦的多。因此，为了避免求解函数的导数，选择使用差商近似代替微商： 于是，牛顿迭代公式转化为： 下面研究弦截法的几何意义： 经过点及点亮点做函数的割线，其点斜式方程为：，其零点为。把X用表示即得到迭代格式。它又成为割线法，需要两个初始值，割线与X轴交点的横坐标就是新的近似值。如图所示： 下面两个定理为弦割法收敛定理： 定理1：设在其零点的邻域内有二阶连续导数，，则当时，由弦截法迭代公式产生的序列收敛于，且收敛的阶为1.618。 定理2：设函数在区间上二阶连续可导，且满足下述三点： （1）； （2）对任意的，所用的计算机型号为acer（宏碁）Aspire V5-471G 33214G50Mabb。下面对本款笔记本电脑的配置进行简单介绍。 2.1 处理器 CPU系列：英特尔 酷睿i3 代系列 CPU型号：Intel 酷睿i3 3217U CPU主频：1.8GHz 总线规格：DMI 5 GT/s 三级缓存：3MB 核心类型：Ivy Bridge 核心/线程数：双核心/四线程 制程工艺：22nm 指令集：AVX，64bit 功耗：17W 2.2 存储设备 内存容量：4GB 内存类型：DDR3 插槽数量：2xSO-DIMM 最大内存容量：32GB 硬盘容量：500GB 硬盘描述：5400转 光驱类型：DVD刻录机 设计类型：光驱内置 光驱描述：支持DVD SuperMulti双层刻录 2.3 显卡 显卡类型：入门级独立显卡 显卡芯片：NVIDIA GeForce GT 620M 显存容量：1GB 显存类型：DDR3 显存位宽：128bit 流处理器数量：96 DirectX：11 屏幕尺寸：14英寸 屏幕比例：16:9 屏幕分辨率：1366x768 背光技术：LED背光 2.5 操作系统 操作系统：Windows 7 旗舰版 第三章 算法的MATLAB实现及结果分析 本章主要阐述了利用数学软件MATLAB编写第一章所叙述的五种求解非线性方程的方法，并求解的零解。 3.1 二分法 二分法的matlab程序代码： function [k,c,err,yc]=erfen(f,a,b,emg1,emg2) %二分法求解非线性方程 %f为所需