热传导方程
热传导方程简介
热传导方程（Heat Equation）是描述热量在介质中传递的偏微分方程，其基本形式为：
[ \frac{\partial u}{\partial t} = \alpha \nabla^2 u ]
其中 ( u(x,t) ) 表示温度分布，( \alpha ) 为热扩散系数，( \nabla^2 ) 是拉普拉斯算子。
有限差分法求解
在MATLAB中，可通过有限差分法（FDM）离散求解热传导方程。以一维情况为例：
离散化方程
空间离散：( x_i = i\Delta x )，时间离散：( t_n = n\Delta t )。
采用显式欧拉格式：
[ u_i^{n+1} = u_i^n + \frac{\alpha \Delta t}{(\Delta x)^2} (u_{i+1}^n - 2u_i^n + u_{i-1}^n) ]
MATLAB代码实现
% 参数设置L=1;% 空间长度T=0.1;% 总时间alpha=0.01;% 热扩散系数Nx=50;% 空间网格数Nt=1000;% 时间步数dx=L/Nx;% 空间步长dt=T/Nt;% 时间步长r=alpha*dt/dx^2;% 初始条件（如高斯脉冲）x=linspace(0,L,Nx+1);u0=exp(-100*(x-0.5).^2)';% 显式迭代求解u=u0;forn=1:Nt u_new=u;fori=2:Nxu_new(i)=u(i)+r*(u(i+1)-2*u(i)+u(i-1));endu=u_new;end% 可视化plot(x,u0,'r--',x,u,'b-');legend('Initial','Final');xlabel('x');ylabel('Temperature');使用PDE Toolbox求解
MATLAB的PDE Toolbox提供了更便捷的偏微分方程求解功能：
- 定义几何模型
model=createpde();geometryFromEdges(model,@squareg);- 设置方程参数
specifyCoefficients(model,'m',0,'d',1,'c',alpha,'a',0,'f',0);- 初始条件与边界条件
setInitialConditions(model,u0);% u0为初始温度分布applyBoundaryCondition(model,'dirichlet','Edge',1:4,'u',0);- 求解与可视化
result=solvepde(model,T);pdeplot(model,'XYData',result.NodalSolution);隐式方法（Crank-Nicolson）
为提高稳定性，可采用隐式Crank-Nicolson格式：
离散方程
[ -\frac{r}{2}u_{i-1}^{n+1} + (1+r)u_i^{n+1} - \frac{r}{2}u_{i+1}^{n+1} = \frac{r}{2}u_{i-1}^n + (1-r)u_i^n + \frac{r}{2}u_{i+1}^n ]
MATLAB实现
需构建三对角矩阵并调用tridiag求解器，或直接使用稀疏矩阵求解：
A=gallery('tridiag',Nx-1,-r/2,1+r,-r/2);B=gallery('tridiag',Nx-1,r/2,1-r,r/2);forn=1:Ntu(2:Nx)=A\(B*u(2:Nx));end注意事项
- 稳定性条件：显式方法需满足 ( \alpha \Delta t / (\Delta x)^2 \leq 0.5 )。
- 边界处理：Dirichlet边界直接赋值，Neumann边界需特殊处理。
- 高维扩展：二维/三维情况下需修改离散格式为对应的五点或七点差分格式。