求动力学模型的特征矩阵、特征方程【matlab工具集_03】

最新推荐文章于 2022-07-20 21:09:46 发布

歪卜巴比

最新推荐文章于 2022-07-20 21:09:46 发布

阅读量1.4k

点赞数

分类专栏： matlab(基础)使用技巧文章标签： matlab 线性代数特征值分解动力学

本文链接：https://blog.csdn.net/l1234x/article/details/120161661

版权

这篇博客介绍了如何利用matlab对SIVRS传染病动力学模型进行特征矩阵和特征方程的计算。通过matlab程序得出结果，并提到特征矩阵基于雅可比矩阵构造，使用collect函数整合同类项。同时提到了求解雅可比矩阵的函数参考。

摘要由CSDN通过智能技术生成

这里我们以如下SIVRS传染病动力学模型为例,其动力学方程如下所示：

(模型的详细介绍可参考几种经典病毒动力学模型【基于matlab的动力学模型学习笔记_3】_L1234X的博客-CSDN博客)

计算此动力学模型的matlab程序为：

%求特征方程
syms x1 x2 x3 x4 y%定义特征值
x=[x1,x2,x3,x4]
N=1;
a=0.6;%易感者转为感染者系数（感染系数）
b=0.8;%感染者转为免疫者系数（康复系数）
c=0.1;%感染者内部变异 
d=0.1;%感染者外部变异   
e=0.05;%死亡系数 
h=0.2;%免疫者转化为易感者概率
g=0.05;%初始出生率
k=2;%网络平均度（每个个体平均接触其它个体数）
f1=[g*N-a*k*x1*x2+h*x4-e*x1,a*k*x1*x2-d*k*x2*x3-(c+b+e)*x2,d*k*x2*x3+c*x2-e*x3,b*x2-h*x4-e*x4]%定义函数，以矩阵形式展示
j=jacobian(f1,x)%计算雅可比矩阵
n=[1 1 1 1];%构造单位向量 
D=y*diag(n)-j%构造得特征矩阵
m=det(D)%计算特征矩阵的值
k=collect(m,y)%得特征方程