实验五 误差分析

《实验五 误差分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验五绪论－－误差分析【实验目的】1、了解数值计算中的误差种类，及避免误差危害的几种手段，2、深刻体会”数学上恒等，数值上不一定恒等”的含义3、为本课程的学习准备良好的数值思想【实验内容】1、误差的来源与分类2、数值计算中避免误差危害的若干方法3、数值实验举例4、根据要求，完成实验报告中的内容【实验指导】1)误差的来源与分类误差的来源是多方面的,通常误差主要由以下4个方面的因素引起：（1）模型误差（ModelingError）――把实际问题向数学问题转化的过程中，忽略了一些对问题影响不是很大的因素，我们称这种忽略了的因素为模型误差；（2）观测误差（MeasurementError

2、）――在一般的数学模型中，往往含有比较多的参数，而这些参数的值一般都需要通过观测得到，而观测得到的结果由于受到观测设备、观测方法等因素的影响往往都有误差，我们称这种由于观测引起的误差为观测误差。（3）截断误差（TruncationError）――当我们不能得到数学模型的精确解时，通常要用数值方法求它的近似解，其近似解与精确解之间的误差称为截断误差。例如：在计算机上直接使用公式计算时，会出现无穷过程的计算，不能在有限时间内得到需要的结果，因此，通常需要将上述无穷过程近似为有穷过程：，由此可以得到近似的计算结果，这样用数值方法中的有穷过程替代数学模型中的无限过程时，就会产生上述截断误

3、差。截断误差又称为方法误差。（4）舍入误差（RoundoffError）――由于计算机的字长有限，在使用计算机进行数据处理时，计算机表示的数据或计算结果会与原始数据或理论上的计算结果有差异，这种误差就是舍入误差。比如说，在计算机上表示时，只能表示成的形式，这里与的误差就是舍入误差。由于误差是不可避免的，我们只能尽可能的减少它对计算结果的影响。在计算方法这门课程中，我们主要关心如何减少截断误差与舍入误差对计算结果的影响。2)数值计算中避免误差危害的若干方法（1）选择稳定的算法――稳定算法是指后面一次运算产生的误差不能把前面一次运算所产生的误差扩大，这样就可以保证在运算过程中计算结果

4、的误差能够控制，保证计算结果的可靠性；（2）避免相近两数相减――由于两个相近数相减，会丢失有效数字，从而会增大相对误差。避免相近数相减需要结合具体的问题采取不同的方法或技巧。例如：当时，利用等式可以避免相近数相减。（2）避免分母过小或用绝对值较大的数作乘数――分母过小或乘数过大会导致计算机数据溢出，从而影响计算结果；（3）避免大数“吃”小数；（4）简化计算，减少计算步骤――一般来说运算次数减少，则计算过程中的积累误差有可能下降，这样就可以达到降低误差的目的。1)数值实验举例（1）下述三个表达式在数学上是恒等的，试在字长为的计算机上，分别利用这三个表达式计算其在区间上一些点处的函数

5、值，并比较计算结果，并说明理由。解：为了能够模拟字长为的计算机上的数值计算我们先编写函数digit（x,m），其功能是将向量x表示成字长为m的规格化浮点向量。函数digit（x,m）的Matlab程序如下：functiony=digit(x,m)%Thisfunctionisusedtoroundxtowards%anearestnormalizedscientificm-digitnumber.%Forexample,%digit(12.345,3)=0.123*10^2;%digit(12.345,4)=0.1235*10^2;%digit(0.012345,3)=0.123

6、*10^{-1}.%Input:%-xisavectorinR^{k}.%-misthegivennumberofsignificant%decimaldigitsofcomputer.%Sep.,26,2007byXuMinghua.k=max(size(x));y=x;%initializethevalueofy.fori=1:kifx(i)<0sign=-1;elsesign=1;endx(i)=abs(x(i));p=0;ifx(i)<0.1&x(i)>epswhilex(i)<0.1x(i)=x(i)*10;p=p-1;endendifx(i)>=1whilex(i)>

7、=1x(i)=x(i)/10;p=p+1;endendy(i)=round(x(i)*10^m)/10^m;y(i)=sign*y(i)*10^p;endreturn为了便于看出上述三种表达式在数值计算上的差异，我们设计算机的字长，并将分别取为：2e-12e-22e-32e-42e-5pi/2pi3*pi/22*pi利用下述主程序可以得到如下计算结果：xABC0.20000.50000.50000.50000.020000.50000.50000.002000.50000.500