MATLAB中sir问题

如何实现以下函数图像
在 Heaviside 函数中，本文设置在 t=20 时加入宣传和防控工作。按照 Heaviside 函数的定义[12]，当 0<t<20 时，还没有实行宣传与防控工作，则宣传与防控工作的效用为 0，即 Heaviside 函数值为 0，即 T=0；当 t>20 时，由于已经实行宣传与防控工作，产生了宣传与防控工作的效用，即 Heav－iside 函数值为 1，即 T=1。此模型设置的参数与 SIR 传染病模型仿真的第一组参数是相同的，并且在 t=20 时，加入宣传防控系数 a=0.2 作用于 SIR 传染病模型，病人人数快速下降，以此与 SIR 传染病模型形成对比，更加符合生活中的实际情况，并更有利于实际应用。

用差分方程，设置好传播系数和时间间隔就行