(数值分析实验一误差的影响.docVIP

内江师范学院 数 值 分 析 实 验 报 告 册 编制 张莉 审定 牟廉明 专业： 班级： 级 班 学号： 姓名： 数学与信息科学学院 2013年9月 说 明 一、学生在做实验之前必须要准备实验，主要包括预习与本次实验相关的理论知识，熟练与本次实验相关的软件操作，收集整理相关的实验参考资料，要求学生在做实验时能带上充足的参考资料若准备不充分，则学生不得参加本次实验，不得书写实验报告； 要求学生要认真做实验，主要是指不得迟到、早退和旷课，在做实验过程中要严格遵守实验室规章制度，认真完成实验内容，极积主动地向实验教师提问等； 学生要认真工整地书写实验报告，实验报告的内容要紧扣实验的要求和目的，不得抄袭他人的实验报告； 实验名称： 实验一 误差的影响 指导教师： 吴开腾 张莉 实验时数： 4 实验设备：安装了Matlab、C++、VF软件的计算机 实验日期：2013年 9月 25日 实验地点： 第五教学楼北902 实验目的： 掌握舍入误差的概念，理解数值稳定性。 实验准备： 在开始本实验之前，请回顾教科书的相关内容； 需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。 实验内容及要求 A题的两个实根，，试编写当时求两个实根的算法程序，并用该程序计算二次方程的两个根，使它至少具有4位有效数字。 B题 舍入误差在数值计算中是一个很重要的概念，在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果，因此，选取数值稳定性的算法，在数值计算中是十分重要的。 对计算定积分 分别采用下面两个递推公式进行计算，并比较实验结果分析出哪个算法是稳定，并给出具体原因。 递推公式（1）； 递推公式（2）。 C题 序列满足递推关系，若取作近似计算，问计算到时误差有多大？这个计算过程稳定吗？ D题 用下列两种方法计算的近似值，问哪种方法能提供较好的近似值？ （1）； （2）。 说明：实验过程应包括对问题的简要分析、求解方法、求解步骤、程序及其必要的图表等内容。 实验过程： A题实验过程 问题分析 本文是针对一元二次方程的解的问题，在本问中的是根据判断方程的关系来进行求解。满足时继续运用一元二次方程的求根公式时会出现大数吃小数的情况，这样就会有一个解就是错误的，但是其中有一个数是正确的。这样在解决此问题时就需要找到那个正确的解，然后运用维达定理求得另一个解。这样就可以避免大数吃小数的情况。 问题求解方法 此问题中的求解方法，就是找到一个正确的解，然后运用维达定理得到另一个解。 寻找正确的解的方法：根据系数的正负可以判断出用求根公式得到的正确解。 、，正确的解为 、，正确的解为 在用以上方法求得正确的解后，运用维达定理求解得另一个解。 求解流程和代码 求解流程示意图如下： 结果及解释 根据以上的程序输入：[1 -56 1]，可以得到其误差为0.173125402902485*e-014，可以得到其精度在有效数字为4位的范围内，则可知此算法是有效的。其所给的方程的解用此算法得到的结果为0.0178628407335566、55.9821371592664。 B题的实验过程 舍入误差在数值计算中是一个很重要的概念，在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果，因此，选取数值稳定性的算法，在数值计算中是十分重要的。 对定积分分别递推公式递推公式进行计算，比较实验结果分析出哪个算法是稳定，给出原因。 、分析 、（其中i=1、2表示利用递推公式一、二的算法），计算出利用两个算法的相对误差值。 最后，画出当n分别取1至20的整数时，利用递推公式一、二得到的结果与定积分理论值的差值的散点图。 四、实验求解结果 利用Matlab10.0软件编程计算。于是得到以下结果： 表1 n的取值 利用积分得到理论值 利用公式一得到的实际值 利用公式二得到的实际值 公式一与积分的 公式二与积分的 1 0.367879 0.367879 0.367879 0 0 2 0.264241 0.264241 0.264241 0 0 3 0.207277 0.207277 0.207277 8.32667E-17 2.77556E-