matlab 真假随机,Matlab 随机数是否随机的问题

博主在MATLAB中进行随机采样模拟,发现实际不符合理论值的百分比远高于预期。问题在于,模拟结果中实际不符合的百分比比理论值高10~20倍。理论值为0.076%,但实际达到1~2%。博主怀疑可能是由于MATLAB使用伪随机数导致的误差。
摘要由CSDN通过智能技术生成

谢老师,您好,我今天研究了一个关于随机采样的问题,本来以为很简单,但是matlab做出来的效果和我预期的相差很大。请指教。

问题描述,对于一个含有n个样本的样本空间,作k * n 次随机采样,k 是多少时可以保证有99.9%以上的概率含有n个样本至少出现一次?

理论出现所有元素的概率p0, p0 = 1- (1 - 1/n) ^(k*n); 当p0 > 99.9%时解出k。

模拟程序如下:

%对于随机取样的情况下需要多大的样本量可以保证99%以上的情况取到所有的样本clear;clc;

numOfElements = 20; %样本容量

factor = 7 ; %放大系数

targetMatchPct =1-(1-1./numOfElements).^(numOfElements*factor); %理论匹配系数

spaceSize = factor*numOfElements;%测试空间大小

numOfTests = spaceSize*10; %模拟次数

fprintf('Size of sample space %d, factor %d\n',spaceSize,factor);

fprintf('Number of simulations %d\n',numOfTests);

res = zeros(1,numOfTests);

for t = 1:numOfTests

randSet = randi(numOfElements,[spaceSize,1]); %产生随机样本

res(t) = numel(unique(randSet)); %随机样本中包含的总样本

end

numOfUnMatch = sum(res

unmatchPct = numOfUnMatch / numOfTests * 100.0; %实际不符合的百分比(这个值要比理论值大10~20倍)

fprintf('Unmatched percent %3.3f%%(practise) vs %3.3f%%(theorical)\n',unmatchPct,(1-targetMatchPct)*100.0);

但是,实际模拟结果,实际不符合的百分比要比理论值大10~20倍。理论值应该是0.076%,而实际值是1~2%。

请问:

1. 我的理论和模拟对马?

2. 如果正确的话,这个问题是怎样引起的?是否是因为matlab是伪随机数引起的?

非常感谢!

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值