传染病SIS模型matlab代码

以下是一个简单的 SIS 模型的 MATLAB 代码示例：

% 初始化参数
beta = 0.005; % 传染率
gamma = 0.2; % 恢复率
N = 1000; % 总人口数
I0 = 10; % 初始感染人数
S0 = N - I0; % 初始易感人数
T = 100; % 模拟时间
tspan = 0:T;
y0 = [S0, I0];

% 定义SIS模型的ODE方程
f = @(t, y) [-beta * y(1) * y(2) / N + gamma * y(2), beta * y(1) * y(2) / N - gamma * y(2)];

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

% 绘制SIS模型的结果
figure;
hold on;
plot(t, y(:, 1), 'b-', 'LineWidth', 2);
plot(t, y(:, 2), 'r-', 'LineWidth', 2);
legend('易感者', '感染者');
xlabel('时间');
ylabel('人数');
title('SIS模型');

在这个示例中，我们使用ode45函数求解SIS模型的ODE方程，然后绘制易感者和感染者的人数随时间的变化情况。