建模案例飛行管理問(wèn)題

1、建模案例飛行管理問(wèn)題一個(gè)飛行管理問(wèn)題一個(gè)飛行管理問(wèn)題 1.1 問(wèn)題描述問(wèn)題描述 在約在約10000米高空的某邊長(zhǎng)為米高空的某邊長(zhǎng)為160km的正方形區(qū)域內(nèi)，的正方形區(qū)域內(nèi)，經(jīng)常有若干架飛機(jī)做水平飛行，區(qū)域內(nèi)每架飛機(jī)的位經(jīng)常有若干架飛機(jī)做水平飛行，區(qū)域內(nèi)每架飛機(jī)的位置和速度向量均由計(jì)算機(jī)記錄其數(shù)據(jù)，以便進(jìn)行飛行置和速度向量均由計(jì)算機(jī)記錄其數(shù)據(jù)，以便進(jìn)行飛行管理，當(dāng)一架欲進(jìn)入該區(qū)域的飛機(jī)到達(dá)區(qū)域邊緣時(shí)，管理，當(dāng)一架欲進(jìn)入該區(qū)域的飛機(jī)到達(dá)區(qū)域邊緣時(shí)，記錄其數(shù)據(jù)后，要立即計(jì)算并判斷是否會(huì)與區(qū)域內(nèi)的記錄其數(shù)據(jù)后，要立即計(jì)算并判斷是否會(huì)與區(qū)域內(nèi)的飛機(jī)發(fā)生碰撞。飛機(jī)發(fā)生碰撞。 現(xiàn)假定條件如下：現(xiàn)假定條件如下

2、： 不碰撞的標(biāo)準(zhǔn)為任意兩架飛機(jī)的距離大于不碰撞的標(biāo)準(zhǔn)為任意兩架飛機(jī)的距離大于8km； 飛機(jī)飛行方向角調(diào)整幅度不應(yīng)超過(guò)飛機(jī)飛行方向角調(diào)整幅度不應(yīng)超過(guò)30 ； 所有飛機(jī)飛行速度均為所有飛機(jī)飛行速度均為800kmh； 進(jìn)入該區(qū)域的飛機(jī)在到達(dá)該區(qū)域邊緣時(shí)，與進(jìn)入該區(qū)域的飛機(jī)在到達(dá)該區(qū)域邊緣時(shí)，與區(qū)域內(nèi)飛機(jī)的距離應(yīng)在區(qū)域內(nèi)飛機(jī)的距離應(yīng)在60km以上；以上； 最多需考慮最多需考慮6架飛機(jī)；架飛機(jī)； 不必考慮飛機(jī)離開(kāi)此區(qū)域后的情況。不必考慮飛機(jī)離開(kāi)此區(qū)域后的情況。 請(qǐng)你對(duì)這個(gè)請(qǐng)你對(duì)這個(gè)避免碰撞避免碰撞的飛行管理問(wèn)題建立數(shù)學(xué)模型，列出計(jì)算的飛行管理問(wèn)題建立數(shù)學(xué)模型，列出計(jì)算步驟，對(duì)以下數(shù)據(jù)進(jìn)行計(jì)算（方向角誤差

3、不超過(guò)步驟，對(duì)以下數(shù)據(jù)進(jìn)行計(jì)算（方向角誤差不超過(guò)0.01 ），要求飛），要求飛機(jī)飛行方向角調(diào)整的幅度機(jī)飛行方向角調(diào)整的幅度盡量小盡量小。 該區(qū)域四個(gè)定點(diǎn)的坐標(biāo)為（該區(qū)域四個(gè)定點(diǎn)的坐標(biāo)為（0，0）、（）、（160，0）、（）、（160，160）、（）、（0，160）。記錄數(shù)據(jù)見(jiàn)表）。記錄數(shù)據(jù)見(jiàn)表21。 表表21 飛機(jī)位置和方向角記錄數(shù)據(jù)飛機(jī)位置和方向角記錄數(shù)據(jù)飛機(jī)編號(hào)橫坐標(biāo)縱坐標(biāo)方向角（飛機(jī)編號(hào)橫坐標(biāo)縱坐標(biāo)方向角（ ）飛機(jī)編號(hào)橫坐標(biāo)）飛機(jī)編號(hào)橫坐標(biāo)x縱坐標(biāo)縱坐標(biāo)y方向角（方向角（ ）飛機(jī)編號(hào)橫坐標(biāo)縱坐標(biāo)方向角（） 飛機(jī)編號(hào)橫坐標(biāo)x縱坐標(biāo)y方向角（）115014024341455015928585

