5800曲线任意里程中边桩坐标正反算

CASIO fx-5800P线元法坐标正反算程序V2.0说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。

该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。

主程序名：ABCYT第1步Deg：Fix 3：10→DimZ第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂K0+＂?S：Prog＂ABCYTSJ＂：If P=0：Then 10^(45)→P：IfEnd：If R=0：Then 10^(45)→R：IfEnd第3步1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If W=1：Then Goto 1：Else Goto2：IfEnd第4步Lbl 1：＂E＂?Z90=N：Abs(S-O)→W：Prog ＂ABCYTZ＂第5步Cls：＂F=＂：Locate 3,1,F°：＂X=＂：Locate 3,2,X：＂Y=＂：Locate 3,3,Y◢第6步Prog＂ABCYTSQX＂：Cls：＂H=＂：Locate 3,1,H◢第7步1→W：90→N：Goto 3第8步Lbl 2：?X：?Y：X→I：Y→J：Prog＂ABCYTF＂：O+W→S第9步Cls：＂K=＂：Locate 3,1,S：＂E=＂：Locate 3,2,Z◢第10步2→W：Goto 3正算子程序名：ABCYTZ第1步0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L第2步1-L→F：1-K→M第3步U+W×(A×cos(G+Q×E×K×W×(C+K×W×D))+B×cos(G+Q×E×L×W×(C+L×W×D))+B×cos(G+Q×E×F×W×(C+F×W×D))+A×cos(G+Q×E×M×W×(C+M×W×D)))→X第4步V+W×(A×sin(G+Q×E×K×W×(C+K×W×D))+B×sin(G+Q×E×L×W×(C+L×W×D))+B×sin(G+Q×E×F×W×(C+F×W×D))+A×sin(G+Q×E×M×W×(C+M×W×D)))→Y第5步G+Q×E×W×(C+W×D)→F：F+ N→Z[1]第6步X+Z×cos（Z[1]）→X：Y+Z×sin（Z[1]）→Y反算子程序名：ABCYTF第1步Lbl 2：（S-O）→W：0→Z：Prog ＂ABCYTZ＂：F-90→Z[9]：(J-Y)×cos(Z[9])-(I-X)×sin(Z[9])→Z[10]第2步If Abs(Z[10])>0.001：Then S+Z[10]→S：Goto 2：Else Goto 1：IfEnd第3步Lbl 1：(Y-J)÷sin(Z[9])→Z数据库名：ABCYTSJ第1步Goto1第2步Lbl 1：If S<7586.707 Or S>13346.96：Then Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：IF S<7946.707：Then 98°56′56″→G：7586.707→o：3378605.445→U：453648.704→V：0→P：4500→R：360→H：1→Q：Return：IfEnd第5步Lbl 1：IF S<11766.03：Then 101°14′26″→G：7946.707→o：3378544.714→U：454003.518→V：4500→P：4500→R：3819.323→H：1→Q：Return：IfEnd第6步Lbl 1：IF S<12126.03：Then 149°52′11″→G：11766.03→o：3376389.890→U：457018.324→V：4500→P：0→R：360→H：1→Q：Return：IfEnd第7步Lbl 1：IF S<13346.96：Then 152°09′41.68″→G：12126.03→o：3376073.846→U：457190.654→V：0→P：0→R：1220.93→H：0→Q：Return：IfEnd第n步……………………………………………………数据输入说明：第1步Goto1第2步Lbl 1：If S<本条线路起点里程S>本条线路止点里程Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：If S<本线元止点里程：Then线元起点切线方位角→G：线元起点桩号→O：线元起点坐标X→U：线元起点坐标Y→V：线元起点半径（直线为0、曲线为半径）→P：线元止点半径（直线为0、曲线为半径）→R：线元长度→H：线元转向（左转为-1、右转为1、直线为0）→Q：Return：IfEnd第n步Lbl 1：每增加一行则为增加一个线元要素。

