在尝试理解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个代表,所以我更新了答案
答案 0 :(得分:4)
再次感谢大家,我想我会更新这个,因为我现在了解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];
提供的最终解决方案是任何(bsxfun(@ eq,dxcp,X'),2),它相当于: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.如果上面的行改为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()更快(平均)。
答案 1 :(得分:0)
您不希望在any(x == y)
测试中抛出每个元素组合,您希望测试dxcp
中的每个元素,看它是否存在于X
中。所以这是短版本,也不需要转置。矢量化也应该比bsxfun
快一点。
arrayfun(@(el) any(el == X), dxcp)
结果是
ans =
0 1 0 0 0 0