资源描述:
《第6章 概率统计方法模型》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、第6章 概率统计方法模型本章详细介绍用概率统计方法建模的基本思路,结合实际的案例,指出如何用随机变量和概率分布来描述随机不确定事件,说明求解概率统计类模型的一般过程,并指出该类数学模型在社会调查、影响因素分析、发展趋势模拟等方面的广泛应用。§6.1 概率模型与MonteCarlo模拟6.1.1 概率模型(1)传染病随机模型假设人群总的规模为n,在总人群中,病人的数量为m,健康人的数量为s,即满足n=m+s。在人们的日常生活中,任意两人之间(包括健康人和病人)接触的概率相同,每人平均与k个人接触。当健
2、康人和病人接触时,被传染的概率为p。在以上假设的参数中,m和s通常是已知的,k和p可以通过专家的经验和统计数据获得。我们分析的目的是寻找健康人群中每天平均被感染的人数与已知参数之间的关系,以及初始参数对传染病的扩散速度和流行趋势的影响。首先以每一名健康人为研究对象,探讨其每天被感染的概率,而每一名健康人被一名指定病人接触并传染的概率等于每名健康人与指定传染者接触的概率乘以接触时感染的概率。记人群中任意两人接触的概率为q,则对每一名健康人来说,其每天接触的人数服从二项分布,分布函数为,(6.1.1)这
3、个分布的期望为k,即,进而。这样,一名健康人被一名指定病人接触并感染的概率为。进一步,对人群中的每一名健康人来说,其每天不被感染的概率为,被感染的概率为.所以,对人群中的所有健康人来说,每天被感染的人数服从二项分布,分布函数为,(6.1.2)每天被感染的人数期望为,标准差为。为了得到简明的结果,对进行近似计算,由于通常人群的总数,且根据Talyor展开,得,(6.1.3)因此,.(6.1.4)通过式(6.1.4)可以看出平均每天被感染的人数与s、m、p和k之间的关系。进而可以度量平均每天被感染人数的
4、相对误差即(6.1.5)159由式(6.1.4)可以看出,对于健康人群来说,每天平均被感染的人数与人群中每人每天平均接触的人数k,健康人与病人接触时被感染的概率p成正比。当n,p,k都确定的情况下,时,也就是在整个人群中,病人和健康人的数量各占一半时,每天被感染的人数达到最大。为了对传染病的传染过程有一个直观的了解,假设一个人口总量n=10000的人群,在日常生活中,平均每人每天接触的人数k=18,健康人与病人接触时被感染的概率p=10%,对于不同的m,平均每天被感染人数与相对误差的变化趋势如图6.
5、1.1和图6.1.2所示。可见被感染人数随着病人数量的增大而增大,直到病人数量占总人群数量的一半时达到最大,随后呈下降趋势。随着病人人口的增加每天被感染人数的相对误差一直呈减少趋势,尤为明显的是病人数量增长的前期,相对误差急剧减少。图6.1.1平均每天被感染人数的趋势图图6.1.2平均每天被感染人数的相对误差趋势R编程如下:crb<-function(m,n=10000,p=0.1,k=18)#函数{u<-(m*(n-m)*p*k)/(n-1);u}m<-1:10000plot(1:10000,cr
6、b(m),xlab="m",ylab="平均每天被的传染人数",type="l",col="blue")crb1<-function(s,n=10000,p=0.1,k=18)#相对误差函数{miugama<-((n-1-m*p*k)/((n-m)*m*p*k))^0.5;miugama}m<-1:6000plot(1:6000,crb1(m),xlab="m",ylab="相对误差",type="l",col="red")(2)企鹅繁殖模型企鹅的繁殖过程是一个典型的随机不确定模型。首先,每只母企鹅
7、下蛋的数量是随机的,服从泊松分布,其次,每个企鹅蛋是否可以成功孵化也是不确定的。针对这一问题,我们在合理假设的基础上,建立概率模型,求企鹅后代个数的期望值。据统计,企鹅生蛋的个数是服从参数为的泊松分布,即,而每个生蛋能发育成企鹅的概率为p,且每个生蛋能否发育成企鹅是彼此独立的随机事件。令代表企鹅后代的个数,且有。可见取非负的整数值0,1,2,…,对于的概率,利用全概率公式159(6.1.6)注意到每个生蛋发育成小企鹅是相互独立的,且发育成小企鹅的概率为p,因此,实际上反映了有k个生蛋,每个生蛋独立发
8、育,恰好发育成r个企鹅的概率。显然,它是一个伯努利试验,因而(6.1.7)于是有(6.1.8)得到企鹅后代的个数服从参数为的泊松分布,从而企鹅后代个数的期望值为。也说明了企鹅后代的个数与生蛋的个数以及发育成功的概率成正比。6.1.2 MonteCarlo模拟MonteCarlo(蒙特卡洛)模拟,也称为统计模拟方法。该方法是上世纪40年代。MonteCarlo模拟是由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的非常重要的数值计算方法。通常蒙特