tpca.m

来自「一个关于数据聚类和模式识别的程序,在生物化学,化学中因该都可以用到.希望对大家有」· M 代码 · 共 19 行

clear all

data = rand(30, 300);
dim_n = 15;

tic
[new_data, eigvec, eigval] = pca(data, dim_n);
toc

tic
[new_data1, eigvec1, eigval1] = pca2(data, dim_n);
toc

max(max(abs(eigval-eigval1)))
% Note that eigvec and eigvec1 may have an angle of
% 180 degrees between corresponding eigenvectors.
max(max(abs(abs(eigvec)-abs(eigvec1))))
max(max(abs(abs(new_data)-abs(new_data1))))