飞行管理 程序

LLian

谁有95年建模B题飞行管理的相关程序啊？？？急求！！！！！！

LLian

百度上有相关论文，但是没有程序也白搭，哪位大神有啊？？？？

LLian

这个约束条件可以写成什么样的程序呢？急求！！！

LLian

大神们！！！百度里给的那种答案，飞行管理模型中这个怎么求啊？只要给个最终结果就好了！！！

LLian

大神~~~~~~们~~~~~

angler啊

附1 的计算（lingo编程）:
        model:
        sets:
        plane/1..6/:x0,y0;
        link(plane,plane):alpha,sin2;
        endsets
        @for(link(n,m)|n#ne#m:
        sin2(n,m)=64/((x0(n)-x0(m))^2+(y0(n)-y0(m))^2);
        );
        @for(link(n,m)|n#ne#m:
        (@sin(alpha*3.14159265/180.0))^2=sin2;
        );
        data:
        x0=150,85,150,145,130,0;
        y0=140,85,155,50,150,0;
        enddata
        end
       
        ALPHA=
       
        ALPHA( 1, 1)        0.000000
        ALPHA( 1, 2)        5.391190
        ALPHA( 1, 3)        392.2310
        ALPHA( 1, 4)        5.091816
        ALPHA( 1, 5)        159.0366
        ALPHA( 1, 6)        2.234507
        ALPHA( 2, 1)        5.391190
        ALPHA( 2, 2)        0.000000
        ALPHA( 2, 3)        4.804024
        ALPHA( 2, 4)        6.613460
        ALPHA( 2, 5)        5.807866
        ALPHA( 2, 6)        3.815925
        ALPHA( 3, 1)        392.2310
        ALPHA( 3, 2)        4.804024
        ALPHA( 3, 3)        0.000000
        ALPHA( 3, 4)        4.364672
        ALPHA( 3, 5)        202.8337
        ALPHA( 3, 6)        2.125539
        ALPHA( 4, 1)        5.091816
        ALPHA( 4, 2)        6.613460
        ALPHA( 4, 3)        4.364672
        ALPHA( 4, 4)        0.000000
        ALPHA( 4, 5)        4.537692
        ALPHA( 4, 6)        2.989819
        ALPHA( 5, 1)        159.0366
        ALPHA( 5, 2)        5.807866
        ALPHA( 5, 3)        202.8337
        ALPHA( 5, 4)        4.537692
        ALPHA( 5, 5)        0.000000
        ALPHA( 5, 6)        2.309841
        ALPHA( 6, 1)        2.234507
        ALPHA( 6, 2)        3.815925
        ALPHA( 6, 3)        2.125539
        ALPHA( 6, 4)        2.989819
        ALPHA( 6, 5)        2.309841
        ALPHA( 6, 6)        0.000000
       

附2 的计算（matlab编程）:
        a=[150,85,150,145,130,0];b=[140,85,155,50,150,0];
        x=a+b*i;
        c=[243,236,220.5,159,230,52]*pi/180;
        v=exp(i*c);
        k=length(a);
        for n=1:k
           for m=1:k
              cita(n,m)=(angle(v(n)-v(m))-angle(x(m)-x(n)))*180/pi;
           end
        end
        cita
       
        cita =
       
                 0  109.2636 -128.2500   24.1798 -186.9349   14.4749
          109.2636         0  -88.8711  -42.2436  -92.3048    9.0000
          231.7500  271.1289         0   12.4763  301.2138    0.3108
           24.1798  -42.2436   12.4763         0    5.9692   -3.5256
          173.0651  267.6952  -58.7862    5.9692         0    1.9144
           14.4749    9.0000    0.3108   -3.5256    1.9144         0
       
附3 最优解的计算lingo程序：
        因与在lingo编程里不易区分，因此将表示为beita
        model:
        sets:
        plane/1..6/:cita;
        link(plane,plane):alpha,beta;
        endsets
        min=@sum(planeabs(cita));
        @for(plane(n):
                    @bnd(-30,cita(n),30);
                     );
        @for(link(n,m)|n#ne#m:
                      @abs(beta(n,m)+0.5*cita(n)+0.5*cita(m))
              >alpha(n,m);
            );
        data:
        alpha=0.000000,5.391190,32.230953,5.091816,20.963361,2.234507,
        5.391190,0.000000,4.8040024,6.813460,5.807866,3.815925,
        32.230953,4.804024,0.000000,4.364672,22.833654,2.125539,
        5.091816,6.613460,4.363673,0.000000,4.537692,2.989819,      20.963361,5.807866,22.833654,4.537692,0.000000,2.309841,
        2.234507,3.815925,2.125539,2.989819,2.309841,0.000000;
        beta=0.000000 109.263642 -128.250000 24.179830 173.065051 13.474934
        109.263642 0.000000 -88.871096 -42.243563 -92.304847 9.000000
        -128.250000 -88.87096 0.000000 12.476311 -58.786243 0.310809
        24.179830 -42.243563 12.476311 0.000000 5.969234 -3.525606
        174.065051 -92.304846 -58.786244 5.969234 0.000000 1.914383
        14.474934 9.000000 0.310809 -3.525606 1.913383 0.000000;
        enddata
        end
        结果
                               CITA( 1)        0.000000           
                               CITA( 2)        0.000000            
                               CITA( 3)        2.836544         
                               CITA( 4)        0.000000            
                               CITA( 5)        0.000000            
                       CITA( 6)       0.7929160   

LLian

angler啊 发表于 2011-8-9 19:48
附1 的计算（lingo编程）:
        model:
        sets:

大神，你太强悍了！！！

LLian

angler啊 发表于 2011-8-9 19:48
附1 的计算（lingo编程）:
        model:
        sets:

大神，还在不？我想请教一下，任意两架飞机相对飞行速度与飞机质点连线夹角有程序么？有结果么？？？还有刚才程序都是求什么的？能简单说明下么？？

angler啊

上面程序的ALPHA，你自己看下程序。只是给你提示，每个人设的都不一样，你要自己编自己的公式呢。

LLian

angler啊 发表于 2011-8-9 22:39
上面程序的ALPHA，你自己看下程序。只是给你提示，每个人设的都不一样，你要自己编自己的公式呢。

O(∩_∩)O谢谢~~~