《我的5800测量程序》由会员分享，可在线阅读，更多相关《我的5800测量程序（3页珍藏版）》请在金锄头文库上搜索。
1、今天终于完成了缓和曲线上由坐标反算里程的程式，误差小于 0.0001m。在 5800 上运行，大约需 要 1 分钟左右。此程式逐步逼近法，所以费时。下面我把程式写出来。 Prog “INIT”-在这个程序中输入曲线参数。 “dka”?A:”ca”:?B:”X0”?C:”Y0”?D:”Ra”?E:”Rb”?F:”dkb”?G #Ra、Rb 分别为起点和终点的曲率，左曲线为负，右曲线为正，直线为 0 #ca 为起点方位角 Prog”FFXPS”-复方辛普森公式，分段数为 9If KG:then G-AO:else K-AO:ifend:(F-E)(G-A)P:E+POQ:#Q 为 i 点的曲率 B
2、+（Q+E）*O*90Z:prog”K”:prog”K05”:#Z 为 i 点的方位角C+H6*（cos（B）+2S+4U+cos（Z） ）X:#B 为起始点方位角 D+H6*(sin(B) +2T+4V+sin(Z)Y:#H 为分为 N 小段后每段长度=（dkb-dka）/9 If L0:then prog”BZJS”:ifend #L 为边桩的宽度，左边桩为负，右边桩为正，0 为中桩 Prog”K”-子程序，求出 S 和 T0S:0T:O9H:1I:#O=dki-dka Lbi 1:A+IHJ:E+P(J-A)X:B+(X+E) *(J-A)*90Y:#X 为 k 点曲率，Y 为 k 点方
3、位角 S+cos(Y)S:T+sin(Y)T:Isz I:#Isz I 是将变量 I 加 1 If I180:then Abs(L-K-360)M:else Abs(L-K)M:ifend:Icos(M)S:A+SK: Lbi 1:0L:Prog”FFXPS”:pol(W-X,R-Y): If J180: then Abs(J-Z-360)M: else Abs(J-Z)M:ifend:Icos(M)S:K+SK: If Abs(S)0.0001: then Goto 1: ifend:#0.0001 决定了所求里程的误差大小。If J-Z180 OR J-Z360+JR:JR: Lbi 1:
4、”IX”?P:”IY”?Q:pol(P-K,Q-L)S: J360+JU:JU:if Abs(U-R)0: then “YP”:W Else “ZP”:W Ifend:Goto 1Prog”XYTOSZ”-坐标反算里程（圆曲线） “Xs”?A:”Ys”?B:”Ca”?C:”FX”?F:”bj”?D:”dka”?E: A+Dcos(C+90F)M:B+Dsin(C+90F)N: Lbi 1:”X”?X:”Y”?Y:pol(X-M,Y-N):J360+JJ:C+90F+180O: J+180P:if P-O180:then Abs(P-O-360)W:else Abs(P-O)W:ifend:D*
5、W*/180S:”DK=”:E+S “BZK”:(D-I)*FProg”SQXJS”-竖曲线计算程序C-IM:C+JN:if KA and KM: then L-(C-K)*FH:ifend: If KM and KC: then L-(C-K)*FH:H+(sqrt(R2+(K-M)2)-R)*E/Abs(E)H:ifend: If KC and KN: then L+(K-C)*GH:H+(sqrt(R2+(N-K)2)-R)*E/Abs(E)H:ifend: if KNand KB: then L+(K-C)*GH:ifend: A-dks 起始里程 B-dke 终点里程 C-dkj 变坡点里程 R-曲线半径 E-变坡点外失距，凸曲线为负，凹曲线为正 F-第一坡度 G-第二坡度（依里程前进方向） I-第一切线长 J-第二切线长 M-竖曲线开始里程 N-竖曲线结束里程 L-变坡点高程 K-待求点里程 H-待求点高程R坡度2坡度1dkeNEdkjMdksprog”QXYS”:Lbi 1:”DKI”?K:Prog”SQXJS”:”H”:H Goto 1 具体应用时，可以把每一段曲线要素写成一人子程序，然后要算哪一段在这个程序中调用就行了， 程序看起来会很短而且不复杂。
