基于遗传算法的投影寻踪代码，提供大家下载，方便查阅-【研学堂】【代码】投影寻踪代码，请验用！！

function Qa=Project_Pursuit(X,a,Alpha)

　 输入参数列表

　X　　　 本指标矩阵，n×p的矩阵，每一行为一个样本，

　　　　　Xij表示第i行第j列指标，X是否已经归一化均可

　a　　　 投影向量，1×p的矩阵，元素取值范围-1～1，要求其元素平方和等于1

　Alpha　 窗口半径系数，典型取值0.1

　 输出参数列表

　Qa　　　投影指标函数

　 第零步：对a的预处理

b=sqrt(sum(a.^2))

a=a/b

　 第一步：归一化处理

[n,p]=size(X)

x=zeros(n,p)

Xjmax=max(X)

Xjmin=min(X)

for i=1:n

　　x(i,:)=(X(i,:)-Xjmin)./(Xjmax-Xjmin)

end

　 第二步：构造投影指标值

Z=zeros(n,1)

for i=1:n

　　Z(i)=sum(a.*x(i,:))

end

　 第三步：计算投影指标函数

　 计算类间类内矩阵散度

meanZ=mean(Z)

Sa=0

for k=1:n

　　sa=(Z(i)-meanZ(i)).^2

　　Sa=Sa+sa

　　Sa=sqrt(Sa/n)

end

R=Alpha*Sa　窗口半径

Da=0

for k=1:n

　　rik=abs(Z(i)-Z(k))

　　if R>rik

　　　　Da=Da+rik　　　　　　　　　　　　　 Da=Da+R-rik

　　end

end

Qa=Sa*Da