如何在Julia中生成任意指定秩的随机矩阵

Julia标准库没有直接提供生成任意指定秩随机矩阵的封装函数，但基于基础线性代数规则，几行代码就能实现需求，不需要额外安装第三方包。

核心逻辑 #

任意一个尺寸为m×n、秩为k的矩阵，都可以分解为两个满秩矩阵的乘积：一个m×k维的列满秩矩阵A，和一个k×n维的行满秩矩阵B，二者相乘得到的M = A*B秩必然为k。只要A、B的元素从连续随机分布中采样，生成过程中出现秩亏的概率为0，结果完全符合要求。

快速示例：生成3×3秩为2的随机矩阵 #

直接用内置的随机数生成函数即可实现：

using LinearAlgebra

# 生成两个满秩因子矩阵
A = randn(3, 2)
B = randn(2, 3)
# 乘积即为目标秩2矩阵
M = A * B

# 校验秩，结果恒为2
rank(M)

通用封装 #

你可以把逻辑封装成通用函数，支持任意尺寸、任意合法秩的随机矩阵生成，默认生成标准正态分布元素的矩阵，不需要额外依赖第三方包：

using LinearAlgebra, Random

function rand_rank(m::Int, n::Int, k::Int; rng=GLOBAL_RNG, randfn=randn)
    @assert 0 ≤ k ≤ min(m, n) "指定的秩k必须满足0 ≤ k ≤ min(矩阵行数, 矩阵列数)"
    k == 0 && return zeros(m, n)
    A = randfn(rng, m, k)
    B = randfn(rng, k, n)
    return A * B
end

调用示例 #

# 对应需求：生成3×3秩为2的随机矩阵（标准正态分布元素）
M2 = rand_rank(3, 3, 2)
rank(M2) # 返回2

# 生成5×7、秩为4的[0,1]均匀分布随机矩阵
M4 = rand_rank(5, 7, 4, randfn=rand)
rank(M4) # 返回4

如果需要生成整数元素、稀疏结构的指定秩随机矩阵，只需要调整传入的randfn采样函数即可，核心乘积逻辑不变。

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