|
本帖最后由 愚工688 于 2023-2-21 00:41 编辑
vfbpgyfk
修正的依据是什么?根据什么确定修正值?有修正公式吗?
------------------------------------------------------------------------------------------------------------
这是依据相对误差的定义推导出来的。
相对误差 Δ=计算值/真值-1;变形一下,1+ Δ=计算值/真值;
于是 真值= 计算值/(1+ Δ),
这就是真值与计算值之间的关系式。
如果已知一个偶数的计算值与相对误差,无疑就能够得到真值。(当然这样是无意义的,因为要得到相对误差,就必须知道真值)
但是在一个区间内,如果我们知道了其中的所有偶数的素数对的平均相对误差δ,那么我们用这个平均相对误差δ来替代每个偶数的真实的相对误差 Δ,尽管两者之间有一定的小差异,但是我们可以肯定通过这样的修正,能够使得计算值的精度得到很大的提升;
由于通过一些样本区间的偶数的素数对连乘式计算值的相对误差是呈现缓慢趋于0.20附近的,
因此可以用一个小样本区域的相对误差平均值δ来替代更大范围内的偶数的真实的相对误差 Δ,以提高素数对数量的计算精度:
我们还可以用比较大些的偶数区域的相对误差平均值δ来代替相对比较小一些的偶数区域的各个偶数的真实的相对误差 Δ,以此得到略微小于素对真值的具有比较高精度的下界素数对计算值。
这就是我从相对误差定义推理出来的高精度素数对的计算式子与比较高精度的素数对下界计算值的数理方法。
部分区域的偶数素数对计算值的相对误差统计计算数据:
M= 99999950 - 100000046 : n= 49 μ= .119 σχ= .001 δ(min)= .117 δ(max)= .122
1亿-100亿的样本的统计计算数据:
(标准偏差的通用符号为σχx ,μ-样本平均值)
100000000 - 100000098 : n=50 μ= .1192 σχ= .0013 δ(min)= .1156 δ(max)= .1224
1000000000 - 1000000098 : n= 50 μ= .1368 σχ= .0004 δ(min)= .1356 δ(max)= .138
2000000000 - 2000000098 : n= 50 μ= .1406 σχ= .0003 δ(min)= .1399 δ(max)= .141
3000000000 - 3000000098 : n= 50 μ= .1431 σχ= .0002 δ(min)= .1425 δ(max)= .1435
4000000000 - 4000000098 : n= 50 μ= .1449 σχ= .0003 δ(min)= .1441 δ(max)= .1456
5000000000 - 5000000098 : n= 50 μ= .1462 σχ= .0003 δ(min)= .1456 δ(max)= .1468
5999999990 - 6000000088 : n= 50 μ= .1471 σχ= .0002 δ(min)= .1466 δ(max)= .1474
8000000000 - 8000000050 : n= 26 μ= .1486 σχ= .0002 δ(min)= .1481 δ(max)= .1490
10000000000 -100亿098 : n= 50 μ= .1494 σχ= .0002 δ(min)= .1491 δ(max)= .1497
15000000000 -150亿098 :n= 50 μ= .15159 σχ= .00014 δ(min)= .1511 δ(max)= .15185
素数对下界计算值的计算实例:
G(110000000000) = 180801081;
inf( 110000000000 )≈ 180550355.5 , Δ≈-0.001387 ,infS( 110000000000 )= 121871489.95 ,
G(110000000002) = 122052830;
inf( 110000000002 )≈ 121871490 , Δ≈-0.001486 ,infS( 110000000002 )= 121871489.95 ,
G(110000000004) = 250274235;
inf( 110000000004 )≈ 249916814.3 , Δ≈-0.001428 ,infS( 110000000004 )= 121871489.95 ,
G(110000000006) = 133138114;
inf( 110000000006 )≈ 132950716.3 , Δ≈-0.001408 ,infS( 110000000006 )= 121871489.95 ,
G(110000000008) = 129058444;
inf( 110000000008 )≈ 128868117.6 , Δ≈-0.001475 ,infS( 110000000008 )= 121871489.96 ,
G(110000000010) = 325654239;
inf( 110000000010 )≈ 325204309 , Δ≈-0.001382 ,infS( 110000000010 )= 121871489.96 ,
G(110000000012) = 156839107;
inf( 110000000012 )≈ 156621995.1 , Δ≈-0.001384 ,infS( 110000000012 )= 121871489.96 ,
G(110000000014) = 122060507;
inf( 110000000014 )≈ 121884990.7 , Δ≈-0.001438 ,infS( 110000000014 )= 121871489.96 ,
G(110000000016) = 244091411;
inf( 110000000016 )≈ 243742979.9 , Δ≈-0.001427 ,infS( 110000000016 )= 121871489.97 ,
G(110000000018) = 122058317;
inf( 110000000018 )≈ 121890323.5 , Δ≈-0.001376 ,infS( 110000000018 )= 121871489.97 ,
G(110000000020) = 165628934;
inf( 110000000020 )≈ 165382515.2 , Δ≈-0.001488 ,infS( 110000000020 )= 121871489.97 ,
G(110000000022) = 271221025;
inf( 110000000022 )≈ 270825533.3 , Δ≈-0.001458 ,infS( 110000000022 )= 121871489.97 ,
计算式:
inf( 110000000000 ) = 1/(1+ .162 )*( 110000000000 /2 -2)*p(m) ≈ 180550355.5 ,
inf( 110000000002 ) = 1/(1+ .162 )*( 110000000002 /2 -2)*p(m) ≈ 121871490 ,
inf( 110000000004 ) = 1/(1+ .162 )*( 110000000004 /2 -2)*p(m) ≈ 249916814.3 ,
inf( 110000000006 ) = 1/(1+ .162 )*( 110000000006 /2 -2)*p(m) ≈ 132950716.3 ,
inf( 110000000008 ) = 1/(1+ .162 )*( 110000000008 /2 -2)*p(m) ≈ 128868117.6 ,
inf( 110000000010 ) = 1/(1+ .162 )*( 110000000010 /2 -2)*p(m) ≈ 325204309 ,
inf( 110000000012 ) = 1/(1+ .162 )*( 110000000012 /2 -2)*p(m) ≈ 156621995.1 ,
inf( 110000000014 ) = 1/(1+ .162 )*( 110000000014 /2 -2)*p(m) ≈ 121884990.7 ,
inf( 110000000016 ) = 1/(1+ .162 )*( 110000000016 /2 -2)*p(m) ≈ 243742979.9 ,
inf( 110000000018 ) = 1/(1+ .162 )*( 110000000018 /2 -2)*p(m) ≈ 121890323.5 ,
inf( 110000000020 ) = 1/(1+ .162 )*( 110000000020 /2 -2)*p(m) ≈ 165382515.2 ,
inf( 110000000022 ) = 1/(1+ .162 )*( 110000000022 /2 -2)*p(m) ≈ 270825533.3 ,
注:
只有在样本的相对误差统计计算中的标准偏差σχ比较小的情况下才能够用修正系数对一个比较大范围内的偶数的素数对进行高精度的计算,若标准偏差σχ比较大,说明区域偶数的素对计算值的相对误差的波动比较大,那么修正系数就可能只能适合部分偶数而不适合另外的一些偶数,就无从达到全部偶数的高精度的计算。
|
|