交點(diǎn)法與線元法

1、本人一直以來想找一個相結(jié)合的坐標(biāo)正反算程序，在網(wǎng)上找了很久，沒能 找到一個較為滿意的，有幸在測量空間看到大歪哥的«Casio5800交點(diǎn)法程序 與線元法（積木法）匝道坐標(biāo)正反算放樣程序，根據(jù)歪哥意見需要的自行 修改結(jié)合XY框架自己修改為數(shù)據(jù)庫反算程序等，本人不才，采用最笨的辦法 將兩個程序綜合了一下，使之能既能進(jìn)行交點(diǎn)法正反算，又能進(jìn)行線元法正反 算。在此特別感謝大歪哥！將程序發(fā)上來，愿與大家一同交流學(xué)習(xí)歡迎大家吐 口水，只要能進(jìn)步就行！程序由一個主程序ZBZFS和8個子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）構(gòu)成，運(yùn)行時只需運(yùn)行主程序即可！本程序適用于單交點(diǎn)對稱型

2、、不對稱型、無緩和曲線單圓曲線型一個交點(diǎn) 范圍內(nèi)（含交點(diǎn)前后有直線段時）的曲線要素核對和坐標(biāo)正反算，手工輸入要 素，對設(shè)計(jì)圖紙的“直線、曲線轉(zhuǎn)角表中交點(diǎn)數(shù)據(jù)進(jìn)行復(fù)核驗(yàn)證，并能對單一 線元進(jìn)行坐標(biāo)正反算。1主程序名：ZBZFS （功能:進(jìn)入計(jì)算主程序）65TDimziZDeg:Fix3lZ” l.JDZFS 2. ZHADAO ZFS ” ? I:TZ61:” l.ZHONGSHU JS2.JS ” ? IlZIf 1=1: Then Gotol: Else Goto2:lfEndiZLbl l:lfZ61=l:Then Prog " JDYSElse Cis:” KO ” ?A:”

3、KN ” ?L:” XO ” ?U:” YO ” ?V :” FO ” ?W :” RO ” ?P:” RN ” ?Q:” ZX:-l,+l,0 ” ?G:IfEndiZLbl 2:Prog ” JS ”2子程序名：JS （功能:選擇正算或反算模式）Cis:° XC " ?H:” YC ” ?ZlZ3/20RJ Or 0 To K = ?p?WHO =>GOS1=Z = ?trpro: XY>=只 X+TCOS(M+B)xkxlz Y+Tsin(M+B)YlZ 36OFrac(M+36oT36OOM0PO_(X,H.-<Z -36OFrac(J+36o)-

4、l-36o4zJ024Z9 Prog = X-5B.7Goto 2lZcFr=l.zs 2FS = ?rrll2HVGOto 3lZLb- lnyr_f zalulHhen = JD ZS Kx+xxx = ?K-Prog - 4=E_se = ZHADAO ZS Kx+xxx = ?K -一fEndlZLb- YCpr904zs?cprLbl3:Cis:lfZ61=l:Then ” JDFS KN+ ” ?K:” X " ?C:” Y ” ?D:Prog ” 4 ” ：Else Cis:” ZHADAO FS ”:” X " ?C:” Y ” ?D:IfEndiZLbl4

