matlab解偏微分方程

代码为2018年全国数学建模竞赛a题第一问matlab源程序，可动态生成三层隔热服距离与温度的关系图，以及三层隔热服的温度分布图。主要内容建立在一维非稳态热传导以及偏微分方程的解法程序

rho  = [  300   862  74.2  1.18];        % [kg/m^3] 密度**
cp   = [ 1377  2100  1726  1005];        % [J/kg/K] 比热
k    = [0.082 0.370 0.045 0.028];        % [W/m/K ] 热传导率
L    = [  0.6   6.0   3.6   5.0];        % [mm    ] 厚度
xL   = cumsum（L）;                           % [mm    ] 材料层接面位置，即图1中的 [xI xII xIII xIV]
mm2m = 1e-3;                                % 将 mm 转为 m 的系数

dx = 0.2;                                   % [mm    ] 空间网格尺寸
xmax = sum（L）;                              % [mm    ] 模拟区域尺寸
x = 0:dx:xmax;                              % [mm    ] 网格节点坐标
nx = length（x）;                             % 节点数目
ind = round（cumsum（L）/dx）;                  % 不同材料层间的交接点序号

dt = 4e-4;                                  % [s     ] 时间步长
tmax = 90 * 60;                             % [s     ] 最