Change matrix elements in a loop?

2015 11 6
1 个回答
回答已采纳
更新时间：2015 11 6
3 次查看（30 天）

0 个投票

I am trying to figure out how to use a for loop to change the elements in a matrix.
I want to do this ten times.
For example, if row 1 of a 6x3 matrix is (a,a,0), I want to make a loop that generates row 2 as (2a,2a,0). I am basically trying to multiply the initial row to make all 6 rows increase by a multiple of a.
Thank you in advance for any suggestions!

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

 采纳的回答

Star Strider
Star Strider 2015-11-6

0 个投票

No loop necessary. Use bsxfun:
M = randi(9, 6, 3);
Mm = bsxfun(@times, M, [1:size(M,1)]')
This multiplies every row of ‘M’ by the column vector [1 2 ... 6]'.

8 个评论
显示 6更早的评论 隐藏 6更早的评论

lsutiger1
lsutiger1 2015-11-6
Thank you for your response. This doesn't seem to be working, but I think that's because I left out an important piece of information as to why I need a loop: my matrix is actually 10x3x4. I have four different matrices that I need to do this to, and they all have different initial rows.
Star Strider
Star Strider 2015-11-6
My pleasure.
That definitely does change your Question. Assuming that you want to do this for every one of the four ‘pages’ of your matrix, loop through them and do the multiplication on each with the loop:
M = randi(9, 10, 3, 4); % Create Matrix
Mm = zeros(size(M)); % Preallocate
for k1 = 1:size(M,3) % Iterate Over Pages
Mp = M(:,:,k1); % Select Page
Mm(:,:,k1) = bsxfun(@times, Mp, [1:size(Mp,1)]'); % Multiply Page
end
Mp % Last Page (Optional)
Mm(:,:,k1) % Multiplied Last Page (Optional)
The two ‘Optional’ lines simply verify that the code did what we want it to. They can be deleted.
lsutiger1
lsutiger1 2015-11-6
Great, that helps a lot! Thank you! If I may ask one more question.. what would I do if I had a similar 10x3x4 matrix set-up, but it was all zeroes, except for the first row of each 'page.'
For example, if row 1 was (1,0,1), what would I need to do to make row 2 (2,0,2), row 3 (3,0,3), etc.?
I tried to do the following:
if ii=1:n
M = bxsfun(@times, X(:,:,ii), [1:size(X,1)]')
end
Star Strider
Star Strider 2015-11-6
My pleasure!
It should work the same way, without modification. Anything multiplied by zero is zero, so the zeros would remain and only the non-zero values would change.
I assume you define ‘n’ in your loop as size(X,3).
lsutiger1
lsutiger1 2015-11-6
The error that I am getting says:
Error using bsxfun Non-singleton dimensions of the two input arrays must match each other
Here is my entire code thus far; maybe this will help.
% Create matricies for atomic coordinates in unit cell
a = input('What is the lattice constant? \n a = ');
n = input('How many molecular orientations are required? \n n = ');
d = input('How many unit cells in each direction are desired? \n d = ');
type(:,:,n) = zeros(d*a,3);
% Generate initial molecular positions in unit cell
for ii=1:n
disp(['What is the position of molecule ' num2str(ii) ' in terms of a?']);
x = input('x-direction: ');
type(1,1,ii) = x;
y = input('y-direction: ');
type(1,2,ii) = y;
z = input('z-direction: ');
type(1,3,ii) = z;
end
% Generate lattic from initial molecular positions
for ii=1:d % create 10 unit cells in each direction
lattice = bsxfun(@times, type(:,:,ii), [1:size(type,3)]');
end
Star Strider
Star Strider 2015-11-6
You changed my code! You’re not creating the vector along the correct dimension.
Use this instead:
lattice = bsxfun(@times, type(:,:,ii), [1:size(type,1)]');
Use the row size (dimension 1) not the page size (dimension 3). That’s why it’s not working. Changing it back to my original code works.
lsutiger1
lsutiger1 2015-11-6
Yes, I changed it because I was getting that error; neither the 1 or 3 seems to make a difference. This is why I am confused.
Star Strider
Star Strider 2015-11-6
I can’t run your code because I have no idea what you are doing and the inputs should be. Please attach a .mat file with the appropriate size and values for ‘type’. Use the ‘paperclip’ icon, and complete both the ‘Choose file’ and ‘Attach file’ steps.
My code does what you want it to, as you requested, but I can’t troubleshoot it with your code without the appropriate data.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by