searching for a better approach for solving a 3-D matrix

Question

Udit Srivastava 2016-10-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/306957-searching-for-a-better-approach-for-solving-a-3-d-matrix

关闭： MATLAB Answer Bot 2021-8-20

在 MATLAB Online 中打开

Hey guys,

I need to solve the following equation.

u(i,j,k+1)=(2*dt/3).*(u(i,j,k).*(u(i+1,j,k)-u(i-1,j,k))./(2*dr)+u(i,j,k-1).*(u(i,j+1,k-1)-u(i,j-1,k-1))./(2*dz));

where dt and dr are some known constants.

The approach I can think of is given below:

for i=2:n, for j=2:n, for k=2:n 
u(i,j,k+1)=(2*dt/3).*(u(i,j,k).*(u(i+1,j,k)-u(i-1,j,k))./(2*dr)+u(i,j,k-1).*(u(i,j+1,k-1)-u(i,j-1,k-1))./(2*dz)); 
end, end, end

Is there any better approach to solve the same equation?

Please help me with this problem and suggest me a better approach, for my approach is not good enough for very long equations including 3 subscripted independent variables i, j, k where each denote r, z and time(t) components respectively, working in cylindrical coordinate system.

Thank you.

Udit Srivastava.

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

此问题已关闭。

Answer 1

KSSV 2016-10-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/306957-searching-for-a-better-approach-for-solving-a-3-d-matrix#answer_238734

在 MATLAB Online 中打开

How about?

i = 2:n ;
j = 2:n ;
for k=1:n
    u(i,j,k+1)=(2*dt/3).*(u(i,j,k).*(u(i+1,j,k)-u(i-1,j,k))./(2*dr)+u(i,j,k-1).*(u(i,j+1,k-1)-u(i,j-1,k-1))./(2*dz)); 
end

2 个评论
显示无隐藏无

Udit Srivastava 2016-10-13

编辑：Udit Srivastava 2016-10-13

Sir, I actually want to remove the 'k' term from my equation and adjust it somewhere else so that the equation works like a 2-D matrix, making it look less cumbersome, like putting an outer loop on 'k' or something like that. Is it anyway possible?