matlab 代码通过有限差分法模拟了二维非稳态导热问题，并将初始时刻、最终时刻的温度分布以及温度场随时间的变化进行了可视化展示
clc;clearall;closeall;xlength=1;%x方向长度ylength=1;%y方向长度a=1e-4;%导温系数nx=15;%x方向网格数目ny=15;%y方向网格数目deltax=xlength/nx;%x方向单个网格长度deltay=ylength/ny;%y方向单个网格长度t=1600;%总时间nt=150;%时间步长deltat=t/nt;%时间步%给网格点编号，定义边