MATLAB: How does bsxfun actually work

bsxfunmatrix manipulation

I tried to multiply two matrices using bsxfun. Here is a sample code.
 a = [0,1];
 b = [1,2];
 q = bsxfun(@times, a, b);
Once I run this code, I get the output
 [0 2]
So I tried the following combinations
 q1 = bsxfun(@times, a, b');
 q2 = bsxfun(@times, a', b);
 q3 = bsxfun(@times, a', b');
and all gave me a valid output and none of these generate a multiplication error.
The help says the function uses "singleton expansion" which, I suppose, means it inserts a singleton dimension in either of the two matrices to make the multiplication possible.
Can someone explain what exactly "singleton expansion" means and how is it used in this function?

Best Answer

"Can someone explain what exactly "singleton expansion" means and how is it used in this function?"
bsxfun (and many operators that support this since R2016b) expands any singleton dimension to match the dimension of the other array. For example, lets say you have two arrays with the following sizes:
size(A) = [4,1,5]
size(B) = [4,9,1]
Note how the second dimension of A and the third of B are both equal to 1 (i.e. they are singletons). Then using bsxfun (or the new operators) to perform any binary operation will expand those singleton dimensions to match the size of that dimension of the other array, giving an output with size:
size(out) = [4,9,5]
Thus the singleton dimensions have been expanded. Note that it is called singleton expansion because it only applies to dimensions with size 1, and is NOT (currently) generalized to other non-matching sizes. For example, this would be an error:
size(A) = [3,1,5]
size(B) = [4,9,1]
because the first dimensions are NOT singleton and have different sizes.
Related Question