我需要生成具有以下属性的随机数。
我尝试了几种语法，但没有任何效果，例如
r=1+8.*rand(100,1);
这给了我一个 1-9 之间的随机数，但它不是整数（例如 5.607 或 4.391），并且每次我计算平均值时它都会变化。
【问题讨论】：
您或许可以根据 Matlab 的randi 函数定义满足您要求的函数。但要小心，很容易定义不产生随机数的随机数生成器的函数。
另一种方法可能适合 - 创建一个概率分布来满足您的要求。在这种情况下，您需要一个由 9 个浮点数组成的向量，它们总和为 1，并且分别表示第 i 个整数出现的概率。例如，分布可能由以下向量描述：
[0.1 0.1 0.1 0.1 0.2 0.1 0.1 0.1 0.1]
这些将区间 [0,1] 分成 9 个部分。然后，取你最喜欢的 rng，它生成 [0,1) 范围内的浮点数并生成一个数字，假设它是 0.45。顺着 0 到 1 的区间读取，发现在第 5 个区间，所以返回整数 5。
显然，我懒得给你一个向量，它给出了 6 作为分布的平均值，但这对你来说应该不难弄清楚。
【讨论】：
这是一个带有循环的算法，通过根据当前迭代的平均值将随机数从向量的一半重新生成到另一个向量，从而达到所需的平均值 xmean（具有所需的精度 xeps）。通过我的测试，它很快就达到了平均值。
n = 100;
xmean = 6;
xmin = 1;
xmax = 9;
xeps = 0.01;
x = randi([xmin xmax],n,1);
while abs(xmean - mean(x)) >= xeps
if xmean > mean(x)
x(find(x < xmean,1)) = randi([xmean xmax]);
elseif xmean < mean(x)
x(find(x > xmean,1)) = randi([xmin xmean]);
end
end
x 是您需要的输出。
【讨论】：
randi（生成随机整数）首先出现在 MATLAB 2008b 版本中。可能你的年龄更大。看到这个SO answer。
您可以使用randi 获取随机整数
【讨论】：
您可以使用 floor 将随机数截断为仅整数值：
r = 1 + floor(9 * rand(100,1));
获得指定的平均值有点棘手；这取决于你追求什么样的分布。
【讨论】：
如果分布并不重要并且您只关心平均值，那么有一个特别简单的函数可以做到这一点：
function x=myrand
x=6;
end
【讨论】：
6。
在您设计随机数生成器之前，您需要指定它应该从中提取的分布。您只是部分地做到了这一点：即，您指定它从 [1,9] 中的整数中提取，并且它具有您希望能够指定的平均值。这仍然留下了无限的分布可供选择。您希望您的发行版具有哪些其他属性？
编辑以下评论：概率分布中任何有限样本的平均值——所谓的样本平均值——只会接近分布的平均值。没有办法解决这个问题。
如上所述，在域 [1,9] 中的整数上的最简单（在最大熵意义上）分布是指数分布：即，
p = @(n,x)(exp(-x*n)./sum(exp(-x*(1:9))));
参数x 确定分布均值。对应的累积分布为
c = cumsum(p(1:9,x));
要从分布p 中抽取，您可以从 [0,1] 中抽取一个随机数，然后找出它位于 c 的哪个子区间：即
samp = arrayfun(@(y)find(y<c,1),rand(n,m));
将返回一个从p 提取的整数数组[n,m]。
【讨论】：