spe.m

王小平《遗传算法——理论、应用与软件实现》随书光盘

function [e]=spe(tx,vx,lv)
%
% Calculates the SPE for a given data sets
%

%
% Calculate the principal components
%
%[sc lds ssq res q tsq]=pca(tx,0,0,lv);
[sc lds ssq]=pcomp(tx);
lds=lds(:,1:lv);


%
% Calculate the predicted training data set
%
spt=tx*lds;
vxt=spt/lds;

%
% Calculate the predicted validation data set
%
sp1=vx*lds;
vx1=sp1/lds;

vxa=[vxt ; vx1];
x=[tx ; vx];

[D L]=size(vxa);

for i = 1 : D
	e(i,1)=sum((x(i,:)-vxa(i,:)).^2);
end