5800计算器程序辛普森公式程序.docx

”XLZB”：坐标计算主程序

V=19:Fixm↓（与预留的自由变量数有关！）

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

Q=0↓

LbI1↓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 ：Pause0:”Y=” ：Y▲(X”X=”：Pause0:Y”Y=” ：

▲)

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

Goto1

(注释：括号内写法为Casio4800P计算器显示方式)

“ZB”:坐标计算子程序V=19:Fixm↓（与预留的自由变量数有关！）W=21↓（与预留的自由变量数有关！）

｛K｝:↓LbI2↓F=0：

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

段号；方框内数字39.903为断链值，短链输“+”，长链输“-”;如遇断链，依此增加即可！）

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

K≥Z[W+3a]=K≤E=V=W－1：≠＞W=W+1:Goto2△△↓K≥24200=K=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就可以满足要求)

LbI0↓

F=0：

T=6=F=39.903：△（注释：方框内数字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=”ENDOFTHEPROGRAM” ▲

△

Goto0

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

曲线元数据库设置表

名

称

第一段

第二段

……

第i段

……

第a-1段

第a段

曲线元长度

21

22

……

20+i

…