Tc=qg/ha+Ta; %地表温度 T0f1(1)=Tc; T0f2(1)=Tc; for j=1:num_layers lamdpor(j) = lamds(j) .* (1 - phi(j)) + phi(j) .* lamdaf(j); dtt=dz*qg/lamdpor(j); %竖向网格温差 if(j==1) j00=2; j01=nz/4+1; elseif(j==2) j00=nz/4+2; j01=2*nz/4+1; elseif(j==3) j00=2*nz/4+2; j01=3*nz/4+1; else j00=3*nz/4+2; j01=nz+1; end for i=j00:j01 Tc=Tc+dtt; if(i<=nj) T0f1(i)=Tc; T0f2(i)=Tc; end end end优化代码

可以将外层循环和内层循环合并，使用向量化方式进行优化。具体代码如下：

lamdpor = lamds .* (1 - phi) + phi .* lamdaf; % 计算lamdpor
dzqg = dz * qg; % 计算dz*qg
Tc = cumsum(dzqg ./ lamdpor) + Ta; % 计算温度，加上Ta
T0f1(1:nj) = Tc(j00:j01); % 将温度赋值给T0f1
T0f2(1:nj) = Tc(j00:j01); % 将温度赋值给T0f2

这样就可以将多重循环优化为一次矩阵运算，大大提高了代码执行效率。注意要将地表温度Tc初始化为Ta，然后加上dz*qg/ha，最后再加上Ta。