主题:求二维热传导方程的ADI差分格式的MATLAB程序
求二维热传导方程的ADI差分格式的MATLAB程序
差份格式是这样的
u(1:11,1:11,1)=0,r=0.01
(1+2r)u1(i,j)-ru1(i+1,j)-ru1(i-1,j)=(1-2r)u(i,j,k)+ru(i+1,j,k)+ru(i,j,k)
(1+2r)u(i,j,k+1)-ru(i,j,k+1)-r(i,j,k+1)=(1-2r)u1(i,j,)+ru1(i+1,j,)+ru1(i-1,j)
u1(i,j)是k层和k+1层中间的过渡变量
请问这个要怎么用FOR循环实现
差份格式是这样的
u(1:11,1:11,1)=0,r=0.01
(1+2r)u1(i,j)-ru1(i+1,j)-ru1(i-1,j)=(1-2r)u(i,j,k)+ru(i+1,j,k)+ru(i,j,k)
(1+2r)u(i,j,k+1)-ru(i,j,k+1)-r(i,j,k+1)=(1-2r)u1(i,j,)+ru1(i+1,j,)+ru1(i-1,j)
u1(i,j)是k层和k+1层中间的过渡变量
请问这个要怎么用FOR循环实现