Matlab 为SIR疾病模型绘制多条线

我需要关于为不同参数绘制多条线的帮助，例如\beta为不同的值

（\beta=0.1、0.2、0.4等）

，图形将用于

dy_2/dt

，这意味着

仅绘制（T，y_2（：，2）

或

绘制（T，y_2（：，3））

然后，针对不同的\beta值，绘制时间与感染的曲线图。我尝试了以下代码：

函数euleristenga（）
全部清除；clc；
ρ=0.7；λ=0.4；β=0.9；θ=0.1；
θ=0.2；θ=0.1；α=0.33；ω=0.4；磅平方英寸=0.3；
R0=ρ*β/（θ1+λ）
选项=odeset（'RelTol'，1e-4，'AbsTol'，[1e-4 1e-4 1e-4 1e-4 1e-4]）；
[T，Y]=ode45（@SIRmodel[0365]，[0.9,0.7,0.5,0.15,0.35]，期权）；
图（T，Y（：，1），'g'，T，Y（：，2），'r'，T，Y（：，3），'k'，T，Y（：，4），'b'，T，Y（：，5），'m'，'Linewidth'，2）；
xlabel（'Time'，'FontSize'，14）；
ylabel（‘人口比例’、‘FontSize’、14）
图例（‘易感’、‘感染’、‘ART上的艾滋病’、‘可检测VL’、‘不可检测VL’）
标题（“人口与时间比例的确定性解决方案，'FontSize'，14）；
网格化
开箱
函数dy=SIRmodel（t，y）
dy=零（5,1）；
dy（1）=rho*Lambda-（β*y（2）+Lambdaα*（y（3）+y（4）+y（5））*y（1）；
dy（2）=（1-rho）*Lambda+（β*y（1）-theta1 Lambda+alpha*（y（3）+y（4）+y（5））*y（2）；
dy（3）=θ1*y（2）+θ2*y（4）+θ3*y（5）-（α+ω+λα*（y（3）+y（4）+y（5））*y（3）；
dy（4）=psi*ω*y（3）-（α+θ2+λα*（y（3）+y（4）+y（5））*y（4）；
dy（5）=（1-psi）*omega*y（3）-（αθ3+Lambdaα*（y（3）+y（4）+y（5））*y（5）；
结束
结束

我希望图表如下图所示：

我需要改变\beta的值以不同的速率变化

这是我的代码生成的数字

---

能否添加您迄今为止尝试过的代码，以便我们更好地了解您的问题所在？我会将所有线的值存储在一个矩阵中，并在循环中逐个绘制。通过指定将图形的hold属性设置为“on”，可以在图形上绘制多条线。您可以使用命令

按住

@Thombou执行此操作。我已附上代码，请您也附上此代码生成的图形。这将给出一个完全重载的绘图。您有5个不同函数的图形，您希望将此函数集的5个变体添加到绘图中。有5个绘图，一个f是有意义的或每个函数及其变体。