将 bsxfun 与匿名函数一起使用

在尝试理解 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 次，所以我更新了答案

再次感谢大家，我想我会更新这个，因为我现在了解 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 = [ 找到第一个元素12 3 6 10 20]； X' =​​ [2;5;9;18];

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

重复此操作，直到 dxcp 的所有值与 X' 的第一行进行比较。按照此逻辑，第二行中的第一个元素使用以下之间的比较进行计算： dxcp = [12 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 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() 更快（平均）。