punish.m

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

function y=punish(x,r)
x1=x(1);
x2=x(2);
x3=x(3);
y=1000-x1^2-2*x2^2-x3^2-x1*x2-x1*x3+r*((x1^2+x2^2+x3^2-25)^2+(8*x1+14*x2+7*x3-56)^2);