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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

那你就应该直接将你的代码贴出来,让大家帮你分析.

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 请问这个差分格式要怎么编

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。