传染病模型

传染病模型中的符号表示


SI模型(艾滋传染模型)


matlab 复制代码
%% 直接求微分方程的解析解
dsolve('Dx1 = -0.1 * x1 * x2 / 1000', 'Dx2 = 0.1 * x1 * x2 / 1000','x1(0) = 999, x2(0) = 1', 't');

%% 根据S + I = N做一个化简

x1 = dsolve('Dx1 = -0.1 * x1 * (1000 - x1) / 1000', 'x1(0) = 999', 't');
x2 = 1000 - x1;

figure(1)
fplot(x1, [0 200], 'b')
hold on
fplot(x2, [0 200], 'r')
legend('易感者', '感染者')

%% 这道题目可以求出解析解,但是后面的大部分题目都是求不出解析解这里我们再熟悉一遍数值解的求解过程

clc; clear
global TOTAL_N
TOTAL_N = 1000;
i0 = 1;
s0 = 999;
[t, x] = ode45('fun1', [1:200], [s0, i0]);
plot(t, x(:, 1), 'r-*');
hold on
plot(t, x(:, 2), 'b-+')'
legend('易感者', '感染者')

matlab 复制代码
function dx = fun1(t, x)
    global TOTAL_N
    beta = 0.1;
    dx = zeros(2, 1);
    dx(1) = - beta * x(1) * x(2) / TOTAL_N;
    dx(2) = beta * x(1) * x(2) / TOTAL_N;
end

SIS模型(普通传染病模型 )


maatlab 复制代码
%% 
clc;clear
global TOTAL_N
TOTAL_N = 1000;
global alfa;
alfa = 0.06;
global beta
beta = 0.1;
i0 = 1;
s0 = 999;
[t, x] = ode45('fun1', [1:500], [s0, i0]);
plot(t, x(:, 1), 'r-*');
hold on
plot(t, x(:, 2), 'b-+')'
legend('易感者', '感染者')

SIR模型


这里关于总人数毫无疑问 N = S + I + R但是在传染过程中,由于康复者已经有抗体且不会再被感染,所以这里有效人群就不能再把R计算在内了,N' = S + I

matlab 复制代码
%%
clc; clear all;
N = 1000;
i0 = 1;
s0 = 999;
r0 = 0;
[t, x] = ode45('fun1', [1:500], [s0 i0 0]);
x = round(x);
figure(1);
plot(t, x(:, 1), 'r-');
hold on 
plot(t, x(:, 2), 'b-');
hold on
plot(t, x(:, 3), 'g-');
hold on
legend('易感者S','感染者I','康复者R')
matlab 复制代码
function dx = fun1(t, x)
   beta = 0.1;
   gamma = 0.02;
   % x(1)表示S,x(2)表示I, x(3)表示R
   dx = zeros(3,1);
   C = x(1) + x(2);
   dx(1) = - beta * x(1) * x(2) / C;
   dx(2) =  beta * x(1) * x(2) / C - gamma * x(2);
   dx(3) = gamma * x(2);
end

对SIR模型的拓展

matlab 复制代码
function dx = fun1(t, x)
    beta = 0.1;
   gamma = 0.02;
   if t > 100
        gamma = gamma * 10;
   end
   % x(1)表示S,x(2)表示I, x(3)表示R
   dx = zeros(3,1);
   C = x(1) + x(2);
   dx(1) = - beta * x(1) * x(2) / C;
   dx(2) =  beta * x(1) * x(2) / C - gamma * x(2);
   dx(3) = gamma * x(2);
end

拓展2:考虑死亡率

SIRS模型


相关推荐
数小模.13 小时前
我与数学建模之顺遂!
数学建模
明月看潮生14 小时前
青少年编程与数学 02-014 高中数学知识点 06课题、数学建模与探究活动
数学建模·青少年编程·数据科学·高中数学·编程与数学
HR Zhou1 天前
群体智能优化算法-算术优化算法(Arithmetic Optimization Algorithm, AOA,含Matlab源代码)
人工智能·算法·数学建模·matlab·优化·智能优化算法
青橘MATLAB学习1 天前
钢管下料问题:基于Matlab的优化求解与实践
开发语言·数学建模·matlab·钢管下料
JINX的诅咒4 天前
CORDIC算法:三角函数的硬件加速革命——从数学原理到FPGA实现的超高效计算方案
算法·数学建模·fpga开发·架构·信号处理·硬件加速器
2501_906801205 天前
BY组态-低代码web可视化组件
前端·物联网·低代码·数学建模·前端框架
2501_906801485 天前
BY组态-低代码web可视化组件
前端·物联网·低代码·数学建模·编辑器·web
电科_银尘6 天前
【Matlab】-- 基于MATLAB的美赛常用多种算法
算法·数学建模·matlab
烟锁池塘柳07 天前
【数学建模】动态规划算法(Dynamic Programming,简称DP)详解与应用
算法·数学建模·动态规划
HR Zhou8 天前
群体智能优化算法-大猩猩部落优化算法(Gorilla Troops Optimizer, GTO,含Matlab源代码)
算法·机器学习·数学建模·matlab·群体智能优化