曲线任意里程中边桩坐标正反算5800p好用.doc

曲線任意里程中邊樁坐標正反算5800pTYQXJS(主程序名)LbI 0:4→Dimz:Cls1.SZ→XY2.XY→SZ?→N:ClsXO?U: YO?V: SO?O:FO?G: LS?H: RO?P:RN?R:?Q1÷P→C：(P-R)÷(2HPR) →D：180÷π→EN=2=Goto 2LbI 1:Cls:?SS=0=Goto 0??(注:當S輸入0時,程序會由LbI 0開始重新運行)?Z:Abs(S-O)→WProgSUB1:ClsXS=:X◢YS=:Y◢FS=:F-90 ?DMS◢Goto 1LbI 2:Cls:?XX=0=Goto 0??(注:當X輸入0時,程序由會LbI 0開始重新運行)?Y:X→I:Y→JProgSUB2:O+W→S:ClsS=:S◢Z=:Z◢Goto 2注:紅色部可以不要,亦可以要,它的作用是將S或X輸入0時,程序會重新由頭開始運行,(即是可以計算下段線元)SUB1(子程序名)0.1739274226→A:0.3260725774→B0.0694318442→K:0.3300094782→L1-L→F:1-K→MG+QEKW(C+KWD)→Z[1]G+QELW(C+LWD)→Z[2]G+QEFW(C+FWD)→Z[3]G+QEMW(C+MWD)→Z[4]A×cos(Z[1])→X??(剩號可省略)X+Bcos(Z[2])→XX+Bcos(Z[3])→XX+Acos(Z[4])→XU+WX→XAsin(Z[1])→YY+Bsin(Z[2])→YY+Bsin(Z[3])→YY+Asin(Z[4])→YV+WY→YG+QEW(C+WD)+90→FX+Zcos(F)→XY+Zsin(F)→YSUB2(子程序名)G-90→T(Y-V)cos(T)-(X-U)sin(T)→WAbs(W)→W:0→ZLbI 4:ProgSUB1T+QEW(C+WD)→L(J-Y)cos(L)-(I-X)sin(L)→ZIf Abs(Z)＜1E-6:Then 0→Z :ProgSUB1(J-Y)÷sin(F)→ZElse W+Z→W:Goto 4:IfEnd注:1E-6即是10的負6次方,亦即是等於0.000001輸入時按SHIFT---log---(-)---6亦可直接將If Abs(Z)＜1E-6寫成If Abs(Z)＜0.000001注意:本程序是根據yhhf編寫的4800P改寫而成,以下是源程序及算例大家有問題亦可在以下網址查看及向yshf發問/celiang/article_view.asp?id=380author=yshf曲线任意里程中边桩坐标正反算(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：RRN：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 22. 正算子程序(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)))