大象群落的穩(wěn)定發(fā)展

摘要

本文根據(jù)非洲某國的國家公園近兩年內(nèi)從公園運出的大象的大致年齡和性

別的統(tǒng)計情況，探討大象的合理的存活率并推測當(dāng)前的年齡結(jié)構(gòu)，針對不同情況

給出如何進(jìn)行避孕注射以達(dá)到控制大象數(shù)量的目的。

首先，充分利用給出的近兩年來運出的大象的數(shù)量與性別統(tǒng)計表，分析近兩年來

的大象群落的情況，建立一個線性方程組的數(shù)學(xué)模型，通過求解方程組得到年齡

在2歲到60歲之間的大象的總數(shù)，并且求出了存活率為：98.9718%;因為假設(shè)

公園內(nèi)2歲到60歲之間的大象占總大象的比例等于運出的2歲到60歲之間的

大象占總移出大象的比例，所以通過一些比例之間的關(guān)系得到這個大象群落的當(dāng)

的的年齡結(jié)構(gòu)（見表1）。

然后，建立一個按年齡分組的種群增長的差分方程模型，運用第一問求出的各年

齡段大象的存活率以與繁殖率，求解當(dāng)前大象群落對應(yīng)的Leslie矩陣的特征根,

發(fā)現(xiàn)該特征根大于1,根據(jù)Leslie矩陣的穩(wěn)定性理論知道：如果不進(jìn)行避孕注射

該大象種群將無限增長（如果環(huán)境允許）；據(jù)此，利用Leslie矩陣穩(wěn)定的充要條

件求出應(yīng)該保持多大的繁殖率才能使種群保持穩(wěn)定，求解的主要思路是：特征根

取為1、把繁殖率當(dāng)成未知數(shù)，將此時的各年齡段的存活率代入方程VI,求解這

個以繁殖率為未知數(shù)的方程可以得到要使種群保持穩(wěn)定繁殖率的取值；根據(jù)需要

避孕掉母象所生的幼象的數(shù)目等于注射避孕藥使得母象沒有繁殖幼象的數(shù)目這

一條件建立一個方程，最后求得每年注射避孕藥的母象頭數(shù)為：1393（頭）。

最后，假設(shè)被轉(zhuǎn)移的大象只考慮處于1—60歲之間，這樣可以認(rèn)為轉(zhuǎn)移后的大象

看成每年多死了這么多頭大象，即意味著死亡率將增加，存活率將減少；仍然按

照解決第二問的模型，只需將此時不同的各年齡段大象的存活率代入那個以繁殖

率為未知數(shù)的方程（方程VI）,求出應(yīng)該保持多大的繁殖率才能使種群保持穩(wěn)定。

考慮到求解的數(shù)據(jù)比較多，采取計算機模擬的方法來確定移出大象后所需要進(jìn)行

避孕的母象頭數(shù)（見表2）,為了檢驗計算機模擬的正確性，用理論去驗證。模

擬的思路方法見計算機模擬流程圖一圖2。

關(guān)鍵詞：關(guān)鍵字：線性方程組、差分方程模型、Leslie矩陣、計算機模擬

問題重述

1.位于非洲某國的國家公園中棲息著近11000頭大象。管理者要求有一個健康

穩(wěn)定的環(huán)境以便維持這個1100。頭大象的穩(wěn)定群落。管理者逐年統(tǒng)計了大象

的數(shù)量，發(fā)現(xiàn)在過去的2。年中，整個大象群經(jīng)過一些偷獵槍殺以與轉(zhuǎn)移到外

地還能保持在1100。頭的數(shù)量，而其中每年大約有近600頭到800頭是被

轉(zhuǎn)移的。由于近年來，偷獵被禁止，而且每年要轉(zhuǎn)移這些大象也比較困難，現(xiàn)

決定采取避孕注射法以維持大象數(shù)量的平衡。我們已知此公園近兩年內(nèi)從這

