CASIO fx5800P CYT 综合曲线坐标正反算程序V3.1.doc

CASIO 卡西欧计算器程序制作：沈阳毅适用：CASIO fx-5800PQ Q：4 4 3 9 8 3 2 7 3 原创 CASIO fx-5800P CYT 综合曲线坐标正反算程序V3.0说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。包含：交点前直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、交点后直线段的线路中、边桩坐标及切线方位角,该程序可计算不对称缓和曲线设计并能直接计算出极坐标放样数据,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。主程序名：CYT第1步 50DimZ：Deg：Fix 6第2步 0=>XS：1=>XY：-1=>DK：?Q：DK?Z第3步 QJ：If J=0：Then Prog CYTXS：Else If J <0：Then Prog CYTFS：Else If J>0：Then Lbl 1：W?P：?Q： Prog CYTSJ：Prog CYTYS：Prog CYTZS：IfEnd第4步 Cls：F=：Locate 3,1,Z27°：X=：Locate 3,2,X：Y=：Locate 3,3,Y第5步 Z(X)?U：Z(Y)?V第6步 Pol(X-U,Y-V)：If J<0：Then J+360J：IfEnd：Prog CYTSQX第7步 Cls：F=：Locate 3,1,J°：S=：Locate 3,2,I：H=：Locate 3,3,X第8步 DK?Z：Goto 1曲线要素子程序名：CYTYS第1步 S(2)÷24÷R-S(4)÷2688÷R(3) + S(6)÷506880÷R(5)Z5：T(2)÷24÷R-T(4)÷2688÷R(3)+ T(6)÷506880÷R(5)Z6第2步 S÷2-S(3)÷240÷R(2)Z7：T÷2-T(3)÷240÷R(2)Z8 第3步 90×S÷÷RZ9：90×T÷÷RZ10 第4步 (Abs(A)-Z9-Z10)×R×÷180Z2 第5步 (R+Z5)×tan(Abs(A)÷2)+(Z6-Z5)÷sin(Abs(A)+Z7Z11 第6步 (R+Z6)×tan(Abs(A)÷2)+(Z5-Z6)÷sin(Abs(A)+Z8Z22 第7步 (R+Z5)÷cos(Abs(A)÷2)-RZ1 第8步 Z2+S+TL：D-Z11Z12：Z12+SZ13：Z13+Z2Z15：Z12+LZ16 第9步 F+180Z17：If Z17>360：Then Z17-360Z17：IfEnd：N-Z11×cos(F)Z18：E-Z11×sin(F)Z19 第10步 F+AZ20：If Z20>360：Then Z20-360Z20：IfEnd：N+Z22×cos(Z20)Z21：E+Z22×sin(Z20)Z23显示曲线要素子程序名：CYTXS第1步 Prog CYTSJ：Prog CYTYS第2步 Cls：T1=：Locate 4,1,Z11：T2=：Locate 4,2,Z22：E=：Locate 4,3,Z1：L=：Locate 4,4,L第3步 Cls：ZH=：Locate 4,1,Z12：HY=：Locate 4,2,Z13：YH=：Locate 4,3,Z15：HZ=：Locate 4,4,Z16第4步 Prog CYT正算子程序名：CYTZS第1步 If ZZ12：Then Goto I：Else If ZZ13：Then Goto E：Else If ZZ15：Then Goto F：IfEnd第2步 If ZZ16：Then Goto G：Else If ZZ16：Then Goto H：IfEnd第3步 Lb1 I：Z12-ZG：Z18+G×cos(Z17)+P×cos(F+Q)X：Z19+G×sin(Z17)+P×sin(F+Q)Y：FZ27：Goto 3第4步 Lb1 E：Z-Z12W：If A<0：Then -1H：Else 1H：IfEnd：(Z-Z12)(2)÷2÷R÷S×180÷Z26第5步 F+Z26×HZ27：W-W(5)÷40÷R(2)÷S(2)Z28：W(3)÷6÷R÷S-W(7)÷336÷R(3)÷s(3)Z29第6步 Z18+Z28×cos(F)-H×Z29×sin(F)+P×cos(Z27+Q)X：Z19+Z28×sin(F)+H×Z29×cos(F)+P×sin(Z27+Q)Y：Goto 3第7步 Lb1 F：Z-Z12K：If A<0：Then -1H：Else 1H：IfEnd：(Z-Z12)-0.5×S)÷R×180÷Z32第8步 F+Z32×HZ33：R×sin(Z32)+Z7Z34：R×(1-cos(Z32)+Z5Z35第9步 Z18+Z34×cos(F)-H×Z35×sin(F)+P×cos(Z33+Q)X：Z19+Z34×sin(F)+H×Z35×cos(F)+P×sin(Z33+Q)Y：Z33Z27：Goto 3第10步 Lb1 G：Z16-ZM：If A>0：Then -1H：Else 1H：IfEnd：Z20+180Z42：If Z42>360：Then Z42-360Z42：IfEnd第11步 M(2)÷2÷R÷T×180÷Z38：Z20+Z38×HZ39：M-M(5)÷40÷R(2)÷T(2)Z40：M(3)÷6÷R÷T-T(4)÷336÷R(3)Z41第12步 Z21+Z40×cos(Z42)-H×Z41×sin(Z42)+P×cos(Z39+Q)X： Z23+Z40×sin(Z42)+H×Z41×cos(Z42)+P×sin(Z39+Q)Y：Z39Z27：Goto 3第13步 Lb1 H：If ZZ16：Then Z-Z16Z43：IfEnd第14步 Z21+Z43×cos(Z20)+P×cos(Z20+Q)X：Z23+Z43×sin(Z20)+P×sin(Z20+Q)Y：Z20Z27第15步 Lbl 3：Return反算子程序名：CYTFS第1步 Lbl 3：X?U：Y?V：0P：0Q：Prog CYTSJ：Prog CYTYS：Lbl 2：Prog CYTZS第2步 Z27-90Z14：(V-Y)×cos(Z14)-(U-X)×sin(Z14)Z25第3步 If Abs(Z25)>0.001：Then Z+Z25Z：Goto 2：Else Goto 1：IfEnd第4步 Lbl 1：(Y-V)÷sin(Z14)Z30第5步 Cls：X=：Locate 4,1,U：Y=：Locate 4,2,V：K=：Locate 4,3,Z：S=：Locate 4,4,Z30第6步 DK?Z：Goto 3平曲线数据库名：CYTSJ第1步 Goto 1第2步 Lbl 1：If Z<1500.000 Or Z>15000：Then Cls：Locate 2,2,PQX：Locate 4,3,CHAOXIAN：Locate 10,4,Stop第3步 Stop：IfEnd第4步 Lb1 1：If Z<5031.443：Then 3379475.002N：450979.99E：4758.99274754061D：8000R：170S：170T：2°4108.03A：113°5320.39F：Return：IfEnd第5步 Lb1 1：If Z<7152.556：Then 3378819.247N：452290.956E：6224.74386459545D：5000R：330S：330T：-17°3732.80A：116°3428.42F：Return：IfEnd第6步 Lb1 1：If Z<15000：Then 3378226.731N：456053.721E：10021.3590194243D：4500R：360S：360T：53°1246.06A：98°5655.62F：Return：IfEnd第n步 注意：数据库中第2步止点里程15000和最后一步中的止点里程15000一致，否则计算最后一段直线会出现错误。数据输入说明：第1步 Goto 1第2步 Lbl 1：If Z<本条线路起点里程Z>本条线路止点里程 Cls：Locate 2,2,PQX：Locate 4,3,CHAOXIAN：Locate 10,4,Stop第3步 Stop：IfEnd第4步 Lbl 1：If Z<本交点HZ里程或下个交点ZH里程：Then本交点坐标XN：本交点坐标YE：交点里程D：曲线半径R：第一缓和曲线长S：第二缓和曲线长T：转角A：第一直线方位角F：Return：IfEnd第n步 Lbl 1：每增加一行则为增加一个交点曲线要素。竖曲线子程序名：CYTSQX第1步 ZC：Prog CYTSQXSJ：AZ1：BZ2：CZ3：If (Z1-Z2)<0：Then 1Z8：Else -1Z8：IfEnd第2步 R×Abs(Z1÷100-Z2÷100)L：L÷2T：T(2)÷2÷RE：K-TZ4：Z4+LZ5第3步 If Z3Z4：Then H-(K-Z3)×Z1÷100X：Goto 3：Else If Z3Z4 And Z3<K：Then H+(Z3-K)×Z1÷100Z6：(Z3-Z4)(2)÷2÷RZ7：Z6+Z7×Z8X：Goto 3：Else If Z3=K：Then H+E×Z8X：Goto 3：Else If Z3>K And Z3Z5：Then H+(Z3-K)×Z2÷100Z6：(Z5-Z3)(2)÷2÷RZ7：Z6+Z7×Z8X：Goto 3：Else If Z3Z5：Then H+(Z3-K)×Z2÷100X：IfEnd：Lbl 3：Return竖曲线数据库：CYTSQXSJ第1步 Goto 1第2步 Lbl 1：If C<5000.000 Or C>12000：Then Cls：Locate 2,2,SQX：Locate 4,3,CHAOXIAN：Locate 10,4,Stop第3步 Stop：IfEnd第4步 Lb1 1：If C<5900：Then 5800K：323.527H：20000R：-0.5800A：0.4200B：Return：IfEnd第5步 Lb1 1：If C<8982：Then 8910K：336.589H：20000R：0.4200A：-0.3000B：Return：IfEnd第6步 Lb1 1：If C<12000：Then 11187K：329.758H：20000R：-0.3000A：-0.5500B：Return：IfEnd第n步 注意：数据库中第2步止点里程12000和最后一步中的止点里程12000一致，否则计算最后一段直坡会出现错误。数据输入说明：第1步 Goto 1第2步 Lbl 1：If C<本条线路竖曲线起点里程C>本条线路竖曲线止点里程n Cls：Locate 2,2,SQX：Locate 4,3,CHAOXIAN：Locate 10,4,Stop第3步 Stop：IfEnd第4步 Lbl 1：If C<本条竖曲线圆直桩号或止点桩号：Then 变坡点桩号K：变坡点高程H：竖曲线半径R：第一纵坡（百分数单位）A：第二纵坡（百分数单位）B：Return：IfEnd第n步 Lbl 1：每增加一行则为增加一个竖曲线要素。程序演示步骤：坐标正算运行主程序CYT0=>XS：1=>XY：-1=>DKQ? 1 (输入0为显示曲线要素、输入1为坐标正算、输入-1为坐标反算)DK?= 10000(计算里程)W?= 5(边桩偏距：左为负，右为正)?= 90(边桩偏角：正交为90°，斜交自拟)FWJ= 127°2302.26(计算切线方位角)X= 3377702.695(计算结果坐标X)Y= 455855.4893(计算结果坐标Y)Z(X)?= 33777000.567(测站坐标X)Z(Y)?= 455950.689(测站坐标Y)F= 271°1649.05(测站至计算点方位角)S= 95.223518(测站至计算点距离)H= 333.319（计算竖曲线高程）说明：按EXE键输入里程、偏距、偏角可进行连续正算。坐标反算运行主程序CYT0=>XS：1=>XY：-1=>DKQ? -1 (输入0为显示曲线要素、输入1为坐标正算、输入-1为坐标反算)DK? 10000（输入近似里程）X? 3377700.567 （输入反算点坐标X）Y? 455950.689 （输入反算点坐标Y）X= 3377700.567（回显计算结果X）Y= 455950.689 （回显计算结果Y）K= 10076.0649 （反算结果里程）S=-51.760382 (反算结果偏距：负为偏左,正为偏右)说明：按EXE键输入近似里程和坐标X、Y可进行连续反算。显示曲线要素运行主程序CYT0=>XS：1=>XY：-1=>DKQ? 0 (输入0为显示曲线要素、输入1为坐标正算、输入-1为坐标反算)DK? 5000（输入里程调用曲线要素）T1= 272.525270(第一切线长)T2= 272.525270(第二切线长)E= 2.348042(外失距)L= 544.975457(曲线长)ZH= 4486.467430(直缓桩号)HY= 4656.467430(缓圆桩号)YH= 4861.442887(圆缓桩号)HZ= 5031.442887(缓直桩号)五大桩对照表ZH=04+486.467,X=3379585.366,Y=450730.8111,FWJ=113°5320.39HY=04+656.467,X=3379515.972,Y=450886.0019,FWJ=114°2951.94QZ=04+758.955,X=3379472.878,Y=450978.9889,FWJ=115°1354.40HY=04+861.443,X=3379428.597,Y=451071.4162,FWJ=115°5756.86HZ=05+031.443,X=3379353.085,Y=451223.7237,FWJ=116°3428.42程序输入说明： 输入程序时每输入完成1步则按EXE换行，需认真、仔细，子程序输入方式和主程序相同，若字体小请按住CTRL键然后再滑动鼠标滚动条进行缩放页面大小。本程序中多数指令在FUNCTION中（2：COMPLX、3：PROG）查找。Deg和Fix 按SHIFT再按MODE SETUP查找。程序数据库可采用“曲线坐标计算程序VBA”自动生成fx-5800P数据库。其他说明须知：当计算里程小于或大于数据库中指定的桩号时会提示超限显示。如：输入里程1000会提示如下错误：提示1： PQXCHAOXIAN Stop该错误表示输入的里程超出了平曲线数据库中指定的最小或最大里程。提示2：SQXCHAOXIAN Stop该错误表示输入的里程超出了竖曲线数据库中指定的最小或最大里程。更新记录：2010-10-01：程序首次发布,共125步骤,释义详细。2010-12-01：V1.0 调整了首次步骤多的问题,精简到了50步,修改了主程序显示内容及相关布置。2010-12-06：V2.0增加了坐标反算功能。2012-12-17：V3.0增加了竖曲线高程计算，修正了程序运行方式。2013-06-30：V3.1修正了缓和曲线坐标计算遇小半径出现误差现象（修改程序字体加粗部分）。 如有雷同纯属巧合