call function for each combination of rows in A and columns in B

Question

danielle sisserman 2020-12-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/705373-call-function-for-each-combination-of-rows-in-a-and-columns-in-b

编辑： Bruno Luong 2020-12-30

Xtrain = [1 2; 3 4; 5 6];
Xtest = [7 8 ; 9 10];
kernel = @(x1,x2) norm(x1 - x2);
m = size(Xtrain, 1);
n = size(Xtest,1);
kernels = kernel(Xtest, Xtrain);

I want "kernels" to be the following matrix:

kernels = [ kernel( [7 8] ,[ 1 2] ) kernel([7 8] ,[3 4] ) kernel([7 8] ,[5 6] )

kernel([9 10] , [1 2] ) kernel([9 10] ,[3 4]) kernel( [9 10] ,[5 6] )]

How can I call the kernel function properly so that I will get back this matrix?

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ameer Hamza 2020-12-30

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/705373-call-function-for-each-combination-of-rows-in-a-and-columns-in-b#answer_587343

在 MATLAB Online 中打开

Try this

Xtrain = [1 2; 3 4; 5 6];
Xtest = [7 8 ; 9 10];
kernel = @(x1,x2) norm(x1 - x2);
m = size(Xtrain, 1);
n = size(Xtest, 1);
XtrainC = mat2cell(Xtrain, ones(m,1), 2);
XtestC = mat2cell(Xtest, ones(n,1), 2);
[R, C] = ndgrid(1:n, 1:m);
kernels = arrayfun(@(r, c) kernel(XtestC{r}, XtrainC{c}), R, C)

2 个评论
显示无隐藏无

danielle sisserman 2020-12-30

Beautiful! Thank you so much.

Ameer Hamza 2020-12-30

I am glad to be of help!

请先登录，再进行评论。

Answer 2

Bruno Luong 2020-12-30

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/705373-call-function-for-each-combination-of-rows-in-a-and-columns-in-b#answer_587393

编辑：Bruno Luong 2020-12-30

在 MATLAB Online 中打开

Simple for-loop is 3-4 times faster than fancy MATLAB pseudo vectorization, and much more readable

Xtrain = rand(1000,2);
Xtest = rand(1000,2);
kernel = @(x1,x2) norm(x1 - x2);
m = size(Xtest,1);
n = size(Xtrain,1);
tic
kernels = zeros(m,n);
for i=1:m
    for j=1:n
        kernels(i,j) = kernel(Xtest(i,:),Xtrain(j,:));
    end
end
toc % Elapsed time is 1.389989 seconds.
tic
XtrainC = mat2cell(Xtrain, ones(m,1), 2);
XtestC = mat2cell(Xtest, ones(n,1), 2);
[R, C] = ndgrid(1:n, 1:m);
kernels = arrayfun(@(r, c) kernel(XtestC{r}, XtrainC{c}), R, C);
toc % Elapsed time is 5.821724 seconds.