Index exceed matrix dimensions.

2 次查看（过去 30 天）
Tarun Pandey
Tarun Pandey 2020-11-20
评论： Tarun Pandey 2020-11-20
% code for PTPB
Cv=0.5; %in m2/day
U0=100; %initial consolidation
U1s=0; %excess pore water pressure at top
U2s=0; %excess pore waler press at bottom
t=30; %in days
z=10; %total thickness in m
n=10; %no.of layers
dz=z/n; %length of interval
dt=1;%length of time interval
lambda= Cv*dt/dz^2;
m=t/dt;
U1=zeros(n,m);
U2=zeros(n,m);
U=zeros(n,m);
U1(2:n-1,1)=100;
U2(2:n-1,1)=100;
j=1;
for k=dt:dt:29*dt
j=j+1;
for i=2:n-1
U1(i,j)=(1-lambda)*U1(i,j-1)/(1+lambda)+(lambda)*U1(i+1,j-1)/(1+lambda)+(lambda)*U1(i-1,j)/(1+lambda);
end
end
for k=dt:dt:29*dt
j=j+1;
for i=n-1:-1:2
U2(i,j)=(1-lambda)*U2(i,j-1)/(1+lambda)+(lambda)*U2(i-1,j-1)/(1+lambda)+(lambda)*U2(i+1,j)/(1+lambda);
end
end
U=0.5*(U1+U2);

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

回答（1 个）

Alan Stevens
Alan Stevens 2020-11-20
Do you just need
j = 1;
before the last
for k=dt:dt:29*dt
loop?
  6 个评论
显示 4更早的评论
Alan Stevens
Alan Stevens 2020-11-20
Or more like this
figure
plot(0:dt:29*dt,U1),grid

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by