MATLAB循环中如何为数组开头补零以构建土壤厚度矩阵？

核心思路是根据当前循环次数，在缩短后的数组开头补对应数量的零，保证每行长度始终等于num_steps。具体修改如下：

修改后的完整循环代码 #

delta_soil_initial = ones(1, num_steps) + delta - 1; % 初始土层厚度
delta_soil_matrix = zeros(num_steps, num_steps); % 初始化存储矩阵

for i = 1:num_steps
    % 更新土层厚度
    delta_soil_initial = delta_soil_initial - f_z .* t_i;
    
    % --- 关键补零操作 ---
    % 计算需要补的零的个数：第i次循环需补i-1个零
    num_zeros = i - 1;
    % 拼接零数组和当前的delta_soil_initial，凑满num_steps长度
    delta_soil_padded = [zeros(1, num_zeros), delta_soil_initial];
    % 存入矩阵对应行
    delta_soil_matrix(i, :) = delta_soil_padded;
    
    % 截断数组，为下一次循环做准备
    f_z(end) = [];
    delta_soil_initial(1) = [];
end

逻辑说明 #

• 第1次循环（i=1）：无需补零，直接存入原长度的数组（刚好等于num_steps）
• 第2次循环（i=2）：补1个零，拼接后长度为1 + (num_steps-1) = num_steps
• 第k次循环（i=k）：补k-1个零，拼接后长度为(k-1) + (num_steps - (k-1)) = num_steps，完美匹配矩阵每行长度

关于circshift无效的原因 #

circshift是循环移位操作，只会调整现有元素的位置，不会新增元素补零，因此无法解决数组长度不足的问题，直接补零拼接才是最适配的方案。

内容的提问来源于stack exchange，提问作者ChemNerd