道路直线段中边桩坐标放样计算CASIO fx-5800P程序090405

CASIO5800计算器中边桩坐标计算程序一．GSGL（主程序）1.“K0”? W:“X0”? U:“Y0”? V:“R”?R: “LH”?L:“A0”?A:“B”?B:“R=1,L=-1”?N:“S-BACK"?E:“T-BACK"?F:2.Lb1 0:3.“KI”? K: E-F S:4.(R+L2÷(24R)-L^(4)÷(2688R^(3)))×tan(B÷2)+L÷2-L^(3)÷(240R2)T:5.S-T H:K-W J:J＜0 =〉Goto 0:6.If H>J: Then U+Jcos(A)X:V+Jsin(A)Y:7.“XI=”:X “YI=”:Y8.A G:Prog “SUB”: Goto 0: IfEnd:9.L=0 =〉Goto 1:10.K-H-W Z: Z-Z^(5)÷(90R2L2)C:11.30 Z2÷(πLR)D: Z＞L =〉Goto 1:12.U+Hcos(A)+Ccos(A+ND)X:V+Hsin(A)+Csin(A+ND)Y:13.“XI=”:X“YI=”:Y14.A+180NZ2÷(2πLR)G:Prog “SUB”: Goto 0:15.Lb1 1:16.πRB÷180-L C: K-H-W-L M: 90M÷(πR)I: 2Rsin(I)O:17.L-L^(3)÷(90R2)P:30LN÷(πR)Q:A+3Q+NI D: M>C =〉Goto 2:18.U+Hcos(A)+Pcos(A+Q)+Ocos(D)X: V+Hsin(A)+Psin(A+Q)+Osin(D)Y:19.“XI=”:X “YI=”:Y20.A+N(180L÷2÷π÷R+180M÷(πR))G: Prog “SUB”: Goto 0:21.Lb1 2:22.L=0 =〉Goto 0:23.A+NB G:24.H-J+L+πRB÷180P:P<-.003 =〉Goto 0:25.P-P^(5)÷(90R2L2)C: 30P2÷(πRL)D:A+NB-180-ND Q:26.U+Scos(A)+Tcos(G)+Ccos(Q)X: V+Ssin(A)+Tsin(G)+Csin(Q)Y:27.“XI=”:X “YI=”:Y28.A+NB-180NP2÷(2πRL)G: Prog “SUB” :Goto 0:二．SUB（子程序）1．Lb1 0:2. 0I：0O：0M:3. “Y=1，N=0”? I：I≠1 =〉Goto 2：4. G+90C：G-90D：“R=1，L=-1”? O：“S-LONG”? M：5. If O=1：Then X+Mcos(C)P : Y+Msin(C)Q：6. “XR=”：P “YR=”：Q7. Else X+Mcos(D)P: Y+Msin(D)Q：8. “XL=”：P “YL=”：Q IfEnd:9. Goto 0：10. Lb1 2说明1）K0：起点的桩号2）X0：起点的X坐标3）Y0：起点的Y坐标4）KI：需计算的桩号5）R：圆曲线半径6）LH：缓和曲线段长度7）A0：起点方位角（角度）8）B：偏角（角度）9）R=1，L=-1：表示路线左右偏，右1，左-110）S-BACK：起点至交点长度11）T-BACK：一般取012）Y=1，N=0：是否计算边桩，是1，否013）R=1，L=-1：左右边桩，右1，左-114）S-LONG：表示中桩至边桩距离。

工程测量放样CASIO fx-5800P 直线正反算坐标计算程序作者：张素辉来源：《城市建设理论研究》2013年第32期【摘要】介绍工程测量放样中CASIO fx-5800P直线正反算计算程序，主要应用于结构物的坐标计算。

【关键词】坐标；计算；CASIO；程序中图分类号：P123 文献标识码：A【引言】对于一个工程测量技术人员来说，能够掌握和运用测量基础知识，正确、熟练地使用计算器（如CASIO fx-5800P）编程和运用，才能在野外方便快捷的计算、复核大地坐标，并能结合全站仪准确、方便的放样。

在此谨以我在施工测量中编辑、调试、实用的CASIO fx-5800P 直线正反算程序奉送读者。

一、主程序名：ZXZFSDeg：Fix 3：20→DimZ↙（设置程序参数）＂DK(QD)＂?A（输入起点桩号）：＂X(QD)＂？B（输入起点坐标X）：＂Y(QD)＂？C （输入起点坐标X）：＂X(ZD)＂?U（输入终点坐标X）：＂Y(ZD)＂?V↙（输入终点坐标Y）“DJ=”：POL(U-B,V-C)：I → M◢(显示直线平距)“FWJ=”：if JGoto0↙ElseJ → N:N▲DMS◢ (显示直线方位角)Goto0 ↙(直线参数设置完成)Lbl 0↙N→F：＂1=>XY.2=>DK＂：＂J＂？W↙（条件选择：输入数字1为正算，输入数字2为反算）If W=1：Then Goto 1：Else If W=2：Then Goto 2：IfEnd↙Lb1 1↙（正算）＂JSDK＂？E（输入待求点的桩号）：＂PJ＂?G（输入待求点偏距）：＂PA＂？H↙（输入待求点偏角左-右+）＂X=＂：B+(E-A) ×cos(F)+G×cos(F+H)→X◢＂Y=＂：C+(E-A) ×sin(F)+G×sin(F+H)→Y◢（计算出待测点的坐标X,Y）＂CZD(X)＂？S ↙（设置全站仪测站坐标X）IfS=0：ThenGoto 1：Else （当测站坐标X输0时，不计算放样数据）＂CZD(Y)＂？T：↙（设置全站仪测站坐标Y）＂FY-PJ=＂：Pol(X-S，Y-T)：I→M◢＂FY-FWJ=＂：If JGoto1↙Else J→Z[1]：Z[1] ▲DMS◢（计算出放样平距及放样坐标方位角）Goto 1↙（正算循环）Lb1 2↙（反算）＂X=＂？K：＂Y=＂？L↙（输入全站仪测量点的坐标X,Y）Pol(K-B，L-C)：I→D：If JIf FSin(Ｏ)×D→Q↙＂DK=＂：A+Q÷tan(Ｏ)→R◢(反算测量点的桩号)＂PJ=＂：sin(Ｏ)×D→Q：If FGoto 2↙Else -1×Q→Q◢（反算出测量点距离直线的偏距：负值为左偏，正值为右偏） Goto2↙ ·（反算循环）备注：1、◢为显示指令。

