




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
x12.5服从常见连续型分布的随机数的产生2.5.1服从均匀分布随机数的生成2.5.2服从指数分布随机数的生成2.5.3服从正态分布随机数的生成2.5.4服从卡方分布随机数的生成x2例1
产生服从区间(a,b)上的均匀分布的随机数。当Y服从(0,1)上的均匀分布,则X服从(a,b)上的均匀分布。2.5.1服从均匀分布随机数的生成np.random.uniform(1,10,10)x3x4
2.5.2服从指数分布随机数的生成np.random.exponential(0.2,20)x5x6
2.5.3服从标准正态分布随机数的生成x7
x8
算法np.random.standard_normal(1000)x9importrandomimportmathimportmatplotlib.pyplotaspltimportnumpyasnpx=[]y=[]Pi=3.1415926defNormRandom(n):for_inrange(n):r1=1.0*random.random()r2=1.0*random.random()x.append(math.sqrt(-2.0*math.log(r1))*math.cos(2.0*Pi*r2))y.append(math.sqrt(-2.0*math.log(r1))*math.sin(2.0*Pi*r2))returnx,yn=1000X,Y=NormRandom(n)plt.hist(X,bins=10,cumulative=False,normed=True)x=np.arange(-5,5,0.01)plt.plot(x,(1/np.sqrt(2*Pi))*np.exp(-(x**2)/2),color='r')x10卡方分布是n个独立的标准正态分布的平方之和,可以先生成n个独立的标准正态分布;再将其平方,最后求和得到卡方分布。2.5.4服从卡方分布随机数的生成思路1x11结论1结论2若X服从Ga(n/2,1),则Y=2X服从自由度为n的卡方分布当n为偶数时,k=n/2是整数,若X服从Ga(k,1),则Y=2X服从自由度为n=2k的卡方分布。当n为奇数时,k=(n-1)/2整数,若X服从Ga(k,1),则令Y=2X服从自由度为2k的卡方分布,Z为标准正态分布的平方,即自由度为1的卡方分布,则W=Y+Z服从自由度为n的卡方分布。且独立,则Y服从Ga(n,1);思路2x12(3)当n为偶数时,k=n/2,抽样公式为X=2Y;n为奇数时,k=(n-1)/2,取Z为标准正态分布的平方,抽样公式为X=2Y+Z。(1)当n为偶数时,k=n/2;当n为奇数时,k=(n-1)/2。产生k个随机数(2)令,则Y服从Ga(k,1);算法np.random.chisquare(12,50)x13importrandomimportmathimportmatplotlib.pyplotaspltimportnumpyasnpPi=3.1415926x=[]defKaFangRandom(Freedom,n):ifFreedom<1orn<1:returnFalse#自由度n的一半k=int(Freedom/2)foriinrange(n):u=1.0x14foriinrange(n):u=1.0forjinrange(k):m=1.0*random.random()#k个U[0,1]随机数之积u*=1.0*m#y=-lnuy=-1*math.log(u)#若n为偶数,x=2yifFreedom%2==0:x.append(2.0*y)#若n为奇数,x=2y+zx15else:r1=1.0*random.random()r2=1.0*random.random()#u服从N(0,1)u=math.sqrt(-1*2.0*math.log(r1))*math.cos(2*Pi*r2)#x=2y+z,z=u2x.append(2.0*y+u*u)returnxx16#随机数的个数n=200#自由度Freedom=13print("产生服从卡方分布的随机数为:")s=KaFangRandom(Freedom,n)print(s[0:20])plt.hist(s,bins=10,cumulative=False,normed=True)x=np.arange(0,40
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- (高清版)DB54∕T 0463-2025 羊肚菌种植技术规程
- 重症自身免疫性脑炎监测与治疗
- 少儿体育运动会活动方案
- 少先队春游活动方案
- 小班规则性游戏活动方案
- 工会活动美食品尝活动方案
- 展厅春节活动方案
- 小跳蛙线上活动方案
- 小梦室外活动策划方案
- 左邻右舍沟通日活动方案
- 安全生产奖罚管理制度
- 2025年全省民政行业职业技能大赛(孤残儿童护理员)备考试题库(含答案)
- 南京鼓楼医院合作协议书
- DB32/T 3375-2018公共场所母乳哺育设施建设指南
- Machine-Cmk-设备能力指数Cmk分析表
- 规培指导教师考试试题及答案
- 科研助理合同协议书
- 离散余弦变换教学课件
- 药学三基培训课件
- 2025-2030中国寿险行业市场现状供需分析及投资评估规划分析研究报告
- 2025-2030医疗美容产业市场深度调研及发展趋势与投资前景预测研究报告
评论
0/150
提交评论