万能全线坐标计算程序

合集下载

5800---万能程序

万能程序包含所有曲线要素的万能测量程序（6种模式）扩展变量操作（15个）：Defm 15 , O字母0数字←┚J-SQX（竖曲线数据输入）“J-SQX”:{NHUVMQP}:Z[9]=N“SJD”:Z[10]=H“JDZ”:Z[11]=U“I1”:Z[12]=V“I2”:Z[13]=M“R”:Z[14]=Q“QD”:Z[15]=P“ZD”:“T O J-JSMS”J-PQX (平面数据输入，自行切换到J-JSMS)A“JD” B“JDX” C“JDY” F“FJ” O“ZJ：Z-，Y+” RE“LS1”K“LS2”：E<1=>E=1E-9⊿K<1=>K=1E-9⊿Z[1]=EE÷24R-E∧4÷2688RRR:Z[2]=E÷2-EEE÷240RR:X=(EE-KK)÷24R÷sin Abs O :Z[3]“T1”=(R+Z[1])tan(Abs O÷2)+Z[2]-X◢Z[4]“T2” =(R+KK÷24R-K∧4÷2688RRR)tan(Abs O÷2)+K÷2-KKK÷240RR+X◢L=Abs OπR÷180+(E+K)÷2◢J=tan-1((R+Z[1])÷(Z[3]-Z[2]):X“E”=(R+Z[1])÷sin J-R◢X=A-Z[3]:Y=X+E:E<1=>X“ZY”◢≠=>X“ZH”◢Y“HY”◢⊿Y“QZ”=X+(L-K-E)÷2+E◢Y=X+L-K:X=X+L:K<1=>X“YZ”◢≠=>Y“YH”◢X“HZ”◢⊿Prog“J-JSMS” ←┚J-JSMS（放样模式主程序）“1-ZS,2-F1,3-FS”:“4-F2,5-Z0,6-DM”:Lb1 0:{Z}:Z“MS”≤2=>Goto 1: ≠=> Z≤4=> Goto2:≠=> Goto 3⊿⊿←┚Lb1 1 :{PDW}:PD“BZ”W“BJ”: Prog“JP”:X“X=”◢(4850输入“X=”：X◢)Y“Y=” ◢(4850输入“Y=”：Y◢)Z<2=> Goto 1:≠=>Prog“JS”:Prog“JC”:H“Z”=H+G◢Goto 1⊿←┚Lb1 2 :{XYW}:XYW“BJ”: Prog“JF”: P“P=”◢(4850输入“P=”：P◢)D“BZ=”◢(4850输入“BZ=”：D◢)Z<4=>Goto2:≠=>{Q}:Q“ZP”:Prog“JS”:Prog“JC”:Prog“JB”:Goto2⊿←┚Lb1 3 :{P}:P: Prog“JS”:H“Z0=”◢(4850输入“Z0=”：H◢)Z<6=>Goto3:≠=>{D}:S“CGZ”:T“JKZ”:Z[8]=Q“LH”:D“BZ”:Prog“JC”:J “JK=”◢(4850输入“JK=”：J◢)G“Z”=H+G◢(4850输入“Z=”：G=H+G◢)Goto 0⊿JS（竖曲线计算子程序）P<Z[14]=>Prog“J-S QX” ⊿P>Z[15]=>Prog“SJ-SQX”⊿←┚N=Z[9]:U=Z[11]:V=Z[12]:H=Z[13]:G=HAbs(U-V)÷200:P<N=>I=U:M=N-G:M>P=>M=P⊿≠=>I=V:M=N+G:M<P=>M=P⊿⊿J=(P-M)2÷2H:U-V>0=>J=-J⊿H=Z[10]+(P-N)×I ÷100+JJP（平曲线正算子程序）Fixm←┚Lb1 1 :J=F:X= B-Z[3]cos F :Y= C-Z[3]sin F :G=E←┚P≤A-Z[3]=>I=A-Z[3]-P:M=-I:N=0:H=F+W:Goto 5←┚≠=>P≤A-Z[3]+E=>I=P-A+Z[3]:H=90II÷REπ:O<0=>H=-H⊿H=H+W+F:G oto 3←┚≠=>P≤A-Z[3]+L-K=>I=P-A+Z[3]-E:Goto 4: ≠=>Goto 2⊿⊿←┚ Lb1 2 :X= B+Z[4]cos(F+O :Y= C+Z[4]sin(F+O :J=F+O+180:G=K←┚ P≤A-Z[3]+L=>I=A-Z[3]+L-P:H=90II÷RKπ:O>0=>H=-H⊿H=H+J+W+180 :Goto 3 ←┚≠=>I=P-A+Z[3]-L:M=-I:N=0:H=J+W+180:Goto 5 ←┚Lb1 3 :M=I-I∧5÷40RRGG:N=III÷6GR-I∧7÷336RRRGGG:Goto 5 ←┚ Lb1 4 :H=(E+2I)×90÷πR:M=R sin H+Z[2]:N=R(1-cosH)+Z[1]:O<0=>H=-H⊿H=J+H+W: Goto 5←┚Lb1 5 :P≤A-Z[3]+L-K=> O<0=>N=-N⊿≠=>O>0=>N=-N⊿⊿Goto 6←┚ Lb1 6 :X =X+M cos J-N sin J+D cos H:Y =Y+M sin J+N cos J+D sin HJF（平曲线反算子程序）Fixm←┚U=X:V=Y:D=0:J=F-W:P=A+(Y-C)cos J-(X-B)sin J←┚Lb1 1:Prog“JP”:J=H-180:I=(V-Y)cos J-(U-X)sin J:AbsI<1E-4=>Goto 2:≠=>P=P+I:Goto 1⊿←┚Lb1 2: D=(V-Y)÷s inHJC（超高加宽计算子程序）Fixm:Z[8]<1=>Z[8]=1E-8⊿I=E:M=A-Z[3]:E<1=>I=Z[8]:M=M-I⊿G=K:N=A-Z[3]+L:K<1=>G=Z[8]:N= N+G⊿←┚P≤0.5(N+M=>U=(P-M)÷I:≠=>U=(N-P)÷G⊿U<0=>U=0⊿U>1=>U=1⊿J=UT 高次抛物线加宽J=T(1-3UU+2UUU ) : D O<0=>J=0⊿I=1:Prog“JD”:I=US:G=I+UN-N:U=AbsD-V:U<0=>U=O⊿U>M+J+Z[5]=>U=M+J+Z[5]⊿Z[7]=1=>Goto 1:≠=>Goto 2←┚Lb1 1:I=G:D O>0=>I=-I:Abs I<N=>I=-N⊿⊿G=UI:N=I:Goto 3←┚Lb1 2:N<I=>N=I⊿G=N(M-U)+Z[5]Z[6]:N=-N:D O<0=>G=G+2IU:N=I⊿Goto 3←┚Lb1 3: U= M+J+V+Z[5]:Abs D≥U-0.01=>G=G-Z[5](N+Z[6])⊿JB（路基边桩放样子程序）H=H+G-Q:H<0=>I=2:H“W=”=Abs H◢≠=>I=3 :H“T=”◢(4850输入“T=”：H◢)⊿Prog“JD”:G=H÷V:G<1=>I=MG:≠=>G<2=>I=M+N(G-1):≠=>I=M+N+Z[6](G-2) ⊿⊿I“BP”=Abs D-VI-Z[7]Int G-U-Z[5] ◢G=V Frac G:G<1.5=>G“PT”◢⊿JD（断面数据储存子程序）I=1=>M=27:N=0.002:Z[5]=0.5:Z[6]=0.002:V=1.5:Z[7]=1:≠=> I=2=> Z[5]=0.8: M=0.75: N= 0.75: Z[6]=1.00: V=8: Z[7]=2.0 :≠=> Z[5]=0.3:M=1.50: N=1.75: Z[6]=2.00 : V=8: Z[7]=2.0⊿⊿ 符号说明：←┚◢⊿≠=> => ≤≥一、程序简介本套程序是共有2个主程序，7个子程序。

高速MC公路测量CASIO4800&4850万能坐标计算程序(完整版)

高速公路测量CASIO4800&4850万能坐标计算程序(完整版）程序特点：真正的全线贯通坐标正反计算、任意斜角计算！！！程序中加入测站点，真正的实现了“坐标法”与“极坐标法”两种放样方法的同时显示的功能，使得放样操作方法选择时更加灵活！！！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，避免频繁调用子程序，提高运算速度。

一、程序：ZBJSW“1.ZS 2.FS 3.SZ”:W=1=>Z[2]=0:V=0：Goto 1 ΔW=2=> Goto 4ΔW=3=> O “KOU LING”:O≠123456=>O=0: “OUT”◢Goto CΔO=0: V=0:Z[1]=0:Goto 0←┘Lbi 0←┘”N0.”：Z[1]+1 ◢Z[1]=0=>{ABCREFGUKO}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G “LS”:U“G”:K“X（00）”: O“Y（00）”: Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]= E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘Lbi A←┘Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“R0”：E “RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1：Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0←┘Lbi 1←┘{DZT }:D：Z：T“RJ”：Z[2]=0:Goto 2←┘Lbi 2←┘V≠1=>Z[2]>Z[1] =>GoToCΔΔD≤Z[Z[2]×8+9]=> A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]: U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2←┘Lbi 3←┘W=3 =>N=8：≠P＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝C+(NPQ+2UR)NS：L=1←┘X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR)LS,L,1,(N-1)))+ZCos（J+ T）←┘Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR)LS,L,1,(N-1)))+Z Sin（J+T）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO CΔW=3=>GOTO AΔZ=0=>“X（Z）=”:X:Pause 0: “Y（Z）=”:Y◢Pol（（X－K），（Y－O））←┘“S（Z）=”:I ◢J<0=> J=J+360Δ“F（Z）=”: J→DMS◢Goto 1ΔZ<0=>“X（L）=”:X:Pause 0: “Y（L）=”:Y◢Pol（（X－K），（Y－O））←┘fx4850①“S（L）=”:I ◢J<0=> J=J+360Δ“F（L）=”: J→DMS◢Goto 1ΔZ>0=>“X（R）=”:X:Pause 0: “Y（R）=”:Y ◢Pol（（X－K），（Y－O））←┘“S（R）=”:J ◢J<0=> J=J+360Δ“F（R）=”: J→DMS◢Goto 1 ←┘Z=0=> X “X（Z）=”◢Y “Y（Z）=”◢Pol（（X－K），（Y－O））←┘I“S（Z）=”◢J<0=> J=J+360ΔJ“F（Z）=”◢Goto 1ΔZ<0=> X “X（L）=”◢Y “Y（L）”◢Pol（（X－K），（Y－O））←┘fx4800②I“S（L）=”◢J<0=> J=J+360ΔJ“F（L）=”◢Goto 1ΔZ>0=> X “X（R）=”◢Y “Y（R）=”◢Pol（（X－K），（Y－O））←┘I“S（R）=”◢J<0=> J=J+360ΔJ“F（R）=”◢Goto 1 ←┘Lbi 4←┘{MH} ：M“X”：H“Y”：Z[2]=0：GOTO 5←┘Lbi 5←┘V=1：D= Z[Z[2]×8+9]：Z=0：T=90：GOTO 2←┘Lbi 6←┘K=（（H -B）Cos（C-90）-（M-A）Sin（C-90））×（（H -Y）Cos（J-90）-（M-X）Sin（J-90））：K≤0=> Goto 7ΔIsz Z[2]：Goto5←┘Lbi 7←┘D=F+Abs((H -B)Cos(C-90)-(M-A)Sin（C-90）)：D>G=> Isz Z[2]: Goto5ΔGoto 8←┘Lbi 8←┘V=2 ：GOTO 3←┘Lbi 9 ←┘K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：Abs K<（1÷E）^3=>Goto BΔD=D+K ：GOTO 8←┘Lbi B←┘V=3 ：Z=0：Goto 3←┘Lbi C←┘Z=(H-Y) ÷Sin(J+90)：“D”:D:Pause 0: “Z”: Z◢4850输出(Z=(H-Y) ÷Sin(J+90)：D“D”◢Z “Z”◢4800输出)GOTO 4←┘Lbi C←┘二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。

