CASIOfx4850万能坐标计算程序[Machao2006.10.23].doc

CASIO fx4850万能坐标计算程序[Machao 2006.10.23]
一、程序内容
程序：Mc
 W“1.js   2.SZ”:W=1=>Goto 1 ΔW=2=> O“KOU  LING”:O≠123456=>O=0“OUT” ◢Goto 4 ΔO=0:
Defm 83:Z[81]=1:Goto 0←┘
Lbi 0←┘
{ABCDEFGU}:A“X0”:B“Y0”:C“F0’:D“R0”:E“RN’:F“D0”:G“LS”:U“G” ←┘
Z[Z[81]]=A:Z[Z[81]+10]=B:Z[Z[81]+20]=C:Z[Z[81]+30]=1÷D:Z[Z[81]+40]=1÷E:Z[Z[81]+50]=F: Z[Z[81]+60]=F+G: Z[Z[81]+70]=U←┘
A=0=> Dsz Z[81]:Goto 1ΔIsz Z[81]:Goto 0←┘
Lbi 1←┘
{HS}:H“D’:S“Z”:Z[82]=1:Goto 2←┘
Lbi 2←┘
Z[81]>Z[82] =>”OUT” ◢:Goto5←┘
H≤Z[Z[82]+60]=>A=Z[Z[82]]:B=Z[Z[82]+10]:Z[Z[82]+20]=C:D=Z[Z[82]+30]: E=Z[Z[82]+40]: F=Z[Z[81]+50]: G=Z[Z[81]+60]: U=Z[Z[81]+70]: Goto3:ΔIsz Z[82]:Goto 2
Lbi 3←┘
P＝U（E-D）÷Abs（G-F）：Q＝Abs（H-F）：I＝PQ：J＝C+90 Q（I+2UD）/π： J<0=>J=J+360ΔM＝C+45 Q（I÷4+2UD）÷2π：N＝C+135Q（3I÷4+2UD）÷2π：K＝C+45Q（I÷2+ 2UD）÷π                                             ←┘
S=0=>Z=0：Goto4ΔS＜0=>Z＝－1：Goto4ΔZ＝1：Goto4←┘
Lbi 4
X＝A+Q(Cos C+4（Cos M+Cos N）+2Cos K+Cos J) ÷12 +ZSCos（J+90Z) ←┘     
Y＝B+Q(Sin C+4（Sin M+Sin N）+2Sin K+Sin J) ÷12 +ZSSin（J+90Z)   ←┘
S=0=>“X”:X:Pause 0: “Y” :Y◢ Goto 1Δ
S<0=>“XL”:X:Pause 0: “YL”:Y◢ Goto 1Δ
S>0=>“XR”:X:Pause 0: “YR”:Y ◢ Goto 1 ←┘
Lbi5←┘

    二、 变量及说明
X0：Y0：F0――――曲线元起点X、Y坐标及起点正切线方位角
R0：RN――――――曲线元起点及终点半径
D0：LS：Q―――曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）
D――――――――曲线元中待求点桩号
Z―――――――――计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0
X：Y
XL:YL
XR:YR ―――――――待求点中边桩桩号的X，Y坐标
扩充变量：Z[Z[81]]: Z[Z[81]+10]: Z[Z[81]+20]: Z[Z[81]+30]: Z[Z[81]+40]: Z[Z[81]+50]: Z[Z[81]+60]: Z[Z[81]+70]: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:Q
        (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
        (2) 当所求点位于中线时，Z=0，坐标显示X  Y；当位于中线左侧时，Z取负值,坐标显示XL  YL,；当位于中线右侧时，Z取正值,坐标显示XR  YR。
        (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
        (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。
         (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45
次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
         (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的
值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
本程序可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行正算。
本程序是对扩充变量的应用，实现了真正意义上的的全线贯通。程序分为两部分：1为计算2为设置。首先对内存变量进行扩充，增加了83个扩充变量，为Z[1]~Z[83]，为防止误操作在输入正确口令123456时方可进行设置。本程序设计为十段线元要素，如有需要可以根据内存容量进行增加。在执行程序中将各个线元要素按规律输入到扩充内存变量中，在坐标计算时首先判断其在哪一线元内，并把其线元要素调出进行计算。在进行计算时直接输入桩号、宽度就可以进行全线坐标计算。
    本程序也可在4800上运行，只需将输出部分进行修改。
 四、 示例
郑石高速路面9标主线线元要素及相关数据如下：
X0             Y0             F0             R0          RN         D0        LS       Q
744450.244  510454.431       220º41 º31.8 º   7000        7000    116075.889   3386.293    1
742505.351  507722.676       248º24 º33.6 º   1045                 1045     119462.182    1980.825    0
741776.461  505880.832       248º24 º33.6 º   8500        8500    121443.007   1517.999    1
741346.546  504427.086      258º38 º30.1 º   1045                 1045     122961.007    1985.977    0

桩号          宽度                X              Y 
117000         0              743791.308    509807.481
118000        -20             743164.306    509027.467
119000         20             742707.550    508137.933
119462.182      0             742505.349    507722.673
120000       -15.253          742293.266     507228.206
121000       15.260           741953.665    506287.142
121443.007      0             741776.461    505880.833
122000       -10.20           741578.902    505359.701
122961.007    50.355          741395.915    504417.169
123650         0             741210.853    503751.587