利用CASIO fx-5800P计算常规道路上中桩和外移桩坐标
代科
【期刊名称】《中国高新科技》
【年(卷),期】2018(000)014
【摘要】常规道路曲线由直线段、缓圆段、圆曲线段和圆缓段组成,该程序可计算此类常规道路曲线中桩坐标以及外移桩坐标,其最大的特点就是在计算外移桩坐标
时引入了相对于中桩两个方向上的偏移量。

文章首先对程序中曲线参数说明进行了分析,并阐述了程序中曲线参数说明,然后对主程序L-COORDINATE、子程序L-PARAMETER进行探讨,最后通过实例进行计算,减少了在计算线路外移桩时的工作量,提高了工作效率。

【总页数】4页(P115-118)
【作者】代科
【作者单位】中铁北京工程局集团有限公司,北京100000
【正文语种】中文
【中图分类】TU198
【相关文献】
1.CASIO fx-5800p计算器在公路互通立交匝道坐标计算中的应用 [J], 张鹏飞;周庆来;曾晟
2.CASIO fx-4500 计算器计算公路曲线要素及公路中桩、边桩坐标的综合程序 [J], 杨宗炉
3.利用CASIO fx-5800P计算常规道路上中桩和外移桩坐标 [J], 代科
4.CASIO fx-5800P公路坐标计算程序的设计与实践 [J], 向继平
5.CASIO fx-4800p编程计算路线中桩及边桩坐标 [J], 齐兴强
因版权原因，仅展示原文概要，查看原文内容请购买。

CASIO fx-4800P、fx-5800P型计算器用于线路施工曲线中线点坐标的计算程序中铁十局三建公司工程技术部摘要：本文介绍了CASIO fx-4800P 、fx-5800P型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。

本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。

关键词：曲线线路施工测量计算程序1．概述过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。

对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。

随着电子技术进步和经济发展，测量仪器和测量方法的不断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。

但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。

根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。

前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。

而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。

为了能够快速准确地为全站仪提供测设数据，发挥全站仪快速测设的特点，提高测量工作效率，应采用可编程的计算器，编制计算程序。

本文主要介绍应用CASIO fx-4800P型计算器的计算程序，供公司测量同行们参照使用。

2．计算程序QXZBJS（文件名：曲线坐标计算fx-4800P）Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+M ◢G=RAπ/180：“S”S=G+L◢LbiA：｛C，V｝：C“CSDLC=”：V“HXPJ=”：D=C-K：D≤L＝＞I=D-D5/（40R2L2）：U=D3/(6RL)-D7/（336R3L3):J=√（I2+U2）：Goto1：≠＞D≤G＝＞O=90（2D-L）/( Rπ)：I=RsinO+M：U=R（1-cosO）+P：J=√（I2+U2）：Goto2：≠＞D=S-（C-K）： = D-D5/（40R2L2）：Z[2]=D3/6RL-D7/（336R3L3）：I=T+（T-Z[1]）cosA-Z[2]sinA：U=（T-Z[1]）sinA+Z[2]cosA：J=√（I2+U2）：Goto3：Lbi1:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-30D2/ （RLπ）：H=F-90D2/ （RLπ）：≠＞Q=F+30D2/ （RLπ）：H=F+90D2/（ RLπ）⊿ Goto4：Lbi2:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I)：H=F-O：≠＞Q= F+tng-1(U/I)：H=F+O⊿Goto4：Lbi3:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I):H=F-(A-90(S-(C-K))2/ （RLπ）): ≠＞Q= F+tng-1(U/I):H=F+(A-90(S-(C-K))2/ （RLπ）):⊿ Goto4：Lbi4:B=90+H:H＜0＝＞H“QXFWJ”=B+360◢≠＞H≥360＝＞H“QXFWJ”=H-360 ◢≠＞H“QXFWJ”=H◢⊿Goto5：Lbi5: X“CSD:X”=JcosQ+N+VcosB◢ Y“CSD:Y”=JsinQ+E+VsinB◢GotoA3．程序说明3.1 输入已知变量R—圆曲线半径，显示R？L—缓和曲线长，显示L？A—曲线转向角，显示A？E—直缓点纵坐标，显示ZH：X=？N—直缓点横坐标，显示ZH：Y=？F—第一切线方位角，即ZH至JD的方位角，显示FWJ=？K—直缓点里程，显示ZH：LC=？3.2 计算待求量T—切线长度，显示T= …S—曲线全长，显示S= …Z[3]—外矢距，即JD到QZ的距离，显示E0= …3.3 输入待求变量K—输入待求（测设）点的里程，显示LC=？V—横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。

