如何生成满足差值和为3000的truncated normal分布随机数序列

原代码问题分析 #

• 循环逻辑错误：每次循环仅固定取预生成随机数数组的前两个元素计算差值，c的数值永远固定，不会随着循环累加，且单个差值最大不超过30，永远无法达到3000的阈值，会陷入死循环
• 未实现目标序列生成：没有将选中的随机数存入存储目标序列的d变量
• 差值计算逻辑缺失：没有按需求累加多组随机数的差值

正确实现思路 #

按常规需求「生成截断正态随机数序列，累加序列相邻元素的差值，直至累加总和≥3000」实现，若你需要其他差值计算规则，调整对应求和逻辑即可。

可运行代码 #

# 未安装依赖包先运行：install.packages("truncnorm")
library(truncnorm)
# 预生成截断正态分布随机数，调大生成量避免不够用
truncnorm_rn <- rtruncnorm(10000, a = 0, b = 30, mean = 0, sd = 1)

sum_diff <- 0 # 差值总和
target_seq <- c() # 存储符合条件的序列
# 先放入第一个随机数作为序列起始
target_seq[1] <- truncnorm_rn[1]
current_index <- 2 # 下一个待取用随机数的下标

while(TRUE) {
  # 取当前下标对应随机数加入序列
  next_num <- truncnorm_rn[current_index]
  target_seq <- c(target_seq, next_num)
  # 累加当前相邻元素的差值，不需要绝对值可去掉abs()
  sum_diff <- sum_diff + abs(target_seq[length(target_seq)] - target_seq[length(target_seq)-1])
  # 达到阈值直接跳出循环
  if (sum_diff >= 3000) {
    break
  }
  current_index <- current_index + 1
  # 防护逻辑：预生成随机数用完时自动补充，避免下标越界
  if (current_index > length(truncnorm_rn)) {
    truncnorm_rn <- c(truncnorm_rn, rtruncnorm(10000, a = 0, b = 30, mean = 0, sd = 1))
  }
}

# 结果验证
print(paste0("差值总和：", sum_diff))
print(paste0("目标序列长度：", length(target_seq)))

内容的提问来源于stack exchange，提问作者Brad Yun