数值分析上机实验报告一.doc

实验报告一 题目： 数值运算中误差分析的方法与原则 摘要：在我们的日常生活与学习中，很多具体问题抽象成数学模型都可以解决，而求解这些数学模型就要用到数值分析，本实验讨论的是数值分析中的误差。 前言：（目的和意义） 掌握误差来源，会对误差进行分析，了解简化计算步骤的基本原理和应用。 数学原理： 误差会随着计算步骤的增加而积累，计算步骤越多，误差越大。 为了减小数值计算结果的误差，应该尽量减少计算步骤，并对误差做好分析与处理。 程序设计一： （1）计算 并估计误差。 本实验采用Matlab的M文件编写，程序如下： I=1-exp(-1); n=input(请输入n的值); format long for N=1:n I=1-n*I; end I 当n=17时，I= -4.769577843020550e+020 程序设计二： （2）计算多项式( )并计算与的值 本实验采用Matlab的M文件编写，程序如下： n=input(请出入n的值); x=input(请出入x的值); a=3; p=3; for i=1:n a=2*a+3; b=x^i; p=p+a*b; end p 计算结果：= 1.099478611479765e+213 =600 结果分析和讨论： 计算时，要防止大数“吃”掉小数 要避免除数绝对值远远小于被除数绝对值； 要避免两相近数相减； 注意简化计算步骤，减少运算次数。同样一个计算问题，若能减少运算次数，不但可以节省计算时间，还能减小舍入误差。例如上述第二题，如果要直接计算的值再逐项相加，那么一共要做次乘法和n次加法。若采用以上程序中用到的算法那么只要用2n次乘法和n次加法就可以计算出的值。计算时还应当避免的就是在保留有效数字时竟可能多的保留例如上面第一题当n=17时I=-4+020，理论上所有的值应该是非负的，但是结果却出现了负值。因此在得出结果之后做好误差分析与处理也是非常重要的。 西安工程大学数值计算方法实验报告 .2.