matlab编程传染病模型拟合（5个代码

这里给出一个 SIR 模型的传染病拟合代码，包含了数据的导入、模型的设定、参数的拟合和结果的可视化：

% 导入数据
data = readmatrix('data.csv');

% 设置模型
function dy = SIR(t,y,beta,gamma)
    dy = zeros(3,1);
    dy(1) = -beta*y(1)*y(2);
    dy(2) = beta*y(1)*y(2) - gamma*y(2);
    dy(3) = gamma*y(2);
end

% 拟合参数
params0 = [0.1, 0.1];
lb = [0, 0];
ub = [1, 1];
params = lsqcurvefit(@SIR, params0, data(:,1), data(:,2:4), lb, ub);

% 可视化结果
tspan = data(:,1);
y0 = data(1,2:4);
[t,y] = ode45(@(t,y) SIR(t,y,params(1),params(2)), tspan, y0);

figure
plot(t,y(:,1),'r-',t,y(:,2),'g-',t,y(:,3),'b-',data(:,1),data(:,2),'ro',data(:,1),data(:,3),'go',data(:,1),data(:,4),'bo')
legend('Susceptible','Infected','Recovered')
xlabel('Time (days)')
ylabel('Number of individuals')

其中，数据文件 `data.csv` 包含了每天的感染人数、治愈人数和死亡人数，格式如下：

Day,Infected,Recovered,Dead
1,100,0,0
2,150,10,0
3,200,20,5
4,250,30,10
...