個地區(qū)運出的大象的大致年齡和性別的統(tǒng)計。根據(jù)這些信息我們需要解決以

下問題：

2.探討年齡在2歲到60歲之間的象的合理的存活率的模型，推測這個大象群落

的當(dāng)前的年齡結(jié)構(gòu)。

3.估計每年有多少母象要注射避孕藥，可以使象群固定在1100。頭左右。這里

不免有些不確定性，是否能估計這種不確定性的影響。

假如每年轉(zhuǎn)移50至300頭象到別處，那么上面的避孕措施將可以有怎樣的改

變？

問題假設(shè)

L假設(shè)大象的性別比近似認(rèn)為1:1,并且采用措施維持這個性別比；

2.假設(shè)母象可以懷孕的年齡為11歲一60歲、最高年齡為70歲，70歲的死亡率

為100%,并且61—70歲的大象的頭數(shù)呈線性遞減；

3.假設(shè)大象在各年齡段中的分布率不變，即年齡結(jié)構(gòu)不變，并采用各種措施維持

這一結(jié)構(gòu)；

4、假設(shè)被轉(zhuǎn)移的大象只考慮處于1一60歲之間，轉(zhuǎn)移后的大象看成每年多死了

這么多頭大象；

5.假設(shè)。歲大象能夠活到1歲的比例為75%;

符號說明

□：表示一年中大象的頭數(shù)（i=0表示0歲大象的頭數(shù)，i=1表示1-60歲大象頭數(shù),

i=2表示61-70歲大象的頭數(shù)）；

□：表示存活率（口表示0歲大象的存活率，□表示1一60歲大象的存活率，□表示

61歲一70歲大象的存活率）；

□：表示時段k第i年齡組的大象數(shù)量；

□:第i年齡組每個（母象）個體在1個時段內(nèi)平均繁殖的數(shù)量；

□:第i年齡組的存活率；

□：Leslie矩陣；

U：LI矩陣的那個唯一止特征根;

□:表示移出大象的頭數(shù)；

問題分析

對于問題一，利用給出的近兩年來運出的大象的數(shù)量與性別統(tǒng)計表，可以分析近

兩年來的大象群落的情況，比如移出的各個年齡段的大象占移出的總的大象的頭

數(shù)的比例是多少，還可以根據(jù)兩年移出大象后大象總數(shù)都是11000來建立方程,

用于求解存活率。

對于問題二，因為考慮的是公園在未來很長一段時間的大象種群控制問題，所以

可以建立一個按年齡分組的種群增長的差分方程模型，根據(jù)差分方程的Leslie

矩陣的特征根，結(jié)合Leslie矩陣的穩(wěn)定性理論對當(dāng)前大象種群的情況進(jìn)行分析。

為了保持大象種群的穩(wěn)定，必須使得Leslie矩陣的最大特征根為1,而這樣，特

征根取為1.把繁殖率當(dāng)成未知數(shù)，將此時的各年齡段的存活率代入方程特征方

程，求解這個以繁殖率為未知數(shù)的方程可以得到要使種群保持穩(wěn)定繁殖率的取

值；根據(jù)需要避孕掉母象所生的幼象的數(shù)目等于注射避孕藥使得母象沒有繁殖幼

象的數(shù)目這一條件建立方程來求解應(yīng)該對多少頭母象進(jìn)行避孕。

對于問題三，由于假設(shè)被轉(zhuǎn)移的大象只考慮處于1—60歲之間，故可以認(rèn)為轉(zhuǎn)移

后的大象看成每年多死了這么多頭大象，即意味著死亡率將增加，存活率將減

少。按照解決第二問的模型，只需將此時不同的各年齡段大象的存活率代替原來

的存活率，就可以求出此時應(yīng)該保持多大的繁殖率才能使種群保持穩(wěn)定。為了方

便，可用采用計算機模擬的方法來確定移出的大象在哪個年齡段，考慮到計算機

