概率论与数学建模.doc

^` 概率论与数学建模 基础知识部分 一、概率论： 1、概率：刻化某一事件在一次试验中发生的可能性大小的数。 注：事件指随机事件（可重复、可预测、结果明确） 例如抛骰子，抛一枚硬币。 2、常见的随机变量：X （1）离散型： 泊松分布： 实际应用：时间t内到达的次数； （小概率事件）一本书中一页中的印刷错误数； 某地区在一天内邮件遗失的信件数； 某一天内医院的急症病人数； 某一地区一个时间间隔内发生交通事故的次数； 一个时间间隔内某种放射性物质发出的经过计数器的α粒子数等等…… （2）连续型： 指数分布： 其中为常数 ，记为 特点:无记忆性。即是 一个元件已经使用了s小时，在此情形下，它总共能使用至少s+t小时的概率，与开始使用时算起它至少能使用t小时的概率相等，即元件对已使用过s小时无记忆。 实际应用：（可靠性理论、排队论）许多“等待时间”都服从指数分布；一些没有明显“衰老”迹象的机械元器件（如半导体元件）的寿命也可也用指数分布来描述…… 正态分布： 记为 标准正态分布： 正态分布标准化：若,则，标准化的目的在于能够方便查阅书后的标准正态分布表。 ““原则： ““原则被实际工作者发现,工业生产上用的控制图和一些产品质量指数都是根据3σ原则制定。 3、随机变量的特征数（数字特征）： 均值（期望）： 方差： 中心极限定理：是独立同分布的随机变量序列，且 则有： 模型一、轧钢中的浪费模型： 问题：将粗大的钢坯制成合格的钢材需要两道工序：粗轧（热轧），形成刚才的雏形；精轧（冷轧），得到规定长度的成品材料。由于受到环境、技术等因素的影响，得到钢材的长度是随机的，大体上呈正态分布，其均值可以通过调整轧机设定，而均方差是由设备的精度决定，不能随意改变。如果粗轧后的钢材长度大于规定长度，精轧时要把多余的部分切除，造成浪费；而如果粗轧后的钢材长度小于规定长度，则造成整根钢材浪费。如何调整轧机使得最终的浪费最小。 （1） 问题概述：成品材料的规定长度已知为，粗轧后的钢材长度的标准差为，粗轧后的钢材长度的均值，使得当轧钢机调整到m进行粗轧，然后通过精轧以得到成品材时总的浪费最少。 （2） 问题分析：精轧后的钢材长度记为，的均值记为m，的方差为，按照题意，。概率密度函数记为f（x），当成品钢材的规定长度给定后，记的概率为， =p（）。在轧钢过程中产生的浪费由两种情况构成：若，则浪费量为；若，则浪费量为。注意到当很大时，的可能性增加，浪费量同时增加；而当很小时，的可能性增加，浪费量也增加，因此需要确定一个合适的使得总的浪费量最小。 （3） 模型建立与求解： 这是一个优化模型，建模的关键是选择合适的目标函数，并用 ，，m把目标函数表示出来。根据前面的分析，粗轧一根钢材平均浪费长度为： 利用，，和 由（1）得：W=m- 以W为目标函数是否合适？ 由于轧钢的最终产品是成品材，浪费的多少不应以每粗轧一根钢材的平均浪费量为标准，而应该用每得到一根成品材浪费的平均长度来衡量。因此目标函数为： 因为是已知的常数，所以目标函数可以等价的取为： 其中， 易见平均每得到一根成品钢材所需要的刚才长度，问题就转化为求m使达到最小。 令则(2)式可表为： 其中： 可用微分法解的极值问题。不难推出最优值Z应满足方程： （*） 记可根据标准正态分布的函数值和制成表格式给出图形。 Z -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 227.0 56.79 18.10 7.260 3.477 1.680 Z 0 0.5 1.0 1.5 2.0 2.5 1.253 0.876 0.656 0.516 0.420 0.355 由上表可得方程（*）的根Z* 注：当给定λ>F（0）=1.253时，方程（*）不止一个根，但是可以证得，只有唯一负根Z*<0，才使取得极小值。 模型二、（美国）一个地区911应急服务中心在过去的一年内平均每月要收到171个房屋火灾电话，基于这个资料的，火灾率被估计为每月171次，下个月收到火灾报警电话只有153个，这表明房屋火灾率实际上实际上是减少了，或是或是它只是一个随机波动？ 分析：Xn——第n-1次和第n次火灾之间的时间（月），X1…，Xn，…是独立的且每一个Xn服从参数为λ的指数分布，λ为报告的房屋火灾率（月），即是：，（Xi>0） 目标：给定λ=171，确定每月收到153次这样的少的电话报警的概率有多大？或者说每月收到153是否属于正常范围内？ 建模：，， 将代入得： （利用3原理）： 若要有95%的把握，则： 若要有99%的把握，则： 选择95%的把握得到： 将=171，n=153代入（1），有： 即： 因此我们的观察值是在正常的变化范围之内 结论：断言火灾报警率降低的证据不充分，它可能是正态随机变量的正常结果。当然，若每月都连续这样低，则需重新评估。 灵敏度分析：当=171代入（1）得： 因为对任何的，区间总会包含1，即在之间都属于正常范围。 对于“每月171次”的假设的敏感性分析。去掉特殊性，假设每月的均值是，我们有一个月的报警电话次数的观测值n=153，代入（1），有： 因为对于任何的之间总会包含1，所以=153属于正常的变化范围。 随机过程与数学建模 基础知识部分 随机过程： 热噪声电压：电子元件或器件由于内部微观粒子的随机热骚动所引起的端电压称为热噪声电压。它在任一时刻t的值是一随机变量，记为V(t),不同时刻对应不同的随机变量，当时间在某区间，譬如在 上推移时，热噪声电压表现为一族随机变量，记为：{V(t),t>=0}。由于热骚动的随机性，在相同条件下每次测量都将产生不同的电压——时间函数。这样，不断的独立的测量就可以得到一族不同的电压——时间函数。 t V1(t) V2（t） V3（t） t t tj tj tj 设T是一无限实数集，我们把依赖于参数的一族（无限多个）随机变量称为随机过程。记为{X(t), }。这时每一个，X(t)是一随机变量，T叫做参数集。 把t看作为时间，称X(t)为时刻t时过程的状态，而X(t)=x或是t=t1时过程处于状态x。对于一切的，X(t)的所有可能的一切值的全体称为随机过程的状态空间。 马尔可夫链及其基本方程： 将时间离散化为n=0，1，2，…对每个n，系统的状态用随机变量Xn表示，设Xn可以取k个离散的值Xn=1，2，…k，且记即状态概率从Xn=i~Xn+1=j的概率记为 即转移概率。如果的取值只取决于Xn的取值及转移概率，而与X1，X2，…，Xn-1的取值无关，则称这种离散状态按照离散的时间的随机转移过程叫做马尔可夫过程。或者说此过程具有马尔可性或无后效性。 注：还可以这样表示 由状态转移的无后效性和全概率公式可以写出马尔可夫链的基本方程为 并且和应满足： （2） 引入状态概率向量和转移概率矩阵 则（1）式可表为： 由此可得 ： （2）式表明转移矩阵P是非负矩阵，且P的行和为1，称为随机矩阵。 说明：对于马尔可夫链模型最基本的问题是构造状态Xn及写出转移矩阵P，一旦有了P，那么给定初始状态概率a（0）就可以用（3）和（4）或计算任意时段n的状态概率a（n） 模型一：人寿保险公司对受保人的健康状况特别关注，他们欢迎年轻力壮的人投保，患病者和高龄人则需付较高的保险金，甚至被拒之门外。人的健康状态随着时间的推移会发生转变，且转变是随机的，保险公司要通过大量数据对状态转变的概率做出估计，才可能制定出不同年龄、不同健康状况的人的保险金和理陪金数额，下面分两种情况进行讨论： （1）健康与疾病： 其中（i，j=1，2） 0.2 0.7 0.3 0.8 n 0 1 2 3 …… 1 0.8 0.78 0.778 7/9 0 0.2 0.22 0.222 2/9 若开始处于疾病状态，即， n 0 1 2 3 …… 0 0.7 0.77 0.777 7/9 1 0.3 0.23 0.223 2/9 更一般的，当时，的趋向与上面两表相同。 结论：当时，趋向于稳定值，与初始状态无关。 （2）健康、疾病、死亡 0.65 0，25 0.18 0.1 0.8 0.02 1 n 0 1 2 …… 50 …… 1 0.8 0.757 0.1293 0 0 0.18 0.189 0.0326 0 0 0.2 0.054 0.8381 1 对于例题中的（1）小问，看出从任意状态出发经过有限次的转移都能达到另外的任意状态，而（2）小问中则不能。 正则链定义：一个有k可状态的马尔可夫链，如果存在正整数N，使从任意的状态i经N次转移都以大于0的概率到达状态j（i，j=1，2，3，…，k）则称为正则链。 Th1.若马尔可夫链的转移矩阵为P，则它是正则链的充要条件是存在正整数N，使>0（指的每一元素大于零）。 Th2.正则链存在唯一的极限状态概率π=（π1，π2，…，πk），使得当时，，π与初始状态概率a(0)无关（π称为稳定概率或稳定状态分布），满足π=π 例如： 由上面方程组可求得， 注：这能够满足我们这样的一种想法：由于随机波动，我们不能期望当系统稳定时状态变量将停留在一个数值上。我们能达到的最好的希望是状态变量的概率分布将趋于一个极限分布。 吸收链：马尔可夫链存在一种状态，系统一旦进入该状态不再会转移到其他状态，并且系统从其他任何状态出发最终都会转移到该状态。 吸收链的定义：转移概率Pij=1的状态i称为吸收状态，如果马尔可夫链至少含一个状态，并且从每一个非吸收状态出发，能以正的概率经过有限次转移到达某个吸收状态，则称这个马尔可夫链为吸收链。 模型二： 一个宠物商店出售容量为15L的水族箱，每个周末商店老板要盘点存货，开出订单。商店的订货策略是如果存货为0，就在这个周末进3个新的15L水族箱。如果，只要商店还保存一个存货，就不在进新的。这个策略是基于商店平均每周销售一个水族箱的事实提出的。这个策略是不是能够保证防止商店缺货时顾客需要水族箱而无货销售的损失？ 分析：商店在每个销售周的开始存货在1个到3个之间，一周销售的个数依赖于供给和需求两方面，需求是每周平均一个，但是是随机波动的。完全在某些周需求大于供给，即使在一周的开始就有3个水族箱的最大库存。我们希望计算需要超过供给的概率。 要解决此问题需给出关于需求的概率特征的假设，假设潜在的购买者在每周以一定的概率随机到达。 因此，在一周内潜在的购买者的数目均值为1的泊松分布。 建模与求解： 符号： 假设： 目标：计算 注：与模型的动态有关，将被用来构成转移状态矩阵P。 已知状态空间设=3 先确定P： 所以，若，则 其余的状态转移概率可以类似的计算，的得： 0.368 0.184 0.632 0.448 0.368 0.264 0.368 0.368 所以 可以看出的概率依赖于n，更具体的说，它依赖于，若=3，则 等等，为了得到关于需求多么经常超过供给的更好的想法，我们需要更多关于的信息，现在我们想找到一个唯一的渐进稳定的概率向量π，故我们有 得： 即当时，近似有 所以，根据全概率公式，有： 即在长时间的运行中，需求将有10%的时间超过供给。