二维稳态热传导问题的MATLAB数值解法与高斯-塞德尔迭代实现

1. 二维稳态热传导问题背景与数学建模

在工程热物理领域，稳态热传导问题是最基础也是最重要的研究课题之一。当我们面对一块方形金属板，已知其四条边界的温度分布，需要计算板内各点的温度值时，这就构成了典型的二维稳态热传导边值问题。

从物理本质来看，热传导遵循傅里叶定律：热流密度与温度梯度成正比。对于各向同性材料，在无内热源的稳态情况下，温度场满足拉普拉斯方程：

∇²T = ∂²T/∂x² + ∂²T/∂y² = 0

这个偏微分方程(PDE)描述了温度在场域内的分布规律。为了数值求解这个连续问题，我们需要进行离散化处理。最常用的方法是有限差分法，将连续的求解区域划分为规则的网格，用差分近似代替微分。

对于方形区域，采用均匀网格划分，设网格间距为Δx=Δy=h。根据中心差分公式，二阶导数可表示为：

∂²T/∂x² ≈ (T(i+1,j) - 2T(i,j) + T(i-1,j))/h²
∂²T/∂y² ≈ (T(i,j+1) - 2T(i,j) + T(i,j-1))/h²

将它们代入拉普拉斯方程，得到离散后的五点差分格式：

(T(i+1,j) + T(i-1,j) + T(i,j+1) + T(i,j-1) - 4T(i,j))/h² = 0

整理后得到每个内节点的温度计算式：

T(i,j) = (T(i+1,j) + T(i-1,j) + T(i,j+1) + T(i,j-1))/4

这表明在稳态情况下，每个内节点的温度等于其四个邻接节点温度的平均值。这个简洁的表达式为后续的迭代求解奠定了基础。

2. 高斯-塞德尔迭代法原理与实现

2.1 迭代法基本思想

对于大型线性方程组，直接解法如高斯消元法会面临计算量大、存储要求高的问题。而迭代法通过逐步逼近的方式，更适合求解这类稀疏系统。高斯-塞德尔迭代是其中经典且高效的一种方法。

与雅可比迭代不同，高斯-塞德尔迭代在计算过程中会立即使用已更新的节点值。具体到我们的热传导问题，当按行优先顺序遍历网格时，计算T(i,j)时已经获得了T(i-1,j)和T(i,j-1)的新值。这种"即时更新"的策略显著提高了收敛速度。

迭代公式表示为：

Tₖ₊₁(i,j) = (Tₖ(i+1,j) + Tₖ₊₁(i-1,j) + Tₖ(i,j+1) + Tₖ₊₁(i,j-1))/4

其中k表