数学建模社区-数学中国
标题:
numpy 实现Lesile 模型
[打印本页]
作者:
2744557306
时间:
2024-3-17 16:43
标题:
numpy 实现Lesile 模型
Leslie 模型是一种经典的人口增长模型，用于描述一个有性别、有年龄结构的人口在不同年龄段的增长和变化。它是由生态学家 Patrick H. Leslie 在 1945 年提出的，被广泛应用于人口生态学和种群生态学领域。
1 `% A4 o# C8 l2 Y# [& @, O% }
Leslie 模型通常假设以下几点：
$ D' U3 n Q& s: z% l) u
: D! L S, ?2 K, D, g
1.人口结构： 人口被划分为不同的年龄组或年龄段。
5 V% }& `# N/ x) g6 V
2.增长率： 每个年龄段的人口增长率可以不同。
5 b) }! W3 T1 d O
3.繁殖率： 不同年龄段的个体具有不同的繁殖率。
- H$ }& V9 ?7 b' S- W/ W
4.死亡率： 不同年龄段的个体具有不同的死亡率。
6 N/ K7 D4 ^# y5 |
Leslie 模型的主要应用包括人口预测、资源管理、生态系统稳定性分析等。通过调整繁殖率和存活率等参数，可以模拟不同情景下人口的增长和变化趋势，为人口政策和资源规划提供科学依据。
+ s) z; r0 |& r2 \
) p9 u9 p% g" h- L. y) `" [
接线来我们来实现Leslie模型
" @- h/ Q* Y6 z7 f, b7 [! K
1.定义数据：
4 O) a2 M u5 k1 H5 i X1 m
; R9 y( A% b. H+ ^5 g4 i/ {
alpha = np.array([0, 4, 3])
^6 a5 y" Y' l1 f, x0 w
beta = np.array([.5, .25])
/ b1 D' h' y! g( j5 u8 f8 n
% a8 {" \: g, n- {4 Y; m( y9 x
L = np.zeros((len(alpha), len(alpha)))
% v$ ~$ |# y }6 @5 p5 ?1 ]) d
L[0, :] = alpha
) _, h% k/ s$ \" l
for index, each in enumerate(beta):
4 Z5 s; a" U& s
L[index+1, index] = each
! U$ O* t( E8 z" ^9 w* y& S" G* C
4 C9 L* k* M2 |" e
x0 = np.array([1000, 1000, 1000])
7 `3 ~6 f0 o6 w7 B8 ~1 _0 i" }$ \5 m
% } t, c: _, S# B) L
- s/ ?, Q ]4 b, y8 ]3 ?
2.定义了初始条件 alpha、beta 和初始状态向量 x0。
, w6 K0 ~ z) R
3.创建了一个转移矩阵 L，其中第一行为 alpha，其余为 beta。
% \5 [; v) b* J. B9 o
) `' g: I1 W( q
. Y; q+ O; x: _! {. N
4.迭代求解：
" ~+ e. R) I4 H, Q
, d% p6 N7 y. J4 V1 b* D
x_list = [x0]
1 R% u6 z6 W+ O" l
for i in range(5):
! \# S: P8 Y" d; F( u0 i
x_list.append(L.dot(x_list[-1].T))
% p8 `; H" C8 m- @7 {
3 g! H, J% O8 i& S( {9 G9 O+ c
2 k$ Z# L* c$ D: K; X
5.使用迭代方法计算了未来 5 年的状态，并将结果存储在 x_list 中。
8 d' h0 o- ]- k" X' O, t
5 J- t# a N, s0 c
0 K9 E. `: [, b% ~
6.打印结果：
. _3 m1 P# }; D- L( l; {. H5 l
% R' B) s+ r# j! _* S
for index, x in enumerate(x_list):
9 _6 z+ C5 F" Q/ g
print('the {}th year: 1y: {}, 2y: {}, 3y: {}'.format(
) N/ ~ T# E& ^" v3 U' I* f
index, x[0], x[1], x[2]))
/ {. j5 _% m0 y
% V# u5 [& }! v% f
/ T+ M r. S6 M2 p4 M8 C
7.打印了每一年的状态向量。
9 [0 o. I6 [2 S3 x
2 q& K% Q. H: w# H* S
" r" b( q) ~: v6 c2 H2 ~" \
8.绘制图形：
; r( t I- Y7 D; E
: x5 P/ p$ O$ w, ]* v% D0 u! E
year_list = np.array(range(len(x_list)))
4 d$ ^/ U2 I* K5 Z
y1,y2,y3 = [],[],[]
! f4 ?6 J, P O) i+ @
for year in year_list:
: u: M' P6 ?9 d) Q; y: N
y1.append(x_list[year][0])
& Z v/ `+ d" F: n- p k) t1 n
y2.append(x_list[year][1])
4 U0 }# Y) F. l) r
y3.append(x_list[year][2])
- e; x2 i u# M0 U
. y5 C& D! V9 ~2 l8 C6 A; z
w = .3
! k" v% T, s( z. p" v/ [
plt.bar(year_list-w,y1, width=w, label='1y')
/ J7 A; T; \+ _4 V! F* I( h
plt.bar(year_list,y2, width=w, label='2y')
" L1 k& x P+ p) s( R( A+ j! v8 _
plt.bar(year_list+w,y3, width=w, label='3y')
6 Z6 w3 J$ {& T. j: x- e. a' Y y1 Q) l
plt.legend()
' k# @; c0 S& T
2 I0 F% G4 ^4 N7 C* p
+ P1 `; d8 j8 K! K: F' @7 y1 _
9.使用 Matplotlib 绘制了每一年不同状态的柱状图，并添加了图例。
( [7 L$ Z& k: m$ ]# `' \* W
' M) F9 e! R7 i: C0 Y: H0 W
这段代码的主要目的是模拟了基于给定的转移矩阵和初始状态的系统状态随时间的演变，并通过图形化展示了模拟结果。
% K& r" H" [' {+ Y3 |3 X
# C: u+ B( @( n
: A1 \ c/ F) ^. h0 }8 F
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5