谁懂kpca,有问题请教以下,谢谢!

clear all
X=randn(15,10);
Xmean=mean(X);
Xdelta=std(X);
a=ones(15,1);
x=(X-a*Xmean)./(a*Xdelta);
for i=1:15,
for j=i:15,
k(i,j) = exp(-norm(x(i, -x(j, )^2/2);
k(j,i)=k(i,j);
end
end
I=ones(15,15)/15;
k1=k-I*k-k*I+I*k*I;
[evecs,evals]=eig(k1);
v(:,i)=v(:,i)/sqrt(d(i,i));
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
Y=randn(15,10);
y=(Y-a*Xmean)./(a*Xdelta);
for m=1:15,
for n=1:15,
k2(i,j)=exp(-(norm(y(m, -x(n, )).^2/2);
end
end
k3=k2-I*k2-k2*I+I*k2*I;
for m=1:15,
t(:,m)=k3*v(:,m);
end
tsquare(m)=t(m,6:15)*inv(d(6:15,6:15))*t(m,6:15)"; 数据挖掘研究院
m=1:15;
tsquare;
plot(m,tsquare(m))% T-2统计图

不知道，你什么 意思

楼主要问什么问题呀

Create By Any-Extract(WL-AE)