有限差分格式
2019-09-03 10:53:50 0 举报有限差分格式
作者其他创作
大纲/内容
t=t+dtu2(0)=0u2(k)=0
N
call chafenfa
a=1+2*s b(0)=-s/(1+2*s) u(1)=u(1)/a
u1(0)=0u(k)=0
t=t+dt
subroutinechafenfa开始
结束
call chasing
a=1+2*s+s*b(n-3)u(n-1)=(s*u(n-2)+u(n-1))/a
后差法
t=0
abs(t0-t)1e-6
subroutine chasing开始
x=a+i*dxu1(i)=sin(pi*x)+x*(1-x)
a=1+2*s+s*b(i-1) b(i)=-s/au(i+1)=(s*u(i)+u(i+1))/a
end
开始
t=t*0.01
前差法
Y
u1(i)=u1(i)+dt*2
u2(i)=s*u1(i-1)+(1-2*s)*u1(i)+s*u1(i+1)+dt*2
x=a+i*dxu0=exp(-pi*pi*t)*sin(pi*x)+x*(1-x)
u(i)=u(i)-u(i+1)*b(i-1)
0 条评论
下一页
