How do I update a matrix in for loop (indexing problem)

Question

DM 2015-8-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/234677-how-do-i-update-a-matrix-in-for-loop-indexing-problem

编辑： Guillaume 2015-8-20

Hello,

I have a matrix A and a vector u. I want to run a for loop to update the vector u. The operation is easy, I will multiple A by u (result is another vector denoted by z). Pick the largest element in z (result is scalar denoted by m). Update u by dividing the vector z with m. I don't know how to index the matrix so I used cell notation but that might not be correct.

     A=[1 2 3 4; 5 6 7 8 ; 9 10 11 12; 13 14 15 16];
     u= ones (4,1);
     for s=1:10
         z{s}= A*u{s};
         m(s)= max(z{s});
         u{s} = z{s}/m(s);
    end

Any suggestions would be very helpful.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Guillaume 2015-8-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/234677-how-do-i-update-a-matrix-in-for-loop-indexing-problem#answer_189854

编辑：Guillaume 2015-8-20

在 MATLAB Online 中打开

If you don't want to keep around the result of each step, then your algorithm is simply:

A = [1 2 3 4; 5 6 7 8 ; 9 10 11 12; 13 14 15 16];
u = ones (4,1);
for s = 1:10
   z = A * u;
   m = max(z);
   u = z / m;
end

If you do want to keep the intermediary results, then you can indeed use cell arrays, like this:

    A = [1 2 3 4; 5 6 7 8 ; 9 10 11 12; 13 14 15 16];
    uinitial = ones (4,1);
    numsteps = 10;
    u = cell(1, numsteps+1); z = cell(1, numsteps); m = cell(1, numsteps);
    u{1} = uinitial;
    for s = 1:numsteps
       z{s} = A * u{s};
       m{s} = max(z{s});
       u{s+1} = z{s} / m{s};
    end

Or you could use arrays which will be easier to inspect

    A = [1 2 3 4; 5 6 7 8 ; 9 10 11 12; 13 14 15 16];
    uinitial = ones (4,1);
    numsteps = 10;
    u = zeros(4, numsteps+1); z = zeros(1, numsteps); m = zeros(1, numsteps);
    u(:, 1) = uinitial;
    for s = 1:numsteps
       z(s) = A * u(:, s);
       m(s) = max(z(s));
       u(:, s+1) = z(s) / m(s);
    end