曲线任意坐标正反算5800p好用

1、曲线任意坐标正反算5800p| TYQXJS主程序名)Lbl 0:4 f Dimz:Cls"1.SZ f XY""2.XY f SZ"?f N:Cls"XO"?U: "YO"?V: "SO"? O"FO"?G: "LS"?H: "RO"?P:"RN"?R:?Q1-Pf C: (P-R) - (2HPR) f D: 180-nf EN=2=>Goto 2LbI 1:Cls:?S:S=0=>Goto 0?Z:

2、Abs(S- C) f WProg"SUB1":Cls"XS=":X 丄"YS=":Y 丄"FS=":F-90? DMMGoto 1Lbl 2:Cls:?XX=0=>Goto 0?Y:Xf I:Y f JProg"SUB2": O+Wf S:Cls"S=":S 丄"Z=":Z 丄Goto 2注:紅色部可以不要，亦可以要，它的作用是將S或 X輸入0時,程序會重新由頭 開始運行，(即是可以計算下段線元)SUB1子程序名)0.1739274226 f A:

3、0.3260725774 f B0.0694318442 f K:0.330009482 f L1-L f F:1-K f MG+QEKW(C+KWDZ1 G+QELW(C+LWfZ2 G+QEFW(C+FWDZ 3 G+QEMW(C+MWD(4 AX cos(Z1) f X X+Bcos(Z2) f X X+Bcos(Z3) f X X+Acos(Z4) f X U+Wf XAsi n( Z1) f YY+Bsi n(Z2) f YY+Bsi n( Z3) f YY+Asi n(Z4) f Y V+Wf YG+QEW(C+WD)+fFX+Zcos(F) f X Y+Zsin(F) f YSU

4、B2子程序名)G-90f T(Y-V)cos(T)-(X-U)s in (T)f WAbs(W)f W:(H ZLbI 4:Prog"SUB1"T+QEW(C+WD)L(J-Y)cos(L)-(l-X)si n(L) f ZIf Abs(Z) v 1E-6:Then 0 f Z :Prog"SUB1":(J-Y) - sin(F) f ZElse W+Zf W:Goto 4:IfE nd数据库If S w终点桩号：Thenlf N:。 U:起点X坐标f A:起点Y坐标f B:起点切线方位角f C:起点曲率f D:终点曲率f E: Goto 0:lf En

5、d 换行冋上下个线兀要素。0000000000000000000000000000000000000000N个Lb1 0 : Return 换行曲率左偏为负右偏为正Return为程序命令线元要素可以用常规程序计算如有多条路数据库另建如Z2 Z3 0 0 0等等需要使用时把主程序 Z1改为Z2依次类推NO？为要求点里程注:1E-6即是10的負6次方，亦即是等於 0.000001輸入時按 SHIFT-log-(-)-6亦可直接將 If Abs(Z) v 1E-6 寫成 If Abs(Z) v 0.000001注意:本程序是根據yhhf編寫的4800P改寫而成，以下是源程序及算例大家有問題亦可在以下

6、網址查看及向yshf發問 ng.n et/celia ng/article_view.asp?id=380&author=yshf曲线任意里程中边桩坐标正反算 (CASIO fx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程一一正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线(完整或非完整型) 的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止 点曲率半径)及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，

7、用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。特别申明：(1).适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm(2) .在引用该核心计算部分时，请 注明来源。二、源程序1. 主程序(TYQXJS)"1.SZ => XY" : "2.XY => SZ" : N： U"X0": V"Y0": OSO" : G"F0": H"LS": P"R0"

8、;:R"RN": Q C=1- P: D=(P-R) - (2HPR): E=18On: N=1=>Goto 1:工 >Goto 2 ZjLbl 1 : SZ : SZ W=Abs(SQ : Prog "SUB1" : X"XS"=XJY"YS"=YJF"FS"=F-90 丄Goto 1 Lbl 2 : XY: XY I=X : J=Y : Prog "SUB2" : S"S"=6WJZ"Z"=Z 丄Goto 22.正算子程

9、序(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(Asi n(G+QEKW(C+KWD)+Bsi n(G+QELW(C+LWD)+Bsi n(G+QEFW(C+FWD)+Asi n( G+QEMW(C+MWD=)G+QEW(C+WD)+90 X=X+ZcosF Y=Y+ZsinF3.反算子程序(S

10、UB2)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-(l-X)sinL: AbsZ<1E-6=>Goto1:工 >W=W+ZGoto 0Lbl 1 : Z=0: Prog "SUB1" : Z=(J-Y) - sinF三、使用说明1、规定(1) 以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时，Q=-1;当线元往右偏时，Q=1;当线元为直线时，Q=(。(2) 当所求点位于中线时，Z=0;当

11、位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半于设计规

12、定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输入部分：1. SZ => XY2. XY = > SZN?选择计算方式，输入1表示进行由里程、边距计算坐标输入2表示由坐标反算里程和边距。X0 ?线元起点的X坐标Y0 ?线元起点的Y坐标SO ?线元起点里程F0 ?线元起点切线方位角LS ?线元长度R0 ?线元起点曲率半径RN ?线元止点曲率半径Q ? 线元左右偏标志（左偏Q=-1，右偏Q=1,直线段Q=0S?正算时所求点的里程Z ?正算时所求点距中线的边距 （左侧取负，值右侧取正值，在中线 上取零）X ?反算时所求点的X坐标Y ?反算时所求点的Y坐标显示部分：X

13、S=xxx正算时，计算得出的所求点的X坐标YS=xxx正算时，计算得出的所求点的Y坐标FS= xxx正算时，所求点对应的中线点的切线方位角S= xxx反算时，计算得出的所求点的里程Z= xxx反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线 +完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程 SO起点坐标XO Y0、起点切线方位角F0、线元长度LS、起点曲率半径 R0止点曲率半径 RN线元左右偏标志Q）如下:SOX0Y0F0LS R0 RNQ500.00019942.83728343.561125 16 31.00269.2561E45 1E450769

14、.25619787.34028563.378125 16 31.0037.4921E45 221.75-1806.74819766.56628594.574120 25 54.07112.779221.75 221.75-1919.52719736.07228701.89391 17 30.6380.285221.75 9579.228-1999.81219744.03828781.65980 40 50.00100.0001E45 1E4501、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785YS=28509

15、.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z= 5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS=19785.25749YS=28575.02270S=780Z=0计算得XS=19781.15561YS=28572.16358S=780Z= 5计算得XS=19777.05373YS=28569.304469YS=28654.130918YS=28652.913798YS=28651.696689118YS=28722.05802得 XS=19736.47687298YS=28722.53168

16、S=870Z=-5S=870Z=0S=870 Z= 5S=940Z=-5.123S=940Z=0YS=28722.35642S=940 Z= 3.009计算得XS=19747.5360计算得XS=19742.6864计算得XS=19737.8368计算得XS=19741.5计算计算得XS=19733.472、反算X=19831.41899.9999974Z= -5 .00018164X=19827.33699.9996493Z= 0.000145136X=19823.25398得 S=699.9999985 Z= 5.000003137X=19785.257490.0000035Z= -5 .000001663X=19781.1556180.0000025Z=- 0.000002979X=19777.0537380.0000016Z= 4.99999578Y=28509.726计算得S=6Y=28506.838计算得S=6Y=28503.95084计算Y=28575.02270计算得S=78Y=28572.16358计算得S=7Y=28569.30446计算得S=7X=19747.536Y=28654.131计算得 S=870.0001137 Z= -4.99941049X=19742.686Y=28652