Index exceed matrix dimensions.

2020 11 20
1 个回答
更新时间：2020 11 20
2 次查看（30 天）

0 个投票

% 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

0 个投票

Do you just need
j = 1;
before the last
for k=dt:dt:29*dt
loop?

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

Tarun Pandey
Tarun Pandey 2020-11-20
can you help me how to plot this?
Alan Stevens
Alan Stevens 2020-11-20
What exactly do you want to plot?
Tarun Pandey
Tarun Pandey 2020-11-20
U vs t (isochrones)
Alan Stevens
Alan Stevens 2020-11-20
编辑：Alan Stevens 2020-11-20
Do you mean something like this at the end of your code
[x, y] = meshgrid(0:dt:29*dt,1:n);
figure
surf(x,y,U1)
xlabel('x'),ylabel('y'),zlabel('U1')
and similarly for U2.
I should probably have used t instead of x!
Alan Stevens
Alan Stevens 2020-11-20
Or more like this
figure
plot(0:dt:29*dt,U1),grid
Tarun Pandey
Tarun Pandey 2020-11-20
It worked thanku sir..

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by