4、236513015023031551552205新進(jìn)入0052說(shuō)明：方向角指飛行方向與x軸正向的夾角。試根據(jù)實(shí)際應(yīng)用背景對(duì)你的模型進(jìn)行評(píng)價(jià)和推廣 * 對(duì)問(wèn)題仔細(xì)閱讀, 首先抓住題目中的關(guān)鍵詞“管理”進(jìn)行聯(lián)想. 抓住諸如“碰撞”、“調(diào)整”、“避免碰撞”、“立即”、“判斷”等等詞語(yǔ). * 聯(lián)系解決問(wèn)題的方案,不加約束繼續(xù)聯(lián)想，再將關(guān)鍵詞搭配起來(lái). 1. 問(wèn)題的前期分析160km160km飛行位置示意圖立即 判斷 碰撞 條件 實(shí)時(shí) 算法 避免 碰撞 調(diào)整 方向角 實(shí)時(shí) 幅度盡量小 相對(duì)距離優(yōu)化問(wèn)題優(yōu)化算法優(yōu)化調(diào)整方案問(wèn)題的初步理解和想法 飛行管理問(wèn)題是優(yōu)化問(wèn)題,在調(diào)整方向角的幅度盡量小的同時(shí)，還必須

5、注意調(diào)整方案及算法的實(shí)時(shí)性.2. 問(wèn)題探究問(wèn)題探究 （1）優(yōu)化問(wèn)題的目標(biāo)函數(shù)為何？）優(yōu)化問(wèn)題的目標(biāo)函數(shù)為何？方向角調(diào)整的盡量小方向角如何表示方向角的概念是什么方向角的平方和目標(biāo)函數(shù) 任意兩架飛機(jī)的距離大于任意兩架飛機(jī)的距離大于8 公里；公里； 飛機(jī)飛行方向角調(diào)整的幅度不應(yīng)超過(guò)飛機(jī)飛行方向角調(diào)整的幅度不應(yīng)超過(guò) 30 ；l（2）優(yōu)化問(wèn)題的約束條件為何？?jī)牲c(diǎn)間距離表示方法判斷避免碰撞的依據(jù)把飛機(jī)視為點(diǎn)總結(jié)：目標(biāo)函數(shù)和約束條件中都含有方向角。（3）分析6架飛機(jī)目前碰撞情況 描點(diǎn)作圖；描點(diǎn)作圖； 分析飛機(jī)在飛行區(qū)域的時(shí)間；分析飛機(jī)在飛行區(qū)域的時(shí)間； 判斷判斷5架飛機(jī)此刻的情況；架飛機(jī)此刻的情況； 判斷新

