Description: 【研学堂】【代码】投影寻踪代码，请验用！！
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
To Search:
File list (Check if you may need any files):
matlab.m
Download