CASIO5800测量程序

1、CASIO fx-4800P型计算器用于线路施工曲线中线点坐标的计算程序中铁十局三建公司工程技术部 摘要：本文介绍了CASIO fx-4800P型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。关键词：曲线 线路 施工测量 计算程序1概述 过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。随着电子技术进步和经济发展，测量仪器和测量方法的不

2、断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。为了能够快速准确地为全站仪提供测设数据，发挥全站仪快速测设的特

3、点，提高测量工作效率，应采用可编程的计算器，编制计算程序。本文主要介绍应用CASIO fx-4800P型计算器的计算程序，供公司测量同行们参照使用。2计算程序QXZBJS（文件名：曲线坐标计算fx-4800P） Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+MG=RA/180：“S”S=G+LLbiA：C，V：C“CSDLC=”：V“HXPJ=”：D=C-K：DLI=D-D5/（40R2L2）：U=D3/(6RL)-D7/（

4、336R3L3):J=（I2+U2）：Goto1：DGO=90（2D-L）/( R)：I=RsinO+M：U=R（1-cosO）+P：J=（I2+U2）：Goto2：D=S-（C-K）： = D-D5/（40R2L2）：Z2=D3/6RL-D7/（336R3L3）：I=T+（T-Z1）cosA-Z2sinA：U=（T-Z1）sinA+Z2cosA：J=（I2+U2）：Goto3：Lbi1:Q:Q“Z=1；Y=2”：Q=1Q=F-30D2/ （RL）：H=F-90D2/ （RL）：Q=F+30D2/ （RL）：H=F+90D2/（ RL） Goto4：Lbi2:Q:Q“Z=1；Y=2”：Q=1Q

