MatLab循环执行时无法向round_values列表写入值且无报错如何解决？

问题诱因 #

• 第一层循环迭代范围写法错误：Matlab 中for 变量 = 单个数值的逻辑仅会让循环执行1次，原代码两层循环都没有指定起始值为1，仅写了上限值，导致仅遍历了chl_data的最后一个元素，其余元素完全没有被处理。
• length(chl_data)取值不符合预期：length函数返回的是数组最大维度的长度，如果chl_data是二维矩阵，该值不等于行总数，当行数小于列数时还会出现行索引越界的风险。
• 空矩阵直接索引赋值的隐含问题：初始round_values = []为空矩阵，直接对指定索引位置赋值时，Matlab会自动给未赋值的前置索引位置补0，会引入多余的0值，不符合预期。

修正方案 #

方案1：修复原循环逻辑 #

chl_data = ncread('L3m_20210901-20210930__420535409_4_GSM-MODVIR_CHL1_MO_00.nc', 'CHL1_mean');
TF = isnan(chl_data);
chl_data(TF)=0;
% 提前初始化和chl_data同尺寸的全0矩阵，避免动态扩展性能问题
round_values = zeros(size(chl_data));
% 正确指定循环遍历的行范围：1到总行数
for kx = 1:size(chl_data,1)
    % 正确指定循环遍历的列范围：1到总列数
    for kz = 1:size(chl_data,2)
        if chl_data(kx,kz)~=0
            round_values(kx,kz) = round(chl_data(kx,kz),1);
        end
    end
end

方案2：更高效的向量化实现（无需循环） #

Matlab原生支持矩阵级别的操作，无需遍历即可完成计算，性能远高于嵌套循环：

chl_data = ncread('L3m_20210901-20210930__420535409_4_GSM-MODVIR_CHL1_MO_00.nc', 'CHL1_mean');
% 替换NaN为0
chl_data(isnan(chl_data)) = 0;
% 初始化全0结果矩阵
round_values = zeros(size(chl_data));
% 筛选所有非0位置，直接批量计算保留1位小数
non_zero_idx = chl_data ~= 0;
round_values(non_zero_idx) = round(chl_data(non_zero_idx), 1);

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