卡西歐計(jì)算器坐標(biāo)的正反算

1、線路中邊樁任意里程坐標(biāo)正、反算改編程序（CASIO fx-4800）可以算任意斜交涵洞軸線的坐標(biāo)，增加T為斜交角度，規(guī)定T為涵軸右側(cè)方向與“線路前進(jìn)方向切線”之間的夾角，當(dāng)涵軸與線路正交時(shí)，T=90,其他操作與原程序一樣；  1. 正算子程序(SUB1)                            

2、                            color=RedA=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+57.2958QKW(1/P+KWD)+Bcos(G+57.2958QLW(1/P+L

3、WD)+Bcos(G+57.2958QFW(1/P+FWD)+Acos(G+57.2958QMW(1/P+MWD)：Y=V+W(Asin(G+57.2958QKW(1/P+KWD)+Bsin(G+57.2958QLW(1/P+LWD)+Bsin(G+57.2958QFW(1/P+FWD)+Asin(G+57.2958QMW(1/P+MWD)：F=G+57.2958QW(1/P+WD)+90：X=X+Zcos(F-90+T)：Y=Y+Zsin(F-90+T)2. 反算子程序(SUB2)  W=Abs(Y-V)cos(G-90)-(X-U)sin(G-90)：Z=0：Lbl 0

4、：Prog "SUB1"：L=(G-90)+57.2958QW(1/P+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：>W=W+Z：Goto 0Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF二.增設(shè)數(shù)據(jù)庫(kù)程序（主程序）     Lb1 4："1.SZ => XY"："2.XY => SZ"：NS：S下一線元起點(diǎn)里程O=本線元起點(diǎn)里程：U=本線元起點(diǎn)X：V=本線元起點(diǎn)Y：G=本

5、線元起算方位角：H=本線元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、-1：Prog“TYQXJS”： Goto0(第一線元數(shù)據(jù)要素)S下一線元起點(diǎn)里程O=本線元起點(diǎn)里程：U=本線元起點(diǎn)X：V=本線元起點(diǎn)Y：G=本線元起算方位角：H=本線元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、-1： Goto0(第二線元數(shù)據(jù)要素)S下一線元起點(diǎn)里程O=本線元起點(diǎn)里程：U=本線元起點(diǎn)X：V=本線元起點(diǎn)Y：G=本線元起算方位角：H=本線元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、-1： Goto0(第三線元數(shù)據(jù)要素)。Goto0(第N-1線元數(shù)據(jù)要素)。Goto0(第N線元數(shù)

6、據(jù)要素)LB1 0：D=(P-R)÷(2HPR)：N=1=>Goto 1：>Goto 2Lbl 1：TZ：Z：W=Abs(S-O)：Prog "SUB1"：X"XS"=XY"YS"=YF"FS"=F-90C“XC”：E“YC”：I=0：J=0：Pol（X-C，Y-E）：I“I=” J0J“J=”+360J“J=” Goto4Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto4 

7、;                                                  

8、;    三. 使用說(shuō)明：  1.把所有相關(guān)的”線元要素“依次輸入”SJK“  2.運(yùn)算時(shí)直接調(diào)用“SJK”運(yùn)行，程序提示輸入里程“S”？時(shí)，正算直接輸入待求點(diǎn)里程，反算輸入所求點(diǎn)“近似“里程  3.如果在正算中途想進(jìn)行反算，重新給”N“？賦值"2"便可4.程序中“XC、YC、”為測(cè)站坐標(biāo)，“”、“”為放樣距離和方位角  5.其它輸入與"yshf"原程序一致（程序中1/P表示1÷P）曲線任意里程中邊樁坐標(biāo)正反算(CASIO fx-480

9、0P計(jì)算器)程序一、程序功能        本程序由一個(gè)主程序(TYQXJS)和兩個(gè)子程正算子程序(SUB1)、反算子程序(SUB2)序構(gòu)成，可以根據(jù)曲線段直線、圓曲線、緩和曲線（完整或非完整型）的線元要素（起點(diǎn)坐標(biāo)、起點(diǎn)里程、起點(diǎn)切線方位角、線元長(zhǎng)度、起點(diǎn)曲率半徑、止點(diǎn)曲率半徑）及里程邊距或坐標(biāo)，對(duì)該曲線段范圍內(nèi)任意里程中邊樁坐標(biāo)進(jìn)行正反算。另外也可以將本程序中核心算法部分的兩個(gè)子程序移植到其它相關(guān)的程序中，用于對(duì)曲線任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序也可以在CASIO fx-4500P計(jì)算器及 CASIO f

10、x-4850P計(jì)算器上運(yùn)行。   特別申明：(1). 適用于弧長(zhǎng)小于2倍半徑的各種線元坐標(biāo)正反算，精度優(yōu)                 于1mm；                二、源程序1.主程序(TYQXJS)"1.SZ => XY&qu

11、ot;："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷：N=1=>Goto 1：>Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=XY"YS"=YF"

12、FS"=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto 2三、使用說(shuō)明1、規(guī)定        (1) 以道路中線的前進(jìn)方向（即里程增大的方向）區(qū)分左右；當(dāng)線元往左偏時(shí)，Q=-1；當(dāng)線元往右偏時(shí)，Q=1；當(dāng)線元為直線時(shí)，Q=0。        (2) 當(dāng)所求點(diǎn)位于中線時(shí)，Z=0；當(dāng)位于中線左鍘時(shí)，

13、Z取負(fù)值；當(dāng)位于中線中線右側(cè)時(shí)，Z取正值。        (3) 當(dāng)線元為直線時(shí)，其起點(diǎn)、止點(diǎn)的曲率半徑為無(wú)窮大，以10的45次代替。        (4) 當(dāng)線元為圓曲線時(shí)，無(wú)論其起點(diǎn)、止點(diǎn)與什么線元相接，其曲率半徑均等于圓弧的半徑。         (5) 當(dāng)線元為完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲

14、率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。         (6) 當(dāng)線元為非完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。2、輸入與顯示說(shuō)明       輸入部分：      

15、1. SZ => XY      2. XY = > SZ       N ?  選擇計(jì)算方式，輸入1表示進(jìn)行由里程、邊距計(jì)算坐標(biāo) ；輸入2表示由坐標(biāo)反算              里程和邊距。     X0 ？線元起點(diǎn)的X坐標(biāo)    &

16、#160; Y0 ？線元起點(diǎn)的Y坐標(biāo)      S0 ？線元起點(diǎn)里程     F0 ？線元起點(diǎn)切線方位角     LS ？線元長(zhǎng)度     R0 ？線元起點(diǎn)曲率半徑     RN ？線元止點(diǎn)曲率半徑     Q ？  線 元左右偏標(biāo)志(左偏Q=-1，右偏Q=1，直線段Q=0）     &#

17、160;S ？ 正算時(shí)所求點(diǎn)的里程      Z ？正算時(shí)所求點(diǎn)距中線的邊距(左側(cè)取負(fù)，值右側(cè)取正值，在中線上取零)      X ？反算時(shí)所求點(diǎn)的X坐標(biāo)      Y ？反算時(shí)所求點(diǎn)的Y坐標(biāo)       顯示部分：       XS=×××  &#

18、160; 正算時(shí)，計(jì)算得出的所求點(diǎn)的X坐標(biāo)       YS=×××    正算時(shí)，計(jì)算得出的所求點(diǎn)的Y坐標(biāo)       FS=×××    正算時(shí)，所求點(diǎn)對(duì)應(yīng)的中線點(diǎn)的切線方位角       S=×××   

19、 反算時(shí)，計(jì)算得出的所求點(diǎn)的里程       Z=×××    反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距四、算例        某匝道的由五段線元（直線+完整緩和曲線+圓曲線+非完整緩和曲線+直線）組成，各段線元的要素（起點(diǎn)里程S0、起點(diǎn)坐標(biāo)X0 Y0、起點(diǎn)切線方位角F0、線元長(zhǎng)度LS、起點(diǎn)曲率半徑R0、止點(diǎn)曲率半徑RN、線 元左右偏標(biāo)志Q）如下：S0   X0

20、60;      Y0                   F0            LS R0 RN     Q500.000    19942.837 &

21、#160;  28343.561     125 16 31.00    269.256     1E45 1E45      0769.256    19787.340    28563.378     125 16 31.00    &

22、#160; 37.492     1E45 221.75     -1806.748    19766.566    28594.574     120 25 54.07     112.779    221.75 221.75     -1919.527 &

23、#160;  19736.072    28701.893       91 17 30.63       80.285    221.75 9579.228   -1999.812    19744.038    28781.659    

24、;   80 40 50.00     100.000    1E45 1E45       0         1、正算          （注意：略去計(jì)算方式及線元要素輸入，請(qǐng)自行根據(jù)所求點(diǎn)所在的線元輸入線元要素）   

25、60;        S=700  Z=-5    計(jì)算得  XS=19831.41785   YS=28509.72590            S=700  Z=0     計(jì)算得  XS=19827.33592&#

26、160;  YS=28506.83837            S=700  Z= 5    計(jì)算得  XS=19823.25398   YS=28503.95084            S=780  Z=-5 

27、60;  計(jì)算得  XS=19785.25749   YS=28575.02270            S=780  Z=0     計(jì)算得  XS=19781.15561   YS=28572.16358         &#

28、160;  S=780  Z= 5    計(jì)算得  XS=19777.05373   YS=28569.30446                       S=870  Z=-5    計(jì)算得&#

29、160; XS=19747.53609   YS=28654.13091            S=870  Z=0     計(jì)算得  XS=19742.68648   YS=28652.91379            S=87

30、0  Z= 5    計(jì)算得  XS=19737.83688   YS=28651.69668                   S=940  Z=-5.123    計(jì)算得  XS=19741. 59118  YS=

31、28722.05802            S=940  Z=0            計(jì)算得  XS=19736.47687   YS=28722.35642           

32、60;S=940  Z= 3.009    計(jì)算得  XS=19733.47298   YS=28722.53168                   2、 反算             X=198

33、31.418   Y=28509.726   計(jì)算得  S=699.9999974  Z= -5 .00018164             X=19827.336   Y=28506.838   計(jì)算得  S=699.9996493  Z= 0.000145136    

34、60;        X=19823.25398   Y=28503.95084    計(jì)算得  S=699.9999985  Z= 5.000003137           X=19785.25749   Y=28575.02270  計(jì)算得  S=780.0000035  Z= -5 .000001663           X=19781.15561   Y=28572.16358   計(jì)算得  S=780.0000025  Z=- 0.000002979