卡西欧5800程序：大地坐标转施工坐标(续)

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；右侧输入“+”值，左侧输入“-”值。

大地坐标与施工坐标相互转换（CASIOfx-5800P）刚刚接触5800计算器，看到大地坐标与施工坐标相互转换是两个分开的。

现在为了学习一下5800程序输入，把它组合成一个程序文件，请大家试用指点。

程序清单:Fix 4:”ALFA=”?→J: ”AXO=”? →C: ”BYO=”? →D↙Lbl 3:”1→AG,2→GA”?Q↙Q=1＝＞Goto 1: Q=2＝＞Goto 2↙Lbl 1:”AA=”?→A:”BA=”?→B: ↙”XG=”:C+ACOS(J)-BSIN(J)▲”YG=”:D+ASIN (J)+BCOS (J)▲Goto 3↙Lbl 2:”XG=”?→X:”YG=”?→Y: ↙”AA=”( X-C) COS(J)+(Y-D)SIN(J)▲”BA=”: ( C-X) SIN (J)+(Y-D)COS (J)▲Goto 3↙说明：ALFA＝大地坐标与施工坐标夹角。

AXO，BYO＝施工坐标系的原点的大地坐标AA，BA＝施工坐标X,Y＝大地坐标例：运行程序显示输入ALFA＝? 大地坐标与施工坐标夹角13'0'0'AXO＝？施工坐标系的原点的大地坐标X2870340.477BYO＝？施工坐标系的原点的大地坐标Y610780.2221→AG,2→GA?1为施工转换大地，2为大地转换施工1AA=?要转换成大地坐标点的施工坐标A54.9BA=?要转换成大地坐标点的施工坐标B82.5计算结果XG=2870375.411YG=610872.9573回车1→AG,2→GA?1为施工转换大地，2为大地转换施工2XG=? 要转换成施工坐标点的大地坐标X 2870375.411YG=? 要转换成施工坐标点的大地坐标Y 610872.9573计算结果A=54.8995B=82.5001（欢迎参与右下角评论）公办中专（中央广播电视中等专业学校）常年滚动招生，滚动注册：专业齐全（中专），建筑类专业（可考二级建造师、安全员、造价员必备学历、征兵入伍可用）药剂类专业（可考执业药师）（注：费用低、出证快，教育部直属中专，可做学历认证，国家学信网永久查询。

公路曲线大地坐标施工放样程序（CASIO fx-58000）一、主程序28:""?:""?:""?:""?:.5^(3)/(240):""?:"":(/(24R )-G ^(4)/(2688R ^(3)))ta n (0.5)T "":/180+G L""?:"":[3]"1,2"?:""?:""?:""?:10:"Z B J SD im Z X C Z Y C C L S G R R G G R Q A A T Q R G L R A Z H O H Z O Z Y Z H J D N P J DE MF F L b K π→====-→==++→=→==+→=====▲▲▲"?:","?:""?:0:0:0:P r "":P r "":0:K V N W V B B K O S IfS T h e n S XY U o g C o g D G o to IfE n d ==-+=-→↵≤→→→↵[][]IfS L :T h e n T +C O S (A )(T +S -L )X :S i n (A )(T +S -L )Y :A U :P ro g "C ":P ro g "D ":G o to 0:IfE n d IfS L -G :T h e n L -S S :P ro g "B ":X Z1:YZ2:T -Z [2]S in (A )+C o s (A )(T -Z [1])X :Z [2]C o s (A )+S in (A )(T -Z [1])Y :L -S S :A -U U :P ro g "C ":P ro g "D ":G o ≥→→→↵≥→→→→→→→2222to 0:IfE n d IfS G :T h e n 90(2S -G )/(R )U :Q +R S in (U )X :G/(24R )-G ^(4)/(2688R ^(3))+R (1-C o s (U ))Y :P ro g "C ":P ro g "D ":G o to :If E n d If S G :T h e n K -0S :S -S ^(5)/(40RG)X :S ^(3)/(6R G )-S ^(7)/(336R ^(3)G ^(3))Y :90S /(R G )ππ↵≥→→→↵≤→→→U :P ro g "C ":P ro g "D ":G o to 0:I f E n d →↵二、子程序 BSC"X=Dta注：”/”表示÷三、程序变量说明以程序运行时计算器显示为准：XC?EC?提示输入架仪点的坐标LS?提示输入测设曲线缓和曲线长R?提示输入测设曲线半径A?提示输入测设曲线转角（以度分秒输入）T= L=显示计算的曲线元素（检验LS、R、A是否输入正确）ZH? 提示输入测设曲线直缓（圆）桩桩号（对比设计资料可检查直缓（直圆）桩桩号输入是否正确Y1,Z2?判断曲线转向，右转输入1，左转输入2JDN?JDE?提示输入交点的N和E坐标F?提示输入曲线计算起始边的坐标方位角（以度分秒格式输入）K?提示输入欲放样点的桩号V?提示输入欲放样点相对中点的法线方向的距离，位于曲线内侧输入负值，外侧为正值，测设中桩时输0B?提示输入欲放样点相对中点的切线方向的距离，位于曲线前进方向输入正值，后退方向输入负值。