线路通用计算坐标程序(5800)

线路通用计算坐标程序(5800)程序符号说明：1、DK:放样点里程；2、B ：（B1）放样半宽（平行于S 角非切线边在里程点左右侧的宽度）；3、S （α1）：放样点与其对应的中线点连线和该中线点切线的夹角（左负右正，即放样点在线路左侧S 为负，在右侧为正）4、C ：(B2)垂直于S 角非切线边前后宽度（前正后负）；5、D （α2）切线（或平行于切线的直线）与S 角非切线边的垂直线的夹角（相对于切线逆时针偏转为负，顺时针偏转为正）；显示数据及其意义：输入J ，M ，0，R ，T ，F ，L ，A,I,E,K,B,S,C,D 值后显示第1个显示 X= 放样点X 坐标第2个显示 Y= 放样点Y 坐标第3个显示 Q= 如果Q=1,则只输入放样点里程、偏距、偏角；否则显示输入曲线要素DK-C 表示放样里程K 值超出本程序的计算范围值值XLTY20→DimZ ↵(扩展变量)Lbl 0:“JDLC=”?J (交点里程)：“XJ=”?M (交点X 坐标)：“YJ=”?O (交点Y 坐标)：“R=”?R (半径)：“T=”?T (切线长)：“L=”?F (曲线长)：“LS=”?L (缓和曲线长)：“A1=”?A (第一切线方位角)：“A2=”?I (第二切线方位角)：“E =”?E (转角符号)：Goto1↵Lbl1: “DK =”?K: “B1=”?B: “α1=”?S: “B2=”?C: “α2=”?D:J-T →Z (ZH 点里程)：Z+L →Z[1](HY 点里程)：Z+F →Z[3](HZ 点里程)：Z[3]-L →Z[2](YH 点里程)：M+TCos(A+180) →Z[4]( ZH 点X 坐标)：O+TSin(A+180) →Z[5]( ZH 点Y 坐标)：M+TCos(I) →Z[6]( HZ 点X 坐标)：O+TSin(I)→Z[7]( HZ 点Y 坐标)：32:40L L G R-→3030[4]()/()[8]LLZ C os A E G C os Z RRππ++→( HY点X 坐标): 3030[5]()/()[9]LLZ Sin A EG C os Z RRππ++→( HY点Y 坐标):Goto2↵Lbl 2:If K ≤Z:Then Goto A:Else If K ≤Z[1]:Then Goto B: If K ≤Z[2]:Then Goto C:If K ≤Z[3]:Then Goto D:Else Goto1:IfEnd: IfEnd: IfEnd: IfEnd ↵ Lbl A: M+(T+Z-K)Cos(A+180)+BCos(A+S)+CCos(A+D) →X ( 前直线X 坐标)：O+(T+Z-K)Sin(A+180)+B Sin (A+S)+C Sin (A+D) →Y ( 前直线Y 坐标)：GotoE ↵Lbl B:252()():40LK Z K Z P R ---→22230()30()90()[4]()/()()S K Z K Z K Z Z C os A E P C os BC os A E RL RL RL πππ+---++++290()()D K Z C C os A E XRLπ+-++→( 第一缓和曲线X 坐标):22230()30()90()[5]()/()()S K Z K Z K Z Z Sin A E P Cos BSin A E RL RL RL πππ+---++++290()()D K Z C Sin A E YRLπ+-++→( 第一缓和曲线Y 坐标):GotoE ↵LblC:90()90()90(22)[8]2()()()S K Z L K Z K Z L Z RSin C os A BC os A E R R Rπππ+-----++++90(22)()D K Z L CCos A EX Rπ+--++→( 圆曲线X 坐标):90()90()90(22)[9]2()()()S K Z L K Z K Z L Z RSin Sin A BSin A E R R Rπππ+-----++++90(22)()D K Z L CSin A EY Rπ+--++→( 圆曲线Y 坐标): GotoE ↵LblD:252([3])([3]):40LZ K Z K P R ---→22230([3])30([3])90([3])[6](180)/()()S Z K Z K Z K Z C os I EP C os EBC os I ERLRLRLπππ+---++-+-290([3])()D Z K C C os I E XRLπ+-+-→( 第二缓和曲线线X 坐标):22230([3])30([3])90([3])[7](180)/()()S Z K Z K Z K Z Sin I EP C os EBSin I ERLRLRLπππ+---++-+-290([3])()Sin D Z K C I E YRLπ+-+-→( 第二缓和曲线线Y 坐标): GotoE ↵Lbl E: “X=”：X ◢“Y=”：Y ◢ Goto Q ↵Lbl 9: “DK-C ” ◢Goto Q ↵Lbl Q: “Q=”?Q:If Q=1:Then Goto 1:Else Goto 0:IfEnd ↵。