模擬的不確定性，可以對模擬結(jié)果進(jìn)行檢驗。

探討大象的存活率和當(dāng)前大象的年齡結(jié)構(gòu)

下面將根據(jù)給出的近兩年來運出的大象的數(shù)量與性別統(tǒng)計表，分析近兩年來的大

象群落的情況，建立一個線性方程組數(shù)學(xué)模型，通過求解方程組得到年齡在2歲

到60歲之間的大象的存活率，并給出大象各年齡所占的比例，進(jìn)而得到這個大

象群落的當(dāng)前的年齡結(jié)構(gòu)。

1.線性方程組模型的建立

(1)首先，計算一年中大象的頭數(shù)。

大象群是山0歲，1—60歲，61歲一70歲組成，且穩(wěn)定在

11000頭。設(shè)。歲的頭數(shù)為X0,1—60歲大象頭數(shù)為XI,61歲一70歲大象頭

數(shù)為X2。所以得到第一個方程：

X0+X1+X2=11000

(D

(2)其次，考慮到前一年大象的總數(shù)等于前兩年存活卜來的大象加上新生

的幼兒再減去運出的大象數(shù)。

設(shè)。歲大象的存活率為，1—60歲大象的存活率為，61歲

-70歲大象的存活率為。則經(jīng)過一年后，新生的大象存活下來的頭數(shù)為

X0;1到60歲的大象存活下來的頭數(shù)為XI;61歲一70歲的大象能

存活下來的頭數(shù)為X2，因此得到第二個方程：

(Xox+Xixpi+X2xp?)+X0-622=l1000

(II)

聯(lián)立(I)、(II)得到方程組：

Xo+X,+X2=llOOO

'X。XPo+X,XPl+XXP2+X-622=11000

2o

2.模型的求解

根據(jù)近兩年來運出的大象的數(shù)量與性別統(tǒng)計表，得到如下分析結(jié)果：

(1)計算0歲的大象頭數(shù)

由表中統(tǒng)計，1歲一10歲的大象占1歲一60歲的大象比例為：

(67/622+169/876)/2=15.032%

所以得到：11歲―60歲能生小象的母象占1歲—60歲的大象比例

為：

(1-15.032%)x0.5=42.48%

因為能生小象的母象每3.5年生一頭小象，且雙胞胎的機會為

1.35%,相當(dāng)于每年生0.2896頭，所以。歲的大象占1歲一60歲的大象比例

為：

0.4248x0.2896=0.12303

這樣0歲的大象共有：=0.12303(III)

(2)計算61歲一70歲的大象頭數(shù)

從表中計算運出的59歲的大象占運出的總大象比率為：

(14/622+22/876)/2=0.0238

由于運出的大象都是1歲一60歲的，所以0.0238也可看為59歲的大象占

1—60歲的大象的頭數(shù)比例，得到60歲的大象占的比例為0.0238,由假設(shè)可

以知道：

61歲一70歲的大象頭數(shù)為：

X2=1/2X10X0.0238XP1XX1（IV）

61歲——70歲的大象經(jīng)過一年能存活下來的頭數(shù)為:

X2xp2=（1/2）x9x0.0238xp,xX,（V）

（3）、將（in）、（V）和（iv）兩個式子代入上面方程組（*）得：

（）.12303xX,+Xl+（l/2）x10x0.0238x/；（xX,=iKXX）

0.12303xX1x/%+X]xn+（1/2）x9x0.0238xnxX1+0.12303xX「622=11000

又由假設(shè)知道，0歲大象的存活率為=75%代入上述方程組,

然后用Mathematica解之得：

px=0.989719

X,=8864.85

再依次將、代入（HI）、（V）和（IV）求得：

X2=1044.07

Xo=1090.66

所以，。歲大象的總頭數(shù)為1091（頭）；1-60歲的大象的存活率為

98.9719%,總頭數(shù)為8865（頭）；61歲―70歲的大象頭數(shù)為1091（頭）。

