5800线路坐标正反算程序.docVIP

曲線任意里程中边桩正反算5800p TYQXJS(主程序名) LbI 0:4→Dimz:Cls 1.SZ→XY 2.XY→SZ?→N:Cls XO?U: YO?V: SO?O:FO?G: LS?H: RO?P:RN?R:?Q 1÷P→C：(P-R)÷(2HPR) →D：180÷π→E N=2=Goto 2 LbI 1:Cls:?S ：?Z:Abs(S-O)→W ProgSUB1:Cls XS=:X◢ YS=:Y◢ FS=:F-90 ?DMS◢ Goto 1 LbI 2:Cls:?X ：?Y:X→I:Y→J ProgSUB2:O+W→S:Cls S=:S◢ Z=:Z◢ Goto 2 SUB1(子程序名) 0.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+ZcosF →X： Y+ZsinF→ Y SUB2(子程序名) G-90→T (Y-V)cos(T)-(X-U)sin(T)→W Abs(W)→W:0→Z LbI 4:ProgSUB1 T+QEW(C+WD)→L (J-Y)cos(L)-(I-X)sin(L)→Z If Abs(Z)＜1E-6:Then 0→Z :ProgSUB1 (J-Y)÷sin(F)→Z Else W+Z→W:Goto 4:IfEnd 注:1E-6即是10的負6次方,亦即是等於0.000001 曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 一、程序功能 本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序( SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。 特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优 于1mm； (2). 在引用该核心计算部分时，请注明来源。 二、源程序 1.主程序(TYQXJS) 1.SZ = XY：2.XY = SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：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：XXS=X◢ YYS=Y◢ FFS=F-90◢ Goto 1←┘ Lbl 2：{XY}：XY：I=X：J=Y：Prog SUB2：SS=O+W◢ ZZ=Z◢ Goto 2 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：AbsZ1E-6=Goto1：≠W=W+Z：Goto 0Δ←┘ Lbl 1：Z=0：Prog SUB1：Z=(J-Y)÷sinF 三、使用说明 1、规定 (1) 以道路中线的前进方向（即里