5、:Prog ” XY-A ” IZ(D-Y)sin(M)+(C-X)cos(M)HiZIf Abs(H)>X10-3 :Then K+HTK:Goto 4:lfEndlZ(D-Y)vcos(M)TlZ3TO:Prog ” XY-B ” :Goto 3iZ3子程序名:XY-A (功能：坐標(biāo)計(jì)算程序)5TN:G(Q-l-P-l)-rAbs (L-A) TF:Abs (K-A) WNTR:90RETS:W+(FNR+2GPl)NSTM:lTElZU+Rv6x(Cos (W)+Cos (M) +4(Cos (W+(E+0.5)FR+2GP-l)x(E+0.5)S),E,0,(N-l) +2J (

6、Cos (W+(EFR+2GP-1)ES,E,1,(N-1) TXlZV+R4-6x(sin (W)+sin (M) +4(sin (W+(E+0.5)FR+2GPl)x(E+05)S),EQ(N：L) +2J (sin (W+(EFR+2GP-1)ES,E,1,(N-1) TYlZ4子程序名：XY-B (功能：顯示正算或反算結(jié)果)If O=2:TheniZCis:” Kxxx=":“ Z= M :n X 二 n ” Y= ” :Locate 6,1, K :Locate 4, 2, T :Locate 4,3, X :Locate 4,4, Y，lfT=O :Then Cis :n

7、 QF(Z)= n :Locate 8,1/ M:MVDMSIfEndlZCis :” Kxxx=":” S= ”：Locate 6,1, K :Locate 4, 2,1 :” F= ”：J:JTDMSIfEndiZIf O=3:Then u X= ”H Y= H :” Kxxx=":“ Z=":Locate 4,1,C:Locate 4,2,D :Locate 6,3,K :Locate 4,4,T，IfEnd:CIsiZ5子程序名:4 （功能：將交點(diǎn)參數(shù)轉(zhuǎn)為線元計(jì)算參數(shù)）Lbll:IFZ48<0:Then -1Z62:Else:lTZ62:IfEndL

8、bl 2:lf K>Z57:Then Z57TA:Z TL:Z23TU:Z24TV:Z31W :10A45P:10A45Q :OTG:IfEnd IZLbl3:lfK>Zl:Then Z1A :Z TL:Z19TU :Z20TV:Z29W:10A45P:Z46TQ:Z62TG:IfEndiZLbl4:lfK>Z2:ThenZ2A:Z4TL:Z25TU:Z26TV:Z32TW:Z46TP:Z46TQ:Z62TG:IfEndiZLbl5:fK>Z4:ThenZ4->A:Z5TL:Z27 TU:Z28TV:Z33TW :Z46TP:10M5TQ:Z62TG:IfEndi

9、ZLbl6: lfK>Z5:Then Z5A :Z5+1000L:Z21TU:Z22TV:Z30W:10A45P :10A45Q :OTG :IfEndiZ6子程序名：JDYS (功能：輸入交點(diǎn)要素、顯示交點(diǎn)要素及主點(diǎn)坐標(biāo))Cis:” BP ” ?H:HTZ5 7：° K (JD) ” ?K:10/20KZ41:” X (JD) ” ?X:XTZ42:u Y (JD) ”？Y:YTZ43:” LSI u ?B:BTZ44:” LS2 ” ?C:CTZ45:?R:RTZ46:” (ZH) FWJ° ” ?M:MTZ4 7:” a (Z-,Y+) ° ” ?0:

10、OTZ48:Z4 刀+Z48TZ49:Prog " 1 ":Prog ” 2 ” ZCis:u Tl= ” :u T2= ” :° LY=":Locate 4,1, Z50:Locate 42 Z51:Locate 4,3,Z52:Locate 4,4, Z53RCis:” E= ” :Locate 7,1, Z54Cis:” K (QD)= ” :n x 二 n ” Y= ” :” FWJ= ” Locate 7,1,Z57:Locate7,2,Z23:Locate 7,3, Z24:Locate 7,4,Z31，Cis:” K (ZH)= ” :n

11、x 二 n n Y= n :” FWJ= ” :Locate7,l,Zl:Locate 7,2,Z19:Locate 7,3, Z20:Locate 7A Z29Cis:” K (HY)= ” ：n X 二 n ” Y= ” :” FWJ= ” :Locate 7,1,Z2:Locate 7,2,Z25:Locate 7,3, Z26:Locate 7,4,Z32，Cis:” K (QZ)= ” :Locate 7,1,Z3Cis:” K (YH)= ” :n X 二 n H Y= n :” FWJ= ” :Locate 7,1/Z4:Locate 7,2,Z2刀：Locate 7,3/ Z2

12、8:Locate 7,4, Z33，Cis :” K (HZ)二”：“ X= “ :u Y= “ :” FWJ= ” :Locate 7,1,Z5:Locate 7,2,Z21:Locate 7,3/ Z22:Locate 7,4, Z30，7子程序名：1 (功能:計(jì)算交點(diǎn)要素)lfZ48<0:Then -1->Z55:Else 1Z55:IfEnd :Z55*Z48Z56lZZ442v24vZ46-Z44A(4)4-2688vZ46 A(3) TZ62Z45 2 4-24vZ46- Z45A(4) 4-2688vZ46 A(3) TZ刀/Z44v2-Z44A(3 片 24OWZ4

13、62TZ8ZZ45v2-Z45A(3 片 24OWZ462TZ9ZZ8+( Z46+Z7-( Z46+Z6)cos(Z56)vsin(Z56)Z50lZZ9+( Z46+Z6-( Z46+Z7)cos(Z56)vsin(Z56)Z51iZZ46* Z56n-rl80+( Z44+ Z45)豐2TZ52 iZZ46* Z56n-rl80-( Z44+ Z45)豐2TZ53IZ(Z46+(Z6+Z 7)v2)vcos(Z56v2)-Z46->Z54lZZ41-Z50->ZllZlZZl+Z44->Z2lZlZZ2+Z53-r2Z3lZZl+Z52-Z45Z4lZZ4+Z45-&g

14、t;Z5lZ8子程序名：2 (功能：計(jì)算主點(diǎn)坐標(biāo)及切線方位角）Z42-Z50cos(Z4 刀)TZ19:(直緩坐標(biāo))Z43-Z50sin(Z4 刀)TZ20 lZZ4刀TZ :360Frac(Z+360片360TZ29(方位角)Z42+Z51cos(Z49)Z21:(緩直坐標(biāo))Z43+Z51sin(Z49)Z22 lZZ49TZ:360Frac(Z+360片360TZ30(方位角)Zl-Z5刀TLZ (HTZ5刀為前直線起點(diǎn)樁號)Z42-( Z50+L)cos(Z47)TZ23 1/(前直線起點(diǎn)坐標(biāo))Z43-( Z50+L)sin(Z47)->Z24 lZZ4刀TZ :360Frac(Z

15、+360片360TZ31 iZ (方位角)Z44Z12:Z44Z13:Prog ” 3 ” ZZ4-ZlL:90(2L-Z44)vZ46-rnZlllZZ46sin(Zll)+Z8Z14:Z46(l-cos(Zll)+Z 6Z15 iZZ19+Z14cos(Z47)-Z55Z15sin(Z47)Z27iZ (圓緩點(diǎn)坐標(biāo))Z20+Z14sin(Z47)+Z55Z15cos(Z47)->Z28lZZ4 刀+Z55Z11Z:360Frac(Z+360片360TZ33 iZ (方位角)Z2-ZlL:90(2L-Z44)vZ46-rnZ58lZZ46sin(Z58)+Z8Z14:Z46(l-cos

16、(Z58)+Z6Z15 iZZ19+Z14cos(Z47)-Z55Z15sin(Z47)Z25iZ (緩圓點(diǎn)坐標(biāo))Z20+Z14sin(Z47)+Z55Z15cos(Z47)->Z26lZZ47+Z55Z58TZ:360Frac(Z+360片360TZ32 iZ (方位角)9子程序名：3 (主點(diǎn)坐標(biāo)計(jì)算輔助程序)lfZ12=0:ThenOZ14:0Z15:ElseiZZ12-Z12A(5) v404-(Z46*Z13)2+ Z12A(9) 4-3456v(Z46*Z13)A(4) -Z14lZZ12A(3)v6-r(Z46*Z13)-Z12A(7)v336-r(Z46*Z13)A+Z12

17、八(ll)v42240v(Z46*Z13)A(5) TZ15lZIfEndiZ程序說明：1、進(jìn)入程序：1.JDZFS 2. ZHADAO ZFS?選1為交點(diǎn)法正反算(以后操作均 為交點(diǎn)法計(jì)算)，選2為線元法正反算(以后操作均為線元法計(jì)算)2、ZHONGSHUJS2.JS?選1重輸參數(shù)，選2直接進(jìn)入交點(diǎn)法或線元法正反 算(參數(shù)為已輸過的參數(shù))3、參數(shù)輸入：一、交點(diǎn)法已知數(shù)據(jù)輸入：BP?上一交點(diǎn)ZH樁號K (JD) ?交點(diǎn)樁號X (JD) ?交點(diǎn)X坐標(biāo)Y(JD) ?交點(diǎn)Y坐標(biāo)LSI ?第一緩和曲線xxLS2 ?第二緩和曲線xxR ?圓曲線半徑(ZH) FWJ°?交點(diǎn)前(即前交點(diǎn)至本交點(diǎn)也

18、即ZH點(diǎn))的正切線方位角a (乙Y+) ?本交點(diǎn)處線路轉(zhuǎn)角(左轉(zhuǎn)為負(fù)，右轉(zhuǎn)為正，度分秒輸入)交點(diǎn)法計(jì)算要素顯示：T"第一切線長T2二第二切線長L二曲線總長LY二圓曲線長E二曲線外距K(ZH)=直緩點(diǎn)樁號K(HY)=緩圓點(diǎn)樁號K(QZ)=曲中點(diǎn)樁號K(YH)=圓緩點(diǎn)樁號K(HZ)二緩直點(diǎn)樁號二、線元法已知數(shù)據(jù)輸入：KO? KN? RO? RN?FO?XO? YO?ZX?分別為線元起點(diǎn)樁號、終點(diǎn)樁號、起點(diǎn)半 徑、終點(diǎn)半徑、起點(diǎn)切線方位角、起點(diǎn)X坐標(biāo)、起點(diǎn)Y坐標(biāo)、線元轉(zhuǎn)向。4、XC?YC?輸入置鏡點(diǎn)即測站的X,Y坐標(biāo)5、1.ZS2.FS選1正算選2反算程序線元判斷原則：(1) 以道路中線的前進(jìn)方向(即里程增大的方向)區(qū)分左右；(2) 當(dāng)所求點(diǎn)位于中線時，Z=0；當(dāng)位于中線左側(cè)時，Z取負(fù)值；當(dāng)位于中線 右側(cè)時，Z取正值。(3) 當(dāng)線元為直線時，其起點(diǎn)、終點(diǎn)的曲率半徑為無窮大，以10的45次代 替輸入xl0x45o（4）當(dāng)線元為圓曲線時，無論其起點(diǎn)、終點(diǎn)與什么線元相接，其曲率半徑均 等于圓弧的半徑。（