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;