資源描述：
《718幾何物理問題作業(yè)》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、1.CUMCM-1995A：—個飛行管理問題functionpicture(XfYfa)X=[50851501451300];Y=[14085155501500];a=[243236220.515923052];figureaxis([0,160,0,160])fori=l:length(X)holdonx=0:0.1:160;y=tan(a(i)*pi/180)*(x-X(i))+Y(i);ifa(i)<=180plot(x,yzr')%紅色虛線為向上飛的飛機航線elseplot(x,y,'-b')%藍(lán)色直線為向下飛的飛機航線endendholdonfori=l:length(X)plo
2、t(X(i)zY(i)z^k')%黑色的號為飛機初始位置endLingo:MODEL:TITLE飛行管理問題的非線性規(guī)劃模型;SETS:Plane/1..6/:xO,yO,citaO,citalzdcita;!citaO^示初始角度，cital為調(diào)整后的角度,d_cita為調(diào)整的角度;link(plane,plane)
3、&1#LT#&2:b,c;ENDSETSDATA:!輸入飛機初始朋標(biāo)和角度（角度制）；xOy0citaO=1150140243885852361150155220.5114550159113015023000052max_ci=30;T_max=0.283;V=800;EN
4、DDATAINIT:dcita=000000;ENDINIT@for(plane:cital-citaO=d_cita);!模型中二次函數(shù)系數(shù)；@for(link(i,j):b(izj)=-2*(xO(i)-xO(j))*@sin((cital(i)+cital(j))*3.14159265/360)+2*(y0(i)-yO(j))*@cos((cital(i)+cital(j))*3.14159265/360);c(i,j)=(xO(i)-xO(j))A2+(yO(i)-yO(j))A2-64;)；!避免碰撞的條件；!右端點非負(fù)；@for(link(i,j):[Right](2*V*T_
5、max*@sin((cital(i)-cital(j))*3.14159265/360))人2+b(i,j)*(2*V*T_max*@sin((cital(i)-cital(j))*3.14159265/360))+c(i,j)>0);!最小點非負(fù)；@for(link(i,j):[Minimum]@if(b(i,j)#lt#O#and#-b(i,j)/4/V/@sin((cital(i)-cital(j))*3.14159265/360)#gt#O#and#-b(i,j)/4/V/@sin((cital(i)-cital(j))*3.14159265/360)#lt#T_max,b(i,j
6、)'2-4*c(i,j),-l)<0);!@for(link(i,j):@if(b(i,j)#lt#O,b(i,j)"2-4*c(i,j),-1)<0);@for(link:@free(b));!調(diào)整角度上下限，單位為角度；@for(plane:@bnd(-max_cita,d_cimax_cita));[obj]MIN=@SUM(plane:(dcita)A2);![obj]MIN=@SUM(plane:@abs(d_cita));ENDLocaloptimalsolutionfound.Objectivevalue:Infeasibilities:Totalsolveriteratio
7、ns:1728.2760?3637979E-1133ModelTitle:飛行管理問題的非線性規(guī)劃模型VariableValueReducedCostD_CITA(1)-28.543290.000000D_CITA(2)0.0000000.000000D_CITA(3)-3.6819650.000000D_CITA(4)-30.000004.450516D_CITA(5)0.0000000.000000D_CITA(6)0.0000000.000000這個結(jié)果得到的是-個局部極小點，調(diào)整角度較大?可以看到，在0.01°的誤差要求下，需要調(diào)整第1,3,4三架飛機的角度，分別調(diào)整-2&54°,
8、-3.68°,-30.00°?調(diào)整量的平方和為1728.276?MODEL:!飛行管理模型；SETS:Plane/1..6/:x0,y0,d_cita;!d.cita為調(diào)整的角度；一link(plane,plane)
9、&1#LT#&2:alpha,beta;EXDSETSDATA:xOyO=15014085851501551455013015000;beta=109.263642-128.25000024.179830-
