标签:xk python sp Leslie Ls print X0 Pinv 莱斯利
莱斯利Leslie种群模型 python
sympy
# 莱斯利Leslie种群模型
import numpy as np
import sympy as sp
X0 = np.array([500, 1000, 500])
L = np.array([[0, 4, 3], [0.5, 0, 0], [0, 0.25, 0]])
X1 = L @ X0
X2 = L @ X1
X3 = L @ X2
Ls = sp.Matrix([[0, 4, 3], [sp.Rational(1, 2), 0, 0], [0, sp.Rational(1, 4), 0]]) # 符号矩阵
lamda = sp.var('lamda')
p = Ls.charpoly(lamda) # 计算特征多项式 #re-A的行列式
# w11 = Ls.eigenvals()
# w22 = Ls.eigenvects()
w1 = sp.roots(p) # 特征值
w2 = Ls.eigenvals() # 直接计算特征值
v = Ls.eigenvects()
# print("特征值", w2)
# print(w1)
# print('特征向量', v)
P, D = Ls.diagonalize() # 相似对角化
Pinv = P.inv() # 求逆矩阵
Pinv = sp.simplify(Pinv)
cc = Pinv @ X0
# print(P)
# print(cc[0])
# print(w1)
# print(v)
k = sp.var('k', positive=True, integer=True)
xk = P @ (D ** k) @ Pinv @ sp.Matrix(X0)
# print(xk)
# print(xk[0])
s = sp.simplify(xk[0])
print(s.subs(k, 2).n()) # # 计算第k个时期的种群数量
标签:xk,python,sp,Leslie,Ls,print,X0,Pinv,莱斯利 来源: https://blog.csdn.net/weixin_51020254/article/details/122784983
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。