曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序.doc

《曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、曲线任意里程中边桩坐标正反算(CASIOfx-5800P计算器)程序一、程序功能      本程序由一个主程序(TYQXJS)和三个子程（正算子程序(SUB1)、反算子程序(SUB2)，数据库子程序（DAT03））构成，可以根据曲线段—直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行夹角为90°的正反算。程序改进：1、增加桩号控制语句，当计算桩号超出本线元自动从数据库读取线元数据，或者提示手工输入。      

2、    2、手工输入下一线元数据时，程序将上一线元终点桩号、坐标、方位角自动代入下一线元作为线元起点数据进行计算。      特别申明：(1).适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；        (2).本程序采用yshf老师4800P程序TYQXJS中核心算法，在此表示感谢，转载请注明！        (3).程序还有不足之处，（如：坐标反算不能全线贯通，只能算本线元内，超出本线元之后需重新读取数据库，然后重新输入反算点坐标），请高手指点！二、程序清单：TYQXJS(曲线坐标计算主程序)“1.SJK，2.SHURU”

3、?N:IfN=1:ThenGoto1:ElseGoto2:IfEndLbl1:“Z”?S:Prog“DAT03”:Goto3Lbl2:“X0”?U:“Y0”?V:“Z0”?O:“F0”?G:“LS”?H:“R0”?P:“RN”?R:?Q:(P-R)÷(2HPR)→D:180÷兀→ELbl3:“1.ZS,2.FS”?C:IfC=1:ThenGoto4:ElseGoto5:IfEndLbl4:“Z”?SIfS＜OOrS＞O+H:ThenIfN=1:ThenProg“DAT03”:Goto4:ElseGoto6:IfEnd:IfEnd“W:Z-,Y

4、+”?ZAbs(S-O)→WProg“SUB1”“FWJ=”:F-90◢DMS◢“X=”:X◢“Y=”:Y◢Goto4Lbl5:?X:?Y:X=I:Y=J:Prog“SUB2”IfW＜0OrW＞H:ThenIfN=1:ThenO+W→S:Prog“DAT03”:Goto5:ElseGoto6:IfEnd:IfEnd“Z=”:O+W→S◢“W=”:Z◢Lbl6:IfS＜OOrW＜0:Then“OUT-OF-RANGE”◢Goto2:IfEndH→W:0→Z:Prog“SUB1”:X→U:Y→V:F-90→G:O+H→O:“LS”?H:“R0”?

5、P:“RN”?R:?Q:(P-R)÷(2HPR)→D:180÷兀→EIfC=1:ThenGoto4:ElseGoto5:IfEndSUB1(正算子程序)0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→MU+W(Acos(G+QEKW(1/P+KWD))+Bcos(G+QELW(1/P+LWD))+Bcos(G+QEFW(1/P+FWD))+Acos(G+QEMW(1/P+MWD)))→XV+W(Asin(G+QEKW(1/P+KWD))+Bsin(G+

6、QELW(1/P+LWD))+Bsin(G+QEFW(1/P+FWD))+Asin(G+QEMW(1/P+MWD)))→YG+QEW(1/P+WD)+90→FX+ZcosF→XY+ZsinF→YReturnSUB2(反算子程序)G-90→T：Abs((Y-V)cosT-(X-U)sinT)→W：0→ZLbl0:Prog"SUB1"T+QEW(C+WD)→L:(J-Y)cosL-(I-X)sinL→Z:IfAbsZ<1×10-4:ThenGoto1:ElseW+Z→W：Goto0:IfEndLbl1:0→Z:Prog"SUB1":(J-Y)÷s

7、inF→ZReturn说明：1、×10-4：×10就是最下面一排键“.”键右边的一个键。-4为负4    2、程序中的字母O都已经加粗，区别于数字零DAT03(数据库子程序)IfS＜路线起点桩号OrS＞路线终点桩号:Then2→N:“X0”?U:“Y0”?V:“Z0”?O:“F0”?G:“LS”?H:“R0”?P:“RN”?R:?QGoto1:IfEndIfS＜线元终点桩号:Then线元起点X坐标→U:线元起点Y坐标→V:线元起点桩号→O:线元起点方位角→G:线元长度→H：线元起点半径→P：线元终点半径→R:曲线偏向(1,0,-1)→Q:Go

8、to1:IfEndIfS＜线元终点桩号:Then线元起点X坐标→U:线元起点Y坐标→V:线元起点桩号→O:线元起点方位角→G:线元长度→H：线元起点半径→P：线元终