一维非稳态导热模型matlab,一维非稳态导热方程求解(附Matlab程序).pdf

使用差分方法求解下面的热传导方程

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 i1, j i , j i1, 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, j1 i, j

t i j k k

上面两式带入原热传导方程

T T T 2T T

i, j1 i, j i1, j i, j i1, j

k h2

2

4 k

令r ，化简上式的

h2

T (12r)T r(T T )

i, j1 i, j i1, j i1, j

如下图：

T

i, j1

rT (12r)T rT

i1, j i, j i1, j

t