MATLAB 索引超出矩阵维度。

开始运行程序没问题，后面在今天用的时候突然提示“索引超出矩阵维度。”我代码和数据集都没换过。。。不知道原因，求大神解一下

function [acc,G,Cls] = GFK(X_src,Y_src,X_tar,Y_tar,dim)
    Ps = pca(X_src);
    Pt = pca(X_tar);
    G = GFK_core([Ps,null(Ps')], Pt(:,1:dim));
    [Cls, acc] = my_kernel_knn(G, X_src, Y_src, X_tar, Y_tar);
end
function G = GFK_core(Q,Pt)

    N = size(Q,2); % 
    dim = size(Pt,2);

    % compute the principal angles
    QPt = Q' * Pt;
    [V1,V2,V,Gam,Sig] = gsvd(QPt(1:dim,:), QPt(dim+1:end,:));
    V2 = -V2;
    theta = real(acos(diag(Gam))); % theta is real in theory. Imaginary part is due to the computation issue.

    % compute the geodesic flow kernel
    eps = 1e-20;
    B1 = 0.5.*diag(1+sin(2*theta)./2./max(theta,eps));
    B2 = 0.5.*diag((-1+cos(2*theta))./2./max(theta,eps));
    B3 = B2;
    B4 = 0.5.*diag(1-sin(2*theta)./2./max(theta,eps));
    G = Q * [V1, zeros(dim,N-dim); zeros(N-dim,dim), V2] ...
        * [B1,B2,zeros(dim,N-2*dim);B3,B4,zeros(dim,N-2*dim);zeros(N-2*dim,N)]...
        * [V1, zeros(dim,N-dim); zeros(N-dim,dim), V2]' * Q';
end

我重新去下了原文件，还是提示索引超出矩阵维度。
此外还进行了断点测试，根本进不了GFK_core的函数里面

刚刚解决了，MATLAB里面我加了drtoolbox的降维包，并且修改了pathdef.m的文件，然后就不能用的，原因尚不清楚

既然是矩阵维度超出，那么问题就在G = GFK_core([Ps,null(Ps')], Pt(:,1:dim));里，
很可能是Pt(:,1:dim)出错，其中dim的值超过了Pt中的最大个数。
可以下一个断点，看看到这一步时Pt的值和dim的值是否有问题，
实在不行就将dim的值设为1，看看这样还会不会报错。