5800计算器程序辛普森公式程序(通用)

1．”XLZB”：坐标计算主程序

↓（与预留的自由变量数有关！）

｛O｝∶O“1.ZX 2.BX”↓

Q=0↓

LbI 1↓

O≠1=>｛QG｝:Q”L”:G”XJJD”:△（显示“L”输入边桩宽度，均为正值，显示“XJJD”输入边桩测点与线路中线上点的连线与线路中线上点的切线的夹角，如右法线为90度，左法线为270度或者-90度，其余同理。）

Prog”ZB”:X=X+Rec(Q,T+G):Y=Y+J:

”X=”：X ：Pause 0 :”Y=”：Y▲(X”X=”：Pause 0 : Y”Y=”：▲)

”T=”：T→DMS▲(T”T=”▲)

Goto 1

(注释：括号内写法为Casio4800P计算器显示方式) 2．“ZB”:坐标计算子程序

↓（与预留的自由变量数有关！）

（与预留的自由变量数有关！）

｛K｝:↓

LbI 2↓

F=0：

△（注释：方框内数字27[第七段]为断链的

段号；方框内数字39.903为断链值，短链输“+”，长链输“-”;

如遇断链，依此增加即可！）

E= Z[W+3a]+Z[W]+F：

K≥Z[W+3a] =>K≤E =>V=W－1：≠＞W=W+1: Goto 2△△↓

K－F：△（方框为断链点的后桩号，即去向桩号，如遇断链，依此增加即可！）

S=K－Z[V+3a+1]:P=(Z[V+2a+1] －Z[V+a+1])÷Z[V+1]:

I=P*S:C=Z[V+4a+1]:D=Z[V+a+1]:

J=C+(I+2D)*S*90/π:L=1：

X=Z[V+5a+1]+S*(CosC+CosJ+4∑（Cos (C+((L+0.5)I/B+2D)×(L+0.5)*90*S/B/π),L,0,(B-1)）+2∑（Cos (C+(L*I/B+2D)*L*90*S/B/π),L,1,(B-1)))/6/B:

Y=Z[V+6a+1]+S*(SinC+SinJ+4∑（Sin (C+((L+0.5)I/B+2D)×(L+0.5)*90*S/B/π),L,0,(B-1)）+2∑（Sin (C+(L*I/B+2D)*L*90*S/B/π),L,1,(B-1)))/6/B:

T=J:

３．“SJK.DAT”:数据库主程序

Mcl:Defm(20+7a):Z[21]=?:……Z[21+3a]=?:Z[21+4a]=?:Z[21+5a]=?:Z[ 21+6a]=?:Prog”SJK”:(输入程序时，Defm后须换成相应数据，a为曲线元段数)

4．“SJK”：数据库计算子程序

a=?:V=20:T=0:B=?(a:曲线元段数；B:等分数，最少为2，一般

为4就可以满足要求)

LbI 0↓

F=0：

（注释：方框内数字6[第七段减一]为断链的段号减去1所得；方框内数字39.903为断链值，短链输“+”，长链输“-”;如遇断链，依此增加即可！）

E= Z[V+3a+1]+Z[V+1]+F:

Z[V+3a+2]=E：

S=Z[V+1] :

I=Z[V+2a+1] －Z[V+a+1]:

C=Z[V+4a+1]: D=Z[V+a+1]:

J= C+(I+2D)*S*90/π:

Z[V+4a+2]=J:L=1：

Z[V+5a+2]=Z[V+5a+1]+S*(CosC+CosJ+4∑（Cos (C+((L+0.5)I/B+2D)×(L+0.5)*90*S/B/π),L,0,(B-1)）+2∑（Cos (C+(L*I/B+2D)*L*90*S/B/π),L,1,(B-1)))/6/B:

Z[V+6a+2]=Z[V+6a+1]+S*(SinC+SinJ+4∑（Sin (C+((L+0.5)I/B+2D)×(L+0.5)*90*S/B/π),L,0,(B-1)）+2∑（Sin (C+(L*I/B+2D)*L*90*S/B/π),L,1,(B-1)))/6/B:

V=V+1:T=T+1:T＞a-2 =>”END OF THE PROGRAM”▲

△

Goto 0

为了使程序能够计算连续多个曲线元的坐标，减少计算过程中的数据输入，可利用CASIOfx-4800/4850P计算器提供的扩展变量设置曲线元数据库，假定某一段线路包含有a个曲线元，则可设置曲线元数据库如下表所示：

曲线元数据库设置表

扩展变量Z[1]~Z[20]用作自由变量，可以用来弥补字母变量的不足，此项设置可自由选择！因此，数据库的起点本程序从Z[21]开始，表中Z[21]~Z[21+3a]以及Z[21+4a]、Z[21+5a]、Z[21+6a]需预先置入，表中空白部分由程序自动计算生成。

CASIOfx-4800/4850P计算器程序

源程序

该算法共有四个程序：1.线路坐标计算主程序（XLZB）；2.数据

库生成主程序(SJK.DAT)；3.坐标计算子程序(ZB)；4.数据库计算子程

序(SJK)。

程序使用说明

不管是整体式还是分离式，主线还是匝道，都能用同一个计算器就可进行各路段相关坐标计算，只是分别采用了各自路

段的数据库，如果再结合施工放样正反算程序及高程程序、隧

道及路基开挖断面程序，则可形成一个完善实用的施工测量程

序包，特别适合于公路、铁路现场施工测量人员使用。

在计算某一段线路时，首先必须运行相应路段数据库生成主程序(SJK.DAT)，然后运行线路坐标计算主程序（XLZB），按

照程序提示输入是计算中线还是边线，若是计算边线除需输入

桩号外，还需输入宽度及斜交角度。因此，根据线路实际情况

可编制相应路段数据库生成主程序即可完成相应路段坐标计算

工作。

算例

本文提供聂让编著《全站仪与高等级公路测量》306~308页算例的计算结果如下：