误差与范数概念附习题与答案.docx

PAGE PAGE 1 / 6 第一章 误差与范数 误差的来源 例 1.1.1 用差商  f (a) f (a h) h f (a) 求  f (x)  ln x 在 x  3 处导数的近似值 . 取 h 0.1 ， h 0.000 1 ， h =0.000 000 000 000 001 和 h =0.000 000 000 000 000 1 分别用 MATLAB软件计算，取十五位数字计算 . 解 在 MATLAB工作窗口输入下面程序 a=3;h=0.1;y=log(a+h)-log(a);yx=y/h 运行后得 yx = 0.32789822822991 将此程序中 h 改为 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 算出的结果反而毫无价值 . 例 1.1.2 分别求方程组 AX b 在下列情况时的解，其中 A 2 2 1 1 . 1 1.01 （1） b ; （ 2） b . 2 2.01 解 （1） 首先将方程组输入程序 AX b 化为同解方程 X A 1b ，然后在 MATLAB工作窗口 b=[2,2];A=[1,1;1,1.01]; X=A\b 2 2 运行后输出当 b 时， AX 2 2 b 的解为 X ； 0 1 （2）同理可得，当 b 时， AX 2.01 b 的解为 X . 1 例 1.1.3 计算 e的近似值 . 解 泰勒级数 x x 2 x 3 x 4 x n 取 x 1 ，得 e 1 x 2 ! 3! 4 ! n! ( x ) , e 1 1 1 1 1 2 ! 3! 4 ! 1 . （ 1.2 ） n! 这是一个无限过程，计算机无法求到精确值 .只能在（ 1.2 ）取有限项时计算，再估计误差 .如果取有限项 1 sn (1) 1 1 2 ! 1 1 1 3 ! 4 ! n! 作为 e的值必然会有误差，根据泰勒余项定理可知其截断误差为 e s n (1) e (0 ( n 1)! 1) . 如果取（ 1.2 ）的前九项，输入程序 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 s8  (1) e (8 1)! 3 8.267196 9! 10 - 6 ( 0  1), 所以 e 的近似值 e s8 (1) 1 1 1 1 2 ! 3! 1 1 4 ! 5 ! 1 1 6! 7 ! 1 2.718 28. 8 ! 1.2 误差和有效数字 例 1.2.1 取 2.718 28 作为 e的四舍五入近似值时，求其绝对误差和相对误差 . 解 在 MATLAB 工作窗口输入程序 juewu=exp(1)-2.71828 运行后输出结果为 juewu = 1.828 459 045 505 326e-006 例 1.2.2 计算 2 0 解 因为被积函数 sin x sin x x d x 的近似值，并确定其绝对误差和相对误差 . x 的原函数不是初等函数，故用泰勒级数求之 . sin x x 2 1 x 4 x６ x８ ( x ) , （ 1.5 ） x ３ ! ５ ! ７ ! ９ !  24x x x６ 2 4 这是一个无限过程，计算机无法求到精确值 .可用（ 1.5 ）的前四项 1 代 替被积函数  sin x x ，得 ３ ! ５ ! ７ ! 3 5 7 2 sin x 2 x 2 x 4 x６ ( ) ( ) ( ) 2 2 2 ? y d x ( 1 )d x = = y . 0 x 0 ３ ! ５ ! ７ ! 2 3 ３ ! 5 ５ ! 7 ７ ! 根据泰勒余项定理和交错级数收敛性的判别定理，得到绝对误差 R y y? ( ) 92 = WU， 9 在 MA TLAB 命令窗口输入计算程序如下： syms x 9 9! f=1-x^2/(1*2*3)+x^4/(1*