|
楼主 |
发表于 2020-8-10 14:03
|
显示全部楼层
本帖最后由 愚工688 于 2020-8-10 06:08 编辑
使用计算式:Xi(M)=t2*c1*M/(logM)^2 进行计算,式中的误差动态修正系数
t2=1.358-log(M)^(0.5)*.05484;
来修正哈-李公式与真值之间的误差偏差,能保持在一定范围内达到比较高的计算精度。
如果使用连乘式进行计算,则我是按照一定大小区域的偶数 取一个相对误差修正系数μ值,计算值乘以1/(1+μ)进行修正的。当然也能够如同t2那样用一个解析式来近似计算误差修正系数μ值,只是没有做。
查表法也是数学上面常用的方法,诸如以前使用的四位数学用表中就有按照尾值查表取误差修正值的。
例:连乘式的相对误差修正:
Sp( 20200809010 *) = 1/(1+ .1533 )*( 20200809010 /2 -2)*p(m) ≈ 34854492.9 , k(m)= 1.346801
Sp( 20200809012 *) = 1/(1+ .1533 )*( 20200809012 /2 -2)*p(m) ≈ 54259744.4 , k(m)= 2.096633
Sp( 20200809014 *) = 1/(1+ .1533 )*( 20200809014 /2 -2)*p(m) ≈ 25879461 , k(m)= 1
Sp( 20200809016 *) = 1/(1+ .1533 )*( 20200809016 /2 -2)*p(m) ≈ 25879461 , k(m)= 1
Sp( 20200809018 *) = 1/(1+ .1533 )*( 20200809018 /2 -2)*p(m) ≈ 56763031.4 , k(m)= 2.193362
G(20200809000) = 71772475;
Sp( 20200809000 *)≈ 71740428.5 , jdz =sp(m)/s(m) ≈0.999553
G(20200809002) = 28764472;
Sp( 20200809002 *)≈ 28754956.7 , jdz =sp(m)/s(m) ≈0.999669
G(20200809004) = 26288489;
Sp( 20200809004 *)≈ 26277606.6 , jdz =sp(m)/s(m) ≈0.999586
G(20200809006) = 65790461;
Sp( 20200809006 *)≈ 65764277.4 , jdz =sp(m)/s(m) ≈0.999602
G(20200809008) = 26532564;
Sp( 20200809008 *)≈ 26520461.2 , jdz =sp(m)/s(m) ≈0.999544
G(20200809010) = 34868882;
Sp( 20200809010 *)≈ 34854492.9 , jdz =sp(m)/s(m) ≈0.999587
G(20200809012) = 54278231;
Sp( 20200809012 *)≈ 54259744.4 , jdz =sp(m)/s(m) ≈0.999659
G(20200809014) = 25889107;
Sp( 20200809014 *)≈ 25879461 , jdz =sp(m)/s(m) ≈ 0.999627
G(20200809016) = 25890731;
Sp( 20200809016 *)≈ 25879461 , jdz =sp(m)/s(m) ≈ 0.999565
G(20200809018) = 56789952;
Sp( 20200809018 *)≈ 56763031.4 , jdz =sp(m)/s(m) ≈ 0.999526
|
|