6、進(jìn)入飛機(jī)與其他判斷新進(jìn)入飛機(jī)與其他5架飛機(jī)碰撞的情況。架飛機(jī)碰撞的情況。（4）求解方法）求解方法 特殊到一般：先考慮特殊到一般：先考慮2架，然考慮架，然考慮3架架 優(yōu)化問(wèn)題為非線性規(guī)劃問(wèn)題，編程求解優(yōu)化問(wèn)題為非線性規(guī)劃問(wèn)題，編程求解?？偨Y(jié)：初等算法和高等算法都可。創(chuàng)新之處：算法。模型一及求解模型一及求解 模型建立模型建立 這個(gè)問(wèn)題顯然是一個(gè)優(yōu)化問(wèn)題。設(shè)第這個(gè)問(wèn)題顯然是一個(gè)優(yōu)化問(wèn)題。設(shè)第i 架飛架飛機(jī)在調(diào)整時(shí)的方向角為機(jī)在調(diào)整時(shí)的方向角為 （題目中已給（題目中已給出），調(diào)整后的方向?yàn)槌觯?，調(diào)整后的方向?yàn)?，題目中就是要求飛機(jī)飛行方向角調(diào)整的幅度題目中就是要求飛機(jī)飛行方向角調(diào)整的幅度盡量小，因此有

7、化的目的函數(shù)可以是：盡量小，因此有化的目的函數(shù)可以是： （1）為了建立這個(gè)問(wèn)題的優(yōu)化模型，只需要明確約束條件就可以了。一個(gè)簡(jiǎn)為了建立這個(gè)問(wèn)題的優(yōu)化模型，只需要明確約束條件就可以了。一個(gè)簡(jiǎn)單的約束是飛機(jī)飛行方向角調(diào)整的幅度不應(yīng)超過(guò)單的約束是飛機(jī)飛行方向角調(diào)整的幅度不應(yīng)超過(guò)30 ，即，即 （2） 題目中要求進(jìn)入該區(qū)域的飛機(jī)在到達(dá)該區(qū)域邊緣時(shí)，與區(qū)域內(nèi)的飛題目中要求進(jìn)入該區(qū)域的飛機(jī)在到達(dá)該區(qū)域邊緣時(shí)，與區(qū)域內(nèi)的飛機(jī)的距離應(yīng)在機(jī)的距離應(yīng)在60km以上。這個(gè)條件是個(gè)初始條件，很容易驗(yàn)證目以上。這個(gè)條件是個(gè)初始條件，很容易驗(yàn)證目前所給的數(shù)據(jù)是滿足的，因此本模型中可以不予考慮。剩下的關(guān)鍵前所給的數(shù)據(jù)是滿足的

8、，因此本模型中可以不予考慮。剩下的關(guān)鍵是是 要滿足題目中描述的任意兩架位于該區(qū)域內(nèi)的飛機(jī)的距離應(yīng)該大要滿足題目中描述的任意兩架位于該區(qū)域內(nèi)的飛機(jī)的距離應(yīng)該大于于8km。但這個(gè)問(wèn)題的難點(diǎn)在于飛機(jī)是動(dòng)態(tài)的，這個(gè)約束不好直。但這個(gè)問(wèn)題的難點(diǎn)在于飛機(jī)是動(dòng)態(tài)的，這個(gè)約束不好直接描述，為此我們首先需要描述每架飛機(jī)的飛行軌跡。接描述，為此我們首先需要描述每架飛機(jī)的飛行軌跡。 記飛機(jī)飛行速率為（記飛機(jī)飛行速率為（800kmh），以當(dāng)前），以當(dāng)前時(shí)刻為時(shí)刻為0時(shí)刻。設(shè)第時(shí)刻。設(shè)第 架飛機(jī)在調(diào)整時(shí)的位架飛機(jī)在調(diào)整時(shí)的位置坐標(biāo)為置坐標(biāo)為 （已知條件），時(shí)刻的位（已知條件），時(shí)刻的位置坐標(biāo)為置坐標(biāo)為 ，則，則 （3）

