SIR传染病模型matlab代码

以下是基于SIR模型的传染病模型的Matlab代码，供参考：

% 设置参数
beta = 0.5; % 感染概率
gamma = 0.1; % 每天康复概率
N = 1000; % 总人口数
I0 = 1; % 初始感染人数
S0 = N - I0; % 初始易感人数
R0 = 0; % 初始康复人数
tspan = [0 200]; % 时间跨度
y0 = [S0 I0 R0]; % 初始状态

% 定义ODE函数
f = @(t,y) [-beta*y(1)*y(2)/N; beta*y(1)*y(2)/N-gamma*y(2); gamma*y(2)];

% 解ODE方程
[t,y] = ode45(f,tspan,y0);

% 绘制图像
plot(t,y(:,1),'b',t,y(:,2),'r',t,y(:,3),'g');
title('基于SIR模型的传染病模型');
xlabel('时间');
ylabel('人数');
legend('易感人数','感染人数','康复人数');