坐标反算c语言程序,道路坐标正反算无限个拐点计算器测量程序(新20070825)

《道路坐标正反算无限个拐点计算器测量程序(新20070825)》由会员分享，可在线阅读，更多相关《道路坐标正反算无限个拐点计算器测量程序(新20070825)(10页珍藏版)》请在人人文库网上搜索。

1、道路坐标正反算无限个拐点计算器测量程序(新20070825)道路坐标正反算4800一、程序功能及原理1. 功能说明：本程序由一个主程序(TYQXjs)和两个子程一一正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段 直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、 起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标，对该曲 线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植 到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在 CASIO fx-4500P计算器及CASIO。

2、 fx-4850P 计算器上运行。2 .计算原理：利用 Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。二、源程序1. 主程序(TYQXjs)Defm 2Lbl 1 : N=1 : N : 1.SZ = XY : 2.XY = SZ : N : Prog ” DAT1 ” :UX0: VY0 : OS0 : GF0: HLS : PR0 : RRN : Q : C=1-P: D=(P-R) p2HPR) : E=180*n: N=1=Goto 1 :工 Goto 2 T=90 : Z=0 : SZT : SZ:T” ANG : W=Abs(S-O) : Prog S。

3、UB1 : XS=:X=X 丄YS=:Y=Y 丄Goto 3 F=F+360:工 F360=F=F-360 ” DIST=” :丄” FW= :F tDMS 丄 Goto 12. 正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:M=0.5:F=1-LX=U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Ncos(G+QEMW(C+MWD)+Bcos(G+QE(1-L)W(C+(1-L)WD)+Acos(G+QE(1-K)W(C+(1-K)。

4、WD): Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Nsin(G+QEMW(C+MWD)+Bsin(G+QE (1-L)W(C+(1-L)WD)+Asin(G+QE(1-K)W(C+(1-K)WD): F=G+QEW(C+WD)+T : X=X+ZcosFY=Y+ZsinF3. 反算子程序(SUB2)T=G-90 : W=Abs(Y-V)cosT-(X-U)sinT) : Z=0 : Lbl 0 : Prog SUB1 : L=T+QEW(C+WD):Z=(J-Y)cosL-(I-X)sinL : AbsZGoto1 :工 W=W+Z : Goto。

5、 0 Lbl 1 : Z=0 : Prog SUB1 : Z=(J-Y) -sinF : Z 工 0=Z=Z -AbsZ 2 2Z=Z ? (I-X)+( J-Y)4. 曲线元要素数据库：DAT1S:S :S*=S v*=U=* ：v=* ：O=* ：G=* ：H=* p=* R=*:Q=* I5. 曲线元要素数据库：DAT2注：如有多个曲线元要素继续添加入数据库DAT1中，每一段(直线、缓和曲线、圆曲线等)线元曲线要素为：起点X坐标U、起点Y坐标V、起点里程0、起点切线方位角 G、线元长度H、起点曲率半径P、 止点曲率半径R、偏角标志Q (左偏Q=-1、右偏Q=1、直线Q=0 )。线元曲线各。

6、点的坐标、起点切线方位角均可用专门的excel程序计算，该excel程序还可以计算岀所需要只的其它坐标。附后专门的excel程序计算。1、规定(1) 以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时， Q=-1 ;当线元往右偏时，Q=1 ;当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0 ;当位于中线左铡时，Z取负值；当位于中线中线右 侧时，Z取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆 弧的半径(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径。

7、为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的 值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等 于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。(7) 曲线元要素数据库(DAT1 )可根据线型不同分为各个线元段输入到DAT1中，即分为直线段、缓和曲线、圆曲线等。2、输入与显示说明提示：计算前可手工输入测站坐标经Z1，Z2，反算前需对该反算线段内，正算一次任意一桩。

8、号，以确定反算起始数据。输入部分：1. SZ = XY2. XY = SZN?选择计算方式，输入 1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。S0 ?线元起点里程S ?正算时所求点的里程S ? EXEZ ?正算时所求点距中线的边距 (左侧取负，值右侧取正值，在中线上取零)ANG ?正算边桩时左右边桩连线与线路中线的右交角X ?反算时所求点的X坐标Y ?反算时所求点的 Y坐标X0 ?线元起点的X坐标Y0 ?线元起点的Y坐标 显示部分:XS=x xx正算时，计算得岀的所求点的X坐标YS=x xx正算时，计算得岀的所求点的Y坐标DIST=xxx正算时，放样距离FW=xxx正算时，。

9、放样方位角X=xxx反算时，测量点的Y=xxx反算时，测量点的放样方位角S=x xx反算时，计算得岀的所求点的里程Z=xxx反算时，计算得岀的所求点的边距(负为左边，正号为右边四、算例S=940 Z=0计算得 XS=19736.47687YS=28722.35642S=940 Z= 3.009计算得 XS=19733.47298YS=28722.531682、反算X=19831.418X=19827.336Y=28509.726Y=28506.838计算得 S=699.9999974 Z= -5 .00018164计算得 S=699.9996493 Z= 0.000145136X=19823.。

10、25398Y=28503.95084计算得 S=699.9999985 Z= 5.000003137X=19785.25749X=19781.15561X=19777.05373Y=28575.02270Y=28572.16358Y=28569.30446计算得 S=780.0000035 Z= -5 .000001663计算得 S=780.0000025 Z=- 0.000002979计算得 S=780.0000016 Z= 4.99999578X=19747.536X=19742.686X=19737.837Y=28654.131Y=28652.914Y=28651.697计算得 S=87。

11、0.0001137 Z= -4.99941049计算得 S=870.0003175 Z=- 0.00041814计算得 S=870.0002748 Z= 4.999808656X=19741.5912Y=28722.0580计算得 S=939.9999786 Z= -5.123024937X=19736.4769Y=28722.3564计算得 S=939.9999862 Z=- 0.000027710某匝道的由五段线元(直线 +完整缓和曲线+圆曲线+非完整缓和曲线+直线)组 成，各段线元的要素(起点里程SO、起点坐标X0 Y0、起点切线方位角 F0、线元长度LS、起点曲率半径 R0、止点曲率半。

12、径 RN、线 元左右偏标志 Q)如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.561125 16 31.00269.2561E451E450769.25619787.34028563.378125 16 31.0037.4921E45221.75-1806.74819766.56628594.574120 25 54.07112.779221.75221.75-1919.52719736.07228701.89391 17 30.6380.285221.759579.228-1999.81219744.03828781.65980 40 50.00100.00。

13、01E451E4501、正算(注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素)S=700Z=-5计算得XS=19831.41785YS=28509.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z= 5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS=19785.25749YS=28575.02270S=780Z=0计算得XS=19781.15561YS=28572.16358S=780Z= 5计算得XS=19777.05373YS=28569.30446S=870Z=-5计。

14、算得XS=19747.53609YS=28654.13091S=870Z=0计算得XS=19742.68648YS=28652.91379S=870Z= 5计算得XS=19737.83688YS=28651.69668S=940 Z=-5.123计算得 XS=19741.59118 YS=28722.05802X=19733.4730Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694所附后的专门的excel程序，可以计算各曲线要素的桩号、中桩坐标、边桩、桥涵等。用法：(1)主要在控制台一栏进行操作；(2)先输入高程、平面要素数据，可单独输入平面要素。输入数。

15、据只需输入黄色区域就可以了，然后 到控制台启动计算程序并刷新和保存。(3)计算中桩时，在控制台启动计算程序，输入起始桩号和计算间距，然后列表后，可以在逐桩坐标看到各桩号中桩坐标和该桩号的切线方位角( 注意：用此项功能可以计算得 ZH、HY、 YH、HZ点的中桩坐标和方位角)。以下是无用的参考资料曲线任意里程中边桩坐标正反算及放样fx-4850程序(第五次修改)本程序是在yshf及hangniu佃73两位师傅的成果上作的一次改动，修改内容主要有一下几点：1、在变量符号上修改，目的是更能与工地实际 结合；2、把原主程序分为两个，更能灵活调用；3、加入放样程序，做到坐标计算到放样一体化；4、使得整个。

16、测量放样过程更加简单。附件(点击下载):GAUSSLE坐标正反算fx 4850程序源程序1.正算主程序GSZSIX0:SY0 : OK0 : GF0 : HKN:PR0 : RRN : Q Q-Z +Y) ” ：D=(P-R) 2(H-0)PR):KL” L(-Z +Y) ” ：M” ANG(YJJ)” =90:(注：此处若不给M赋值，则可计算斜交点)J=Abs(K-O): ProgSUB1:” FWJ= :F=FM:” X=” ：X=X:Pause0” Y=” :Y=YJ2.反算主程序GSFSXY : Z2=X : Z3 =Y :IX0:SY0 : OK0 : GF0 : HKN:PR0 :。

17、 RRN : Q Q-Z +Y) ” :D=(P-R) 2(H-0)PR):J=Abs(Y-S)cos(G-90 (X-I)sin(G-90):L=0 : M M(YJJ) ” =90:Lbl 0: Prog SUB1:L=(Z3-Y)cos(G- 90+QJ(1 十 P+JD)X 180 円不2 -X)sin(G-90+QJ(1 十 P +JD) X 180:-%)AbsLGoto1:半 J=J+L Goto 0 -1Lbl 1 : L=0 : Prog SUB1 : L=(Z3-Y) sinF:” K=” ：K=O+J:Pause0: ” L=” 1L=L3.正算子程序(SUB1)Defm。

18、 4:A=0.1184634425:B=0.2393143352:Z4=0.28444 44444:C=0.0469100770:E=0.2307653449:Z1=0.5:X=l+J(Acos(G+QCJ(1 十 P+CJD)X 180 十 n )+Bcos( G+QEJ(1 十 P+EJD)X 180 十 n )+Z4cos(G+QZ1J(1 十 P+Z1JD) X 180 十 n )+Bcos(G3)Q(1 P+(1- E)JD) X 180 十 n )+Acos(G(HQC)J(1 P+(1-C)JD) X 180 十 n ):Y=S+J(Asin(G+QCJ(1 十 P+CJD)X 。

19、180 十 n )+Bsin( G+QEJ(1 十 P+EJD)X 180 十 n )+Z4sin(G+QZ1J( 1 十 P+Z1JD) X 180 十 n )+Bsin(G&Q(1 P+(1- E)JD) X 180 十 n )+Asin(G+Q-C)J(1 P+(1-C)JD) X 180 十 n )F=G+QJ(1 十 P+JD) X 180 十 n +MX=X+LcosF : Y=Y+LsinF4.曲线兀要素数据库：DAT-MK 0=KVH=*：S=*：O=*:G=* H=* p=* R=* Q=*K O=KVH=I=*:S=*：O=*:G=* H=* p=* R=* Q=*K O=。

20、KVH=I=*:S=*：O=*:G=* H=* p=* R=* Q=*K O=KVH=I=*:S=*:O=*:G=* H=* p=* R=* Q=*K O=KC H=I=*:S=*:O=*:G=* H=* p=* R=* Q=*K O=KC H=I=* ：s=* ：o=* ：G=* ：H=* p=* R=* Q=*(注：如有多个曲线元要素继续添加入数据库DAT-M 中)5、M线(坐标正算)组合程序 MG-ZBProg ” DAT” ：Prog ” GSZS 6、 M线(坐标计算放样)组合程序 MG-FY Prog ” MGZB :Prog ” LTKZD” : Prog ” FY”7、M线(坐。

21、标反算)组合程序 M-GSFBProg ” DAT” :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点通用公式计算线路中边桩坐标并计算放样数据。利用待求点至线元起点切线作垂线，逐次迭代趋 近原理反算里程及边距。10/10。