chenmh_pi.m

来自「Frequently used algorithms for numerical」· M 代码 · 共 10 行

function index = Chenmh_PI(V_all,U);
T = length(V_all);
[p,r] = size(V_all{1});
index = zeros(T,1);
for i = 1:T
    V = V_all{i};
    index(i) = trace(V'*U(:,r+1:end)*U(:,r+1:end)'*V)/...
               trace(V'*U(:,1:r)*U(:,1:r)'*V);
end