《数值计算方法》试验报告册要点.doc

《数值计算方法》 实 验 报 告 册 姓名: 学号: 班级: 教师: 安徽农业大学理学院 应用数学系 学年 第 学期 目 录 目 录 i 实验报告范例 1 实验一 5 实验二 7 实验三 12 实验四 15 实验五 17 实验六 19 实验报告范例 实验名称 误差传播与算法稳定性 姓 名 学 号 班 级 指导教师 实验日期 成 绩 实验目的 理解数值计算稳定性的概念; 了解数值计算方法的必要性; 体会数值计算的收敛性与收敛速度。 实验题目 计算 实验原理 由，知 则 可得递推关系 1． ， 2． ， 下面分别以1,2递推关系求解： 方案一 ， 当时， 递推公式为 方案二 ， 当时 则 即 取递推初值 递推公式为 实验内容 将递推公式(1)中的初值，得 将递推公式(2)中的初值，得 实验结果 计算结果如下表所示 n n(1) I(2) 0 1 2 3 4 5 6 7 8 9 10 0.095310 0.046900 0.031000 0.023333 0.016667 0.033333 -0.166667 1.809524 0.095310 0.046898 0.031018 0.023153 0.018465 0.015353 0.013138 0.011481 0.010188 0.009232 0.008678 实验分析 由递推公式(1)知当时,应当为精确解,递推公式的每一步都没有误差的取舍,但计算结果,出现负值.由此看出,当较大时,用递推公式(1)中的近似是不准确的. 主要原因为初值不是精确值,设有误差,由递推公式(1)知 则有 误差随的增大而迅速增加,增加到得倍.由此可见,递推公式计算的误差不仅取决于初值的误差,公式的精确性,还有赖于误差的传递即计算的稳定性. 由递推公式(2)知,为估计值,并不精确,有,而由得 误差随递推公式逐步缩小.综上所述,在递推计算中,数值计算方法是非常重要的,误差估计、误差传播及递推计算的稳定性都会直接影响递推结果. 评阅意见 签名： 评阅日期： 附表 程序代码 说明： 具体实验题目与实验内容可自行根据实验指导书自行拟定； 报告填写用“宋体”（小四）格式字体； 实验报告完成后，以学生的“实验序号+姓名+学号”作为该word文件名保存，例如“张三”学号为，则本次实验报告的保存文件名为：“实验X张三.doc”； 在规定的时间内，学生将本报告通过电子邮件提交给授课教师，邮件的主题为：实验X张三。 算法编程语言可自选,程序代码可直接复制于实验报告附表八中,也可将可执行文件连同将实验报告压缩为rar格式文件一同提交。 实验一 实验名称 非线性方程求根 姓 名 学 号 班 级 指导教师 实验日期 成 绩 实验目的 理解非线性方程求根的基本算法； 掌握相应数值算法的程序编写； 探讨迭代法及初始值与迭代收敛性的关系。 实验题目 用迭代法求方程的根； 用Newton法求方程在附近的根； 求方程的全部实根，。 实验原理 实验内容 实验结果 实验分析 评阅意见 签名： 评阅日期： 附表 程序代码 问题2 (参考Matlab代码) function s=fun2(x) s=x^3-x-1; function y=dfun2(x) y=3*x^2-1; function [k,xk,yk,wucha]=newtonqx(x0,err,ddmax) x(1)=x0; for i=1: ddmax x(i+1)=x(i)-fun2(x(i))/dfun2(x(i)); wucha=abs(x(i+1)-x(i)); i=i+1; xk=x(i);yk=fun2(x(i)); if wuchaerr k=i-1; xk=x(i); return; end end if iddmax disp(请注意：迭代次数超过设定的最大值ddmax。) k=i-1; xk=x(i); end 实验二 实验名称 线性方程组的数值解法 姓 名 学 号 班 级 指导教师 实验日期 成 绩 实验目