generate random numbers in a interval from specify distribution

发表于 | 更新于 | 分类于 | 评论数：

问题

如何利用一个分布（高斯分布），生成指定区间内的随机数？
比如在[-100,100]区间内，利用生成分布产生200个随机数，这个其实是叫截断正态分布。
思路：
我的想法是先使用正态分布随机生成一个数，然后判断这个数是否是满足区间内，如果不满足，就生成下一个。

代码实现

C

1
2
3
4
```

### python
``` python

matlab

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
```

### R
``` R
rand_interval = function(min, max, number)
{
  sd = (as.numeric(max)-as.numeric(min))/3
  count = 0
  res = 1:number
  while(TRUE)
  {
    x = rnorm(1, mean=0,sd=sd)
    if(x >= min && x <= max)
    {
      count = count+1
      res[count] = x
    }
  }
  return (res)
}

min = -100
max = 100
number = 20
x1 = rand_interval(min, max, number)
x2 = rand_interval(min, max, number)
x1
x2

参考文献

1.https://www.ilovematlab.cn/thread-302139-1-1.html
2.https://stats.stackexchange.com/questions/113230/generate-random-numbers-following-a-distribution-within-an-interval