如何在Julia中生成元素为模1复数的任意大小随机数组？

模为1的复数可以表示为 ( e^{i\theta} )（其中( \theta ) 是[0, 2π)范围内的随机角度），基于这个数学性质，Julia里有几种简洁的实现方式：

1. 直接广播生成（推荐，效率最高） #

利用随机角度直接转换为复数，这是最直接高效的方式：

# 生成3×5的模为1的随机复数数组
m, n = 3, 5
unit_complex_arr = exp.(im .* 2π .* rand(m, n))

• rand(m, n) 生成[0,1)区间的均匀随机数矩阵
• 乘以2π将范围映射到[0, 2π)，得到随机角度
• 通过exp(im * θ)转换为模为1的复数，广播操作.确保对数组每个元素生效

你也可以封装成复用函数：

function rand_unit_complex(dims...)
    exp.(im .* 2π .* rand(dims...))
end

# 调用示例：生成4维数组
arr = rand_unit_complex(2, 2, 2, 2)

2. 归一化随机复数（备选，效率稍低） #

先生成普通随机复数，再对每个元素做归一化处理：

# 生成2×2的随机复数数组后归一化
arr = rand(ComplexF64, 2, 2)
arr ./= abs.(arr)

这种方法需要额外计算每个元素的模长再做除法，效率不如第一种方法，但代码直观易懂。

验证结果 #

可以通过计算数组元素的模长来验证：

# 检查所有元素的模是否接近1（浮点误差会有极小偏差）
println(abs.(unit_complex_arr))

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