si传染病模型matlab代码

以下是SI模型的MATLAB代码示例：

% 设置模型参数
beta = 0.5; % 传染率
N = 1000; % 总人口数
I0 = 1; % 初始感染人数
T = 100; % 模拟时间

% 初始化模型状态
S = N - I0; % 初始易感人数
I = I0; % 初始感染人数

% 运行模拟
for t = 1:T
    % 计算新感染人数
    newI = round(beta * S * I / N);
    
    % 更新状态变量
    S = S - newI;
    I = I + newI;
    
    % 绘制图形
    plot(t, S, 'bo', t, I, 'ro');
    xlabel('时间');
    ylabel('人数');
    legend('易感者', '感染者');
    drawnow;
end

该代码使用SI模型来模拟传染病的传播。其中，beta是传染率，N是总人口数，I0是初始感染人数，T是模拟时间。代码中使用for循环来迭代模型状态，每次计算新感染人数，更新状态变量，并绘制图形。最后，运行模拟并输出结果。