Running for loop on an equation

Question

aliza mustafa 2022-10-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1815590-running-for-loop-on-an-equation

评论： aliza mustafa 2022-10-1

采纳的回答： Torsten

Hi everyone,

I have the equation:

(1)

I have:

and

I have to compute

for 0, 1, 2, ... 10.

I did it this way in MATLAB:

S = [1/2, 1/3, 1/2, 1/3, 1/3, 1/2;

1/2, 0, 0, 0, 0, 0;

0, 1/3, 1/2, 0, 0, 0;

0, 0, 0, 0, 1/3, 0;

0, 1/3, 0, 1/3, 0, 0;

0, 0, 0, 1/3, 1/3, 1/2];

u_0 = [0, 0, 1, 0, 0 ,0]';

for i = 0:10

u_k = S.*u_i;

I am stuck here. How can I write the equation (1) in matlab? Any help will be highly appreciated. Thanks alot in advance.

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

请先登录，再进行评论。

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

Answer 1

Torsten 2022-10-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1815590-running-for-loop-on-an-equation#answer_1064675

编辑：Torsten 2022-10-1

在 MATLAB Online 中打开

S = [1/2, 1/3, 1/2, 1/3, 1/3, 1/2;
    1/2, 0, 0, 0, 0, 0;
    0, 1/3, 1/2, 0, 0, 0;
    0, 0, 0, 0, 1/3, 0;
    0, 1/3, 0, 1/3, 0, 0;
    0, 0, 0, 1/3, 1/3, 1/2];
u_0 = [0, 0, 1, 0, 0 ,0]';
u = zeros(6,10);
u(:,1) = u_0;
for i = 2:11
    u(:,i) = S*u(:,i-1);
end
u
u = 6×11
         0    0.5000    0.5000    0.4583    0.4444    0.4433    0.4441    0.4444    0.4445    0.4445    0.4444
         0         0    0.2500    0.2500    0.2292    0.2222    0.2216    0.2220    0.2222    0.2222    0.2222
    1.0000    0.5000    0.2500    0.2083    0.1875    0.1701    0.1591    0.1535    0.1507    0.1494    0.1488
         0         0         0         0    0.0278    0.0278    0.0285    0.0278    0.0278    0.0278    0.0278
         0         0         0    0.0833    0.0833    0.0856    0.0833    0.0834    0.0833    0.0833    0.0833
         0         0         0         0    0.0278    0.0509    0.0633    0.0689    0.0715    0.0728    0.0734
%u./vecnorm(u)
%null(S-eye(6))

3 个评论
显示 1更早的评论隐藏 1更早的评论

Torsten 2022-10-1

To run it for 1,2,...,10, you must run the loop from 2 to 11. Think about it.

I corrected it above in my code.

aliza mustafa 2022-10-1

Thank you so much for your guidance.

请先登录，再进行评论。

Answer 2

Davide Masiello 2022-10-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1815590-running-for-loop-on-an-equation#answer_1064685

编辑：Davide Masiello 2022-10-1

在 MATLAB Online 中打开

I suspect you want to do something like this

S   = [ 1/2, 1/3, 1/2, 1/3, 1/3, 1/2;...
        1/2, 0, 0, 0, 0, 0;...
        0, 1/3, 1/2, 0, 0, 0;...
        0, 0, 0, 0, 1/3, 0;...
        0, 1/3, 0, 1/3, 0, 0;...
        0, 0, 0, 1/3, 1/3, 1/2;...
        ];
u       = zeros(6,11);
u(:,1)  = [0, 0, 1, 0, 0 ,0]';
for col = 2:11
    u(:,col) = S*u(:,col-1);
end
u
u = 6×11
         0    0.5000    0.5000    0.4583    0.4444    0.4433    0.4441    0.4444    0.4445    0.4445    0.4444
         0         0    0.2500    0.2500    0.2292    0.2222    0.2216    0.2220    0.2222    0.2222    0.2222
    1.0000    0.5000    0.2500    0.2083    0.1875    0.1701    0.1591    0.1535    0.1507    0.1494    0.1488
         0         0         0         0    0.0278    0.0278    0.0285    0.0278    0.0278    0.0278    0.0278
         0         0         0    0.0833    0.0833    0.0856    0.0833    0.0834    0.0833    0.0833    0.0833
         0         0         0         0    0.0278    0.0509    0.0633    0.0689    0.0715    0.0728    0.0734