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的函数里面图片说明图片说明

2个回答

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

qq_43327800
qq_43327800 是因为drtoolbox 里面的PCA函数和MATLAB自带PCA函数冲突了,正解
大约一年之前 回复
qq_43327800
qq_43327800 修改回来就可以用了。。。主要是用了包以后pca返还的不是每个分量的系数,而是变成了降维、、
大约一年之前 回复

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

qq_43327800
qq_43327800 还是感谢了
大约一年之前 回复
qq_43327800
qq_43327800 刚刚解决了,MATLAB里面我加了drtoolbox的降维包,并且修改了pathdef.m的文件,然后就不能用的,原因尚不清楚
大约一年之前 回复
Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
立即提问