把0—70歲的大象分為八個年齡段，由假設(shè)知道，各個年齡段占總

數(shù)可以用各個年齡段移出的頭數(shù)除以移出的總頭數(shù)來衡量。下面以1—1。年齡段

的大象頭數(shù)計算為例：

前一年總共移出622頭，其中1—10歲移出為67頭；前兩年總共

移出876頭，其中1—10歲移出169頭。故1—10年齡段的大象頭數(shù)可以這樣

計算：

X“=8865x[（2+繇/2]=1333（頭）

其他的年齡段用同樣的方法計算，得到如下表（附餅形圖）：

表1(大象年齡結(jié)構(gòu))

年齡01-1011-2021-3031-4041-5051-6061-70

頭數(shù)10911333177710691255188815441044

嬲101216101117149

圖1(大象年齡結(jié)構(gòu)餅圖)

3.結(jié)果分析

(1)由結(jié)果可以知道，2-60歲大象的存活率為98.9718%,這與題目給出

的大于95%是相一致的，所以可以認(rèn)為結(jié)果是合理的；

(2)從圖1可以看出，各個年齡段的大象所占的比例基本上是一樣的，21

-30歲和41-50歲的大象比例相對比較大，因為這段大象正處于年齡的黃金

時期。由此，可以認(rèn)為求出的大象年齡也是合理的。

估計每年注射避孕藥的母象頭數(shù)

為了估計每年注射避孕藥的母象頭數(shù)，首先建立一個按年齡分組的種群增長的差

分方程模型；然后用Leslie矩陣穩(wěn)定的充要條件分析如果不進(jìn)行避孕注射種群的

增K情況；最后仍然利用Leslie矩陣穩(wěn)定的充要條件求出應(yīng)該保持多大的繁殖率

才能使種群保持穩(wěn)定，進(jìn)而利用一個方程求出每年注射避孕藥的母象頭數(shù)。

1.按年齡分組的種群增長的差分方程模型的建立

記時段k第i年齡組的大象數(shù)量為，k=0,1,2,2,…n,第i年齡組的

繁殖率為，即第i年齡組每個(母象)個體在1個時段內(nèi)平均繁殖的數(shù)量，第i

年齡的存活率為，我們這里假設(shè)和不隨時段k變化，在穩(wěn)定的環(huán)境F這個

假設(shè)是合理的。和可由統(tǒng)計資料獲得。的變化規(guī)律由以下的基本事實得到:

時段k+1第1年齡組種群數(shù)量是時段k各年齡組繁殖數(shù)量之和，即

用伙+1)=才2若(幻

/=1

時段k+1笫i+1年齡組的種群數(shù)量是時段k第i年齡組存活下來的數(shù)量，即

XM(A+1)=s/j(k),z=1,2...n-\

記時段k種群按年齡組的分布向量為

x(k)=[%(k),%伏)%⑹F

由繁殖率e和存活率士構(gòu)成的矩陣為

飛a…為九b\\.........do*…％)“70

So0...000...............00...00

0...000...............00...00

00...000...............00...00

00…§900...............00...00

00...0鳥00.........00...00

L=??????...00...............00...00

??????...000.........00...00

00...000...............00...00

00...000..............$600...00

00...000...............0%...00

00...000...............00...00

00...000.........00…$690

根據(jù)Leslie矩陣的性質(zhì)可以得到如下的定理:

定理1：矩陣有唯一的正特征根，且它是單重根的，對應(yīng)正特征向量

L矩陣的其他n-1個特征根人都是滿足同工4叢=2,3,...,〃

該定理表明L矩陣的特征方程

=

4"一(4++bn_xSxS2...Sn_2X+)0

只有一個正根，并且易知，

2.如果不進(jìn)行避孕注射種群的增長情況

(1)建立Leslie矩陣

首先，由第一問的求解知道，0歲的大象的存活率為075;1-60歲大象的存

