曲线任意里程中边桩坐标正反算及放样fx-4850程序.doc

曲线任意里程中边桩坐标正反算及放样fx-4850程序（第五次修改）作者：李正云 本程序是在yshf及hangniu1973两位师傅的成果上作的一次改动，修改内容主要有一下几点：1、在变量符号上修改，目的是更能与工地实际结合；2、把原主程序分为两个，更能灵活调用；3、加入放样程序，做到坐标计算到放样一体化；4、使得整个测量放样过程更加简单。GAUSSLE坐标正反算fx－4850程序源程序1.正算主程序??GSZSIX0：SY0：OK0：GF0：HKN：PR0：RRN：Q”Q(-Z??+Y)” :D=(P-R)÷(2(H-O)PR)：KL”L(-Z??+Y)” :M”ANG(YJJ)”=90:（注：此处若不给M赋值，则可计算斜交点）J=Abs(K-O)：ProgSUB1：”FWJ=”:F=F-M:”X=”:X=X:Pause0: ”Y=”:Y=Y◢2. 反算主程序??GSFSXY：Z[2]=X：Z[3] =Y：IX0：SY0：OK0：GF0：HKN：PR0：RRN：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÷π)：AbsL1E-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+LsinF4. 曲线元要素数据库：DAT-MK≥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-ZBProg”DAT-M”:Prog”GSZS”6、M线（坐标计算－放样）组合程序??MG-FYProg”MG-ZB”:Prog”LTKZD”: Prog”FY”7、M线（坐标反算）组合程序??M-GSFBProg”DAT-M”:Prog”GSFS”说明：一、程序功能及原理1.功能说明：本程序由两个主程序——正算主程序(GSZS)、反算主程序(GSFS)和两个子程——正算子程序(SUB1)、线元数据库(DAT-M)构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本程