曲線任意里程中邊樁坐標正反算.docVIP

?曲線任意里程中邊樁坐標正反算5800p TYQXjs(主程序名)LbI?0:4Dimz: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:XI:Y→JProgSUB2:O+W→S:ClsS=:S◢Z=:Z◢Goto?2注:S=0=Goto?0和X=0=Goto0?可以不要,亦可以要,它的作用是將S或者X輸入0時,程序會重新由頭開始運行,(即是可以計算下段線元)SUB1(子程序名)0.1739274226A: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-90T(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?0Z：?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注意:本程序是根據yshf編寫的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=XYYS=Y◢FFS=F-90◢Goto?1←┘Lbl?2：{XY}：XY：I=X：J=Y：Prog?SUB2：SS=O+WZZ=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