2009-1-1(遗传)

第七章 概率统计模型 第一讲 遗传模型（2课时）

教学目的

掌握概率统计模型。

教学内容

运用概率统计模型方法建立遗传模型及随机存储模型并求解。

模型Ⅰ 遗传模型

1．模型背景与问题提出

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

问题：某植物园中一种植物的基因型为AA，AB和BB.现计划采用AA型植物与每种基因型植物相结合的方案培育植物后代，试预测，若干年后，这种植物的任一代的三种基因型分布情况。

2．模型假设

（1）按问题分析，后代从上一代亲体中继承基因A或B是等可能的,即有双亲体基因型的所有可能结合使其后代形成每种基因型的概率分布情况如表1。

下一代基因型（n代） AA AB BB

(2) 以an,bn和cn分别表示第n代植物中基因型为AA，AB和BB的植物总数的百分率，

上一代父-母基因型（n-1代） AA-AA 1 0 0 AA-AB 1/2 1/2 0 AA-BB 0 1 0 表1 AB-AB 1/4 1/2 1/4 AB-BB 0 1/2 1/2 BB-BB 0 0 1 x(n)表示第n代植物的基因型分布，即有

x(n)?an?????bn?, n?0,1,2,? (1) ?c??n?特别当n=0时，x(0)?(a0,b0,c0)T表示植物基因型的初始分布(培育开始时所选取各种基因型分布)，显然有a0?b0?c0?1.

3．模型建立

注意到原问题是采用AA型与每种基因型相结合，因此这里只考虑遗传分布表的前三列。 首先考虑第n代中的AA型，按上表所给数据，第n代AA型所占百分率为

an?1?an?1?1?bn?1?0?cn?1 2即第n-1代的AA与AA型结合全部进入第n代的AA型，第n-1代的AB型与AA型结合只有一半进入第n代AA型，第n-1代的BB型与AA型结合没有一个成为AA型而进入第n代AA型，故有

1an?an?1?bn?1 (2)

2同理，第n代的AB型和BB型所占有比率分别为

1bn?bn?1?cn?1 (3)

2cn?0 (4)

将(2)、(3)、(4) 式联立，并用矩阵形式表示，得到

x(n)?Mx(n?1),

其中

(n?1,2,?) (5)

?11/20???M??01/21?

?000???利用(5)进行递推，便可获得第n代基因型分布的数学模型

x(n)?Mx(n?1)?M2x(n?2)???Mnx(0) (6)

(6)式明确表示了历代基因型分布均可由初始分布x4．模型求解

这里的关键是计算M.为计算简便，将M对角化，即求出可逆阵P，使PMP??，即有

n?1(0)与矩阵M确定。

M?P?P?1

从而可计算 M?P?P (n?1,2,?)

nn?1其中?为对角阵，其对角元素为M的特征值，P为M的特征值所对应的特征向量。分别为

?1?1, ?2?，?3?0

12?1??1??1???????p1??0?,p2???1?,p3???2?

?0??0??1???????故有

?1?1??11????1?,P??0?1?2??P?1 ???2???00??1?0????即得

1??111??1??11??????1?Mn??0?1?2??0?1?2?? n2???00???1001?0???????11???11?n1?n?1?22??11???0 nn?1??22?000??????1?11?n2?an?????1??bn???0n2?c??0?n?0???1?12n?112n?10????a0??????b0?

????c0????于是 x(n)或写为

1n1n?1?a?1?()b?()c00?n22?1n1n?1? ?bn?()b0?()c022??c?0n??由上式可见，当n??时，有

an?1,bn?0,cn?0

即当繁殖代数很大时，所培育出的植物基本上呈现的是AA型，AB型的极少，BB型不存在。

5．模型分析

（1）完全类似地，可以选用AB型和BB型植物与每一个其它基因型植物相结合从而给出类似的结果.特别是将具有相同基因植物相结合，并利用前表的第1、4、6列数据使用类似模型及解法而得到以下结果：

11an?a0?b0,,bn?0,cn?c0?b0

