用Matlab写SIRP传染病模型

SIRP模型是针对人口传染病的流行病学模型，其中S代表易感人群，I代表感染人群，R代表恢复或死亡人群，P代表人口总数。以下是用Matlab实现SIRP模型的代码：

% 设置初始参数
N = 1000000; % 总人口数
I0 = 1; % 初始感染人数
R0 = 0; % 初始恢复人数
S0 = N - I0 - R0; % 初始易感人数
P = N; % 总人口数

beta = 0.2; % 感染率
gamma = 0.1; % 恢复率
delta = 0.01; % 死亡率

% 定义时间范围
tspan = [0 100];

% 定义初始状态
y0 = [S0; I0; R0; 0];

% 定义ODE方程
f = @(t, y) [-beta * y(1) * y(2)/P; beta * y(1) * y(2)/P - gamma * y(2) - delta * y(2); gamma * y(2); delta * y(2)];

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

% 绘制图形
figure;
plot(t, y(:,1), '-b', 'LineWidth', 2);
hold on;
plot(t, y(:,2), '-r', 'LineWidth', 2);
hold on;
plot(t, y(:,3), '-g', 'LineWidth', 2);
hold on;
plot(t, y(:,4), '-m', 'LineWidth', 2);
legend('易感人群', '感染人群', '恢复人群', '死亡人群');
title('SIRP模型');
xlabel('时间');
ylabel('人数');

在这个代码中，我们首先设定了总人口数、初始感染人数和恢复人数，从而得到初始易感人数。接下来，我们设置了感染率、恢复率和死亡率。我们还需要定义时间范围和初始状态，并使用ode45函数求解ODE。最后，我们绘制了易感人群、感染人群、恢复人群和死亡人群随时间变化的曲线。