飞行管理问题的优化模型.doc

下载文档 降价啦

9
0
约9.11千字
约 12页
2017-07-08 发布于湖北
举报
版权申诉
保障服务

飞行管理问题的优化模型.doc

1、本文档共12页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

飞行管理问题的优化模型

MATLAB: %约束条件函数 function [c,ceq]=aircon(delta) x0=[150, 85, 150, 145, 130, 0]; y0=[140, 85, 155, 50, 150, 0]; alpha0=[243,236,220.5,159,230,52]; %第一组数据; %x0=[60, 70, 80, 50, 40, 0]; %y0=[100, 100, 100, 100, 100, 40]; %alpha0=[270, 270, 270, 270, 270, 0]; %第二组数据; %x0=[60, 60, 60, 60, 60, 0]; %y0=[80, 70, 60, 90, 100, 60]; %alpha0=[180, 180, 180, 180, 180, 0]; %第三组数据; %x0=[0, 55, 90, 40, 80, 0]; %y0=[70, 5, 60, 130, 5, 60]; %alpha0=[0, 90, 180, 270, 180, 0]; alpha0=alpha0*pi/180; v=800; co=cos(alpha0+delta); si=sin(alpha0+delta); for i=2:6 for j=1:i-1 b(i,j)=(x0(i)-x0(j))*(co(i)-co(j))+(y0(i)-y0(j))*(si(i)-si(j)); a(i,j)=v*((co(i)-co(j))^2+(si(i)-si(j))^2); t(i,j)=-b(i,j)/a(i,j); if t(i,j)0 d(i,j)=1000; else d(i,j)=(x0(i)-x0(j)+v*t(i,j)*(co(i)-co(j)))^2+(y0(i)-y0(j)+v*t(i,j)*(si(i)-si(j)))^2; end end end c=64-[d(2,1), d(3,1:2), d(4,1:3), d(5,1:4), d(6,1:5)]; ceq=[ ]; /***************分割符****************************/ %目标函数 function f=airfun(delta) f=delta*delta; %f=sum(abs(delta)); /***************分割符****************************/ %计算程序 clear; clc lb=-pi/18*ones(1,6); ub=pi/18*ones(1,6); options=optimset(LargeScale,off); n=20; fval=100; for i=1:n delta0=10*(rand(1,6)-0.5); [dt,feval]=fmincon(@airfun,delta0,[ ],[ ], [ ], [ ],lb,ub,@aircon,options); if (fevalfval) fval=feval; dl=dt; end end dl=dl*180/pi, f=dl*dl %f=sum(abs(dl)) clc x=[60, 70, 80, 50, 40, 0]; y=[100,100,100,100,100,40]; scatter(x,y,30,r,filled); for i=1:6 text(x(i)+2, y(i), num2str(i)); end c=[5 4 2 1 6 3]; axis([-10, 195, -10, 170]); grid on; hold on; plot([0,160,160,0,0],[0,0,160,160,0],b); zt=[270,270,270,270,270,0]; %zt=[270, 270, 270, 270, 270, 0]+[6.04 5.56 5.08 -6.15 -5.67 4.78]; %zt=[270, 270, 270, 270, 270, 0]+ [-7.7172 5.5435 0.3477 -1.0448 6.0843 -3.4276]; zt1=zt*pi/180; vt=1; dx=vt*cos(zt1); dy=vt*sin(zt1); for vt=1:160 x1=x+dx; y1=y+dy; scatter(x1,y1,11,c