莱斯利Leslie种群模型 python sympy
2022/2/4 20:12:50
本文主要是介绍莱斯利Leslie种群模型 python sympy,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
莱斯利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个时期的种群数量
这篇关于莱斯利Leslie种群模型 python sympy的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-08有遇到过吗?同样的规则 Excel 中 比Python 结果大
- 2024-03-30开始python成长之路
- 2024-03-29python optparse
- 2024-03-29python map 函数
- 2024-03-20invalid format specifier python
- 2024-03-18pool.map python
- 2024-03-18threads in python
- 2024-03-14python Ai 应用开发基础训练,字符串,字典,文件
- 2024-03-13id3 algorithm python
- 2024-03-13sum array elements python