compare_int.m

来自「复合形法matlab优化程序」· M 代码 · 共 39 行

function xh=compare_int(x0)
%x0(:,1)=[0.25 0.5]';x0(:,2)=[0 1]';
%x0(:,3)=[1 0]';x0(:,4)=[0.48 0.55]';
%以矩阵形式输入初始四个顶点x0(2,4)
%%找出最好点--最差点
n=2;k=1;
syms x1 x2;
f=x1^2+2*x2^2-2*x1^2*x2^2;%原函数
for i=1:4
f0(i)=subs(f,{x1,x2},x0(:,i));
end
a=0;
for i=1:4
    if a<f0(i)
        a=f0(i);
        xh(:,1)=x0(:,i);
    end
end
b=0;
for i=1:4
    if b<f0(i)&f0(i)<a;
        b=f0(i);
        xh(:,2)=x0(:,i);
    end
end
c=0;
for i=1:4
    if c<f0(i)&f0(i)<b;
        c=f0(i);
        xh(:,3)=x0(:,i);
    end
end
d=0;
for i=1:4
    if d<f0(i)&f0(i)<c;
        d=f0(i);
        xh(:,4)=x0(:,i);
    end
end