主程序坐标正反算程序名称：001XYLB1 A↙Cls : Fix 4 : 30→Dim Z :"v=1,v=0"?Z ↙If Z=1:Then Goto B:IfEnd↙(Z=1进入里程点坐标正算）If Z=0:Then Prog"FSLCBZ ": Goto A:IfEnd↙(Z=0进入反算里程边桩）Lb1 B ↙“K=”?K ：（计算里程）Prog"SJK1 ":Prog”ZBFY” :Goto A↙子程序数据库 SJK1IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN 本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：1→O（注：左偏曲线输入-1→O,右偏曲线输入1→O）: 半径→R : 曲线偏角→A:第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Prog”XLZBJSCX”:Return: IFEND↙IF…………Prog”XLZBJSCX”:Return:IFEND（曲线段分段输入）↙补充直线段输入如下(单独直线段）IF K<本段直线终点里程AND K≥本段直线起点里程:THEN 本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E: Z[16]+ (K- Z[3])*COS E→Z[18]:Z[17]+ (K- Z[3])*SIN E→Z[19] :E →T:Return:IFEND子程序坐标正算，名称：XLZBJSCXLB1 2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙（M1）Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙（M2）Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙（P1）Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙（P2）π*A*R/180+0.5*( Z[6]+ Z[7])→Z[25] ↙（曲线总长）90* Z[6]/(R*π) →Z[14] ↙（第一缓和曲线总偏角）90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略）Z[8]＋（R+Z[10])TAN(A/2)-(Z[10]-Z[11] )/SIN A→Z[12]↙(切线T1)Z[9]＋（R+Z[11])TAN(A/2)+(Z[10]-Z[11] )/SIN A→Z[13]↙(切线T2)B+ Z[12]*COS (E+180)→ Z[16] ↙（ZH点X）C+ Z[12]*SIN(E+180)→ Z[17] ↙（ZH点Y）Z[1]- Z[25]→Z[3] ↙(ZH点里程)Z[3]+ Z[6]→Z[4] ↙(HY点里程)Z[1]- Z[7]→Z[5] ↙(YH点里程)GOTO 3 ↙LB1 3 ↙(判断里程点与曲线关系)IF K≤Z[3] AND K> Z[2] : THEN GOTO 4 : IFEND ↙IF K≤Z[4] AND K> Z[3] : THEN GOTO 5 : IFEND ↙IF K≤Z[5] AND K> Z[4] : THEN GOTO 6 : IFEND ↙IF K≤Z[1] AND K> Z[5] : THEN GOTO 7 : IFEND ↙LB1 4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →Z[23] : 0→Z[24] : E→T : GOTO 8↙LB1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →Z[23]↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Z[24] ↙90*H^2/( R*π* Z[6]) →T ↙IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 8 ↙LB1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*SIN T+ Z[8]→Z[23] ↙R*(1-COS T)+ Z[10]→Z[24] ↙IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 8 ↙LB1 7 ↙（第二缓和曲线独立坐标）Z[1] -K →H ↙H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]COS A+ Z[12]-U*COS A-V*S IN A→Z[23] ↙Z[13]*SIN A-U*SIN A+V*COS A→Z[24]↙IF O>0 ：THEN F-T→T : T<0=>360+T→T : ELSE F+T →T : IFEND ↙GOTO 8 ↙LLb1 8IF O<0 : THEN - Z[24]→Z[24] : IFEND ↙Z[16]+Z[23]*COS E-Z[24]*SIN E→Z[18] ↙Z[17]+Z[23]*SIN E+Z[24]*COS E→Z[19] ↙ReTurn↙子程序反算里程边桩，名称：FSLCBZLb1 1↙"POINT=X"?X:" POINT=Y"?Y↙（输入任意测点的XY坐标）:“K=”?K （试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）↙Lb1 2↙Prog"SJK1 ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin W)→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK1":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else IF Q<S: THEN K+Q→K: Goto 2 :Else IF Q>S :THEN K-Q→K:Goto 2:Ifend:Ifend:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"W(-Z,Y+)=":I◢(偏距）（由于该程序不能准确判断边桩左右方向，暂作修改，取消左右边判定）"K=":K+S→K◢（里程）Return↙子程序坐标放样：ZBFYLB1 0 ↙"P(N)="?M:"P(E)="?N ↙“a=”?P◢(输入边桩与线路夹角，左-右+)“W=”?D◢（输入边桩距）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢（边桩X）“YP=”: Z[21] ◢（边桩Y）Pol(Z[20]-M,Z[21]-N):"HD=":I◢（边桩放样距）J<0=>J+360→J:"HAR=":J→DMS◢（边桩放样方位角）Return↙直线短链要将该直线分成二条直线来对待编制数据库文件即可，曲线长链要单独建判定语句子程序并然后编制二个数据库（按来向里程和去向里程分别推出曲线起终点里程，以此编制曲线段的数据库）。

FILE非常好用的卡西欧5800p程序一、坐标反算(XY--FD)Deg：Fix 3 ↓(换行)Lbl 0↓〝X0=〞?A:〝Y0=〞?B:〝XN=〞?C:〝YN=〞?D↓(换行)Pol（（C-A）,(D-B)）J＜0＝＞J+360→J↓(换行)〝JL=〞：I◢〝FWJ=〞：J：▶DMSGoto 0二、坐标正算(DF--XY)Deg：Fix 3 ↓(换行)〝X0=〞?A；〝Y0=〞?B：〝L=〞?L：〝C=〞?C↓(换行)A+Lcos(C)→X↓B+Lsin(C)→Y↓〝X=〞：X◢〝Y=〞：Y◢三、直线段中桩、边桩坐标计算（该程序从起点往终点计算）(ZX)1、Deg：Fix 3 ↓(换行) 设置小数保留位数2、Lbl 1↓程序标签3、〝X0〞?X：〝Y0=〞?Y：〝K0=〞?E：〝FWJ=〞?J↓赋值起点坐标、桩号、走向方位角4、?K：?L：?A↓输入所求的桩号、边距、相对于主线的方位角（如果只计算中桩坐标，边距及方位角皆为0）5、X+(K-E)cos(J)+Lcos(J+A)→C↓计算坐标X6〝XP=〞：C◢显示坐标X7、Y+(K-E)sin(J)+Lsin(J+A)→D↓计算坐标Y8、〝YP=〞：D◢显示坐标Y9、Goto 1↓转向程序标签循环执行计算四、圆曲线坐标计算程序（该程序从起点往终点计算）(YQX)1、Deg：Fix 3 ↓(换行)2、Lbl 0：3、〝K0=〞?K：?X：?Y：?A：?R：4、〝L=-1，R=+1〞：?G：5、〝ZS=0，FS≠0〞：?S：6、S≠0 =＞Goto 1↓7、〝Ki〞?Z↓8、90(Z-K)÷(∏R)→B↓9、2RsinB→D↓10、〝Xi=〞：X+Dcos（A+GB）→P◢11、〝Yi=〞：Y+Dsin（A+GB）→Q◢12、A+180G(Z-K)÷(∏R)→0↓13、〝ZBJ〞？E↓14、〝XiL=〞:P+Ecos(0-90) ◢15、〝YiL=〞:Q+Esin(0-90) ◢16. 〝YBJ〞？F↓17、“XiR=”:P+Fcos(O+90)◢18、“YiR=”:Q+Fsin(O+90)◢19、Goto 0↓20.、Lbl 1：21、“XC”?U: “YC”?V↓22、X+Rcos(A+90G)→M↓23、Y+Rsin(A+90G)→N↓24、Pol(U-M,V-N)→H↓25、J＜0 =＞J+360→J↓26.、“PJ=”:R-H→L◢27、“XCCZ=”: M+Rcos(J)→T◢28.、“YCCZ=”: N+Rsin(J)→W◢29 、Pol(T-X,W-Y) →Z↓30.、“CZZH=”: K+sin-1(Z÷2÷R)πR÷90◢31.、Goto 0说明：1.设置小数保留三位。

线路测量坐标正、反算计算原理及卡西欧fx-5800P程序说明一、计算原理在各测量书中对于坐标正算的相关计算式均有说明，故在此不做详解，仅对正算过程中需要用到的原理及公式做一汇总。

对于坐标反算，虽然都采用无限趋近原理进行计算，但计算方式各有差别，本文仅对其中一种自认为相对简单易懂并便于操作的原理进行详解。

1.1 坐标转换[1]如图1，设X P、Y P为P点在国家控制网坐标系中的坐标；x P、y P为P点在工程独立控制网坐标系中的坐标。

X O、Y O为工程独立坐标系原点o在国家坐标系中的坐标，Δα为两坐标系纵坐标轴的交角。

如果一条边在国家坐标系中的坐标方位角为A，而在工程独立坐标系中的坐标方位角为α，则：∆α=A−α(1-1)当由工程独立坐标系坐标换算至国家坐标系坐标时，换算公式为：X=x cos∆α−y sin∆α+X O(1-2)Y=x sin∆α+y cos∆α+Y O}当由国家坐标系坐标换算至工程独立坐标系坐标时，也可使用式(1-2)，此时应将X、Y与x、y互换，且∆α=α−A。

