二维稳态导热微分方程的数值求解matlab程序，包含温度边界、热流边界、对流换热边界三种情况。《传热学》、《数值传热学》、《工程热力学》等课程上级作业。

数值传热的二维稳态导热实验。根据确定的边界条件，利用C语言编写程序对研究区域进行网格划分，获取各网格点的温度值，再利用MATLAB进行数值模拟，确定温度场分布云图

 二维稳态导热的数值计算主要采用了热平衡法。用差分法建立节点的热平衡方程， 将节 点所在的单元体的四个方向传递的热流密度，内热源在单元体产生的热流密度， 根据能量守 恒的原则建立方程，可以得到每一个节点的离散化代数方程。进行数值计算的方法是： 先设定初值， 在根据初值对每一个节点进行迭代可以求得节点 的值。再将初值与新值进行比较，判断迭代的敛散性。比较常用的迭代方法有两种： Gauss-Seidel 法和 Jacobi 法。 Gaus-Seidel 法每次迭代计算， 均是使用节点温度的最新值。 Jacobi 迭代法每次迭代计算均用上一次迭代计算出的值。 对于一个代数方程组， 若选用的迭代方式 不合适有可能导致迭代过程发散， 而对于常物性导热问题组成差分方程组， 每一个方程都选 用导出方程的中心节点温度作为迭代变量则迭代一定收敛。