单交点单圆曲线中边桩坐标计算程序(PM4-3)(2) 主程序——PM4-3，占用内存1070字节。

"SINGLE CIRCLE CURVE"显示程序标题1"METHOD OF COORDINATE PM4-3"显示坐标法标题2Deg:ClrStat:FreqOn:Fix 3基本设置14DimZ定义额外变量"JD MILEAGE PEG(m)="?Z输入以m为单位的交点桩号"JD X(m)="?U输入交点x坐标"JD Y(m)="?V输入交点y坐标"ZD X(m)="?M输入转点x坐标"ZD Y(m)="?N输入转点y坐标Pol(U-M,V-N):Cls计算ZD JD的边长与方位角If J<0:Then J+360A:Else J A:IfEnd判断ZD JD的方位角A"TURNING ANGLE -L,+R(Deg)="?Q输入交点路线转角,左偏为负,右偏为正If Q<0:Then -1Z[10]:Else 1Z[10]:IfEnd确定转角数值Q Z[4]:Abs(Q)D保存转角及符号"R(m)="?R输入圆曲线半径Rtan(0.5D)T计算切线长RDπ÷180L计算曲线长R(cos(0.5D)-1-1)E计算外距2T-L J计算切曲差"T(m)=":T显示切线长"L(m)=":L显示曲线长"E(m)=":E显示外距"J(m)=":J显示切曲差Z-T Z[1]:Z[1]+0.5L Z[2]:Z[2]+0.5L Z[3]计算ZY，QZ，YZ点桩号Z[1]List X[1]存储ZY点的桩号到统计串列U-Tcos(A)List Y[1]:V-Tsin(A)List Freq[1]计算并存储ZY点的测量坐标到统计串列"ZY PEG(m)=":List X[1]显示ZY点桩号及坐标"ZY-X(m)=":List Y[1]"ZY-Y(m)=":List Freq[1]"WL(m),0 NO="?B输入路面左半幅宽度If B>0:Then List Y[1]+Bcos(A-90)Z[11]计算并显示ZY点的左边桩坐标List Freq[1]+Bsin(A-90)Z[12]"XL(m)=":Z[11]"YL(m)=":Z[12]IfEnd"WR(m),0 NO="?C输入路面右半幅宽度If C>0:Then List Y[1]+Ccos(A+90)Z[13]计算并显示ZY点的右边桩坐标List Freq[1]+Csin(A+90)Z[14]"XR(m)=":Z[13]"YR(m)=":Z[14]IfEndZ[2]List X[2]:2K存储ZY点的桩号到统计串列Prog "SUB4-31"调子程序计算并显示QZ点坐标Z[3]List X[3]:3K存储YZ点桩号Prog "SUB4-31"调子程序计算并显示YZ点坐标If R<30:Then 5I:Else If R<60:Then 10I:Else 20I:IfEnd:IfEnd确定整桩间距"INT DIST(m)="?O输入整桩间距If O>0:Then O I:IfEnd输入的整桩间距大于零时，使用输入值"ONLY CALC +PEG(1)"?G输入1为只计算加桩点的坐标If G=1:Then Goto 1:IfEndK+1K计数变量计数Int(Z[1]÷I)I+I List X[K]计算并存储从ZY点开始的第一个整桩号Prog "SUB4-31"调子程序计算并显示坐标Do计算ZY点至YZ点K+1K计数变量计数List X[K-1]+I List X[K]计算并存储整桩号Prog "SUB4-31"调子程序计算并显示坐标LpWhile List X[K]+I<Z[3]没计算到YZ点时继续循环Lbl 1Do"+PEG(m)="?F输入加桩号If F List X[1]:Then Break:IfEnd加桩号小于ZH点桩号时结束程序运行If F List X[3]:Then Break:IfEnd加桩号大于HZ点桩号时结束程序运行K+1K计数变量计数F List X[K]存储加桩号到统计串列Prog "SUB4-31"调子程序计算并显示偏角弦长LpWhile F>0没有计算到YZ点时继续循环"PM4-3END"(3) 子程序——SUB4-31，占用内存498字节。

