diff --git a/SIRC.m b/SIRC.m index ba0d994987a94f7dca5d232d967639af0e9b76b4..7443c2425e015d086bde4bde50b80018fc3e6333 100644 --- a/SIRC.m +++ b/SIRC.m @@ -32,5 +32,20 @@ legend('S', 'I', 'R', 'C'); %% Add noise +noise = wgn(size(t,1), 4, 0, 42); +y_noise = y + noise; + +figure(2) +plot(t, y_noise); +xlabel('Time(days)'); +ylabel('Number of individuals'); +legend('S', 'I', 'R', 'C'); +[t_noisy, y_noisy] = ode45(@(t, y) noisy_deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q), [0 time], y0); + +figure(3) +plot(t_noisy, y_noisy); +xlabel('Time(days)'); +ylabel('Number of individuals'); +legend('S', 'I', 'R', 'C');