bsxfun(minus) vs normal minus

5 次查看（过去 30 天）
DIMITRIOS THEODOROPOULOS
回答： Rik 2018-7-3
i have X=eye(3) and A=magic(3) What is the difference between Result1=A-X and the Result2 with this loop
for i=1:3
Result2=bsxfun(@minus,A,X(i,:));
end
  2 个评论
James Tursa
James Tursa 2018-7-3
编辑：James Tursa 2018-7-3
Run it and see. For one, your loop overwrites Result2 with each iteration, so you are not even doing the same calculations and thus you shouldn't expect them to match. And you don't define Y (was this supposed to be X?). What are you really trying to compare?
DIMITRIOS THEODOROPOULOS
DIMITRIOS THEODOROPOULOS 2018-7-3
Ι dont compare anything special.I just try to understand how bsxfun works.To be honest i still dont understand. For example in the first loop what operations are executed??I expect from each element of the first row of A to substract each element of the first row of X.

请先登录，再进行评论。

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

采纳的回答

Rik
Rik 2018-7-3
Have you read the documentation of bsxfun?
Input arrays, specified as scalars, vectors, matrices,
or multidimensional arrays. Inputs A and B must have
compatible sizes. For more information, see Compatible
Array Sizes for Basic Operations. Whenever a dimension
of A or B is singleton (equal to one), bsxfun virtually
replicates the array along that dimension to match the
other array. In the case where a dimension of A or B is
singleton, and the corresponding dimension in the other
array is zero, bsxfun virtually diminishes the
singleton dimension to zero.
This means the following:
A=[1 2];% 1 by 2
B=[3;4];% 2 by 1
%A will be replicated along the first dimension to make it match B
%B will be replicated along the second dimension to make it match A
%then an element-wise operation can be performed:
C1=bsxfun(@minus,A,B);
C2=repmat(A,size(B,1),1)-repmat(B,1,size(A,2));

更多回答（1 个）

Guillaume
Guillaume 2018-7-3
In your loop, which as James pointed out, wouldn't do anything useful since it overwrites Result2 at each step, for each i,
bsxfun(@minus, A, X(i, :))
is exactly equivalent in term of result to
A - repmat(X(i, :), size(A, 1), 1)
but uses much less memory since it doesn't actually replicate the X row.
Note that since R2016b, which introduced implicit expansion, this is also the same as
A - X(i, :)
Before R2016b, the above would have resulted in an error.

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by