1.2 坐标方位角关系计算1.2.1 正、反坐标方位角[2]一条直线的坐标方位角与直线的前进方向有关，沿直线前进方向的坐标方位角称为正坐标方位角，与其相反方向的坐标方位角称为反坐标方位角。

如图2，由于轴子午线之间是互相平行的，因此同一直线的正、反坐标方位角相差180°，即：α正=α反±180°(2-1)当α反<180°时，取“+”号；当α反>180°时，取“-”号。

1.2.2 坐标方位角的推算[3] 1.2.2.1 转折角为右角如图3(a)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β右为该两边所夹的右角，则：α23=α12±180°−β右=α21−β右 (2-2)1.2.2.2 转折角为左角如图3(b)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β左为该两边所夹的右角，则：α23=α12+β左±180°=α21+β左(2-3) 无论用右角还是左角推算，如遇出现负数的情形，应加上360°。

CASIO5800计算器编程在河口村水库高边坡开挖施工测量放样中的应用摘要:河口村水库大坝趾板高边坡岩石开挖最大高差约为100米，施工放样工作非常艰巨，为保证高边坡开挖后趾板体型，高边坡开挖需要很高的精度。

测量放线实施中CAD及可编程计算器的结合运用，既减少了计算，又提高了测量精度。

测量数据通过EXCEL处理，在CAD上的成图，又效验了测量成果的准确性。

