When assignment is allowed?
显示 更早的评论
The question concerns the requirement of the sizes of the assignment between two arrays.
To my great surprise this is allowed
lhs = zeros(3,4,5);
rhs = ones(4,5);
lhs(1,:,:) = rhs;
I would thought the last statement would throw an error
"Unable to perform assignment because the size of the left side is 1-by-4-by-5 and the size of the right side is 4-by-5."
but obviously it's not.
Can someone points me to the official document that explains such assigment is allowed and how?
2 个评论
The documentation on how assignment works is nonexistent as far as I can tell.
Will be interested to see if there is a doc-based answer, or if this is a case of "Matlab knows what you want"
Also works for other dimensions, so the leading unity dimension is not some sort of special case
lhs = zeros(3,4,5);
rhs = ones(3,5);
try
lhs(:,2,:) = rhs;
disp('success');
catch
disp('error');
end
the cyclist
2023-8-6
This would have taken me by surprise as well. I did not find this behavior in the documentation, but I did find it in this blog post from @Loren Shure. It is in the section "Retaining Array Shape During Assignment".
采纳的回答
The relevant documentation is at,
Singleton dimensions are ignored when determining if the right and left hand sides have matching dimensions.
x=1:8
x(1:3)=[10;20;30] %assign 3x1 into a 1x3 data region
8 个评论
Paul
2023-8-6
Hi Matt J,
If that doc section answers Bruno's question, it's too bad that such a basic rule on assignment, which is used all day and every day in Matlab, is buried deep in some function related to programming with classes.
Also, I don't see how your example is related. The right hand side of the expression doesn't use indexing, but indexing seems to be the entire point of the linked doc section. I thought your example was more of an historical vagary from first implementation of Matlab.
Interestingly enough, the examples portion of that doc page doesn't include any examples with an indexed variable on the right hand side.
Hi Paul,
The indexing used to generate the right hand side is superfluous, though I agree that may add confusion to the documentation. The indexing of B shown on the right hand side in the documentation is just a way of telling the reader what the dimensions of the assigned data block will be. Bruno's example also doesn't have any indexing on the right hand side.
"Bruno's example ...."
Which was one reason why I wasn't sure that that doc section was on point. But Bruno accepted the Answer so I guess he thinks it is.
As far as assigning 3x1 into a 1x3, I'm still thinking that's some sort of historical special case. In the following code, the 12-element rhs could be strung out into a column vector and assigned into the 12 elements of the lhs in column order, which looks like what happens in the 1-dimensional case. But it doesn't work in two dimensions. I'm not suggesting that it should, but it does look like an inconsistency.
rhs = zeros(3,4);
lhs = ones(4,3);
try
lhs(:,:) = rhs(:,:); % try to assign 12 elements into 12 elements
disp('success');
catch
disp('error');
end
I also don't understand how the subsassgn works in this case, because that page doesn't explain how it handles nonscalar subscripts that are paired up with scalar subscripts.
x = rand(1,3);
y = rand(3,1);
x(1,1:3) = y(1:3,1)
I'm really quite disappointed (and have been for a while) that the doc doesn't have a page dedicated to how assignment works, which needs to cover more than just indexing, e.g., assignment involving different data types.
I also don't understand how the subsassgn works in this case, because that page doesn't explain how it handles nonscalar subscripts
Doesn't it? From the doc: "The order and length of all nonscalar subscripts specified for A matches the order and length of nonscalar subscripts specified for B".
the 12-element rhs could be strung out into a column vector and assigned into the 12 elements of the lhs in column order, which looks like what happens in the 1-dimensional case
No, subsasgn is not stringing anything out. It's just ignoring singleton dimensions.
lhs=nan(4,3);
rhs=rand(1,1,1,1,1,1,1,4,1,3);
lhs(:,:)=rhs
Bruno Luong
2023-8-7
编辑：Bruno Luong
2023-8-7
Bruno Luong
2023-8-7
编辑：Bruno Luong
2023-8-7
Rik
2023-8-7
It really does help for 3D images. The mri.mat example dataset is [row,col,color,page] (even though it is grayscale), but most other places will use [row,col,page,color]. This sytax allow subindexing and swapping conventions with relative easy, without having to use permute or squeeze every second line
Bruno Luong
2023-8-7
编辑：Bruno Luong
2023-8-7
更多回答（0 个）
类别
在 帮助中心 和 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 
