arrays - 将 bsxfun 与匿名函数一起使用

Question

在尝试理解 bsxfun 函数后，我尝试在脚本中实现它以避免循环。我试图检查数组中的每个单独元素是否包含在一个矩阵中，返回一个与分别包含 1 和 0 的初始数组大小相同的矩阵。我创建的匿名函数是:

myfunction = @(x,y) (sum(any(x == y)));

x 是一个矩阵，它将包含每句话的“接受值”。 y 是输入数组。到目前为止，我已尝试以这种方式使用 bsxfun 函数:

dummyvar = bsxfun(myfunction,dxcp,X)

我知道 myfunction 等于匿名函数的句柄，bsxfun 可以用来完成这个我只是不明白以下错误的原因:

Non-singleton dimensions of the two input arrays must match each other. 

我正在使用以下测试数据:

dxcp = [1 2 3 6 10 20];
X = [2 5 9 18];

希望输出是:

dummyvar = [1,0,0,0]

干杯，NZBRU。

编辑:达到 15 个代表所以我更新了答案

Answer 1

再次感谢大家，我想我会更新这个，因为我现在了解 Divakar 提供的解决方案是如何工作的。这可能会阻止其他阅读过我最初的问题并对 bsxfun() 的工作原理感到困惑的人感到困惑，我认为把它写出来也有助于我更好地理解它。

注意:以下可能不正确，我只是试图通过查看这个案例来了解该功能的运作方式。

bsxfun 函数的输入是 dxcp 和 X 转置。使用的函数句柄是 @eq，因此比较了每个元素。

%%// Given data
dxcp = [1 2 3 6 10 20];
X = [2 5 9 18];

以下代码:

bsxfun(@eq,dxcp,X')

将第一个输入变量 dxcp 的每个值与 X' 的每一行进行比较。下面的矩阵是这个的输出:

dummyvar =

 0     1     0     0     0     0
 0     0     0     0     0     0
 0     0     0     0     0     0
 0     0     0     0     0     0

通过比较1和2找到第一个元素 dxcp = [1 2 3 6 10 20]; X' =​​ [2;5;9;18];

通过比较 2 和 2 dxcp = [1 2 3 6 10 20] 找到第一行的下一个； X' =​​ [2;5;9;18];

重复此操作，直到 dxcp 的所有值都与 X' 的第一行进行比较。按照这个逻辑，第二行中的第一个元素使用以下比较计算:dxcp = [1 2 3 6 10 20]; X' =​​ [2;5;9;18];

提供的最终解决方案是 any(bsxfun(@eq,dxcp,X'),2) 等同于:any(dummyvar,2)。 http://nf.nci.org.au/facilities/software/Matlab/techdoc/ref/any.html似乎很好地详细解释了 any 功能。基本上，说:

A = [1,2;0,0;0,1]

如果运行以下代码:

result = any(A,2)

然后函数 any 将检查每一行是否包含一个或多个非零元素，如果是则返回 1。这个例子的结果是:

result = [1;0;1];

因为第二个输入参数等于 2。如果将上面的行更改为 result = any(A,1) 那么它将检查每一列。

使用这个逻辑，

result = any(A,2)

用于获取最终结果。

1
0
0
0

如果需要可以将其转置为等于

[1,0,0,0]

性能-运行以下代码后:

tic
dummyvar = ~any(bsxfun(@eq,dxcp,X'),2)' 
toc

发现时长为:

Elapsed time is 0.000085 seconds.

下面的替代方案:

tic
arrayfun(@(el) any(el == dxcp),X)
toc

使用 arrayfun() 函数(将函数应用于数组的每个元素)导致运行时间为:

Elapsed time is 0.000260 seconds.

^以上运行时间是每次运行 5 次的平均值，这意味着在这种情况下 bsxfun() 更快(平均)。