|
即使笔误的不算,你的几个 相对误差值都在10%以上,是远不如连乘式的误差水平的。
10万以上附近的偶数,连乘式计算值的相对误差大多数在0.01~0.05之间 。
我计算一些连续的10万以上的偶数看看计算值的相对误差怎么样?
M= 100018 S(m)= 635 S1(m)= 627 Sp(m)≈ 630.4 δ(m)≈-.0073 K(m)= 1.0244
M= 100020 S(m)= 1602 S1(m)= 1585 Sp(m)≈ 1641 δ(m)≈ .0244 K(m)= 2.6667
M= 100022 S(m)= 674 S1(m)= 664 Sp(m)≈ 671.3 δ(m)≈-.0039 K(m)= 1.0909
M= 100024 S(m)= 599 S1(m)= 592 Sp(m)≈ 615.4 δ(m)≈ .0274 K(m)= 1
M= 100026 S(m)= 1232 S1(m)= 1218 Sp(m)≈ 1230.9 δ(m)≈-.0009 K(m)= 2
M= 100028 S(m)= 627 S1(m)= 618 Sp(m)≈ 656.5 δ(m)≈ .047 K(m)= 1.0667
M= 100030 S(m)= 972 S1(m)= 961 Sp(m)≈ 984.7 δ(m)≈ .0131 K(m)= 1.6
M= 100032 S(m)= 1212 S1(m)= 1194 Sp(m)≈ 1230.9 δ(m)≈ .0156 K(m)= 2
M= 100034 S(m)= 670 S1(m)= 661 Sp(m)≈ 683.9 δ(m)≈ .0207 K(m)= 1.1111
M= 100036 S(m)= 594 S1(m)= 587 Sp(m)≈ 624.8 δ(m)≈ .0518 K(m)= 1.0151
M= 100038 S(m)= 1191 S1(m)= 1177 Sp(m)≈ 1231 δ(m)≈ .0336 K(m)= 2
M= 100040 S(m)= 815 S1(m)= 807 Sp(m)≈ 856 δ(m)≈ .0503 K(m)= 1.3907
M= 100042 S(m)= 604 S1(m)= 598 Sp(m)≈ 615.5 δ(m)≈ .0191 K(m)= 1
M= 100044 S(m)= 1475 S1(m)= 1460 Sp(m)≈ 1477.3 δ(m)≈ .0016 K(m)= 2.4
M= 100046 S(m)= 614 S1(m)= 608 Sp(m)≈ 615.5 δ(m)≈ .0025 K(m)= 1
M= 100048 S(m)= 658 S1(m)= 652 Sp(m)≈ 690.7 δ(m)≈ .0497 K(m)= 1.1221
M= 100050 S(m)= 1724 S1(m)= 1705 Sp(m)≈ 1783.4 δ(m)≈ .0344 K(m)= 2.8971
我曾经作过的10万-11万之间偶数素对计算值的相对误差统计:
M=[ 100002 , 102000 ] r= 317 n= 1000 μ= .017 σχ= .016 δmin =-.0338 δmax = .0804
M=[ 102002 , 104000 ] r= 317 n= 1000 μ= .019 σχ= .017 δmin =-.0381 δmax = .0777
M=[ 104002 , 106000 ] r= 317 n= 1000 μ= .023 σχ= .016 δmin =-.0244 δmax = .0791
M=[ 106002 , 108000 ] r= 317 n= 1000 μ= .028 σχ= .016 δmin =-.0376 δmax = .0906
M=[ 108002 , 110000 ] r= 331 n= 1000 μ= .03 σχ= .016 δmin =-.0214 δmax = .0798
-------------------------------------------------------------------------------------------------------------------------------------
M=[ 100002 , 110000 ] r= 331 n= 5000 μ= .0233 σχ= .017 δmin=-.0381 δmax = .0906
可以看到,最大的相对误差值也小于10%;相对误差中间值 0.025左右。因此说你的方法的误差水平不如连乘式的相对误差水平,是有实际依据的。
这也就是为什么连乘式要达到高精度计算值的偶数对象,必须在大偶数区域,因为其时各个偶数的素对计算值的相对误差的波动性比较小,才能进行相对误差的修正。
|
|