Mathematica程序 計(jì)算兩線段

1、Mathematica程序計(jì)算兩線段之間的最短距離（*輸入四個(gè)點(diǎn)坐標(biāo)*）（*將坐標(biāo)輸成7/2形式可以得到準(zhǔn)確解，輸成3.5只能得到有一定精確度 的數(shù)值解.*）（*如果兩線段分別為AB,CD,按照以下規(guī)則輸入坐 標(biāo)：xA,yA,xB,yB,xC,yC,xD,yD或者xA,yA,xB,yB,xC,yC,xD,yD或者兩種混和輸入，必須保證按 照順序輸入*）zuobiao = 0,0, 13,9, 6,6, 0, 8;zuobiao = PartitionFlattenzuobiao, 2;（*下面進(jìn)行坐標(biāo)旋轉(zhuǎn)，目的是讓兩條直線均不與坐標(biāo)軸平行，省去以后 很多分情況的麻煩，一勞永逸*）xuanzj

2、= 0, Pi/6, Pi/4;brr = ;err = ;juli = ;frr = ;chuizuf =.;Fori = 1, i = Lengthxuanzj, i+,Forj = 1, j =0 & (brr31 - brr51)*(brr51- brr41) = 0,Print-兩線段相交，最短距離為0;xiangjiao = 1,xiangjiao = 0;(*下面先定義一個(gè)函數(shù)chuizu用來計(jì)算點(diǎn)到直線的垂足坐標(biāo)*)chuizuzuobiao_:=Modulearr,arr = Flattenzuobiao;(-(-arr1 arr3A2 + 2 arr1 arr3 arr5-

3、arr1 arr5A2 - arr3 arr2 arr4 +arr5 arr2 arr4 - arr5 arr4A2 +arr3 arr2 arr6 - arr5 arr2 arr6 +arr3 arr4 arr6 + arr5 arr4 arr6-arr3 arr6A2)/(arr3A2 - 2 arr3 arr5 +arr5A2 + arr4A2 - 2 arr4 arr6 +arr6A2), -(-arr1 arr3 arr4 +arr1 arr5 arr4 + arr3 arr5 arr4-arr5A2 arr4 - arr2 arr4A2 +arr1 arr3 arr6 - arr3

4、A2 arr6-arr1 arr5 arr6 + arr3 arr5 arr6 +2 arr2 arr4 arr6-arr2 arr6A2)/(arr3A2 - 2 arr3 arr5 + arr5A2 + arr4A2 - 2 arr4 arr6 + arr6A2)；(*下面是兩線段線不相交時(shí)，計(jì)算各個(gè)端點(diǎn)到所對(duì)線段的最短距離9crr = brr1, brr3, brr4, brr2, brr3,brr4, brr3, brr1, brr2, brr4, brr1, brr2;Ifxiangjiao = 0,Fori = 1, i =0 & (crri31 - chuizuf1)/(crri

5、31- crri21) = 0,frr = Appendfrr, crri1, chuizuf;juli =Appendjuli, SimplifyNormcrri1 - chuizuf(*A勺垂足 F 在CD內(nèi)，取線段AF距離*);If(crri 21 - chuizuf1)/(crri21- crri31) = 0,frr = Appendfrr, crri1, crri 2 ;juli =Appendjuli,SimplifyNormcrri1 - crri2(*A的垂足F在C 的 外側(cè)，取線段AC距離*);crri31) =0 & (crri31 - chuizuf1)/(crri31

6、- crri21) 0,frr = Appendfrr, crri1, crri3;juli =Appendjuli,SimplifyNormcrri1 - crri3(*A勺垂足F 在D 的 外側(cè)，取線段AD距離*);Print準(zhǔn)確解;PrintSortjuli, Less1;Print數(shù)值解;PrintSortjuli / N1;(*下面是作圖*)Forj = 1, j = Lengthbrr, j+,err = Appenderr, brrj1*CosAlpha +brrj2*SinAlpha, -brrj1*SinAlpha +brrj2*CosAlpha;err = Partitionerr, 2;Iffrr != ,grr = TableNormfrri1 - frri 2 / N, i, 1, Lengthfrr;Fori = 1,