|
因为素对用连乘式计算式,随着偶数的增大,区域偶数的计算值相对误差的均值会缓慢增大,因此修正值 μ 也将缓慢增大。
μ=0.21 是计算全部偶数的素对下界值时使用,因为小偶数时连乘积式的相对误差的波动比较大些;个别偶数的计算值偏离0位比较大,要使得下界式成立而修正值 μ必须达0.21;
在偶数趋大到十万亿以上时,连乘积式的相对误差值基本集中在0.175以上;
在偶数趋大到五十万亿以上时,连乘积式的相对误差值基本集中在0.18左右;
在偶数趋大到二百万亿以上时,连乘积式的相对误差值基本集中在0.182左右;
这是几乎是目前我的计算机能够计算的偶数极限了。
当然计算的偶数范围过大,必然有些区域偶数的计算值精度略有降低,可以用内插法增加一些 μ值以缩小偶数的计算范围;
用比较大偶数区域的 μ值代入计算相对比较小些区域的偶数的素对数量,得到的必然是略小于真值的下界素对计算值,而比较 μ=0.21则会有精度方面的提高。
(用 μ=0.21计算几百亿达的偶数的下界值,通常精度在95%以上。)
至于【修正系数=1/1.065 ? 笔误?范围 】—— 计算30万左右的偶数可以使用。但是比较小偶数的相对误差波动大些,达不到高精度计算。
比如:
G(280000) = 2225;inf( 280000 )≈ 2196.0 , Δ≈-0.0130,infS( 280000 )= 1372.48 , k(m)= 1.6
G(280002) = 2899;inf( 280002 )≈ 2875.7 , Δ≈-0.0080,infS( 280002 )= 1372.49 , k(m)= 2.09524
G(280004) = 1395;inf( 280004 )≈ 1372.5 , Δ≈-0.0161,infS( 280004 )= 1372.50 , k(m)= 1
G(280006) = 1392;inf( 280006 )≈ 1379.8 , Δ≈-0.0088,infS( 280006 )= 1372.51 , k(m)= 1.00529
G(280008) = 2804;inf( 280008 )≈ 2745.0 , Δ≈-0.0210,infS( 280008 )= 1372.52 , k(m)= 2
G(280010) = 1856;inf( 280010 )≈ 1830.0 , Δ≈-0.0140,infS( 280010 )= 1372.53 , k(m)= 1.33333
G(280012) = 1398;inf( 280012 )≈ 1372.5 , Δ≈-0.0182,infS( 280012 )= 1372.54 , k(m)= 1
G(280014) = 3455;inf( 280014 )≈ 3382.1 , Δ≈-0.0211,infS( 280014 )= 1372.55 , k(m)= 2.46411
G(280016) = 1658;inf( 280016 )≈ 1606.9 , Δ≈-0.0308,infS( 280016 )= 1372.56 , k(m)= 1.17073
G(280018) = 1381;inf( 280018 )≈ 1372.6 , Δ≈-0.0061,infS( 280018 )= 1372.57 , k(m)= 1
可以清楚的看到,区域下界计算值 infS(m) 随偶数增大而线性增大的特征。
计算式:
inf( 280000 ) = 1/(1+ .065 )*( 280000 /2 -2)*p(m) ≈ 2196 , k(m)= 1.6
inf( 280002 ) = 1/(1+ .065 )*( 280002 /2 -2)*p(m) ≈ 2875.7 , k(m)= 2.09524
inf( 280004 ) = 1/(1+ .065 )*( 280004 /2 -2)*p(m) ≈ 1372.5 , k(m)= 1
inf( 280006 ) = 1/(1+ .065 )*( 280006 /2 -2)*p(m) ≈ 1379.8 , k(m)= 1.00529
inf( 280008 ) = 1/(1+ .065 )*( 280008 /2 -2)*p(m) ≈ 2745 , k(m)= 2
inf( 280010 ) = 1/(1+ .065 )*( 280010 /2 -2)*p(m) ≈ 1830 , k(m)= 1.33333
inf( 280012 ) = 1/(1+ .065 )*( 280012 /2 -2)*p(m) ≈ 1372.5 , k(m)= 1
inf( 280014 ) = 1/(1+ .065 )*( 280014 /2 -2)*p(m) ≈ 3382.1 , k(m)= 2.46411
inf( 280016 ) = 1/(1+ .065 )*( 280016 /2 -2)*p(m) ≈ 1606.9 , k(m)= 1.17073
inf( 280018 ) = 1/(1+ .065 )*( 280018 /2 -2)*p(m) ≈ 1372.6 , k(m)= 1
|
|