Principal component analysis

function [finalData]=pca2(t,n,target)%t is the original data set.and n is the number of features you need
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
temp = t;
result = cov(temp);%calculate the covariance matrix of the temp' matrix
[vect value] = eig(result);% calculate the eignvector and the eignvalue
EigArray = eig(result);%calculate the eignvalue
eign=EigArray';
e_length = length(eign);%compute the length of the eignvalue vector

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%the following loops are to sort the eignvector and record the subscript
%corresponding to each eignvalue
for j=1:1:e_length
    val=-inf;
    for i=1:1:e_length
        if val<eign(i)
            val=eign(i);
            e_sub(j)=i;
        end  
    end
    eign(e_sub(j))=-inf;
end
featureVector=[vect(:,e_sub(1))];
% disp('feature is');
% featureVector
for m=2:1:n
    featureVector=[featureVector,vect(:,e_sub(m))];

end
    disp('feature is');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%calculate the meanVector of the temp'
meanVector = temp;
for i=1:1:100
    mean = 1/42*sum(meanVector(:,i));
    for j=1:1:42
        meanVector(j,i)=meanVector(j,i)-mean;
    end  
end
 
finalData=featureVector'*meanVector';%get the final data which is processed use PCA method

end