马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
大家帮我看看这个程序要怎么改才能实现这个差分格式,刚接触matlab,这个差分格式中的u1(i,j)和是第k层和第k+1层中间的过度变量,一直处理不好,大家帮帮忙啊!
差分格式
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)
u(1:(m+1),1:(n+1),1)=0
r=0.01;
u1=u
for k=2:k
for j=2:m
for i=2:n
u1(i,j)=(r*u1(i+1,j)+r*u1(i-1,j)+r*u(i,j+1,k-1)+(1-2*r)*u(i,j,k-1)+r*u(i,j-1,k-1))/(1+2*r);
u2(i,j,k)=(r*u2(i,j+1,k)+r*u2(i,j-1,k)+r*u1(i+1,j)+r*u1(i-1,j)+(1-2*r)*u1(i,j))/(1+2*r);
end
end
u=u2;
end
[ 本帖最后由 sigma665 于 2007-12-14 12:28 编辑 ] |