a題飛行管理問題

1、1995年A題飛行管理問題A題 一個飛行管理問題 在約10,000米高空的某邊長160公里的正方形區(qū)域內(nèi)，經(jīng)常有若干架飛機作水平飛行.區(qū)域內(nèi)每架飛機的位置和速度向量均由計算機記錄其數(shù)據(jù)，以便進行飛行管理.當一架欲進入該區(qū)域的飛機到達區(qū)域邊緣時，記錄其數(shù)據(jù)后，要立即計算并判斷是否會與區(qū)域內(nèi)的飛機發(fā)生碰撞.如果會碰撞，則應(yīng)計算如何調(diào)整各架（包括新進入的）飛機飛行的方向角.以避免碰撞.現(xiàn)假定條件如下：1） 不碰撞的標準為任意兩架飛機的距離大于8公里2） 飛機飛行方向角調(diào)整的幅度不應(yīng)超過30度3） 所有飛機飛行速度均為每小時800公里4） 進入該區(qū)域的飛機在到達區(qū)域邊緣時，與區(qū)域內(nèi)飛機的距離應(yīng)在60公

2、里以上5） 最多需考慮6架飛機6） 不必考慮飛機離開此區(qū)域后的狀況. 請你對這個避免碰撞的飛行管理問題建立數(shù)學模型.列出計算步驟，對以下數(shù)據(jù)進行計算（方向角誤差不超過0.01度）.要求飛機飛行方向角調(diào)整的幅度盡量小. 設(shè)該區(qū)域4個頂點的座標為（0,0）,(160,0),(160,160),(0,160)記錄數(shù)據(jù)為： 飛機編號 橫座標X 縱座標Y 方向角（度） 1 150 140 243 2 85 85 236 3 150 155 220.5 4 145 50 159 5 130 150 230 新進入 0 0 52注：方向角指飛行方向與X軸正向的夾角. 試根據(jù)實際應(yīng)用背景對你的模型進行評價與推

3、廣.參考解答1 問題分析 根據(jù)題目的條件，可將飛機飛行的空域視為二維平面xoy中的一個正方形區(qū)域，頂點（0,0）,(160,0),(160,160),(0,160).各架飛機的飛行方向角為飛行方向與x軸正向夾角（轉(zhuǎn)角）.根據(jù)兩飛機不碰撞的標準為二者距離大于8km，可將每架飛機視為一個以飛機坐標點為圓心、以4km為半徑的圓狀物體（每架飛機在空域中的狀態(tài)由圓心的位置矢量和飛行速度矢量確定）.這樣兩架飛機是否碰撞就化為兩圓在運行過程中是否相交的問題.兩圓是否相交只要討論它們的相對運動即可.2.模型假設(shè)(1) 飛機進入?yún)^(qū)域邊緣時,立即作出計算,每架飛機按照計算后的指示立即作方向角改變;(2) 每架飛機

4、在整個過程中至多改變一次方向;(3) 忽略飛機轉(zhuǎn)向的影響(轉(zhuǎn)彎半徑和轉(zhuǎn)彎時間的影響);(4) 新飛機進入空域時,已在空域內(nèi)部飛行的飛機的飛行方向已調(diào)合適,不會碰撞;(5) 對每架飛機方向角的相同調(diào)整量的滿意程度是一樣的.3.模型的建立(1) 圓狀模型. 由前面的分析將飛機作為圓狀模型進行研究.兩圓不相交,則表明不會發(fā)生碰撞事故;若兩圓相交,則表明會發(fā)生碰撞事故.為了研究兩飛機相撞問題,采用相對速度作為研究對象,因為飛機是否相撞的關(guān)鍵是相對速度,圖10-3給出任意兩架飛機間的關(guān)系其中符號含義如下:i,j-第i,第j架飛機的圓心;aij-第i架飛機與第j架飛機的碰撞角,是兩圓的切線交角中指向圓的那

5、個角的一 半,aij=aji;ij-第i架飛機相對于第j架飛機的相對飛行速度;lij-第i架飛機與第j架飛機的圓心距;ij-第i架飛機對于第j架飛機的相對速度與兩架飛機圓心連線的交角.規(guī)定以第i架 飛機為原點,ij連線從i指向j為正方向,逆時針旋轉(zhuǎn)為正角,順時針旋轉(zhuǎn)為負角; AB,CD為兩圓的公切線,mi/AB,ni/CD.另外再引入記號:i-第i架飛機的飛行方向與直角坐標xoy中x軸正向的夾角(轉(zhuǎn)角);xi-第i架飛機在坐標xoy中的位置矢量;i-第i架飛機的飛行速度矢量. 由圖10-3中的關(guān)系得到兩飛機不相撞(兩圓不相交)的充要條件是|ij |>aij.當|ij |aij時,則通過調(diào)