本文就结合河口村水库大坝标的趾板边坡开挖测量实践，对在高边坡开挖施工的CASIO5800计算器编程及测量数据校核做一个简要分析。

关键词:测量放样CASIO5800计算器编程1、工程概况河口村水库混凝土面板堆石坝最大坝高122.5m，坝顶高程288.5m，防浪墙高1.2m，坝顶长度530.0m，坝顶宽9.0m，上游坝坡1：1.5，下游坝坡1：1.5。

大坝趾板从X1～X12，一共为10段，其中X1～X10的趾板边坡均为岩石边坡，需要进行爆破开挖。

垂直于趾板”X”线馒头组岩层为1:0.75，其余岩层开挖为1:0.5，马道宽2.0m。

以右岸趾板X2～X3段为例，边坡顶部高程为EL263，开挖至底部EL165高程，从EL240开始至EL180，每隔20米设2米宽马道，EL220以上部分岩石开挖边坡为1:0.75，EL220高程以下岩石开挖边坡为1:0.5。

图1右岸趾板开挖平面图表1设计提供趾板控制点X线控制点坐标（左右坝肩开挖蓝图提供数据）控制点 Y坐标X坐标高程趾板桩号趾板坡降备注X1 377147.2090 3897387.6940 286.0000 0X2 377206.4590 3897297.3950 226.5000 108.0022 1.8152 X1-X2 X3 377263.5880 3897198.0380 166.0000 222.6126 1.8944 X2-X32、编程放样在了解施工内容及图纸之后，为便于我们能够通过数据直观的了解地物间相对位置，同时减少放样过程中的计算器输入，我们首先根据设计提供的测量控制点建立施工坐标系，建立以趾板“X”线为轴线的趾板施工坐标系，在以后的趾板开挖及浇筑施工过程中都将使用趾板施工坐标系。