9、 如果要嚴(yán)格表示兩架位于該區(qū)域內(nèi)的飛機(jī)的如果要嚴(yán)格表示兩架位于該區(qū)域內(nèi)的飛機(jī)的距離應(yīng)大于距離應(yīng)大于8km，則需要考慮每架飛機(jī)在區(qū)，則需要考慮每架飛機(jī)在區(qū)域內(nèi)的飛行時(shí)間的長(zhǎng)度。記域內(nèi)的飛行時(shí)間的長(zhǎng)度。記Ti 為第為第 架飛架飛機(jī)飛出區(qū)域的時(shí)間，即機(jī)飛出區(qū)域的時(shí)間，即 （4） 記記 時(shí)刻第時(shí)刻第 架飛機(jī)與第架飛機(jī)與第 架飛機(jī)的距離架飛機(jī)的距離為為 ，并記，并記 ，這時(shí)在區(qū)，這時(shí)在區(qū)域內(nèi)飛機(jī)不相撞的約束條件就變成了域內(nèi)飛機(jī)不相撞的約束條件就變成了 （5） 其中其中 （6） 此外，經(jīng)過(guò)計(jì)算可以得到此外，經(jīng)過(guò)計(jì)算可以得到 （7） （8） （9） （10） 所以所以 是一個(gè)關(guān)于是一個(gè)關(guān)于 t 的二次函數(shù)，

10、表示的的二次函數(shù)，表示的是一條開(kāi)口向上的拋物線。當(dāng)是一條開(kāi)口向上的拋物線。當(dāng) 即即 （記為（記為 ）時(shí)，）時(shí)， 函數(shù)取最小函數(shù)取最小值值 。注意到。注意到 （初始時(shí)刻不相（初始時(shí)刻不相撞），如果撞），如果 （即（即 ）則此時(shí)約束條）則此時(shí)約束條件（件（5）一定成立，所以）一定成立，所以 如果如果 且且 ，只要在右端點(diǎn)的函數(shù)值，只要在右端點(diǎn)的函數(shù)值非負(fù)即可，即非負(fù)即可，即 （11） 如果如果 且且 ，只需要，只需要 求最小值求最小值 即可，即即可，即 （12） 實(shí)際上，約束（實(shí)際上，約束（11）表示的是）表示的是 在右端點(diǎn)在右端點(diǎn)的函數(shù)值非負(fù)，這個(gè)約束在（的函數(shù)值非負(fù)，這個(gè)約束在（12）的條件）

11、的條件下也是自然成立的，所以可以是對(duì)約束下也是自然成立的，所以可以是對(duì)約束（11）不再附加且的條件。）不再附加且的條件。 ( )ijft( )ijft 于是我們的模型就是于是我們的模型就是 （13） （14） （15） （16） 模型求解模型求解上面這是一個(gè)非線性規(guī)劃模型，雖然是嚴(yán)格滿足題目要求的模型，但上面這是一個(gè)非線性規(guī)劃模型，雖然是嚴(yán)格滿足題目要求的模型，但得到的模型邏輯關(guān)系比較復(fù)雜，約束（得到的模型邏輯關(guān)系比較復(fù)雜，約束（16）是在一定條件下才成）是在一定條件下才成立的約束，而且其中的計(jì)算式（立的約束，而且其中的計(jì)算式（4）也含有相當(dāng)復(fù)雜的關(guān)系式，）也含有相當(dāng)復(fù)雜的關(guān)系式，使用使用LI

12、NGO軟件不太容易將模型很方便的輸入，因?yàn)檫壿嬏庈浖惶菀讓⒛Ｐ秃芊奖愕妮斎?，因?yàn)檫壿嬏幚聿皇抢聿皇荓INGO的優(yōu)勢(shì)所在。即使想辦法把這個(gè)模型輸入到的優(yōu)勢(shì)所在。即使想辦法把這個(gè)模型輸入到LINGO，也，也不一定能求出好的解（筆者嘗試過(guò)，但是不一定能求出好的解（筆者嘗試過(guò)，但是LINGO運(yùn)行時(shí)有時(shí)會(huì)出運(yùn)行時(shí)有時(shí)會(huì)出現(xiàn)系統(tǒng)內(nèi)部錯(cuò)誤，可能是系統(tǒng)有問(wèn)題，無(wú)法繼續(xù)求解）。而且，現(xiàn)系統(tǒng)內(nèi)部錯(cuò)誤，可能是系統(tǒng)有問(wèn)題，無(wú)法繼續(xù)求解）。而且，在實(shí)時(shí)飛行調(diào)度中顯然需要快速求解，所以下面我們想辦法簡(jiǎn)在實(shí)時(shí)飛行調(diào)度中顯然需要快速求解，所以下面我們想辦法簡(jiǎn)化模型。化模型。 這個(gè)模型麻煩之處就在于，要求嚴(yán)格表示兩這個(gè)模型

