CASIO fx-4500 坐标正算程序

用CASIO fx—4500计算器
计算道路平曲线中桩及其边桩的大地坐标

一：已知条件：
　　缓和曲线长：L1 ；平曲线半径：R ；导线边方位角：A ；交点里程桩号：J ；边桩横距：S ；直缓点、缓直点的大地坐标：X0、Y0。
二：计算程序：
ZBJS (程序名)
L1 Lbi 0
L2 E“L1”：R：A：B：J：S
L3 T=Rtan（Abs B/2）+E2/（24R）×tan　
（Abs B/2）+E/2-E xY3/（240R2）▲
L4 L=Abs（RBπ/180）+Ｅ▲
L5 W=28.6479Ｅ/Ｒ：B<0 => W= -WΔ
L6 Z“ZH”=J-T▲H“HY”=Z+Ｅ▲
G“QZ”=Z+L/2▲D“YH”=Z+L-Ｅ▲
I“HZ”=Z+L▲
L7 Lbi 1:{N}:NGoto　0△
N>I=>Goto　0△
L8 N≤H =>F=N-Z:GotO 2△
L9 N≤D =>F=N-H:Goto 3△
L10 N >D =>F=I-N:Goto 2△
L11 Lbi 2:K=W×(F/E)2/3:C=F：Goto 4　
L12 Lbi 3:K=90F/(πR)):
　　C＝2Rsin K:B<0 => K= -K△Goto 5
L13 Lbi 4:N>D = Goto 6△ U“X0”：V“Y0”
L14 X=U+C cos（A+K）▲Y=V+C sin（A+K）▲
L15 X[3]“X1”=X+S Cos（A+3K+90）▲
Y[3]“Y1”=Y+S Sin（A+3K+90）▲
L16 X[5]“X2”=X+S Cos（A+3K-90）▲
Y[5]“Y2”=Y+S Sin（A+3K-90）▲
L17 X[7]=X：Y[7]=Y：Goto 1
L18 Lbi 5：U=X[7]：V=Y[7]
L19 X = U+C Cos（A+W+K）▲
Y = V+C Sin（A+W+K）▲
L20 X[3]“X1”=X+S Cos（A+W+2K+90）▲
Y[3]“Y1”=Y+S Sin（A+W+2K+90）▲
L21 X[5]“X2”=X+S Cos（A+W+2K-90）▲
Y[5]“Y2”=Y+S Sin（A+W+2K-90）▲
Goto 1
L22 Lbi 6：P“X0”：Q“Y0”
L23 X= P+C Cos（A+B+180-K）▲
Y= Q+C Sin（A+B+180-K）▲
L24 X[3]“X1”=X+S Cos（A+B+90-3K）▲
Y[3]“Y1”=Y+S Sin（A+B+90-3K）▲
L25 X[5]“X2”=X+S Cos（A+B-90-3K）▲
Y[5]“Y2”=X+S Cos（A+B-90-3K）▲
Goto 1
三、符号标识：
L1：缓和曲线；R：平曲线半径；A：导线边方位角；B：转角；J：交点里程；S：边桩横距；T：切线长；L：曲线长；ZH：直缓点里程；HY：缓圆点里程；QZ：曲中点里程；YH：圆缓点里程；HZ：缓直点里程；N：曲线上任一点里程桩号；X0、Y0：直缓点、缓
直点大地坐标；X、Y曲线上任一点大地坐标；
X1、Y1、X2、Y2：对应点边桩大地坐标。
四、算例：
某弯道平曲线缓和曲线长200m，半径1800m，起始导线边方位角87о09′06″，左转角15ο26′20″，交点里程102K+044.126，直缓点大地坐标：（38310.860,40822.807）,缓直点大地坐标(38435.931,41493.211)。边桩横距30m。
运行本程序，根据计算器提示逐一输入数据，可显示：
1:曲线元素：切线长T=344.107m;
曲线长L=685.027m
2:曲线主点里程:ZH=101+700.019;
HY=101+900.019;QZ=102+042.533;
YH=102+185.046;HZ=102+385.046