曲线任意里程中边桩坐标正反算及放样fx
曲线任意里程中边桩坐标正反算及放样fx-4850程序
本程序是在yshf及hangniu1973两位师傅的成果上作的一次改动，修改内容主要有一下几点：
1、在变量符号上修改，目的是更能与工地实际结合；
2、把原主程序分为两个，更能灵活调用；
3、加入放样程序，做到坐标计算到放样一体化；
4、使得整个测量放样过程更加简单。
GAUSSLE坐标正反算fx－4850程序
源程序
1.正算主程序 GSZS
I"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” :
D=(P-R)÷(2(H-O)PR)：
KL”L(-Z +Y)” :M”ANG(YJJ)”=90:（注：此处若不给M赋值，则可计算斜交点）
J=Abs(K-O)：Prog"SUB1"：
”FWJ=”:F=F-M:”X=”:X=X:Pause0: ”Y=”:Y=Y◢
2. 反算主程序 GSFS
XY：Z[2]=X：Z[3] =Y：
I"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” :
D=(P-R)÷(2(H-O)PR):
J=Abs((Y-S)cos(G-90)－(X-I)sin(G-90))：
L=0：M”M(YJJ)”=90:
Lbl 0：Prog "SUB1"：
L=(Z[3]-Y)cos(G-90+QJ(1÷P+JD)×180÷π)-(Z[2]-X)sin(G-90+QJ(1÷P +JD) ×180÷π)：AbsL<1E-6=>Goto1：≠>J=J+L：Goto 0Δ←┘
Lbl 1：L=0：Prog "SUB1"：L=(Z[3]-Y)÷sinF:
”K=”:K=O+J:Pause0:”L=”:L=L◢
3. 正算子程序(SUB1)
Defm 4:
A=0.1184634425:B=0.2393143352:Z[4]=0.2844444444:C=0.0469100770:E=0.2307653449: Z[1]=0.5:
X=I+J(Acos(G+QCJ(1÷P+CJD)×180÷π)+Bcos(G+QEJ(1÷P+EJD)×180÷π)+Z[4]cos(G+QZ[ 1]J(1÷P+Z[1]JD)×180÷π)+Bcos(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Acos(G+Q
(1-C)J(1÷P+(1-C)J D) ×180÷π)):
Y=S+J(Asin(G+QCJ(1÷P+CJD)×180÷π)+Bsin(G+QEJ(1÷P+EJD)×180÷π)+Z[4]sin(G+QZ[1 ]J(1÷P+Z[1]JD)×180÷π)+Bsin(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Asin(G+Q
(1-C)J(1÷P+(1-C)JD) ×180÷π))：
F=G+QJ(1÷P+JD) ×180÷π+M：X=X+LcosF：Y=Y+LsinF
4. 曲线元要素数据库：DAT-M
K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
K≥O=>K≤H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
……………………………
K≥O=>K≤H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘
(注：如有多个曲线元要素继续添加入数据库DA T-M中)
5、M线（坐标正算）组合程序 MG-ZB
Prog”DAT-M”:Prog”GSZS”
6、M线（坐标计算－放样）组合程序 MG-FY
Prog”MG-ZB”:Prog”LTKZD”: Prog”FY”
7、M线（坐标反算）组合程序 M-GSFB
Prog”DAT-M”:Prog”GSFS”
说明：
一、程序功能及原理
1.功能说明：
本程序由两个主程序——正算主程序(GSZS)、反算主程序(GSFS)和两个子程——正算子程序(SUB1)、线元数据库(DA T-M)构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本程序可以在CASIO fx-4800P计算器及CASIO fx-4850P计算器上运行。由于加入了数据库(DAT-M)，可实现坐标正反算的全线贯通。
组合程序5可实现M线的正算贯通，组合程序7可实现M线的反算贯通，组合程序6可实现坐标计算到放样一体化。
2．计算原理：
利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。
利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。
二、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，L=0；当位于中线左侧时，L取负值；当位于中线右侧时，L取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
(7)曲线元要素数据库（DAT-M）可根据线型不同分为各个线元段输入到DAT-M中，即分为直线段、缓和曲线、圆曲线等。
(8)正算时可仅输入里程和边距及右交角可实现全线计算，但反算时只能通过首先输入里程K值读取数据库DA T-M，计算器自动将里程K所在线元数据赋给反算主程序GSFS 进行试算，试算出的里程和边距须带入正算主程序GSZS中计算坐标，若坐标吻合则反算正确。
2、输入与显示说明
（1）输入部分：
X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
K0 ？线元起点里程
F0 ？线元起点切线方位角
KN ？线元终点里程
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
Q ？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0)
K ？正算时所求点的里程
L ？正算时所求点距中线的边距(左侧取负值，右侧取正值，在中线上取零) ANG？正算边桩时左右边桩连线与线路中线的右交角
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
M ? 斜交右角
线元要素数据库中K≥O=>K＜H=>中的O和H分别为该段线元起点里程和终点里程A、B、Z[4] 是Gauss-Legendre求积公式中的插值系数
C 、E、Z[1] 是Gauss-Legendre求积公式中的求积节点
（2）显示部分：
X=×××正算时，计算得出的所求点的X坐标
Y=×××正算时，计算得出的所求点的Y坐标
K=×××反算时，计算得出的所求点的里程
L=×××反算时，计算得出的所求点的边距
三、算例
某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：
S0 X0 Y0 F0 LS R0 RN Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0 769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1 806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1 919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1 999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 （注：该算例中线元要素Ls为程序修改前须输入的线元长度，程序修改后改为输入线元终点里程KN）
放样程序 FY
U”XJ”:V”YJ”:W”XH”:Z”YH”:XY:I=W-U:J=Z-V:Pol(I,J):J<0=>J=J+360Δ
“HSJL=”:I◢
O”HF”=J:X=0=>U=U+LCos(J+P”μ”):
V=V+LSin(J+P”μ”):
”X=”:U:Pause0:”Y=”:V:≠>I=X-U:J=Y-V: Pol(I,J):J<0=>J=J+360:Δ
H”QF”=J:J=H-O:J<0=>J=J+360:Δ
“ANG=”:J→DMS （该处不输Pause0语句亦可实现“ANG=”和“QSJL=”的同时出现）
“QSJL=”:I
凌铁大桥控制点数据库 LTKZD
A=21:B=22:C=23:D=24:E=25:F=26:G=27:H=28:I=29←┘
T “ZJDH”:←┘
T=11=>U=****:V=****⊿←┘
T=12=>U=****:V=****⊿←┘
T=13=>U=****:V=****⊿←┘
T=15=>U=****:V=****⊿←┘
T= 4 =>U=****:V=****⊿←┘
T= A=>U=****:V=****⊿←┘
T= B=>U=****:V=****⊿←┘
T= C=>U=****:V=****⊿←┘
……………………………
(注：如有多个控制点继续添加入数据库LTKZD中)
N “HSDH”:←┘
N=11=>W=****:Z=****⊿←┘
N=12=>W=****:Z=****⊿←┘
N=13=>W=****:Z=****⊿←┘
N=15=>W=****:Z=****⊿←┘
N= 4 =>W=****:Z=****⊿←┘
N= A=>W=****:Z=****⊿←┘
N= B=>W=****:Z=****⊿←┘
N= C=>W=****:Z=****⊿←┘
……………………………
(注：如有多个控制点继续添加入数据库LTKZD中)
组合程式 LTFY
Prog ”LTKZD”: Prog ”FY”
说明：
一、程序功能及原理
1、功能说明：
FY程序有两个功能，当前视点坐标X输入为零时，程式执行转点功能，即首先通过全站仪测定任意转点距置镜点距离L及置镜点与转点连线和置镜点与后视点连线的夹角P”μ”并输入计算器，即可计算出任意转点坐标；另一功能是坐标放样，即输入放样点X、Y坐标，从而计算出与置镜点距离“QSJL=”以及全站仪归零放样拨角“ANG=”（前视方位角－后视方位角）。
凌铁大桥控制点数据库LTKZD功能是与FY程序配合后，可实现坐标放样时仅输入置镜点、后视点的编号，从而减少其数据的大量输入。
2、计算原理：
利用坐标计算和方位角的基本知识以及计算器极坐标转换功能键POL
二、使用说明
1、规定
当前视点坐标X输入为零时，程式执行转点功能
2、输入与显示说明
（1）输入部分：
XJ ？置镜点X坐标
YJ ？置镜点Y坐标
XH ？后视点X坐标
YH ？后视点Y坐标
X ？放样点X坐标
Y ？放样点Y坐标
L ？实测转点距置镜点距离
P”μ”？实测置镜点与转点连线和置镜点与后视点连线的夹角
T ？置镜点点号ZJDH
N ？后视点点号HSDH
11、12、13、15、4、A、B、C、D、E、F、G、H、I分别为凌铁大桥控制点点号，计算输入时按以上数字或字母输入计算器即可
(2)显示部分
HSJL=**** 后视距离（置镜点与后视点之间距离）
ANG= **** 全站仪归零放样拨角
QSJL=**** 前视距离（置镜点与放样点之间距离）
X= **** 转点X坐标
Y= **** 转点Y坐标
三、其它说明
1、若需要编制另一线路控制点程序，仅需新建另一数据库子程序LTKZD，新建时仅需改变上文中红色显示数据即可。
2、程式中11~15、及A~I均为控制点编号，程式开头首先给A、B赋值：A=21:B=22等，目的是为了在程式运行要求输入点号时能直接输入字母A、B、C等。
坐标反算
在现场工作中，以往我们都是已知某点的里程及边距，来计算出该点的坐标，但有时我们如果能在测得某点坐标后，计算出该点的里程和距线路中心的距离(在这里我姑且称之为坐标反算)的话，将会帮我们大大减轻野外工作量，提高我们的工作效率。例如：路基填了几层后要精确检查一下路基是否够宽，那么按照我们以往的做法，就是要先将线路中心线放出来，然后用尺拉一下路基宽度，与其在此高程的设计宽度作比较，这样做对高填方而言极不方便。或者是先按所测高程,计算的宽度放出路基边桩，再与所填边线作比较。以上两种方法现场工作量都比较大。较为简便的方法是，我们可以测一下已填路基边线上任一点的三维坐标，然后将其反算求出该点的里程，及其距中线的距离(即所填宽度),由计算出的里程，可算出该里程的路面设计高程，再有所测高程，可计算出该点的设计宽度，两宽度作比较即可。同样在桥面铺装施工时，我们也无须再像以往那样，先放出某点再测其高程，然后与设计高程比较计，算出该点铺装厚度，而可以沿桥面外边线随意布点，测其三维坐标，计算出其里程及到中心线的距离，便可由其里程及距中心距离，计算出该点的设计高程，与其测得高程作比较得出应铺厚度。这样便大大减轻了外业工作强度(由放出点后再测其高程，变为测任意点高程)，而内业计算量与常规相当。另外在临时增加桥涵时，也常用到此方法来计算变更桥涵的中心里程(斜交或正交均可).如目前我标段就存在很多临时变更涵洞,按以往我们的方法是先估计该处大概里程,然后放出所估计里程的中心桩,再用皮尺量出所要增加涵洞处与该中心桩的距离,以此来推算出涵洞的中心里程,这一过程即繁琐又不准确。而目前我们采用的方法是用全站仪测得跨路基现有水沟两端的沟底坐标,计算出其与路基的夹角,按所测坐标及此夹角就可以准确、快速地反算出水沟中心所对应的线路中心里程了。我们在日常测量工作中的很多方面，也会用到这一方法来减轻野外工作量。在目前我标段的S334分离式立交桥的架设过程中,也同样用到了此方法.支座安装好后,对支座中心位置检及高程查无误后开始梁板架设,但是尽管测量控制放样符合规范要求,可是因为其它方面的各种原因可能会使梁板出现偏位高程也可能会出现偏差,那么对现在这种问题该如何检查呢?其实方法是一样的,首先我们可以用全站仪测得架设好后梁的边板外边缘任一点的三维坐标,由此坐标反算出该点所对应的中心里程和距中心的距离,就可以和设计图纸上的距离作比较来检查其是否存在偏位,该点的设计高程也可以由反算所得的中心里程和距中心的边距算出,与所测得的实际高程作一下比较也就可以了.那么通过以上讨论问题归结到了一点，那就是如何在测得任一点坐标后，计算出其所对应的线路中心里程，及其到线路中心的距离（或是斜交的长度）呢？解决此类问题，对目前一些测量软件来说早已不成问题，但是在现场工作中我们用的更多、更方便的还是计算器，那么能否用我们常用的4800或5800计算器编程，来计算此类问
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。 2、线元数据可自动从数据库调用，也可手工输入。 3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。 4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。 5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。第三次及以后试算才要求输入近似里程。 6、程序代码规范简洁，便于阅读、理解。 完整程序清单: ZFS %正反算主程序 B=.1739274226:C=.5-B: Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS": ≠>U=1=>Prog"FS":≠>Goto 1
ZS %正算子程序 {K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB" FS %反算子程序 {KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2Δ M=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1 ZZ %高斯法中桩子程序(4节点) Prog"XL":M=K-L:O=(P-R)÷2PQR: D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D: I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1: S=X+M(BcosD+CcosE+CcosG+BcosH: T=Y+M(BsinD+CsinE+CsinG+BsinH WY %外移点计算子程序 Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1 WZ %位置显示子程序 "KJ":K:Pause 1:J◢ ZB %坐标显示子程序 "XY":S:Pause 1:T◢ YC %异常处理子程序 U=1=>K=L:U=2Δ U=3=>K=M:U=4Δ
道路坐标计算公式
曲线坐标计算 1、曲线要素计算 （1）缓和曲线常数计算 内移距R l 24/p 2 s = 切垂距 23 s 240/2/m R l l s -= 缓和曲线角R l R l s s πβ/902/0??== （2）曲线要素计算 切线长 m R T ++=2/tan )p (α 曲线长 ?+=?-+=180/]180/)2([20απβαπR l R l L s s 外矢距 R R E -+=)]2/cos(/)p [(0α 切曲差 L T q -=2 2、主要点的里程推算
s s s S l YH HZ )/22l -(L QZ YH )/22l -(L HY QZ l +=+=+=+=-=ZH HY T JD ZH 检核： HZ T JD =-+q 3、方位角计算 根据已知JD1和JD2的坐标计算出 21JD JD -α 偏角βαα±=--211JD JD JD ZH ?±-=-18011JD ZH ZH JD αα 4、计算直线中桩坐标 (1)计算ZH 点坐标： ZH JD JD ZH ZH JD JD ZH T y y T x x --?+=?+=1111sin cos αα (2)计算HZ 点坐标： 2 11211cos cos JD JD JD HZ JD JD JD HZ T y y T x x --?+=?+=αα (3)计算直线上任意点中桩坐标 待求点到JD1的距离为i L 2 112 11sin cos -JD JD i JD i JD JD i JD i i L y y L x x HZ T L --?+=?+=+=αα里程 待求点里程 5、计算缓和曲线中桩坐标 (1)第一缓和曲线上任意点中桩坐标 在切线坐标系中的坐标为： s i s i Rl l y Rl l l x 6/)(40/3 25=-= ZH 到所求点方位角：
5800简单全线坐标计算程序
5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N)
卡西欧5800公路坐标正反算程序教程文件
卡西欧5800公路坐标正反算程序
目录 一、坐标正算基本公式 (02) 二、坐标反算原理 (04) 三、高程数据库录入变换 (05) 四、计算器程序 (07) 01、ZBZS(坐标正算) (07) 02、ZBFS（坐标反算） (08) 03、GCJF(高程积分) (09) 04、PJFY（坡脚放样） (10) 05、JFCX（积分程序） (11) 06、ZBFY（坐标放样） (11) 07、DT（递推） (12) 08、HP（横坡） (13) 09、LK（路宽） (14) 10、SJK1（平面数据库） (14) 11、SJK2（纵面数据库） (14) 12、SJK3（左路宽度数据库） (15) 13、SJK4（右路宽度数据库） (15) 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库） (16) 16、SJK7（左上边坡数据库） (17) 17、SJK8（右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序
一、正算所涉及的计算公式 X R i d X αβB d Y d l d βI 图表 1 在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标 点。设： A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为 ρA ,A 点的里程为L A ,B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的 里程为L I 。I 点的切线角为β。 由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。故此任 意点的曲率为; C L R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为; A B A i A B A i L L L L --?-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。设：
道路逐桩坐标计算
道路逐桩坐标计算[可读数据库文件(.mdb)，或读文本文件(.txt或.dat) 附件(点击下载)： ;;; by yshf ;;;道路逐桩坐标计算[可读数据库文件(.mdb)，或读文本文件(.txt或.dat)] ;;;1. 根据“道路设计参数文件”[.txt或.dat(文本文件)， ;;; 或者.mdb(Access 2000 数据库)] ”中的平面曲线线元参数、 ;;; 道路纵断面参数成批地计算所求点坐标和相应中线点的设计高程， ;;; 并在Auto CAD中绘制出逐桩坐标表。 ;;; ;;;2. 必须将下载的文件“zbjgchjsb1.fas”存到“E:\\算例文件夹”中， ;;; 如存入其它地方，则程序不会进行计算。 ;;; ;;;3. 运行环境为：Auto CAD 2000以上版，Access 2000以上版数据库。 ;;; ;;;4. 计算前，先准备数据： ;;; (一)平面曲线 ;;; 平面曲线按线元法将各线元要素录入到Access 2000以上版数据 ;;; 库的“道路平面曲线线元参数表”中，或者录入到文本文件(.txt或 ;;; .dat)。当曲线左偏时，其线元长度输入负值；右偏及直线时其线元 ;;; 长度输入正值。 ;;; 起点切线方位是以度.分分秒秒的形式录入的，例如57°09′13.32″ ;;; 录入为57.091332。 ;;; ;;; (二)平曲线曲率半径约定如下： ;;; (1).当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次;;; 代替。 ;;; (2).当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半;;; 径均等于圆弧的半径。 ;;; (3).当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，;;; 以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直;;; 线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半;;; 径等于圆曲线的半径。 ;;; (4) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计;;; 规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接;;; 时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的;;; 半径。 ;;; ;;; (三)竖曲线 ;;; 竖曲线按变坡点里程、变坡点高程、竖曲线半径的方式录入到 ;;; “道路纵断面参数表”中，在变坡点未设有竖曲线的，其竖曲线半径 ;;; 输入0。 ;;; ;;; (四)注意事项
线元法万能坐标计算程序
线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.wendangku.net/doc/2714171227.html,/：本论文仅供学习交流使用，本站仅作合理转载，原作者可来邮要求删除论 文。 摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。利用CASIO fx-9750GⅡ编程函数计算器强大的内存（可诸存63000个字符）和编程功能，编写各种计算程序，能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料，同时也使我们有更多的时间去挑战更富有创造性的工作。 关键词：坐标放线线元测量程序 1、前言 本程序采用Gauss-Legendre（高斯-勒让德）五节点公式作内核,计算速度（太约2秒）适中，计算精度很高。在此之前，本人曾用过以下公式作内核：①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点，⑦高斯-勒让德5节点，经过测试③计算最快，⑦代码稍长但计算速度只比③⑥稍慢，精度最高，可满足线元长小于1/2πD 的所有线形的精度要求。⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55m m,3/4πD时偏差为31.63mm，πD时偏差为968mm，偏差按半径倍数增大，如线元长大于1/2πD（1/2圆周长）时，可将其拆分二个或多个线元单位，以确计算保精度。 2、程序特点 事先将所有的平曲线交点的线元要素诸存到计算器内，测量时只输桩号、边距等程序会自动寻找各类要素，一气呵成地完成施工测量任务，中途不需人工转换各类要素数据，本程序可诸存几百条线路的要素数据，计算时可按需选择线路编号进行测量。测量时不需查阅及携带图纸，仅一台CASIO fx-9750GⅡ编程函数计算器即可。 本程序含一个主程序:3XYF，五个子程序:GL（公式内核）、QD（线路选择）、XL（线路要素判断）、GF（坐标反算）、File 1 （要素存放的串列工作簿）。可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正反算。 3、计算公式及原理 如图：BC 间为一曲线元，曲线元上任一点的曲率随至B 点的弧长作线性变化。设起点B 的曲率为KA ，终点C 的曲率为KB ，R 为曲线半径。±表示曲线元的偏向，当曲线元左偏时取负号，当曲线元右偏时取正号,直线段以1的45次方代替（即半径无穷大）。 式中：αΑ＝起始方位角l ＝p 点到B的距离lS=曲线总长αp＝p 点切线方位角 R1=R5=0.118463442528095 ,R2 = R4 = 0.239314335249683 , R3 = 0.28444444444444 V1=1-V5= 0.046910070 ,V 2= 1-V4 = 1 0.2307653449 V3= 0.5 利用上面公式及CASIO fx-9750GⅡ编程函数计算器可编写下列计算程序。 4、程序清单 （1）、3XYF（主程序） "1→XY2→FS"?→V:V=1=>Goto 1:V=2=>Goto 2↙（选择计算功能） Lbl 1:File 1:”XLn”?→S:Prog “QD”↙（选择线路）
关于坐标正反算的应用
关于坐标正反算的应用
关于坐标正反算应用备注 1、图纸上如果单纯只注明曲线的半径，就说明这段曲线是缓和曲线； 2、如果有注明曲线的所有要素，就说明这段曲线是圆曲线； 3、如果在图纸的下方有注明曲线的各要素，而且将切线分成几段，就说明这段曲线是缓和曲线带圆曲线； 4、在同一条直线同一方向上任何点的方位角都是相同的。 5、在计算方位角时，两个坐标输入次序先后不同时，得出的方位角不同，但反算距离是一样的。
○1关于坐标正反算的应用 （先点击解析交会和工具；曲线的转角=转向角也是偏角） 一、以知两点坐标，求距离方位角？称为反算 例： 测站点坐标待定点坐标 1、点击坐标正反算→点击坐标反算→输入起点坐标X（1234.5678），Y （8765.4321）→再输入终点坐标X(1293.7422),Y(8870.2181)。 2、点击计算，得出反算方位角60.324509(即至待定点方位角60o32′45.09")。得出反算距离120.339999(即至待定点距离)。 下表是按上面算式计算的结果数据
二、已知一个点坐标，至待定点距离(120.339999)，坐标点至待定点方位角 (60.324509,即60o32′45.09")。求待定点坐标？称坐标正算 例： 已知坐标点待定点坐标 1、点击坐标正反算→点击坐标正算→输入起点坐标X（1234.5678）→ Y(8765.4321)。 2、再输入已知方位角（60.324509）,输入已知距离(120.339999)。 3、点击计算，得出待定点坐标结果：X=1293.742201 Y=8870.218099
圆曲线、缓和曲线反推里程和偏距公式
中铁八局三公司-陈超 ——若有疑问请至Email:ceo@https://www.wendangku.net/doc/2714171227.html, 圆曲线、缓和曲线反推里程和偏距计算 已知数据：交点坐标JD(X)=3378226.731,交点坐标JD(Y)=456053.721,交点里程JD(DK)=DK10+021.359,圆曲线半径R=4500,缓和曲线长LS1=360, LS2=360,转角αy=53°12’46.1”,第一直线方位角:98°56’55.62” 方法：由线路外任意一点反推该点里程和偏距通称为坐标反算或者投影归算。 曲线上任意一点法线垂距计算公式：90)-sin()-(-)90-cos()(i i i αα**-=''x x y y d p i p P ’点距中桩距离公式：)90-sin(/)(p p p p y y D α'-= 例题：假设以DK10+500左5m 作为测点未知点位,其坐标为：xp ’=3377385.2777,yp ’=456241.6099,将通过这个坐标反推该点里程和偏距。 解： 1、假定这个坐标的近似里程为：DK10+000,计算坐标其切线方位角为x i =3377706.6675, y i =455858.5250,αi =127°23′2.26” 则p ’距DK10+000法线垂距为： ∵90) -sin()-(-)90-cos()(i i i αα**-=''x x y y d p i p ∴d1=(456241.6099-455858.525)*cos(127°23’2.26”-90)-(3377385.2777-3377706.6675) *sin(127°23’2.26”-90) =499.526359m 2、根据1计算的垂距d1,则d1+近似里程DK10+000得DK10+499.526359,计算坐标其切线方位角为xi=3377381.994, yi=456237.8102,αi=133°44’38.86” 则p ’距DK10+499.526359法线垂距为： ∵90)-sin()-(-)90-cos()(i i i αα**-=''x x y y d p i p ∴d2=(456241.6099-456237.8102)*cos(133°44’38.86”-90)-(3377385.2777-3377381.994) *sin(133°44’38.86”-90) =0.474557m 3、根据2计算的垂距d2,则d2+ DK10+499.526359得DK10+500.00092,计算坐标其切线方位角为xi=3377381.665, yi=456238.153,αi=133°45’00.61” 则p ’距DK10+500.00092法线垂距为： ∵90)-sin()-(-)90-cos()(i i i αα**-=''x x y y d p i p ∴d3=(456241.6099-456238.153)*cos(133°45’00.61”-90)-(3377385.2777-3377381.665) *sin(133°45’00.61”-90) =-0.001104m d3由于小于垂距限差3mm,故认为DK10+500.00092为p ’点对应的中桩桩号,计算坐标其切线方位角为xi=3377381.6653, yi=456238.1530,αi=133°45’00.61” 则p ’距中桩距离为： ∵) 90-sin(/)(p p p p y y D α'-= ∴(456238.1530-456241.6099)/sin(133°45’00.61”-90) =-4.999022m 因此,p ’点位置所对应在曲线上的桩号为DK10+500.00092左侧4.999022米。
坐标计算方法
旋转坐标系法求缓和曲线坐标 1、旋转坐标系原理 1.1旋转公式 1cos 1sin 1sin 1cos x x y y x y αααα =-=+ 对于测量坐标系逆时针旋转为α取正值，顺时针为负。例如：原坐标系中的()1,1点，坐标系旋转45 °后，在目标坐标系为(。 1cos 451sin 4501sin 451cos 45x y =*?-*?==*?+*?=
2、利用旋转坐标计算缓和曲线任意点的坐标原理 利用缓和曲线坐标公式求 5913 48 16 3711 2610 14034565990401633642240l l l x l A A A l l l y A A A =-+-=-+ 然后旋转坐标轴，γ为方位角，把原坐标系逆时针旋转方位角。 1cos 1sin 1sin 1cos x x y y x y γγγγ =-=+ 3、用旋转坐标系法求曲线坐标 已知： ①缓和曲线上任一点离ZH 点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：0l ④过ZH 点的切线方位角：γ ⑤转向角系数：K (1或－1)左转为-1右转为1 计算过程： 3.1、求直缓点ZH 的坐标 3.1.1缓和曲线要素
A =2 03 00 2242240()tan 2 l p R l l m R T R p q α = =- =++ 00cos sin z z x x T y y T γγ =-=- 3.1.2求第一缓和曲线上任意点在原坐标系中的坐标 5913 4816 3711 2610 14034565990401() 633642240l l l x l A A A l l l y K A A A =-+- =-+ 左转为K=-1右转为K=1，因为右转时y1为正，左转时y1为负 3.1.3旋转坐标系 1cos 1sin 1sin 1cos z z x x x y y y x y γγγγ =+-=++ 3.2、求圆曲线上任意点的坐标 3.2.1求圆曲线上任意点在原坐标系上的坐标
fx-4800P缓和曲线和圆曲线坐标正反算程序(正确版)
缓和曲线和圆曲线坐标正反算程序主程序“TYQXJS” Lb1 0↙→（EXE） {NUVOGHPRQ}:“1.SZ=>XY”: “2.xy=>SZ”: N:U“QDX”:V“QDY”:O“QDLC”:G“FWJ”: H“LS”:P“RO”:R“RN”:Q“ZP=-1,YP=+1,ZZ=0” :C=1÷P:D=(P-R)÷(2HPR):E=180÷π:N=1=>Goto 1:≠=>Goto 2◣↙Lb1 1:{SZ}:SZ:W= Abs(S-O): Prog“1”: X“XS”=X◢ Y“YS”=Y◢ F“FS”=F-90◢ Goto 3↙ Lb1 2:{XY}:XY:I=X:J=Y: Prog“2”:S“S”=O+W◢ Z“Z”=Z◢ Goto 3↙ Lb1 3↙ {DE}:E“QX-JJ,Z-1,Y+1”:D“BZ-JL”↙ F=F+E↙ X=X+D Cos F◢ Y=Y+D Sin F◢ {DE}:D“BZ-JJ”:E“JJ,Z-1,Y+1”↙ F=F+E↙
X=X+D Cos F◢ Y=Y+D Sin F◢ Goto 0 子程序1：“1” A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1 -L:M=1-K:X=U+W(ACos(G+QEKW(C+KWD))+BCos(G+QELW(C+LWD))+BCos(G+Q EFW(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+QE W(C+WD)+90:X=X+Z Cos F:Y=Y+Z Sin F 子程序2: “2” T=G-90:W=Abs((Y-V)Cos T-(X-U)Sin T):Z=0: Lb1 0:Prog“1”:L=T+QEW(C+WD):Z=(J-Y)Cos L-(I-X)Sin L:Abs Z<1E-6=>Goto 1: ≠=>W=W+Z:Goto 0◣↙(E为：4800P键盘的EXE键) Lb1 1:Z=0Prog“1”:Z=(J-Y)÷Sin F◣ 注：第一缓和曲线起点半径输入无穷大（10 45），终点输入圆曲线半径；第二缓和曲线起点半径输入圆曲线半径，终点半径输入无穷大（10 45）；圆曲线输入给出的起点和终点半径；直线段则都输入无穷大（10 45）。 其中起点切线方位角用此程序计算。 ◣(代替空心) 坐标正算 程序“ZBZS” Lb1 0↙ X“HSX”:Y“HSY”:U“CZX”:V“CZY”↙ X-U≥0=>Goto 1: ≠=>Goto 2◣↙
通过逐桩坐标计算曲线要素
通过逐桩坐标表推算曲线要素（CAD篇） 摘要：现在从事工程行业的都流行使用AutoCAD进行绘制图形，为了更好的利用这个绘图工具来绘制线路曲线要素，本文将讲解如何通过设计院提供的逐桩坐标表推算未知曲线要素。 关键词：AutoCAD 技巧曲线要素 说明：AutoCAD已经成为国际上广为流行的绘图工具。具有良好的用户界面，通过交互 菜单或命令行方式便可以进行各种操作。它的多文档设计环境，让非计算机专业人员也能很快地学会使用。在不断实践的过程中更好地掌握它的各种应用和开发技巧，从而不断提高工作效率。 如何提高CAD速率？ 通常在开始绘图的时候一些人由于对工具命令不熟悉直接使用工具栏等查找命令，这样对制图的效率会大打折扣从而导致绘图的速率缓慢，提高制图的方法需要掌握CAD的快捷命令，孰能生巧的记住，然后择优选用其中的一些常用的绘图命令，把繁琐的长命令转化为简单的命令使用，其次需要多练习绘图的方式与方法才会提高绘图水平。 推算原理： 通过逐桩坐标表（含曲线五大桩）然后利用Excel生成展点命令在AutoCAD中进行坐标展点，再通过工具或命令绘制进行查询曲线长、切线长、外失距、交点坐标、交点里程、曲线半径、方位角、转角等。 准备工作： 1、逐桩坐标表X、Y（含曲线五大桩） 2、AutoCAD绘图软件 演示版本为：AutoCAD 2007 示例文件：某高速铁路逐桩坐标表 演示范围：DK07+586.707~DK12+126.03（由于该交点属于大转角则演示明显）
操作流程：坐标展点→绘制半径→绘制切线长→查询方位角→查询转角→查询交点坐标→查询交点里程→查询外失距→绘制缓和曲线。（请注意逐桩坐标表中所提供的ZH、HY、QZ、YH、HZ等说明） 准备操作如下： 1、打开“逐桩坐标表”并复制（里程桩号、坐标X、坐标Y）数据到“曲线坐标计算程序VBA 4.6”的“交点法正算”表格中，效果图如下： 逐桩坐标表见（本文附件）下载地址附后！
圆曲线缓和曲线计算公式
圆曲线缓和曲线计算公式
圆曲线缓和曲线计算公式 2011-09-13 15:19:36| 分类：默认分类|字号订阅 第九章道路工程测量(圆曲线缓和曲线计算公式) 学习园地2010-07-29 13:10:53阅读706评论0 字号：大中小订阅 [教程]第九章道路工程测量(圆曲线缓和曲线计算公式)未知2009-12-09 19:04:30 广州交通技术学院第九章道路工程测量(road engineering survey) 内容：理解线路勘测设计阶段的主要测量工作（初测控制测量、带状地形图测绘、中线测设和纵横断面测量）；掌握路线交点、转点、转角、里程桩的概念和测设方法；掌握圆曲线的要素计算和主点测设方法；掌握圆曲线的切线支距法和偏角法的计算公式和测设方法；了解虚交的概念和处理方法；掌握缓和曲线的要素计算和主点测设方法；理解缓和曲线的切线支距法和偏角法的
计算公式和测设方法；掌握路线纵断面的基平、中平测量和横断面测量方；了解全站仪中线测设和断面测量方法。 重点：圆曲线、缓和曲线的要素计算和主点测设方法；切线支距法和偏角法的计算公式和测设方法；路线纵断面的基平、中平测量和横断面测量方法 难点：缓和曲线的要素计算和主点测设方法；缓和曲线的切线支距法和偏角法的计算公式和测设方法。 § 9.1 交点转点转角及里程桩的测设一、道路工程测量概述 分为：路线勘测设计测量(route reconnaissance and design survey) 和道路施工测量(road construction survey) 。（一）勘测设计测量(route reconnaissance and design survey) 分为：初测(preliminary survey) 和定测(location survey) 1、初测内容：控制测量(control survey) 、测带状地形图(topographical map of a zone) 和纵断面图(profile) 、收集沿线地质水文资
坐标计算方法
已知：①缓和曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当计算第二缓和曲线上的点坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与计算第一缓和曲线时相反 xZ，yZ为点HZ的坐标 切线角计算公式：
已知：①圆曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当只知道HZ点的坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与知道ZH点坐标时相反
xZ，yZ为点HZ的坐标 三、曲线要素计算公式 公式中各符号说明： l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）l1——第一缓和曲线长度 l2——第二缓和曲线长度 l0——对应的缓和曲线长度 R——圆曲线半径 R1——曲线起点处的半径 R2——曲线终点处的半径
P1——曲线起点处的曲率 P2——曲线终点处的曲率 α——曲线转角值 四、竖曲线上高程计算 已知：①第一坡度：i1(上坡为“＋”，下坡为“－”) ②第二坡度：i2(上坡为“＋”，下坡为“－”) ③变坡点桩号：SZ ④变坡点高程：HZ ⑤竖曲线的切线长度：T ⑥待求点桩号：S 计算过程： 五、超高缓和过渡段的横坡计算
卡西欧5800公路坐标正反算程序文件
0年11月 14日 目 录 一、坐标正算基本公式...............................................................02 二、坐标反算原理...............................................................04 三、高程数据库录入变换.........................................................05 四、计算器程序...............................................................07 01、ZBZS(坐标正算)...............................................................07 02、ZBFS （坐标反算）...............................................................08 03、GCJF(高程积分)...............................................................09 04、PJFY （坡脚放样）...............................................................10 05、JFCX （积分程序）...............................................................11 06、ZBFY （坐标放样）...............................................................11 07、DT （递推）...............................................................12 08、HP （横坡）...............................................................13 09、LK （路宽）...............................................................14 10、SJK1（平面数据库）......................................................14 11、SJK2（纵面数据库）......................................................14 12、SJK3（左路宽度数据库）......................................................15 13、SJK4（右路宽度数据库）......................................................15 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库）......................................................16 16、SJK7（左上边坡数据库）......................................................17 17、SJK8（右上边坡数据库）......................................................18 五、后记 (19)
曲线坐标计算
曲线坐标计算 一、 圆曲线 圆曲线要素：α---------------曲线转向角 R---------------曲线半径 根据α及R 可以求出以下要素： T----------------切线长 L----------------曲线长 E----------------外矢距 q----------------切曲差（两切线长与曲线全长之差） 各要素的计算公式为： ??=180π αR L (弧长) )12(sec -=αR E （sec α=cos α的倒数） 圆曲线主点里程：ZY=J D －T QZ=ZY ＋L ／2 或 QZ=JD －q /2 YZ=QZ ＋L ／2 或 YZ=JD ＋T －q JD=QZ ＋q ／2（校核用） 1、基本知识 ◆ 里程：由线路起点算起，沿线路中线到该中线桩的距离。 ◆ 表示方法：DK26＋284.56。 “+”号前为公里数，即26km ，“+”后为米数，即284.56m 。
CK ——表示初测导线的里程。 DK ——表示定测中线的里程。 Ｋ——表示竣工后的连续里程。 铁路和公路计算方法略有不同。 2、曲线点坐标计算（偏角法或弦切角法） 已知条件：起点、终点及各交点的坐标。 1）计算ZY、YZ点坐标 通用公式： 2）计算曲线点坐标 ①计算坐标方位角 i 点为曲线上任意一点。 li 为i 点与ZY点里程之差。 弧长所对的圆心角 弦切角 弦的方位角 当曲线左转时用“-”，右转时用“+”。 ②计算弦长
③计算曲线点坐标 此时的已知数据为： ZY（x ZY，y ZY）、 ZY- i、C。 根据坐标正算原理： 切线支距法这种方法是以曲线起点ZY或终点YZ为坐标原点，以切线为X轴，以过原点的半径为Y轴，则圆曲线上任意一点的切线支距坐标可通过以下公式求得： 利用坐标平移和旋转，该点在大地平面直角坐标系中的坐标可由以下公式求得： 式中：α为ZY(YZ)点沿线路前进方向的切线方位角。当起点为ZY 时，“±”取“＋”，X0=X(ZY), Y0=Y(ZY), 曲线为左偏时应以y i=-y i代入；当起点为YZ时，“±”取“-”，X0=X(YZ), Y0=Y(YZ), 曲线为左偏时应以y i=-y i代入； 注：1、同弧所对的圆周角等于圆心角的一半 2、切线性质圆的切线与过切点的半径相垂直 3、弦切角定理弦切角等于它所夹弧上的圆周角 4、弧长公式 由L/πR=n°/180°得L=n°πR/ 180°=nπR/180 二、缓和曲线（回旋线） 缓和曲线主要有以下几类： A：对称完整缓和曲线（基本形）------切线长、ls1与ls2都相等。B: 非对称完整缓和曲线---------------切线长、ls1与ls2都不相等
缓和曲线上任意点坐标计算程序
第一缓和曲线加圆曲线上任意点坐标计算程序：L1：U=U"Ｘ0"：V=V"Y0"：F"FANG"=F：E=E"LEFT-1"：LbI 0 L2：｛B｝：｛D｝：｛P｝ L3：L=AbS（B-A"ZHD"） L 4： L5：X=L-LX Y5/（40R2S2） L6：Y= LX Y3/（6RS）- LX Y7/（336RX Y3SX Y3）：G=90L2/（∏RS） L GOtO 2 7： L8：LbI 1 L9：L=L-S L10：O=90S/（∏R）+90L/（∏R） L11：M=2（Rsin（90L/∏/R）） L12：X=S-SX Y3/（40R2）+Mcos O L13：Y=S2/（6R）+MsinO:G=90S/（∏R）+180L/（∏R） L GOtO 2 14： L15：LbI 2 L16：W=tan-1（Y/X）:Q=√（X2+Y2） L 17： L E=1=>G=-G 18： L19：X[1]=U+Qcos（F+W）+Dcos（F+G+P）◢ L20：Y[1]=V+Qsin（W+F）+Dsin（F+G+P）◢ L21：GOtO 0
注、○1、XO—为起点X坐标 EXE ○2、YO—为起点Y坐标 EXE ○3、F？—方位角 EXE ○4、LEFT-1？—左偏取1右偏取0 EXE ○5、B？—所求坐标点里程（起点输0时为到起点长度）EXE ○6、ZHD？—为直缓点里程或直圆点里程（起点可以输0）EXE ○7、S？—缓和曲线长、圆曲线时输为0 EXE ○8、R？—半径EXE ○9、D？—中桩到边桩长度EXE ○10、P？—左右方向与中线切线交角、法线方向时左-90右+90 EXE ○11、上述每一步输完后必须确认、结果显示字后转到B进行循环操作。
fx-5800p全线坐标正反算带高程计算程序(线元法)
曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序在CASIO fx-5800P计算器运行。 2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXJS)(A) Deg:fix 3 119→DimZ “INPUT(0) Or DATA(Else)”?I Lbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?N If N=1 Or N=5:Then Goto 1 Else If N=2 Or N=3 Or N=4:Then Goto 2 Else Goto 3 IfEnd:IfEnd Lbl 1:“K(m)=”?S If S<0:Then Goto 0:IfEnd “JL(m)=”?Z If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd If I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEnd S-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEnd Prog “SUB-ZS”:Prog “SUB-GC” If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F?DMS◢ IfEnd:IfEnd