热传导方程有限差分法实现matlab,热传导方程有限差分法的MATLAB实现-史策

《热传导方程有限差分法的MATLAB实现-史策》由会员分享，可在线阅读，更多相关《热传导方程有限差分法的MATLAB实现-史策(4页珍藏版)》请在读根文库上搜索。

1、理论物理与应用物理学研究 收稿日期 ：2009-04-20 作者简介 ： 史 策 (1986-)，男 ，陕西兴平市人 ，西安建筑科技大学理学院硕士研究生 ，研究方向为微分方程数值解法 。 2009 年 7 月 咸阳师范学院学报 Jul.2009 第 24 卷 第 4 期 Journal of Xianyang Normal University Vol.24 No.4 近些年来 ，求解热传导方程的数值方法 1 取得进 展 ，特别是有限差分区域分解算法 2 ，此类算法的特 点是在内边界处设计不同于整体的格式 ， 将全局的 隐式计算化为局部的分段隐式计算 。 使人从感觉上 认为这样得到的解会比全局。

2、隐式得到的解的精度 差 ，但大量的数值实验表明事实正好相反 ，用区域分 解算法求得的解的精度更好 。 MATLAB 具有强大的图形绘制功能 3 ，为科学 计算和图形处理提供了很大的方便 。 用户只须指定 绘图方式 ,并提供充足的绘图数据 ,用很少的程序指 令就可得到直观 、形象的图形结果 。 因此 ，近些年 来 ， 越来越多的人开始使用 MATLAB 来求解热传 导方程 4，5 。 借助 MATLAB 的数值计算和图形处理 技术 6 ，我们可以绘制出热传导方程数值解的二维 、三 维图形 ，从而可以更好地理解热传导方程解的意义 。 一维热传导方程 坠u 坠t =a 2 坠 2 u 坠x 2 ， 。

3、是最简单的偏微分方程之一 ， 其定解问题的数值解 法主要有有限元法和有限差分法等 ， 对于有限元法 来说 ，适用处理复杂区域 、精度可选 ；缺点在于内存 和计算量巨大 ， 不易编程实现 。 对于有限差分法来 说 ，虽然比较直观 、理论也比较成熟 、精度可选 ；但是 不规则区域处理繁琐 ， 网格生成可以使有限差分方 法 7 (FDM)应用于不规则区域 ，但是对区域的连续性 等要求较严 。 适用 FDM 的好处在易于编程 ，易于并 行 。 鉴于以上情况 ，本文考虑以下边界值问题 ： 坠u 坠t =a 2 坠 2 u 坠x 2 ，0l,t0 u| x=0 =0,u| x=l =0, t0 u| x=。

4、0 =sin x l ! ,0 x # % % $ % % clc; formatshorte a=input(请输入系 数 a的 值 ：); l=input(请输入长 度 l的 值 ：); M=input(请输入将区 间 0,l等分的个 数 M：); ot=input(请输入时间增 量 ot的 值 ：); n=input(请输入运行次 数 n的 值 ：); ox=l/M;x0=zeros(