如何在Julia中生成含orthonormal columns的任意维度随机矩阵？

当然有，最常用且高效的方法是通过QR分解实现，步骤如下：

• 首先导入Julia的线性代数工具包：using LinearAlgebra
• 生成指定维度的随机高斯矩阵（要求行数m≥列数n，保证能生成n个线性无关的列向量）：randn(m, n)
• 对该矩阵做QR分解，提取结果中的Q矩阵，它就是列正交归一的目标矩阵

示例代码 #

using LinearAlgebra

# 生成5行3列的列正交归一随机矩阵
m, n = 5, 3
Q = qr(randn(m, n)).Q

# 验证正交性：Q'Q 应近似于3×3单位矩阵
println(Q'Q)

补充说明 #

• 如果需要完整的m×m正交方阵（而非仅前n列），可以在QR分解时加上full=true参数：qr(randn(m, n), full=true).Q，再截取前n列即可
• QR分解是数值稳定的方法，比手动实现Gram-Schmidt正交化更可靠，适合生成大规模矩阵

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