直线段中边桩坐标计算程序

1、5、直线段中、边桩坐标计算程序路线直线段的起点里程为ZA=K21+372.446，已知A点的坐标（xA=1874.132，yA =4975.213）和B点的坐标（xB=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、桩坐标为：xj=xA+(Zj-ZA)cosaAByj=yA+(Zj-ZA)sinaAB左边桩坐标为：xjL=xj+wLcos(aAB-90°)yjL=yj+wLsin(aAB-90°)右边桩坐标为：xjR=xj+wRcos(aAB+90°)yjR=yj+wRsin(aAB+90°)式中，wL为路线左半宽幅，wR为路线右半宽幅。（2）程序流程（3）程序代码程序中的标识符规定如下：Z：A点桩号A：A点的x坐标B：A点的y坐标C：B点的x坐标D：B点的y坐标J：AB的方位角L：AB的水平距离K：计数变量O：整桩间距list1：A点桩号list2：中桩x坐标li

3、st3：中桩y坐标list4：左边桩x坐标list5：左边桩y坐标list6：右边桩x坐标list7：右边桩y坐标1“ZBZZB”显示程序标题2Deg:Fix3基本设置3"A dian ZH(m)="?Z输入以m 为单位的A 点桩号4"XA(m)="?A："YA(m)="?B 输入A点的x坐标、y坐标5"XB(m)="?C："YB(m)="?D 输入B点的x坐标、y坐标6Pol(C-A,D-B) 计算AB 的方位角与平距7List ans1J: list ans2L方位角与平距值赋给J、L8If

4、 J<0:Then J+360R:Else JR:IfEnd9"BEAR AB(DMS)=":R显示AB 的方位角10"DIST AB(m)=":L显示AB 的水平距离111K:ZList 1K 计数变量K；存储A 点桩号12"A DIAN ZH(m)=":List 1K显示A 点桩号13Prog "SUB10"调子程序计算A 点中边桩点坐标14If L>0:Then K+1K15Z+LList 1K 计算B 点桩号16"B DIAN ZH(m)=":List 1K显示B 点桩号17

5、Prog "SUB10"调子程序计算B 点中边桩点坐标18IfEnd19"ZHENG ZJJ(m)="?O 输入整桩间距20K+1K:Int(Z÷O)O+OList 1K 计算并存储从A 点开始的第一个整桩号21"ZHENG ZH(m)=":List 1K显示整桩号22Prog "SUB10"调子程序计算中边桩点坐标23Do24K+1K: List 1K-1+OList 1K 计算并存储整桩号25" ZHENG ZH (m)=":List 1K显示整桩号26Prog "SUB

6、10"调子程序计算并显示坐标27If C>0:Then If List 1K+O-ZL:Then 1Y:Else 0Y:IfEnd整桩号<直线端点桩号时循环,否则结束28IfEnd29LpWhile Y=130END子程序SUB101A+(List 1K-Z)cos (R)List 2K 保存中桩点x坐标到串列2B+(List 1K-Z)sin (R)List 3K保存中桩点y坐标到串列3"X(m)=":List 2K显示中桩点x坐标4"Y(m)=":List 3K显示中桩点y坐标525EE为路线左半宽幅6If E>0:The

7、n List 2K+Ecos(R-90)list4K计算左边桩x坐标，并保存到list4K7List 3K+Esin(R-90)list5K计算左边桩y坐标，并保存到list5K8"XL(m)=": list4K:"YL(m)=": list5K显示左边桩xy坐标9IfEnd1020FF为路线右半宽幅11If F>0:Then List 2K+Fcos(R+90)list6K 计算右边桩x坐标，并保存到list6K12List 3K+Fsin(R+90)list7K计算右边桩y坐标，并保存到list7K13"XR(m)=": l

8、ist6K:"YR(m)=": list7K显示右边桩xy坐标14IfEnd15Return（4）程序运行报告按“EXE”键 屏幕显示 A DIAN ZH(M)=？输入求点桩号372. 446 屏幕显示 XA(M)=?输入求点桩号1874.132 屏幕显示 YA(M)=?输入求点桩号4975.213 屏幕显示 XB(M)=?输入求点桩号1316.429 屏幕显示 YB(M)=?输入求点桩号5139.789 屏幕显示 JULI AB(M)=581.479按“EXE键 屏幕显示 FWJ AB(Deg)=163°3331.95按“EXE键 屏幕显示 A DIAN ZH=372.446按“EXE键 屏幕显示 X(M)=1874.132按“EXE键 屏幕显示 Y(M)=4975.213按“EXE键 屏幕显示 XL(M)=1881.208按“EXE键