eg10_4.m

来自「数学建模中的飞行管理问题的matlab程序,分了三个m文件」· M 代码 · 共 18 行

clear;
clc;
deltaini=[0,0,0,0,0,0];
vlb=-pi/6*ones(1,6);
vub=pi/6*ones(1,6);
options=optimset('LargeScale','off');
    [dl,fval]=fmincon(@airfun,deltaini,[],[],[],[],vlb,vub,@airfunco,options);
n=10
for i=1:n-1
    deltaini=[(rand(1,6)-0.5)*0.1];
    [dt,feval]=fmincon(@airfun,deltaini,[],[],[],[],vlb,vub,@airfunco,options);
    if feval<fval
        fval=feval;
        dl=dt;
    end
end
dl=dl*180/pi
fval=fval*(180/pi)^2