13、麻煩之處就在于，要求嚴(yán)格表示兩架飛機(jī)的飛行距離應(yīng)大于架飛機(jī)的飛行距離應(yīng)大于8km，所以需要考，所以需要考慮每架飛機(jī)在區(qū)域內(nèi)的飛行時(shí)間的長(zhǎng)度，比慮每架飛機(jī)在區(qū)域內(nèi)的飛行時(shí)間的長(zhǎng)度，比較繁瑣。注意到區(qū)域?qū)蔷€的長(zhǎng)度只較繁瑣。注意到區(qū)域?qū)蔷€的長(zhǎng)度只有有 ，任何一架飛機(jī)在所考慮的區(qū)域內(nèi)，任何一架飛機(jī)在所考慮的區(qū)域內(nèi)停留的時(shí)間不會(huì)超過(guò)停留的時(shí)間不會(huì)超過(guò) 。因此這里我們簡(jiǎn)化一下問(wèn)題；因此這里我們簡(jiǎn)化一下問(wèn)題； 不再單獨(dú)考慮每架飛機(jī)在區(qū)域內(nèi)停留的時(shí)間，不再單獨(dú)考慮每架飛機(jī)在區(qū)域內(nèi)停留的時(shí)間，而是以最大時(shí)間而是以最大時(shí)間 （這是已經(jīng)是一個(gè)常數(shù)）（這是已經(jīng)是一個(gè)常數(shù)）代替之，此時(shí)所有代替之，此時(shí)所有 ，這實(shí)際

14、上強(qiáng)化，這實(shí)際上強(qiáng)化了問(wèn)題的要求，即考慮了有些飛機(jī)可能已經(jīng)了問(wèn)題的要求，即考慮了有些飛機(jī)可能已經(jīng)飛出區(qū)域，但仍不允許兩架飛機(jī)的距離小于飛出區(qū)域，但仍不允許兩架飛機(jī)的距離小于8km。 程序：程序：MODEL:TITLE 飛行管理問(wèn)題的非線性規(guī)劃模型飛行管理問(wèn)題的非線性規(guī)劃模型;SETS:Plane/1.6/: x0, y0, cita0, cita1, d_cita;! cita0表示初始角度，表示初始角度，cita1為調(diào)整后的角度為調(diào)整后的角度,d_cita為調(diào)整的為調(diào)整的角度角度;link(plane, plane)|&1 #LT# &2: b,c;ENDSETSDATA:1

