qq_43327800
qq_43327800
采纳率14.3%
2019-01-05 18:39 阅读 12.7k

MATLAB 索引超出矩阵维度。

10

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

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条回答 默认 最新

  • 已采纳
    qq_43327800 qq_43327800 2019-01-06 10:29

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

    点赞 2 评论 复制链接分享
  • qq_34540145 qq_34540145 2019-01-06 09:43

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

    点赞 评论 复制链接分享

相关推荐