活率為0.989718;根據(jù)假設(shè)61-70歲大象頭數(shù)是線性遞減的，而且到7。歲

所有的大象都死完了，所以很容易求出存活率為(1-0.1)

=0.90.989718=0.8907;11一60歲大象的繁殖率為0.1448。

然后根據(jù)上面的矩陣L建立起如果不進(jìn)行避孕注射種群的增長的Leslie矩陣如下

所示：

■00...00.14480.1448........0.14480...00'

0.750...000........00...00

00.9897...000........00...00

00...000........00...0()

00...0.989700........00...00

00...00.98970........00...00

^71x71=?????????000.9897........00...00這是一

.????????000........00...00

00...000........00...00

00...000........0.98970...00

00...000........00.8907...00

00...000........00...00

00...000........00...0.89070

個71，71的矩陣。

(2)討論的特征根，分析種群增長規(guī)律

用Matlab軟件求得特征根為R=1.0481,根據(jù)定理1知道，如果不進(jìn)行避孕注

射，該大象種群將無限增長下去(如果環(huán)境允許)，所以要進(jìn)行避孕注射。

3.求出每年注射避孕藥的母象頭數(shù)

根據(jù)Leslie矩陣的性質(zhì)知道，要保持種群穩(wěn)定，必須使得特征根r=l,

即使得下面式子成立：

=

%+睜0+力2sosI+???+S2,,'Sn-2^nS0SiS2,,'Sn-\1(VI)

具體到本題來就是使得如下成立：

)O24y

bo-S-S(1H-S-l-SH-...H-S)=1

解這個方程求出要保持大象種群的穩(wěn)定，繁殖率應(yīng)該為=0.0377

保持大象種群數(shù)量不變的繁殖率b0與沒采取避孕時的繁殖率b有一定的差

距，所以需要避孕掉具有(b-b0)繁殖率母象所生的幼象。假設(shè)每年要避孕

頭大象，由于一次注射可以使得一頭成熟的母象在兩年內(nèi)不會受孕，所以每年實

際上共有2頭大象處于避孕期。

這樣根據(jù)需要避孕掉具有(b-b0)繁殖率母象所生的幼象的數(shù)目等

于注射避孕藥使得母象沒有繁殖幼象的數(shù)目這個條件得到一個方程：

(1/2)x85%xX]x(0.1448-0.0377)=2n0xO.l448

解之得n°=1393

所以每年注射避孕藥的母象頭數(shù)為：1393(頭)

4.分析不確定因素的影響

(1)最初一兩年避孕母象發(fā)情期增多，與未避孕母象產(chǎn)生競爭求偶的公象,

使部分能懷孕的母象不能懷孕。而避孕的母象每月發(fā)情一次，會擾亂了正常求偶

的母象，這樣會造成未避孕母象的繁殖率出現(xiàn)下降，避孕的母象數(shù)量應(yīng)該減少。

(2)隨著時間的增長，如果持續(xù)使用避孕藥，會使象的年齡結(jié)構(gòu)發(fā)生變化,

象的結(jié)構(gòu)呈老齡化，所以隨著時間的增K,要保證象群的穩(wěn)定，避孕藥的使用量

必定會逐年減少直至禁用。

考慮轉(zhuǎn)移大象時母象的避孕策略

被轉(zhuǎn)移的大象只考慮處于1-60歲之間，轉(zhuǎn)移后的大象看成每年多死了這么多

頭大象，即意味著死亡率將增加，存活率將減少。下面首先通過計算機模擬來確

定移出大象后所需要進(jìn)行避孕的母象頭數(shù)；然后用理論去驗證計算機模擬的正確

性。

1.通過計算機模擬確定需要進(jìn)行避孕的母象頭數(shù)

產(chǎn)生了n(可自己指定)個0~1的隨機數(shù)。具體算法如下頁圖所示。

圖2（計算機模擬流程圖）

初始化：按照各個年齡段

大象所占比重，給每個年

齡段分配一個區(qū)間，每個

