Я пытаюсь смоделировать зараженную популяцию, используя стандартную модель SIR . Я подгоняю данные с помощью функции curve_fit в Python, описанной здесь , но обнаруживаю, что модель сильно зависит от общей совокупности N. В частности, чем ближе значение N к максимальному числу зараженных людей (в данном случае 40,114), тем лучше подходит. Однако это нереально c, и общая численность населения должна составлять порядка 10 000 000.
Вот данные о ежедневном заражении:
[1.0000e+00 1.0000e+00 1.0000e+00 2.0000e+00 4.0000e+00 7.0000e+00
1.1000e+01 1.7000e+01 2.6000e+01 3.6000e+01 5.1000e+01 1.0700e+02
1.7900e+02 2.6500e+02 4.0700e+02 5.2700e+02 7.7200e+02 1.2080e+03
1.8700e+03 2.6130e+03 4.3490e+03 5.7390e+03 7.4170e+03 9.3090e+03
1.1293e+04 1.3753e+04 1.6378e+04 1.9392e+04 2.2952e+04 2.6312e+04
2.8995e+04 3.1798e+04 3.3765e+04 3.6098e+04 3.7752e+04 3.8926e+04
3.9321e+04 4.0114e+04 4.0101e+04]
Кто-нибудь знает, как решить эту проблему? выдавать? Т.е. как добиться лучшего соответствия с более реалистичной c общей численностью населения N?