新版万能坐标计算公式(0-18(带斜交通涵放样程序)

说明书
本公式采用积木法计算直、缓、圆曲线，使用者只要将曲线要素按提示录入到曲线要素表中，此路线的任意点坐标都能计算出来，对于工程测量的内业资料员整理资料是非常有用
的。

注意：曲线要素表编辑好后，请不要随便打开并修改，因此造成的损失本人将不付任何责
任。

作者QQ：254418166
本程序使用方法：请在曲线要素表中将每段线元的起点桩号，终点桩号，起点X坐标，起点Y坐标，起点桩号的路线方位角，起点桩号的半径，终点桩号的半径，线元判断值，一一输入到曲线要素表中，这样曲线要素表就开好了，然后想算坐标就在坐标计算工作簿中输入里程桩号及距中桩距离，就能算出相应的坐标。

有的设计图纸是给的交点坐标及桩号，使用者要将此交点及桩号转化成直直线元段，直缓线元段，缓圆线元段，圆圆线元段，然后按里程桩号输入进曲线要素表，就可以计算任意里程桩号的中桩和边桩任意坐标。

有什么不懂的请咨询设计本程的作者。

断值，圆向
左偏请输入
-1,向右偏
请输入1,直
本模板的里程桩号是从K0+000开始，到K18+000结束,半径如果是1E+45表示是直线段，表示半径无穷大。

坐标检查区域是对本段曲线元终点和下段曲线元的起点坐标复核，误差很大的时候请使用都仔细检查数据录入情况。

红色区域为公式区域，请不要去修改，绿色区域是使用者录入曲线要素区域。

坐标计算工作簿的第２已经隐藏，里面的作用是对你录入的里程桩号进行判断，若你要计算的里程桩号超出曲线要素表的范围，将会变成红色格子，
以示提醒。

密码：447290900。

坐标万能计算公式

万能坐标计算公式X=起点x+(待求点桩号-起点桩号)*【0.1184634425*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.046910077+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.046910077^2)+0.2393143352*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.2307653449+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.2307653449^2)+0.2844444444*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.5+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.5^2)+0.2393143352*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.7692346551+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.7692346551^2)+0.1184634425*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.953089923+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.953089923^2))+边距*COS(偏角+起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)】Y=起点Y+(待求点桩号-起点桩号)*【0.1184634425*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.046910077+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.046910077^2)+0.2393143352*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.2307653449+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.2307653449^2)+0.2844444444*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.5+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.5^2)+0.2393143352*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.7692346551+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.7692346551^2)+0.1184634425*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.953089923+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.953089923^2))+边距*SIN(偏角+起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)】切线方位角A=起始方位角+(终点曲率*(待求点桩号-起点桩号) +0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号) ^2)*180/3.14159265。

