{"id":14276,"url":"\/distributions\/14276\/click?bit=1&hash=721b78297d313f451e61a17537482715c74771bae8c8ce438ed30c5ac3bb4196","title":"\u0418\u043d\u0432\u0435\u0441\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0432 \u043b\u044e\u0431\u043e\u0439 \u0442\u043e\u0432\u0430\u0440 \u0438\u043b\u0438 \u0443\u0441\u043b\u0443\u0433\u0443 \u0431\u0435\u0437 \u0431\u0438\u0440\u0436\u0438","buttonText":"","imageUuid":""}

Прогноз пандемии коронавируса. Математическая модель

​Прогноз развития эпидемии в России coronavirussia.online

Как и все сознательные люди, я сижу дома и пытаюсь себя всячески развлечь. Но со всех утюгов только и слышно, что про коронавирус.

Решил я узнать, когда же закончится эта канитель. Покопался в интернете и на то время ничего конкретного ни нашел.

На своей работе я уже давно ничего не программировал, а руки-то чешутся. Вот и придумал написать сайт, который соберет данные из открытых источников, красиво их отобразит их, да и еще сделает прогноз, что будет дальше.

Согласно статистике, развитие эпидемий в мире происходит по нормальному закону распределения Гаусса. Его и взял за основу.

Для начала подключился к открытым источникам и научился выкачивать данные раз в 15 минут. Отобразил их красиво с помощью Chart.js и перешел к прогнозу. Взял для работы функцию нормального распределения

y(t) = b · e -c·(t-a)2

Исходные данные ti - ежедневные отметки времени xi - число активных случаев болезни в этот день

Взял натуральный логарифм от y(t) и от данных xi

ui = ln(xi)

u(t) = ln(y) = ln(b)-c·t2 + 2c·a·t-c·a2

Сделаем замены: А=-с; B=2a·c; C=ln(b)-c·a2; Получаем

u(t) = A·t2+B·t+C

Чтобы найти А, В и С - использую метод наименьших квадратов, согласно которому сумма по i от 1 до n квадратов разностей

Σ(xi-y(ti))2 -> min

должна стремится к минимуму.

Назовем эту сумму Ф(А,В,С)

Продифференцируем dФ по dA, dB и dC. Производные при это должны быть равны нулю.

Получаем систему уравнений

A·Σti4+B·Σti3+C·Σti2 = Σui·ti2

A·Σti3+B·Σti2+C·Σti = Σui·ti

A·Σti2+B·Σti+C·n = Σui

Решаем эту систему по методу Крамера

А = ∆A/∆

B = ∆B/∆

C = ∆C/∆

Находим нужные коэффициенты a,b,c для исходной функции

c = -A a = B/2·c b = eC+c·a2

И, вуаля, готово! Программируем и рисуем красивые графики и узнаем когда это кончится!

Так как данные обновляются с появлением новых случаев, и мой прогноз обновляется тоже.

Чем больше данных, тем точнее прогноз. На данный момент, реальные данные экспертов в Италии говорят о пике эпидемии, и кривая Гаусса показывает тоже самое

​Прогноз развития эпидемии в Италии coronavirussia.online

Будьте здоровы!

0
64 комментария
Написать комментарий...
Lev MU

Профессор из Израиля изобрел способ победить коронавирус за 2 месяца
 Так утверждает профессор Амнон Шаашуа, генеральный директор компании Mobileye и один из ведущих деятелей израильского хай-тека.

31.03.20

Шаашуа вместе с двумя коллегами, разработал математическую модель, с помощью которой можно, по его словам, изолировать группы риска и позволить большинству населения страны постепенно выйти из режима ЧП и вернуться на рабочие места.

Ответить
Развернуть ветку
max krasnoff

Бля, изобретение )))))) Это же просто элементарно и очевидно как дважды два - изолировать только группы риска, а всех остальных не трогать и не заебывать карантинами.

Ответить
Развернуть ветку
Lev MU

Бля, не бля, а к нему уже обратились правительства многих стран, с предложением претворить мат. модель в реальность. Губернатор одного из штатов Америки согласен внедрить эту модель. Это у вас там сидя возле печки всё нипочем, всё то вы знаете , во всем с ходу разбираетесь..хе-хе!! Это тот самый профессор, который внедрил в авто 3-ий  глаз. Такие как вы всёзнальщики тоже чертыхались, эк там изобретение мол невиданное. А это Мобилайн и он продал своё изобретение Интел за 15,3 миллиарда долларов. Вот тебе и эк...

Ответить
Развернуть ветку
61 комментарий
Раскрывать всегда