数学建模飞行管理问题最简单易懂方法matlab代码-Toy模板网

这篇具有很好参考价值的文章主要介绍了数学建模飞行管理问题最简单易懂方法matlab代码。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

clc
clear
x=[150;85;150;145;130;0];
y=[140;85;155;50;150;0];
z=zeros(6,1);

oo=100;
minn=100;

for o1=243
    z(1)=o1;
    for o2=236
        z(2)=o2;
        for o3=223.5-1:0.01:223.5+1
            z(3)=o3;
            for o4=159
                z(4)=o4;
                for o5=230
                    z(5)=o5;
                    for o6=53-1:0.01:53+1
                        z(6)=o6;k=0;
                        for i=1:5
                            for j=i+1:6

                                v=800;

                                if z(i)<=90
                                    t1=min((160-x(i))/(v*cosd(z(i))),(160-y(i))/(v*sind(z(i))));
                                elseif z(i)>90&&z(i)<=180
                                    t1=min((0-x(i))/(v*cosd(z(i))),(160-y(i))/(v*sind(z(i))));
                                elseif z(i)>180&&z(i)<=270
                                    t1=min((0-x(i))/(v*cosd(z(i))),(0-y(i))/(v*sind(z(i))));
                                else
                                    t1=min((160-x(i))/(v*cosd(z(i))),(0-y(i))/(v*sind(z(i))));
                                end
                                if z(j)<=90
                                    t2=min((160-x(j))/(v*cosd(z(j))),(160-y(j))/(v*sind(z(j))));
                                elseif z(j)>90&&z(j)<=180
                                    t2=min((0-x(j))/(v*cosd(z(j))),(160-y(j))/(v*sind(z(j))));
                                elseif z(j)>180&&z(j)<=270
                                    t2=min((0-x(j))/(v*cosd(z(j))),(0-y(j))/(v*sind(z(j))));
                                else
                                    t2=min((160-x(j))/(v*cosd(z(j))),(0-y(j))/(v*sind(z(j))));
                                end
                                t=min(t1,t2);
                                aa=v^2*((cosd(z(i))-cosd(z(j)))^2+(sind(z(i))-sind(z(j)))^2);
                                bb=v*((x(i)-x(j))*(cosd(z(i))-cosd(z(j)))+(y(i)-y(j))*(sind(z(i))-sind(z(j))));
                                cc=(x(i)-x(j))^2+(y(i)-y(j))^2;
                                if bb>=0
                                    d=cc;
                                elseif -bb>0&&-bb<aa*t
                                    d=cc-bb^2/aa;
                                else
                                    d=aa*t^2+2*bb*t+cc;
                                end

                                if d>64
                                    l=1;
                                else
                                    l=-1;
                                end
                                k=k+l;
                                if k==15
                                    oo=abs(o1-243)+abs(o2-236)+abs(o3-220.5)+abs(o4-159)+abs(o5-230)+abs(o6-52);
                                    k=0;
                                end
                                if k==15
                                    oo=abs(o1-243)+abs(o2-236)+abs(o3-220.5)+abs(o4-159)+abs(o5-230)+abs(o6-52);
                                    k=0;

                                end
                                if(oo<minn)
                                    minn=oo;


                                end
                                if oo<=3.63
                                
                                 disp(o1); disp(o2);disp(o3);disp(o4);disp(o5);disp(o6);
                                end
                            end
                        end
                    end
                end
            end
        end

    end
end

代码中六个角度for循环步长自己改改即可

在约 10,000 米高空的某边长 160 公里的正方形区域内，经常有若干架飞机作水平飞行。

区域内每架飞机的位置和速度向量均由计算机记录其数据以便进行飞行管理。当一架欲

进入该区域的飞机到达区域边缘时，记录其数据后，要立即计算并判断是否会与区域内

的其它飞机发生相撞。如果发生相撞则应计算如何调整各架（包括新进入）的飞机的

飞行方向角，以避免碰撞。现假设条件如下：

1 不相撞的标准为任意两架飞机的距离大于 8 公里

2 飞机飞行方向角调整的幅度不应超过 30 度

3 所有飞机的飞行速度均为每小时 800 公里

4 进入该区域的飞机在到达区域边缘时与区域内飞机的距离应在 60 公里以上

5 最多需考虑 6 架飞机

6 不必考虑飞机离开此区域后的情况。