实验指导书3_误差..docVIP

实验指导书_误差 一、误差的来源 在运用计算方法解决实际问题的过程中，会出现各种各样的误差，必须注重误差分析.否则，一种合理的计算也可能得出错误的结果. 例1 用差商求在处导数的近似值. （1）取和，用手工计算，取五位数字计算； （2）取，，=0.000 000 000 000 001和=0.000 000 000 000 000 1分别用MATLAB软件计算，取十五位数字计算； （3）比较以上的运算结果，说明是否越小则计算结果越准确. 解 根据导数定义，可以用差商求在处导数的近似值 . 从理论上讲，越小则计算结果越准确.如果用手工计算，取五位数字计算.当时，得 0.328 0， 与导数的精确值比较，这项计算还是可取的.但是当时，得 ， 算出的结果反而毫无价值. 如果应用MATLAB软件计算，取十五位数字计算，结果就完全不同了.在MATLAB工作窗口输入下面程序 format long g; a=3;h=0.1;y=log(a+h)-log(a);yx=y/h 运行后得 yx = 0.32789822822991 将此程序中改为0.000 1，运行后得 yx = 0.33332777790385 后者比前者好.再取h = 0.000 000 000 000 001，运行后得 yx = 0.44408920985006 不如前者好.取h = 0.000 000 000 000 000 1，运行后得 yx = 0 算出的结果反而毫无价值. 例2 计算的近似值. 解 泰勒级数 , 取，得 . （1.1） 这是一个无限过程，计算机无法求到精确值.只能在（1.1）取有限项时计算，再估计误差.如果取有限项 作为的值必然会有误差，根据泰勒余项定理可知其截断误差为 . 如果取（1.）， n=8; s=1;S =1; for k=1:n s=s*k; S=S+1/s, end s, S, R=3/(s*(n+1)) 或 S1=1+1+1/2+1/(1*2*3)+1/(1*2*3*4)+1/(1*2*3*4*5)+1/(1*2*3*4*5*6)+1/(1*2*3*4*5*6*7)+1/(1*2*3*4*5*6*7*8), R1=3/(1*2*3*4*5*6*7*8*9) 运行后结果 S = R = 2.71827876984127 8.267195767195768e-006 因为截断误差为 所以e的近似值2.718 28. 二、误差和有效数字 例1 取作为的四舍五入近似值时，求其绝对误差. 解 在MATLAB工作窗口输入程序 jueduwucha=exp(1)-2.71828 运行后输出结果为 jueduwucha = 1.828 459 045 505 326e-006 例2 计算d 的近似值，并确定其绝对误差和相对误差. 解 因为被积函数的原函数不是初等函数，故用泰勒级数求之. , （1.1） 这是一个无限过程，计算机无法求到精确值.可用（1.1）的前四项代替被积函数，得 d)d==. 根据泰勒余项定理和交错级数收敛性的判别定理，得到绝对误差 = WU， MATLAB命令窗口输入计算程序如下： syms x f=1-x^2/(1*2*3)+x^4/(1*2*3*4*5)-x^6/(1*2*3*4*5*6*7) y=int(f,x,0,pi/2),y1=double(y) y11=pi/2-(pi/2)^3/(3*3*2)+(pi/2)^5/(5*5*4*3*2)-(pi/2)^7/(7*7*6*5*4*3*2) infd=double(inf) ,inf=int(sin(x)/x,x,0,pi/2) WU =(pi/2)^9/(9*9*8*7*6*5*4*3*2), R =infd-y11 因为运行后输出结果为： 1.370 762 168 154 49，=1.370 744 664 189 38， 1.750 396 510 491 47e-005， WU= 1.782 679 830 970 664e-005. 所以，的绝对误差为，故d.的相对误差为 0.007 3%. 三、数值计算中应注意的问题 一个数值问题的解决，往往要经过多次运算.每一步运算都可能产生误差，在反复多次计算的过程中，必然产生误差的传播和积累.显然，当误差积累偏大时，会使计算结果失真.因此，在每一步计算中，都应该防止产生误差升级的现象. 例1 求数的近似值. 解 （1）直接用MATLAB命令 x=(7^15)*(sqrt(1+8^(-19))-1) 运行后输出结果 x = 0 问题出现在两个相近的数与相减时，计算