5、=F-tng-1(U/I)：H=F-O：Q= F+tng-1(U/I)：H=F+OGoto4：Lbi3:Q:Q“Z=1；Y=2”：Q=1Q=F-tng-1(U/I):H=F-(A-90(S-(C-K)2/ （RL）): Q= F+tng-1(U/I):H=F+(A-90(S-(C-K)2/ （RL）): Goto4： Lbi4:B=90+H:H0H“QXFWJ”=B+360 H360H“QXFWJ”=H-360 H“QXFWJ”=HGoto5：Lbi5: X“CSD:X”=JcosQ+N+VcosB Y“CSD:Y”=JsinQ+E+VsinBGotoA 3程序说明3.1 输入已知变量 R圆曲

6、线半径，显示R？ L缓和曲线长，显示L？ A曲线转向角，显示A？ E直缓点纵坐标，显示ZH：X=？ N直缓点横坐标，显示ZH：Y=？ F第一切线方位角，即ZH至JD的方位角，显示FWJ=？ K直缓点里程，显示ZH：LC=？ 3.2 计算待求量 T切线长度，显示T= S曲线全长，显示S= Z3外矢距，即JD到QZ的距离，显示E0= 3.3 输入待求变量 K输入待求（测设）点的里程，显示LC=？ V横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。3.4 计算结果和说明X待求（测设）点的纵坐标，显示X= Y待求（测设）点的横坐标，显示Y= 该程序既可对

7、曲线内中线任意点的坐标进行计算，又可以对中线外移桩的坐标进行计算。计算原理如下：先建立整个曲线的坐标系统，确定曲线起点（ZH）的坐标，以ZH到JD的第一切线方位角为计算方位角。 在程序计算过程中，当屏幕显示Z=1；Y=2？的提示时可根据曲线的转向，从下列两种情况中选择一种予以回答。若曲线为左转，则输入数值1；若曲线为右转，则输入数值2。 本程序式输入时，凡遇到、两种符号时，不应将其分开，应视为一个整体的符号。 由于变量的增加，程序式中使用到Z1、Z2、Z3三个计算量，此时应将计算器的内存变量扩大，另行增加36个变量即可。在计算器上的操作过程为：按SHIFT键再按Defm键，然后按数字36即可。

8、3.5 本程序在实际测量工作中的应用 本程序不仅可以计算曲线线路上任意中线点的坐标，亦可以计算线路左、右侧的法线偏距点的坐标，用于曲线中桩和边桩放样。 可以计算曲线线路中心点的切线方位角和法线方位角，该方位角可以用于桥梁墩、台（亦为曲线上桥梁墩、台截面纵、横轴线的方位角）定位。其切线方位角在计算器上显示为“QXFWJ= ”。 本程序具有多方面的计算功能，可根据仪器状况和测设对象以及工作需要而灵活应用，可在测设现场计算，减轻内业的计算时间，提高工作效率。4计算举例 结合某线路的曲线资料，计算曲线线路上中桩或边桩的坐标。4.1已知要素 R=1500， L=150， AZ=30° 主点里程

9、：ZHDK8+150.000 HYDK8+300.000 YHDK8+935.398 HZDK9+085.398 4.2 建立平面控制坐标系 建立独立坐标系统，若以ZH为坐标原点，以ZHJD的切线为X轴的正方向，则ZH点的坐标可设为：X=8150.000，Y=1000.000；切线的初始方位角F“FWJ= ”可设为0°0000。现以计算圆曲线内DK8+400的中线点或左侧15米边桩点的坐标为例。4.3 计算操作步骤 1输入曲线半径R=1500按EXE键 2输入缓和曲线长L=150按EXE键 3输入曲线转向角A=30°按EXE键 4输入ZH点X坐标ZH：X=8150.000按

10、EXE键 5. 输入ZH点Y坐标ZH：Y=1000.000按EXE键 6输入切线方位角FWJ=0°0000按EXE键 7输入ZH点里程ZH：LC=8150按EXE键 显示T=477.0849919（切线长）按EXE键；显示S=935.398164（曲线长）；按EXE键显示E0=53.56126046（外矢距）；按EXE键8输入测设点里程LC=8400；按EXE键9输入中线点或左侧边桩的距离，如计算中桩点坐标输入PJU=0，如需计算左侧15米边桩的坐标则输入PJU=-15米；按EXE键10屏幕显示Z=1；Y=2？提示，本例为左偏曲线故输入1，按EXE键,显示FXFWJ=83.31549

11、239，按SHIFT键再按°，即DK8+400中桩点的法线方位角为83°1855.77按EXE键11. 显示CSD：X=8399.5970289，按EXE键； 显示CSD：Y=989.1782961，按EXE键。若需计算DK8+400左侧15米边桩点的坐标，按EXE键，屏幕提示LC=8400；再按EXE键，提示PJU=？，输入-15米，再按EXE键，显示Z=1；Y=2？输入数字1按EXE键，显示QXFWJ=83.31549239，按EXE键；显示CSD：X=8397.850997，按EXE键；显示CSD：Y=974.2802637，按EXE键。12若需要再计算其他线路中线点

12、或边桩点的坐标，重复811步骤操作，不必再输入各已知数据，仅输入需要计算点的里程或横向偏距及Z=1；Y=2？的参数即可。ZBFS(坐标正算)A”X1”:B”Y1”:C”X2”:D”Y2”:Fixm：Pol（C-A,D-B）:I”S=” J<0J=J+360J”FWJ” 测站点到后视点的距离及方位角X1: Y1为测站点坐标X2: Y2为后视点坐标S=为测站点到后视点的距离FWJ为测站点到后视点的方位角ZBZS(坐标正算)A”X1”:S:F:D”FWJ”: C”X2”=A+SCOS(F) D”Y2”=B+SSIN(F) X1: Y1为测站点坐标X2: Y2为后视点坐标S=为测站点到后视点的距

13、离FWJ为测站点到后视点的方位角JDZB（文件名：曲线坐标计算CASIO fx-5800P）“JDX1”?A:“JDY1”?B: “JDX2”?C:“JDY2”?D:Pol(C-A,D-B):R：L：A：JF：“ZH：LC=”?K：L2/24R-L4/2688R3P:L/2-L3/240R2M：（R+P）*tng(A/2)+MT: “T=”:TRA/180G： G+LS: “S=”:SC+TCOS(F+180) N: “ZH.X=”:ND+TSIN(F+180) E: “ZH.Y=”:E LbiA：C，V：C“DK?+”：V“PJ”：C-KD：If DL:Then D-D5/（40R2L2)

14、I: D3/6RL-D7/（336R3L3）U:（I2+U2）J：Goto1：IfEnd：IF DG:Then O=90（2D-L）/ （R）： RsinO+MI：R（1-cosO）+PU：（I2+U2）J：Goto2：Else S-（C-K）D:D-D5/（40R2L2）Z1： D3/6RL-D7/（336R3L3）Z2：T+（T-Z1）cosA-Z2sinAI：（T-Z1）sinA+Z2cosAU：（I2+U2）J：Goto3：IfEndLbi1:“Z=1；Y=2”？X：If X=1 Then F-30D2/ （RL）Q：F-90D2/ （RL）H：Else F+30D2/ （RL）Q：F

15、+90D2/（ RL）H:IfEnd: Goto4：Lbi2: X“Z=1；Y=2”：If X=1 Then F-tng-1(U/I) Q： F-OH：Else F+tng-1(U/I) Q：F+OH: IfEnd:Goto4：Lbi3: X“Z=1；Y=2”：If X=1 Then F-tng-1(U/I) Q:F-(A-90(S-(C-K)2/ （RL）) H: Else F+tng-1(U/I) Q: F+(A-90(S-(C-K)2/ （RL）) H: IfEnd: Goto4： Lbi4: 90+HB: If H0 Then B+360H: IfEnd: If H360 Then H

16、-360H : Else HH: IfEnd:Goto5：Lbi5:“X=”:JcosQ+N+VcosBX “Y=”:JsinQ+E+VsinB“FWJ=”:HGotoA 3程序说明3.1 输入已知变量“JDX1” 交点1X坐标“JDY1” 交点1Y坐标 “JDX2” 交点2X坐标“JDY2” 交点2Y坐标 R圆曲线半径，显示R？ L缓和曲线长，显示L？ A曲线转向角，显示A？ E直缓点纵坐标，显示ZH：X= N直缓点横坐标，显示ZH：Y= F第一切线方位角，即ZH至JD的方位角 K直缓点里程，显示ZH：LC=？ 3.2 计算待求量 T切线长度，显示T= S曲线全长，显示S= 3.3 输入待求

17、变量 K输入待求（测设）点的里程，显示DK?+？ V横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。3.4 计算结果和说明X待求（测设）点的纵坐标，显示X= Y待求（测设）点的横坐标，显示Y= 该程序既可对曲线内中线任意点的坐标进行计算，又可以对中线外移桩的坐标进行计算。计算原理如下：先建立整个曲线的坐标系统，确定曲线起点（ZH）的坐标，以ZH到JD的第一切线方位角为计算方位角。 在程序计算过程中，当屏幕显示Z=1；Y=2？的提示时可根据曲线的转向，从下列两种情况中选择一种予以回答。若曲线为左转，则输入数值1；若曲线为右转，则输入数值2(非1也对

18、)。3.5 本程序在实际测量工作中的应用 本程序不仅可以计算曲线线路上任意中线点的坐标，亦可以计算线路左、右侧的法线偏距点的坐标，用于曲线中桩和边桩放样。 可以计算曲线线路中心点的切线方位角和法线方位角，该方位角可以用于桥梁墩、台（亦为曲线上桥梁墩、台截面纵、横轴线的方位角）定位。其切线方位角在计算器上显示为“FWJ= ”。 本程序具有多方面的计算功能，可根据仪器状况和测设对象以及工作需要而灵活应用，可在测设现场计算，减轻内业的计算时间，提高工作效率。4计算举例 结合某线路的曲线资料，计算曲线线路上中桩或边桩的坐标。4.1已知要素 “JDX1” 交点1X坐标“JDY1” 交点1Y坐标 “JDX

19、2” 交点2X坐标“JDY2” 交点2Y坐标 R圆曲线半径，显示R？ L缓和曲线长，显示L？ A曲线转向角，显示A？ E直缓点纵坐标，显示ZH：X= N直缓点横坐标，显示ZH：Y= F第一切线方位角，即ZH至JD的方位角 K直缓点里程，显示ZH：LC=？ “JDX1” =2920660.136“JDY1” =475778.751 “JDX2” =2919624.99“JDY2” =477022.978R=4500， L=340， AZ=26°2423”(左转) 主点里程：ZH.LC=DK348+692.554 4.3 计算操作步骤 1输入交点1的X坐标 “JDX1” =2920660

20、.136按EXE键 2输入交点1的Y坐标 “JDY1” =475778.751按EXE键3输入交点2的X坐标 “JDX2” =2919624.990按EXE键 4输入交点2的Y坐标 “JDY2” =477022.978按EXE键5. 输入曲线半径R=4500按EXE键 6输入缓和曲线长L=340按EXE键 7输入曲线转向角A=26°2423”按EXE键 8输入交点2的ZH点里程ZH：LC=348692.554按EXE键 显示交点2的T=1225.972（切线长）按EXE键；显示交点2的S=2413.952（曲线长）；按EXE键E直缓点纵坐标，显示交点2的ZH：X=2920409.07

21、4 N直缓点横坐标，显示交点2的ZH：Y=476080.52279输入交点2的测设点里程DK?+?=349032.554；按EXE键10输入交点2中线点或左侧边桩的距离，如计算中桩点坐标输入PJU=0，如需计算左侧15米边桩的坐标则输入PJU=-15米；按EXE键11屏幕显示Z=1；Y=2？提示，本例为左偏曲线故输入1，按EXE键12. 显示X=2920194.946，按EXE键； 显示Y=476344.5952，按EXE键。显示QXFWJ=127°3540.43”， 即DK349+032.554中桩点的切线方位角为127°3540.43”按EXE键若需计算DK349+03

22、2.554左侧15米边桩点的坐标，按EXE键，屏幕提示DK?+?= 349032.554；再按EXE键，提示PJU=？，输入-15米，再按EXE键，显示Z=1；Y=2？输入数字1按EXE键， 显示X=2920206.831，按EXE键；显示476353.746，按EXE键。显示QXFWJ=127°3540.43”，按EXE键；13若需要再计算其他线路中线点或边桩点的坐标，重复811步骤操作，不必再输入各已知数据，仅输入需要计算点的里程或横向偏距及Z=1；Y=2？的参数即可。竖曲线程序fx-4800p一、程序名：SQX二、程序:RDNCI1：FI2：T=0.5Abs(F-C)RE=T2

23、/(2R)VTU=1,AO=-1LbI1：L：LD=>Goto2L-DT=>H=M+ (L-D)F -V(T-L +D)2/(2R)：Goto3：H=M+ (L-D)F：Goto3H=M （L-D）F：Goto3LbI2：D-LT=>H=M-（D-L）C -V（T-D+ L）2/(2R)：Goto3：H=M-（D-L）C：Goto3LbI3：HGoto1三、程序执行1、CI1 前进方向输入变坡点后段坡率2、FI2 输入变坡点前段坡率以路线前进方向3、R 输入竖曲线半径。计算出T竖曲线切线长，E竖曲线外矢距。4、L 输入所求点的桩号（里程）5、D 输入变坡点桩号（里程）6、M

24、输入变坡点高程。7、H 以路线计算出所求点高程8、VTU=1,AO=-1 凸曲线时输入1，凹曲线时输入-1四、说明1、计算凸曲线，所求点桩号（里程）不小于后一个凹曲线临界（终点）桩号，不大于前一个凹曲线临界（起点）桩号。2、计算凹曲线，所求点桩号（里程）不小于后一个凸曲线临界（终点）桩号，不大于前一个凸曲线临界（起点）桩号。曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序    一、程序功能    本程序由一个主程序(TYQXJS)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线

25、、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。    特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；        

26、;      (2). 在引用该核心计算部分时，请注明来源。    二、源程序1.主程序(TYQXJS)"1.SZ => XY"："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

27、7;：N=1=>Goto 1：>Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=XY"YS"=YF"FS"=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto 22. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K

28、：X=U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD)：Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Bsin(G+QEFW(C+FWD)+Asin(G+QEMW(C+MWD)：F=G+QEW(C+WD)+90：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y

29、)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：>W=W+Z：Goto 0Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF三、使用说明1、规定   (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。   (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。   (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的4

30、5次代替。   (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。   (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。   (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线

31、的半径。2、输入与显示说明   输入部分：   1. SZ => XY   2. XY = > SZ    N ?  选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。   X0 ？线元起点的X坐标    Y0 ？线元起点的Y坐标    S0 ？线元起点里程   F0 ？线元起点切线方位角   LS ？线元长度   R0 ？线元起点曲率半径&

32、#160;  RN ？线元止点曲率半径   Q ？  线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）   S ？ 正算时所求点的里程   Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)   X ？反算时所求点的X坐标   Y ？反算时所求点的Y坐标   显示部分：   XS=×××    正算时，计算得出的所求点的X坐标  

33、; YS=×××    正算时，计算得出的所求点的Y坐标   FS=×××    正算时，所求点对应的中线点的切线方位角   S=×××    反算时，计算得出的所求点的里程   Z=×××    反算时，计算得出的所求点的边距高程源公式|编程计算器卡西欧4800 2008-10-11 01:02:

34、42  作者：LEE.CI  来源：网易博客  浏览次数：165  文字大小：【大】【中】【小】 超高及高程（源程序）L1: R:L“LS”:Z“ZH”:H“HZ”:V“RS”:F“I1”:G“I2”:O“BPH”:N“BP”:T=Vabs(GF)÷2:X“X(R+1,L-1)”:P“P(R+,L-1)”L2: Lbl 6 : SL3: S+TN<0 D=O+(SN)F: SN<0 D=O+(SN)FU(SN+T)2÷2÷V: SNT D=O+(SN)GU(TS+N)2&#

35、247;2÷V: D=O+(SN)G D“ H=” L4: R< E= : M= Goto 1: R< E= : M= Goto 1:  L5: Lbl 1 : M L C=L: Goto 2: C=M  L6: Lbl 2 : P X >0 Goto 3 : Goto 4  L7: Lbl 4 : S Z I=0.02 : Goto 5: S<Z+C I=0.02(0.02EPX)÷C×(SZ) : Goto 5: S<HC I=EPX : Goto 5: S<H I=EPX+(0.02EPX)&

36、#247;C×(SH+C) : Goto 5: I=0.02 : Goto 5L8: Lbl 3 : Z1=0.04÷(EPX+0.02)×C : S Z + Z1 I=0.02 : Goto 5: S<Z+C I=0.02+(EPX0.02)÷(CZ1)×(SZZ1) : Goto 5: S<HC I=EPX : Goto 5: S<HZ1 I=EPX(EPX0.02)÷(CZ1)×(SH+C) : Goto 5: I=0.02 : Goto 5L9: Lbl 5 : B : W“HB”=DIBL10:

37、Goto 6说明： R:平曲线半径 L“LS”:缓和曲线长 P：偏转，Z=1,Y=1X“XL”：线路（Z,Y）左线1,右线1 Z“ZH”: 直缓点里程H“HZ” ：缓直点里程 V“RS”：竖曲线半径 F“I 1” ：前坡堵（带符号）G“I2”： 后坡度值（带符号） O“BPH”：变坡点高程 N“BP”： 变坡点里程T： 切线长 E：最大超高值 C ：渐变段长度（LC） I： 横坡值 M：最大渐变长度（Lmax） B ：距中心距离 W“HB”：边部高程 字四、QXZB(曲线坐标) A”ZH(ZH)=:B”X(ZH)=:C”Y(ZH)= ”X(HZ)=:E”Y(HZ)=:F”FWJ-1=:G”FW

38、J-2=:H”R=:I”L1=:J”L=:K”L2=:L”ZJ(Z=1,Y=-1)= O=180/pi (pi为圆周率，电脑中找不到此符号以此代替) LbI A MN M”ZHUANG HAO= P=M-A MM>(A+I+J+K)?GOTO A N”ZHUANG JU(Z+,Y-)= M>(A+I+J)?GOTO D M>(A+I)?GOTO C LbI B(第一缓和曲线) X=B+(P-(P5÷40÷H2÷I2)COSF+(P3÷6÷H÷I-P7÷336÷H3÷I3)LSINF Y=

39、C+(P-(P5÷40÷H2÷I2)SINF-(P3÷6÷H÷I-P7÷336÷H3÷I3)LCOSF X=X+NCOS(F-(P2÷2÷H÷I)OL-90) Y=Y+NSIN(F-(P2÷2÷H÷I)OL-90) GOTO A LbI C(圆曲线) X=B+(HSIN(P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2)COSF-(H(1-COS(P-I)÷H+I&

40、#247;2÷H)O)+ I2÷24÷H)(-1)LSINF Y=C+(HSIN(P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2)SINF+(H(1-COS(P-I)÷H+I÷2÷H)O)+I2÷24÷H)(-1)LCOSF X=X+NCOS(F-(P-I)÷H+I÷2÷H)OL-90) Y=Y+NSIN(F-(P-I)÷H+I÷2÷H)OL-90) GOTO A LbI D(第二缓

41、和曲线) Q=A+I+J+K-M X=D-(Q-Q5÷40÷H2÷K2)COSG+(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LSING Y=E-(Q-Q5÷40÷H2÷K2)SING-(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LCOSG X=X+NCOS(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90) Y=Y

42、+NSIN(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90) GOTO A 变量说明： 以上变量依次为： ZH(ZH)=直缓点桩号，无缓和曲线则为直圆点； X(ZH)=直缓点X坐标；无缓和曲线则为直圆点； Y(ZH)=直缓点Y坐标，无缓和曲线则为直圆点； X(HZ)=缓直点X坐标，无缓和曲线则为圆直点； Y(HZ)=缓直点Y坐标，无缓和曲线则为圆直点； FWJ-1=第一切线方位角； FWJ-2=第二切线方位角； R=半径； L1=第一缓和曲线长； L=圆曲线长； L2=第二缓和曲线长； ZJ(Z

43、=1,Y=-1)=转角形式(左=1，右=-1) ZHUANG HAO=所求点桩号 ZHUANG JU(Z+,Y-)=边桩桩距，左边桩为正，右为负，中桩0 我已经用过了，没发现什么问题，请大家提意见。 1、极坐标法放样Prog：FYLb1 0：A“X0”：B“Y0”：I=0：J=0：Pol（C“XA”-A），（D“YA”-B）：J<0=>G“FW- OA”=J+360L“L0”=IGoto 1：> G“FW O-A”=JL“L0”=ILb1 1：EQ：E“Xi”：Q“Yi”：Pol（E-A），（Q-B）：J<0=>J=J+360：Goto 2：> Goto 2

44、Lb1 2： F“FW-OB”=JL=I0=F-G：O<0=>O“BJ”=O+360Goto 3：> O “BJ” Lb1 3：P=O-180Goto 1注：a、输入：（X0、Y0）、（XA、YA）测站点坐标、后视点坐标 Xi、Yi 放样点坐标b、输出：FW-OA测站至后视边方位角、L0后视边长 FW-OB测站至放样点方位角、L放样边长BJ后视边置零，放样点顺时针拨角P偏角（+为右偏、-为左偏）本值用于计算路线偏角2、公路竖曲线高程计算程序Prog:SQXLbl A:A“+（-）i1”:B“+（-）i2” W=（B-A）÷100：R：T=Abs（RW）÷2

45、：L=T*2：E=T2÷（2R）：K“JD K+”：G“JD H”：C=K-T：D=K+T：Lbl 0：J“Ki+”：J<0=>Goto 1：> Goto 2Lb1 1：“Out QX1”：H=G-(K-J)A÷100Goto 5Lb1 2：J>D=>Goto 4 W<0=>F=-1W>0=>F=1J>K=>Goto 3H=G-（K-J）A÷100+F(J-C)2÷（2R）Goto 5Lb1 3：H=G+（J-K）B÷100+F（D-J）2÷（2R）Goto 5Lb1

46、4：“OUT QX2”：H=G+（J-K）B÷100Goto 5Lb1 5：M“DHi”：H=H+M注：a、公式：L=|R(i2-i1)| 、T=L÷2、E=T2÷（2R）、h=l2÷(2R)b、功能：已知前后坡度%、竖曲线半径，计算各桩高程。c、输入：i1、i2前、后坡度%（上坡+、下坡-） JD K+、 JD H、R变坡点桩号、高程、竖曲线半径 Ki+计算桩号 DHi计算点与设计线高差（高+、低-）d、输出：H待计算点高程3、直线坐标计算Prog：ZXLb1 0：M“K0”：A“X1”：B“Y1”：C“X2”：D“Y2”：I=0：J=0：Pol（C-

47、A），（D-B）：L=IJ<0=>F”FW”=360+JGoto 1：>F”FW”=J Lb1 1：KZQ:K“Ki”：Z“Lb”：Q=90：Q“ANG”：L=K-M：Rec（L，F）：X=A+I：Y=B+J：Rec（Z，F+Q）：O=I：P=J：X=X+OY=Y+P Goto 1注：a、公式：X=X0+L*cos(Fw)、Y=Y0+L*sin(Fw)b、功能：已知直线起终点坐标，计算中间桩中边桩坐标c、输入： K0、（X1、Y1）、（X2、Y2）起点桩号、起点坐标、终点坐标Ki、Lb、Ang 待求桩号、边桩距（左-、右+）、边桩与中线右夹角（90）d、输出：L、FW直线长度

48、、起点至终点方位角 X、Y计算点坐标a多边形面积计算（一）Prog：AreaLb1 0：S=0：L=0：E“X0”：F“Y0”：A=E：B=FLb1 1：CD：C“Xi： XiX0，Yi=Y0 END”：D“Yi”：M=（C-E）2+（D-F）2）：S=S+（AD-CB）：L=L+（C-A）2+（D-B）2）：A=C：B=D：S=0Goto 1：Lb1 3：M0=Goto 1：S“Area”=Abs S÷2L注：a、公式：Area =（XiYi+1-Xi+1Yi）÷2b、功能：已知多边形各转点坐标，计算多边形面积。c、输入：X0、Y0起点坐标 Xi、Yi各点坐标（按顺序输入

49、、最后要闭合即最后一点要输入起点坐标）d、输出：Area面积 L周长4、对称标准曲线中边桩坐标计算Prog：QXLb1 A: M“JDK+”:A“PJ=>R（+），L（-）”:R：W“LS”: A0=>G=1:Goto E: >G=-1:Lb1 E:Defm 10:A=Abs A:B=90W÷R÷：P=W2÷（6R）-W4÷（336R3）-R（1-cosB）:Q=W-W3÷(40R2)-RsinB：T=（R+P）tan（A÷2）+Q：L=RA÷180+W:E=(R+P)÷cos(A÷2)

50、-R:N=M-T:Z5=N+W:Z=N+L÷2:C=N+L-W:D=N+L:A=AG:U:U“ZD-2,XY-1”=2=>TELN“ZH”Z5“HY”Z“QZ”C“YH”D“HZ”Goto A：> Goto B Lbl B：F“FW”：V“JD-X”：H“JD-Y”：I=0：J=0：Rec（T，F+180）：Z1“ZH X”=I+V：Z2“ZH Y”=J+H：Rec（T，F+A）：Z3“HZ X”=I+V：Z4“HZ Y”=J+H：L=0Lbl 0：E=90：KSE：K“JS K+”：S“LB”：E“ANG”：S=0=>U=1：Goto C：>U=S÷

51、;（Abs S）Lbl C：S=Abs S：W=0=> Goto D：>K<N=> Goto 1：>K<Z5 => Goto 2：>K<D=> Goto 4：> Goto 5Lbl D：K<N=>Goto 1：>K<D=>Goto 3：> Goto 5Lbl 1：Rec（M-K，F+180）：X=V+I：Y=J+H：Rec（S，F+180-（180-E）U）：X“Xi”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 2：Z=K-N：O=90Z2÷R÷W÷：

52、X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：Lbl 6：Rec（X，F）：X=Z1+I：Y=Z2+J：Rec（Z，F+90G）：X=X+I：Y=Y+J：Rec（S，F+OG+EU）：X“Xi”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 3：Z=K-Z5：O=（K-N-W÷2）*180÷÷R：X=Q+Rsin（B

53、+180Z÷÷R）：Goto 6Lbl 4：Z=D-K：O=90Z2÷R÷W÷：X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：Rec（X，F+A+180）：X=Z3+I：Y=Z4+J：Rec（Z，F+A+180-90G）：X=X+I：Y=Y+J：Rec（S，F+A+180-OG-（180-E）U）：X“X

54、i”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 5：Rec（K-D+T），（F+A）：X=V+I：Y=H+J：Rec（S，F+A+EU）：X“Xi”=X+IY“Yi“=Y+J S=SU： Goto 7Lbl 7：Goto 0注：输入：JDK+、PJ、R、LS、交点桩号、路线偏角（右+、左-）、半径、缓和曲线长“ZD-2,XY-1”要显示主点要素输入2，否则进入坐标计算FW、（JD X、JD Y）起始边方位角（即ZHJD方位角）、交点坐标JS K+、LB、Ang待计算桩号、边桩距（右+、左-）、边桩与路中线右夹角（默认90度）输出：T、E、L、ZH、HY、YH、HZ、曲线要素、 Xi

55、,Yi待求点坐标5、公路任何曲线通用中边桩坐标正反算（复化辛卜生公式也称积木法）Prog：HQXDefm 4V“V=1=>XY： 2=>Ki”Lbl 0：KLWLbl 4：Q“one=0：1、2” =0=> Goto 1Q=1=>Prog“ZD”（线元要素数据文件名） Q=3=>Prog“要素表文件名” （注：本句是为了全线及多条线路计算用，如仅单个可省略本句）Lbl 1：A"K1"：B"K2"：C"FW0"：D"R1"：E"R2"：F"X0"：

56、G"Y0"D0=>I=1÷D：=>I=DE0=>J=1÷E：=>J=EAbsD+AbsE=0=>M=2：=>M=16V=2=>L=0：W=90K“Ki”：L“L-BZ”：W"ANG"：N=0：Z1=0：Z2=0：Z3=0：Z4=0 进入坐标迭代计算Lbl 2：N=N+1：H=2(K-A)÷M：R=NH÷2+A：R=C+180÷×(I+(J-I)÷2÷(B-A)×(R-A)×(R-A)：Int(N÷2)=N

57、÷2=>Z1=Z1+cosR：Z2=Z2+sinR：=>Z3=Z3+cosR：Z4=Z4+ sinRN=M=>Goto 3：>Goto 2Lbl 3：X=F+H÷6×(cosC+4Z3+2Z1-cosR)+Lcos(R+W)：Y=G+H÷6×(sinC+4Z4+2Z2-sinR)+Lsin(R+W)：V=2=>Goto 6X"Xi="Y"Yi="R"FW"=R-360Intg(R÷360) Q=0=>KLW：Goto 1 ：> Got

58、o 0Lbl 6： Pol(T"XF"-X,U"YF"-Y) 进入桩号求算O=Icos(J-R)：AbsO0.0001=>O"DL"K=K+O O"LP"=Isin(J-R)TU：Goto 6：=>K=K+O：L=0：Goto 4线元要素数据文件格式：K>起点桩号=> A"K1"=？：B"K2"=？：C"FW0"=？：D"R1"=？：E"R2"=？：F"X0"=？：G"Y0"=？：K> B=>A=B：B=？（下曲线要素）注：a、功能：利用公路通用复化辛卜生积分公式进行中边桩坐