From 178f62b439c3c61af1225f9abe6bd12e2734d5c9 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rivny=C3=A1k=20T=C3=ADmea?= <rivnyak.timea@hallgato.ppke.hu> Date: Tue, 1 Jun 2021 21:01:06 +0000 Subject: [PATCH] Add noise to the system --- SIRC.m | 15 +++++++++++++++ 1 file changed, 15 insertions(+) diff --git a/SIRC.m b/SIRC.m index ba0d994..7443c24 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'); -- GitLab