function_u.m

来自「无线传感器网络数学优化模型程序」· M 代码 · 共 17 行

function y=function_u(x,a,k,m)
global dimension
global popsize
T=zeros(dimension,popsize);
A=x;
B=k*(A-a).^m;
B1=k*(-A-a).^m;
C=x>a;
C1=find(C);
D=length(C1);
T(C1(1:D))=B(C1(1:D));
C=x<-a;
C1=find(C);
D=length(C1);
T(C1(1:D))=B1(C1(1:D));
y=sum(T);