回 帖 发 新 帖 刷新版面

主题:求二维热传导方程的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循环实现

回复列表 (共1个回复)

沙发

倒数第3行错了,应该是这样的

(1+2r)u(i,j,k+1)-ru(i,j+1,k+1)-r(i,j-1,k+1)=(1-2r)u1(i,j,)+ru1(i+1,j,)+ru1(i-1,j)

我来回复

您尚未登录,请登录后再回复。点此登录或注册