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

%计算A矩阵
A(1,1)=-2*dt/(rho(1)*c(1)*dx(1))*(h(1) lam(1)/dx(1));
A(1,2)=2*dt/(rho(1)*c(1)*dx(1)*dx(1))*lam(1);

for n=2:LEN1-1
A(n,n-1)=F(1);A(n,n)=-2*F(1);A(n,n 1)=F(1);
end
A(LEN1,LEN1-1)=2*D(1)*(lam(1)/dx(1));