matlab - 对矩阵的行调用通用函数

Question

我想有效地计算通用核函数的核矩阵MATLAB。这意味着我需要计算 X 的每一行 x 的 k(x,y)以及 Y 的每一行 y。这是一些 matlab 代码，可以计算我想要的内容喜欢，但速度相当慢，

function K=compute_kernel( k_func, X, Y )
    m = size(X,1);
    n = size(Y,1);
    K = zeros(m,n);
    for i = 1:m
        for j = 1:n
            K(i,j) = k_func(X(i,:)', Y(j,:)');
        end
    end
end

有没有其他方法可以解决这个问题，例如一些 bsxfun 变体对从 X 和 Y 获取的每一行调用一个函数？

Answer 1

pdist2(X,Y, dist_func) 不幸的是计算的是 dist_func(X, Y(i,:))，而不是 dist_func(X(i, :)，Y(i,:))。所以我需要的实际功能是，

function K=compute_kernel( k_func, X, Y )
  % Woohoo! Efficient way to compute kernel
  size(X)
  size(Y)
  m = size(X,1);
  n = size(Y,1);
  for i = 1:n
    K(:,i) = pdist2( Y(i,:), X, k_func);
end

它不如仅使用 pdist2 那么好，但仍然比之前的情况高效得多。