Casio5800交点法程序（输入版）（歪哥哥2009版）本程序由一个主程序JD和三个子程序（JDA、JDB、JDC）构成，运行时只需运行主程序即可！本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标计算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并为线元法程序提供起点坐标起点切线方位角等数据！当然本程序也可单独逐交点输入进行放样计算用！鉴于5800计算器的空间和以上所述本程序的主要目的，故此程序不修改为数据库版本！需要的自行修改结合XY框架自己修改为数据库反算程序等！主程序名：5.JD24→Dimz↙Cls :＂XC＂?U :＂YC＂?V :＂K（JD）＂?K :＂X（JD）＂?X :＂Y（JD）＂?Y :＂LS1＂?B :＂LS2＂?C : ?R :＂（ZH）FWJ°＂?M : ＂α（Z-,Y+）°＂?O : M+O→N :Prog ＂JDA＂↙Cls :＂T1=＂:＂T2=＂:＂L=＂:＂L Y=＂: Locate 4,1,S : Locate 4,2,T : Locate 4,3,L : Locate 4,4,Q◢Cls :＂E=＂:＂K（ZH）=＂: Locate 7,1,E : Locate 7,2,Z[1]◢Cls : ＂K（HY）=＂:＂K（QZ）=＂:＂K（YH）=＂:＂K（HZ）=＂: Locate 7,1, Z[2] : Locate 7,2, Z[3] : Locate 7,3, Z[4] : Locate 7,4, Z[5]◢LbI 0 : ＂K×+×××＂?P : ＂Z＂?D : If D≠0 :Then ＂RJ＂?H : IfEnd : Prog ＂JDB＂↙If D＜0 :Then Cls : ＂X(L)=＂:＂Y(L)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(L)=＂: Locate 6,1,I : ＂F(L)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙If D=0 :Then Cls : ＂X(Z)=＂:＂Y(Z)=＂: Locate 6,1,F : Locate 6,2,G : ＂QXFWJ(Z)=＂: Z▼DMS◢Pol(F-U,G-V : Cls : ＂S(Z)=＂: Locate 6,1,I : ＂F(Z)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙If D＞0 :Then Cls : ＂X(R)=＂:＂Y(R)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(R)=＂: Locate 6,1,I : ＂F(R)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙子程序1名：JDAIf O＜0 :Then -1→W : Else 1→W : IfEnd : WO→A ↙B2 ÷24÷R-B^(4)÷2688÷R ^(3) →Z[6] ↙C2 ÷24÷R-C^(4)÷2688÷R ^(3) →Z[7] ↙B÷2-B^(3)÷240÷R2 →Z[8] ↙C÷2-C^(3)÷240÷R2 →Z[9] ↙Z[8]+((R+Z[7]-(R+Z[6])cos(A))÷sin(A))→S↙Z[9]+((R+Z[6]-(R+Z[7])cos(A))÷sin(A))→T↙RAπ÷180+(B+C) ÷2→L↙RAπ÷180-(B+C) ÷2→Q↙(R+(Z[6]+Z[7])÷2)÷cos(A÷2)-R→E↙K-S→Z[1] ↙↙Z[1]+B→Z[2] ↙↙Z[2]+Q÷2→Z[3]↙Z[1]+L-C→Z[4]↙Z[4]+C→Z[5]↙子程序2名：JDBX-Scos(M)→Z[19]:Y-Ssin(M)→Z[20]↙X+Tcos(N)→Z[21]:Y+Tsin(N)→Z[22]↙If P＞Z[1]:Then Goto 1 :IfEnd↙Z[1]-P→L↙X-(S+L)cos(M)+Dcos(Z+H)→F↙Y-(S+L)sin(M)+Dsin(Z+H)→G↙M→Z : Goto 5↙LbI 1 : If P＞Z[2]:Then Goto 2 :IfEnd↙P-Z[1]→L:L→Z[12]:B→Z[13]:Prog＂JDC＂↙Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+H)→F↙Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+H)→G↙M+90WL2 ÷(BRπ)→Z↙Goto 5↙LbI 2 : If P＞Z[4]:Then Goto 3 :IfEnd↙P-Z[1]→L:90(2L-B)÷R÷π→Z[11]↙Rsin(Z[11])+Z[8]→Z[14]:R(1-cos(Z[11]))+Z[6]→Z[15]↙Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+H)→F↙Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+H)→G↙M+WZ[11]→Z↙Goto 5↙LbI 3 : If P＞Z[5]:Then Goto 4 :IfEnd↙Z[5]-P→L:L→Z[12]:C→Z[13]:Prog＂JDC＂↙Z[21]-Z[14]cos(N)-WZ[15]sin(N)+Dcos(Z+H)→F↙Z[22]-Z[14]sin(N)+WZ[15]cos(N)+Dsin(Z+H)→G↙N-90WL2 ÷(CRπ)→Z↙Goto 5↙LbI 4 : P-Z[5]→L↙X+(T+L)cos(N)+Dcos(Z+H)→F↙Y+(T+L)sin(N)+Dsin(Z+H)→G↙N→Z↙Goto 5↙LbI 5 : 360Frac((Z+360)÷360→Z↙子程序3名：JDCIf Z[12]=0 :Then 0→Z[14]: 0→Z[15]:Else↙Z[12]- Z[12]^(5)÷40÷(RZ[13])2+ Z[12]^(9)÷3456÷(RZ[13])^(4) →Z[14]↙Z[12]^(3)÷6÷(RZ[13])-Z[12]^(7)÷336÷(RZ[13])^(3)+ Z[12]^(11) ÷42240÷(RZ[13])^(5)→Z[15] ↙IfEnd↙程序说明：已知数据输入：XC ? 测站X坐标YC ? 测站Y坐标K（JD）?交点桩号X（JD）?交点X坐标Y（JD）?交点Y坐标LS1 ?第一缓和曲线长度LS2 ?第二缓和曲线长度R ? 圆曲线半径（ZH）FWJ°?交点前（即前交点至本交点也即ZH点）的正切线方位角α（Z-,Y+）?本交点处线路转角（左转为负，右转为正，度分秒输入）K×+×××? 待求桩号Z ?待求桩号距中距离（左负值，右正值，中为0）RJ ?斜交右角（线路切线前进方向与边桩右侧夹角）计算结果显示：T1=第一切线长T2=第二切线长L=曲线总长L Y=圆曲线长E=曲线外距K(ZH)=直缓点桩号K(HY)=缓圆点桩号K(QZ)=曲中点桩号K(YH)=圆缓点桩号K(HZ)=缓直点桩号X= Y=待求点的坐标（其中：L-左Z-中R-右）QXFWJ(Z)=待求点的中桩切线方位角（当求中桩坐标时显示）S= F=测站至待求点的水平距离、方位角（其中L-左Z-中R-右）（歪哥哥郑重声明：代码无误，输入注意！使用者使用不当所造成的一切损失概不负责！！）。

C A S I O f x-5800P计算与道路坐标放样计算
一、前言
本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。

原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。

改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是：
1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦；
2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）；
3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）；
4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现
程序的特点：
1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写；
2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；
3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。

程序代码
计算器操作正算说明记住要输入正算1
C匝道计算结果。

FX-5800p计算器公路测量常用程序一、程序功能本程序由6个主程序、5个次子程序及5个参数子程序组成。

主要用于公路测量中坐标正反算，设计任意点高程及横坡计算,桥涵放样，路基开挖口及填方坡脚线放样。

程序坐标计算适应于任何线型.二、源程序1.主程序1：一般放样反算程序(①正算坐标、放样点至置仪点方位角及距离；②反算桩号及距中距离)程序名:1ZD-XYLb1 0:Norm 2F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量)Z[1]=90（与路线右边夹角）Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2Lb1 1:Ｆix 3:＂Ｘ=＂：Locate 6,4,Ｘ◢＂Ｙ＝＂：Locate 6,4,Ｙ◢Ｐrog＂3JS”:Goto 0:Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢Ｇoto 0２．主程序２：高程序横坡程序(设计任意点高程及横坡)程序名：2GCLbI 0:Norm 2“KM”?Z:?D:Prog”H”:Fix3:”H=”:Locate 6,4,H◢“I=”: Locate 6,4,I◢Goto 03.主程序3：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：3JSX：Y：1268．123→K(置仪点X坐标)2243．545→L（置仪点Y坐标，都是手工输入,也可以建导线点数据库子程序,个人认为太麻烦）Y-L→E：X-K→F：Pol(F,E):IF J<0:ThenJ+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示)Fix 3:”S=”:Locate 6,4,I◢4．主程序4：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序3中输入置仪点坐标后计算放样点至置仪点方位角及距离)程序名：4JH-XYLbI 0:Norm 290→Z[1](涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量)1→F:Prog”THB”:?L:Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ=＂：Locate 6,4,Ｘ◢＂Ｙ＝＂：Locate 6,4,Ｙ◢Ｐrog＂3JS”:Goto 0:5．主程序5:路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名：5FBXLbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog“THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:GotoH:Else L-0.75→D:Goto H:IfEndLbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5] →G:Goto W:Else Z[5] →G:Goto T:LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else If G>Z[10]:Then Goto 2:Else Goto 3:IfEnd:LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z:LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:Goto Z:LbI 3:L+Z[12]+G×Z[7]:Goto z:LbI T:L+0.5→N:If G>Z[17]:Then (N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto Z:Else (N+G×Z[15])→S:Goto z:LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:”L0=”:L Locate 6,4,T◢＂Ｄ＝＂：Locate 6,4,Ｄ◢“TW=”: Locate 6,4,Z[5]◢Goto 06．主程序6：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至L)程序名：6GDProg “G0”Z-C→E:(B-A)×E/S+A→L:L:7．坐标计算次程序（ＴＨＢ）程序名：THB18→DimZ：＂ＫＭ＂?Z：Prog ＂Ｘ０＂1÷P→Ｃ： (P-R)÷(2HPR) →Ｓ：180÷π→Ｅ：Ｆ=1=>Goto 1：Ｆ=2=>Goto 2←┘Lbl 1：?Ｄ： Abs(Ｚ－O) →W：Prog "Ａ"：X：Ｙ：Ｇoto 3LbI 2:Ｘ：Ｙ：Ｘ→I：Ｙ→J：Prog ＂Ｂ＂：Ｏ＋Ｗ→Z：Ｄ→Ｄ：Ｇoto 3LbI 3:IF Ｆ＝１Then Ｘ：Ｙ：Else Ｚ：Ｄ8. 正算子程序(Ａ)程序名：A0.1184634425→A： 0.2393143352→B：0.2844444444→N 0.046910077→K：0.2307653449→L：0.5→M： U+W(Acos(G+QEKW(C+KWＳ))+Bcos(G+QELW(C+LWＳ))+Ncos(G+QEMW(C+MWＳ))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X：V+W(Asin(G+QEKW(C+KWＳ))+Bsin(G+QELW(C+LWＳ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WＳ))) →Y：G+QEW(C+WＳ)+Ｚ［１］→Ｚ［２］：X+Ｄcos（Ｚ［２］）→Ｘ： Y+Ｄsin（Ｚ［２］）→Ｙ9. 反算子程序(Ｂ)程序名：BG-90→Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) →Ｗ：0→Ｄ：Lbl 0：Prog "Ａ"： T+QEW(C+WＳ) →Ｌ： (J-Y)cos(L)-(I-X)sin(L)→Ｄ：IF Abs(Ｄ)<0.01:Then Goto1：Else W+Ｄ→W:Goto 0←┘Lbl 1:0→D：Prog "Ａ"：(J-Y)÷sin(Ｚ［２］) →D:10．高程计算子程序（H）程序名：HProg “S0”:R:T:C:G:I:C-T→F:Z-F→L:C+T→E:G-TI→Q:If T=O:Then Q+LI→H:Goto 0:Else IfZ<F:Then Q+LI→H:Goto 0:Else If Z≤E:Then Q+LI+L2÷2÷R→H:Goto 0:LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I:H:I:11.高程超高计算程序（I）程序名：IProg”I0”:W=1=> Goto 0:W=2=>Goto 1:LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd: LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×(((3((Z-C)÷L)2-2((Z-C)÷L)∧(3))×(N-M))+M)→V:Goto 2:IfEnd:LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:12．数据子程序(附后示例)①程序名：X0（坐标计算要素程序）If Z≥25900 And Z≤26615.555:Then 25900→O:11587.421→U:1847.983→V:101。