四大主点的坐标和方位角计算程序：LU-XIANClrstat:20→Dimz:Fix6:Lib0:"NO.J(2)="?S:"R="?R:"L(H)="?L:Prog"FANG- WEI":J+180→J:J-O→K:K<0＝＞K+360→K:K-180→F:If>0:ThenL(1):ElseR(1): IfEnd←L(1)＝＞-1→W:R(1)＝＞1→W:J+W(I*180/(2πR))→z[6]:"B(HY)=":(z[6]-180)▲DME◢I*180/(2πR)→B:I/2-I^(3)/(240*R^(2)→G:I^(2)/(24R)→P:"T=":(R+P)*tan(Abs(F)/2)+G→z[7]→◢ "L(Y)=":πR( Abs(F)-2B)/180→z[8]◢ "NO.ZH=":S-z[7]→z[9]◢ "NO.HZ=":z[9]+2I+z[8]◢ "X-ZH=":C+z[7]*cos(J)→z[1]◢ "Y-ZH=":D+z[7]*sin(J)→z[2]◢ "X-HZ=":C+z[7]*cos(O)→z[3]◢ "Y-ZH=":D+z[7]*sin(O)→z[4]◢I-I^(3)/(40R^(2))+I^(5)/(3456*R^(4))-I^(7)/(599040R^(6)→z[5]:"X-HY=":z[1]+z[5]/cos(30I/(πR))*cos(J-180+30WI/(πR))◢"Y-HY=":z[2]+z[5]/cos(30I/(Rπ))*sin(J-180+30WI/(πR))◢"B(Y)=":O+(-W)*(180I/(2πR)◢ "X-YH=":z[3]+z[5]/cos(30I/(πR))*cos(O+180-30WI/(πR))◢ "Y-YH=":z[4]+z[5]/cos(30I/(πR))*sin(O+180-30WI/(πR))◢ Goto 0 ←NO.JD(2),R,L(H)分别为：交点里程，半径，缓和曲线长;T,L(Y),NO.ZH,NO.HZ, X-ZH,Y-ZH,X-HZ,Y-HZ,X-HY,Y-HY,X-YH,Y-YH,B(YH) 分别为:切线长，圆曲线长，直缓点里程，缓直点里程，四大主点的坐标和圆缓点方位角方位角计算程序：FANG-WEIClrstat:Deg:"X1="?A:"Y1="?B:"X2="?C:"Y2="?D:"X3="?E:"Y3="?F←C-A→X:D-B→Y:Arg(X+Yi)→J:J<0＝＞J+360→J:E-C→M:F-D→U:Arg(M+Ui)→O:O<0＝＞O+360→O:"J=":J▲DME◢←"O=":O▲DME◢X1,Y1,X2,Y2,X3,Y3:分别为三个交点的坐标。