使用差分方法求解下面的热传导方程
2
T (x, t) a T (x, t) (0 x 1, 0 t 0.2, a 1)
t xx
初值条件:T ( x,0) 4 x 4 x2 ;
T (0, t) 0
边值条件: ;
T (1, t) 0
使用差分公式
T ( x h, t ) 2T ( x , t ) T ( x h, t ) T 2T T
i j i j i j 2 i1, j i , j i1, j
T ( x , t ) O(h )
xx i j h2 h2
T ( x , t k) T ( x , t ) T T
T ( x , t ) i j i j O(k) i, j1 i, j
t i j k k
上面两式带入原热传导方程
T T T 2T T
i, j1 i, j i1, j i, j i1, j
k h2
2
4 k
令r ,化简上式的
h2
T (12r)T r(T T )
i, j1 i, j i1, j i1, j
如下图:
T
i, j1
rT (12r)T rT
i1, j i, j i1, j
t