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

%以下是重复出现的数，先计算方便书写
F=lam./c./rho.*dt./dx./dx;
B=[dx(1)*h(1)/lam(1),dx(4)*h(2)/lam(4)];
D=[dt/(dx(1)*rho(1)*c(1) dx(2)*rho(2)*c(2)),dt/(dx(3)*rho(3)*c(3) dx(2)*rho(2)*c(2)),dt/(dx(3)*rho(3)*c(3) dx(4)*rho(4)*c(4))];

%计算矩阵A
A=zeros(LEN4);
%外界与第一段接触点的位置
A(1,1)=2*F(1)*(1 B(1));A(1,2)=-2*F(1);
for n=2:LEN1-1