曲线任意坐标正反算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 nEN=2=>Goto 2LbI 1:Cls:?S:S=0=>Goto 0?Z

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

3、f A:0.3260725774 f B0.0694318442 f K:0.330009482 f L1-L f F:1-K f MG+QEKW(C+KWD)f Z1G+QELW(C+LWD)f Z2G+QEFW(C+FWD)f Z3G+QEMW(C+MWD) f Z4A Xcos(Z1) f XX+Bcos(Z2) f XX+Bcos(Z3) f XX+Acos(Z4) f XU+WX f XAsin(Z1) f YY+Bsin(Z2) f YY+Bsin(Z3) f YY+Asi n(Z4) f YV+WY f YG+QEW(C+WD)+90f FX+Zcos(F) f XY+Zsin(

4、F) f YSUB2 (子程序名)G-90 f T(Y-V)cos(T)-(X-U)s in (T)f WAbs(W) f W:0 f ZLbl 4:Prog"SUB1"T+QEW(C+WD) f L(J-Y)cos(L)-(I-X)si n(L) f ZIf Abs(Z) v 1E-6:Then 0 f Z :Prog"SUB1":(J-Y) *sin(F) f ZElse W+Z f W:Goto 4:lfEnd数据库If SW终点桩号：Then1 f N: °°f U:起点X坐标f A:起点Y坐标f B:起点切线方位角f C:起

5、点曲率f D:终点曲率f E: Goto 0:If End换行| ij I . |o oooooooooo000000000000000000000000000000N个Lb1 0: Return 换行曲率左偏为负右偏为正Return为程序命令线元要素可以用常规程序计算如有多条路数据库另建如Z2 Z3 °。等等需要使用时把主程序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編寫

6、的4800P改寫而成，以下是源程序及算例大家有問題亦可在以下網址查看及向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" : O"S0" : G"F

8、0" : H"LS":P"R0" : R"RN" : Q : C=1 讦：D=（P-R） *（2HPR） : E=180 *n:N=1=>Goto 1 :工Goto2 A Lbl 1 : SZ: SZ: W=Abs(S- O): Prog "SUB1" : X"XS"=X 丄Y"YS"=Y 丄F"FS"=F-90 丄Goto 1 Lbl 2 : XY: XY: I=X : J=Y : Prog "SUB2" : S&qu

9、ot;S"= O+W 丄Z"Z"=Z 丄Goto 22. 正算子程序(SUB1)A=0.1739274226: B=0.3260725774: K=0.0694318442: L=0.3300094782F=1-L :M=1-K : X=U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bco s(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+M

10、WD): F=G+QEW(C+WD)+90 : X=X+ZcosF : Y=Y+ZsinF3. 反算子程序(SUB2)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+Z : Goto 0 ALbl 1 : Z=0 : Prog "SUB1" : Z=(J-Y) FinF三、使用说明1、规定(1) 以道路中线的前进方向(即里程增大的方

11、向)区分左右；当线元往左偏时，Q=-1 ;当线元往右偏时，Q=1 ;当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0 ;当位于中线左铡时，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=0)S ?正算时所

13、求点的里程Z ?正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零 ）X ?反算时所求点的X坐标Y ?反算时所求点的Y坐标显示部分：XS= xx>（正算时，计算得出的所求点的X坐标YS= xx>（正算时，计算得出的所求点的Y坐标FS= xxx:正算时，所求点对应的中线点的切线方位角S= xxx反算时，计算得出的所求点的里程Z= xxx反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线） 组成，各段线元的要素（起点里程 SO、起点坐标X0 Y0、起点切线方位角F0、线 元长度LS、起点曲率半径 R0、止点曲率半径 R

14、N、线 元左右偏标志Q）如下：SO X0Y0F0LS R0 RNQ500.00019942.83728343.561125 16 31.00269.2561E451E450769.25619787.34028563.378125 16 31.0037.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.81219744.03828781.659804050.00100.0001E451E

15、4501、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785YS=28509.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.30446S=870Z=-5计算得XS=1974

16、7.53609YS=28654.13091S=870Z=0计算得XS=19742.68648YS=28652.91379S=870Z=5计算得XS=19737.83688YS=28651.69668S=940 Z=-5.123计算得 XS=19741.59118 YS=28722.05802S=940Z=0计算得XS=19736.47687YS=28722.35642S=940Z=3.009计算得XS=19733.47298YS=28722.531682、反算X=19831.418Y=28509.726计算得 S=699.9999974 Z=-5 .00018164X=19827.336Y=2

17、8506.838计算得 S=699.9996493 Z=0.000145136X=19823.25398Y=28503.95084计算得1 S=699.9999985 Z=5.000003137X=19785.25749Y=28575.02270计算得S=780.0000035 Z=-5 .000001663X=19781.15561Y=28572.16358计算得S=780.0000025 Z=-0000002979X=19777.05373Y=28569.30446计算得S=780.0000016 Z=4.99999578X=19747.536Y=28654.131计算得 S=870.0001137 Z=-4.99941049X=19742.686