matlab二阶迎风差分格式,热传导方程几种差分格式的MATLAB数值解法比较

第2 5卷

第 2期

沈

阳

化

工

大

学

学

报

V0 . 5 N . 1 o2 2

J n. 01l u 2

2 011 0 . 6

J OURNAL HENYANG OF S UNI VERS TY I OF CHEM I CAL TECHNOLoGY

文章编号： 29 0 5—29 ( 0 10 0 7 0 18 2 1 )2— 19— 4

热传导方程几种差分格式的 MA L T AB数值解法比较

冯立伟

(阳化工大学数理系，宁沈阳 10 4 )沈辽 1 12摘要：对求解热传导方程的几种差分格式进行分析， 并讨论使用 MAT A L B编程求解偏微分方

程的方法.编制几种差分格式的 MA L B程序， TA使用算例进行数值实验，不同网格比情况下，在比

较几种算法的优劣. 关键词：热传导方程； MA L B; C akN cl n离散 T A rn— ios o

中图分类号： 0 4 .2 2 18 文献标识码： A

许多工程问题需要研究热量在物体内部的

为空间步长与时间步长， 2族平行直线 x 用 j=

传导情况或某种物质在液体中的扩散情况，因此 j, ht=k矩形区域[,]X[,]分割成矩 r将 0T 0L研究热传导问题特别是非稳态热传导问题十分形网格. 重要 . l热传导方程已有多种求解格式 J E前 . 显式格式： MAT AB是目前最流行、用最广泛的科学和 L应 r工程计算软件 . T A MA L B基于矩阵运算，有强具 +。 大的数值运算能