jaykangta 发表于 2007-12-6 17:18

请问这个差分格式要怎么编

u(1:11,1:11,1)=0,r=0.01,k=1:300
(1+2r)u'(i,j,k+1/2)-ru'(i+1,j,k+1/2)-ru'(i-1,j,k+1/2)=(1-2r)u(i,j,k)+ru(i,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)u'(i,j,k+1/2)+ru'(i+1,j,k+1/2)+ru'(i-1,j,k+1/2)

xjzuo 发表于 2007-12-6 17:23

请将问题讲清楚一些，例如贴一下公式。
另：作for循环应该就行了．

jaykangta 发表于 2007-12-6 17:38

方程是
Ut=Uxx+Uyy
U(x,y,0)=0
差分格式是ADI交替隐式
u(1:11,1:11,1)=0,r=0.01,k=1:300
(1+2r)u'(i,j,k+1/2)-ru'(i+1,j,k+1/2)-ru'(i-1,j,k+1/2)=(1-2r)u(i,j,k)+ru(i,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)u'(i,j,k+1/2)+ru'(i+1,j,k+1/2)+ru'(i-1,j,k+1/2)
我用FOR循环编了一下,可是它根本不循环

xjzuo 发表于 2007-12-6 17:43

那你就应该直接将你的代码贴出来,让大家帮你分析.
页: [1]
查看完整版本: 请问这个差分格式要怎么编