100.0009765362300000.000970703123000-0.00501757812500-2.047021484375021 表2误差序列{x-n,{xrpm}和{xm-9m} Xa-Fn Xe-De 0006000m00 0000m000 0000m0mm 0.003000000000000 .003000000000000 0.003000000000000 0.001500000000000 0.00400000000000 0.007300000000000 0.000750000000000 0.005250000000000 0.015750000000000 0.000375000000000 0.005625000000000 0.0318750000m0000 0.0001875 0.00581250000 0.0639375000 0.0000375000000 0.005906250000000 0.1279%8750000001 0.000046875000000 0.005953125000000 0.255984375000003 0.0000m3437500000 0.005976562500000 0.511992187500005 0.00001171875000 0.0059%828125000 1.0239609375001 10 0.000005859375000 0.005994140625000 2.047998046875021 绘制:图】误差序列xn:图2误差序列{r-p:图3误差序列x-gm} 程序代码: %The propagation oferrors %Three ways to compute the sequence(1/2) The original erors are.6.0.0030.003 AP-;q=x x11;0)0.994p(1lp20.497:q11q(20.497 for=2:11 x(0-1/2*x(-1) =l2*i-1上 end for i=3:11 p0-32*p(11/2*p-2 q(0=5/2*g(i-1g(i2X end10 0.000976562500000 0.000970703125000 -0.005017578125000 -2.047021484375021 表 2 误差序列{xn-rn},{ xn-pn}和{ xn-qn} n xn-rn xn-pn xn-qn 0 0.0060000000 0.0000000000 0.0000000000 1 0.003000000000000 0.003000000000000 0.003000000000000 2 0.001500000000000 0.004500000000000 0.007500000000000 3 0.000750000000000 0.005250000000000 0.015750000000000 4 0.000375000000000 0.005625000000000 0.031875000000000 5 0.000187500000000 0.005812500000000 0.063937500000001 6 0.000093750000000 0.005906250000000 0.127968750000001 7 0.000046875000000 0.005953125000000 0.255984375000003 8 0.000023437500000 0.005976562500000 0.511992187500005 9 0.000011718750000 0.005988281250000 1.023996093750010 10 0.000005859375000 0.005994140625000 2.047998046875021 绘制:图 1 误差序列{xn-rn};图 2 误差序列{xn-pn};图 3 误差序列{xn-qn}。 程序代码: % The propagation of errors % Three ways to compute the sequence {1/2^n} % The original errors are 0.006,0.003,0.003 x=zeros(11,1);r=x;p=x;q=x; x(1)=1;r(1)=0.994;p(1)=1;p(2)=0.497;q(1)=1;q(2)=0.497; for i=2:11 x(i)=1/2*x(i-1); r(i)=1/2*r(i-1); end for i=3:11 p(i)=3/2*p(i-1)-1/2*p(i-2); q(i)=5/2*q(i-1)-q(i-2); end