function f=constrains(x)
h=[];
n=length(x);
mm=(n-2)/3;
i=1:mm;j=(mm+1):(n-2);
h(i)=x(i);
h(j)=0.0015-abs(x(j));
c=0;
b=0.002:0.00008:0.01;
for w1=b
c=c+1;d=0;
for w2=b
d=d+1;
g(c,d)=P1fun(x,w1,w2);
end
end
h(n-1)=x(n-1)-max(g(:));
cc=0;
bb=0.008:0.0002:0.04;
for w3=bb
cc=cc+1;
gg(cc)=P2fun(x,w3);
end
h(n)=x(n)-max(gg(:));
f=h;
这是我在进行约束优化时的一个约束函数,其中后两个约束我想不到不用循环就能实现的方法.请院长看一下有没有方法化简程序.尽量避免循环出现.(我看了您关于矢量编程的帖子后仍然没有思路)