5直线段中、边桩坐标计算程序.docVIP

5、直线段中、边桩坐标计算程序 路线直线段起点里程ZA=K21+372.446，xA=1874.132，yA =4975.213xB=1316.429，yB=5139.789路面左半幅宽WL=25m，右半幅宽WR=20m，整桩间距取50m。 要求： 1、写出数学模型 2、编写流程 3、程序代码 4、程序运行报告 5、程序运行结果 =================================================== （1）数学模型 图1 直线段中、边桩坐标计算程序计算原理 在图1中，已知路线直线段起点A 的桩号ZA及其平面坐标(xA, yA)，端点的已知数据为B点的平面坐标(xB,yB)。设直线上任意点j的桩号为Z j，则j点的中桩坐标为： 左边桩坐标为： 右边桩坐标为： 式中，为路线左半宽幅，为路线右半宽幅。 （2）程序流程 （3）程序代码 程序中的标识符规定如下： Z：A点桩号 A：A点的x坐标 B：A点的y坐标 C：B点的x坐标 D：B点的y坐标 J：A→B的方位角 L：A→B的水平距离 K：计数变量 O：整桩间距 list1：A点桩号 list2：中桩x坐标 list3：中桩y坐标 list4：左边桩x坐标 list5：左边桩y坐标 list6：右边桩x坐标 list7：右边桩y坐标 1 “ZBZZB” 显示程序标题 Deg:Fix3 基本设置 A dian ZH(m)=?Z 输入以m 为单位的A 点桩号 XA(m)=?A：YA(m)=?B 输入A点的x坐标、y坐标 XB(m)=?C：Y(m)=?D 输入B点的x坐标、y坐标 Pol(C-A,D-B) 计算A→B 的方位角与平距 List ans[1]J: list ans[2]L 方位角与平距值赋给J、L If J0:Then J+360R:Else JR:IfEnd 9 BEAR AB(DMS)=:R 显示A→B 的方位角 DIST AB(m)=:L 显示A→B 的水平距离 1K:ZList 1[K] 计数变量K；存储A 点桩号 A DIAN ZH(m)=:List 1[K] 显示A 点桩号 Prog SUB10 调子程序计算A 点中边桩点坐标 If L0:Then K+1K 15 Z+LList 1[K] 计算B 点桩号 B DIAN ZH(m)=:List 1[K] 显示B 点桩号 Prog SUB10 调子程序计算B 点中边桩点坐标 IfEnd 19 ZHENG ZJJ(m)=?O 输入整桩间距 K+1K:Int(Z÷O)O+OList 1[K] 计算并存储从A 点开始的第一个整桩号 ZHENG ZH(m)=:List 1[K] 显示整桩号 Prog SUB10 调子程序计算中边桩点坐标 Do 24 K+1K: List 1[K-1]+OList 1[K] 计算并存储整桩号 ZHENG ZH (m)=:List 1[K] 显示整桩号 Prog SUB10 调子程序计算并显示坐标 If C0:Then If List 1[K]+O-Z≤L:Then 1Y:Else 0Y:IfEnd 整桩号直线端点桩号时循环,否则结束 IfEnd 29 LpWhile Y=1 30 END 子程序——SUB10 1 A+(List 1[K]-Z)cos (R)List 2[K] 保存中桩点坐标到串列 B+(List 1[K]-Z)sin (R)List 3[K] 保存中桩点坐标到串列 X(m)=:List 2[K] 显示中桩点坐标 Y(m)=:List 3[K] 显示中桩点坐标 E E为路线左半宽幅 If E0:Then List 2[K]+Ecos(R-90)list4[K] 计算左边桩x坐标，并保存到list4[K] List 3[K]+Esin(R-90)list5[K] 计算左边桩y坐标，并保存到list5[K] XL(m)=: list4[K]:YL(m)=: list5[K] 显示左边桩坐标 IfEnd 10 20F F为路线右半宽幅 If F0:Then List 2[K]+Fcos(R+90)list6[K] 计算右边桩x坐标，并保存到list6[K] List 3[K]+Fsin(R+90)list7[K] 计算右边桩y坐标，并保存到list7[K] XR(m)=: list6[K]:YR(m)=: list7[K] 显示右边桩坐标 IfEnd 15 Return （4）程序运行报告 按“EXE”键 屏幕显示 A DIAN ZH(M)=？ 输入求点桩号372. 446 屏幕显示 XA(M)=? 输入求点桩号1874.