数值计算(二分法简单迭代法Newton迭代法弦截法(割线法双点弦法)).docVIP

本科生实验报告 实验课程 数值计算方法 学院名称 信息科学与技术学院 专业名称 计算机科学与技术 学生姓名 学生学号 指导教师 实验地点 实验成绩 二〇 一六六 实验一 非线性方程求根 1.1问题描述 实验Newton迭代法、弦截法(割线法、双点弦法)，求 x5-3x3+x-1= 0 在区间 [-8,8]上的全部实根，误差限为10-6。 要求：讨论求解的全过程，对所用算法的局部收敛性，优缺点等作分析及比较， 第2章 算法思想 2.1二分法 思想：的单调有根区间内，将有根区间二分寻找的解 步骤： 1.取中点(x0+x1)/2 2.若f(mid)=0,则mid为方程的，否则比较与两端的符号若f(x0) 异号，则根[x0,mid]之间，否则在[mid,x1] 3并重复上述步骤，达到精度要求，则mid方程的解 2.2 简单迭代法 思想：迭代法是一种逼近的方法，是固定公式校正的，使之逐步精确最后得到精度要求的结果。 迭代公式x)，迭代公式必须是收敛。 2.x1,x1=f(x0). 3.判断|x1-x0|是否满足精度要求，如不满足则重复上述步骤。 4输出x1,方程的 2.3 Newton迭代法 思想：设r是的根，选取作为r的初始近似值，过点做曲线的切线L，L的方程为，求出L与x轴交点的横坐标，称x1为r的一次近似值。过点做曲线 的切线，并求该切线与x轴交点的横坐标，称为r的二次近似值。重复以上过程，得r的近似值序列，其中，称为r的 次近似值原函数的导数’(x);构造牛顿迭代公式 2.,若f’(x0)=0,退出计算，否则继续向下迭代。 3.若|x1-x0|满足精度要求，x1即为方程的近似解。 思想:为加速收敛，两个端点都在变动的弦，用牛顿迭代公式的导f’(x)。 步骤： 1.构造双点弦法公式 2x2=x1-f(x1)(x1-x0)/f(x1)-f(x0); 3.判断f(x2)是否满足精度要求若没有则继续，否则输x2.x2即为方程的近似解。 第3章 测试结果及分析 测试结果 函数图像 函数 Y=x5-3x3+x-1 二分法(表1-1，1-2，1-3) [-1.6,-1.3] k xk k xk k xk 0 -1.45 5 -1.50156 10 -1.50493 1 -1.525 6 -1.50391 11 -1.505 2 -1.4875 7 -1.50508 12 -1.50504 3 -1.50625 8 -1.50449 13 -1.50506 4 -1.49688 9 -1.50479 14 -1.50507 表1-1 区间[-1.2,-0.9] k xk k xk k xk 0 -1.05 5 -0.998437 10 -1.00005 1 -0.975 6 -1.00078 11 -0.999976 2 -1.0125 7 -0.999609 12 -1.00001 3 -0.99375 8 -1.0002 13 -0.999994 4 -1.00312 9 -0.999902 14 -1 表1-2 区间[1.5,1.8] k xk k xk k xk 0 1.65 7 1.69102 14 1.69029 1 1.725 8 1.69043 15 1.69029 2 1.6875 9 1.69014 16 1.69029 3 1.70625 10 1.69028 17 1.69028 4 1.69687 11 1.69036 18 1.69028 5 1.69219 12 1.69032 6 1.68984 13 1.6903 表1-3 简单迭代法 初值-1.5 k xk k xk k xk 1 -1.5 7 -1.50435 13 -1.50493 2 -1.50217 8 -1.50453 14 -1.50497 3 -1.50287 9 -1.50466 15 1.50499 4 -1.50341 10 -1.50476 16 -1.50501