MATLAB热力学方程仿真：应用实例

K=0.1;%定义k温度传导系数

%%定义的温度
T=rand(1,100)*10;

N=200;                  %迭代次数
for i=1:N
    ddT=circshift(T,-1)-T;%边缘假设不变
    ddT=circshift(ddT,-1)-ddT;
    ddT=ddT(1:end-2)
    T(2:end-1)=T(2:end-1)+ddT*K

    c=T

    a=heatmap(T)
    colormap('hot')
    drawnow
    pause(0.3);
    delete(a);
end