sir传染病模型的MATLAB实现

有无大佬可以详细讲述一下MATLAB怎么跑sir传染病模型啊。数模小白，不知道关于传染病的数据应该插在哪行代码里面。。。谢谢各位大佬了

我这里有个对SEIR模型简单的注释

clear;clc;

N = 10000;                                                                  %人口总数
E = 0;                                                                      %潜伏者
I = 1;                                                                      %传染者
S = N - I;                                                                  %易感者
R = 0;                                                                      %康复者

r = 20;                                                                     %感染者接触易感者的人数
B = 0.03;                                                                   %传染概率
a = 0.1;                                                                    %潜伏者转化为感染者概率
y = 0.1;                                                                    %康复概率

T= 1:140; %时间140天
%进行迭代                                       
for idx = 1:length(T)-1
    S(idx+1) = S(idx) - r*B*S(idx)*I(idx)/N;
    E(idx+1) = E(idx) + r*B*S(idx)*I(idx)/N-a*E(idx);
    I(idx+1) = I(idx) + a*E(idx) - y*I(idx);
    R(idx+1) = R(idx) + y*I(idx);
end
%画图
figure
plot(T,S,T,E,T,I,T,R);
grid on;
xlabel('天');ylabel('人数')
legend('易感者','潜伏者','传染者','康复者')
title('第1版基础seir模型')

我专栏有些到挺多关于传染病模型的文章，希望对你有帮助，望采纳
如果还有不懂的可以私信