实验五误差分析

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

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

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

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

5、上，分别利用这三个表达式计算其在区间 【，2二上一些点处的函数值，并比较计算结果，并说明理由。2 2 2A = (1-cosx)/x ,B=(sinx)/x) /(1 cosx), C = 2(sin(x/2)/x)解：为了能够模拟字长为m的计算机上的数值计算我们先编写函数digit ( x, m)，其功能是将向量x表示成字长为 m的规格化浮点向量。函数digit( x, m)的Matlab程序如下: fun ctio n y=digit(x,m)% This fun cti on is used to round x towards% a n earest no rmalized scie

6、ntific m-digit nu mber.% For example,% digit(12.345,3)=0.123*1A2;% digit(12.345,4)=.1235*1A2;% digit(0.012345,3)=.123*1A-1.% In put:%- x is a vector in RAk.%- m is the give n nu mber of sig ni fica nt%decimal digits of computer.% Sep., 26, 2007 by Xu Min ghua. k=max(size(x);y=x; % in itialize the va

7、lue of y.for i=1:kif x(i)0sig n=-1;elsesig n=1;endx(i)=abs(x(i);p=;if x(i) eps while x(i) =1while x(i)=1 x(i)=x(i)/10;P=P+1；endendy(i)=rou nd(x(i)*10Am)/10Am; y(i)=sig n*y(i)*10y；endreturn为了便于看出上述三种表达式在数值计算上的差异，我们设计算机的字长 m = 2，并将x 分别取为：2e-12e-22e-32e-42e-5 pi/2 pi 3*pi/2 2*pi利用下述主程序可以得到如下计算结果：xABC0.

8、20000.50000.50000.50000.020000.50000.50000.002000.50000.50000.000200.50000.50000.000000.50000.50001.60000.40000.41000.41003.10000.21000.23000.20004.70000.04600.04400.05106.300000.00000.0002从上述计算结果可以知道：数学上恒等，数值上不一定恒等。结合防止误差的危害手段 思考造成上述差异的原因。上述计算的主程序如下：% Main program for Example 1% Show the differe ne

9、e of the followi ng 3 expressi ons% 1. A=(1-eosx)/xA2,% 2. B=(s in x/x)A2/(1+eosx),% 3. C=2(si n( x/2)/x)A2% in numerical eomputation.elem=2;x=digit(2e-1 2e-2 2e-3 2e-4 2e-5 pi/2 pi 3*pi/2 2*pi,m); a=(1-digit(cos(x),m)./digit(x.A2,m);a=digit(a,m); b=digit(digit(digit(si n( x),m)./x),m).A2,m)./digit(

10、1+cos(x),m);b=digit(b,m);c=2*digit(digit(si n(digit(x./2,m),m)./x),m).A2; c=digit(c,m);x, a, b, c(2)编写程序按照给定的顺序，在给定字长的计算机上分别计算S =12345 + 6 +6 + +61000S =，2 丄.丄%000 12345其中-i =1:100,为区间0, 0.5上的随机数，在此基础上演示数值计算过程中大数吃小 数的现象。解：下述程序先在许可的计算机精度下计算S和5，然后模拟字长为 m的计算机按照上述给定顺序计算S和$，为了能够看到大数吃小数的现象，字长取为m = 5。程序与结果

11、分别如下：Matlab 程序：% Compute $S 仁 A+B(1,1)+ cdots +B(1, n)$and $S2=B(1,1)+cdots+B(1, n)+A$, and% show big nu mber eat ing small nu mber.% Assume the computer is an m-digit computer.clc;m=5;% in itialize the dadaA=12345;n=1000;B=digit(0.5*digit(ra nd(1, n),m),m);% Compute $S 仁 A+B(1,1)+ cdots +B(1, n)$an

12、d $S2=B(1,1)+cdots+B(1, n)+A$,% in mach ine precisi on.S仁A;for i=1: nS1=S1+B(1,i);endS2=0;for i=1: nS2=S2+B(1,i);endS2=S2+A;fprintf(Numerical results computed in machine precision n);fprintf(S1= %8.4e,S2 = %8.4e , S1-S2= %8.4e n, S1, S2, S1-S2 );% Compute $S 仁 A+B(1,1)+ cdots +B(1, n)$and $S2=B(1,1)

13、+cdots+B(1, n)+A$,% in an m-digit nu mber computer.C1=A;for i=1: n6= digit(C1+B(1,i),m);endfprin tf(n);fprin tf(Numerical results computed in a computer with m-digit nu mber n)fprin tf(S1= %8.4e, The error of S1 is %8.4e n, C1, C1-S1);C2=0;for i=1: nC2=digit(C2+B(1,i),m);endC2=C2+A;fprin tf(S2= %8.4

14、e, The error of S2 is %8.4e n, C2, C2-S2);fprin tf(S1-S2= %8.4e n, C1-C2);运行上述程序可得结果如下：Numerical results computed in mach ine precisi onS仁 1.2596e+004, S2 = 1.2596e+004 , S1-S2= -9.0949e-012Numerical results computed in a computer with m-digit nu mberS仁 1.2345e+004, The error of S1 is -2.5088e+002S2= 1.2596e+004, The error of S2 is -1.0320e-001S1-S2= -2.5078e+002从上述结果可以看出如果利用现有的机器精度来计算S和S,则两者的差异不大，但如果在字长为5的计算机上计算s和s，