5800P曲线坐标计算编程使用(非标准曲线).docVIP

曲线任意里程中边桩坐标正反算(CASIO fx-800P计算器)程序 曲线任意里程中边桩坐标正反算(CASIO fx-800P计算器)程序一、程序功能????????本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-0P计算器上运行。?? 特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；????????????????????(2). 在引用该核心计算部分时，请注明来源。二、源程序主程序(QXJS)1.SZ = XY：2.XY = SZN:X0?U: Y0?V: S0?O: FWJ0?G: LS?H:R0?P?RS?R: PX?Q： 1÷PC： (P-R)÷(2HPR)D： 180÷πE：N=1:Then Goto 1：Goto 2:IfEndLbl 1: KS?Z:ES?S：Abs(-O) →W：Prog :XS= :X→X◢YS=: Y→Y◢FWJS=:(F-90)DMS◢ Prog GCGoto1 Lbl 2： ?X:YS?Y： XI： YJ：Prog :O+W→Z KS=:Z→Z◢ ES=: S→S◢ FWJS=:(F-90)DMS◢ Prog GCGoto 2 1.全线计算主程序() 1.SZ = XY：2.XY = SZ：?N: Lbl 0：?Z:ES?S：Prog :1÷P→C:(P-R)÷(2HPR) →D:180÷π→E：N=1:Then Goto 1: Else Goto 2:IfEndLbl 1:Abs(Z-O) →W：Prog :XS= :X→X◢YS=: Y→Y◢FWJS=:(F-90)DMS◢ Prog GCGoto 0 Lbl 2： ?X:YS?Y： XI： YJ：Prog :O+W→Z KS=:Z→Z◢ ES=: S→S◢ FWJS=:(F-90)DMS◢ Prog GCGoto 22. 正算子程序()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+Scos(F→X:Y+Ssin(F→Y Return3. 反算子程序()G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T))→W:0→S：Lbl 0：Prog :T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→S:If Abs(S)10^(-6):Then Goto1:Else S+W→W:Goto 0:IfEndLbl 1:0→S：Prog :(J-Y)÷sin(F)→S Return数据库：XXX Lbl 1:If Z≤下一起点里程: Then本线路起点里程→O：U：V:本线路方位角→G：H：P：R:±1→Q (左偏Q=-1，右偏Q=1，直线段Q=0）O：U：V:本线路方位角→G：H：P：R:±1→Q :IfEnd: IfEnd Return 四、数据库 1、线路平面数据库 数据名：“XXX”可根据自己须要更改 Lbl 1:If Z≤80(计算段终点里程) : Then 0→O (计算段起点里程): 62450.391→U(起点里程纵坐标): 56853.282→V (起点里程横坐标): 220°31′33″→G(计算段起点里程切线方位角): 80→H(缓和线): 10^45→P曲率半曲率半曲率半曲率半 2、数据库数据输入时应注意缓和曲线段、曲线段、两曲线中有一直线段为一计算段，详上面数据库内示意数据，带下划画线的为一计算段，逐个往下编直至整条线路编完，但注意的是CASIO fx-4850P计算器 GC=: E+(A+B)V+C→K◢ Return 6．道路放样高程主程序（DLCP）道路抄平 注：该