718幾何物理問(wèn)題作業(yè)

《718幾何物理問(wèn)題作業(yè)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、1.CUMCM-1995A：—個(gè)飛行管理問(wèn)題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')%紅色虛線為向上飛的飛機(jī)航線elseplot(x,y,'-b')%藍(lán)色直線為向下飛的飛機(jī)航線endendholdonfori=l:length(X)plo

2、t(X(i)zY(i)z^k')%黑色的號(hào)為飛機(jī)初始位置endLingo:MODEL:TITLE飛行管理問(wèn)題的非線性規(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:!輸入飛機(jī)初始朋標(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;)；!避免碰撞的條件；!右端點(diǎn)非負(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);!最小點(diǎn)非負(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:飛行管理問(wèn)題的非線性規(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這個(gè)結(jié)果得到的是-個(gè)局部極小點(diǎn)，調(diào)整角度較大?可以看到，在0.01°的誤差要求下，需要調(diào)整第1,3,4三架飛機(jī)的角度，分別調(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-