4800放样程序

1 XY边长及方位角反算Prog"A":L=J:Prog"B"◢Fix 3:I"D="◢Lbl 1:{S}:V"X"=X+Rec(S, L ◢U"Y"=Y+J ◢Goto 1输入：X1,Y1,X2,Y2.输出: L°=,D=.2 XY-1输入里程桩号计算中间点坐标Prog"A":L=J:Prog"B"◢Fix 3:I"D="◢{Z}: Z:Lbl 1:{G}:V"X"=X+Rec(Abs(G-Z, L ◢U"Y"=Y+J ◢Z-起始点里程桩号Goto 1Ｇ－任一点里程桩号3 JD两直线交点坐标计算Prog"A":Prog"C":L"L1":Prog"F":E=W+K:{L}:L"L2":Prog"F":C=V+K:S=tan E:T=tan C:G=S-T:A=A+0.0001U"X"=(SA-MT-B+N)/G◢Z"Y"=S(U-A)+B◢Pol(U-M,Z-N:I"D="输入X1,Y1,X2,Y2,X3,Y3.X4,Y4,L1,L2,输出X=,Y=,D=,4 ZBFD坐标放点Prog"A":Lbl 3:Fix 3:{MN}:Pol(M"X3"-A,N"Y3"-B:I"D="◢J=J-W+180:J<0=>J=J+360L=J:Prog"B"◢Goto 3输入:X1,Y1,X2,Y2,X3,Y3.输出:D,L°5 Z自由导线推算E=0:1F=0:Prog"A":Fixm:C=A:D=B:J[F]=J:Lbl 4:{LH}:L"L0":Prog"F":I[F]=J[F]+K-180:H[E+F]=H"D":Prog"R":C"X="◢D"Y="◢J[F]=I[F]:Goto 46 ZDX支导线计算X1Y1,X2Y2,为已知点坐标Prog"A":L°-角度用小数输入Lbl 4:输一组L°、D.得一个点的坐标{LS}:L:Prog"F":M"X3"=A+Rec(S,W+K+180◢N"Y3"=B+J◢Goto 4输入:X1,Y1,X2,Y2,L,S.输出:X3,Y3.L°—观测角度以下按小数点输入7 DX单一导线简易平差计算Prog"A":Prog"C":J=W:Ｉ=V:C=A-M:D=B-N:F=E"N":I[F]=J:Lbl 1:F=F-1:{GH}:I[F]=J[F]+G"L°"-180:I[F]<0=>I[F]=I[F]+360:≠>I[F]>360=>I[F]=I[F]-360H[E+F]=H"D":F>1=>Goto 1{G}:G=I[F]+G"L°"-180:G<0=>G=G+360:≠>G>360=>G=G-360I=G-I:W=12″[S]=620.790G"W"=602FracI ◢相对中误差:fS/[S]=1/88119fX=-0.007fY=0.000G=0:Lbl 2:I[F]=I[F]-I(E-F)/E:Prog"R":说明:F=F+1:N-观测角数,L°-左角按度分秒输入G=H+G:W-角度闭合差,fX,fY 坐标闭合差F<E=>Goto 2若用外(右)角计算闭合差程序中之I=G-I-360G"[S]="◢若用内(左)角计算闭合差程序中之I=G-II"fX"=C ◢本程序只能计算八个未知点,多于八个时每增加一个C=A:未知点,必需扩展两个储存器.A"fY"=D ◢D=B:B"fS/[S]"=Int(G/ √ (I I+AA ◢Lbl 3: F=F －1:Prog"R":C"X"=C-I H/G ◢D"Y"=D-AH/G ◢F>1=>Goto 38 WDX 无定向导线Prog"A":D=I:C=0:M=0:N=0:P:Lbl 0:P=P-1:{LS}:L"B°":Prog"F":F=K:C=C+F:C=0=>Goto 1C=C+180:Goto 1X1Y1,X2Y2,-已知点坐标Lbl 1:P-边数,B°-观测角,第一个角输入0,Rec(S,C:S-边长,Xf,Yf-增量累加数,记录侍用Fix 4:输入一组Xf,Yf,求得一组X Y,M"Xf"=M+I ◢用左角,度以下按小数输入.N"Yf"=N+J ◢P>0=>Goto 0Pol(-M,-N:R=W+180-J:K=D/I:V=Cos R:点号观测角边长Xf Yf X Y E=Sin R:Lbl 2:{HO}:H"Xf":O"Yf":Fix 3:T"X"=X+K(HV-OE ◢Z"Y"=Y+K(HE+OV ◢Goto 2B 为观测角,第一个角输入04119.5814110.7005102.2204969.5125029.9914887.630142.64164.8800-66.203843.11280.0000-211.5145-303.145023 0°00′00″ 42 28 06 277 33 42 A 1164.88313.274329.1104426.9529 QJD前交点(正弦公式)Prog"A":Prog"J":S=I Sin E/Sin(E+K:Prog"O"A°、B°-观测角,X1、Y1、X2、Y2-正切公式Prog"A":Prog"J":S=tan E:T=tan K:G=S －1＋T -1:M"X3"=(X/T+A/S-Y+B)/G ◢N"Y3"=(Y/T+B/S+X-A)/G 输入:X1,Y1.X2,Y2,A°,B°,输出:X3Y310 FX 方向交会Prog"A":Prog"J":S=tan E:T=tan K:G=S-T:M"X3"=(SX-AT-Y+B)/G ◢N"Y3"=S(M-X)+YX1Y1,X2Y2,为起算点坐标.A°,B°-方位角.X,Y--交点坐标11 CJ 侧方交会点Prog"A":Prog"J":S=I sin(E+K)/Sin K:{K}:K=0=>M"X"=A+Rec(S,W+E-180◢≠N"Y"=B+JX1Y1,X1'Y1',X2Y2--已知点坐标,A°B°-观测角,左右分开.左K=0,右K=1A=63° 33′06″B=92° 20′51″从左右两组算出坐标取中数为XYA ′=109° 59′39″B ′=50° 57′06″12 HJ后交点计算Prog"A":M"X3":X1=992.692X2=1012.499N"Y3":Prog"J":S=tan E:T=tan K:H=(Y-B)/S+X-A:F=(X-A)/S-Y+B:Q=(N-B)/T-M+A:V=(M-A)/T+N-B:K=(H+Q)/(F+V:X3=1063.248D=(H-FK)/(1+K2:Y4=999.997Y3=1000.000O"X4"=A+D◢X1Y1,X2Y2,X3Y3---已知点坐标P"Y4"=B+KD A°、B°----观测角X4Y4--测点坐标12HJ—1Prog〝A〞︰D=I︰（另一种后Pol（A－M〝X3〞，B－N〝Y3〞︰交公式）S=I︰G=J︰Prog〝J〞︰Z=180－（W－G＋E＋K）÷2︰Q=tan－1（D Sin C÷S Sin E︰H=tan－1（tan Z÷tan（45＋Q︰R=Z＋H︰T=D Sin（E＋R）÷Sin E︰O〝X〞=X＋T Cos（W＋R◢P〝Y〞=Y＋T Sin（J＋R13 LDHJ两点后交Prog"A": L"A°":Prog"F":C=K:{L}:X2=857.025X1=944.202L"B°":Prog"F":O=K:{L}:L"E°":Prog"F":E=K:{L}:L"F°":Prog"F":F=K:Y3=139.053N=Sin O/Sin(O+F:P=√(S2+M2-2SM Cos(C-O:观测角A°B°E°F°度以下按小数点输入.T=√(D2+N2-2DN Cos(F-E:K"D1"=2I"D0"/(P+T◢S"S1"=KS◢D"S2"=KD◢M"S3"=KM◢N"S4"=KN◢H=Cos-1((I2+S2-M2)/2IS:L=H:Prog"B":Z"H°="◢G=Cos-1((I2+N2-D2)/2IN:L=G:Prog"B":Z"G°="◢Fix 3:Q"X3"=A+Rec(S,W+H+180◢R"Y3"=B+J◢Q"X4"=A+Rec(D,W+E-F-L◢R"Y4"=B+JProg"A":Lbl 0:Prog"W":Fix 2:S=D(Sin K)2◢H=S/tan K+G◢{L}:Prog"F":输入:X1,Y1,X2,Y2,H0-起点高,T仪高,Prog"O"◢D-视距,S-平距,Z°—天顶距,P标高Goto 0输出:X3,Y3,S平距15 JZB极坐标法测碎部点坐标、高程Prog"A":G=E"H0"+T:Lbl 8:{LSV}:Prog"F":Prog"O"◢H=G+V"h"-P◢Goto 8输入:X1,Y1,X2,Y2, H0,T,L°，S，h求:X,Y,H.X3=273.29916 SD 垂距,垂点坐标计算Y3=208.502Prog"A":O=Sin J:P=CosJ:E=M"X3"-X:F=N"Y3"-Y:S=EP+FO ◢D=EO-FP ◢Y=252.157V"X"=X+Rec(S,W ◢U"Y"=Y+J注:D 为正时,3点在1-2之左侧;S 为正,3在1-2之线上;为负则在2--117 MM 坐标法计算多边形面积N: P=0: V=0:A"X":B"Y":E=A: F=B:Fix 3:Lbl 0:{XY}:X: Y:N=4(边数)Lbl 1:Pol(X-A,Y-B: I"S"◢P=P+(Y-B)(X+A)/2:A=X:B=Y:V=V+1:V ≠N=>Goto 0P"m 2"=Abs P ◢Q"mu"=1.5 m P18断面图封闭图形面积计算(高程相当于X,距离相当于Y,用相对坐标计算,)N:P=0:V=0:Y=0:A"H":B"D":E=A:69.8F=B: Fix 3:Lbl 0:N=V+1=>H=E: Y=F: Goto 1｛HD ｝:H:Y=Y+D:Lbl 1:m =32.795A=H:B=Y:N=多边形边数,V=V+1:H=转折点高程V ≠N=> Goto 0D=相邻两点距离P"m 2"=Abs P第一点D 输入0。

1.主程序(ZCX1)"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=>Goto 1：≠>Goto 2Δ换行Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FS"=F-90◢Prog“ZDS”2.Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢Prog“ZDS”2. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+QEKW (C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD)) +Acos(G+QEMW(C+MWD)))：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+ Asin(G+QEMW(C+MWD)))：F=G+QEW(C+WD)+90：X=X+ZcosF：Y=Y+ZsinF换行3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ换行Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF换行Lbl 0：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙Lbl 1：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙。

1.)TYQXJS(程序名)Defm 2:“1.SZ=>XY”: “1.XY=>SZ”:N:N=1=>Goto 1≠>Goto 2⊿←Lb1 1:｛SZ｝:SZ←S≤每段线终点桩号=>J=1：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=2：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=3：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=4：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=5：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=6：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=7：Prog“DAT1”: Goto 3⊿←Lb1 3:W=Abs(S-O): Prog“SUB1”:X“XS”=X▲Y“YS”=Y▲Z[5]=X：Z[6]=Y：F“FS”=F-90▲Lb12:Prog“DAT1”: ｛XY｝:XY:I=X:J=Y:Prog“SUB2”:S“S”=O+W▲Z“Z”=Z▲Goto 22.)SUB1(程序名)A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1-L:Z[1]=1-K:X=U+W(A COS(G+QEKW(C+KWD))+B COS(G+QELW(C+LWD))+B COS(G+QEFW(C+FWD))+A COS(G+QE Z[1]W(C+Z[1]WD))):Y=V+W(A SIN(G+QEKW(C+KWD))+B SIN(G+QELW(C+L WD))+B SIN(G+QEFW(C+FWD))+A SIN(G+QEZ[1]W(C+Z[1]WD))):M=90:F=G+QEW(C+WD)+M:X=X+Z COS F:Y=Y+Z SIN F←3.)SUB2(程序名)M=90:T=G-M:W=Abs((Y-V)COS T-(X-U)SIN T):Z=0:Lb10:Prog“SUB1”:L=T+QEW(C+WD):Z= (J-Y)COS L-(I-X) SIN L: Abs Z＜1E-6=> Goto 1:≠>W=W+Z: Goto 0⊿←Lb1 1:Z=0: Prog“SUB1”:Z=(J-Y)÷SIN F4.)DAT1(程序名)J=1=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=2=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=3=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=4=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=5=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=6=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=7=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿……………………………………………………U“X0”：V“Y0”：O“S0”：G“F0”：H“LS”：P“R0”：R“RN”：Q：C=1÷P：D=（P-R）÷（2HPR）：E=180÷∏←说明：1．Q-------以道路中线前进方向，直线输0，左拐曲线输-1，右拐曲线输12．P------ 起点曲率R----- 终点曲S-----桩号Z----边桩距离左为负，右为正，0为中线坐标，N输入1时为正算，输入2时为反算；M为正交或斜交，正交输90，斜交输设计数M在程序内容里改3．当线段为直线时，起点、止点的曲率半径输入1E45；当线段为圆曲线时，起点、止点的曲率半径均输入圆的半径；4．当线段为完整缓和曲线时，起点与直线相接时，起点曲率半径输1E45，与圆曲线相接时，起点曲线半径等于圆曲线半径；止点与直线相接时，止点曲率半径输1E45，与圆曲线相接时，止点曲率半径等于圆曲线半径5．当线段为非完整缓和曲线时，起点与直线相接时，起点曲率半径等于设计规定值，与圆曲线相接时，曲线半径等于圆曲线半径；止点与直线相接时，终点曲率半径等于设计规定值，与圆曲线相接时，曲线半径等于圆曲线半径6．本程序可以连续计算整条路的坐标，但要把整条路的要素输入计算机，也可以计算一部分。

4800竖曲线计算放样程序SJBGJS(主程序)Lb1 1:｛KB X｝:Prog“YSJG”:N=-V－W:N>0＝＞N=-1:≠＞S=0.14L÷（0.07+0.035I）：Y=（0.07+0.035I）÷LC=Abs（K-Q）:K>Q＝＞C<T＝＞D=（K-Q+T）2÷（2RLb1 2:C<T＝＞D=（K-Q-T）2÷（2RLb1 3:K<Q＝＞P=V÷100: ≠＞P=W÷＝＞Prog“123”NormProg“BBBB”:Goto 1123K<J-L＝＞Prog“11”: ≠＞K<J＝＞Prog“22”: ≠＞K≤○＝＞Prog “33”: ≠＞K≤○+L＝＞Prog“22”: ≠＞Prog“11”33G=G+0.01I（MB）11G=G-Abs B×0.0222M=1＝＞Prog“44”: ≠＞Prog“55”44K≤J-L+S＝＞Prog“E”: ≠＞K<J＝＞Prog“C”: ≠＞K<○+L-S＝＞Prog“C”: ≠＞Prog“F”55K≤J-L+S＝＞Prog“A”: ≠＞K<J＝＞Prog“C”: ≠＞K<J＝＞Prog “C”: ≠＞K<○+L-S＝＞Prog“C”: ≠＞Prog“B”EB≤0＝＞G=G+0.02B: ≠＞G=G+B((K-J+L)×Y-0.07)÷FB≤0＝＞G=G+0.02B: ≠＞G=G+B((○+L-K)×Y-0.07)÷CK<J＝＞G=G((K-J+L)×Y-0.07)÷3.5×MB: ≠＞G=G+((○+L-K)×Y-0.07)÷3.5×AB>0＝＞G=G-0.02B: ≠＞G=G-B((K-J+L)×Y-0.07)÷BB>0＝＞G=G-0.02B: ≠＞G=G-B((○+L-K)×Y-0.07)÷BBBBLb1 0:｛Z｝H=G-Z:H>0＝＞Prog“BBB”: ≠＞Prog“SBB”Fix 3:C=C+0.0001:“LL=”:Pause 0BBF“P11”U“P22”H≥-10＝＞C=Abs B+1.9-（H+0.03）F: ≠＞C=Abs B+1.9+9.97F+1-（H+10.03）H<8＝＞C=Abs B+1.5H+0.3: ≠＞H≤16＝＞C=Abs B+13.433+（H-8）×1.5+0.3: ≠＞Abs B+26.865+（H-16）×ACB>0＝＞G=G-○×B: ≠＞G=G-B((K-J+L)×Y-0)÷BCB>0＝＞G=G-○×B: ≠＞G=G-B((○+L-K)×Y-0)÷SBBF“P11”U“P22”H≥-10＝＞C=Abs B+1.9-（H+0.03）F: ≠＞H≥-20＝＞C=Abs B+1.9+9.97F+1-（H+10.03）U：≠＞H≥-30＝＞C=Abs B+1.9+9.97F+9.97U+2-（H+20.03）U: ≠＞H≥-40＝＞………SJGK≤30231.33＝＞Q=30350:A=311.580:R=15000:T=105.061:V=-0.967: W=+2.368:J=30049.62:○=30141.33:L=90:M=+1:I=7:Goto 1: ≠＞K ≤30594.22＝＞Q=30350:A=311.580:R=15000:T=105.061:V=-0.967:W =+2.368:J=30311.33:○=30391.27:L=80:M=+1:I=7:Goto 1: ≠＞K≤31123.63＝＞Q=30350:A=311.580:R=15000:T=105.061:V=-0.967:W=+ 2.368:J=30674.22:○=30768.06:L=80:M=+1:I=7:Goto 1: ≠＞K≤31 123.63＝＞………Lb1 1Q----为该段竖曲线的变坡点桩号A----变坡点高程R----竖曲线半径R----竖曲线切线长E----竖曲线外距P1---变坡点前一个纵坡坡率P2---变坡点后一个纵坡坡率N----竖曲线纵要素因子，凸为-1，凹为+1 K0---直线直线段边桩放样辅助程序BBBBBB（主程序）Defm 3Lb1 0:｛G ,Z｝H=G-Z:H>0＝＞Prog“BBBBB”: ≠＞Prog“BBBB”Goto 0BBBBBH<12＝＞L=9+1.5H: ≠＞H≤20＝＞L=9+8×1.5+1.5+（H-8）×1.75: H ≤20＝＞L=9+8×1.5+1.5+8×1.75+3+（H-16）×2:BBBBC“P1”:D“P2”:S“P3”H>-12＝＞L=11.25-（H+0.045）C: ≠＞H>-20＝＞L=11.25+9.955 C+1.5-（H+10+0.045）D: ≠＞L=11.25+9.955 C+9.955 D+3-（H+20+0.045）说明：G---路肩设计高Z---实测高L---理论宽度H---填挖值。

4800公路测量放样计算程序CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F －Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。

卡西欧Fx-4800P计算器测量程序一、坐标放样源程序ZBFYU“Xc”:M“Yc”:O“XH”:Q“YH”:I=O-U:J=Q-M: Pol(I，J):J＜0＝＞J=J+360:≠＞J=J:⊿T“HF” =J◢I“HL”◢X=0=＞W“XQ” =U+Lcos(T+P“a”)◢S“YQ” =M+Lsin(T+P)◢≠＞I=X-U:J=Y-M: Pol(I，J):J＜0＝＞J=J+360:≠＞J=J⊿I“QL=”◢H“QF” =J◢J=H-T:J＜0=＞J=J+360◢≠＞J=J使用说明:XH—后视点X坐标 YH—后视点Y坐标XC—置镜点X坐标 YC—置镜点Y坐标HF—后视方位角 QF—前视方位角HL—后视距离 QL—前视距离X、Y—前视点坐标 J—仪器后视归零后拨角当X=0时运行坐标反算QL—测设点到置镜点的距离а—测设点与后视点夹角XQ—测设点的X坐标 YQ—测设点的Y坐标二、复合辛卜生公式计算曲线元坐标源程序XBS—XY(已知里程求坐标)A“QD”:Z“ZD”:R“R0”:V“RN”:D“X0”:E“Y0”:F“F0”:K:L“B”:P“а”:J:Prog “SIM1”:W＜0＝＞W=W+360:≠＞W=W⊿W“FWI=”◢X=X◢Y=Y◢XBS—L、B(已知坐标求里程和边距)A“QD”:Z“ZD”:R“R0”:V“RN”:D“X0”:E“Y0”:F“F0”:K”K”(输入大致里程)I=X:B=Y:K=A+Abs((Y-E)cos(F+P“a”-180)+(X-D)sin(F+P-180):Lbl 1:L=0:Prog“SIM1”:C=(B-Y)cos(W+P-180)- (I-X)sin(W+P-180): AbsC＜10∧-4＝＞Goto2: ≠＞K=K+C∶Goto1:⊿Lbl2:L=0:Prog“SIM1”:C=(B-Y)÷sin(W+P):K=K◢C“B=”◢子程序SIM1T=(J÷V-J÷R)÷(Z-A)H=(K-A)÷NS=90H÷πW=F+((TNH+J÷R)+J÷R)NS: G=1:X=D+H÷6×(cosF+cosW+4×Σ(cos(F+((G+0.5)HT+2J÷R)(G+0.5)S),G,0,(N-1))+2Σ(cos(F+(THG+2J÷R)GS,G,1, (N-1)))+Lcos(W+P):Y=E+H÷6×(sinF+sinW+4×Σ(sin(F+((G+0.5)HT+2J÷R)(G+0.5)S,G,0,(N-1))+2×Σ(sin(F+(THG+2J÷R)GS,G,1,(N-1)))+Lsin(W+P)使用说明:QD—起点里程 ZD—终点里程R0—起点半径 RN—终点半径J—曲线转向 (左转为-1 直线为0 右转为+1)X0、Y0—线元起点坐标 F0—起点方位角K—待求点里程 N—计算精度(可取4～12的整数)а—边桩与中线右交角 B—边桩到中桩的距离，左“-”右“+”X、Y —待求点坐标 K—待求点里程三、数据库的编制以每个线元为一段先付参数值，例K≥542.759＝＞K≤686.148＝＞A”起点里程”=542.7590:Z终点里程=163.5:R起点半径=10^45: V终点半径=10^45:D起点坐标=22236.147:E起点坐标=32145.669:F起点切线方位角=81°22′42″:J＝0:⊿⊿K≥987.324＝＞K＜10234.021＝＞A=……J=－1⊿⊿……Prog“XBS—XY”四、纵断高程计算Lbi 0R:T:D”K(JD)”:H”H(JD)”:U”I1”:V”I2”:｛K｝V＞U=＞G=1:≠＞V＜U=＞G=-1⊿⊿K≤D-T=＞L=D-KP=H-LU/100≠＞K≤D=＞L=K-D+TP=H-U(D-K)/100+GL2/2R≠＞K≤D+T=＞L=D+T-KP=H+V(K-D)/100+GL2/2R≠＞K＞D+T=＞L=K-DP=H+LV/100⊿⊿⊿⊿P”H=”▲Goto 0H:高程、 K：变坡点里程、H：变坡点高程、I1、I2前后坡度五、中轴旋转、超高计算: (1)IR,(2)IL线路左转输2，W（Z）：中分带宽度，K：桩号，H（Z）：中桩设计高程，W（L）左侧宽度，W（R）：右侧宽度Lbi 1G”(1)IR,(2)IL”:D”K(ZH)”:S”LS1”:H”LS2”:Q”LY”:W”W(Z)”:I”I1”:J”IB”｛KVAE｝:K:V”h(Z)”:A”W(L)”:E”W(R)”G=2=＞G=-1⊿｛KPAB｝:K:P”H(Z)”:A”W(L)”:B”W(R)”Prog”XZ”:U”H(L)=”▲V”H(R)=”▲Goto11、子程序：XZG=1=＞U=A：V=B：≠＞G=-1=＞U=B：V=A⊿⊿K≤D=＞R=P-I（U-W/2）/100：E=P-I（V-W/2）/100：≠＞K≤D+S=＞Prog”ZXA”: ≠＞K≤D+S+Q=＞R=P+J(U-W/2)/100:E=P-J(V-W/2)/100: ≠＞K≤D+S+Q+H=＞Prog”ZXB”: ≠＞K＞D+S+Q+H=＞R=P-I(U-W/2)/100:E=P-I(V-W/2)/100⊿⊿⊿⊿⊿G=1=＞U=R:V=E: ≠＞G=-1=＞U=E:V=R⊿⊿2、子程序： ZXAL=K-DR=P-I(U-W/2)/100+L(U-W/2)(I+J)/100SL≤2IS/(I+J) =＞E=P-I(V-W/2)/100: ≠＞L＞2IS/(I+J) =＞E=P-LJ(V-W/2)/100S ⊿⊿3、子程序： ZXBL=D+S+Q+H-KR=P-I(U-W/2)/100+L(U-W/2)(I+J)/100HL≤2IH/(I+J) =＞E=P-I(V-W/2)/100: ≠＞L＞2IH/(I+J) =＞E=P-IJ(V-W/2)/100H ⊿⊿六：边轴旋转超高计算：L：平曲线长度，I1正常横坡，IB超高横坡。

曲线计算公式、程序及算例广西建设职业技术学院测量教研室李向民2007年4月一、圆曲线中桩与边桩坐标计算公式1．根据已知的中线偏角、半径和交点桩号计算切线长T 、曲线长L 及外距E ，以及曲线主点的桩号切线长 2tan αR T = 曲线长180παR L =外 距 )12(sec2cos-=-=ααR R R E切曲差（超距） D = 2T －L 式中α以度为单位。

交点的桩号已由中线丈量得到，根据交点的桩号和曲线测设元素，可计算出各主点的桩号，ZY 桩号 =JD 桩号－TQZ 桩号 = ZY 桩号＋2LYZ 桩号 = QZ 桩号＋2L为了避免计算中的错误，可用下式进行计算检核：JD 桩号 = YZ 桩号－T ＋D2．根据已知的两个交点的坐标，计算曲线主点的坐标和细部桩点的坐标（1）计算直圆点（ZY ）坐标根据JD 1和JD 2的坐标(x 1 , y 1)、(x 2 , y 2)，用坐标反算公式计算第一条切线的方位角α2-1和ZY 点坐标（x ZY , y ZY ），212112arctanx x y y --=-α122122sin cos --+=+=ααT y y T x x ZY ZY（2）计算圆心坐标因ZY 点至圆心方向与切线方向垂直，其方位角为αZY-O =α2-1 ± 90°（左偏角时用“+”，右偏角时用“－”）则圆心坐标（x o , y o ）为ozy ZY o o zy ZY o R y y R x x --+=+=ααsin cos（3）计算圆心至各细部点（和主点）的方位角设ZY 点至曲线上某细部里程桩点的弧长为l i ，其所对应的圆心角i β按下式计算得到：πβ180⋅=R l i i 则圆心至各细部点的方位角αi 为αi =（αZY-O ＋180°）±i β （左偏角时用“－”，右偏角时用“+”）（4）计算各细部点（和主点）的坐标根据圆心至细部点的方位角和半径，可计算细部点坐标io i i o i R y y R x x ααsin cos +=+=（5）计算各桩点边桩坐标设左边桩与中桩之间的垂直距离为D 左，右边桩与中桩之间的垂直距离为D 右，则左边桩的坐标为：io i i o i D R y y D R x x ααsin )(cos )(左左左左±+=±+=（左偏角时用“－”，右偏角时用“+”）右边桩的坐标为：io i i o i D R y y D R x x ααsin )(cos )(右右右右±+=±+=（左偏角时用“+”，右偏角时用“－”）注：可用全站仪按极坐标法测设线路主点和细部中桩，以及边桩；二、圆曲线计算程序（CASIO fx-4800P ）“YQS ”文件（主程序，计算圆曲线、直线的中桩坐标、边桩坐标）{}{}{}[][][][][][][][]{}{}2▲)90sin()(""▲)90cos()("":1:10:)sin )((:)cos )((:,211▲sin )(""▲cos )("")(180)(,11180sin cos 90sin cos 12""4:212""3:""2:180""1:)2tan(:,31:3:20""Pr :"":"":"":"",,,,2211Goto H H C B AbsH N Y V C B AbsH M X U C C H B K I P N B K I D M HK H K Lb Goto H H Q HZ R N Y V Q HZ R M X U OZ C Q R L O T I K L H K Lb J C J R F N J R E M Z B J B T P F B T D E Z Z YZ Z Z Z QZ Z T I ZY Z RA L Z A R T ZA Z A Lb Goto Goto R R OB og I Y P X D Y Y X X I P D Y X -=++=++=-=≠⇒=⇒≥-+=-+=-=++=++=⨯+=÷=--=+=+=+=⨯-=+=+=+=÷+=-=÷=÷=≠⇒⇒=ππ“OB ”文件 （子程序，计算两点间的坐标方位角）BB Lb Goto W B W B P Y Lb Goto W B Goto D X D X P Y W Lb Goto B B P Y Lb Goto Goto D X =+=≠⇒=⇒>-+=≠⇒⇒>--÷-==≠⇒=⇒>-≠⇒⇒=----:414:360tan :tan 0:314:180tan :30:)()(:214:270:900:11:2:10111注释:R 为圆曲线半径（当算直线时输入0） I 为2JD 桩号1X ,1Y 是1JD 的XY 坐标2X ,2Y 是2JD 的XY 坐标A 为偏角Z ：右偏角时Z=1，左偏角Z=－1 K 为待放点桩号H 为边桩距：当为路线左边桩时符号为正，右边桩时符号为负，求中桩时输入0。

5800程序使用说明书一、匝道放样程序（ZDJS)1.屏幕显示“K?”，输入待算点里程；2.屏幕显示“X=”，计算出待算点坐标X；3.屏幕显示“Y=”，计算出待算点坐标Y；4.屏幕显示“W?”，输入夹角；5.屏幕显示“L?”，输入左边距；6.屏幕显示“I?”，输入右边距；7.屏幕显示“XL=”，计算出待算点左边距坐标X；8.屏幕显示“YL=”，计算出待算点左边距坐标Y；9.屏幕显示“XR=”，计算出待算点右边距坐标X；10.屏幕显示“YR=”，计算出待算点右边距坐标Y。

（继续按“确认”将返回步骤1）二、匝道子程序（S1)If K≤304.758（终点桩号）:Then000.000（起点桩号）→Q:304.758（终点桩号）→Z:2762282.155（起点X坐标）→A:If K≤终点桩号：Then起点桩号→Q:终点桩号→Z：起点X坐标→A:起点Y坐标→B：起点切线方位角→F：起点曲率1/R1→D:终点曲率1/R2→R: （如是直线均输0）Goto 0:If End换行下一个线元要素输入同上。

Lb1 0：Return ↙三、坐标正反算（2-ZBFS、4-ZBZS）（一）、坐标反算1.屏幕显示“X1=?”，输入置镜点（已知点）X坐标；2.屏幕显示“Y1=?”，输入置镜点（已知点）Y坐标；3.屏幕显示“X2=?”，输入后视点（已知点）X坐标；4.屏幕显示“Y2=?”，输入后视点（已知点）Y坐标；5.屏幕显示“D=”，计算出两点之间的平距；6.屏幕显示“FWJ=”，计算出两点之间的坐标方位角；7.返回步骤3，屏幕显示“X2=?”，输入后视点X坐标。

（二）、坐标正算1.屏幕显示“X1=?”，输入置镜点（已知点）X坐标；2.屏幕显示“Y1=?”，输入置镜点（已知点）Y坐标；3.屏幕显示“D=？”，输入已知点与待算点之间的平距；4.屏幕显示“FWJ=？”，输入已知点与待算点之间的坐标方位角；5.屏幕显示“X2=”，计算出待算点X坐标；6.屏幕显示“Y2=”，计算出待算点Y坐标；7.返回步骤3，屏幕显示“D=？”，输入已知点与待算点之间的平距。

已知任意点的里程、半宽，计算该点平面坐标的程序前言1、本计算程序适用范围仅限于某ZH点至下一相邻ZH点间的线路范围。

2、使用本程序前需进行导线点坐标计算并求出各个曲线的曲线要素如第一切线方位角、转向角、曲线半径、缓和曲线长度、切线长度、曲线长度等。

3、使用过程中要注意长短链对里程传递的影响关系，否则将导致计算结果错误。

4、往计算器中输入程序符号时将⊿改为。

5、本程序仅供参考,实际使用时请先进行验证。

6、本程序适用的计算器型号为卡西欧4800/4850。

“wangDongping” (程序名称)Defm 13↙XYVWRLZTC：Fixm↙Z[4]=0.5×L－L3÷(240×R2)↙Z[2]=Z[4]+(R+L2÷24÷R)Abs tan(T÷2)：G=Z＋Ans：Fix3：G＂JDDK=＂◢Z[9]=R Abs Tπ÷180＋L：Z[9]＂QXL=＂◢Z[13]=Z＋Ans：Z[13]＂HZDK=＂◢D=C＋T：D＜0⇒D=D＋3600⊿D＞3600⇒D=D－3600⊿D◢Z[11]＂HZX＂=V＋Z[2]×(CosC＋CosD)◢Z[12]＂HZY＂=W＋Z[2]×(SinC＋SinD)◢Norm：Prog＂WEIXONG2＂(程序“WaangDongping”运行到此时，自动转到程序“WangDangping2”进行下一步运行)“WangDongping2”（程序名称）Lb1 0：{K,S,Q,E,F,N,U,H}：P=Q＋Abs(E－F)×N：K≤Z⇒Goto1⊿Z＜K⇒K≤Z＋L⇒Goto2⊿Z＋L＜K⇒K≤Z＋Z[9]÷2⇒Goto3⊿Z＋Z[9]÷2＜K⇒K＜Z[13]－L⇒Goto4⊿Z[13]－L≤K⇒K＜Z[13]⇒Goto5⊿ K≥Z[13]⇒Goto6⊿Goto0↙Lb1 1：Z[1]=V＋(Z－K)×cos(C＋1800)＋Pcos(C＋S)＋Ucos(C＋H)：Z[3]=W＋(Z －K)×sin(C＋1800)＋Psin(C＋S)+Usin(C+H)：Goto7⊿Goto0↙Lb1 2：Z[5]=K－Z－(K－Z)5÷(40R2L2)：Z[6]=（K－Z）3÷（6RL）：A=T÷AbsT：Z[7]=(K－Z)2×180÷(6πRL)：Z[1]=V＋√((Z[5])2＋(Z[6])2)×cos(C＋A×Z[7])＋Pcos(C＋3×A×Z[7]＋S)＋Ucos(C＋3×A×Z[7]＋H)：Z[3]=W＋√((Z[5])2＋(Z[6])2)×sin(C＋A×Z[7])＋Psin(C＋3×A×Z[7]＋S)＋Usin(C＋3×A×Z[7]＋H)：Goto7⊿Goto0↙Lb1 3：B=L×180÷(2πR)：Z[8]=B＋(K－Z－L)×180÷(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(K－Z)－L)×180÷(2πR)：A=T÷AbsT：Z[1]=V+ √((Z[5])2+(Z[6])2)×cos(C+A×Z[7])+Pcos(C+A×Z[10]+S)＋Ucos(C＋A×Z[10]＋H)：Z[3]=W+√((Z[5])2＋(Z[6])2)×sin(C+A×Z[7])+Psin(C+A×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 4：B=L×180÷(2πR)：Z[8]=B+(Z[13]－K－L)×180(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(Z[13]－K)－L)×180÷(2πR)：A=－T÷AbsT：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+A×Z[10]+S)＋Ucos(D＋A×Z[10]＋H)：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+A ×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 5：Z[5]=Z[13]－K－(Z[13]－K) 5÷(40R2L2)：Z[6]=(Z[13]－K) 3÷(6RL)：A=－T÷AbsT：Z[7]=(Z[13]－K)2×180÷(6πRL)：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+3×A×Z[7]+S)+Ucos（D+3×A ×Z[7]+H）：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+3×A×Z[7]+S)＋Usin(D+3×A×Z[7]+H)：Goto7⊿Goto0↙Lb1 6：Z[1]=Z[11]+(K－Z[13])×CosD+Pcos(D+S)＋Ucos(D+H)：Z[3]=Z[12]+(K －Z[13])×sinD+Psin(D+S)＋Usin(D+H)：Goto7⊿Goto0↙Lb1 7：M=Z[1]－X：O=Z[3]－Y：Pol(M,O)：J＜0 J=J+3600⊿Fix3：J◢I◢P◢E－F◢M+X◢O+Y◢Norm：Goto0↙坐标程序符号说明1、X，Y：测站坐标2、V，W：ZH点X,Y坐标3、R：曲线半径4、L：缓和曲线长度5、Z：ZH点里程6、T：曲线转向角（左负右正）7、C：第一切线方位角8、K：放样点里程9、Q：放样半宽（平行于S角非切线边在里程点左右侧的宽度，若已知实际半宽P，则Q=P；若放路基边桩则输入变坡点处的半宽P）10、S：放样点与其对应的中线点连线和该中线点切线的夹角（左负右正，即放样点在线路左侧S为负，同理在右侧为正）11、F：放样点地面标高12、N：路基边坡坡度13、E：变坡点设计标高14、P=Q+Abs（E－F）×N（Abs为绝对值），若P为已知，则Q的值直接输入P的值即可，E、F、N均为零。

一、QXFY 辛甫森公式放样程序1. U“X0”：V“Y0”2. A“XA”：B“Y A”：C“CA”：D“1÷RA”：E“1÷RB”：F“KA”：G“KB”3. Lb1 1：{HLR}：H“KI”：L“JJ”：R“DZ+Y-”4. H>G Goto 15. H<F Goto 16. P=(E-D) ÷Abs(G-F)：Q=Abs(H –F)：Z=P×Q：T=D+Z7. W=C+（Z+2D）Q×90÷π8. M=C+(Z÷4+2D)Q×22.5÷π：N=C+（3Z÷4+2D）Q×67.5÷π：K=C+(Z÷2+2D)Q×45÷π9. X=A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W-L)10. Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinW) ÷12+ Rsin (W-L)11. Pol((X-U)，（Y-V）)J≤0J=J+360J“FWJ”12. I“S”13. Goto 1程序显示说明：须输入参数：X0:置镜点X (对应变量U)Y0:置镜点Y (对应变量V)XA:曲线起算点X (对应变量A)Y A:曲线起算点Y (对应变量B)CA：曲线起算点切线方位角(对应变量C) 1÷RA：1÷半径，即起算点曲率，右偏为正左偏为负(对应变量D)1÷RB：1÷半径，即曲终点曲率，右偏为正左偏为负(对应变量E)KA:起算点里程(对应变量F)KB:曲终点里程(对应变量G)KI:待求点里程(对应变量H)JJ:夹角(与前进方向切线方位角夹角，向右为正)(对应变量L)D“Z+Y－”：偏距，左正右负(对应变量R)计算结果: W: 待求点切线方位角(对应变量W) X:计算点x (对应变量X)Y:计算点y (对应变量Y)FWJ:放样方位角(对应变量J)S:放样距离(对应变量I)K BX 0Y 0本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R ）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩及斜交或正交的边桩。