gusucode.com > 很好用的投影寻踪程序 matlab开发源码 > Monkey_PPE/Target.m

    %subfunction of RAGA_PPC
function y=Target(x,a)
[m,n]=size(x);
for i=1:m
    s1=0;
    for j=1:n
       s1=s1+a(j)*x(i,j);
    end
    z(i)=s1;
end
%求z的标准差Sz
Sz=std(z);
%计算z的局部密度Dz
R=0.1*Sz;
s3=0;
for i=1:m
    for j=1:m
        r=abs(z(i)-z(j));
        t=R-r;
        if t>=0
            u=1;
        else
            u=0;
        end
        s3=s3+t*u;
    end
end
Dz=s3;
%计算目标函数值Q
y=Sz*Dz;