Rechenvorschrift

Aufbauend auf den Annahmen und Bezeichnungen zuvor können wir nun das SIR-Modell einführen: Wir betrachten weiterhin die Zeitschritte und verwenden die folgenden Schreibweisen:

Um das SIR-System modellieren zu können, führen wir zusammenfassend die folgenden Parameter ein:

Anzahl der Individuen der gesamten Population
Anteil der (anfangs) infizierten Individuen
Kontaktrate
Genesungsrate

Insbesondere die Kontaktrate sowie die Genesungsrate werden wir verwenden, um die folgenden Abhängigkeit zwischen den Gruppen zu modellieren:

  1. Für die Krankheit empfänglich Individuen (S) können sich an bereits infizierte Individuen (I) anstecken.
  2. Infizierte Individuen (I) können genesen und sind anschließend immun (R).

Diese beiden Abhängigkeit modellieren wir unter Verwendung der eingeführten Parameter folgendermaßen:

Dieses System mag auf den ersten Blick recht komplex wirken und wir sollten uns ein wenig Zeit nehmen, es genauer zu verstehen. Zunächst ist es wichtig zu verstehen, dass

die Anzahl der Individuen ist, die sich pro Zeitschritt mit der Krankheit infizieren. Wichtig dabei ist insbesondere der folgende Anteil:

Dann dadurch wird die Anzahl der überhaupt möglichen Kontakte zwischen infizierten und für die Krankheit empfänglichen Individuen modelliert. Diese Anzahl wird (analog zum einfachen Modell mit der Sterberate) sowohl mit der Kontaktrate als auch mit der Schrittweite multipliziert.

Zusammenfassend ist die Anzahl der für die Krankheit empfänglichen Individuen zum Zeitpunkt gleich der Anzahl der für die Krankheit empfänglichen Individuen zum Zeitpunkt abzüglich des Terms, der die Anzahl der Individuen beschreibt, die sich im letzten Zeitschritt mit der Krankheit infiziert haben:

Einfacher nachvollziehen lässt sich, dass

die Anzahl der (infizierten) Individuen ist, die pro Zeitschritt genesen. Entsprechend lässt sich

verstehen und schließlich auch

für die Veränderung der infizierten Individuen.

Anfangsbedingungen

Um eine Simulation unter Verwendung der Rechenvorschrift durchzuführen, müssen wir Anfangsbedingungen definieren. Hierzu verwenden wir den Parameter , der mit einem Wert zwischen 0 und 1 den Anteil der (anfangs) infizierten Individuen beschreibt. Damit definieren wir

als Anfangsbedingungen, die wir in den folgenden Simulationen verwenden werden.

Aufgabe

Mach dich mit der Rechenvorschrift aus diesem Abschnitt vertraut und zeige, dass die Anzahl der Individuen insgesamt zeitlich konstant bleibt. Berechne dazu

und berücksichtige die Anfangsbedingungen.

Quiz
Simulation