2018A高温作业服设计 (显式 隐式差分代码)(9)

T=37*ones(LEN4,1);%初始温度37°
Tend=ones(1,t/dt-1);%Tend 矩阵记录每次迭代皮肤层温度
A_=A^-1;
%开始迭代，T为每次迭代后各微元的温度

for i=1:t/dt-1
T=A_*(T C);
Tend(1,i)=T(end,1);
end
%Tm记录每一秒人体温度
Tm=zeros(t,1);