9860线元法隧道测量(坐标正反计算、高程)计算程序
CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序
1. A（此为主程序）
Lbl 0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V：
If V=1：Then GOTO 1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ：
If V=3：Then GOTO 3 ：Else GOTO 0 ：IfEnd：Lbl 3：〝ZH=H〝？→H ：〝SDZF=Z〝？→Z：Prog〝ZGCZCX〝：GOTO 0：Lbl 1 ：〝ZH=L〝？→L：
If L>173000 And L<174661.96：Then GOTO 4 ：Else GOTO 0 ：IfEnd ：
Lbl 4：L→L：〝SDZF=Q〝？→Q：〝XLZJ,-Z+Y=Q〝：Q+0.125→Q：
Prog 〝ZBQXYS〝：〝JSJD=J〝：90→J▲Prog 〝ZSZB〝：
〝X=〝：X ▲ 〝Y=〝：Y▲〝FWJ=O〝：O▼DMS▲
L→H：Q-0.125→Z：Prog 〝ZGCZCX〝：GOTO 0：
Lbl 2：〝XO=M〝?→M：〝YO=R〝?→R：173300→L：
If M>3845505.273 And M<3846506.099
And R>499371.832 And R<500352.224 ：Then GOTO 5：Else GOTO 2：
IfEnd ：Lbl 5：0→Q：0→J：
Prog 〝ZBFS〝：〝LC=L〝：L ▲〝JL=Q〝：Q▲ 〝SDZJ，-Z+Y=Q〝：Q-0.125→Q ▲
L→H：Q→Z：Prog 〝ZGCZCX〝：GOTO 0
2.正算坐标ZBZS
( L-S ) / 4→H：90/π→F：HHF（1/T-1/I）/（K-S）→U：2HF/ I→D：C+4D+16 U→Ｏ：
O+J→P ：C+ D+ U→E：C+2D+4U→W：C+3D+9U→G：
A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+Qcos P→X ：
B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P→Y
3.反算坐标：ZBFS
Lbl 0：Prog 〝ZBQXYS〝：Prog 〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X) sinZ→P ：L+P→L：If Abs P≥0.001：Then GOTO 0 ：Else GOTO 1：IfEnd ：
Lbl 1： (R-Y)cosO-(M-X) sinO→Q
4. 曲线元要素数据库：ZBQXYS
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：
If L≥S And L＜K ：Th en **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：
……………………………
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ┘
(注：如有多个曲线元要素继续添加入数据库ZBQXYS中)
5 高程计算主程序ZGCZCX (后有修改说明)
Lbi1 ：〝SCGC=X〝？→X：〝R=M〝：5.98→M：〝CS=N〝：1.603→N：Prog＂GC SJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ：Then 0→K：GOTO 2：IfEnd ：If H>B-T And H If H≥B+T ：Then 0→K：D→C：GOTO 2：IfEnd ： Lbi 2 ：〝XLZG=G〝：A-（B-H）C-K^2÷2R→G▲If Z≥0： Then If X>100 Then 〝YO1XGC=J〝：G+N→J▲ 〝YGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲ 〝YKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲ 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ 〝YGCCQ,C+,Q-=J〝：X-J→J▲Else 〝YO1XGC=J〝：G+N→J▲ 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ IfEnd ： Else 〝ZO1XGC=J〝：G+N→J▲If X>100 ：Then 〝ZGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲ 〝ZKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲ 〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ 〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else 〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J ▲IfEnd ：IfEnd 6高程计算主程序子程序：GCSJK If H>起点桩号 And H≤第一个竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd： If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd： If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd ………………………(继续添加要素) 说明：第一部分坐标部分（1、2、3、4） V=1进入坐标正算 V=2进入坐标反算V=3进入单独的高程计算 当V不等于1、2、3时，则返回程序，要求再次输入V值。 变量说明： S……..起点里程A…..起点X坐标B….起点Y坐标C….起点方位角 I ….….起点半径K…..终点里程T .…终点半径L….计算点里程 J…输入左右（左-，零，右+）角度Q….中桩到边桩的距离 〝SDZJ-Z+Y〝为隧中支距，左〝-〝右〝+〝 坐标正算中输入的Z值为线路支距 坐标反算中得出的Z值也为线路支距 注意：1：计算中桩坐标J D,JL 等于零。 2：这里的方位角是指〝主点作为起点〝的方位角 例如：（ZY , YZ , ZH ,QZ, HY , YZ）每个点都是起点所以特别注意 3、规定 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， 半径＝半径*-1；当线元往右偏时，半径=半径；当线元为直线时，半径=以10的45次代替。 (2) 当所求点位于中线时，Q,J=0；当位于中线左侧时，J取负值；当位于中线右侧时，J取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 4、输入与显示说明 V=1的时候 JL ？正算时所求点距中线的边距(在中线上取零) J ？边桩时左右边桩连线与线路中线的交角线元左右偏标志 (左偏J=-1*角度，右偏J=角度，直线段J=0) 角度 X=×××正算时，计算得出的所求点的X坐标 Y=×××正算时，计算得出的所求点的Y坐标 FWJ=*** 计算得出所求点的方位角°′″ V=2的时候 XO=×××所求点的X坐标 YO=×××所求点的Y坐标 LC=**** 求点所对应的里程 JL=**** 求点到所对应的里程的距离 -5就是左边相反就是右边 想知道是否正确V=1 看看坐标是否一致 注意：验算的时候不需要改动JL ？ J ？[/center]◢ 说明：第二部分高程部分（5高程计算主程序：SQXBG 6高程计算主程序子程序：S JK） 1、A：交坡点高程 2、B：交坡点桩号 3、C：前坡度（上坡为正，下坡为负） 4、D：坡度（上坡为正，下坡为负） 5、R：竖曲线半径 6、H：待求点桩号 7、E：下一竖曲起点桩号 8、X〝SCGC〝测点实测高程 9、M〝R〝 O1点半径（此处如须计算开挖线加上二衬厚度、沉降量、其他结构层的厚度即可） 10、N〝CS〝线路中桩设计高程至O1点的距离 11、J〝YO1XGC〝 O1点（圆心点）高程 12、J〝YSBSJGC〝根据桩号及支距（正反算程序中）计算得到的测点（设计）高程J〝YGCFSKD〝反算中根据高程反算得到的支距 J〝YKDCQ,+C,-Q〝实测高程计算得支距与坐标反算得支距进行对比,得到超欠挖 J〝YGCCQ,C+,Q-〝实测高程与坐标反算得到高程进行对比,得到超欠挖 13、J〝ZO1XGC〝 O1点（圆心点）高程 14、J〝YSBSJGC〝根据桩号及支距（正反算程序中）计算得到的测点（设计）高程J〝ZGCFSKD〝反算中根据高程反算得到的支距 J〝ZKDCQ,+C,-Q〝实测高程计算得支距与坐标反算得支距进行对比,得到超欠挖 J〝ZGCCQ,C+,Q-〝实测高程与坐标反算得到高程进行对比,得到超欠挖 （本人施测的是一公路隧道的上行线，线路中心线至隧道中心线的距离为0.125米，线路中心线在隧道中心线的左侧。Z〝SDZJ-Z+Y〝为隧中支距，Z+0.125为线路支距。如线路中线在隧道中线的右侧，则为Z-0.125为线路支距。其中Z在左侧取-，右侧取+。） SQXBG 修改说明此式中X为实测的该点高程，V=1和3时，请将X输为0，当V=2时，输入为实测的高程数据。一般在测量时，将拱部点测得三维后，反算得到里程和支距， 即用实测高程计算得到该高程点的理论支距，再用理论支距与坐标反算的支距进行对比，超欠挖自明。同时用反算的支距也可计算出该点的理论高程，用实测的高程与计算的理论高程对比，自可得到拱部是高度是否超欠挖。 Lbi1 ：〝SCGC=X〝？→X：：〝R=M〝？→M：〝CS=N〝？→N：Prog＂GCSJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ：Then 0→K：GOTO 2：IfEnd ：If H>B-T And H If H≥B+T ：Then 0→K：D→C：GOTO 2：IfEnd ：Lbi 2 ： 〝XLZG=G〝：A-（B-H）C-K^2÷2R→G▲If Z≥0： 〝YO1XGC=J〝：G+N→J▲Then If X>100 Then 〝YGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲ 〝YKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲ 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ 〝YGCCQ,C+,Q-=J〝：X-J→J▲Else 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J：IfEnd ：▲ Else ：〝ZO1XGC=J〝：G+N→J▲If X>100 ：Then 〝ZGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲ 〝ZKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲ 〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ 〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else 〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J ▲IfEnd ：IfEnd 如在4850中请将X〝X=〝▲ Y〝Y=〝▲改为〝X=〝：X◢〝Y=〝：Y◢ 将L〝LC=〝▲Q〝JL=〝▲改为〝LC=〝：L◢〝JL=〝：Q◢ 将G〝XLZG〝=A-（B-H）C-K2÷2R▲改为〝XLZG=〝：G=A-（B-H）C-K2÷2R▲将J〝YO1XGC〝=G+N▲改为J〝YO1XGC=〝：G+N▲ 将J〝YKDCQ,+C,-Q〝= J- Abs (Z+0.125) ▲改为J〝YKDCQ,+C,-Q=〝：J- Abs (Z +0.125) ▲ 将J〝YSBSJGC〝=G+N+√(M^2-(Z+0.125)^2) ▲ 改为J〝YSBSJGC〝: G+N+√(M^2-(Z+0.125)^2) 将J〝YGCCQ,C+,Q-〝=X-J ▲改为J〝YGCCQ,C+,Q-〝: X-J 左侧同理。 查看楼主其它问题 →回复内容只查看楼主与本人回复按支持人数排序 沙发：lbs326 楼主10-4-4 0:29 算例0403 1、平曲线要素 交点号JD37 交点桩号SK174+133.90 交点坐标X=3845932.100 交点坐标Y=499588. 222 右偏角45°54′41″前方位角297°18′17″半径=1200 LS1=LS2=150 ZH点桩号SK173+550.394 HY点桩号SK173+700.394 YH点桩号SK174+511.96 HZ点桩号SK174+661.960 A1=A2=424.264 T1=T2=583.596 SK173+300 QD点坐标X=3845549.528 Y=500329.288 C=297°18′17″ SK173+550.394 ZH点坐标X=3845664.390 Y=500106.7933 C=297°18′17″ SK173+700.394 HY点坐标X=3845735.948 Y=499974.9917 C=300°53′08.5″ SK174+106.177 QZ点坐标X=3845998.621 Y=499668.235 C=320°15′37.47″ K174+511.960 YH点桩号X=3846348.183 Y=499465.9908 C=339°38′06.45″ ZBQXYS If L≥173300 And L＜173550.394：Then 173300→S：3845549.528→A：500329.288→B：297°18′17″→C：1E45→I：173550.394→K：1E45→T：IfEnd： If L≥173550.394 And L＜173700.394：Then 173550.394→S：3845664.390→A：50010 6.7933→B：297°18′17″→C：1E45→I：173700.394→K：1200→T ：IfEnd： If L≥173700.394 And L＜174106.177：Then 173700.394→S：3845735.948→A：49997 4.9917→B：00°53′08.55″→C：1200→I：174106.177→K：1200→T： IfEnd： If L≥174106.177 And L＜174511.960：Then 174106.177→S：3845998.621→A ：4996 68.235→B：320°15′37.47″→C：1200→I：174511.96→K：1200→T ：IfEnd： If L≥174511.96 And L＜174661.960：Then 174511.960→S：3846348.183→A：499465. 9908→B：339°38′06.45″→C：1200→I：174661.96→K：1E45→T ： IfEnd 注：1E45，意即无穷大。 2、竖曲线要素 第一变坡点桩号SK173+350 变坡点高程=1332.762 竖曲线半径=24498.033 前坡度-0.761% 后坡度-2.431% 竖曲线起点：SK173+145.441 竖曲线止点：SK173+554.559 第二变坡点桩号SK174+395.249 变坡点高程=1307.352 竖曲线半径=56000 前坡度-2.431% 后坡度-1.690% 竖曲线起点：SK174+187.949 竖曲线止点：SK174+602.909 第三变坡点的竖曲线起点桩号：SK174+900。 GCSJK If H>173.145.441 And H≤174187.949 ：Then 1332.762→A：173350→B：0.00761→C：-0.02431→D：24498.033→R：IfEnd ： If H>174187.949 And H≤174900 Then 1307.352→A：174395.429→B：-0.02431→C：-0.0169→D：56000→R：IfEnd 3、实例计算如下 其中：X”SCGC”输入任意实测高程值M”R”=5.98N”CS”=1.602 V=1 L输入173806 JL输入0.125 SDZJ,-Z+Y后显示0 X=3845794.175 Y=499886.936 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 V=1 L输入173806 JL输入-4.875 SDZJ,-Z+Y后显示-5 X=3845790.126 Y=499884.002 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1326.559 V=1 L输入173806 JL输入5.125 SDZJ,-Z+Y后显示5 X=3845798.224 Y=499889.869 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”ZGCFSKD”= Abs√(M^2-(X-（G+N）)^2) ▲J”ZKDCQ,+C,-Q”=J-Abs(Z+0.125) ▲J”ZSBSJGC”=G+N+√(M^2-(Z+0.125)^2) ▲J”ZGCCQ,C+,Q-”=X-J 均为相应值 V=2 L输入173806 X0=3845794.175 Y=499886.9355 LC=173805.9996 JL=0.125210 SDZJ,-Z+Y=0.000210 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 J”ZGCFSKD”J”ZKDCQ,+C,-Q”J”ZGCCQ,C+,Q-” 均为相应值 V=2 L输入173806 X0=3845790.126 Y=499884.002 LC=173805.9996 JL=-4.875065 SDZJ,-Z+Y=-5.000065 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”ZGCFSKD”J”ZKDCQ,+C,-Q”J”ZGCCQ,C+,Q-” 均为相应值 V=2 L输入173806 X0=3845798.224 Y=499889.869 LC=173806.0004 JL=-5.124899 SDZJ,-Z+Y=4.999899 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”ZGCFSKD”J”ZKDCQ,+C,-Q”J”ZGCCQ,C+,Q-” 均为相应值 V=3 ZH输入桩号173806 SDZF,-Z+Y =0 XLZJ输入0.125 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 V=3 ZH输入桩号173806 SDZF,-Z+Y =-5 XLZJ输入-4.875 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 V=3 ZH输入173806 SDZF,-Z+Y =5 XLZJ输入5.125 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 线元法万能曲线正反算简介 我的线元法是把线形分为直线和曲线，直线就不用说了，起止点桩号，坐标和方位角就可以算了；曲线最基本的组合：是由一段缓和曲线+一段圆曲线组成，任意复杂的曲线都可以分解成缓和曲线+圆曲线或者其中之一就可以。 分析最复杂的曲线可以看到: 一般复杂线形由Ls1 ,R1,Ls2, R2组成，相邻的Ls1+R1,一般满足A*A=Ls1*R1,这就是一个线元法单元，即使不满足也可以作为一个线元： 当Ls1= Ls2，且R1= R2时，为单曲线 当Ls1≠ Ls2，或者R1≠R2时，为复合曲线 当Ls1= Ls2=0时，线性为圆曲线， 当圆曲线长度为0时，线性为缓和曲线+缓和曲线， 当A*A≠Ls1*R1时，为卵形曲线，需要计算虚拟起点坐标 综合以上线形，本程序正反算计算全部可以处理。结合目前流行的线元法，本程序也可以，分为缓和曲线和圆曲线录入，方法是一样的，所不同的是起点要注意，复杂曲线，是两边向中间定义数据库，缓和曲线永远是ZH点或HZ点为起点。 曲线要素说明（有9个）： 1、起点桩号：（一般为ZH点或HZ点，或ZY点或YZ点，或者卵形公切点GQ） 2~3、起点坐标：（X，Y） 4、起点方位角：FWJ 114°15′24.33″写成：114.152433 5、线性特征：直线，左偏，右偏；三个选一个 6、终点桩号：如果起点为ZH点，终点一边为YH点，QZ点，HY点，都可以，一般为YH点，缓和曲线+圆曲线。如果缓和曲线Ls=0，就是YZ点；大小不一定按路线顺序，如果起点为HZ点，终点根据缓和曲线+圆曲线的特点，和上个线元对接上就可以了。 7、缓和曲线长度Ls： 8、圆曲线半径R: 9、回旋参数A: 一般满足A*A=Ls1*R1，不满足条件的是卵形曲线。 可以处理任意数量断链。 操作流程：1、先编辑线元数据，保存后推出。 2、如果有线元断链的输以下线元断链数据 3、打开线元万能曲线计算单点计算就可以了。 目前，已有一个例子文件在里面，在安装文件目录下“ \dmfx4.0\demo\左线”，有个CAD文件，里面有校核数据，可以看到本软件处理的逐桩表和要素表，可以验证软件的数据，任意数据坐标反算可以得到桩号和距中，任意输入桩号和距中可以正算得坐标。 授权版用户，可以通过运行交点文件编辑，保存后，退出；打开线元法数据编辑，浏览正在使用的主项目文件，就可以看到一个线元数据，点击这个文件确定，保存退出。就完成交点法数据转换线元法数据过程。 CASIO fx-4800P、fx-5800P型计算器用于线路施工 曲线中线点坐标的计算程序 中铁十局三建公司工程技术部 摘要：本文介绍了CASIO fx-4800P 、fx-5800P型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。关键词：曲线线路施工测量计算程序 1．概述 过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。随着电子技术进步和经济发展，测量仪器和测量方法的不断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。为了能够快速准确地为全站仪提供测设 数据，发挥全站仪快速测设的特点，提高测量工作效率，应采用可编程的计算器，编制计算程序。本文主要介绍应用CASIO fx-4800P型计算器的计算程序，供公司测量同行们参照使用。 2．计算程序 QXZBJS（文件名：曲线坐标计算fx-4800P） Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+M ◢G=RAπ/180：“S”S=G+L◢ LbiA：｛C，V｝：C“CSDLC=”：V“HXPJ=”：D=C-K：D≤L＝＞I=D-D5/（40R2L2）：U=D3/(6RL)-D7/（336R3L3):J=√（I2+U2）：Goto1：≠＞D≤G＝＞O=90（2D-L）/( Rπ)：I=RsinO+M：U=R（1-cosO）+P：J=√（I2+U2）：Goto2：≠＞D=S-（C-K）： = D-D5/（40R2L2）：Z[2]=D3/6RL-D7/（336R3L3）：I=T+（T-Z[1]）cosA-Z[2]sinA：U=（T-Z[1]）sinA+Z[2]cosA：J=√（I2+U2）：Goto3：Lbi1:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-30D2/ （RLπ）：H=F-90D2/ （RLπ）：≠＞Q=F+30D2/ （RLπ）：H=F+90D2/（ RLπ）⊿ Goto4：Lbi2:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I)：H=F-O：≠＞Q= F+tng-1(U/I)：H=F+O⊿Goto4：Lbi3:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I):H=F-(A-90(S-(C-K))2/ （RLπ）): ≠＞Q= F+tng-1(U/I):H=F+(A-90(S-(C-K))2/ （RLπ）):⊿ Goto4： Lbi4:B=90+H:H＜0＝＞H“QXFWJ”=B+360◢≠＞H≥360＝＞H“QXFWJ”=H-360 ◢≠＞H“QXFWJ”=H◢⊿Goto5： Lbi5: X“CSD:X”=JcosQ+N+VcosB◢ Y“CSD:Y”=JsinQ+E+VsinB◢ 经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。 2、线元数据可自动从数据库调用，也可手工输入。 3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。 4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。 5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。第三次及以后试算才要求输入近似里程。 6、程序代码规范简洁，便于阅读、理解。 完整程序清单: ZFS %正反算主程序 B=.1739274226:C=.5-B: Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS": ≠>U=1=>Prog"FS":≠>Goto 1 ZS %正算子程序 {K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB" FS %反算子程序 {KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2Δ M=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1 ZZ %高斯法中桩子程序(4节点) Prog"XL":M=K-L:O=(P-R)÷2PQR: D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D: I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1: S=X+M(BcosD+CcosE+CcosG+BcosH: T=Y+M(BsinD+CsinE+CsinG+BsinH WY %外移点计算子程序 Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1 WZ %位置显示子程序 "KJ":K:Pause 1:J◢ ZB %坐标显示子程序 "XY":S:Pause 1:T◢ YC %异常处理子程序 U=1=>K=L:U=2Δ U=3=>K=M:U=4Δ 曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序在CASIO fx-5800P计算器运行。 2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXJS)(A) Deg:fix 3 119→DimZ “INPUT(0) Or DATA(Else)”?I Lbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?N If N=1 Or N=5:Then Goto 1 Else If N=2 Or N=3 Or N=4:Then Goto 2 Else Goto 3 IfEnd:IfEnd Lbl 1:“K(m)=”?S If S<0:Then Goto 0:IfEnd “JL(m)=”?Z If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd If I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEnd S-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEnd Prog “SUB-ZS”:Prog “SUB-GC” If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F?DMS◢ IfEnd:IfEnd 线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.360docs.net/doc/9411297819.html,/：本论文仅供学习交流使用，本站仅作合理转载，原作者可来邮要求删除论 文。 摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。利用CASIO fx-9750GⅡ编程函数计算器强大的内存（可诸存63000个字符）和编程功能，编写各种计算程序，能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料，同时也使我们有更多的时间去挑战更富有创造性的工作。 关键词：坐标放线线元测量程序 1、前言 本程序采用Gauss-Legendre（高斯-勒让德）五节点公式作内核,计算速度（太约2秒）适中，计算精度很高。在此之前，本人曾用过以下公式作内核：①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点，⑦高斯-勒让德5节点，经过测试③计算最快，⑦代码稍长但计算速度只比③⑥稍慢，精度最高，可满足线元长小于1/2πD 的所有线形的精度要求。⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55m m,3/4πD时偏差为31.63mm，πD时偏差为968mm，偏差按半径倍数增大，如线元长大于1/2πD（1/2圆周长）时，可将其拆分二个或多个线元单位，以确计算保精度。 2、程序特点 事先将所有的平曲线交点的线元要素诸存到计算器内，测量时只输桩号、边距等程序会自动寻找各类要素，一气呵成地完成施工测量任务，中途不需人工转换各类要素数据，本程序可诸存几百条线路的要素数据，计算时可按需选择线路编号进行测量。测量时不需查阅及携带图纸，仅一台CASIO fx-9750GⅡ编程函数计算器即可。 本程序含一个主程序:3XYF，五个子程序:GL（公式内核）、QD（线路选择）、XL（线路要素判断）、GF（坐标反算）、File 1 （要素存放的串列工作簿）。可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正反算。 3、计算公式及原理 如图：BC 间为一曲线元，曲线元上任一点的曲率随至B 点的弧长作线性变化。设起点B 的曲率为KA ，终点C 的曲率为KB ，R 为曲线半径。±表示曲线元的偏向，当曲线元左偏时取负号，当曲线元右偏时取正号,直线段以1的45次方代替（即半径无穷大）。 式中：αΑ＝起始方位角l ＝p 点到B的距离lS=曲线总长αp＝p 点切线方位角 R1=R5=0.118463442528095 ,R2 = R4 = 0.239314335249683 , R3 = 0.28444444444444 V1=1-V5= 0.046910070 ,V 2= 1-V4 = 1 0.2307653449 V3= 0.5 利用上面公式及CASIO fx-9750GⅡ编程函数计算器可编写下列计算程序。 4、程序清单 （1）、3XYF（主程序） "1→XY2→FS"?→V:V=1=>Goto 1:V=2=>Goto 2↙（选择计算功能） Lbl 1:File 1:”XLn”?→S:Prog “QD”↙（选择线路） CASIO 5800计算器测量计算程序 上上月做这个东西的时候没仔细检查，有好几处输错了的地方，今天把它修改过来。 简要介绍： 1. 新版程序把线元法和交点法已经集成在一个模块中了，用户只需修改JD 程序和ZA程序中的数据部分即可,其余不需作任何的改动。 2. 因为每条路高程计算不尽相同，且比较复杂，现在可利用PC机EXCEL计算好打印成表格带到工地上使用，所以本版程序未对线路高程序进行专门的编程计算，而是利用统计计算模式中来输入桩号（第一列X）及左、右高程（第二、三列Y，Freq），这种输入数据的方式最为直观，易发现错误，也易修改，输入完毕后运行S程序对数据按桩号进行排序，在程序中通过调用GG程序来进行内插计算，SG=-1得左标高，SG=1得右标高（若SG输入0,则可进行一般的线性内插计算）。 3. 在JD程序和XY程序中，先将一个计算单元的数据置入矩阵F中（1行8列或1行9列），这样程序可读性极好。 4．相比原CASIO4850程序操作习惯，作了一点小小的改动，测站坐标存在Z[10]，N中，X坐标原存在M中容易被误操作修改，而设计标高存在M中，这样易于修改，因为CASIO5800没有IN，OUT功能，很不方便。 4. 程序利用Z[2]变量值来判断是采用交点法还是线元法模型计算，Z[2]＝0为线元法，否则为交点法。 一、PQX程序：计算中边桩坐标及近似的桩号反算，在运行模式直接调用。 ①Z[10]→S:”XO”?S:S→Z[10]:”YO”?N:Prog “AU” ②Lbl 2:?L:Prog “Z”:Prog “E”:1n→O:90→S ③Lbl 4:”JJ”?S:”YC”?O:SO=0 =>Goto 2…原来lbl 后没有标号４的。 ④O=-1 =>Goto 6 ⑤“X,Y”:R+OCos(Z+S)→X▲U+OSin(Z+S)→Y▲Prog “D”:Goto 4 ⑥Lbl 6:Z[7]→X:Z[8]→Y:”XF”?X:”YF”?Y:XY=0 =>Goto 4 ⑦X→Z[7]:Y→Z[8]:Pol(X-R,Y-U+1p):Z+S-J→J:”YC,DL,L”:ICos(J)→O▲ISin( J)→I▲L+I▲Goto 6 二、P程序：在程序中提供一个自由运算的模式。 ①Lbl 1:”TMP”?I:If I≠0:Then “RST”:I▲Goto 1:IfEnd 二、LYC程序：进行桩号反算及边坡放样，在运行模式直接调用。 ①Prog “AU” ②Lbl 1:Z[7]→X: Z[8]→Y: Z[6]→S: ”XF”?X :X→Z[7]:”YF”?Y:Y→Z[8]: ”ZF”?S: S→Z[6] ③Lbl 2:Prog “Z”:Y=U =>Y+1p→Y ④Pol(X-R,Y-U):J-Z→J:Isin(J)→O:Icos(J)→I ⑤If Abs(I)≤0.1:Then Prog “E”:”L,YC”:L+I→L▲O▲Goto 3:IfEnd ⑥If Z[9]≠0:Then Pol(Z[9]-SO,I):πJZ[9]÷180→I:IfEnd ⑦”DL”:I▲L+I→L:Goto 2 ⑧Lbl 3: Z[6]→S:If S=0:Then Goto 1:IfEnd ⑧M→Z ⑨Lbl 4:”SG”?Z:Z→M:If Abs(Z)=1: Then Prog “GG”:Y→Z:If X=1:Then CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序 悬赏分：+62 作者：刘工 2010-4-3 1. A（此为主程序） Lbl 0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V： If V=1：Then GOTO 1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ： If V=3：Then GOTO 3 ：Else GOTO 0 ：IfEnd：Lbl 3：〝ZH=H〝？→H ：〝SDZF=Z〝？→Z：Prog〝ZGCZCX〝：GOTO 0：Lbl 1 ：〝ZH=L〝？→L： If L>173000 And L<174661.96：Then GOTO 4 ：Else GOTO 0 ：IfEnd ： Lbl 4：L→L：〝SDZF=Q〝？→Q：〝XLZJ,-Z+Y=Q〝：Q+0.125→Q： Prog 〝ZBQXYS〝：〝JSJD=J〝：90→J▲Prog 〝ZSZB〝： 〝X=〝：X ▲ 〝Y=〝：Y▲〝FWJ=O〝：O▼DMS▲ L→H：Q-0.125→Z：Prog 〝ZGCZCX〝：GOTO 0： Lbl 2：〝XO=M〝?→M：〝YO=R〝?→R：173300→L： If M>3845505.273 And M<3846506.099 And R>499371.832 And R<500352.224 ：Then GOTO 5：Else GOTO 2： IfEnd ：Lbl 5：0→Q：0→J： Prog 〝ZBFS〝：〝LC=L〝：L ▲〝JL=Q〝：Q▲ 〝SDZJ，-Z+Y=Q〝：Q-0.125→Q ▲ L→H：Q→Z：Prog 〝ZGCZCX〝：GOTO 0 2.正算坐标ZBZS ( L-S ) / 4→H：90/π→F：HHF（1/T-1/I）/（K-S）→U：2HF/ I→D：C+4D+16 U→Ｏ： O+J→P ：C+ D+ U→E：C+2D+4U→W：C+3D+9U→G： A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+Qcos P→X ： B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P→Y 3.反算坐标：ZBFS Lbl 0：Prog 〝ZBQXYS〝：Prog 〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X) sinZ→P ： 3、交点法、线元法坐标计算 坐标计算是根据图纸中“直线及曲线转角一览表”提供的数据计算道路中桩坐标，然后和图纸提供的“逐桩坐标表”比对，如果一样则说明输入平曲线参数输入正确，可以计算边桩坐标和其他结构物坐标了；如果中桩坐标不一样，一般是平曲线参数输入有误，需要重新检查输入，另一种结果是图纸有错，这种情况少见，但不代表没有。“直线及曲线转角一览表”和“逐桩坐标表”见附件1、附件2。 线元法是以路线的起点坐标、方位角、起终点桩号等节点元素来计算出要求的坐标；交点法是以路线的交点要素和路线的主要要素来求得坐标。 ①交点法 交点：路线的转折点，路线改变方向是相邻两直线的延长线相交的点。用JD表示， 有些图 纸上用 IP表示。 看下图： 交 点是针对曲线的（包含圆曲线和缓和曲线），一段曲线就有一个交点。交点参数有：坐标（X,Y）、交点桩号、转角值、圆曲线半径R、缓和曲线长度。 教学提供软件（轻松测量、双心软件、测量工具）交点法曲线要素输入说明： 1、QD起点坐标： 起点坐标必须在直线段上，或填写前一交点的坐标。 2、JD交点曲线要素： （1）交点桩号 （2）交点坐标（X,Y） （3）曲线半径R 始点的话，起始里程有时候需要校正，当然，并不是每个图纸给出的起点里程都需要校正，大多数图纸的起点里程已经被设计院校正过，我们输入平曲线的时候需要验证一下。如果我们按照图纸给出的起点里程输入，发现后面的交点里程都和图纸相差一个相同的值，这就表明我们输入的起点里程需要校正。 起始点里程正常输入，第二、三个交点输入完成后，检查第二个交点的切线长和交点 里程是否和图纸一样，如果切线长正确，交点里程不正确，说明起点里程需要校正，将第二个交点的里程与正确里程的差值，应用到起点里程中，从而使第二个交点里程和后面交点的里程与图纸吻合。 注意：交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。对于非普通的三单元曲线，交点法不适用。非普通的三单曲线例如下页的JD18及JD19处的平曲线， 的输入是否正确，有的图纸给的方位角数据较少，需要每隔几个线元才能检验方位角。 公用模块： Option Explicit Public Const PI = 3.14159265358979 '已知A、B两点坐标计算方位角，JSFWJ的中文意思是计算方位角 Public Function JSFWJ(xa As Double, ya As Double, xb As Double, yb As Double) As Double '已知A、B两点坐标计算方位角函数过程Dim vx As Double, vy As Double vx = xb - xa: vy = yb - ya '如果A、B两点坐标相同，出现提示对话框 If vx = 0 And vy = 0 Then MsgBox "您选择的是同一个点！", vbOKOnly + vbExclamation, "提示信息" JSFWJ = 999999999# End If '计算方位角的值 If vx = 0 And vy > 0 Then '与y轴正半轴平行 JSFWJ = RadianToAngle(PI / 2#) ElseIf vx = 0 And vy < 0 Then '与y轴负半轴平行 JSFWJ = RadianToAngle(PI * 3# / 2#) ElseIf vy = 0 And vx > 0 Then '与x轴正半轴平行 JSFWJ = RadianToAngle(0) ElseIf vy = 0 And vx < 0 Then '与x轴负半轴平行 JSFWJ = RadianToAngle(PI) ElseIf vx > 0 And vy > 0 Then '第一象限 JSFWJ = RadianToAngle(Atn(vy / vx)) ElseIf vx < 0 And vy > 0 Then '第二象限 JSFWJ = RadianToAngle(Atn(vy / vx) + PI) ElseIf vx < 0 And vy < 0 Then '第三象限 JSFWJ = RadianToAngle(Atn(vy / vx) + PI) ElseIf vx > 0 And vy < 0 Then '第四象限 JSFWJ = RadianToAngle(Atn(vy / vx) + 2 * PI) End If End Function '已知A、B两点坐标计算距离，JSJLS的中文意思是计算距离S Public Function JSJLS(xa As Double, ya As Double, xb As Double, yb As Double) As Double Dim vx As Double, vy As Double vx = xb - xa: vy = yb - ya '如果A、B两点坐标相同，出现提示对话框 If vx = 0 And vy = 0 Then MsgBox "您选择的是同一个点！", vbOKOnly + vbExclamation, "提示信息" JSJLS = 99999999# End If '计算距离 JSJLS = Sqr(vx * vx + vy * vy) End Function '弧度化角度 Public Function RadianToAngle(ByVal alfa As Double) As Double Dim alfa1 As Double, alfa2 As Double alfa = alfa * 180# / PI CASIO fx5800p全线高程计算程序 GAOCHEN 主程序 Lbl 1 “KM=，<0,Stop”:?K：K<0=>Stop:“PY=”？L：Prog”GK” C-D→E:Abs(RE/2)→T：R(Abs(E)/E)→R If K≤B-T:Then 0→H:Else:If K≥B+T Then 0→H:D→C:Else K-B+T→H:Ifend:Ifend A-(B-K)C-H2/(2R)→G：Cls “KM=”：Locate 4，1，K：Locate 10，1，“PY=”：Locate 13，1，L：Fix 3 “H=”：Locate 4，2，G Prog “PODU”：（E-B）/（D-A）（K-A）+B→I：(F-C)/（D-A）（K-A）+C→J “HL=”：G+IL→X：Locate 4，3，X：Locate 11，3，“I=”：Locate 13，3，I*100 “HR=”：G+JL→Y：Locate 4，4，Y：Locate 11，4，“I=”：Locate 13，4，J*100◢Cls：Norm 2：“BM+HS≤0,Goto 1”？Z：Z≤0=> Goto 1：Cls (输入视线高) “KM=”：Locate 4，1，K：Locate 10，1，“PY=”：Locate 13，1，L：Fix 3 “QSM=”: Locate 6，2，Z-G (显示中桩读数) “QSL=”: Locate 6，3，Z-X (显示左桩读数) “QSR=”: Locate 6，4，Z-Y◢(显示右桩读数) Norm 2：Cls：Goto1 （后面可加已知视线高计算读数部分，不想计算读数则视线高输入0或负数如不想显示麻烦，可将Locate语句去掉） 以下两个子程序不需运行，只是两个独立的数据库赋值程序，字母重复不影响计算结果 GK 数据库子程序 If K≤第二曲线起点桩号：Then 第一曲线交点高程→A：第一曲线交点桩号→B：第一曲线前坡→C：第一曲线后坡→D：第一曲线半径→R：Return：Ifend …………….（有几个变坡点编几个If语句） PODU 计算坡度子程序 If K≤第一变（非变）坡段终点：Then 第一曲线起点桩号→A：第一曲线起点左坡→B：第一曲线起点右坡→C：第一曲终点桩号→D：第一曲终点左坡→E：第一曲终点右坡→F：Return：Ifend ………………（每一个超高变化线元一个If语句） 结果显示： KM=0000.000 PY=0.000 H= 00.000 HL=00.000 I=-1.5 HR=00.000 I=-1.5 KM=0000.000 PY=0.000 QSM= 00.000 QSL= 00.000 QSR= 00.000 Casio5800 计算器数据库型全线高程计算程序(更新2) 一、主程序：2H-SZY Lbl 0: HS ?U: BM-H ?X: SXG= : U+X丄输入后视及水准点设计高程显示视线高程 Lbl 1: “QS”?P: “K X x+xx乂' ?K : Prog “ S.Z”/ C 十100—C: D - 100—D: R Abs(D-C) - 2—T / If D>C :Then 1 —W:Else1—W:lfEnd / If K 在测量岗位工作已经有三个月到时间了，三个月的时间学习和收获了许多，现对这三个月的工作学习做一下总结。 测量工作内容主要有以下两个方面：测量放线（坐标计算），高程控制。 一、测量放线 测量放线到主要技术包括坐标计算和仪器使用。坐标计算包括直线段坐标计算和曲线段坐标计算。 1、直线段坐标计算。直线坐标计算分为中桩坐标计算和边桩坐标计算。 1）中桩坐标计算。根据公式 ααsin ,cos d Y Y d X X +=+=起中起中 d — 所求点到起点距离； α— 该直线坐标方位角。在此顺带详细介绍一下坐标方位角到计算方法： （1）坐标方位角的计算 AB AB A B A B AB x y x x y y ??=--=arctan arctan α当 R y x R y x R y x R y x -360,0,0180,0,0-180,0,0;,0,0?=?+?=??>?αααα；； （2）坐标方位角的推算 ， ， 218021*********βαβααβαβαα-?+=-=+?+=+=B B AB BA B 由此推出：βαα±?+=180后前（“左”→“+”， “右”→“-”），计算中，若α值大于360°，应减去360°；若小于0°，则加上360°。 2）边桩坐标计算 应用公式 )90sin(90cos(?±+=?±+=ααl y y l x x 中边中边）， 进行边桩坐标到计算。北客站为直线车站，坐标计算较简单，现以位于机场线第二段底板的变电所夹层东北角C 点为例进行计算： 以机场线右线为基准来计算中、边桩坐标。已知起点坐标A （22264.4009，11553.2031），终点坐标B （22180.2655，11279.0739），起点里程为YDK0+255.275,C 点里程为YDK0+286.075,偏距为15.33m ，则由以上公式计算C 点坐标： α=arctan(（11279.0739-11553.2031）/（22180.2655-22264.4009）)+180°=252.938°， =中x 22264.4009+（286.075-255.275）*cos252.938°=22255.3640 =中y 11553.2031+（286.075-255.275）*sin252.938°=11523.7586 =c x +15.33*cos （252.938°+90°）=22270.0193 = c y +15.33*sin （252.938°+90°）=11519.2606，则可求出C （22270.0193，11519.2606）。 2、曲线段坐标计算 1）不带缓和曲线的圆曲线中、边桩坐标计算 北 中 x 中 y CASIO5800计算器 公路测量计算程序 程序设计:魏加训 2009.2.28 Casio 5800计算器数据库型万能坐标正反算计算程序 一、主程序：1XY Lbl 0:“1.ZS 2.FS”?→V↙ If V=1: Then “CZ X” ?H: “CZ Y” ?T:Goto 1: Else If V=2: Then Goto 2 :IfEnd: IfEnd ↙ Lbl 1: “K××+×××”?D:?Z: “RJ”?G↙ Prog “P.Z”↙(注:计算另一线路时修改替换此处和FS子程序中的P.Z为对应线路的数据库名称即可) Prog “ZS” ↙ If Z<0: Then Cls:“X（L）=”: “Y（L）=”: Locate 6，1，X : Locate 6，2，Y◢ Pol（X-H，Y-T）: Cls: “S（L）=”: Locate 6，1，I : "F（L）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ If Z=0: Then Cls:“X（Z）=”: “Y（Z）=”: Locate 6，1，X : Locate 6，2，Y : “QXFWJ （Z）=”: 360Frac（（J+360）÷360▼DMS◢ Pol（X-H，Y-T）: Cls: “S（Z）=”: Locate 6，1，I : "F（Z）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ If Z >0: Then Cls:“X（R）=”: “Y（R）=”: Locate 6，1，X : Locate 6,2，Y◢ Pol（X-H，Y-T）: Cls: “S（R）=”: Locate 6，1，I : "F（R）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ Lbl 2: 0→Z:0→G:”X” ?M:”Y” ?I:Prog “FS”: Cls:“K=”：“Z=”：Locate 4，1，D : Locate 4，2，Z◢ Goto 2↙ 二、正算子程序：ZS 5→N: U(E-1-R-1)÷Abs（K-F）→P: Abs（D-F）÷N→Q: 90Q÷π→S: (注:此处5→N是控制计算精度可修改的,一般取值为4～6即可) C+(NPQ+2UR-1)NS→J:1→L↙ A+Q÷6×(Cos (C)+Cos (J) +4∑（Cos (C+((L+0.5)PQ+2UR-1)×(L+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR-1)LS,L,1,(N-1))）+ZCos（J+G) →X : B+Q÷6×(Sin(C)+Sin( J) +4∑（Sin (C+((L+0.5)PQ+2UR-1)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR-1)LS,L,1,(N-1))）+Z Sin（J+G)→Y ↙ 三、反算子程序：FS Lbl 0:Prog “P.Z”:Prog “ZS”↙ (注:计算另一线路时修改替换此处和1XY主程序中的P.Z为对应线路的数据库名称即可) (I-Y)sin(J)+(M-X) cos(J)→ P :D+P→ D ↙ If Abs(P)≥0.001:Then Goto 0 : Else Goto 1 : IfEnd↙ Lbl 1: (I-Y)cos(J)-(M-X) sin(J) →Z 一个excle 模板的制作 在当今社会，excle的使用已经是越来越来频繁了，几乎涉及所有的行业，路桥施工也 不例外。我在某路桥公司曾经负责过某项目部的测量工作。大家都知道，测量最主要的就 是计算了，如坐标、高程、横坡度等。我现在给大家推荐一款我自己编制的关于测量计算 的excel模板。 首先我会跟大家介绍一下模板的作用，然后再一一讲解此模板的制作过程。 首先给大家看一下此模板的界面如下： 也许大家咋一看，切~ 这算啥，我也会做这张表格，实在是太简单了。不错，如果仅 仅是靠手动输入这样子的数字，也许只要懂一点点excle 的人都会制作出这张表格吧。不过，这张表格并不是你表面所看到的仅仅是几个数字而已，其内在的公式才是它的亮点。也许 这样讲大家还不是很清楚，我继续给大家截个图，看看它里面的公式是什么。 大家注意到上面的公式了吗，并不是仅仅是输入数字就完事的，它是一个自定义函数 zbx（），那么后面的都是一样吗？完全正确，后面的都是自定义函数，它们分别是zby（）、sqx（）、hpz（）、hpy（）。也许大家会问，恩，是不错，但是有什么用呢？那让我先给大家 简述一下这个自定义函数的用法。竟然是一个函数，那么它就必须要有一个自变量，这几个 函数的自变量又是什么呢？其实这个模板里面所有函数的自变量只有一个，就是桩号。什么 意思？就是只要你给出任意一个桩号，都能得到其对应的坐标、中桩高程和横坡度。假设我 们要K38+000~K38+200 段落内每隔20M 一个断面所有点的坐标、中桩高程、以及左右横坡。我就用这个模板给大家演示一下（此模板暂时数据只针对黄祁高速公路六标项目部）。 先在桩号那一列把K38+000~K38+200 输入进去，可不要真的把字母“K”和加号 竖曲线全线高程计算程序（数据库模式-精简版） 主程序：ＫSＨ Lb1 0 :＂K=＂? Ｋ:＂D=＂?D:＂H=＂?Ｈ Lb1 1:Prog＂ＫSＨ０＂ If Abs(Ｓ-Ｋ)≥Abs(Ｔ):Then Ｐ+Ｉ(Ｋ-Ｓ)→Ｕ：IfEnd If Abs(Ｓ-Ｋ) < Abs(Ｔ):Then (Ｓ-(Ｔ)-Ｋ)2/(2R)+(P+I(K-S))→Ｕ：IfEnd Lb1 2 ：Ｕ-D-Ｈ→G Locate1,1，＂Ｋ=＂: Locaet3,1，Ｋ: Locaet1,2，＂D=＂: Locaet3,2，D: Locaet1,3，＂Ｈ=＂: Locate3,3，Ｈ: Locaet1,4，＂G=＂: Locaet3,4，G◢ Goto 0 数据库Ｋ：KSＨ０ K< 0 => Stop If K>0:Then 200→S:-12000→R:45.04→P:0.01345→I:62.7→T:If End If K>262.7:Then 520→S:7000→R:46→P:0.003→I: 94.481→T :IfEnd If K> …… 输入提示： If K>前段竖曲线终点里程: Then 本段竖曲线交点里程→S: 曲率半径（凸负凹正）→R: 交点高程→P: 纵坡（上正下负，以小数表示，如0.3%为0.003）→I: 切线长→T : Ifend 符号说明： K=？: 输入待求桩号； D=？: 若计算路面高程输入0；若计算路基高程则输入路面与路基的高差； （如：路面比路基高0.5米，则输入0.5，算出结果就是路基的设计高程。） H=？: 计算设计高程则输入0；若计算高差时；则输入实测高程； G= : 若H输入0，G表示设计高程；若H输入为实测高程，G表示高差（正为填，负为挖） 注：纵坡的设计精度不够会使程序计算出的精度也不够，可以微调纵坡调整（如0.3%的0.003改为0.00301或0.002995），以每段竖曲线的前端来验证，若计算出的高程比设计的高程低，就加大纵坡比，相反则减小纵坡比。 本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《Casio5800交点法程序》与《线元法（积木法）匝道坐标正反算放样程序》，根据歪哥意见“需要的自行修改结合XY框架自己修改为数据库反算程序等”，本人不才，采用最笨的办法将两个程序综合了一下，使之能既能进行交点法正反算，又能进行线元法正反算。在此特别感谢大歪哥！将程序发上来，愿与大家一同交流学习欢迎大家吐口水，只要能进步就行！ 程序由一个主程序ZBZFS和8个子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标正反算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并能对单一线元进行坐标正反算。 1主程序名：ZBZFS（功能：进入计算主程序） 65→Dimz↙ Deg:Fix 3↙ ＂1.JD ZFS 2. ZHADAO ZFS＂? I: I→Z[61]: ＂1.ZHONG SHU JS 2. JS＂? I↙ If I=1: Then Goto1: Else Goto2:IfEnd↙ LbI 1 :If Z[61]=1: Then Prog＂JDYS＂:Else Cls:＂K0＂?A:＂KN＂?L :＂X0＂?U :＂Y0＂?V :＂F0＂?W :＂R0＂?P :＂RN＂?Q:＂ZX:-1,+1,0＂?G:IfEnd↙ LbI 2 :Prog＂JS＂ 2子程序名：JS（功能：选择正算或反算模式） Cls:＂XC＂?H:＂YC＂?Z↙ Cls:＂1.ZS 2.FS＂? I: I=2=>Goto 3↙ LbI 1 : Cls: If Z[61]=1: Then＂JD ZS KX+XXX＂?K :Prog＂4＂: Else ＂ZHADAO ZS KX+XXX＂?K :IfEnd↙ LbI 2: Cls:90→B: Cls:＂RJ Or 0 To K＂?B:B=0 =>Goto 1:＂Z＂?T↙ Prog ＂XY-A＂↙ X+Tcos(M+B)→X↙ Y+Tsin(M+B)→Y↙ 360Frac((M+360)÷360→M↙ Pol(X-H,Y-Z : 360Frac((J+360)÷360→J↙ 2→O: Prog ＂XY-B＂:Goto 2↙ LbI 3 : Cls: If Z[61]=1: Then＂JD FS KN+＂?K:＂X＂?C:＂Y＂?D:Prog＂4＂:Else Cls: ＂ZHADAO FS＂:＂X＂?C:＂Y＂?D:IfEnd↙ LbI 4 :Prog ＂XY-A＂↙ (D-Y)sin(M)+(C-X)cos(M)→H↙ If Abs(H)＞X10-3 :Then K+H→K:Goto 4:IfEnd↙ (D-Y)÷cos(M)→T↙ 3→O: Prog ＂XY-B＂:Goto 3↙ 3子程序名：XY-A（功能：坐标计算程序） 5→N: G(Q-1-P-1)÷Abs（L-A）→F: Abs（K-A）÷N→R: 90R÷π→S:线元法简介
CASIO fx-5800p测量程序
线元法线路坐标正反算程序
fx-5800p全线坐标正反算带高程计算程序(线元法)
线元法万能坐标计算程序
CASIO 5800计算器测量计算程序
[整理]9860G线元法隧道(坐标正反计算、高程)计算程序
交点法线元法坐标计算
《坐标方位角及距离计算小程序》代码——Access实现
CASIO fx5800p全线高程计算程序
5800数据库型全线高程计算程序(2)
测量坐标计算及高程计算
5800坐标高程计算程序
怎样计算高速公路路线坐标及高程
卡西欧5800-竖曲线全线高程计算程序
交点法与线元法