%SEIR 模型
clear;clc;
%参数设置
N=60000000;
I=1;%传染者
R=0;%康复者
D=0;%死亡患者数量
E=0;%潜伏者
S=N-I;%易感染者
r=2;%接触病患的人数
a=0.125;%潜伏者患病概率
B=0.8;%感染概率
y=0.143;%康复概率
B2=0.03;%接触潜伏者
B3=0.859;%转阴率
r2=10;%潜伏者每天接触的人数
k=0.025373;
T=0:200;
% for idx =1:length(T)-1
% S(idx+1)=S(idx)-rBI(idx)S(idx)/N;%易感人数代
% E(idx+1)=E(idx)+rBS(idx)I(idx)/N-aE(idx)%潜伏者人数迭代
% I(idx+1)=I(idx)+aE(idx)-(y+k)I(idx);%患病人数迭代
% R(idx+1)=R(idx)+yI(idx);%康复人代
% D(idx+1)=D(idx)+kI(idx);%死亡患者人数迭代
% end
%修正
for idx =1:length(T)-1
%S(idx+1)=S(idx)-rBI(idx)-r2B2E(idx)+B3E(idx);%易感人数迭代
%E(idx+1)=E(idx)+rBI(idx)+r2B2E(idx)-B3E(idx);%潜伏者人数迭代
%I(idx+1)=I(idx)+aE(idx)-(y+k)I(idx);%患病人数迭代
%R(idx+1)=R(idx)+yI(idx);%康复人数迭代
%D(idx+1)=D(idx)+k*I(idx);%死亡患者人数迭代
end
%plot(T,S,T,E,T,I,T,R,T,D);
plot(T,E,'--r',T,I,'-k');
grid on;
xlabel('日期');
ylabel('人数/人');
%legend('易感者','潜伏者','潜伏者转阳','康复者','死亡者');
title('SEIR模型');
%plot(T,E,T,I,T,R,T,D);
plot(T,E,'--r',T,I,'-k');
grid on;
xlabel('时间/天');
你这里E这个变量是0吧？把它作为y来画图？