如何在Julia中生成和为0、范围在-1到1之间的随机数列表？

基于你简化问题的转换方法 #

你已经把问题简化到生成n个[0,1]区间且和为n/2的随机数，接下来只需要一步线性变换就能得到[-1,1]区间且和为0的数：

• 假设已生成数组x，其中每个x_i ∈ [0,1]且sum(x) = n/2
• 对每个元素执行变换：y_i = 2x_i - 1
• 验证：每个y_i的范围是-1到1，总和sum(y) = 2*sum(x) - n = 2*(n/2) -n = 0，完全符合要求

生成[0,1]且和为n/2的随机数的Julia实现 #

如果还没完成这一步，这里提供两种实用方法：

方法1：归一化调整法 #

function generate_sum_half(n)
    # 生成n个[0,1]随机数
    x = rand(n)
    # 计算缩放比例，将总和调整为n/2
    scale = (n/2) / sum(x)
    return x .* scale
end

这种方法实现简单，分布近似均匀，完全满足遗传算法的需求。

方法2：Dirichlet分布变体（更接近均匀分布） #

如果追求更严格的均匀分布特性，可以借助Dirichlet分布：

using Distributions

function generate_sum_half_dirichlet(n)
    # 生成均匀单纯形上的样本
    dir_sample = rand(Dirichlet(n, 1.0))
    # 缩放至总和为n/2
    return dir_sample .* (n/2)
end

完整流程示例 #

将两步合并，直接生成符合要求的数组：

function generate_centered_random(n)
    x = rand(n)
    scale = (n/2) / sum(x)
    x_scaled = x .* scale
    # 转换到[-1,1]区间且总和为0
    return 2 .* x_scaled .- 1
end

# 测试示例
n = 5
result = generate_centered_random(n)
println("生成数组：", result)
println("范围验证：所有元素∈[-1,1]？", all(-1 .≤ result .≤ 1))
println("总和验证：sum(result)≈0？", isapprox(sum(result), 0; atol=1e-10))

内容的提问来源于stack exchange，提问作者Dustin Miao