產(chǎn)生n個0~1隨機數(shù)

記為向量a（i）

()<a(i)W0.148<a(i)

kl=kl+

0.346<a(i)三

經(jīng)k3=k3+

0.469va(i)k4=k4+

0.6!7<a(i)三堡k5=k5+

0.827<a(i)W1醇k6=k6+

根據(jù)kL?k6求出各年齡段

大象的存活率、運用Leslie

No

i=i+Li=n矩陣穩(wěn)定的條件求出此時

繁殖率應(yīng)該為多少才能保

持大象種群穩(wěn)定，進(jìn)而求

出需避孕的大象頭數(shù)

下面議n=100為例進(jìn)行計算機模擬。

令n=100,進(jìn)行10次計算機模擬，得到當(dāng)運出大象頭數(shù)為100頭時，要

使大象頭數(shù)穩(wěn)定在1100。頭需避孕的大象頭數(shù)如下：

第一次模擬：1164（頭）

第二次模擬：1195（頭）

第三次模擬：1192（頭）

第四次模擬：1167（頭）

第五次模擬：1206（頭）

第六次模擬：1163（頭）

第七次模擬：1207（頭）

第八次模擬：1190（頭）

第九次模擬：1173（頭）

第十次模擬：1178（頭）

10次模擬得到需避孕的大象頭數(shù)的平均數(shù)為1184頭。因此可以認(rèn)為當(dāng)運

出大象頭數(shù)為100頭時，要使大象頭數(shù)穩(wěn)定在11000頭需避孕的大象頭數(shù)為

1184頭。

同理，可以得到題目中要求的當(dāng)運出大象50-300頭要使大象頭數(shù)穩(wěn)定在

11000頭需避孕的大象頭數(shù)分別為：

表2（移出大象與對應(yīng)得避孕母象頭數(shù)）

彩

出

5060708090100110120130140150160170

頭

數(shù)

避12127125123121118116115111109103103101

孕97731244026600

頭

數(shù)

移

出18

190200210220230240250260270280290300

頭0

數(shù)

遽

學(xué)97714

932917860838822764738654653591570

頭0

數(shù)

圖3（運出大象數(shù)目與應(yīng)該避孕母象數(shù)目的關(guān)系）

1400

戴

冰

￡

<

的

S

群

綴

粕

國

名

日

濱

象

的母

避孕

應(yīng)該

0頭而

100

定在1

數(shù)穩(wěn)

象總

使大

為了

看出

可以

趨勢

變化

圖的

從該

數(shù)得

次函

合二

行擬

乘進(jìn)

小二

過最

系。通

數(shù)關(guān)

次函

成二

近似

頭數(shù)

象的

出大

與運

數(shù)目

：

關(guān)系

數(shù)y

象頭

孕母

與避

數(shù)X

象頭

出大

到移

為:

圖像

畫出

）

比圖

系對

合關(guān)

與擬

關(guān)系

數(shù)據(jù)

（原始

圖4

原始關(guān)系與擬合得到關(guān)系對比圖