工程测量坐标正反算通用程序（终极篇）

工程测量坐标正反算通用程序（终极篇）第五篇坐标正反算通用程序(终极篇)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?DMS◢第6行：Goto 0K——计算点的里程BIAN——计算点到中桩的距离（左负右正）α——取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为ZBFS)第1行：”X1=”? C:”Y1=”?D:”K1=”?K第2行：Lbl 0:Prog “A”第3行：Pol(C-N,D-E):Icos(F-J)→S:K+S→K第4行：Abs(S)>0.0001=>Goto 0第5行：”K1=”:K◢第6行：”BIAN=”:Isin(J-F)→Z◢X1——取样点的X坐标Y1——取样点的Y坐标K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号Z——偏距（左负右正）注：在9860或9960中需将第3行替换为Pol(C-N,D-E): List Ans[1]→I ：List Ans[2]→J：Icos(J-F)→S:K+S →K，正反算主程序所有输入赋值多加一赋值符号（→），其他所有除数据库外的程序均保持不变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+S I)×90÷π→F第5行：F<0=>F+360→F: F>360=>F-360→F4. 数据库（命名为A）第1行：K≤175.191=>Stop（超出后显示Done）第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（第一缓和曲线）第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（圆曲线）第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（第二缓和曲线）第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（下一曲线的第一缓和曲线，示例为S型曲线，超出后显示Done）第6行：Lbl 1:Prog “XYF”A——曲线段起点的里程N——曲线段起点的x坐标E——曲线段起点的y坐标F——曲线段起点的坐标方位角P——曲线段起点的曲率(半径倒数，直线为0，左负右正)Q——曲线段终点的曲率(半径倒数，直线为0，左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加两个变量，具体方法参考本程序集中的第1篇辛普生公式的坐标计算通用程序（2）适用于任意线形：直线（0→P、0→Q）、圆曲线（圆半径倒数→P、圆半径倒数→Q）、缓和曲线（0或圆半径倒数→P、圆半径倒数或0→Q）、卵形曲线（接起点圆的半径倒数→P、接终点圆的半径倒数→Q），曲线左转多加一负号。

精典的全线三维坐标计算程序

精典的全线三维坐标计算程序，先传上坐标（CASIO4850）400字节主程序（ZB )01 Prog "DATE 1":L=Z-D:V=1:R=90/3.14BC02 L>L-S=>L=S-L:V=-1:C=U:G=G+HA:T=OΔ03 l<0=>X=L:Y=0:P=0:W=0:Goto1Δ04 LX=$(cos(RX2,0,L:Y=$(sin(RX2,0,L:P=RL2:W=L:Goto1Δ05 LP=RC(2L-C):X=BsinP+C/2-C^3/240B2:Y=B(1-cosP) + C2/24B:W=CΔ06 Lbl 1:X=VX:Y=HY:R=(X+Yi)(cosG+isinG)07 K=E-VT(cosG+isinG):O=R+K:P=G+HVP08 V=5+WQ(1+H)/2C:U=5+WQ(1-H)/2C09 "ZHONG X Y P":O◢ P->DMS◢"YOU X Y":M=O+V(cos(P+90)+isin(P+90◢"ZUO X Y":N=O+U(cos(P-90)+isin(P-90◢数据库（DATE 1)Z<缓直点=>A=xxx:B=xxx:C=xxx:D=xxx:E=xxx:G=xxx:K=xxx:Q=xxx:U=xxx:Goto 2...Lbl 1:U=CLbl 2:T=0.5C-C^3/240B2+(B+U2/24B-(B+C2/24BcosA)/sinAO=0.5U-U^3/240B2+(B+C2/24B-(B+U2/24BcosA)/sinAS=3.14AB/180+0.5(C+U)H=A/AbsA:A=AbsA说明：1.主程序中$为积分符号，注意i为复数的i,02计算第二缓和曲线，03计算直线，04计算缓和曲线，05计算圆曲线。

公路测量中的万能计算程序

公路测量中的万能计算程序使用卡西欧(CASIO—fx4850P)计算器进行编制：文件名：ZHY程序式：Lbl 0:A “X0”B “ Y0”:C “A0”:D “ 1- R0':E “1 + R1'：F “DK0':G “ DK1': Lbl 1:{P}:P “DKP -=(E^ Abs(GF):Q= Abs(P-F): “ J= ” :J= Abs(C+90Q(IQ+2D +n ) ▲ S“X" = A + / (COS ( C n+ 180+(2D+IX)X +2 r ,0,Q :C > 0= > Goto 2:工=> Goto 3A : Lbl2:T “ Y" =B+/ (Sin(C n+ 180+(2D+IX)X -2 r ,0,Q ) : Goto 4: Lbl 4:Lbl 3:T “Y" -=B(Sin(C n+ 180+(2D+IX)X -2 r ,0,Q : Goto 4: Lbl 4:{Z}:Z “BZM=J+R “Q': “ XP=' :V=S+Z COS M ▲“ YP=' :N=T+Z Sin M ▲ Goto 1说明：X0.Y0.A0---起点坐标和起点方位角R0 ------- 起点半径（直线为0）R1 ------- 终点半径（直线为0）1 --------- 线路转向:-1为线路左偏,1为线路右派偏DK0 ----- 起点桩号DK1 ----- 终点桩号DKP ----- 计算点桩号J --------- 计算点方位角Q -------- 所计算边距坐标与线路的任意夹角（右角）•法线角度为90度BZ ------- 中桩到边桩的水平距离，负为左边距，正为右边距•XP.YP------计算点桩号的坐标•方位角的计算公式1. 字母所代表的意义:X1：QD的X坐标屮:QD的Y坐标X2：ZD的X坐标y2：ZD的Y坐标S: QD- ZD的距离a：QD- ZD的方位角2. 计算公式：1）当y- y i>0, X2- x i>0 时："J丄d=3€0°+arctg —一^—2）当y2- y i<0, X2- x i>0 时：F —T3）当X2- x i<0 时：v '-二、平曲线转角点偏角计算公式1.字母所代表的意义：a i: QD- JD的方位角a 2: JD〜ZD的方位角3：JD处的偏角2.计算公式：3 = a 2- a i （负值为左偏、正值为右偏）三、平曲线直缓、缓直点的坐标计算公式1.字母所代表的意义：U: JD的X坐标V: JD的Y坐标A:方位角（ZH- JD）T:曲线的切线长,D: JD偏角，左偏为-、右偏为+ 2.计算公式:直缓(直圆)点的国家坐标： X' =U+Tcos(A+180 ) Y' =V+Tsi n(A+180°)缓直(圆直)点的国家坐标：X " =U+Tcos(A+D) Y' =V+Tsi n( A+D)四、平曲线上任意点的坐标计算公式 1. 字母所代表的意义： P:所求点的桩号B:所求边桩〜中桩距离，左-、右+ M 左偏-1，右偏+1 C: JD 桩号 D: JD 偏角 L s ：缓和曲线长 A:方位角(ZH- JD ) U: JD 的X 坐标 V: JD 的Y 坐标I=C-T :直缓桩号 J=I+L :缓圆桩号K=H+L 缓直桩号T :曲线的切线长,--- _厶 ISO圆缓桩号(f—色聲負 r 蛍审」■ _「 (f —氐+ £6乌+『I ((F—CUI KI6十汕06」滾十『|「汕O(Q6+/\AIAI+<U一S 8+E A H FO u一S O+(O8L巴—>|紀-汕(CXI(o 6+<u一S 8+E A H泸(06+<S 008-T X U丈->剖卷曰 (。

坐标计算程序及坐标计算公式

坐标计算程序及坐标计算公式一、坐标计算公式直线段：中桩公式：x=x1+(z-c)*cosay=y1+(z-c)*sina边桩公式：x=x1+(z-c)*cosa±d*cos(a-90)y=y1+(z-c)*sina±d*sin(a-90)说明： x1——起点x坐标，y1——起点y坐标，z——计算点桩号， c——起点桩号，a——方位角，d——距中桩距离。

“±”左边桩为“+”，右边桩为“-”。

二、方位角计算公式1、直线段方位角图纸提供。

2、若为单一圆曲线时，起点用直线段方位角图纸提供。

3、若为缓和曲线时：第一段缓和曲线方位角为直线段方位角图纸提供；第二段圆曲线起点方位角a1=a0±Ls*180/2∏r。

第三段缓和曲线方位角a2=a1±L*360/2∏r。

（a0为直线段方位角，Ls为缓和曲线长度，L为圆曲线长度，r为圆曲线半径，“±”右偏角为“+”、左偏角为“-”。

）三、5800计算器坐标程序坐标程序由1个主程序、5个子程序和1个数据库组成进入计算器编程模式（5：PROG—1：NEW新建程序名—3：EDIT），输写程序。

1、主程序adminFix 3（回车换行）Lb1 0:150→DimZ:“ZX→0,A→1,B1→2,B2→3,C→4,D→5,CR→6”?N:N→Z[149]:Prog “DAT2”：“ZS=>1，FS=>2”?N:If N=1:Then Goto 1:Else N=2=>Goto 2:IfEnd:Goto 0:Lb1 1:?S:“ANG=”?M:?Z:S=0=>Goto 0:0→N（回车换行）Lb1 5:Isz N:If S≤Z[8N+2]+Z[8N+4]:Then N→J:Prog“DAT1”:Else Goto5:IfEnd（回车换行）Abs(S-0)→W:Prog“SUB1”：“XS=”：X◢“YS=” ：Y◢“FWI=”：F-M→F：If F≤360:Then F→F:Else F-360→F:IfEnd:F►DMS◢Goto4（回车换行）Lb1 2:?X:?Y:“ANG=”？M：M→Z[148]：If M＜90:Then 180-M→M:IfEnd（回车换行）X-Z[4]:Y-Z[5]:X=0=>Goto 0:0→N（回车换行）Lb1 A:Isz N:N→Z[150]:Z[8N+3]-M→A:Z[8(N+1)+3]-M→B:Prog “ZX1”:If Z[6]×Z[7]≤0:Then N→J:Prog “DAT1”:Goto B:IfEnd（回车换行）Z[8N+3]+M→A:Z[8(N+1)+3]+M→B:Prog “ZX1”:If Z[6]×Z[7] ≤0:Then N→J:Prog “DAT1”:Goto B:Else Goto A: IfEnd（回车换行）Lb1 B:Prog “SUB2” （回车换行）Z[150]→N:0+W→S:If S>Z[8N+2]+Z[8N+4]+.001:Then Goto A: IfEnd（回车换行）If N>13:Then 0→N: Goto A: IfEnd（回车换行）If Z[148]>90: Then S+2Zsin(M-90) →S: IfEnd:“S=”:S◢“Z=”: Z◢“OK→2,NO→1”?N:If N=1: Then Z[150]→N: Goto A:Else Goto 2: IfEnd（回车换行）Lb1 4 :0→J:“DIST=”：Pol （X－Z[1]，Y－Z[2]）→I◢ J-F：If F<0：Then F+360→F：IfEnd：“FW=”： F►DMS◢ Goto 12、子程序DAT1Z[8J]→U：Z[8J+1]→V：Z[8J+2]→O：Z[8J+3]→G：Z[8J+4]→G：Z[8J+4]→H：Z[8J+5]→P：Z[8J+6]→R：Z[8J+7]→Q：1÷P→C：（P－R）÷（2HPR）→D：（180÷∏）→E：Return3、子程序ss“A”？→A：If A<0:Then 10A◢Else 9A◢IfEnd:Ans×1.05（回车）4、子程序SUB1.1184634425→A：.2393143352→B：.28444444444→N：.046910077→K：.2307653449→L：.5→Z[3]:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C +Z[3]WD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD)))→X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Nsin(G+QEZ[3]W(C+Z[3 ]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD))) →Y：G+QEW(C+WD)+M→F：X+Zcos（F）→X：Y+Zsin（F）→Y: Return5、子程序SUB2G-M→T：Abs（（Y-V）cos（T）-（X-U）sin（T））→W：0→Z：Lb1 0：Pros “SUB1” （回车换行）T+QEW（C+WD）→L：（Z[5]-Y）cos（L）-（Z[4]-X）sin（L）→Z：If Abs（Z）<1×10∧(-4):Then Goto 1:Else W+Z→W: Goto 0: IfEnd（回车换行）Lb1 1:0→Z:Prog“SUB1”：（Z[5]-Y）÷sin（F-2M+180）→Z：Return6、子程序ZX1（Z[5]-Z[8N+1]）cos（A）-（Z[4]-z[8N]）sin（A）→Z[6]：（Z[5]-Z[8（N+1）+1]）C5]-Z[8（N+1）+1]）cos（B）-（Z[4]-Z[8（N+1）]sin（B）→Z[7]：Return7、数据库DAT2If N=0:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）If N=1:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）说明：1、所有路线都是从Z[8]开始。

卡西欧 5800程序

卡西欧 5800程序卡西欧FX5800全线贯通万能正、反算程序时间:2021-01-27 23:56:25 来源:本站作者:xdw_815101 我要投稿我要收藏投稿指南FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版 ZHUCHENGXU 主程序 \1正算，输入2反算“NEW=0,OLD ≠0”?Z If Z=0:Then“X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifen d:Q=2=>Goto 2Lbl 1 :“KM=,<0 Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?LLbl Z:Z=1=> Prog“01”：Z=2=> Prog“02” 选择数据库文件，可增加 H-F→X:0.5(E-D)÷(G-F)→NC+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→PA+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→UB+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→V Q=2=>Goto 4：Cls：Fix 3\，1，U：\，2，V：“FWJ=”：P?DMS�� Norm 2：Cls：Goto 1Lbl 2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S“KMDG=”？H ：90→O:0→L：Goto Z （H线路范围内的任意桩号） Lbl 4:Pol（R-U，S-V）：J<0 => J+360→J While abs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J P+360→P:IF P-J>0:then -I→L:else I→L ifendGoto 3: Whileend：H+Icos（J-P）→H:Goto Z Lbl 3:Cls：Fix 3“KM=”: Locate 4，1，H：“PY=”: Locate 4，2，L�� Norm 2：Cls：Goto 201（数据库子程序） 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、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

5800公路线路计算万能程序

5800公路线路计算万能程序本程序以高斯四节点法为计算内核,多条线路正反算,中桩,边桩高程计算,超高计算,超高缓和计算,加宽计算,加宽缓和计算,边坡开口线计算,挡土墙坡脚线计算，考虑了中间绿化带的影响，适用与国家高速公路至乡村四级公路计算和放样，路基路面工程可以直接得出中边桩的设计三维坐标，去掉加宽和超高影响的计算困难，边坡和坡脚线计算可直接在边坡上提取坐标带入程序，经过所有设计因素的综合，得出更改边距并显示修改偏差后的坐标．反算速度明显提高,只需要3秒.程序无错，可放心输入，另本程序可以增加隧道超欠挖计算子程序,非常方便,因为每个隧道的断面数据不一样,所以在此没有明确写出.计算器主程序：ZHU-CHENG-XULbi0:“1，ZS=FS，2ZS，3FS，4XY＝＞SG，5。

”？U:U=1＝＞Prog”ZS-FS”: 进入公路三维程序U=2＝＞Porg”ZS”: 进入坐标正算程序U=3＝＞Porg”FS”: 进入坐标反算程序U=4＝＞Porg”XY ＝＞SG”: 进入大地坐标转施工坐标程序…………Goto0:说明：计算器总的主程序，进入选择各种分支计算程序。

1为公路三维计算，2为普通正算，3为普通反算，4为大地坐标转施工坐标。

选择错误重新选择。

此程序可以不用输入，只为给大家一个思路,可以把计算器所有程序集中到一个主程序内管理.公路三维部分正反算选择程序：ZS-FSDeg: //设置角度模式20→DimZ: //扩展变量“1LZ＝＞XY，2XY＝＞LZ,3BIANPO-FY”？U： //正反算选择，正算选1，反算选2,坡口坡脚选3If U=1: ThenProg”ZS-XH”:IfEnd: //进入正算循环主程序IfU=2:Then Prog”F S-XH”:IfEnd: //进入反算循环主程序IfU=3:ThenProg”BIANPO-XH”:IfEnd: //进入边坡开挖主程序正算循环主体程序ZS-XH“1PT-2SJ”？W： //普通计算和设计边距计算选择“XL-XZ“？U： //选择线路 1～NIf W=1: Then //运行普通计算模式边踞自由输入Lbi0:”L=”?L:”Z=”?Z: //输入桩号和边距Porg”ZS-XH-1“：Goto0:IfEnd:运行正算循环子程序If W=2:Then //运行设计计算模式边距按设计输入Lbi1:”L=”? L： //输入桩号Porg”0。

统一曲线参数输入万能坐标计算程序

统一曲线参数输入万能坐标计算程序姬根杰李永红中铁十七局集团四公司重庆 404100E-mail ：jcglj.ji@【内容摘要】运用复化辛普生公式、结合Excel VBA （Visual Basic for Application ）的编程，对于各种曲线线形,进行统一曲线参数输入,中桩、边桩坐标计算【关键词】曲线编程坐标计算在线路测量的计算工作中，离不开直线、圆曲线、缓和曲线（包括完整型和非完整型或卵型）的线元，对于各种线形的中桩坐标、边桩坐标的计算，真是百家争鸣，方法各异，尤其立交桥，组成立交的基本单元是匝道，匝道的平面线形组合相对比较复杂，计算烦琐，特别是卵形曲线的计算更加抽象难懂;主要有曲直法、解析法、拟合法、积木法、综合法、弦切法等。

有没有一种统一的解决方法呢?答案是肯定的, 运用复化辛普生公式、结合Excel VBA （Visual Basic for Application ）的编程，在我们普通的Excel 办公软件中，对于各种线形的中桩坐标、边桩坐标的计算，实现了真正意义上的的全线贯通。

使用此程序具有以下优点：1)使重复的批量计算任务自动化，减轻测量内外业计算的工作强度。

2)方便快捷，创建满足工程要求的特定报表。

3)便于测量数据的管理和维护。

一、复化辛普生的基本原理及计算思路抛物形公式 (又称辛普生Simpson 公式)： ()[()4()()]62b a b a a bf x dx f a f f b -+≈++⎰ (1)抛物形公式的几何意义：抛物形公式是用抛物线围成的曲边梯形的面积近似代替()f x 围成的曲边梯形的面积。

为了抛物形公式的计算结果满足精度要求，如果我们将积分区间划分成若干个小区间，在各小区间上采用低次的抛物形公式，然后再利用积分的区间可加性，把各区间上的积分加起来，得到新的求积公式，就是复化辛普生Simpson 公式的思想。

图1 抛物形公式示意图在线路中，无论直线、圆曲线、缓和曲线、卵型曲线无论那一种线元，可以看作由若干个小区间组合成的曲线，每个小区间的X、Y坐标增量，利用辛普生公式求积，再利用积分的区间可加性，只需要知道起点A的桩号La、半径Ra、X、Y坐标、方位角Fa以及终点B的桩号Lb、半径Rb，就可以计算AB段线元的任意点坐标和方位角。

5800 万能坐标程序TYQX

1.主程序TYQXLbl 0： 1÷P→C: (P-R)÷(2HPR) →D: If N=1 Then Goto1：Else Goto2：IfEndLbl 1: "Z"?Z： Abs(S-O) →W：Prog "SUB1"： "XS="：X◢"YS="：Y◢360Frac（(F-90+360) ÷360→Z[1]：“QXFWJ=”：Z[1]④DMS◢Prog"FY"ReturnLbl 2："X"？X： "Y"？Y：X→I：Y→J：Prog“SUB2”：“S=” ：O+W→S◢“Z=”：Z→Z◢Return2. 正算子程序SUB10.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F： 1-K→M：U+W(Acos(G+Q（180÷∏） KW(C+KWD))+Bcos(G+Q（180÷∏）LW(C+LWD))+Bcos(G+Q （180÷∏）FW(C+FWD))+Acos(G+Q（180÷∏）MW(C+MWD))) →X： V+W(Asin(G+Q（180÷∏）KW(C+KWD))+Bsin(G+Q（180÷∏）LW(C+LWD))+Bsin(G+Q（180÷∏）FW(C+FWD))+Asin(G+Q（180÷∏）MW(C+MWD))) →Y：G+Q（180÷∏）W(C+WD)+E→F： X+Zcos（F）→X： Y+Zsin（F）→Y：Return2. 反算子程序SUB2G-E→T： Abs((Y-V)cosT-(X-U)sin(T) →W：0→Z：Lbl 0：Prog "SUB1"：T+Q（180÷∏）W(C+WD) →L： (J-Y)cos（L）-(I-X)sin（L）→Z：If Abs（Z）<1E-6：Then Goto1：Else W+Z→W： Goto0 ：IfEnd Lbl 1 0→Z：Prog "SUB1" ：(J-Y)÷sin（F）→Z：Return数据库子程序TYQX-EPLbl 0："1.SZ => XY"："2.XY => SZ"：？ N：“J 。

线元法万能坐标计算程序（适用于CASIOfx-9750GⅡ计算器）

线元法万能坐标计算程序（适用于CASIOfx-9750GⅡ计算器）线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器)摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。

利用CASIO fx-9750G Ⅱ编程函数计算器强大的内存（可诸存63000个字符）和编程功能，编写各种计算程序，能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料，同时也使我们有更多的时间去挑战更富有创造性的工作。

关键词：坐标放线线元测量程序1、前言本程序采用Gauss-Legendre（高斯-勒让德）五节点公式作内核,计算速度（太约2秒）适中，计算精度很高。

在此之前，本人曾用过以下公式作内核：①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点，⑦高斯-勒让德5节点，经过测试③计算最快，⑦代码稍长但计算速度只比③⑥稍慢，精度最高，可满足线元长小于1/2πD的所有线形的精度要求。

⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55mm,3/4πD时偏差为31.63mm，πD 时偏差为968mm，偏差按半径倍数增大，如线元长大于1/2πD（1/2圆周长）时，可将其拆分二个或多个线元单位，以确计算保精度。

2、程序特点事先将所有的平曲线交点的线元要素诸存到计算器内，测量时只输桩号、边距等程序会自动寻找各类要素，一气呵成地完成施工测量任务，中途不需人工转换各类要素数据，本程序可诸存几百条线路的要素数据，计算时可按需选择线路编号进行测量。

测量时不需查阅及携带图纸，仅一台CASIO fx-9750GⅡ编程函数计算器即可。

FX5800计算器公路全线坐标正、反算计算程序

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、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

5800计算器全线坐标计算放样程序

5800计算器全线坐标计算放样程序（修改版）“XLZBJSCX” ◢LB1 0 ↙CLS : FIX 4 : 30→DIM Z ↙“XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y):Pol(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":J→DMS◢(后视方位角)LB1 1 ↙“K=”?K ◢（计算里程）IF K<本段曲线终点里程AND K≥上段曲线终点里程：THEN本段终点里程→Z[1] :上段曲线终点里程→Z[2] ：1→O （注：左偏曲线输入-1→O,右偏曲线输入1→O）: 偏角→A：半径→R : 第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : GOTO 2 : IFEND↙…………（曲线段分段输入）补充直线段输入如下(止点独立直线段)IF K<本段直线终点里程 AND K≥本段直线起点里程:THEN 1→O:本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E:GOTO 4:IFENDLB1 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])→S ↙（曲线总长）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]-S→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] →X : 0→Y : E→T : PROG“TYZBCX” ：GOTO 1 ↙LB1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →X ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Y ↙90*H^2/( R*π* Z[6]) →T ↙IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙PROG“TYZBCX” ：GOTO 1 ↙LB1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*SIN T+ Z[8]→X ↙R*(1-COS T)+ Z[10]→Y ↙IF O>0：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙PROG“TYZBCX” ：GOTO 1 ↙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*SIN A→X ↙Z[13]*SIN A-U*SIN A+V*COS A→Y ↙IF O>0：THEN F-T→T : T<0=>360+T→T : ELSE F+T →T : IFEND ↙PROG“TYZBCX” ：GOTO 1 ↙子程序：“TYZBCX” ↙（统一坐标计算）IF O<0 : THEN -Y→Y : IFEND ↙“QXJ=” :T◢（计算里程点切线方位角，可以不显示）Z[16]+X*COS E-Y*SIN E→Z[18] ↙Z[17]+X*SIN E+Y*COS E→Z[19] ↙“XI=” : Z[18] ◢（中线X）“YI=” : Z[19] ◢（中线Y）Pol(Z[18]-M,Z[19]-N):"DI=":I（中桩放样距）◢J<0=>J+360→J:"FI=":J→DMS◢（中桩放样方位角）“PJ=”?P◢ (输入边桩与线路夹角，左-右+)“PD=”?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):"DP=":I◢（边桩放样距）J<0=>J+360→J:"FP=":J→DMS◢（边桩放样方位角）RETURN↙经演算可以放心使用（短链情况：IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN IF K>短链来向里程：THEN K-短链长度→K:IFEND: 本段终点里程-短链长度→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 K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN IF K>长链来向里程：THEN K+长链长度→K:IFEND: 本段终点里程+长链长度→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↙）2009年5月12日星期二修正使用卢海源。

万能全线坐标计算程序

单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段Mail:ceo@ Site: Time:2010-04-28 14:21坐标计算表单位名称：中铁八局集团有限公司遂渝铁路二线工程指挥部第一项目经理部工程名称：改建铁路遂宁至重庆增建二线Ⅰ标段计算：陈超复核：日期：。

5800计算器全线坐标计算放样程序(交点法)

王建清“ZBZFS” (主程序)LB1 A ↙CLS ：FIX 4：30→DIM Z ：“0=ZS，1=FS”？Z↙（0坐标正算，1坐标反算）IF Z=0：THEN Goto B：IFEND↙If Z=1:THEN Prog“FSLCBZ”:GOTO A: IFEND↙LB1 B↙“K=”?K ：“SJK？0=SJK0 1=SJK1…..”?N：N=0=> Prog“SJK0”：N=1=> Prog “SJK1”：....N=20=> Prog“SJK20”: Prog“ZBFY”: GOTO A↙（K是要求点里程桩号，N为选择线路数据库）“FSLCBZ”(反算里程边桩子程序)LB1 1↙“XK=”?X：“YK=”?Y↙（反算点X,Y坐标）“K=”?K↙（该点大致在线路中桩号）LB1 2“SJK？0=SJK0 1=SJK1…..”?N：N=0=> Prog“SJK0”：N=1=> Prog“SJK1”：.... N=20=> Prog“SJK20”:T-90→W:W<0=>360+W→W：Abs((Y-Z〔19〕)×COS(W) -(X-Z〔18〕)×SIN(W))→SIF S<0.0001：THEN GOTO 4：Else Goto 3：IFEND↙LB1 3↙K+S→K：N=0=> Prog“SJK0”：N=1=> Prog“SJK1”：....N=20=> Prog“SJK 20”: T-90→W:W<0=>360+W→W：Abs((Y-Z〔19〕)×COS(W)-(X-Z〔18〕)×S IN(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〕) ：J<0=>J+360→J：(J-T)<0=>-I→I↙CLS:Locate 1,1,“DP(-Z,+Y)=” :Locate 1,2,I: Locate 1,3,“K=”:K+S→K: L ocate 4,3,K◢（分行分列一次显示反算所得里程和边距数据）Return↙“ZBFY”(坐标显示子程序)LB1 0↙“QXFWJ=”:T▲DMS◢（显示该点切线方位角）CLS:Locate 1,1,“K=” :Locate 4,1,K: Locate 1,2,“X=”: Locate 4,2,Z〔1 8〕: Locate 1,3,“Y=”: Locate 4,3,Z〔19〕◢（分行分列一次显示计算所得里程和X,Y坐标数据）CLS:”B JIAO=”?P:”B D=”?D↙（B为以路线前进方向夹角（一般正交为90度），D为距中线距离左为“-”右为“+”）Z〔18〕+D×COS(T+P)→Z〔20〕↙Z〔19〕+D×SIN(T+P)→Z〔21〕↙CLS:Locate 1,1,“K=” :Locate 4,1,K: Locate 1,2,“D=”: Locate 4,2,D: Locate 1,3,“B X=”: Locate 5,3,Z〔20〕:Locate 1,4,“B Y=”: Locate 5, 4,Z〔21〕◢（分行分列一次显示桩号，距中距离，X，Y数据）“XLZBJSCX”(线路坐标计算程序)LB1 2↙Z〔6〕÷2－Z〔6〕^3÷（240R^2）+ Z〔6〕^5÷(34560×R^4) →Z〔8〕↙Z〔7〕÷2－Z〔7〕^3÷（240R^2）+ Z〔7〕^5÷(34560×R^4) →Z〔9〕↙Z〔6〕^2÷(24×R)－Z〔6〕^4÷（2688×R^3）→Z〔10〕↙Z〔7〕^2÷(24×R)－Z〔7〕^4÷（2688×R^3）→Z〔11〕↙π*A*R/180+0.5*( Z[6]+ Z[7])→S ↙（曲线总长）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] ↙（Z H点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] ↙(Y H点里程)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)-H^13/(599040*R^6* Z[6]^6) +H ^17/(178472640*R^8* Z[6]^8)→Z[23]↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3)+H^11/(42240*R^5* Z[6]^5)-H^15/(967680 0*R^7* Z[6]^7)+H^19/(3530096640*R^9* Z[6]^9)→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)-H^13/(599040*R^6* Z[7]^6) +H ^17/(178472640*R^8* Z[7]^8)→U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3)+H^11/(42240*R^5* Z[7]^5)-H^15/(967680 0*R^7* Z[7]^7)+H^19/(3530096640*R^9* Z[7]^9)→V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]COS A+ Z[12]-U*COS A-V*SIN A→Z[23] ↙（坐标转换）Z[13]*SIN A-U*SIN A+V*COS A→Z[24] ↙（坐标转换）IF O>0 ：THEN T +E→T ：ELSE E-T →T ：T<0=>360+T→T ：IFEND ↙GOTO 8 ↙LB1 8↙（施工坐标转换大地坐标）IF 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↙“SJK.. ...”(数据库子程序)IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN 本段终点里程→Z [1] : 上段曲线终点里程→Z[2] ：1→O（注：左偏曲线输入-1→O,右偏曲线输入1→O）: 偏角→A：半径→R ：: 第一缓和曲线→Z[6] ：第二缓和曲线→Z[7] ：交点X→B ：交点Y→C ：小里程向交点方位角→E ：交点向大里程方位角→F ：Prog”XLZBJSCX”： Return: IFEND↙…………（曲线段分段输入）补充单独直线段输入如下：IF K<本段直线终点里程AND K≥本段直线起点里程:THEN本段直线终点里程→Z[3]:终点X→Z[16]：终点Y→Z[17]：方位角→E：E→T：Z[16]+(K- Z[3]) COSE →Z[18]：Z[17]+(K- Z[3]) SIN E →Z[19]：Return：IFEND。