东南大学_数值_第二章_牛顿迭代法.doc

第二章 非线性方程的解法 ——牛顿迭代法 ****（学号） ****（姓名） 算法与程序题目见教材P56 上机题目20。 一、算法原理 根据题目的要求，是关于用牛顿迭代法法求解方程的通用算法。该法是一种通过斜率迭代的算法，其速度比二分法和简单迭代法都要快。其简单原理如下： 设则存在一个数使得由如下定义的迭代序列收敛到： 对于函数，则其递推规则是 定义序列，则序列收敛到，也可表示为。现简要证明： 对于，写出牛顿迭代公式 该公式可化简为 二、流程图 题目要求于用牛顿迭代法法求解方程的通用算法。其计算过程主要用到迭代，图流程图1所示。 输入各参数 k=1 迭代 计算各误差 误差在允许范围之内 T F break k=k+1 kmaxN 输出 图1 题1关于牛顿迭代法的算法流程图 三、计算代码 核心代码 p1=……; if (errdelta)|(relerrdelta)|(abs(y)epsilon),break; if kmaxN ,go to 1 完整代码 程序1：Newton.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Description : 牛顿迭代法 % Author : panyunqiang % Versoin : 1.0 % Date : 2012-9-21 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [p0, err, k, y]=Newton(p0, delta, epsilon, maxN) % input - p0 is the initial approximation to a zero of f % - delta is the tolerance for p0 % - epsilon is the tolerance for the function values y % - maxN is the maxium number of iterations %output - p0 is the Newton approximation to a zero % - err is the error estimate for p0 % - k is the number of iterations % - y is the function value f(p0) for k=1:maxN %%递归 p1=2*p0^3/(3*p0^2-3); %%计算误差 err=abs(p1-p0); relerr=2*err/(abs(p1)+delta); p0=p1; %%当前求出的根的函数值 y=p0^3/3-p0; %%判断 if (errdelta)|(relerrdelta)|(abs(y)epsilon) break; end end 程序2：Newton_Step.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Description : 寻找题目中关于牛顿迭代法收敛的尽可能大的delta % 搜索步进为step=10^(-6),即精确到小数点后六位 % Author : panyunqiang % Versoin : 1.0 % Date : 2012-9-21 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% format long step=10^(-6); delta=10^-8;epsilon=10^-8;maxN=1000; ps=0.6; [p0, err, k, y]=Newton(ps, delta, epsilon, maxN); while((abs(p0)=epsilon)(p0~=NaN)) ps=ps+step; [p0, err, k, y]=Newton(ps, delta, epsilon, maxN); end ps-step 四、计算结果及分析 运行程序Newton_Step.m，获得Newto