【疫情模型】基於matalb改進的SEIR模型【含Matlab原始碼 667期】
阿新 • • 發佈:2021-06-24
一、簡介
基於SEIR模型改進後的傳染病模型
二、原始碼
clear;clc; %-------------------------------------------------------------------------- % 引數設定 %-------------------------------------------------------------------------- N = 10000; E = 10; I = 0; R = 0; UE = 0; Q = 0; D = 0; S = N-I-E-R-UE-Q-D; R0 = 3; a = 0.5; b1 = 0.05; b2 = 0.95; r1=0.90; e=0.01; o1=0.96; r2=0.10; o2=0.4; y = 0.004; T = 1:100; for idx = 1:length(T)-1 S(idx+1) = S(idx)-R0*a*S(idx)*E(idx)/N(1); E(idx+1) = E(idx)+R0*a*S(idx)*E(idx)/N(1)-b2*E(idx)-b1*E(idx); UE(idx+1) = UE(idx)+b1*E(idx)-r1*UE(idx)-r2*UE(idx); Q(idx+1) = Q(idx)+b2*E(idx)+r1*UE(idx)-e*Q(idx)-o1*Q(idx); I(idx+1) = I(idx)+r2*UE(idx)+e*Q(idx)-o2*I(idx)-y*I(idx); R(idx+1) = R(idx)+o1*Q(idx)+o2*I(idx); D(idx+1) = D(idx)+y*I(idx); end
三、執行結果
四、備註
版本:2014a
完整程式碼或代寫加1564658423