Ich habe das SIR Modell mal durchgerechnet für einen Datensatz. Das Vorgehen war dann wie folgt.
Lösen der Dgl. mit den Parametern \( \beta \) und \( \gamma \) die in dem Dgl.System vorkommen und anschliesend die Parameter so variiert, das die Lösungen des Dgl. Systems die gegebenen Daten im quadratischen Mittel am besten approximieren. Ich habe aber nur Daten der Infizierten dazu herangezogen.
Wie Du siehst, geht es einigermassen. Allerdings hängt die Sachen stark von den Startparmetern ab, mit dem die Optimierung beginnt.