2024年5月8日发(作者：甘文心)
曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序（修
改版）
一、
程序功能及原理
1.功能说明：本程序由一个主程序(XYZ)和两个子程——正算子程序(A)、反
算子程序(B)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或
非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度
、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任
意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个
子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反
算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运
行。
2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计
算放样数据。
二、
源程序
1、
主程序(XY-ZD)
Defm 2:F=1：(F=1正算，F=2反算)Z[1]=90（右边线与路中线右夹角）
Prog”THB”:Fix 3:”X=”:X▲”Y=”:Y▲（F=2时，为“KM=“：Z▲”D=“：D
▲）
2次程序：THB
Z“KM“:A=第一段线终点桩号：B=第二段线终点桩号：C=第三段线终点桩号：E=
第四段线终点桩号：F=第五段线终点桩号：……（有多少段就加多少段;超过十五
段，要另加子程序转过来）
第 1 页
Z
Z
三角）
C=1÷P：S=(P-R)÷(2HPR)：E=180÷π：F=1=>Goto 1：≠>Goto 2Δ←
Lbl 1：D：W=Z-O：Prog "A"：Goto 3
Lbl 2：{XY}：XY：I=X：J=Y：Prog "B"：Z=O+W：
D=D：Goto 3
Lbl 3：F=1=>X:Y :≠>Z:D:
3. 正算子程序(A)
A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.:L=0.230765
3449:M=0.5:X=U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C
+MWS))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))): Y=V+W
(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE
(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS)))：Z[2]=G+QEW(C+WS)+Z[1]
：X=X+Dcos Z[2]：Y=Y+Dsin Z[2]
4. 反算子程序(B)
T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：D=0：Lbl 0：Prog "A"：L=T+QEW(C+WS)
：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘
Lbl 1：Z=0：Prog "A"：Z=(J-Y)÷sinZ[2]
5.曲线元要素数据库：C1~C5…….
第 2 页
2024年5月8日发(作者：甘文心)
曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序（修
改版）
一、
程序功能及原理
1.功能说明：本程序由一个主程序(XYZ)和两个子程——正算子程序(A)、反
算子程序(B)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或
非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度
、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任
意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个
子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反
算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运
行。
2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计
算放样数据。
二、
源程序
1、
主程序(XY-ZD)
Defm 2:F=1：(F=1正算，F=2反算)Z[1]=90（右边线与路中线右夹角）
Prog”THB”:Fix 3:”X=”:X▲”Y=”:Y▲（F=2时，为“KM=“：Z▲”D=“：D
▲）
2次程序：THB
Z“KM“:A=第一段线终点桩号：B=第二段线终点桩号：C=第三段线终点桩号：E=
第四段线终点桩号：F=第五段线终点桩号：……（有多少段就加多少段;超过十五
段，要另加子程序转过来）
第 1 页
Z
Z
三角）
C=1÷P：S=(P-R)÷(2HPR)：E=180÷π：F=1=>Goto 1：≠>Goto 2Δ←
Lbl 1：D：W=Z-O：Prog "A"：Goto 3
Lbl 2：{XY}：XY：I=X：J=Y：Prog "B"：Z=O+W：
D=D：Goto 3
Lbl 3：F=1=>X:Y :≠>Z:D:
3. 正算子程序(A)
A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.:L=0.230765
3449:M=0.5:X=U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C
+MWS))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))): Y=V+W
(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE
(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS)))：Z[2]=G+QEW(C+WS)+Z[1]
：X=X+Dcos Z[2]：Y=Y+Dsin Z[2]
4. 反算子程序(B)
T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：D=0：Lbl 0：Prog "A"：L=T+QEW(C+WS)
：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘
Lbl 1：Z=0：Prog "A"：Z=(J-Y)÷sinZ[2]
5.曲线元要素数据库：C1~C5…….
第 2 页