2021年牛顿迭代法实验报告.docVIP

2021年牛顿迭代法实验报告 2021年牛顿迭代法实验报告 PAGE / NUMPAGES 2021年牛顿迭代法实验报告 材料科学与工程试验教学中心 试验汇报 课程名称 试验名称 牛顿迭代法方程求根 专业 班级 学号 姓名 试验日期 材料科学与工程试验教学中心试验汇报 一．试验目 经过对牛顿迭代法作编程练习与上机运算, 深入体会牛顿迭代法不一样特点 二．试验原理（方法） 给定初始值x0, ε为根许可误差, η为｜f(x) ｜许可误差, N为迭代次数许可值。 假如f′(x0)=0或迭代次数大于N, 则算法失败, 结束; 不然实施② 计算x1=x0-f(x0)/f′(x0)。 若｜x1-x0｜ε或｜f(x1) ｜η, 则输出x1, 程序结束; 不然实施④ 令x0=x1, 转向① 三．试验仪器设备及操作步骤 电脑一台 操作步骤: 打开电脑, 打开VB程序, 编辑程序 四．试验数据及处理（N-S图和程序） 开始 f(x)=x3+x2-3x-3 输入初始值X。, ε 是 是 f′(X。)=0 或迭代次数大于N 否 否 X?=X。- f(X。)/f′(X。) 否 否 X。= X? ∣X?-X。∣ X。= X? 或∣f(x?)∣η 是 是 输出X? 结束 Private Sub Command1_Click() Dim x0, x!, eps!, k x0 = Val(InputBox(x=0, ê?è?3??μ, 0.5)) eps = Val(InputBox(eps=, ê?è?epsμ??μ, 0.000005)) x = x0 - f(x0) / p(x0) k = 1 Print Do While Abs(x - x0) = eps Print k= k, x= x, f(x) f(x) x0 = x k = k + 1 x = x0 - f(x0) / p(x0) Loop End Sub Function f(ByVal x!) As Single f = x ^ 3 + x ^ 2 - 3 * x - 3 End Function Function p(ByVal x!) As Single p = 3 * x ^ 2 + 2 * x - 3 End Function 五．试验结果分析及结论（算例及结果） 六．思索题 成绩 指导老师 日期