圆曲线上坐标计算程序

坐标计算程序及坐标计算公式一、坐标计算公式直线段：中桩公式：x=x1+(z-c)*cosay=y1+(z-c)*sina边桩公式：x=x1+(z-c)*cosa±d*cos(a-90)y=y1+(z-c)*sina±d*sin(a-90)说明： x1——起点x坐标，y1——起点y坐标，z——计算点桩号， c——起点桩号，a——方位角，d——距中桩距离。

“±”左边桩为“+”，右边桩为“-”。

二、方位角计算公式1、直线段方位角图纸提供。

2、若为单一圆曲线时，起点用直线段方位角图纸提供。

3、若为缓和曲线时：第一段缓和曲线方位角为直线段方位角图纸提供；第二段圆曲线起点方位角a1=a0±Ls*180/2∏r。

第三段缓和曲线方位角a2=a1±L*360/2∏r。

（a0为直线段方位角，Ls为缓和曲线长度，L为圆曲线长度，r为圆曲线半径，“±”右偏角为“+”、左偏角为“-”。

）三、5800计算器坐标程序坐标程序由1个主程序、5个子程序和1个数据库组成进入计算器编程模式（5：PROG—1：NEW新建程序名—3：EDIT），输写程序。

1、主程序adminFix 3（回车换行）Lb1 0:150→DimZ:“ZX→0,A→1,B1→2,B2→3,C→4,D→5,CR→6”?N:N→Z[149]:Prog “DAT2”：“ZS=>1，FS=>2”?N:If N=1:Then Goto 1:Else N=2=>Goto 2:IfEnd:Goto 0:Lb1 1:?S:“ANG=”?M:?Z:S=0=>Goto 0:0→N（回车换行）Lb1 5:Isz N:If S≤Z[8N+2]+Z[8N+4]:Then N→J:Prog“DAT1”:Else Goto5:IfEnd（回车换行）Abs(S-0)→W:Prog“SUB1”：“XS=”：X◢“YS=” ：Y◢“FWI=”：F-M→F：If F≤360:Then F→F:Else F-360→F:IfEnd:F►DMS◢Goto4（回车换行）Lb1 2:?X:?Y:“ANG=”？M：M→Z[148]：If M＜90:Then 180-M→M:IfEnd（回车换行）X-Z[4]:Y-Z[5]:X=0=>Goto 0:0→N（回车换行）Lb1 A:Isz N:N→Z[150]:Z[8N+3]-M→A:Z[8(N+1)+3]-M→B:Prog “ZX1”:If Z[6]×Z[7]≤0:Then N→J:Prog “DAT1”:Goto B:IfEnd（回车换行）Z[8N+3]+M→A:Z[8(N+1)+3]+M→B:Prog “ZX1”:If Z[6]×Z[7] ≤0:Then N→J:Prog “DAT1”:Goto B:Else Goto A: IfEnd（回车换行）Lb1 B:Prog “SUB2” （回车换行）Z[150]→N:0+W→S:If S>Z[8N+2]+Z[8N+4]+.001:Then Goto A: IfEnd（回车换行）If N>13:Then 0→N: Goto A: IfEnd（回车换行）If Z[148]>90: Then S+2Zsin(M-90) →S: IfEnd:“S=”:S◢“Z=”: Z◢“OK→2,NO→1”?N:If N=1: Then Z[150]→N: Goto A:Else Goto 2: IfEnd（回车换行）Lb1 4 :0→J:“DIST=”：Pol （X－Z[1]，Y－Z[2]）→I◢ J-F：If F<0：Then F+360→F：IfEnd：“FW=”： F►DMS◢ Goto 12、子程序DAT1Z[8J]→U：Z[8J+1]→V：Z[8J+2]→O：Z[8J+3]→G：Z[8J+4]→G：Z[8J+4]→H：Z[8J+5]→P：Z[8J+6]→R：Z[8J+7]→Q：1÷P→C：（P－R）÷（2HPR）→D：（180÷∏）→E：Return3、子程序ss“A”？→A：If A<0:Then 10A◢Else 9A◢IfEnd:Ans×1.05（回车）4、子程序SUB1.1184634425→A：.2393143352→B：.28444444444→N：.046910077→K：.2307653449→L：.5→Z[3]:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C +Z[3]WD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD)))→X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Nsin(G+QEZ[3]W(C+Z[3 ]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD))) →Y：G+QEW(C+WD)+M→F：X+Zcos（F）→X：Y+Zsin（F）→Y: Return5、子程序SUB2G-M→T：Abs（（Y-V）cos（T）-（X-U）sin（T））→W：0→Z：Lb1 0：Pros “SUB1” （回车换行）T+QEW（C+WD）→L：（Z[5]-Y）cos（L）-（Z[4]-X）sin（L）→Z：If Abs（Z）<1×10∧(-4):Then Goto 1:Else W+Z→W: Goto 0: IfEnd（回车换行）Lb1 1:0→Z:Prog“SUB1”：（Z[5]-Y）÷sin（F-2M+180）→Z：Return6、子程序ZX1（Z[5]-Z[8N+1]）cos（A）-（Z[4]-z[8N]）sin（A）→Z[6]：（Z[5]-Z[8（N+1）+1]）C5]-Z[8（N+1）+1]）cos（B）-（Z[4]-Z[8（N+1）]sin（B）→Z[7]：Return7、数据库DAT2If N=0:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）If N=1:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）说明：1、所有路线都是从Z[8]开始。

