概率统计模型.doc

1、概率统计模型模型 遗传模型1模型背景与问题提出所谓常染色体遗传，是指后代从每个亲体的基因中各继承一个基因从而形成自己的基因型.如果所考虑的遗传特征是由两个基因 A 和 B 控制的，那么就有三种可能的基因型：AA，AB 和 BB。例如，金鱼草是由两个遗传基因决定它开花的颜色，AA 型开红花，AB型的开粉花，而 BB 型的开白花.这里的 AA 型和 AB 型表示了同一外部特征(红色) ，则人们认为基因 A 支配基因 B，也说成基因 B 对于 A 是隐性的。当一个亲体的基因型为 AB，另一个亲体的基因型为 BB，那么后代便可从 BB 型中得到基因 B，从 AB 型中得到 A 或B，且是等可能性地得到

2、。问题：某植物园中一种植物的基因型为 AA，AB 和 BB.现计划采用 AA 型植物与每种基因型植物相结合的方案培育植物后代，试预测，若干年后，这种植物的任一代的三种基因型分布情况。2模型假设（1）按问题分析，后代从上一代亲体中继承基因 A 或 B 是等可能的,即有双亲体基因型的所有可能结合使其后代形成每种基因型的概率分布情况如表 1。上一代父 -母基因型（n-1 代）下一代基因型（n 代） AA-AA AA-AB AA-BB AB-AB AB-BB BB-BBAA 1 1/2 0 1/4 0 0AB 0 1/2 1 1/2 1/2 0BB 0 0 0 1/4 1/2 1表 1(2) 以 和

3、分别表示第 n 代植物中基因型为 AA， AB 和 BB 的植物总数的百分率，nba,c表示第 n 代植物的基因型分布，即有)(nx(1),)(ncbax,210特别当 n=0 时， 表示植物基因型的初始分布(培育开始时所选取各种基Tax),(0)0因型分布)，显然有 .10cb3模型建立注意到原问题是采用 AA 型与每种基因型相结合，因此这里只考虑遗传分布表的前三列。首先考虑第 n 代中的 AA 型，按上表所给数据，第 n 代 AA 型所占百分率为11102nnncba即第 n-1 代的 AA 与 AA 型结合全部进入第 n 代的 AA 型，第 n-1 代的 AB 型与 AA 型结合只有一半

