您所在位置：主题：求助二维热传导方程的matlab解法
zhaofengyancongsen
[专家分：0] 发布于 2010-06-22 13:13:00
你好：
我现在在编写二维热传导方程，为什么老出现错误啊！请求高人指点。万分感谢！
时间不怎到怎么来处理。程序如下：
求二维热传导方程的ADI差分格式的MATLAB程序
差份格式是这样的
u(1:11,1:11,1)=0,r=0.01
(1+r)*u(i,j，k+1/2)-r/2*u(i+1,j,k+1/2)-r/2*u(i-1,j,k+1/2)=(1-r)*u(i,j,k)+r/2*u(i,j-1,k)+r/2*u(i,j+1,k)
(1+r)*u(i,j,k+1)-r/2*u(i,j-1,k+1)-r/2*(i,j+1,k+1)=(1-r)*u(i,j,k+1/2)+r/2*u(i+1,j,k+1/2)+r/2*u(i-1,j,k+1/2)
u(i,j,k+1/2)是k层和k+1层中间的过渡变量
请问这个要怎么用FOR循环实现?在这个过程中我用的是追赶法解的。
我现在在编写二维热传导方程，为什么老出现错误啊！请求高人指点。万分感谢！
时间不怎到怎么来处理。程序如下：
求二维热传导方程的ADI差分格式的MATLAB程序
差份格式是这样的
u(1:11,1:11,1)=0,r=0.01
(1+r)*u(i,j，k+1/2)-r/2*u(i+1,j,k+1/2)-r/2*u(i-1,j,k+1/2)=(1-r)*u(i,j,k)+r/2*u(i,j-1,k)+r/2*u(i,j+1,k)
(1+r)*u(i,j,k+1)-r/2*u(i,j-1,k+1)-r/2*(i,j+1,k+1)=(1-r)*u(i,j,k+1/2)+r/2*u(i+1,j,k+1/2)+r/2*u(i-1,j,k+1/2)
u(i,j,k+1/2)是k层和k+1层中间的过渡变量
请问这个要怎么用FOR循环实现?在这个过程中我用的是追赶法解的。