实验五误差分析.docVIP

实验 绪论－－误差分析 【实验目的】 1、 2、”数学上恒等，数值上不一定恒等”的含义 3、 【实验内容】 1、 2、 3、、根据要求，完成实验报告中的内容 【实验指导】 误差的来源是多方面的, 通常误差主要由以下4个方面的因素引起： 模型误差（Modeling Error）――把实际问题向数学问题转化的过程中，忽略了一些对问题影响不是很大的因素，我们称这种忽略了的因素为模型误差； 观测误差（Measurement Error）――在一般的数学模型中，往往含有比较多的参数，而这些参数的值一般都需要通过观测得到， 截断误差（Truncation Error）――当我们不能得到数学模型的精确解时，，计算时，会出现无穷过程的计算，不能在有限时间内得到需要的结果，因此，通常需要将上述无穷过程近似为有穷过程：，由此可以得到近似的计算结果，这样用数值方法中的有穷过程替代数学模型中的无限过程时，Roundoff Error）――由于计算机的字长有限，在使用计算机进行数据处理时，计算机表示的数据或计算结果会与原始数据或理论上的计算结果有差异，这种误差就是舍入误差。比如说，在计算机上表示时，只能表示成的形式，这里与的误差就是舍入误差。 由于误差是不可避免的，我们只能尽可能的减少它对计算结果的影响。在计算方法这门课程中，我们主要关心如何减少截断误差与舍入误差对计算结果的影响。 数值计算中避免误差危害的若干方法 选择稳定的算法――稳定算法是指后面一次运算产生的误差不能把前面一次运算所产生的误差扩大，这样就可以保证在运算过程中计算结果的误差能够控制，保证计算结果的可靠性； 避免相近两数相减――由于两个相近数相减，会丢失有效数字，从而会增大相对误差。避免相近数相减需要结合具体的问题采取不同的方法或技巧。例如：当时，利用等式 可以避免相近数相减。 避免分母过小或用绝对值较大的数作乘数――分母过小或乘数过大会导致计算机数据溢出，从而影响计算结果； 避免大数“吃”小数； 简化计算，减少计算步骤――一般来说运算次数减少，则计算过程中的积累误差有可能下降，这样就可以达到降低误差的目的。 数值实验举例 （1）下述三个表达式在数学上是恒等的，试在字长为的计算机上，分别利用这三个表达式计算其在区间上一些点处的函数值，并比较计算结果，并说明理由。 解：为了能够模拟字长为的计算机上的数值计算我们先编写函数digit（x, m），其功能是将向量x表示成字长为m的规格化浮点向量。函数digit（x, m）的Matlab程序如下： function y=digit(x,m) % This function is used to round x towards % a nearest normalized scientific m-digit number. % For example, % digit(12.345,3)=0.123*10^2; % digit(12.345,4)=0.1235*10^2; % digit(0.012345,3)=0.123*10^{-1}. % Input: % - x is a vector in R^{k}. % - m is the given number of significant % decimal digits of computer. % Sep., 26, 2007 by Xu Minghua. k=max(size(x)); y=x; % initialize the value of y. for i=1:k if x(i)0 sign=-1; else sign=1; end x(i)=abs(x(i)); p=0; if x(i) 0.1 x(i) eps while x(i) 0.1 x(i)=x(i)*10; p=p-1; end end if x(i) =1 while x(i)=1 x(i)=x(i)/10; p=p+1; end end y(i)=round(x(i)*10^m)/10^m; y(i)=sign*y(i)*10^p; end return 为了便于看出上述三种表达式在数值计算上的差异，我们设计算机的字长，并将分别取为：2e-1 2e-2 2e-3 2e-4 2e-5 pi/2 pi 3*pi/2 2*pi 利用下述主程序可以得到如下计算结果：