4、进入第 n 代 AA 型，第 n-1 代的 BB 型与 AA 型结合没有一个成为 AA 型而进入第 n代 AA 型，故有(2)112nnba同理，第 n 代的 AB 型和 BB 型所占有比率分别为(3)1nncb(4)0将(2)、(3)、(4) 式联立，并用矩阵形式表示，得到(5),)1()(nnMx),2(其中 01/利用(5)进行递推，便可获得第 n 代基因型分布的数学模型(6)0()2()1()( xMxMxnn(6)式明确表示了历代基因型分布均可由初始分布 与矩阵 M 确定。)0(4模型求解这里的关键是计算 .为计算简便，将 M 对角化，即求出可逆阵 P，使 ，n 1即有 1P从而可计

5、算 1PMn),2(其中 为对角阵，其对角元素为 M 的特征值，P 为 M 的特征值所对应的特征向量。分别为，,12031,0,321pp故有1021,21PP即得 1022101nnM00211nn于是 01)( 021cbacbaxnnn或写为 0)21(0nnnccba由上式可见，当 时，有0,1nncba即当繁殖代数很大时，所培育出的植物基本上呈现的是 AA 型， AB 型的极少， BB 型不存在。5模型分析（1）完全类似地，可以选用 AB 型和 BB 型植物与每一个其它基因型植物相结合从而给出类似的结果.特别是将具有相同基因植物相结合，并利用前表的第 1、4、6 列数据使用类似模型及

6、解法而得到以下结果： 002,21bcbann 这就是说，如果用基因型相同的植物培育后代，在极限情形下，后代仅具有基因 AA 与BB，而 AB 消失了。（2）本例巧妙地利用了矩阵来表示概率分布，从而充分利用特征值与特征向量，通过对角化方法解决了矩阵 n 次幂的计算问题，可算得上高等代数方法应用于解决实际的一个范例。模型 随机存储模型1模型背景与问题提出存储问题的数学模型涉及以下的主要经济变量：1需求量：某种物资在单位时间内的需求量，以 D 表示，如年需求量、月需求量、日需求量。需求量有时是常量，而在许多情况下则是随机变量，这时它的变化规律应当是能够掌握的。对需求量进行科学地预测和估计是解决存储

7、问题的重要依据。2批量：为补充存储而供应一批物资的数量称为批量，以 表示.由外部订货供应的Q批量称为订货批量；由内部生产供应的批量称为生产批量。3 货点；为补充存储而发生订货时的存储水平，以 R 表示。4备运期：发生订货的时间与实际收到订货入库的时间的间隔。5存储费：保管存货的费用，包括存储所占用资金的利息、仓库和场地费用、物资的存储损耗费用、物资的税金、保险费用等，以 表示。1C6订货费：为补充存储而订货所支付的费用，包括准备和发出订货单的费用、货物的堆放和装运的费用等，以 K 表示。7缺货损失费：发生需求时，存储不能提供而引起的费用，包括利润的损失、信誉的损失、停工待料的损失以及没有履行交

8、货合同的罚款等，以 表示。2存储费、订货费和缺货损失费构成了库存的总费用，即总费用=存储费+订货费+缺货损失费. 使总费用最小是建立和求解存储模型的主要目标。为实现该目标，需要确定批量和订货点，这就是所谓存储决策.批量与订货点即决策变量.因而存储模型的主要形式有：总费用= （批量）或总费用= （批量，订货点） ，即 F= （ ）或 F= （ ，R） 。ff fQf为了更具体理解随机性存储模型，先来看一个具体实例。2报童问题报童每日早晨从报社以每份报纸 0.30 元的批发价购得当日的日报，然后以每份 0.45 元的零售价售出。若卖不完，则每份报纸的积压损失费为 0.30 元；若不够卖，则缺一份报

9、纸造成潜在损失的缺货损失费为 0.15 元。该报童对以往的销量作了连续一个月的统计，其记录如表 2 所示。表 2 销量 统计那么， 报童每日 应订多少 份报纸，才能使总损失费最小？假定报童每日订报 份，并设当日需求量为 D，则Q日需求量 D 120 130 140 150 160频率P（D） 0.15 0.2 0.3 0.25 0.1当 时，积压损失费为 ；DQ)(30.DQF当 时，缺货损失费为 .15于是可以将报童订报的决策与相应的总费用如表 3 所示120 130 140 150 160F DQ P 0.15 0.2 0.3 0.25 0.1平均损失总费用120 0 1.5 3 4.5

10、6 2.95130 3 0 1.5 3 4.5 2.1140 6 3 0 1.5 3 2.175150 9 6 3 0 1.5 3.6160 12 9 6 3 0 6.15表 3 订报的决策与相应的总费用从表中可看出，当报童每日订报 130 份时，平均损失费用最小，最小损失总费用为 2.1元。下面建立这一报童问题模型的数学解析式，用求极值的方法求解最小损失总费用。设平均总费用为 ，则)(QTF. （7）DQDPP)(15.0)(30.为求使 最小的 值，解下列不等式组：)(.0)()dTF其中 且 ,10|minDQd .16,54,312Sd上式等价于 QDQDddP.0)(15.)(30.

11、即 QDd.015.)()15.03(故 （8）QDdQDP).(3)(亦即).()120(3.0)()120( QPdQP由于 。3,515. P因 此且可以看到，上述结果与通过列表得到的结果是一致的。报童问题是一个离散型问题.若考虑相应的连续型问题，则总费用公式为 QQxdPxxdPTF0 ).()(.0)()(3.)(这里， 为一定时期内销售量的概率密度.为求总费用的最小值，令xP.0)(dTF得 QxP0.15.)()15.3(于是 .3)(*0d问题的关键成为如何从这个积分等式中求出 ，其求法通常用迭代法。*利用求极值的数学方法求解存储模型，这是解决存储问题的主要思路。尤其对于连续型

12、存储模型，用求极值的方法求解模型就显得更为有效和更为重要。存储问题中的随机性主要由以下两个因素产生；第一，对物资的需求量经常发生随机波动；第二，订货的到达时间经常发生随机性的提前或推迟。3不允许缺货情形由于需求量是随机的，所以，可考虑其平均需求量，而且不允许缺货也只是指在一定置信度下的不允许缺货。设 D 为年平均需求，则类似于确定性存储的 模型，可得到相应的最佳批量 如EOQ*Q下：（9）.21*CKDQ这里，K 为一次定购费， 为该种物资一个单位存储一年的费用。为在一定置信度下对不缺货提供安全保证，可将安全库存量加到正常存货中以提供所希望达到的服务水平（即不缺货的概率） 。这时，有. （10

13、）lR式中，R 为订货点， 分别为备运期内的销售量 L 的均值与均方差， 为安全库存系和 数， 为安全库存量。安全库存系数 即为给定置信度 下的上 100 百分位点，其值满足等式1，可通过查概率分布表得到。)(XP因此，订货策略为，当备运期大于零时，若存储量降低到 R，则以 为订货量进行订*Q货。例 6. 设某公司订购一种备件，一次订货费为 60 元，年平均需求量为 500 件，每件年存储费为 40 元，备运期 8 天，备运期中的销售量服从均值为 15、均方差为 2 的正态分布.为使不缺货的概率达到 99.9%且总费用最小，问订货点是多少，每次订多少件？注意到 D=500 件/年，K=60 元

14、， =40 元，则1C件.3940562*Q根据不缺货的概率达到 99.9%，查正态分布表得 =3，订货点为件。215R故订货点为 21 件，每次订货 39 件。4允许缺货情形设 同前， 为单位缺货损失费，并设存储量降到 R 时订货，订货数量为 ，1,CKD2 Q备运期中的需求量 服从密度为 的分布函数 ，则在缺货要补的情况下，订货刚x)(xf)(xF到之前的平均存储量（平均最小存储量）与订货刚到之后的平均存储量（平均最大存储量）分别为，RRdxfQdxf00)()(与则年平均存储量为 .Q2年平均存储费为 .RdxfC01)(年平均订货费为 KD/Q。当备运期中的需求量超过订货点 R 时，就

15、发生缺货，因此，缺货量的均值为.dxfx)(故年平均缺货损失费为.RfQDC)(2于是年总费用 为),(RTFRRdxfxQDCKTF.)(2),(01为求 的最小值，令),(QRTF. RRdxfdxfTF021 0)()(),(可得. （11）RDCQdxf021)(由（11）得.12)(dxfRxK故解得最佳批量 与订货点 满足如下方程组：*Q*R121() (12)()()3RCFDKxfdFR 最佳批量 和订货点 可按以下步骤解出：*Q* 取 ;112CD 将 代入(12)求 R ;1 将 R= 代入 (13)求 ;12Q 将 代入(12)。2重复、一直迭代到收敛为止，最后得到的即为

16、最佳值 和 。*QR设某公司购进某种物资，其年平均需求量为 1000 件，每件的年存储费为元，一次订货费为 10 元，缺货损失费每年每件元，备运期的需求量服从0，200上的均匀分布，试求最佳批量与最佳订货点。注意到1000 件，10 元， 元， 元。由备运期的需求量服从21C20，200上的均匀分布可知，其概率密度为f(x)= , x 0, 200201将其代入（12）式和（13）式。由（12）式，得.2010QdxR即. （14）1/2QR由（13）式得 2201050 RRdx（15）.12取 件，代入（.0）式，得 件。021Q192R将 代入（15）式，得 。1R9.132Q将 代入（14）式，得 。2 7R再将 代入（15）式，得 。.043再将 代入（14）式，得 。3Q19因 与 相等，故得到订货点 R =192, 件。2R*Q