资源描述：
《718几何物理问题作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
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')%蓝色直线为向下飞的飞机航线endendholdonfori=l:length(X)plo
2、t(X(i)zY(i)z^k')%黑色的号为飞机初始位置endLingo:MODEL:TITLE飞行管理问题的非线性规划模型;SETS:Plane/1..6/:xO,yO,citaO,citalzdcita;!citaO^示初始角度，cital为调整后的角度,d_cita为调整的角度;link(plane,plane)
3、&1#LT#&2:b,c;ENDSETSDATA:!输入飞机初始朋标和角度（角度制）；xOy0citaO=1150140243885852361150155220.5114550159113015023000052max_ci=30;T_max=0.283;V=800;EN
4、DDATAINIT:dcita=000000;ENDINIT@for(plane:cital-citaO=d_cita);!模型中二次函数系数；@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;)；!避免碰撞的条件；!右端点非负；@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);!最小点非负；@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));!调整角度上下限，单位为角度；@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:飞行管理问题的非线性规划模型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这个结果得到的是-个局部极小点，调整角度较大•可以看到，在0.01°的误差要求下，需要调整第1,3,4三架飞机的角度，分别调整-2&54°,
8、-3.68°,-30.00°•调整量的平方和为1728.276・MODEL:!飞行管理模型；SETS:Plane/1..6/:x0,y0,d_cita;!d.cita为调整的角度；一link(plane,plane)
9、&1#LT#&2:alpha,beta;EXDSETSDATA:xOyO=15014085851501551455013015000;beta=109.263642-128.25000024.179830-
