基于MATLAB的有限差分法解决二位瞬态导热问题

本程序解决的是应对平板的二维导热问题，可以定义边界温度及平板材料，获得一定时间后的温度云图。它只考虑热传导，没有考虑对流传热及热辐射，自己做个记录吧。

程序代码如下：

clc;
clear;
%%
%本文档为有限差分法计算二位导热问题的实例
%为了计算的方便，设定铝板的边界初始为设定初始的值，且温度在传导过程中保持不变
%导热方程为:  T_deta/t_deta=a(d(dT/dx)/dx+d(dT/dy)/dy)
%author:sealingj
%%
%网格划分 %%%%%%%%%%%%%%%%%%
%            TB
%      -----------------
%     |                 |
%  LB |      T(M,N)     | RB    （L_H)M
%     |                 |
%     |_________________|
%             BB（L_W）N          
%%%%%%%%%%%%%%%%%%%%%%%%%%%
M=30;   %划分网格节点
N=30;   %划分网格节点
T(M,N)=0;  %平板初始温度
T(:,:)=30;
X=T; %设置临时温度场分布
K=moviein(200); %获得图像


%%
%初始数据设置
L_H=0.05;        %板的高度
L_W=0.05;        %板的宽度