在匿名函数中使用bsxfun

Question

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

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

X是将包含“接受值”的矩阵。Y是输入数组。到目前为止，我已经尝试以这种方式使用bsxfun函数：

dummyvar = bsxfun(myfunction,dxcp,X)

我知道我的函数等于匿名函数的句柄，可以使用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个rep，所以我更新了答案

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 2 0；X‘= 2;5;9;18；

比较2和2 dxcp =1 2 3 6 10 2 0；X‘= 2;5;9;18；

这是重复的，直到dxcp的所有值与X‘的第一行相比较。按照这种逻辑，第二行中的第一个元素使用以下比较进行计算: dxcp = 1 2 3 6 10 20；X‘= 2;5;9;18；

最后提供的解决方案是any(@eq，dxcp，X‘)，它相当于: any(dummyvar,2)。http://nf.nci.org.au/facilities/software/Matlab/techdoc/ref/any.html似乎很好地解释了任何函数。基本上，说：

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

如果运行以下代码：

result = any(A,2)

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

result = [1;0;1];

因为第二个输入参数等于2。如果上面的行被更改为结果= 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()更快(平均)。

Answer 2

您不希望每个元素组合都抛到您的any(x == y)测试中，而是希望测试dxcp中的每个元素，看看它是否存在于X中。因此，这是短版本，它也不需要转座子。矢量化也应该比bsxfun快一点。

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

结果是

ans =

     0     1     0     0     0     0