15、501402438585236150155220.5145501591301502300052x0 y0 cita0 =max_cita = 30;T_max = 0.283;V=800;ENDDATAINIT:d_cita = 0 0 0 0 0 0;ENDINITfor(plane: cita1 - cita0 = d_cita);for(link(i,j): b(i,j) = -2*(x0(i) -x0(j)*sin (cita1(i)+cita1(j)*3.14159265/360) +2*(y0(i) -y0(j)*cos (cita1(i)+cita1(j)*3.14159265/

16、360); c(i,j) = (x0(i) -x0(j) 2 + (y0(i) -y0(j) 2 - 64;);! 避免碰撞的條件避免碰撞的條件;! 右端點(diǎn)非負(fù)右端點(diǎn)非負(fù);for(link(i,j): Right (2*V*T_max*sin(cita1(i)-cita1(j)*3.14159265/360)2 + b(i,j)*(2*V*T_max*sin(cita1(i)-cita1(j)*3.14159265/360) + c(i,j) 0);! 最小點(diǎn)非負(fù)最小點(diǎn)非負(fù);for(link(i,j): Minimum if(b(i,j)#lt#0 #and# -b(i,j)/4/V/sin(

17、cita1(i)-cita1(j)*3.14159265/360) #gt#0 #and# -b(i,j)/4/V/sin(cita1(i)-cita1(j)*3.14159265/360) #lt#T_max , b(i,j)2-4*c(i,j),-1) 0);!for(link(i,j): if(b(i,j)#lt#0, b(i,j)2-4*c(i,j), -1) ij(2)由圓狀模型導(dǎo)出的方程討論ij的改變量與第i第j兩架飛機(jī)飛行方向角改變量i, j的關(guān)系由題目條件知|vi|=A=800，可用復(fù)數(shù)表示速度設(shè)第i，j飛機(jī)飛行方向改變前的速度分別為改變后的速度分別為改變前后相對(duì)速度分別為兩者

18、之商的幅角就是ij定理：對(duì)第i，第j兩架飛機(jī)，其相對(duì)速度方向ij的改變量ij等于兩飛機(jī)飛行方向角改變量之和的一半，即模型目標(biāo)函數(shù)：Min其中為各飛機(jī)方向角調(diào)整量的最大值或?yàn)榧s束條件：調(diào)整方向角時(shí)不能超過(guò)300：調(diào)整飛行方向后飛機(jī)不能碰撞：模型為化為線性規(guī)劃模型由于i可正可負(fù)，為使各變量均非負(fù)，引入新變量：模型化為模型求解ij的計(jì)算model:sets:plane/1.6/:x0,y0;link(plane,plane):alpha,sin2;endsetsfor(link(i,j)|i#ne#j:sin2(i,j)=64/(x0(i)-x0(j)2+(y0(i)-y0(j)2););for(l

19、ink(i,j)|i#ne#j:(sin(alpha*3.14159265/180.0)2=sin2;);data:x0=150,85,150,145,130,0;y0=140,85,155,50,150,0;enddataendALPHA( 1, 1) 1.234568 ALPHA( 1, 2) 5.391190 ALPHA( 1, 3) 752.2310 ALPHA( 1, 4) 5.091816 ALPHA( 1, 5) 202X.963 ALPHA( 1, 6) 2.234507 ALPHA( 2, 1) 5.391190 ALPHA( 2, 2) 1.234568 ALPHA( 2,

20、 3) 4.804024 ALPHA( 2, 4) 6.613460 ALPHA( 2, 5) 5.807866 ALPHA( 2, 6) 3.815925 ALPHA( 3, 1) 752.2310 ALPHA( 3, 2) 4.804024 ALPHA( 3, 3) 1.234568 ALPHA( 3, 4) 4.364672 ALPHA( 3, 5) 1102.834 ALPHA( 3, 6) 2.125539 ALPHA( 4, 1) 5.091816ALPHA( 4, 2) 6.613460 ALPHA( 4, 3) 4.364672 ALPHA( 4, 4) 1.234568 AL

21、PHA( 4, 5) 4.537692 ALPHA( 4, 6) 2.989819 ALPHA( 5, 1) 202X.963 ALPHA( 5, 2) 5.807866 ALPHA( 5, 3) 1102.834 ALPHA( 5, 4) 4.537692 ALPHA( 5, 5) 1.234568 ALPHA( 5, 6) 2.309841 ALPHA( 6, 1) 2.234507 ALPHA( 6, 2) 3.815925 ALPHA( 6, 3) 2.125539 ALPHA( 6, 4) 2.989819 ALPHA( 6, 5) 2.309841 ALPHA( 6, 6) 1.2

22、34568ijJ=123456i=10.0000005.39119032.2309535.09181620.9633612.23450725.3911900.0000004.8040246.6134605.8078663.815925332.2309534.8040240.0000004.36467222.8336542.12553945.0918166.6134604.3646720.0000004.5376922.989819520.9633615.80786622.8336544.5376920.0000002.30984162.2345073.8159252.1255392.98981

23、92.3098410.000000整理可得ij的值（單位角度）也可以用MATLAB計(jì)算ij的值x=150,85,150,145,130,0;y=140,85,155,50,150,0;k=length(x);alpha=zeros(k);for i=1:k for j=1:k if i=j alpha(i,j)=0; else alpha(i,j)=(180/3.14159265)*asin(8/sqrt(x(i)-x(j)2+(y(i)-y(j)2); end endendalpha計(jì)算ij的值程序?yàn)橛?jì)算結(jié)果為 0 5.391190237223 5.391190237223 032.2309

24、52672331 4.80402393379720.963360893128 5.80786624342132.230952672331 5.091816448550 0 4.364671899111 4.364671899111 0 22.833654204009 4.53769246240220.963360893128 2.234506736995 4.537692462403 2.989819139045 0 2.309841365405 2.309841365405 0ij的計(jì)算:a=150,85,150,145,130,0;b=140,85,155,50,150,0;x=a+b*i

25、;c=243,236,220.5,159,230,52*pi/180;v=exp(i*c);k=length(a);for i=1:k for j=1:k beita(i,j)=(angle(v(i)-v(j)-angle(x(j)-x(i)*180/pi; endendbeita用matlab程序編寫beita = 0 109.2636 -128.2500 24.1798 -186.9349 14.4749 109.2636 0 -88.8711 -42.2436 -92.3048 9.0000 231.7500 271.1289 0 12.4763 301.2138 0.3108 24.1

26、798 -42.2436 12.4763 0 5.9692 -3.5256 173.0651 267.6952 -58.7862 5.9692 0 1.9144 14.4749 9.0000 0.3108 -3.5256 1.9144 0運(yùn)算結(jié)果最優(yōu)解的計(jì)算用LINGO求解程序如下model:sets:plane/1.6/:cita;link(plane,plane):alpha,beta;endsetsmin=sum(plane:abs(cita);for(plane(i): bnd(-30,cita(i),30); );for(link(i,j)|i#ne#j: abs(beta(i,j)

27、+0.5*cita(i)+0.5*cita(j) alpha(i,j); );data:alpha=0.000000,5.391190,32.230953,5.091816,20.963361,2.234507, 5.391190,0.000000,4.8040024,6.813460,5.807866,3.815925,32.230953,4.804024,0.000000,4.364672,22.833654,2.125539,5.091816,6.613460,4.363673,0.000000,4.537692,2.989819, 20.963361,5.807866,22.83365

28、4,4.537692,0.000000,2.309841,2.234507,3.815925,2.125539,2.989819,2.309841,0.000000; beta=0.000000 109.263642 -128.250000 24.179830 173.065051 13.474934109.263642 0.000000 -88.871096 -42.243563 -92.304847 9.000000-128.250000 -88.87096 0.000000 12.476311 -58.786243 0.31080924.179830 -42.243563 12.4763

29、11 0.000000 5.969234 -3.525606174.065051 -92.304846 -58.786244 5.969234 0.000000 1.91438314.474934 9.000000 0.310809 -3.525606 1.913383 0.000000;enddataend用MATLAB計(jì)算編程如下function f,g=plane(x)alph=0.000000,5.391190,32.230953,5.091816,20.963361,2.234507,5.391190,0.000000,4.8040024,6.813460,5.807866,3.81

30、5925,32.230953,4.804024,0.000000,4.364672,22.833654,2.125539,5.091816,6.613460,4.363673,0.000000,4.537692,2.989819,20.963361,5.807866,22.833654,4.537692,0.000000,2.309841,2.234507,3.815925,2.125539,2.989819,2.309841,0.000000; bet=0.000000 109.263642 -128.250000 24.179830 173.065051 13.474934109.2636

31、42 0.000000 -88.871096 -42.243563 -92.304847 9.000000-128.250000 -88.87096 0.000000 12.476311 -58.786243 0.31080924.179830 -42.243563 12.476311 0.000000 5.969234 -3.525606174.065051 -92.304846 -58.786244 5.969234 0.000000 1.91438314.474934 9.000000 0.310809 -3.525606 1.913383 0.000000;f=abs(x(1)+abs

32、(x(2)+abs(x(3)+abs(x(4)+abs(x(5)+abs(x(6);g(1)= alpha(1,2) -abs(beta(1,2)+0.5*x(1)+0.5*x(2);g(2)= alpha(1,3) -abs(beta(1,3)+0.5*x(1)+0.5*x(3);g(3)= alpha(1,4) -abs(beta(1,4)+0.5*x(1)+0.5*x(4);g(4)= alpha(1,5) -abs(beta(1,5)+0.5*x(1)+0.5*x(5);g(5)= alpha(1,6) -abs(beta(1,6)+0.5*x(1)+0.5*x(6);g(6)= al

33、pha(2,3) -abs(beta(2,3)+0.5*x(2)+0.5*x(3);g(7)= alpha(2,4) -abs(beta(2,4)+0.5*x(2)+0.5*x(4);g(8)= alpha(2,5) -abs(beta(2,5)+0.5*x(2)+0.5*x(5);g(9)= alpha(2,6) -abs(beta(2,6)+0.5*x(2)+0.5*x(6);g(10)= alpha(3,4) -abs(beta(3,4)+0.5*x(3)+0.5*x(4);g(11)= alpha(3,5) -abs(beta(3,5)+0.5*x(3)+0.5*x(5);g(12)=

34、 alpha(3,6) -abs(beta(3,6)+0.5*x(3)+0.5*x(6);g(13)= alpha(4,5) -abs(beta(4,5)+0.5*x(4)+0.5*x(5);g(14)= alpha(4,6) -abs(beta(4,6)+0.5*x(4)+0.5*x(6);g(15)= alpha(5,6) -abs(beta(5,6)+0.5*x(5)+0.5*x(6);執(zhí)行程序x0=0,0,0,0,0,0; v1=-30*ones(1,6); v2=30*ones(1,6);opt=;x= constr(plane,x0,opt,v1,v2)結(jié)果：x = -0.0000

35、0576637983 -0.00000576637983 2.58794980234726 -0.00001243487985 0.00003620473095 最優(yōu)解：模型檢驗(yàn)各飛行方向按此方案調(diào)整后，系統(tǒng)各架飛機(jī)均滿足|ij+( I+ j)/2|ij ,結(jié)果是正確的。模型的評(píng)價(jià)與推廣（1）此模型采用圓狀模型分析碰撞問(wèn)題是合理的，同時(shí)采用相對(duì)速度作為判別標(biāo)準(zhǔn)，既體現(xiàn)了碰撞的本質(zhì)（相對(duì)運(yùn)動(dòng)），又簡(jiǎn)化了模型的計(jì)算；（2）建模中用了適當(dāng)?shù)暮?jiǎn)化，將一個(gè)復(fù)雜的非線性規(guī)劃問(wèn)題簡(jiǎn)化為線性規(guī)劃問(wèn)題，既求到合理的解，又提高了運(yùn)算速度，這對(duì)解決高速飛行的飛機(jī)碰撞問(wèn)題是十分重要的。此模型對(duì)題目所提供的例子計(jì)算得出的結(jié)果是令人滿意的。（3）由對(duì)稱性知模型中的約束個(gè)數(shù)是 （n是飛機(jī)數(shù)），所有約束條件數(shù)是 ，計(jì)算量增加不大。2nC2)7(42nnnCn建模方案二：建模方案二：1.問(wèn)題分析目的：不碰撞手段：調(diào)整飛行方向角要求：調(diào)整的幅度盡量小