CASIOfx—5800P公路坐标计算程序的设计与实践作者：向继平来源：《电脑知识与技术》2018年第33期摘要：CASIO fx-5800P是测绘与土木工程领域广泛使用的计算器，特别适合现场计算。

本文通过比较分析，选取基于线元的Gauss-Legendre数学模型，结合CASIO fx-5800P计算器的特点，通过设计实用的公路坐标计算程序，提高现场计算的工作效率和精度。

关键词：Gauss-Legendre；线元；程序；矩阵；坐标计算中图分类号：TP311 文献标识码：A 文章编号：1009-3044（2018）33-0241-03在公路测量工作中，要经常对线路的中桩、边桩、结构物位置进行实地放样，这就需要在现场实时快速地获取放样点的坐标。

我国现行《公路路线设计规范》规定，公路平面线形由直线、圆曲线和缓和曲线（回旋线）三种要素组成。

在现代高速公路上，有时缓和曲线所占的比例超过了直线和圆曲线，成为平面线形的主要组成部分。

缓和曲线的计算复杂，尤其是非完整缓和曲线，计算器编程计算是解决公路坐标现场计算的有效手段。

1 数学模型的选择公路坐标计算的关键是数学模型的选择。

传统的计算采用交点法，它针对组成公路平面线形的直线、圆曲线和缓和曲线采用不同的数学模型进行计算，计算时要先建立各个计算单元（直线、圆曲线、缓和曲线）的相对坐标系，再根据点的桩号确定其所在的计算单元，采用对应的计算公式求出该点在相对坐标系中的坐标，然后转换为公路统一坐标系的坐标，计算过程复杂。

