第十六章差分方程模型
第十六章 差分方程模型
离散状态转移模型涉及的范围很广,可以用到各种不同的数学工具。下面我们对差分方程作一简单的介绍,下一章我们将介绍马氏链模型。
§1 差分方程
1.1 差分方程简介
规定t只取非负整数。记yt为变量y在t点的取值,则称ytyt1yt为yt的一阶向前差分,简称差分,称2yt(yt)yt1ytyt22yt1yt为yt的二阶差分。类似地,可以定义yt的n阶差分nyt。
由t、yt及yt的差分给出的方程称为yt的差分方程,其中含yt的最高阶差分的阶数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如,二阶差分方程
2ytytyt0也可改写成yt2yt1yt0。
满足一差分方程的序列yt称为差分方程的解。类似于微分方程情况,若解中含有
的独立常数的个数等于差分方程的阶数时,称此解为该差分方程的通解。若解中不含任意常数,则称此解为满足某些初值条件的特解。
称如下形式的差分方程
a0ynta1ynt1anytb(t) (1) 为n阶常系数线性差分方程,其中a0,a1,,an是常数,a00。其对应的齐次方程为
a0ynta1ynt1anyt0 (2)
容易证明,若序列yt(1)与yt(2)均为(2)的解,则ytc1yt(1)c2yt(2)也是方程(2)的解,其中c1,c2为任意常数。若yt(1)是方程(2)的解,yt(2)是方程(1)的解,则
ytyt(1)yt(2)也是方程(1)的解。
方程(1)可用如下的代数方法求其通解: (I)先求解对应的特征方程
a0na1n1a00 (3) (II)根据特征根的不同情况,求齐次方程(2)的通解。
(i)若特征方程(3)有n个互不相同的实根1,,n,则齐次方程(2)的通解为
tc11cntn (c1,,cn为任意常数)
(ii)若是特征方程(3)的k重根,通解中对应于的项为(1ktk1)t,
i(i1,,k)为任意常数。
(iii)若特征方程(3)有单重复根i,通解中对应它们的项为
1tcost2tsint,其中22为的模,arctg
为的幅角。
(iv)若i是特征方程(3)的k重复根,则通解对应于它们的项为
(1ktk1)tcost(k12ktk1)tsint
-182-
i(i1,,2k)为任意常数。
(III)求非齐次方程(1)的一个特解t。若yt为方程(2)的通解,则非齐次方程(1)的通解为tyt。
求非齐次方程(1)的特解一般要用到常数变易法,计算较繁。对特殊形式的b(t)也可使用待定系数法。例如,当b(t)btpk(t),pk(t)为t的k次多项式时可以证明:若b不是特征根,则非齐次方程(1)有形如btqk(t)的特解,qk(t)也是t的k次多项式;若b是r重特征根,则方程(1)有形如bttrqk(t)的特解。进而可利用待定系数法求出qk(t),从而得到方程(1)的一个特解t。
例1 求解两阶差分方程yt2ytt。
解 对应齐次方程的特征方程为10,其特征根为1,2i,对应齐次方程的通解为
2
ytc1cos
2
tc2sin
2
t
11
,b,故原方程的通解22
原方程有形如atb的特解。代入原方程求得a为
11
tt
2222
例2 在信道上传输仅用三个字母a,b,c且长度为n的词,规定有两个a连续出现
c1cos
tc2sin
的词不能传输,试确定这个信道容许传输的词的个数。
解 令h(n)表示容许传输且长度为n的词的个数,n1,2,,通过简单计算可求得:h(1)3,h(2)8。当n3时,若词的第一个字母是b或c,则词可按h(n1)种方式完成;若词的第一个字母是a,则第二个字母是b或c,该词剩下的部分可按h(n2)种方式完成。于是,得差分方程
h(n)2h(n1)2h(n2),(n3,4,)
其特征方程为
2220
特征根
则通解为
113,213
h(n)c1(13)nc2(13)n,(n3,4,)
利用条件h(1)3,h(2)8,求得
h(n)
2323
(13)n
223
(13)n,(n1,2,)
在应用差分方程研究问题时,我们常常需要讨论解的稳定性。对常系数非齐次线性差分方程(1),若不论其对应齐次方程的通解中任意常数c1,,cn如何取值,在t时总有yt0,则称方程(1)的解是稳定的。根据通解的结构不难看出,非齐次方
-183-
程(1)稳定的充要条件为其所有特征根的模均小于1。
1.2 常系数线性差分方程的Z变换解法
常系数线性差分方程采用解析解法比较容易,而且对其解的意义也容易理解,但采用这种解法求解常系数线性非齐次差分方程比较繁琐,通常是采用Z变换,将差分方程变换为代数方程去求解。
设有离散序列x(k),(k0,1,2,),则x(k)的Z变换定义为
X(z)Z[x(k)]
x(k)z
k0
k
(4)
其中z是复变量。显然上式右端的级数收敛域是某个圆的外部。
X(z)的Z反变换记作
x(k)Z1[X(z)]
1.2.1 几个常用离散函数的Z变换 (i)单位冲激函数(k)的Z变换 Z[(k)]
(k)z
k0
k
[1zk]k01
即单位冲激函数的Z变换为1。
(ii)单位阶跃函数U(k)的Z变换
Z[U(k)]即
Z[U(k)]
U(k)z
k0
k
1zk,
k0
z
(|z|1)
z1
(iii)单边指数函数f(k)ak的Z的变换(a为不等于1的正常数)
zkkk
Z[a]az(|z|a)
zak0
1.2.2 Z变换的性质
(i)线性性质
设Z[f1(k)]F1(z),Z[f2(k)]F2(z),则
Z[af1(k)bf2(k)]aF1(z)bF2(z)
其中a,b为常数。收敛域为F1(z)和F2(z)的公共区域。
(ii)平移性
设Z[f(k)]F(z),则
Z[f(k1)]z[F(z)f(0)], Z[f(kN)]z[F(z)
1N
N1k0
f(k)z
N1k1
k
],
Z[f(k1)]z[F(z)f(1)z], Z[f(kN)]z例3 求齐次差分方程
-184-
N
[F(z)f(k)zk]
x(k2)3x(k1)2x(k)0,x(0)0,x(1)1
的解。
解 令Z[x(k)]X(z),对差分方程取Z变换,得 z2X(z)z3zX(z)2X(z)0, X(z)
zzz
,
z23z2z1z2
对上式取z反变换,便得差分方程的解为
x(k)(1)k(2)k。
§2 蛛网模型
2.1 问题提出
在自由竞争的社会中,很多领域会出现循环波动的现象。在经济领域中,可以从自由集市上某种商品的价格变化看到如下现象:在某一时期,商品的上市量大于需求,引起价格下跌,生产者觉得该商品无利可图,转而经营其它商品;一段时间之后,随着产量的下降,带来的供不应求又会导致价格上升,又有很多生产商会进行该商品的生产;随之而来的,又会出现商品过剩,价格下降。在没有外界干扰的情况下,这种现象将会反复出现。
如何从数学的角度来描述上述现象呢? 2.2 模型假设
(i)设k时段商品数量为xk,其价格为yk。这里,把时间离散化为时段,一个时期相当于商品的一个生产周期。
(ii)同一时段的商品的价格取决于该时段商品的数量,把
ykf(xk) (5) 称之为需求函数。出于对自由经济的理解,商品的数量越多,其价格就越低,故可以假设:需求函数为一个单调下降函数。
(iii)下一时段商品数量由上一个时段的商品的价格决定,把
xk1g(yk) (6) 称之为供应函数。由于价格越高可以导致产量越大,故可假设供应函数是一个单调上升的函数。
2.3 模型求解
在同一个坐标系中做出需求函数与供应函数的图形,设两条曲线相交于
则P因为此时x0g(y0),y0f(x0),若某个k,有xkx0,P0(x0,y0),0为平衡点。则可推出
yly0,xlx0,(lk,k1,)
即商品的数量保持在x0,价格保持在y0,不妨设x1x0,下面考虑xk,yk在图上的变化(k1,2,)。如下图所示,当x1给定后,价格y1由f上的P1
-185-
点决定,下一时段的数量x2由g上的P2点决定,y2又可由f上的P3点决定。依此类推,可得一系列的点P1(x1,y1),P2(x2,y1),P3(x2,y2),P4(x3,y2),图上的箭头表示求出Pk的次序,由图知:
k
limPk(x,y)P0(x0,y0),
即市场经济将趋于稳定。
并不是所有的需求函数和供应函数都趋于稳定,若给定的f与g的图形如下图所示,得出的P1,P2,就不趋于P0,此时,市场经济趋向不稳定。
上两图中的折线P1P2,P2P3,P3P4,形似蛛网,故把这种模型称为蛛网模型。在进行市场经济分析中,f取决于消费者对某种商品的需要程度及其消费水平,g取决于生产者的生产、管理等能力。
当已经知道需求函数和供应函数之后,可以根据f和g的性质判断平衡点P0的稳定性。利用结论:当|x1x0|较小时,P0点的稳定性取决于f与g在P0点的斜率,即当
|f'(x0)||g'(y0)| (7) 时,P0点稳定,当
|f'(x0)||g'(y0)| (8)
时,P0点不稳定。
这一结论的直观解释是:需求曲线越平,供应曲线越陡,越有利于经济稳定。 设|f'(x0)|,式得
-186-
1
|g'(y0)|,在P0点附近取f与g的线性近似,由(5),(6)
yky0(xkx0) (9) xk1x0(yky0) (10) 上两式中消去yk,得
xk1xk(1)x0 (11) (11)式对k1,2,均成立,有
xk1xk(1)x0
()xk()2xk1()(1)x0 ()2xk1()3xk2()2(1)x0 ……………………………………………… ()k2x3()k1x2()k2(1)x0 ()k1x2()kx1()k1(1)x0 以上k个式子相加,有
xk()kx1(1)x0[1()()k1] ()x1[1()]x0
k
k
(12)
此为(11)式的解。 若P0是稳定点,则应有:
limxk1x0
k
结合(12)式考虑,P0点稳定的条件是
1 (13)
即
同理,P0点不稳定的条件是
1
1 (14)
即
k
1
此时,limxk1。这与(7),(8)式是一致的。
2.4 模型的修正 在上面模型假设的第(iii)点中引进了供应函数,并且知道g取决于管理者的生产、管理水平。如果生产者的管理水平更高一些,他们在决定该商品生产数量xk1时,不仅考虑了前一时期的价格yk,而且也考虑了价格yk1。为了简化起见,不妨设xk1由
1
(ykyk1)决定,则供应函数可写成 2
-187-
xk1g(ykyk1)
2在P0附近取线性近似,则有 xk1x0
1
2
(ykyk12y0) (15)
由(9)式有
yky0(xkx0) yk1y0(xk1x0) 将上两式代入(15)式,整理得
2xk1xkxk1(1)x0,(k2,3,)
这是一个二阶线性差分方程,其特征方程为
220
经计算,可得其特征根
()28
(16) 1,2
4
结论:若方程的特征根均在单位圆内,即|1|1,|2|1,则P0为稳定点。 当8时,(16)式有两个实根,因 ()28 2,
44
则有|2|2,故此时P0不是稳定点。
当8时,(16)式有两个共轭复根,此时
|1,2
12|8()
244
22
1
2
要使P0为稳定点,只需
2
与(13)式相比,与的范围扩大了。这是由于经营者经营管理水平的提高带来的结果。 §3 商品销售量预测
在利用差分方程建模研究实际问题时,常常需要根据统计数据并用最小二乘法来拟合出差分方程的系数。其系统稳定性讨论要用到代数方程的求根。对问题的进一步研究又常需考虑到随机因素的影响,从而用到相应的概率统计知识。
例4 某商品前5年的销售量见表。现希望根据前5年的统计数据预测第6年起该
从表中可以看出,该商品在前5年相同季节里的销售量呈增长趋势,而在同一年中销售量先增后减,第一季度的销售量最小而第三季度的销售量最大。预测该商品以后的销售情况,根据本例中数据的特征,可以用回归分析方法按季度建立四个经验公式,分别用来预测以后各年同一季度的销售量。例如,如认为第一季度的销售量大体按线性增长,可设销售量yt(1)atb,由
x=[[1:5]',ones(5,1)];y=[11 12 13 15 16]';z=x\y 求得az(1)1.3,bz(2)9.5。
(1)
根据yt(1)1.3t9.5,预测第六年起第一季度的销售量为y617.3,
(t)y718.6,…。由于数据少,用回归分析效果不一定好。
如认为销售量并非逐年等量增长而是按前一年或前几年同期销售量的一定比例增长的,则可建立相应的差分方程模型。仍以第一季度为例,为简单起见不再引入上标,以yt表示第t年第一季度的销售量,建立形式如下的差分公式:
yta1yt1a2 或
yta1yt1a2yt2a3
等等。
上述差分方程中的系数不一定能使所有统计数据吻合,较为合理的办法是用最小二乘法求一组总体吻合较好的数据。以建立二阶差分方程yta1yt1a2yt2a3为例,选取a1,a2,a3使
5
[y
t3
t
(a1yt1a2yt2a3)]2
最小。编写Matlab程序如下: y0=[11 12 13 15 16]';
y=y0(3:5);x=[y0(2:4),y0(1:3),ones(3,1)]; z=x\y
求得a1z(1)1,a2z(2)3,a3z(3)8。即所求二阶差分方程为
ytyt13yt28。
虽然这一差分方程恰好使所有统计数据吻合,但这只是一个巧合。根据这一方程,可迭代求出以后各年第一季度销售量的预测值y621,y719,…等。
上述为预测各年第一季度销售量而建立的二阶差分方程,虽然其系数与前5年第一季度的统计数据完全吻合,但用于预测时预测值与事实不符。凭直觉,第六年估计值明显偏高,第七年销售量预测值甚至小于第六年。稍作分析,不难看出,如分别对每一季度建立一差分方程,则根据统计数据拟合出的系数可能会相差甚大,但对同一种商品,这种差异应当是微小的,故应根据统计数据建立一个共用于各个季度的差分方程。为此,将季度编号为t1,2,,20,令yta1yt4a2或yta1yt4a2yt8a3等,利用全体数据来拟合,求拟合得最好的系数。以二阶差分方程为例,为求a1,a2,a3使得
-189-
Q(a1,a2,a3)
[y
t9
20
t
(a1yt4a2yt8a3)]2
最小,编写Matlab程序如下:
y0=[11 16 25 12 12 18 26 14 13 20 27 15 15 24 30 15 16 25 32 17]'; y=y0(9:20);
x=[y0(5:16),y0(1:12),ones(12,1)]; z=x\y
求得a1z(1)0.8737,a2z(2)0.1941,a3z(3)0.6957,故求得二阶差分方程
yt0.8737yt40.1941yt80.6957,(t21)
根据此式迭代,可求得第六年和第七年第一季度销售量的预测值为 y2117.586,9y2519.1676
还是较为可信的。
§4 遗传模型
随着人类的进化,人们为了揭示生命的奥妙,越来越重视遗传学的研究,特别是遗传特征的逐代传播,引起人们更多的注意。无论是人,还是动植物都会将本身的特征遗传给下一代,这主要是因为后代继承了双亲的基因,形成自己的基因对,基因对将确定后代所表现的特征。下面,我们来研究两种类型的遗传:常染色体遗传和x链遗传。根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因型的分布。
4.1 常染色体遗传模型 常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对,基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因A和a控制的,那么就有三种基因对,记为AA,Aa,aa。例如,金鱼草由两个遗传基因决定花的颜色,基因型是AA的金鱼草开红花,Aa型的开粉红色花,而aa型的开白花。又如人类眼睛的颜色也是通过常染色体遗传控制的。基因型是AA或Aa的人,眼睛为棕色,基因型是aa的人,眼睛为蓝色。这里因为AA和Aa都表示了同一外部特征,我们认为基因A支配基因a,也可以认为基因a对于A来说是隐性的。当一个亲体的基因型为Aa,而另一个亲体的基因型是aa时,那么后代可以从aa型中得到基因a,从Aa型中或得到基因A,或得到基因a。这样,后代基因型为Aa或aa的可能性相等。下面给出双亲
例5 农场的植物园中某种植物的基因型为AA,Aa和aa。农场计划采用AA型的植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任一代的三种基因型分布如何?
(a)假设
令n0,1,2,。
-190-
(i)设an,bn和cn分别表示第n代植物中,基因型为AA,Aa和aa的植物占植物总数的百分率。令x当n0时
x(0)a0
(n)
为第n代植物的基因型分布:
x(n)an
bnb0
Tcn Tc0
表示植物基因的初始分布(即培育开始时的分布),显然有
a0b0c01
(ii)第n代的分布与第n1代的分布之间的关系是通过上面的表格确定的。 (b)建模 根据假设(ii),先考虑第n代中的AA型。由于第n1代的AA型与AA型结合,后代全部是AA型;第n1代的Aa型与AA型结合,后代是AA型的可能性为第n1代的aa型与AA型结合,后代不可能是AA型。因此当n1,2,时
an1an1即
1
;而2
1
bn10cn1 2
1
anan1bn1 (17)
2
类似可推出
bn
1
bn1cn1 (18) 2
cn0 (19)
将(17),(18),(19)式相加,得
anbncnan1bn1cn1
根据假设(i),有
anbncna0b0c01
对于(17),(18),(19)式,我们采用矩阵形式简记为 x其中
(n)
Mx(n1),n1,2, (20)
1M0
0
由(20)式递推,得
1
2120
01 0
x(n)Mx(n)M2x(n2)Mnx(0) (21)
(21)式给出第n代基因型的分布与初始分布的关系。
编写如下Matlab程序: syms n a0 b0 c0
-191-
M=sym('[1,1/2,0;0,1/2,1;0,0,0]'); [p,lamda]=eig(M);
x=p*lamda.^n*p^(-1)*[a0;b0;c0]; x=simple(x) 求得
nn1
11an1b0c0
22
nn111
(22) bnb0c0
22
cn0
1
当n时,0,所以从(22)式得到
2
an1,bn0,cn0
即在极限的情况下,培育的植物都是AA型。
(c)模型的讨论
若在上述问题中,不选用基因AA型的植物与每一植物结合,而是将具有相同基因
(n)n(0)
并且xMx,其中
n
1M0
0
1
41214
00 1
编写如下Matlab程序:
syms n a0 b0 c0
M=sym('[1,1/4,0;0,1/2,0;0,1/4,1]'); [p,lamda]=eig(M);
x=p*lamda.^n*p^(-1)*[a0;b0;c0]; x=simple(x)
求得
-192-
11n1ana0b0
22
n1
(23) bnb0
211n1cnc0b02211
当n时,ana0b0,bn0,cnc0b0。因此,如果用基因型相同
22
的植物培育后代,在极限情况下,后代仅具有基因AA和aa。
4.2 常染色体隐性病模型
现在世界上已经发现的遗传病有将近4000种。在一般情况下,遗传病与特殊的种族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征,不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者的概率。
(a)假设
(i)常染色体遗传的正常基因记为A,不正常基因记为a,并以AA,Aa,aa分别表示正常人,隐性患者,显性患者的基因型。
(ii)设an,bn分别表示第n代中基因型为AA,Aa的人占总人数的百分比,记
x
(n)
an
,n1,2, bn
(iii)为使每个儿童至少有一个正常的父亲或母亲,因此隐性患者必须与正常人结
(b)建模 由假设(iii),从第n1代到第n代基因型分布的变化取决于方程
1
bn1 21
bn0an1bn1
2
(n)(n1)
所以xMx,n1,2,,其中
anan1
-193-
112 M
10
2
(n)n(0)(0)
如果初始分布x已知,那么第n代基因型分布为xMx,n1,2,。
易知
n
1an1b02
n1,2, (24)
n1bb0n2
当n时,an1,bn0,隐性患者逐渐消失。从(24)式中可知
1
bnbn1
2
1
这说明每代隐性患者的概率是前一代隐性概率患者的。
2
(c)模型讨论
研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
bn
bn1
(25) 1bn1
2
下面给出数值的例子:
某地区有10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式可知下一代(大约27年)的隐性患者将减少到5%;如果随机结合,根据(25)式,可以预言下一代人中有9.5%是隐性患者,并且可计算出大约每出生400个黑人孩子,其中有一个是显性患者。
4.3 X链遗传模型
雌性具有两个基因AA,或Aa,或aa。X链遗传是指雄性具有一个基因A或a,
其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与X链遗传有关的近亲繁殖过程。
(a)假设
(i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它们产生的后代中任选两个结成配偶。如此继续下去。
(ii)父体与母体的基因型组成同胞对,同胞对的形式有(A,AA),(A,Aa),
(A,aa),(a,AA),(a,Aa),(a,aa)六种。初始一对雌雄的同胞对,是这六种类型
-194-
(iii)在每一代中,配偶的同胞对也是六种类型之一,并有确定的概率。为计算这些概率,设an,bn,cn,dn,en,fn分别是第n代中配偶的同胞对为(A,AA),(A,Aa),
(A,aa),(a,AA),(a,Aa),(a,aa)型的概率,n0,1,。令
fn,n0,1,
(iv)如果第n1代配偶的同胞对是(A,Aa)型,那么它们的雄性后代将等可能地得到基因A和a,它们的雌性后代的基因型将等可能地是AA或Aa。又由于第n代雌雄结合是随机的,那么第n代配偶的同胞对将等可能地为四种类型(A,AA),(A,Aa),(a,AA),(a,Aa)之一。对于其它类型的同胞对,我们可以进行同样分析,因此有
(n)(n1) xMx,n1,2, (26)
bn
cn
dn
en
T
x(n)an
其中
100M
000
从(26)式中易得
1
4140
0001
1
000411
1044
1
000
4
14100
4
000 001
x(n)Mnx(0),n1,2,
编写如下Matlab程序:
syms n a0 b0 c0 d0 e0 f0
M=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0; 0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1]; M=sym(M);
[p,lamda]=eig(M);
x=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0]; x=simple(x)
由上述程序计算结果可以看出
-195-
2121abcd03030303e00
0(n)
当n时,x,
0
02121bcdef000003333
因此,在极限情况下所有同胞对或者是(A,AA)型,或者是(a,aa)型。如果初始的父母体同胞对是(A,Aa)型,即b01,而a0c0d0e0f00,于是,当n
时
12
x(n) 000033
21
即同胞对是(A,AA)型的概率是,是(a,aa)型的概率是。
33
习 题 十 六
1. (汉诺塔问题)n个大小不同的圆盘依其半径大小依次套在桩A上,大的在下,小的在上。现要将此n个盘移到空桩B或C上,但要求一次只能移动一个盘且移动过程中,始终保持大盘在下,小盘在上。移动过程中桩A也可利用。设移动n个盘的次数为an,试建立关于an的差分方程,并求an的通项公式。
2. 设第一月初有雌雄各一的一对小兔。假定两月后长成成兔,同时(即第三月)开始每月初产雌雄各一的一对小兔,新增小兔也按此规律繁殖。设第n月末共有Fn对兔子,试建立关于Fn的差分方程,并求Fn的通项公式。
3. 在常染色体遗传的问题中,假设植物总是和基因型是Aa的植物结合。求在第n代中,基因型为AA,Aa和aa的植物的百分率,并求当n趋于无穷大时,基因型分布的极限。
T
-196-