如何在Octave中对expand函数的for循环实现矢量化

你原本的循环逻辑是将输入向量统一转为列向量后，逐列重复填充为ncol列的矩阵，这类逐列重复的逻辑完全可以用Octave的原生特性实现矢量化，不需要写显式for循环。

原代码需要修正的冗余内容 #

原函数末尾的f = 1:5属于无关代码，不仅会在调用时额外输出变量内容，还会干扰函数正常逻辑，直接删除即可。

可选矢量化实现方案 #

• 方案1：使用原生repmat函数（兼容所有Octave版本，完全匹配原循环逻辑）
  repmat是Octave专门用于矩阵/向量重复的内置函数，参数格式为repmat(待重复对象, 行重复次数, 列重复次数)：

  function Z = expand(xx,ncol)
      xx = xx(:); % 保留原逻辑：将输入强制转为列向量
      Z = repmat(xx, 1, ncol);
  endfunction
  

  调用expand(1:5, 2)即可得到每行元素相同的2列矩阵，返回结果为：

  ans =
  
     1   1
     2   2
     3   3
     4   4
     5   5
  

• 方案2：利用隐式广播特性（Octave 4.0+ 支持，写法更简洁）
  Octave支持维度匹配时的隐式扩展，不需要显式重复向量即可实现批量赋值：

  function Z = expand(xx,ncol)
      xx = xx(:);
      Z(:, 1:ncol) = xx; % 自动将列向量xx广播填充到所有指定列
  endfunction
  

  这种写法不需要预分配矩阵空间，Octave会自动处理维度匹配，运行效率和原循环实现相比有显著提升。
• 方案3：运算触发广播
  可以通过列向量和行向量的逐元素乘法触发广播，实现相同效果：

  function Z = expand(xx,ncol)
      xx = xx(:);
      Z = xx .* ones(1, ncol);
  endfunction
  

  列向量xx维度为[m, 1]，全1行向量维度为[1, ncol]，逐元素相乘时Octave会自动将两者扩展为[m, ncol]的矩阵，每列内容和xx完全一致。

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