22这就是说，如果用基因型相同的植物培育后代，在极限情形下，后代仅具有基因AA与BB，

而AB消失了。

（2）本例巧妙地利用了矩阵来表示概率分布，从而充分利用特征值与特征向量，通过对角化方法解决了矩阵n次幂的计算问题，可算得上高等代数方法应用于解决实际的一个范例。

模型Ⅱ 随机存储模型

1．模型背景与问题提出

存储问题的数学模型涉及以下的主要经济变量：

1．需求量：某种物资在单位时间内的需求量，以D表示，如年需求量、月需求量、日需求量。需求量有时是常量，而在许多情况下则是随机变量，这时它的变化规律应当是能够掌握的。对需求量进行科学地预测和估计是解决存储问题的重要依据。

2．批量：为补充存储而供应一批物资的数量称为批量，以Q表示.由外部订货供应的批量称为订货批量；由内部生产供应的批量称为生产批量。

3． 货点；为补充存储而发生订货时的存储水平，以R表示。 4．备运期：发生订货的时间与实际收到订货入库的时间的间隔。

5．存储费：保管存货的费用，包括存储所占用资金的利息、仓库和场地费用、物资的存储损耗费用、物资的税金、保险费用等，以C1表示。

6．订货费：为补充存储而订货所支付的费用，包括准备和发出订货单的费用、货物的堆放和装运的费用等，以K表示。

7．缺货损失费：发生需求时，存储不能提供而引起的费用，包括利润的损失、信誉的损失、停工待料的损失以及没有履行交货合同的罚款等，以C2表示。

存储费、订货费和缺货损失费构成了库存的总费用，即总费用=存储费+订货费+缺货损失费. 使总费用最小是建立和求解存储模型的主要目标。为实现该目标，需要确定批量和订货点，这就是所谓存储决策.批量与订货点即决策变量.因而存储模型的主要形式有：总费用=f（批量）或总费用=f（批量，订货点），即F=f（Q）或F=f（Q，R）。

为了更具体理解随机性存储模型，先来看一个具体实例。

2．报童问题

报童每日早晨从报社以每份报纸0.30元的批发价购得当日的日报，然后以每份0.45元的零售价售出。若卖不完，则每份报纸的积压损失费为0.30元；若不够卖，则缺一份报纸造成潜在损失的缺货损失费为0.15元。该报童对以往的销量作了连续一个月的统计，其记录如表2所示。

日需求量D 频率P（D） 120 0.15 130 0.2 140 0.3 150 0.25 160 0.1 表2 销量统计

那么，报童每日应订多少份报纸，才能使总损失费最小？

假定报童每日订报Q份，并设当日需求量为D，则 当Q?D时，积压损失费为F?0.30(Q?D)； 当Q?D时，缺货损失费为F?0.15(D?Q). 于是可以将报童订报的决策与相应的总费用如表3所示

F D Q P 120 130 140 150 160 120 0.15 0 3 6 9 12 130 0.2 1.5 0 3 6 9 140 0.3 3 1.5 0 3 6 150 0.25 4.5 3 1.5 0 3 160 0.1 6 4.5 3 1.5 0 平均损失总费用 2.95 2.1 2.175 3.6 6.15 表3 订报的决策与相应的总费用 从表中可看出，当报童每日订报130份时，平均损失费用最小，最小损失总费用为2.1元。

下面建立这一报童问题模型的数学解析式，用求极值的方法求解最小损失总费用。

设平均总费用为TF(Q)，则

TF(Q)??0.30(Q?D)P(D)??0.15(D?Q)P(D). （7）

D?QD?Q为求使TF(Q)最小的Q值，解下列不等式组：

??TF(Q)?TF(Q?d)?0 ???TF(Q)?TF(Q?d)?0.,130,140,150,160}. 其中 d?min{|Q?D|}?10, 且 Q?d?S?{120Q?D上式等价于

??0.30P(D)??0.15P(D)?0?D?Q?dD?Q?d ?0.30P(D)?0.15P(D)?0.???D?QD?Q?即

本文来源：https://www.bwwdw.com/article/aieg.html