曲线任意里程中边桩坐标正反算及放样fx-4850程序（第五次修改）

概要：算斜交点）J=Abs(K-O)：Prog"SUB1"：”FWJ=”:F=F-M:”X=”:X=X:Pause0: ”Y=”:Y=Y◢2. 反算主程序 GSFSXY：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

曲线任意里程中边桩坐标正反算及放样fx-4850程序（第五次修改）,http://www.gczl8.com
曲线任意里程中边桩坐标正反算及放样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点通用公式计算线路中边桩坐标并计算放样数据。
    利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。

二、使用说明
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值读取数据库DAT-M，计算器自动将里程K所在线元数据赋给反算主程序GSFS进行试算，试算出的里程和边距须带入正算主程序GSZS中计算坐标，若坐标吻合则反算正确。
2、输入与显示说明
   （1）输入部分：
X0 ？线元起点的X坐标
     Y0 ？线元起点的Y坐标

[1] [2] [3]  下一页