精选-上机实验报告非线性方程的数值解法.docVIP

二分法实验 (1) 上机题目：二分法的应用 实验目的 ：熟悉二分法并在计算机上实现 实验要求 ： ①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告， 内容应该包括： 所用的算法语言， 算法步骤 陈述，变量说明，程序清单，输出计算结果，结果分析等等； ③用编好的程序在Ｍ atlab 环境中执行。 算法说明： ①找出 计算 f(x) 在有限根区间 [a, b]端点的值， f(a),f(b) ②计算 计算 f(x) 在区间中点（ a b ）处的值 f( a b ) . 2 2 ③判断 若 f( a b )=0 ，则 a b 即是根，计算过程结束，否则检验若 ( a b )f(a)0 ，则以 2 2 2 b 代替 b,否则以 a b 代替 a.反复执行步骤 ②和步骤 ③，直到区间 [a, b] 长度小于允许误 2 2 差 ，此时中点 a b 即为所求近似根。 2 计算例题 ：求 f(x)=x3- x-1 在[1,1.5] 的零点 . f(1)0,. f(1.5)0,delta= d=10 6 不动点迭代法实验⑵ 上机题目：不动点迭代法的实现 实验目的 ：熟悉迭代法并在计算机上实现 实验要求 ： ①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告， 内容应该包括： 所用的算法语言， 算法步骤 陈述，变量说明，程序清单，输出计算结果，结果分析等等； ③用编好的程序在Ｍ atlab 环境中执行。 算法说明： ①准备 将 f (x) 0 改写为等价的形式 x(x) 。若要求 x* 满足 f (x* ) 0 ，则 x* ( x* ) ；反之亦然。求 f ( x) 的零点，等价于求的不动点，选择一 个初始近似值 x0 代 入 x (x) 右端，得 x1 (x0 ) 。如此反复迭代计算 xk 1 ( xk ), k 0,1, ，此为迭 代公式。 ②不动点的存在性迭代法的收敛性 定理 1 设 (x) C [a,b] 满足： (1) 对任意 x [ a, b] 有 a ( x) b ; (2) 存在正常数 L1，使对任意 x, y [ a,b] 都有 | ( x) ( y) | L | x y | ，则 ( x) 在 [a,b] 上存在唯一不动点 x* 。 定理 2 设 (x) C [a,b] 满足定理 1 中的两个条件，则对任意 x0 [a, b] ，由迭代公 { x } * 式得到的迭代序列 收敛到 ( x) 的不动点 x ，并有误差估计 k | xk x* | 1LkL | x1 x0 | 。 计算例题 : 求 f(x)=x3 -x-1 在[1,1.5] 的零点 . 取 x0=1.5, d=10 6 牛顿法实验⑶ 上机题目 : 实现 Newton法 实验目的 : 编制求单变量非线性方程组的程序 . 实验要求 : ①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告，内容应该包括：所用的算法语言， 算法步骤陈述，变量说明，程序清单，输出计算结果，结果分析等； ③用编好的程序在Ｍ atlab 环境中执行简化 Newton法 . 来解方程。 计算步骤： ①准备 选定初始近似值 x 0 ，计算 f 0 =f(x 0 ), f 0 f (x0 ) 。 ② 迭 代 按 公 式 x 1 x0 f 0 f0 迭 代 一 次 ， 得 新 的 近 似 值 x 1 , 计 算 f 1 f ( x1 ), f1 f ( x1 ) . ③控制 如果 x 1 满足 1 或 f 1 2 ，则终止迭代，以 x1 作为所求的根；否则 x1 x0 ,当 x1 c 转步骤 4，此处 1， 2 是允许误差，而 x1 x0 时。 ,当 x1 c x1 其中 c 是取绝对误差或相对误差的控制常数，一般可取 c=1. ④修改 如果迭代次数达到预先指定的次数 N，或者 f =0，则方程失败； 否则以 （ x1 ， f1 , f 1 ）代替（ x0 , f 0 , f 0 ）转步骤 2 继续迭代。 例题：用牛顿方法解方程 x 3 -x-1=0 ，取迭代初值 x 0 =1.5, d=10 6 . 程序代码： ①二分法： function [x,k,error]=erfenfa(f,a,b,delta) 二分法求方程的根 % = 左边： [ 给定区间内的一个实根 % = 右边： erfenfa( 函数表达式 f=inline(f);  , 迭代次数 , 误差] , 区间左端点，区间右端点，预定精度  ) if f(a)*f(b)0 x=char( k=0; err