- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[实验六随机数问题1
实验六 随机数1.随机数的之间的随机数,通常用表示,即
,
为的概率分布密度,记为。
2.产生区间(0,1)内均匀分布的随机数位数的数字自乘以后,得到一个位数的数字,然后取其中间的是位数,便得到一个新的位数,将新的位数的数字再一次自乘以后,又得到一个位数的数字,然后再取其中间的是位数,便得到第三个新的位数,如此不断进行下去,便得到了一个位字长的随机数序列,而将这个数列中的每个数都除以最大的位整数,便得到了一个分布在之间的伪随机数列。
3.产生区间(0,1)内均匀分布的随机数的公式是:,
。
其中,是乘因子,是模数或。第一式称做以为模数的同余式,即以除后得到的余数记为。,均为之间的整数。从上述构造过程可知(因而对也同样)最多有个相异值, 即这表明()有周期,且。因此,不是真正的随机数序列,但若充分大,在一个周期内能经得起统计检验。周期的长度取决于参数的选择时,。
时,
混合同余法:
递推公式为:,
。
其中,是非负整数。通过适当选取参数可以改善随机数的统计性质。一般取作小于的任意奇数正整数,最好使其与模互素。其他参数的选择与乘同余法相类似, 对随机数的统计检验。 (1)参数检验:子样平均值和理论平均值的差异性检验。 在(0,1)上均匀分布随机变量和的平均值和方差分别为: 如果随机数是X的N个独立观测值,令: 则它们的平均值和方差为: 由中心极限定理知统计量: 渐近服从正态N(0,1)。 给定显著性水平,根据正态分布确定临界值,据此判断与, 与的差异是否显著,从而决定能否把 看成(0,1)上均匀分布的随机变量X的N个独立样本值。 (2)分布均匀性检验: 均匀性检验又称频率检验,它用来检验用随机数(样本值)确定的经验频率和均匀分布频率是否有显著性差异。常用的统计检验方法有检验和累积频率检验(K-S检验)。 (3)独立性检验: 只要是检验随机数中前后各数的统计相关性是否显著。
#include iostream
#include ctime
#include cstdlib
using namespace std;
int main()
{
double random(double,double);
srand(unsigned(time(0)));
for(int icnt = 0; icnt != 10; ++icnt)
cout No. icnt+1 : int(random(0,10)) endl;
return 0;
}
double random(double start, double end)
{
return start+(end-start)*rand()/(RAND_MAX + 1.0);
}
/* 运行结果
* No.1: 3
* No.2: 9
* No.3: 0
* No.4: 9
* No.5: 5
* No.6: 6
* No.7: 9
* No.8: 2
* No.9: 9
* No.10: 6
*/
伪随机数生成程序:
//rand01.c
#include
static unsigned int RAND_SEED;
unsigned int random(void)
伪随机数
{
伪随机数RAND_SEED=(RAND_SEED*123+59)%65536;
return(RAND_SEED);
}
void random_start(void)
{
int temp[2];
movedata(0x0040,0x006c,FP_SEG(temp),FP_OFF(temp),4);
RAND_SEED=temp[0];
}
main()
{
unsigned int i,n;
random_start();
for(i=0;i10;i++)
printf(%u\t,random());
printf(\n);
}
文档评论(0)