mix.m

来自「使用惩罚函数法对机械设计中的问题进行优化处理」· M 代码 · 共 19 行

x0=[3;4;1];
r0=0.2
c=4;
e=0.1;
[xk,yk]=newtontype(x0,r0);
r=r0*c
[xk1,yk1]=newtontype(xk,r);
e1=abs((yk1-yk)/yk);
e2=sqrt((xk1-xk)'*(xk1-xk));
while(e1>e||e2>e)
    xk=xk1;
    yk=yk1;
    r=r*c
    [xk1,yk1]=newtontype(xk,r); 
    e1=abs((yk1-yk)/yk);
    e2=sqrt((xk1-xk)'*(xk1-xk));
end
xbest=xk1;
ymin=mass(xbest);