6、整兩飛機的方向角I, j,使飛機不相撞.(2) 決策目標. 題目要求飛機飛行方向角調(diào)整的幅度盡量小,這個盡量小是針對每架飛機而言,同時也要求整體滿意程度(即對管理層而言,應(yīng)使每架飛機的調(diào)整都盡量的小).因此構(gòu)造目標函數(shù)時,可以認為若對方向角調(diào)整量最大的飛機而言,其調(diào)整量可滿意,則由假設(shè)(5)對其余飛機調(diào)整量均可滿意.即要求每架飛機的調(diào)整量都小于某個數(shù)(0).故可取目標函數(shù)為求其最小值min.(3) 由圓狀模型導出的方程. 首先討論相對飛行速度方向角ij的改變量與第i,第j兩架飛機飛行方向角改變量i, j的關(guān)系.由題目條件知,對第i架飛機|i|=800=A(km).于是可用復數(shù)表示飛機速度i=.

7、設(shè)第i,j兩架飛機飛行方向改變前的速度分別為,改變飛行方向后的速度為, 則飛行方向改變前后的相對速度分別為= = 即與交角相之差為.將其歸納為定理 對第i,j架飛機,其相對速度方向ij的改變量ij等于兩飛機飛行方向角改變量之和的一半. 由題目的要求調(diào)整飛行方向角時不能超過30°即|i|30 , i=1,2,6要保證調(diào)整飛行方向后飛機不碰撞，應(yīng)有|ij+ij |>aij由前面構(gòu)造的目標函數(shù)為min 030總結(jié)以上得如下優(yōu)化模型min （1）s.t. |ij+ij |>aij, (2)|i|, I=1,6 (3)|i|30, I=1,6 (4)030 (5)(4) 線性規(guī)劃模

8、型.將上述優(yōu)化模型進行化簡,可轉(zhuǎn)化為線性規(guī)劃模型.當ij>0時,(10.2)式可化為ij+ij>aij;當ij<0時,(10.2)式可化為ij+ij>-aij;由于i可正負,為使各變量均為正,引入新的變量,使i=- , 030 , 030于是條件(3),( 4)可化為-30-30-這樣,優(yōu)化模型(1)(5)就轉(zhuǎn)化為如下線性規(guī)劃模型min (6)s.t.ij>0時, -+->2ij-2ij , (7)ij<0時, -+-<-2ij-2ij (8)-30 (9)-30 (10)- (11)- (12) , , , 0 (14)其中ij,ij可由題中已

9、知的參數(shù)計算得到,ij=arg(i-j)-arg(j-xi)ij=arcsin(8/|xi-xj|)4.模型求解(1) 記錄各飛機狀態(tài)(位置矢量、速度矢量);(2) 計算任兩架飛機間的參數(shù)ij, ij;(3) 利用計算線性規(guī)劃的軟件(如LINDO)求解(10.6)(10.14).這一 步也可使用單純形法自己編程求解.5.結(jié)果檢驗對題目所給實例進行計算得如下調(diào)整方案1=0 , 2=0, 3=1.814732,4=-0.743155, 5=0, 6=1.814732.各飛行方向角按此方案調(diào)整后,系統(tǒng)各架飛機均滿足|ij|>ij(即不會相撞).其中有些飛機對有|ij|-ij<0.01(0.01°是題目要求的計算精度).如果希望|ij|ij+0.01,只須將模型中的ij用=ij+0.01代替即可.將調(diào)整后各量再代入模型進行計算得min=0即此時無需再改動.經(jīng)模擬程序運行可觀察動態(tài)結(jié)果是正確的.6.模型評價與推廣(1) 此模型采用圓狀模型分析碰撞問題是合理的,同時采用相對速度作為判斷標準,即體現(xiàn)了碰撞的本質(zhì)(相對運動),又簡化了模型的計算;(2) 建模中用了適當?shù)暮喕?將一個復雜的非線性規(guī)劃問題簡化為線性規(guī)劃問題.既求到合理的解,又提高了運算速度.這對于解決高速運行的飛機碰撞問題是十分重要的.此模型對題目所提供的例子計算得出的結(jié)果是令