Problem with for loop within a loop of a loop

1 次查看（过去 30 天）

Mette Marie Bondegaard Petersen 2019-9-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479333-problem-with-for-loop-within-a-loop-of-a-loop

评论： Mette Marie Bondegaard Petersen 2019-9-8

The following code calculates the excess pore water pressure u after consolidation, and generates a m x n matrix, where the thickness of the soil layer i divided into m equal parts and u is calculated for n time steps.

In this simplified exampel, beta is a constant, but I want to impliment specifik values of beta for different equal sized depths indicated by m.

So say I have a beta-vector containing 5 values and a soil layer divided into m=50 equal parts. beta(1) will run from m(1):m(10), beta(2) will run from m(11):m(20) and so forth.

Does anyone have a suggestion on how to impliment the beta-vector?

beta=0.2; n= 10; m= 5
u_1=[60 54 41 29 19 15]';
u(:,1)=u_1;
for j=(1:n)
  for i=(2:m) 
    u(i,j+1)=u(i,j)+beta*(u(i-1,j)+u(i+1,j)-2*u(i,j));
    for i=(m+1) % As we look at a half-closed layer
      u(i,j+1)=u(i,j)+beta(end)*(2*u(i-1,j)-2*u(i,j))
    end
  end
end

(Ultimatly I will end up with lenght(beta)=30, m= 300 and n=>130 and even with beta as a constant it takes some time to run.)

Much appriciated!

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

Guillaume 2019-9-7

在 MATLAB Online 中打开

Note that

for i = (m+1) %why the () ?

is simply

i = m+1;

Mette Marie Bondegaard Petersen 2019-9-8

Understood :)

This is a result of me trying to figure out how to run the loop for different intervals of m, corresponding to different values of beta, as the actual beta is a vector.

请先登录，再进行评论。

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

采纳的回答

Guillaume 2019-9-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479333-problem-with-for-loop-within-a-loop-of-a-loop#answer_390833

在 MATLAB Online 中打开

Probably, the biggest source of slow down is the lack of preallocation of u. As a result, it grows one column at a time, necessating reallocation and copying on each j step.

u = zeros(m+1, n+1);  %preallocation
u(:, 1) = u1;
for j = 1:n   %why the () ?
    for i = 2:m   %so u(1, j) is 0 for all but j = 1? 
        u(i,j+1)=u(i,j)+beta*(u(i-1,j)+u(i+1,j)-2*u(i,j));
    end
    u(m+1,j+1)=u(m+1,j)+beta*(2*u(m,j)-2*u(m+1,j));  %not sure why there was a beta(end) here.
end

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

Guillaume 2019-9-8

在 MATLAB Online 中打开

I completely missed your question about beta.

The simplest thing is to repelem your different beta across the rows. e.g.

beta = [0.2, 0.3, 0.5, 0.7, 0.9];
ubeta = repelem(beta(:), m/numel(beta));

and then you use ubeta(i) in your equation.

Also, I've just noticed that you don't need the i loop:

for j = 1:n
    u(2:end-1, j+1) = u(2:end, j) + ubeta(2:end) .* (u(1:end-2, j) + u(3:end, j) - 2*u(2:end-1, j));
    u(end, j+1) = u(end, j) + 2*ubeta(end) *(u(end-1, j) - u(end, j));
end