matlab求一维热传导方程数值解代码,一维热传导方程的数值解

一维热传导方程的数值解

第3卷第3期

2004年8月

淮阴师范学院学报(自然科学版)JOURNALOFHUAIYINTEACHERSCOLLEGE(NATURALSCIENCEEDITION)　Vol13No13May.2004

一维热传导方程的数值解

徐建良,汤炳书

(连云港高等师范专科学校物理系,江苏连云港　222006)

摘　要:利用差分法对数理方程的多个较复杂的一维热传导问题进行分析,给出了直观的图象.

关键词:差分法;一维热传导方程;边界条件

中图分类号:O551.3　　文献标识码:A)0320210205

,有时很难有解析解,即便解出,其解的形式也通常方法.],探讨利用差分法对数理方程进行数值解的

1　第一类边界条件的一维热传导混合问题的解法设含第一类边界条件的一维热传导混合问题为:

2)ut=auxx+f(x,t)　(0

(0≤t≤∞)u(0,t)=μ1(t),　u(l,t)=μ2(t)

u(x,0)=φ(x)　(0≤x≤l)(1)(2)(3)

这是一个11类非齐次边界条件的一维热传导问题,通常这一类混合问题是很难解的,即便解出,其解也通常是一个无穷级数的形式,对该解的物理意义不能直接讨论,不能给出直观的图象.

1.1　计算方法

为求解方程(1),首先定义函数u(x,t)的时间与空间的网格,将x坐标分成N等份,将t坐标分成M等份.令i表示位置x横轴,j表示时间t纵轴.网格上每个格点对应一个温度值.用中心差分近似代替对空间的偏微分,即

2=2Δx2x(4)

用向前差分近似代替对时间的偏微分,即

ui,j+1-ui,j=Δtt(5)

(5)代入(1)式得以(4)、

ui,j+1-ui,j2ui-1,j-2ui,j+ui+1,j=a+f(i,j)ΔtΔx2(6)

解得

ui,j+1=c(ui-1,j+ui+1,j)+(1-2c)ui,j+Δtf(i,j)(7)

其中

Δta2

c=Δx2

收稿日期:2004202221

作者简介:徐建良(19642),男,江苏武进人,讲师,主要从事物理教学研究.(8)