fx5800任意平曲线极坐标与直角坐标放样,求里程及边距.doc

CASIO-fx5800程序源(任意平曲线极坐标与直角坐标放样,求里程及边距) CASIO-fx5800程序源(任意平曲线极坐标与直角坐标放样,求里程及边距)CASIO-fx5800(任意平曲线极坐标与直角坐标放样,求里程及边距)湖南津市 赵济汉主程序PQX Deg : Fix 3 : “XZ0:YZ1”？A : IfA = 1: Then Goto 1 : IfEnd IfA = 0 : Then “ BS0:XY1:AND2:DK3:L(I)4 ” ？O : IfEnd IfO = 4: ThenGoto 1 : IfEnd IfO = 3: ThenProg“ F.2 ” : IfX= 0 : ThenGoto 1 : IfEnd : IfEnd IfO 1: Then “ X1 ” ？X : “ Y1” ？Y : XZ11: YZ12: “ X2 ” ？P : “ Y2” ？Q : Pol( PX , QY) : IfJ0 : ThenJ + 360J : IfEnd : Cls : “ S12= ” :Locate 6 ,1, I : “ B12= ” : J ?DMS1B : IfEnd IfO = 1: Then “XY(0) 0: 1” ? B: IfEnd Lbl1 : IfA = 1AndO = 3 :ThenProg“ F.2 ” : IfEnd : Prog “ Q.1 ”子程序PPQX （给定综合曲线属性）“ ZQX 0: * 1: * 2 * ” ？C IfC=0: Then 100Z98: n0 Z99: IfEnd IfC=1: Thenm1Z98: n1 Z99: IfEnd IfC=2: Thenm2Z98: n2 Z99: IfEnd *Prog “ PQX ” 其中：mi为第i条综合曲线矩阵变量起始单元，(第1条综合曲线矩阵变量起始单元应自100开始，以便其他程序临时使用100以前的矩阵变量，综合曲线中每一基本单元需8个矩阵存储单元)； ni为第i条综合曲线中基本单元数，每增加一个条综合曲线则增加一行IfC=*Z99: IfEnd语句。还应注特别意，原有程序中的矩阵变量定维语句DimZ应取消，只要原有程序中的矩阵变量定维不大于100，原有程序的运行不受影响。子程序Q . 5(建立数据库)Deg : Fix 5 0I : “ N ” ? W : W Z99: “ZMIN”? V : VZ98Lbl1 : “NO.I= ”: I LblC : V + 8 IG “DK ”？A : AZG“X ”？B : BZG + 1“Y ”？C : CZG +2“AT ”？D : DZG +3“(1Y+1) ”？K : KZG +4“L ”？E : EZG +5“R1 ”？F : FZG +6“R2”？O : OZG +7Lbl2 : “JX1: NO0 : OK2”? M : IfM = 2 : ThenGoto 3 : IfEnd : IfM = 1: ThenGotoC : IfEnd I + 1I : Goto1 Lbl3 : “END ”子程序Q . 6Lbl1 : “S ” ? W : IfW = 0 : ThenGoto 2 : IfEnd IWW : Cls : “S= ” : Locate 5 , 1 , W : “ BP= ” : J ?DMSGoto1 Lbl2 : Cls : “ OK ”子程序Q.1IfX = 0AndO = 3 : ThenGoto2 : IfEnd Lbl1 : Z98 I : Z99 J : “ DKP ” ? S IfS = 1: ThenProg“JH . ZJ ” : Goto1 : IfEnd （加桥台转轴时）IfS = 2: ThenProg“ZH . 2 ” : Goto1 : IfEnd （加涵通转轴时）IfS = 0: ThenGoto2 : IfEnd : IfSZ I : Then “ DKP Z I + 8 J : Then “ DKPMAX ” Goto1 : IfEnd IfO 4AndB 0 : Then “ B ” ? V : V : Prog “ J ” : TU : “()D2 ”? N : Else0 N : IfEnd Prog “ Q.2 ” : Goto1 Lbl2 : “ END ”子程序Q . 2IfO 3: Then “ ” : IfEnd 0I : Z98+ 8 M IfS ZM8 AndS ZM :ThenGoto1 : IfEnd Lbl2 : I +1 I : M + 8 M IfS ZM :ThenGoto1 : ElseGoto2 : IfEnd Lbl1 : IfO = 3: ThenI B :IfEnd : MZ27:IfO 4: ThenProg “Q. 3” : Else“ L(I) , I = ” : I Cls : “ MIN = ” : Locate6 , 1 , ZZ98+8 I : “ MAX = ” : Locate6 ,2 , ZZ98+8 I +8:IB : IfEnd 子程序Q . 3IfO = 3AndU = 0 : ThenBI : IfEnd Z98+8 I G ZGZ ZG +1Z1: ZG +2Z2: ZG +3Z3 Z3:Prog“ J ” : TZ3ZG +4Z4: ZG +5M : 1 ZG +6Z5: 1 ZG +7Z6Z6Z5Z6: SZD IfO = 3: Then0N : IfEnd : Prog“Q . 4 ”子程序Q .4(正算)Lbl0 : 0I : 0E:0F : 0G : 0H Lbl1 : I + 1I : Z3+ Z4( ( I D12) ( Z5+ I D Z624M ) ( 180 ) W : IfI 12: ThenIfFrac ( I2 ) 0 : ThenE + cos ( W ) E : F + sin ( W ) F : ElseG + cos ( W ) G : H + sin ( W ) H : IfEnd : IfEnd IfI 12 : ThenGoto 1：IfEnd Z1+（D36）( cos ( Z3) + cos ( W ) + 2G + 4E ) + N cos ( U + W ) P Z2+（D36）( sin ( Z3) + sin ( W ) + 2H + 4F ) + N sin ( U + W ) Q IfO = 3 : ThenGoto2 : IfEnd IfO 0 : ThenCls : “ XP = ” : Locate 5 , 1 , P : “ YP = ” : Locate 5 , 2 , Q : “ AT = ” : W?DMSIfEnd Z27I : IfZI S = 0AndZI + 1= 0 AndN = 0 : ThenP ZI + 1: QZI + 2: W : Prog “ H ” : TZI +3: IfEnd IfO 1: ThenPol( PZ11, QZ12) : IfJ ZM +8: ThenZM +8S ; IfEndIfU = 0 : ThenProg“Q . 3”: ElseProg“Q . 2”: IfEnd “? ” Lbl1 : ( YQ ) cos （W90）(XP) sin (W90) Z8Z + Z7+ Z8SIfAbs ( Z8) 0.0001: Then Z7+ Z8D : Prog“Q .4”: DZ7: Goto1: IfEnd IfU =1: ThenIfS + 0.01 Z : Then “ DKP ” : S IfB Z991 : Then B1B : ZZ98+ 8B I : ZJ : Else ZI : ZZ98+ 8 Z99J : IfEnd : “ L(I)= ” : B Cls : “ MIN= ” : Locate 6 , 1 , I : “ MAX= ” : Locate 6 , 2 , J : “ END, ?ZH.D(LI)”X Z13: 0X : 3O: Goto3 : IfEnd: IfEnd IfU =0: ThenIfS Z 0.01 : Then “ DK ZZ98+ 8(B+1) + 0.01: Then “ DKMAX , L(I+1) ” Goto2 : IfEnd : IfEnd Cls : “ DKP=” : Locate6 ,1, S : XPP : IfP =0: Then109P : IfEnd : Pol( P , YQ) : Ifsin（WJ）0: Then II : IfEnd : “ D2= ”: Locate 6 , 2, I : Goto2Lbl3子程序 JAns . 36 16 ( Int( Ans ) ) 9Int (100 (Ans Int (Ans ) ) 90T子程序 H. 6 4 Int( Ans )+. 36 Ans + . 004 Int ( ( 60 ( Ans ( Int (Ans ) ) T 为了实现桥台转轴坐标计算，现（2009.04.19）加一个子程序如下：子程序 JH . ZJ (桥台转轴)“?B ”? V : V : Prog“J : W + TT Lbl1 : “?(Y +) D ”? N : N = 0 ? Goto 2 : “?( +) S ”? SPol(N , S ) : Rec( I ,TJ ) Cls : “X = ”: Locate6, 1 , I + P : “Y = ”: Locate6 , 2 , J + Q Goto 1 Lbl2 : “ OK ” 在子程序Q.1 中语句 Lbl1 : Z98 I : Z99 J : “ DKP ” ? S 后加IfS = 1: ThenProg“JH . ZJ ” : Goto1 : IfEnd （加桥台转轴时）就可以了，(2009.04.19以加，见绿色标注语句。桥台转轴直角坐标计算 运行提示符 直角坐标放样令DKP=-1时，转向桥台转轴直角坐标计算提示符说明 以上一测点中桩为中心作偏转? B ? 问转轴偏角（以小数点为度）?(Y + ) D?问偏距（与直角坐标放样时意义一致）。赋值0时转向OK?( + ) S? 问支距（小里程方向为，大里程方向为+）X= 给出点位纵坐标Y= 给出点位横坐标OK本转轴各点计算结束，返回直角坐标放样。子程序 ZH . 2 (涵通转轴)“?B ”? V : V : Prog“J : W + TTLbl1 : “?( +) ? K ”? S : S = 0 ? Goto 2 : “?(Y +) D ”? N N Cos ( T) + P X : N Sin ( T) + Q YCls : “X = ”: Locate6, 1 ,S Cos( W ) + X : “Y = ”: Locate6 , 2 , S Sin(W ) + Y Goto 1 Lbl2 : “ OK ” 在子程序Q.1 中语句IfS = 1: ThenProg“JH . ZJ ” : Goto1 : IfEnd 后加 IfS = 2: ThenProg“ZH . 2 ” : Goto1 : IfEnd （加涵通转轴时)即可。 涵通转轴直角坐标计算中桩直角坐标放样令DKP=2时，提示符示例说明 以上一测点中桩为中心作偏转? B ? 问涵通偏角（以小数点为度）?( + ) K?问里程差（相对于涵通中桩，小里程为、大里程为+），赋值0时转向OK?(Y + ) D?问偏距（与直角坐标放样时意义一致）X= 给出点位纵坐标Y= 给出点位横坐标OK 本转轴各点计算结束，返回直角坐标放样。算例1 例某工程C匝道: 根据算例数据确定综合曲线属性: 令名 C.ZD 令编号0 令（矩阵存起始储单元 ）ZMIN= 100 （综合曲线段数,不含终点）N = 4红色数据应在Q.5程序运行时按提示赋值一 给矩阵存储单元定维:（COMP状态下） 140DimZ 二 编写程序PPQX（PROG状态下给定综合曲线属性）“ C. ZD 0 ” ？ C IfC=0: Then 起始单元100Z98: 段数4 Z99: IfEnd Prog “ PQX ” 三启动Q . 5按下列数据建立矩阵存储表（COMP状态下按提示赋值，详见附件运行提示符） 自然段名 称L0L1 L2 L3L4 DK0 190 355.927472.168561.791X20934.49521066.11921142.859 21076.99320988.267Y89274.17289411.18289552.42789640.28689644.485AT 47.08177 45.09322 93.33281160.09323 185.500271 -11 11 0L190 165.927116.24189.623 0R155005500100 100 0R25500100 100 10 20 0 备注： 起始矩阵存储单元 Z100 终了矩阵存储单元 Z139 。 (139 =100 +8 5-1)四 启动程序PPQX（COMP状态下指定需计算的综合曲线编号） C. ZD 0 ？提示C匝道赋值 0回车（自动转向PQX进入各项计算）为了提高计算速度，现提供Q.4 (5点法加速),以该程序取代原有Q.4,其它无需改变。取代后按指定区间反算时间大约4秒，有兴趣的朋友可以试试。Q.4 (5点法加速)IfZ4= 0 : Then 1Z4: IfEndZ4ZG+6Z9: Z4ZG+7Z10( Z9Z10) （ 2MZ9Z10）Z13: SZ Z14: Z13Z14Z13: 1 Z9Z220.1739274226 Z15: 0.3260725774 Z16: 0.0694318442 Z17: 0.3300094782 Z181Z18Z19: 1Z17Z20: 180 Z21Z3+ Z17Z14(Z22+ Z17Z13) Z21Z17Z3+ Z18Z14(Z22+ Z18Z13) Z21Z18Z3+ Z19Z14(Z22+ Z19Z13) Z21Z19Z3+ Z20Z14(Z22+ Z20Z13) Z21Z20Z3+ Z14 ( Z22+ Z13) Z21WZ1+ Z14 ( Z15cos( Z17) + Z16cos( Z18) + Z16cos( Z19) + Z15cos( Z20) ) +N cos( U + W ) PZ2+ Z14 ( Z15sin( Z17) + Z16sin ( Z18) + Z16sin ( Z19) + Z15sin ( Z20) ) + N sin (U + W) QIfO = 3 : ThenGoto 2 : IfEndIfO 0 : ThenCls : “ XP= ” : Locate 5 , 1 , P : “ YP= ” : Locate 5 , 2 , Q : “ AT= ” : W?DMSIfEndZ27I : IfZI S = 0AndZI+1= 0And N = 0 : ThenPZI+1: Q ZI+2: W : Prog “ H ” : T ZI+3: IfEndIfO 1: ThenPol( PZ11, QZ12) : IfJ MAX提示测点里程超出本综合段终点里程DKP中桩te=mailto:雨林木风 dateTime=2009-03-30T20:07赋值0，边桩赋值1DKP?te=mailto:雨林木风 dateTime=2009-03-30T20:09问测点里程（赋值0时转向程序运行终止）B?te=mailto:雨林木风 dateTime=2009-03-30T20:11问测点偏角（以小数点为度）te=mailto:雨林木风 dateTime=2009-03-30T20:12，中te=mailto:雨林木风 dateTime=2009-03-30T20:08桩不问(+)D2?te=mailto:雨林木风 dateTime=2009-03-30T20:12问测点偏距，te=mailto:雨林木风 dateTime=2009-03-30T20:13中te=mailto:雨林木风 dateTime=2009-03-30T20:09桩不问XP=YP=AT=te=mailto:雨林木风 dateTime=2009-03-30T19:58给出测点te=mailto:雨林木风 dateTime=2009-03-30T19:59纵te=mailto:雨林木风 dateTime=2009-03-30T20:00坐标te=mailto:雨林木风 dateTime=2009-03-30T19:59给出测点te=mailto:雨林木风 dateTime=2009-03-30T20:00横坐标te=mailto:雨林木风 dateTime=2009-03-30T19:59给出测点te=mailto:雨林木风 dateTime=2009-03-30T20:01切线te=mailto:雨林木风 dateTime=2009-03-30T20:00方向te=mailto:雨林木风 dateTime=2009-03-30T20:02ENDte=mailto:雨林木风 dateTime=2009-03-30T20:02序运行程终止 (DKP=0时)te=mailto:雨林木风 dateTime=2009-03-30T20:17两种坐标都te=mailto:雨林木风 dateTime=2009-03-30T20:18需要te=mailto:雨林木风 dateTime=2009-03-30T20:14问计算状态，赋值2时te=mailto:雨林木风 dateTime=2009-03-30T20:17te=mailto:雨林木风 dateTime=2009-03-30T20:17DKP?同上te=mailto:雨林木风 dateTime=2009-03-30T20:18*从略te=mailto:雨林木风 dateTime=2009-03-30T20:18*te=mailto:雨林木风 dateTime=2009-03-30T20:21同时te=mailto:雨林木风 dateTime=2009-03-30T20:19给出以上te=mailto:雨林木风 dateTime=2009-03-30T20:20两种结果te=mailto:雨林木风 dateTime=2009-03-30T20:21ENDte=mailto:雨林木风 dateTime=2009-03-30T20:02程序