基于MATLAB编程语言的一维非稳态导热问题计算机仿真程序设计

1. 研究背景，提出问题

1.1 问题背景

一座房屋有一面平整的巨大墙壁由内饰层、结构层、保温层和外饰层共4层构成，自外向内各层分别为A、B、C、D型材料，各种材料的物性参数及厚度见下表。该房屋室内开有空调，控制室温恒定为27℃，内墙对流换热系数固定为25wm·K。白天时，室外温度与室内及墙壁平衡，均为27℃，自下午5点起，室外温度开始下降。房屋外墙直接暴露在大气中，受环境温度影响，认为墙壁外表面与大气温度相同，以每小时0.9℃的速率均匀下降。求该墙壁到晚上10点的温度分布，以及这段时间内墙壁温度的变化情况。

表1墙体建筑材料的物性参数及厚度

材料	导热率 wm·K	密度 kgm3	比热容 Jkg·K	厚度 m
A	1.5	1800	1200	0.050
B	2.0	2000	920	0.100
C	0.8	1500	1500	0.050
D	1.2	1800	1000	0.050

图1 墙体结构示意图

1.2 问题分析

上述问题是一个实际生活场景中的常见问题，为了便于分析研究，可以将实际问题中的一些要素进行抽象化处理。平整巨大的墙壁的温度分布问题可以抽象为无限大平板的一维导热问题；由于室外温度和墙壁外侧温度是随时间变化的，因而是瞬态导热问题；白天时室内外和墙壁的温度平衡，相当于给出了初始条件；墙壁内侧给出了对流换热系数，室外墙壁给出了温度变化规律，即给出了边界条件。因此，该实际问题经过抽象分析可以得出以下要点：

1）这是一个一维平板瞬态导热问题

2）给出了初始条件，墙壁温度均匀为25℃

3）给出了边界条件，墙壁内侧为第三类边界条件，外侧为第一类边界条件

根据以上信息，可以开始建立模型并编写计算机程序来进行计算了。

2.设计思路

对于工程问题的数值计算方法得出结论，要编写计算机程序来进行数值计算，首先要选择一门简单便捷的计算机语言，要求是要计算精确，易于理解，有良好的交互界面，方便注释和后期修改。综合来看，本文选择MATLAB编程语言来进行程序编写。

MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，非常适合学生作为学习计算机编程的良好工具。

程序设计的基本思路

对于一维瞬态导热问题数值计算程序，程序要尽可能的易于理解，代码尽量短小精悍，同时要注意有良好的数据输入、输出接口，便于数据的后处理。此外，在程序编写时，还要注意程序代码的通用性和易读性，以便后期对于代码进行升级，如有需要时，随时可以通过改写升级为二维或多维非稳态的导热或流动的数值分析程序。这就要求在代码编写时注意变量命