CASIOfx-4500PA计算器程序和操作一、水准测量高程计算程序1、程序名：H-（1）第一行：L1 lbl 0 第二行：L2 Z：A：｛BK｝第三行：L3 H=Z＋A－B◢第四行：L4 V=K-H ◢第五行：L5 Goto 02、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示H-（1）文件名为止：③按EXE键，显示Z？，输入水准点BM后视点高程：如129.919；④按EXE键，显示A？，输入后视水准尺读数：如1.633；⑤按EXE键，显示B？，输入前视水准尺读数：如4.186；⑥按EXE键，显示H=127.366，此结果是K3＋300中桩地面髙程；⑦按EXE键，显示K？，输入K3＋300中桩设计高程127.032；⑧按EXE键，显示V=K-H，-0.334此结果是K3＋300中桩挖出数据；⑨按EXE键，显示B？，连续算下去……（－挖；＋填）。

二、坐标反算程序①1、程序名；1ZBFS第一行：L1 AB：Fixm：｛CD｝第二行：L2 Pol(C-A，D-B）◢第三行: L3 W<0 =〉W=W+360第四行：L4 lntW＋0.01lnt（）＋0.006 Frac（60 Frac W）◢2、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示1ZBFS文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567（也可当测站）；④按EXE键，显示B？，输入第一点y坐标：如814.454（也可当测站）；⑤按EXE键，显示C？，输入第二点x坐标：如406.260 （也可作后视点）；⑥按EXE键，显示D？，输入第二点y坐标：如1029.145（也可作后视点）；⑦按EXE键，显示Pol(C-A，D-B）：218.89476为两点间距离；⑧按EXE键，显示lntW＋0.01lnt（60 Frac W）＋0.006 Frac（60 Frac W）78.45108两点方位角；(即是78°45′11″)三、坐标正算程序1、程序名；2ZBZS第一行：L1 AB：Lbl5：｛ST｝第二行：L2 C=A＋Rce（S，T）◢第三行：L3 D=B＋W◢第四行：L4 Goto52、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示2ZBZS文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567；④按EXE键，显示B？，输入第一点y坐标：如814.454；⑤按EXE键，显示S？，输入两点间距离：如218.89476；⑥按EXE键，显示T？，输入两点方位角：如78°45′10.8″；⑦按EXE键，显示C=A＋Rce（S，T），为第二点x坐标：如406.260；⑧按EXE键，显示D=B＋W，为第二点y坐标：如1029.14499；四、坐标反算程序②1、程序名；ZF第一行：L1 Lbl 0第二行：L2 A：B：｛CD｝：C≤0 G=〉oto2第三行：L3 X=C－A第四行：L4 Y=D－B第五行：L5 S=Pol（X，Y）◢第六行：L6 T=W第七行：L7 T＜0=〉第八行：L8 T=W ◢第九行：L9 Goto 0第十行：L10 Lbl 1第十一行：L11 T″T″=360＋T◢第十二行：L12 Goto 0第十三行：L3 lbl 2第十四行：L14 ｛AB｝第十五行：L15 Goto 02、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示ZF文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567 ；④按EXE键，显示B？，输入第一点y坐标：如814.454；⑤按EXE键，显示C？，输入第二点x坐标：如354.618；⑥按EXE键，显示D？，输入第二点y坐标：如553.341；⑦按EXE键，显示Pol(C-A，D-B）S=Pol（X，Y）：261.2663为两点间距离；⑧按EXE键，显示T=268.0371接着按SHIFT °′″键显T=268°2′13.55″为两点的方位角；⑨按EXE键，显示C？，输入下一点x值；⑩重复以上操作方法（略）五、线路直线段点坐标计算程序(1) ㈠本可计算前切线YZ至ZY点(或后切1、程序名；Z-X,Y(1) 线ZY至YZ)间直线上任点中、边桩坐标。

● 圆曲线方法一：sin (1cos )180i i i i i i x R y R l R ϕϕϕπ⎧⎪=⎪=-⎨⎪︒⎪=⋅⎩——i l 为待定点i P 至起点间的弧长i ϕ为i l 所对的圆心角R 为曲线半径方法二：11802l A R π︒=⋅⋅ 2sin l R A =⋅00cos(/)sin(/)x x l A y y l A =+⋅+-⎧⎨=+⋅+-⎩起点方位角左减右加起点方位角左减右加——00(,)x y 为圆曲线起点坐标方法三：180l A R π︒=⋅ 00cos(/)sin(/)x x R B A y y R B A =+⋅+-⎧⎨=+⋅+-⎩——l 为圆曲线上任意一点距起点距离00(,)x y 为圆曲线圆心坐标B 为圆心到圆曲线起点的方位角，A 为任意点对应的圆心角● 缓和曲线522030406l x l R l ly Rl ⎧=-⎪⎪⎨⎪=⎪⎩——l 为曲线上任一点至起点的曲线长R 曲线半径0l 为缓和曲线全长圆曲线、缓和曲线计算方法1、直线段：先由JD1以及JD2的坐标算出JD1到JD2的方位角，即直线段方位角A ，故可算出HZ 、ZH 坐标及其直线段各点坐标。

2、缓和曲线：以HZ 、ZH 为起点，缓和曲线上任意一点离HZ 、ZH 距离为l ，利用公式522003040()6l x l R l l R ly Rl ⎧=-⎪⎪⎨⎪=⎪⎩为缓和曲线全长,为圆曲线半径算出该点的相对起点坐标，利用arctan y x算出该点相对起点的方位角B ，再根据线路走向及直线段方位角可算出该点的方位角C （顺时针加，逆时针减），用可求出该点相对起点的距离D ，最后用00cos sin x x D C y y D C =+⎧⎨=+⎩可求出该点的坐标。

（00(,)x y 为缓和曲线起点的坐标）3、圆曲线：用上述方法求出圆曲线两端点HY 、YH 坐标，算出HY 到YH 的方位角F ，以及两点间的距离E ，用12arccos ER可算出两端点连线与起点到圆心连线的夹角G ，根据线路走向求出起点到圆心的方位角H （H=F+/-G ），00(,)x y 圆曲线为起点坐标，根据00cos sin x x R H y y R H=+⎧⎨=+⎩，求出圆心坐标。