CASIO万能坐标计算程序.doc

CASIO万能坐标计算程序[MC2007.01.11]? 程序：MCW“1.JS?? 2.SZ”:W=1=Z[2]=0:Goto 1 ΔW=2= O“KOU??LING”:O≠123456=O=0: “OUT” Goto 5ΔO=0:Defm 42:Z[1]=0:Goto 0←┘Lbi 0←┘{ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” ←┘Z[Z[1]*8+3]=A:Z[Z[1]*8+4]=B:Z[Z[1]*8+5]=C:Z[Z[1]*8+6]=1÷R:Z[Z[1]*8+7]=1÷E:Z[Z[1]*8+8]=F: Z[Z[1]*8+9]=F+G: Z[Z[1]*8+10]=U←┘A=0= Z[2]=0:Goto 1ΔIsz Z[1]: Goto 0←┘Lbi 1←┘{DZ}:D：Z：Z[2]=0:Goto 2←┘Lbi 2←┘D≤Z[Z[2]*8+9]=A=Z[Z[2]*8+3]:B=Z[Z[2]*8+4]: C =Z[Z[2]*8+5]:R=Z[2]*8+6]: E=Z[Z[2]*8+7]: F=Z[Z[2]*8+8]: G=Z[2]*8+9]: U=Z[Z[2]*8+10]: Goto3:ΔIsz Z[2]:Goto 2Lbi 3←┘P＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）：I＝PQ：J＝C+90 Q（I+2UR）/π： J0=J=J+360ΔM＝C+45 Q（I÷4+2UR）÷2π：N＝C+135Q（3I÷4+2UR）÷2π：K＝C+45Q（I÷2+ 2UR）÷π:L=C+45Q(I÷8+2UR) ÷4π: S=C+135Q(3I÷8+2UR) ÷4π: T=C+225Q(5I÷8+2UR) ÷4π: H=C+315Q(7I÷8+2UR) ÷4π←┘Lbi 4X＝A+Q(Cos C+4（Cos L+ Cos S +Cos T+Cos H）+2(Cos M + Cos N+Cos K)+Cos J) ÷24+ZCos（J+90)←┘???? Y＝B+Q(Sin C+4（Sin L + Sin S + Sin T+ Sin H）+2(Sin M + Sin N+ Sin K) + Sin J) ÷ 24 +ZSin（J+90)←┘Z=0=“X”:X:Pause 0: “Y” :Y◢ Goto 1ΔZ0=“XL”:X:Pause 0: “YL”:Y◢ Goto 1Δ?????????? fx4850Z0=“XR”:X:Pause 0: “YR”:Y ◢ Goto 1 ←┘Z=0= X “X” ◢Y “Y”◢ Goto 1ΔZ0= X “XL” ◢Y “YL”◢ Goto 1Δ????????????????fx4800Z0= X “XR” ◢Y “YR”◢ Goto 1 ←┘Lbi 5←┘????二、 变量及说明X0：Y0：F0――――曲线元起点X、Y坐标及起点正切线方位角R0：RN――――――曲线元起点及终点半径D0：LS：Q―――曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）D――――――――曲线元中待求点桩号Z―――――――――计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0X：YXL:YLXR:YR ―――――――待求点中边桩桩号的X，Y坐标扩充变量：Z[Z[1]*8+3]: Z[Z[1]*8+4]: Z[Z[1]*8+4]: Z[Z[1]*8+6]: Z[Z[1]*8+7]: Z[Z[1]*8+8]: Z[Z[1]*8+9]: Z[Z[1]*8+10]: 分别为各线元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) 当线元为非完整缓和曲线时