直線段中邊樁坐標(biāo)計(jì)算程序

1、5、直線段中、邊樁坐標(biāo)計(jì)算程序路線直線段的起點(diǎn)里程為ZA=K21+372.446，已知A點(diǎn)的坐標(biāo)（xA=1874.132，yA =4975.213）和B點(diǎn)的坐標(biāo)（xB=1316.429，yB=5139.789），路面左半幅寬為WL=25m，右半幅寬為WR=20m，整樁間距取50m。求直線段的中、邊樁坐標(biāo)。要求：1、寫出數(shù)學(xué)模型2、編寫流程3、程序代碼4、程序運(yùn)行報(bào)告5、程序運(yùn)行結(jié)果=（1）數(shù)學(xué)模型圖1 直線段中、邊樁坐標(biāo)計(jì)算程序計(jì)算原理在圖1中，已知路線直線段起點(diǎn)A 的樁號(hào)ZA及其平面坐標(biāo)(xA, yA)，端點(diǎn)的已知數(shù)據(jù)為B點(diǎn)的平面坐標(biāo)(xB,yB)。設(shè)直線上任意點(diǎn)j的樁號(hào)為Z j，則j點(diǎn)的中

2、樁坐標(biāo)為：xj=xA+(Zj-ZA)cosaAByj=yA+(Zj-ZA)sinaAB左邊樁坐標(biāo)為：xjL=xj+wLcos(aAB-90°)yjL=yj+wLsin(aAB-90°)右邊樁坐標(biāo)為：xjR=xj+wRcos(aAB+90°)yjR=yj+wRsin(aAB+90°)式中，wL為路線左半寬幅，wR為路線右半寬幅。（2）程序流程（3）程序代碼程序中的標(biāo)識(shí)符規(guī)定如下：Z：A點(diǎn)樁號(hào)A：A點(diǎn)的x坐標(biāo)B：A點(diǎn)的y坐標(biāo)C：B點(diǎn)的x坐標(biāo)D：B點(diǎn)的y坐標(biāo)J：AB的方位角L：AB的水平距離K：計(jì)數(shù)變量O：整樁間距l(xiāng)ist1：A點(diǎn)樁號(hào)list2：中樁x坐標(biāo)li

3、st3：中樁y坐標(biāo)list4：左邊樁x坐標(biāo)list5：左邊樁y坐標(biāo)list6：右邊樁x坐標(biāo)list7：右邊樁y坐標(biāo)1“ZBZZB”顯示程序標(biāo)題2Deg:Fix3基本設(shè)置3"A dian ZH(m)="?Z輸入以m 為單位的A 點(diǎn)樁號(hào)4"XA(m)="?A："YA(m)="?B 輸入A點(diǎn)的x坐標(biāo)、y坐標(biāo)5"XB(m)="?C："YB(m)="?D 輸入B點(diǎn)的x坐標(biāo)、y坐標(biāo)6Pol(C-A,D-B) 計(jì)算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 計(jì)數(shù)變量K；存儲(chǔ)A 點(diǎn)樁號(hào)12"A DIAN ZH(m)=":List 1K顯示A 點(diǎn)樁號(hào)13Prog "SUB10"調(diào)子程序計(jì)算A 點(diǎn)中邊樁點(diǎn)坐標(biāo)14If L>0:Then K+1K15Z+LList 1K 計(jì)算B 點(diǎn)樁號(hào)16"B DIAN ZH(m)=":List 1K顯示B 點(diǎn)樁號(hào)17

5、Prog "SUB10"調(diào)子程序計(jì)算B 點(diǎn)中邊樁點(diǎn)坐標(biāo)18IfEnd19"ZHENG ZJJ(m)="?O 輸入整樁間距20K+1K:Int(Z÷O)O+OList 1K 計(jì)算并存儲(chǔ)從A 點(diǎn)開始的第一個(gè)整樁號(hào)21"ZHENG ZH(m)=":List 1K顯示整樁號(hào)22Prog "SUB10"調(diào)子程序計(jì)算中邊樁點(diǎn)坐標(biāo)23Do24K+1K: List 1K-1+OList 1K 計(jì)算并存儲(chǔ)整樁號(hào)25" ZHENG ZH (m)=":List 1K顯示整樁號(hào)26Prog "SUB

6、10"調(diào)子程序計(jì)算并顯示坐標(biāo)27If C>0:Then If List 1K+O-ZL:Then 1Y:Else 0Y:IfEnd整樁號(hào)<直線端點(diǎn)樁號(hào)時(shí)循環(huán),否則結(jié)束28IfEnd29LpWhile Y=130END子程序SUB101A+(List 1K-Z)cos (R)List 2K 保存中樁點(diǎn)x坐標(biāo)到串列2B+(List 1K-Z)sin (R)List 3K保存中樁點(diǎn)y坐標(biāo)到串列3"X(m)=":List 2K顯示中樁點(diǎn)x坐標(biāo)4"Y(m)=":List 3K顯示中樁點(diǎn)y坐標(biāo)525EE為路線左半寬幅6If E>0:The

7、n List 2K+Ecos(R-90)list4K計(jì)算左邊樁x坐標(biāo)，并保存到list4K7List 3K+Esin(R-90)list5K計(jì)算左邊樁y坐標(biāo)，并保存到list5K8"XL(m)=": list4K:"YL(m)=": list5K顯示左邊樁xy坐標(biāo)9IfEnd1020FF為路線右半寬幅11If F>0:Then List 2K+Fcos(R+90)list6K 計(jì)算右邊樁x坐標(biāo)，并保存到list6K12List 3K+Fsin(R+90)list7K計(jì)算右邊樁y坐標(biāo)，并保存到list7K13"XR(m)=": l

8、ist6K:"YR(m)=": list7K顯示右邊樁xy坐標(biāo)14IfEnd15Return（4）程序運(yùn)行報(bào)告按“EXE”鍵 屏幕顯示 A DIAN ZH(M)=？輸入求點(diǎn)樁號(hào)372. 446 屏幕顯示 XA(M)=?輸入求點(diǎn)樁號(hào)1874.132 屏幕顯示 YA(M)=?輸入求點(diǎn)樁號(hào)4975.213 屏幕顯示 XB(M)=?輸入求點(diǎn)樁號(hào)1316.429 屏幕顯示 YB(M)=?輸入求點(diǎn)樁號(hào)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鍵