CASIO5800坐标正反算计算程序CASIO 9860/5800坐标正反算计算程序1. 坐标正算主程序(命名为ZBZS)第1行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B第2行：Prog “A”第3行：”X=”:N+Zcos(F+B)◢第4行：”Y=”:E+Zsin(F+B)◢第5行：”F=”:F◢第6行：Goto 0K——计算点的里程BIAN——计算点到中桩的距离（左负右正）α——取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为ZBFS)第1行：”X1=”? C:”Y1=”?D: ”K1=”?K第2行：Lbl 0:Prog “A”第3行：Pol(C-N,D-E)第4行：List Ans[1]→I第5行：List Ans[2]→J第6行：Icos(F-J)→S:K+S→K第7行：Abs(S)>0.0001=>Goto 0第8行：”K1=”:K◢第9行：”BIAN=”:Isin(J-F)→Z◢X1——取样点的X坐标Y1——取样点的Y坐标K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号Z——偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-A→S：(Q-P)÷L→I第2行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N第3行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E第4行：F+S(2P+SI)×90÷π→F4. 数据库（命名为A）第1行：K≤175.191=>Stop第2行：175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/ 240→Q:70.417→L:K≤A+L =>GoTo 1（第一缓和曲线，圆半径为240）第3行：245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q:72.915→L: K≤A+L =>Goto 1（第圆曲线，半径为240）第4行：318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104→L: K≤A+L =>Goto 1（第二缓和曲线，圆半径为240）第5行：373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:-1/180→Q:67.222→L:K≤A+L=>Goto 1:Stop（第一缓和曲线，圆半径为180）第6行：Lbl 1:Prog “XYF”A——曲线段起点的里程N——曲线段起点的x坐标E——曲线段起点的y坐标F——曲线段起点的坐标方位角P——曲线段起点的曲率(左负右正)Q——曲线段终点的曲率(左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）在9860中，程序中所有公式和部分函数结果均存储在ListAns列表数组中，要想多次调用最好随公式取出结果，并赋给变量。

卡西欧5800P坐标正反算程序超好用5800P坐标正算反算程序U：“XO”起点的X坐标，V：“YO”起点的Y坐标，O：“SO”起点的里程桩号，G：“FO”起点的方位角，H：“LS”计算段线元的长度，P：“RO”线元起点的曲率半径，R：“RN”线元终点曲率半径，“Q ”线元图形信息（直线段Q=0，左偏Q=-1，右偏Q=1）以道路中线前进方向（即大里程方向）区分左右。

DY=ZB:LC：输入1为正算，既输入里程桩号和偏距算坐标，输入2为反算，既输入坐标算里程桩号和偏距。

S:“LC”？S：待求点的里程桩号，Z:“PJ”？Z：待求点的里程偏距（在中线上取零，左侧取负值，右侧取正值）。

1．当线元为直线时，其起点，终点的曲率半径为无穷大，以10的45次方代替。

2．当线元为圆曲线时，无论其起点，终点与什么线元相接，其曲率半径均等于圆曲线的半径。

3．当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替，与圆曲线相接时，曲率半径等于圆曲线半径。

终点与直线相接时，曲率半径为无穷大，以10的45次方代替，与圆曲线相接时，曲率半径等于圆曲线半径。

4．当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值，与圆曲线相接时，曲率半径等于圆曲线半径，终点与直线相接时，曲率半径等于设计规定值，与圆曲线相接时，曲率半径等于圆曲线半径。

1．坐标正反算主程序A-ZBZS3U：“XO ”？U：V：“YO”？V：O：“SO”？O：G：“FO”？G：H：“LS”？H：P：“RO”？P：R：“RN”？R：“Q”？Q：(1÷P)->C:(P-R)÷(2HPR)->D:(180÷Π)->E:N：“DY=ZB:LC”?N:If N=1:Then Goto 1:Else Goto 2:IfEnd:Lbl1:S:“LC”？S：Z:“PJ”？Z：Abs（S-O）->W：Prog“SUB1”：“XS=”：X◢“YS=”：Y◢“FS=”：F-90◢Goto 1:Lbl 2:“X”?X: “Y”?Y:X->I:Y->J:Prog “SUB2”:S: “LC=”:O+W->S◢Z: “PJ=”:Z◢Goto 2:2．坐标正算子程序SUB10.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）+90->F：X+Zcos(F)->X：Y+Zsin（F）->Y：3.坐标反算子程序SUB2G-90->T:Abs((Y-V)cos(T)-(X-U)sin(T))->W:0->Z:Lbl0:Prog“SUB1”:T+QEW(C+WD)->L:(J-Y)cos(L)-(I-X)sin(L)->Z:If Abs(Z)<10^(-3):Then Goto 1:Else W+Z->W:Goto 0:IfEnd:Lbl 1:0->Z:Prog “SUB1”:(J-Y)÷sin(F)->Z:。

卡西欧5800P计算器曲线任意⾥程中边桩坐标正反算卡西欧5800P计算器曲线任意⾥程中边桩坐标正反算⼀、程序功能本程序由⼀个主程序(ZHUCHENXU)和⼏个⼦程序——正算⼦程序(SUB1)、反算⼦程序(SUB2)、数据库(SUB3)、算⽅位⾓程序(FA)、使⽤⽅位⾓算坐标(NE)、边坡放样程序(BIANPO)、隧道圆⼼放样程序（SDYX）构成，可以根据直线、圆曲线、缓和曲线（完整或⾮完整型）的线元要素（起点坐标、起点⾥程、起点切线⽅位⾓、线元长度、起点曲率半径、⽌点曲率半径）及⾥程边距或坐标，对该曲线段范围内任意⾥程中边桩坐标进⾏正反算。

⼆、源程序增加变量16→DimZ1.主程序(ZHUCHENXU)"1.SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO "？N：Lbl 1：”ZHUANGHAO=”?S：Prog "SUB3"：N≥2=>Goto 2：Abs(S-O)→W：”BIANZHU=”?→Z:Prog "SUB1"："XS=”:X→X ◢"YS=”:Y→Y◢"FS=”:(F-90) →F:F►DMS◢Prog "FA"：Goto 1:Lbl 2：”CX=”?X:X→I：”CY=”?Y: Y→J：”DMG=”?→Z[7]：Prog "SUB2"："S=":(O+W) →S◢"Z="Z→Z◢If N=3:ThenProg”SDYX”:IfEnd ：If N=4:Then Prog”BIANPO”:IfEnd:Goto 22. 正算⼦程序(SUB1)0.1739274226 →Z[1]：0.3260725774→B：0.0694318442→K：0.3300094782→L：(1-L) →F：(1-K)→M：U+W(Z[1]cos(G+QEK W(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Z[1]cos(G+QEMW(C+MWD)))→X：V+W(Z[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+ZcosF→X：Y+ZsinF→Y3. 反算⼦程序(SUB2)G-90 →Z[2]：Abs((Y-V)cosZ[2]-(X-U)sinZ[2]) →W：0→Z：Lbl 0：Prog "SUB1"：Z[2]+QEW(C+WD) →L：(J-Y)cosL-(I-X)sinL→Z：If Abs(Z)<1E-6:Then Goto1：Else W+Z→W：Goto 0: IfEndLbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z4. 数据库(SUB3)If S<226100:Then Cls:Stop: IfEndIf S<226255.833:Then 226100→O：99037.9736→U：4282.3590→V：201.3266304→G：1045→P：1045→R：253.543→H：0→Q: Goto 1:IfEnd…………If S>236005.651:Then Cls:Stop: IfEndLbl 1: 1÷P→C：(P-R)÷(2HPR) →D：180÷π→E：5. 算⽅位⾓程序(FA)LbI0：“XC=”? A : “YC=”? T :X：Y0→J：POI（X-A，Y-T）→I：“BC=”：I◢If J＜0 Then 360+J→J：ElseJ→J :IfEnd“FW=”：J►DMS◢ Goto06. 使⽤⽅位⾓算坐标(NE)LbI0：“X=”? X: “Y=”? Y: “L=”? L:“F0”?F：L×Cos（F）+X→N: “N=”N◢L×Sin（F）+Y→E: “E=”E◢Goto07. 边坡放样程序(BIANPO)“QPDJL=”? → Z[3]: “QPDGD=”? → Z[4]: “PG=” ? → Z[5]：2→Z[6]：“P1=” ? → Z[8]：“P2=” ? → Z[9]：“P3=” ? → Z[10] (没有的话就删)：Z[7]-Z[4] → Z[11]：If Z[11]〈0 Then AbsZ[11] →Z[11]: “T=”：Z[11] ◢Else“W=”：Z[11] ◢IfEndZ[11]÷Z[5] →Z[12]：If Z[12]〈1 Then Z[8]×Z[11] →Z[13]：Else Z[8]×Z[5]+Z[9]×Z[5]×（Z[12]-1）→Z[13]：IfEnd:AbsZ-（Z[13]+Z[6]×Int Z[12]+Z[3]）→Z[14]:“BP+C-J=”: Z[14] ◢说明：Z[3]= QPDJL ………………起坡点离中桩距离Z[4]=“QPDGD” ……………起坡点设计标⾼Z[5]=“PG” ………………坡⾼Z[6]=“PTK” ……………… 平台宽度Z[7]=“DMG” ………………实测地⾯⾼度Z[8]=“P1” ………………⼀级边坡坡度Z[9]=“P2” ………………⼆级边坡坡度Z[10]=“P3” ……………… 三级边坡坡度(没有的话就删)Z[14]=“BP+C-J=” ………………边坡超⽋挖宽度（+表超，-表⽋）8. 隧道圆⼼放样程序（SDYX）IF S<229000:THEN 243.997+0.0205×(S-226520)-(S-227318)^2/2/16000+1.503 →Z[10]:“NGC=”：Z[10] ◢ELSE 286.385-0.023×(S-230500)- (S-230593.753)^2/2/165000+1.503→Z[10]:“YLJ=”：Z[10] ◢IfEndZ[10]- Z[7] →Z[15]:“+SG,-JD”：Z[15] ◢√（Z[15]^2+(ABS(Z)-6.8)^2）→Z[16]: “R=”:Z[16] ◢三、使⽤说明1、规定(1) 以道路中线的前进⽅向（即⾥程增⼤的⽅向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。