1300^----------(------------1------------1-----r----------

——原始關(guān)系-

1200"-----擬合圖像

1100

o

100o

90

80

O

2.用理論去驗證計算機模擬的正確性

假設(shè)所轉(zhuǎn)移的大象是有目的的挑選的，即挑選大象時是按照第i(i=l,2,3,4,

5,6)年齡組所占比例進(jìn)行的，這也是符合情理。設(shè)每年被轉(zhuǎn)移的大象共有M頭,

1—60歲的大象共有W頭。

設(shè)第i年齡組占1—60歲的比例為，第i年齡組的存活率，則有

5lf.=(X[xMxWj-Mx叱)/(X]x叫)=(X|xI-M)/X[

則要保持大象群落穩(wěn)定，如第二問的做法有：

2

Z?X50Xx(l+5,4-52+...+)=1

當(dāng)移出大象頭數(shù)M=50時，解得b=0.0455。

又由第三問，把這個方程的換為n、0.0377換為0.0455,求得

n=1291o

所以，現(xiàn)在只需避孕1291頭，由計算機模擬得到的是1297頭，非常地接

近。同樣的道理可以驗證當(dāng)M=60,70,80…的時候也是與計算機模擬很接近的。

由此，可以說明我們用計算機模擬的方法是有效的。

模型的評價和改進(jìn)方法

（1）1.模型的優(yōu)點

（2）本文解決問題的模型都是比較簡單的，但是這并不影響得到的結(jié)果的準(zhǔn)確

性，因為這些簡單的模型都有很強的理論依據(jù)；

（3）在求解第二問的時候，充分利用Leslie矩陣穩(wěn)定性理論來求解應(yīng)該讓多少

母象進(jìn)行避孕注射，這些理論在差分方程中都是經(jīng)典的理論，經(jīng)得起許多

事實的考驗；

（4）第三問的求解中運用了計算機模擬方法來模擬移出大象屬于哪個年齡段,

這樣不僅求解方便、簡潔（只需要把算法程序?qū)懞镁涂梢缘玫浇Y(jié)果），得

到的結(jié)果與實際也更接近；

第三問用計算機模擬得到數(shù)據(jù)后，又用理論去驗證，這樣使得結(jié)果更具有說服

力；

2.模型需要改進(jìn)的地方

（1）因為假設(shè)了大象性別是嚴(yán)格地1：1關(guān)系，而實際中不一定那么地嚴(yán)格是這

樣，所以如果能夠把各個年齡段大象的性別比例分別計算，那么模型的結(jié)

果可能更接近實際；

（2）在進(jìn)行計算機模擬時，最開始的隨機數(shù)的產(chǎn)生個數(shù)只有幾十個，這幾十個

隨機數(shù)不能很好的反映各個年齡段的大象所占的比重，這樣勢必會對結(jié)果

造成一定的誤差;

參考文獻(xiàn)：

L姜啟源數(shù)學(xué)模型（第三版）[M].北京：高等教育出版社2003

2.趙靜數(shù)學(xué)建模與數(shù)學(xué)實驗（第2版）[M].北京：高等教育出版社

2003

3.冏曉陽數(shù)學(xué)實驗與Matlab[M].武漢：華中科技大學(xué)出版社2002

4.鄭諫當(dāng)代數(shù)學(xué)的若干理論與方法[M].上海：華東理工大學(xué)出版社2002

5、李尚志數(shù)學(xué)建模競賽教程[M].江蘇：江蘇教育出版社1996

6.北峰數(shù)模網(wǎng)2006-8-13

附錄：

1.計算機模擬Matlab程序代碼

function[bys,bO]=moni(n)

b=[121610121714];

c=b/sum(b);

b=cumsum(c)';

a=rand(n,1);

kl=0;k2=0;k3=0;k4=0;k5=0;k6=0;

i=l；

fori=l:n

if(a(i)v=b⑴)

kl=kl+l；

elseif(a(i)>b(l)&a(i)<=b(2))

k2=k2+l;

elseif(a⑴>M2)&a(i)v=b⑶)

k3=k3+l;

elseif(a⑴〉b⑶&a(i)<=b⑷)

k4=k4+l;

elseif(a(i)>b(4)&a(i)<=b(5))

k5=k5+l;

else(a(i)>b(5)&a(i)<=b(6))

k6=k6+l;

end

end

sl0=l-(11000*0.12*(l-0.9897)+kl)/(l1000*0.12);

s20=1-(11000*0.16*(1-0,9897)+k2)/(ll000*0.16);

s30=l-(l1000*0.l*(l-0.9897)+k3)/(l1000*0.1);

s40=l-(11000*0.12*(l-0.9897)+k4)/(l1000*0.12);

s50=l-(l1000*0.17*(l-0,9897)+k5)/(l1000*0.17);

s60=l-(l1000*0.14*(l-