任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序

1.加编数据库作为主程序,计算中不必逐项输入"线元要素",提高运算速度，避免现场忙中出错
2.将原来的主程序并入数据库
3.计算直观,人性化
4.正算直接输入里程和边距,反算输入近似里程便可
5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样
6.愿收获与大家共享
7.核心计算程序摘自“yshf”
1. 正算子程序(SUB1)
A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L： M=1-K：X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LWD))+Bcos(G+57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD)))：Y=V+W(Asin(G+57.2958QKW(1/P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW(1/P+MWD)))：F=G+57.2958QW(1/P+ WD)+90：X=X+ZcosF：Y=Y+ZsinF
2. 反算子程序(SUB2)
T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+57.2958QW(1/P+ WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘ Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF
3、.增设数据库程序（ＳＪＫ主程序）
Lb1 4："1.SZ => XY"："2.XY => SZ"：｛NS｝：S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1： Goto0Δ←┘(第一线元数据要素) S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1： Goto0Δ←┘(第二线元数据要素) S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1： Goto0Δ←┘(第三线元数据要素) 。。。。。。。。。。Goto0Δ←┘(第N-1线元数据要素) 。。。。。。。。。。Goto0Δ←┘(第N线元数据要素) LB1 0：D=(P-R)÷(2HPR)：N=1=>Goto 1：≠>Goto 2Δ←┘ Lbl 1：{Z}：Z：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢ Y"YS"=Y◢ F"FS"=F-90◢ C“XC”：E“YC”：I=0：J=0：Pol（X-C，Y-E）：I“I=” ◢ J∠0＝＞J“J=”+360◢ ≠＞J“J=” ◢ ΔGoto4←┘ Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢ Z"Z"=Z◢ Goto4

三、使用说明
1、规定
(1).把所有相关的”线元要素“依次输入”SJK“
(2).运算时直接调用“SJK”运行，程序提示输入里程“S”？时，正算直接输入待求点里程，反算输入所求点“近似“里程
(3).程序中“XC、YC、”为测站坐标，“Ｉ＝”、“Ｊ＝”为放样距离和方位角
(4) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(5) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。
(6) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(7) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。
(8) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(9) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
2、输入与显示说明
输入部分：
1. SZ => XY 2. XY = > SZ N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算 里程和边距。
显示部分：
XS=××× 正算时，计算得出的所求点的X坐标 YS=××× 正算时，计算得出的所求点的Y坐标 FS=××× 正算时，所求点对应的中线点的切线方位角 S=××× 反算时，计算得出的所求点的里程 Z=××× 反算时，计算得出的所求点的边距

其中要注意AbsZ<1E-6=>Goto1中的“E-”是输入“EXP”及“（-）”键。

整理自本站论坛，站长未验证其正确性，请参考学习