线元法将组成公路平面线形的直线、圆曲线和缓和曲线均看作独立的单元——线元。

直线可以看作是曲率半径为∞的曲线，圆曲线为曲率半径为定值的曲线，缓和曲线为曲率半径渐变的曲线，因此，公路中线可以看成是由曲线构成的。

曲线元上任一点的曲率都随弧长作线性变化，基于这一共性推导出的任意点坐标计算积分通式——Gauss-Legendre公式可作为公路坐标计算的通用数学模型，该公式与点所在的位置无关。

5800计算器公路全线坐标正、反算计算程序FX5800全线贯通万能正、反算程序（一体化、超好用、短小、易懂）FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版ZHUCHENGXU主程序"1.ZS,2.FS" ?→Q←┘输入1正算，输入2反算“NEW=0,OLD≠0”?Z←┘IfZ=0:Then“X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2←┘Lbl1 :“KM=,<0 Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?L←┘LblZ:Z=1=> Prog“01”：Z=2=> Prog“02”←┘选择数据库文件，可增加H- F→X:0.5(E-D)÷(G-F)→N←┘C+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→P←┘-A+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→U←┘B+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→V←┘Q=2=>Goto4：Cls：Fix 3←┘"Xn=":Locate4，1，U："Yn=": Locate5，2，V：“FWJ=”：PDMS◢Norm 2：Cls：Goto1←┘Lbl2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S←┘“KMDG=”？H ：90→O:0→L：GotoZ←┘（H线路范围内的任意桩号）Lbl4:Pol（R-U，S-V）：J<0 => J+360→J←┘Whileabs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:If P-J>0:then -I→L:else I→L :Ifend:Goto3: Whileend：H+Icos（J-P）→H:GotoZ←┘Lbl3:Cls：Fix 3←┘“KM=”: Locate4，1，H：“PY=”: Locate4，2，L◢Norm 2：Cls：Goto2←┘01（数据库子程序）If H<=第一曲线终点桩号：then第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend ……………程序说明：1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3…代替输入）；如果在数据库中没有要计算线路的数据，则输入曲线要素X0：曲线起点X坐标；Y0：曲线起点Y坐标；C0:曲线起点方位角；R0-1、 RI-1:曲线起点、终点曲率，直线为0，曲线左偏输入负值，右偏输入正值；SP、 EP：曲线起点桩号，终点桩号；KM：待求点桩号；PJ：正斜交的设定；PY：偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；4、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

