CASIO-4800完整缓和曲线任意中边桩坐标计算(绝对正确实用)
DEFM 4
Lbl 0:{Z,S,O,U}
A:B:C”JDX”:D”JDY”:E”JD”:R:Y”LS”:T”LC”:L:J
Lbl 2:Z【1】=E-T: Z【2】= Z【1】+Y: Z【3】= Z【2】+L: Z【4】= Z【3】+Y:Z≤Z【1】=>GOTO 3: ≠=>Z＞Z【1】=>GOTO 4
Lbl3:Q”X(Z1)”=C+T*COS(A+180) ▲V”Y(Z1)” =D+T*SIN(A+180) ▲K=A: GOTO D
Lbl 4: Z≤Z【2】=>GOTO 5: ≠=>Z＞Z【2】=>GOTO 7
Lbl 5: F=(Z- Z【1】)-( Z- Z【1】)^5/(40R^2Y^2)
I=(Z- Z【1】)^3/(6RY)- (Z- Z【1】)^7/(336R^3Y^3):M=A+30*J*( Z- Z【1】)^2/(∏RY)
Lbl 6: Q”X(H1)”=C+T*COS(A+180)+ ABS(F^2+I^2)^1/2*COS M▲V”Y(H1)” =D+T*SIN(A+180) + ABS(F^2+I^2)^1/2*SIN M▲K=A+90*J( Z- Z【1】)^2/(∏RY): GOTO D
Lbl 7: Z≤Z【3】=>GOTO 8: ≠=>Z＞Z【3】=>GOTO A
Lbl 8: M=A+90*J*Y/(∏R)+90*J*( Z- Z【2】) /(∏R): F=R SIN(180*( Z- Z【2】) /(∏R)):I=R(1-COS(180*( Z- Z【2】) /(∏R))):N= A+30*J*Y/(∏R):P= ABS((Y-Y^3/(40R^2)) ^2+(Y^2/(6R)) ^2) ^1/2
Lbl 9: Q”X(C)”=C+T*COS(A+180)+ ABS(F^2+I^2)^1/2*COS M+P*COS N▲V”Y(H1)” =D+T*SIN(A+180) + ABS(F^2+I^2)^1/2*SIN M+P*SIN N▲K=A+180*J*(Y/2+Z- Z【2】) /(∏R): GOTO D
Lbl A: Z﹤Z【4】=>GOTO B: ≠=>Z≥Z【4】=> Q”X(Z2)”=C+T*COS B▲V”Y(Z2)” =D+T*SIN B▲K=B: GOTO D
Lbl B: F=(Z【4】- Z)-( Z【4】- Z)^5/(40R^2Y^2) :I=( Z【4】- Z)^3/(6RY): M=B-30*J*( Z【4】- Z)^2/(∏RY)
Lbl C: Q”X(H2)”=C+T*COS B- ABS(F^2+I^2)^1/2*COS M▲V”Y(H1)” =D+T*SIN B -ABS(F^2+I^2)^1/2*SIN M▲K=B-90*J(Z【4】- Z)^2/(∏RY)
Lbl D: W”X(L)”=S”LD”*COS(K+180+U)+Q▲H”Y(L)”=S*SIN(K+180+U)+V▲X”X(R)”=O”RD”*COS(K+U)+Q▲H”Y(L)”=O*SIN(K+U)+V▲GOTO 0

曲线任意里程中边桩坐标正反算及放样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)JD) ×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=**⊿⊿←┘
(注：如有多个曲线元要素继续添加入数据库DAT-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)、线元数据库(DAT-M)构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本程序可以在CASIO fx-4800P计算器及 CASIO fx-4850P计算器上运行。由于加入了数据库(DAT-M)，可实现坐标正反算的全线贯通。
组合程序5可实现M线的正算贯通，组合程序7可实现M线的反算贯通，组合程序6可实现坐标计算到放样一体化。
2．计算原理：
利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。
利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。

还有比较多的程序，大家下载就可以知道，觉得合适就下载吧，搂住很辛苦的把他们归拉在在一起的，值得下载

有位工程师编写了一个程序,变为圆曲线和直线正反算,非常适合于水利工程中渠道的测量放样,但不知哪个语句有问题,圆曲线右偏时,计算结果错误,左偏和直线段都是对的,请大师百忙之中能抽出点空修改一下,不胜感谢!!!
CASIO fx-4800P程序编制（里程中边桩计算坐标及反算）
程序编制
一、源程序(SZ-XY-SZ) "1.SZ => XY"："2. XY => SZ"：NU"X0"V"Y0"O"S0"G"F0"Q：Q=0=>R=1E-45：≠> R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙ Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：X"XS" = U+Rec(C,F)+ZcosE◢ Y"YS"=V+J+ZsinE◢ Goto 1←┙ Lbl 2：W=G+90Q：E=U+Rec(R,W)：F=V+J←┙ Lbl 3：{XY}：XY：D=Pol(X-E,Y-F):T=J：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+ DKcosP+ MπRP÷180◢ Z"Z"=DKsinB+MQ(R-D)◢ Goto 3←┙ 二、使用说明 1.由里程、中边桩计算坐标 启动程序后按以下提示输入数据： N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的X坐标） Y0 ？368.166 (直线或圆曲线起点的Y坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？60°54′03″(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且 右偏时Q=1) R ？ 100 (圆曲线的半径) S ？ 197.016 (中线里程) Z ？ -5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0) 显示计算结果 XS=172.7467767 YS=395.2316807 2.由坐标反算里程、中边桩 启动程序后按以下提示输入数据： N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的X坐标） Y0 ？368.166 (直线或圆曲线起点的Y坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？60°54′03″(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且 右偏时Q=1) R ？ 100 (圆曲线的半径) X ？ 172.7467767 (待求里程及中边桩的点的X坐标) Y ？ 395.2316807 (待求里程及中边桩的点的Y坐标) 显示计算结果 S=197.016 Z=-5.251 三、注意事项 1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等 曲线的计算； 2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。
在圆曲线右偏时反算计算错误,哪位高手能修改一下程序?不胜感谢!!!

好东西啊，谢谢楼主共享。

好东西啊，谢谢楼主共享。

好东西啊，谢谢楼主共享。有5800的么？谢谢了。。。。

非常感谢楼主！下载了！

请问大师有没有5800的程序呢？
这程序能否用在5800计算器上呢？

非常感谢楼主！下载了！