以上代码运行有误，请优化，报错信息为：位置 1 处的逻辑索引包含一个在数组边界之外的 true 值。 出错 tianchong (line 21) predicted_values = X(~observed_data_idx, :) * beta;

非常抱歉，给您带来了困扰。错误是由于在构建`predicted_values`时，逻辑索引超出了数组边界导致的。我已经修复了代码，下面是更新后的示例：

% 生成带有缺失值的示例二维数组
rng(1); % 设置随机种子以保持结果的一致性
data = rand(5, 3); % 原始数据
missing_idx = randperm(numel(data), 5); % 生成5个缺失值的索引
data(missing_idx) = NaN; % 将对应索引位置的值设置为NaN

% 创建回归模型并填充缺失值
for j = 1:size(data, 2)
    % 获取当前列的数据和索引
    col_data = data(:, j);
    missing_data_idx = isnan(col_data);
    observed_data_idx = ~missing_data_idx;
    
    % 使用线性回归模型拟合已观测数据
    X = [ones(sum(observed_data_idx), 1), (1:sum(observed_data_idx))'];
    Y = col_data(observed_data_idx);
    beta = X \ Y; % 计算回归系数
    
    % 使用回归模型预测缺失值
    predicted_values = X(missing_data_idx(observed_data_idx), :) * beta;
    
    % 填充缺失值
    col_data(missing_data_idx) = predicted_values;
    
    % 更新原始数据的当前列
    data(:, j) = col_data;
end

disp('填充后的二维数组:');
disp(data);

在修复后的代码中，我更新了`predicted_values`的构建方式，以确保逻辑索引不超出数组边界。我将逻辑索引`missing_data_idx(observed_data_idx)`应用于`X`矩阵，以选择对应的行进行预测。

请注意，这只是一个简单的示例代码，实际的回归填充可能需要更多的数据预处理、模型选择和调参等步骤。此外，根据数据集的特点，可能需要选择其他类型的回归模型。如有其他问题，请随时提问。