CASIO fx—5800P计算器路线坐标计算程序（单个交点）主程序名CALXY （计算中桩、边桩坐标）Lbl 1：“ZZ”? →Z[27] ：“D(-+m)”? →D：If D≠0 ：Then “∠(DMS)”? →V：Ifend』If Z[27]≤Z[17]：Then Z[18]+ (Z[27] —Z[17]) Cos( F)+D Cos( F+V)：“X=”：Ans▲Z[19]+ (Z[27] —Z[17]) Sin( F )+D Sin( F+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[17] And Z[27]≤Z[20]：Then Z[27] —Z[17] →L：Prog “HHXY”：Z[18] +I Cos( F+ JH）+D Cos( F+WH+V)：“X=”：Ans▲：Z[19] +I Sin(F+ JH）+D Sin( F+WH+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[20] And Z[27]≤Z[24] —Z[6]：Then Z[27] —Z[20] →L：Prog “YUXY”：Z[21] +I Cos( Z[23]+ JH）+D Cos( Z[23]+WH+V)：“X=”：Ans▲：Z[22] +I Sin(Z[23]+ JH）+D Sin( Z[23]+WH+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[24] —Z[6] And Z[27]≤Z[24]：Then Z[24] —Z[27] →L：Prog “HHXY”：Z[25] +I Cos( C—JH）—D Cos( C—WH+V)：“X=”：Ans▲：Z[26] +I Sin(C—JH）—D Sin( C—WH+V)：“Y=”：Ans▲：Ifend』If Z[27] ＞Z[24]：Then Z[25] —(Z[27] —Z[24]) Cos( C )—D Cos( C+V)：“X=”：Ans▲Z[26] —(Z[27] —Z[24]) Sin( C —D Sin( C+V)：“Y=”：Ans▲：Ifend』Goto 1子程序名JD （输入曲线参数，计算曲线要素）30→Dim Z：Fix 5：“JD”? →Z[1] ：“XJD”? →Z[2] ：“YJD”? →Z[3]：“F0(DMS)”? →F：“A（-+ DMS）”？→Z[4]：“R”? →R：“LH1”? →Z[5]：“LH2”? →Z[6]：√￣(Z[4]2)→A：Int(Z[4] ÷A)→H』Z[5] →S：Prog “PQ”：Q →Z[7]：P→Z[8]：B →Z[9]：R →S：Prog “PQ”：Q →Z[10]：P→Z[11] ：B →Z[12] 』（R+ Z[8]）Tan(A÷2)+ Z[7] —(Z[8] —Z[11]) ÷Sin(A)→Z[13] ：“T1=”：Ans ▲（R+ Z[11]）Tan(A÷2)+ Z[10] —(Z[8] —Z[11])÷Sin(A)→Z[14] ：“T2=”：Ans▲』√￣((Z[13] —Z[7])2+( R+ Z[8])2) —R →Z[15]：“E=”：Ans▲R( A —Z[9] —Z[12] )÷57.2958+ Z[5] + Z[6] →Z[16]：“L=”：Ans▲』Z[1] —Z[13]→Z[17]：Z[2] —Z[13]Cos( F) →Z[18]：Z[3] —Z[13] Sin( F) →Z[19]：If Z[5]＞0 ：then Z[17] +Z[5]→Z[20]：Z[5] →L：Prog “HHXY”：Z[18] +I Cos( （F+ JH）)→Z[21]：Z[19] +I Sin( （F+ JH）)→Z[22]：（F+ WH）→Z[23] ：Else Z[17]→Z[20]：Z[18]→Z[21]：Z[19]→Z[22]：F→Z[23]：Ifend』Z[17] + Z[16]→Z[24]：Z[2] +Z[14] Cos( （F+ Z[4]）)→Z[25]：Z[3] +Z[14] Sin( （F+ Z[4]）)→Z[26]：F+ Z[4] —180→C』Prog “CALXY”』子程序名PQS÷2—S3÷（240 R2）→Q：S2÷（24 R）—S4÷（2384 R3）→P：28.6479S ÷R→B』子程序名HHXYL—L5÷(40S2 R2)+L9÷(3456 R4S4) →X：L3÷（6S R）—L7÷（336S3 R3）+L11÷（42240 R5S5）→Y：28.6479L2÷(S R) →W：Pol（X,Y）』子程序名YUXY57.2958L÷R →W：RSin(W)→X：R（1—Cos(W）)→Y：Pol（X,Y）』CASIO fx—5800P计算器路线坐标计算程序（单个交点）使用说明首先运行子程序“JD”，输入曲线参数，计算曲线要素。

一．使用说明：1.规定：（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。

（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。

（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。

（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。

（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。

2.输入与显示说明：（1）输入部分：N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐标计算里程边距。

X0？线元起点X坐标。

Y0？线元起点Y坐标。

S0？线元起点里程。

F0？线元起点切线方位角。

LS？线元长度。

R0？线元起点曲率半径。

RN?线元终点曲率半径。

Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）S?正算时所求点里程。

Z？正算时所求点距中线边距。

位于线路左侧输入负值，在右侧输入正值。

AR?斜交角度。

（为线路右角，在0-180度之间）X?反算时所求点X坐标。

Y？反算时所求点Y坐标。

（2）显示部分：XS=计算结果，所求点X坐标。

YS=计算结果，所求点Y坐标。

FS=计算结果，所求点切线方位角。

S=计算结果，反算时所求点里程Z=计算结果，反算时所求点边距。

二．程序清单1（用于任意曲线元的坐标计算）1.主程序(TYQXJS)″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 22.正算子程序(SUB1)0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) →X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y:G+QEW(C+WD) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(SUB2)G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lbl 0:Prog ″SUB1″:T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z:IF Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。

卡西欧CASIO fx-5800 道路坐标放样程序源程序包括：一个主程序和两个子程序主程序：PRO-ROAD行号程序说明1 Deg：Fix3：26→DimZ2 “X（JD1）”?X：“Y(JD1)”?Y：“X（JD0）”?M：“Y(JD0)”?N3 Prog”SUB-QXYS1”4 Pol(X-M,Y-N)：J→Z[20]：J+O→Z[19]5 X-Scos(Z[20]) →Z[21]：Y-Ssin(Z[20]) →Z[22]6 X+Tcos(Z[19]) →Z[23]：Y+Tsin(Z[19]) →Z[24]7 “X0”?U：”Y0”?V8 Lbl 09 “KP”?P：IFP>Z[1]：Then Goto 1 ：IfEnd10 Z[1]-P→L11 “XP=“：X-(S+L)cos(Z[20]) →F12 “YP=“：Y-(S+L)sin(Z[20]) →G13 Z[20] →Z ：Goto 514 Lbl 115 If P>Z[2] ：Then Goto 2 ：IfEnd16 P-Z[1] →L ：L→Z[14]：B→Z[15]：Prog”SUB-HXY”17 “XP=“：Z[21]+Z[16]cos(Z[20]-WZ[17]sin(Z[20]) →F18 “YP=“：Z[22]+Z[16]cos(Z[20]-WZ[17]sin(Z[20]) →G19 Z[20]+90WLL÷(BRπ) →Z20 Goto 521 Lbl 222 If P>Z[4] ：Then Goto 3 ：IfEnd23 P-Z[1] →L ：90(2L-B)÷R÷π→Z[13]24 Rsin(Z[13])+Z[8] →Z[16]：R(1-cos(Z[13]))+Z[6] →Z[17]25 “XP=“：Z[21]+Z[16]cos(Z[20])-WZ[17]sin(Z[20]) →F26 “YP=“：Z[22]+Z[16]sin(Z[20])-WZ[17]cos(Z[20]) →G27 Z[20]+WZ[13] →Z28 Goto 529 Lbl 330 If P>Z[5]：Then Goto 4：IfEnd31 Z[5]-P→L：L→Z[14]：C→Z[15]：Prog”SUB-HXY”32 “XP=“：Z[23]-Z[16]cos(Z[19])－WZ[17]sin(Z[19]) →F33 “YP=“：Z[24]-Z[16]sin(Z[19])＋WZ[17]cos(Z[19]) →G34 Z[19]－90WLL÷(CRπ) →Z35 Goto 536 Lbl 437 P－Z[5] →L38 “XP=“：X＋(T+L)cos(Z[19]) →F39 “YP=“：Y＋(T+L)sin(Z[19]) →G子程序1：SUBQXYS1子程序2：SUB-HXY40 Z[19] →Z41 Lbl 542 If Z<0 ：Then Z+360→Z：IfEnd43 “BP=“：Z DMS44 If U=0：Then Goto 7：Else Pol(F-U,G-V)：IfEnd45 Lbl 646 If J<0：Then J+360→J：IfEnd47 “A1=“：J DMS48 “D1=“：I49 Lbl 750 “D ANGLE”?H：IfH=0：Then Goto 0 ：IfEnd ：?D51 “XB=“：F+Dcos(Z+H) →Z[11]52 “YB=“：G+Dsin(Z+H) →Z[12]53 If U=0：Then Goto 7：Else Pol(Z[11]-U,Z[12]-V)：IfEnd54 Goto 6行号程序说明1 “ANGLE”?O：?R：”LS1”?B：”LS2”?C：”K(JD1)”?K2 If O<0：Then -1→W：Else 1→W：IfEnd：WO→A3 B2÷24÷R－B^(4)÷2688÷R^(3)＋B^(6)÷506880÷R^(5) →Z[6]4 C2÷24÷R－C^(4)÷2688÷R^(3)＋B^(6)÷506880÷R^(5) →Z[7]5 B÷2－B^(3)÷240÷R2＋B^(5)÷34560÷R^(4) →Z[8]6 C÷2－C^(3)÷240÷R2＋C^(5)÷34560÷R^(4) →Z[9]7 “T1=“：Z[8]＋(R＋Z[7]－(R＋Z[6])cos(A))÷sin(A) →S8 “T2=“：Z[9]＋(R＋Z[6]－(R＋Z[7])cos(A))÷sin(A) →T9 “L=“：RAπ÷180＋(B＋C)÷2→L10 “LY=“：L－B－C→Q11 “E=“：(R＋(Z[6]＋Z[7])÷2)÷cos(A÷2)－R→E12 “K(ZH)=“：K－S→Z[1]13 “K(HY)=“：Z[1]＋B→Z[2]14 “K(QZ)=“：Z[1]＋L÷2＋(B－C)÷4→Z[3]15 “K(YH)=“：Z[1]＋L－C→Z[4]16 “K(HZ)=“：Z[4]＋C→Z[5]行号程序说明1 If Z[14]=0：Then 0→Z[16]：0→Z[17]：Else2 Z[14]－Z[14]^(5)÷40÷(RZ[15])2＋Z[14]^(9)÷3456÷(RZ[15])^(4) →Z[16]3 Z[14]^(3)÷6÷(RZ[15])－Z[14]^(7)÷336÷(RZ[15])^(3)＋Z[14]^(11) ÷42240÷(RZ[15])^(5) →Z[17]4 IfEnd。