ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Публикационная активность

(сведения по итогам 2017 г.)
2-летний импакт-фактор РИНЦ: 0,500
2-летний импакт-фактор РИНЦ без самоцитирования: 0,405
Двухлетний импакт-фактор РИНЦ с учетом цитирования из всех
источников: 0,817
5-летний импакт-фактор РИНЦ: 0,319
5-летний импакт-фактор РИНЦ без самоцитирования: 0,264
Суммарное число цитирований журнала в РИНЦ: 6012
Пятилетний индекс Херфиндаля по цитирующим журналам: 404
Индекс Херфиндаля по организациям авторов: 338
Десятилетний индекс Хирша: 17
Место в общем рейтинге SCIENCE INDEX за 2017 год: 527
Место в рейтинге SCIENCE INDEX за 2017 год по тематике "Автоматика. Вычислительная техника": 16

Больше данных по публикационной активности нашего журнале за 2008-2017 гг. на сайте РИНЦ

Вход


Забыли пароль? / Регистрация

Добавить в закладки

Следующий номер на сайте

4
Ожидается:
16 Декабря 2018

Разработка алгоритма и программы управления процессом эпидемии

Working out of algorithm and the program of control of epidemic process
Статья опубликована в выпуске журнала № 3 за 2009 год.[ 17.09.2009 ]
Аннотация:Рассматривается задача поиска оптимального управления эпидемией с помощью вакцинации, карантина и информационно-просветительской программы «Здоровье» в неоднородном сообществе, состоящем из четырех возрастных групп.
Abstract:The task of search of optimum control by epidemic by means of vaccination, quarantine and the informational-educational program «Health» in the inhomogeneous community consisting of four age groups is considered.
Авторы: Андреева Е.А. (natmat68@mail.ru) - Тверской государственный университет, , , доктор физико-математических наук, Овсянникова Н.И. (natmat68@mail.ru) - Поморский государственный университет, г. Архангельск, ,
Ключевые слова: программное управление, алгоритм численного решения, оптимальное управление эпидемией
Keywords: programmed control, algorithm of numerical solution, optimum control of epidemic
Количество просмотров: 5340
Версия для печати
Выпуск в формате PDF (4.21Мб)

Размер шрифта:       Шрифт:

Существует немало программных продуктов, позволяющих прогнозировать процесс эпидемии, строить динамику распространения инфекции в однородном сообществе, зная начальные данные и интенсивность управления. Данный программный продукт в отличие от предыдущих учитывает три различных вида управления, направленных на погашение инфекции сразу в n различных группах.

Рассмотрим динамику управляемого процесса распространения эпидемии в неоднородном сообществе, состоящем из n возрастных групп, в виде следующей системы дифференциальных уравнений:

(1)

, ,

, , ,                               (2)

где  и  – число подверженных заболеванию и инфицированных в i-й группе (i=1, …, n) в момент t;  – функция, характеризующая скорость заражения подверженных заболеванию из i-й группы от инфицированных из j-й группы с вероятностью  ();  – количество людей, восстановивших за единицу времени свое здоровье в i-й группе без воздействия внешних средств: карантина, вакцинации и др. (– среднее время естественного выздоровления);  – коэффициент естественной смертности людей в i-й группе;  – коэффициент смертности от данной инфекции в i-й группе;  – средний показатель рождаемости в i-й группе;  – доля вакцинированных среди подверженных заболеванию в i-й группе;  – доля больных, отправленных на карантин;  – доля людей, подверженных заболеванию, на которых успешно воздействовали с помощью программы «Здоровье».

Целью управления является минимизация затрат на погашение эпидемии, выраженных функционалом:

,                                          (3)

где A – средняя стоимость одного больного для общества в неделю (известная величина, для России составляет примерно 50 долларов [1]); для одного человека из i-й группы: Di – стоимость вакцинации, Ci – стоимость изоляции, Zi – стоимость программы «Здоровье». Если принять A равным одной условной денежной единице, то (3) перепишется в виде:

,                                               (4)

где di – относительная стоимость вакцинации;  – относительная стоимость программы «Здоровье»; ci – относительная стоимость карантина в i-й группе.

Ограничения на функцию управления заданы в следующем виде:

,     (5)

где Аi, Вi, Ci – максимальные нормы управлений в i-й группе, ограниченные техническими и материальными возможностями.

Разобьем отрезок интегрирования [0; T] на q равных частей точками  так, что  . Задача (1)–(5) аппроксимируется следующей дискретной задачей оптимального управления:

,       (6)

 (7)

где  .

                                                       (8)

, ,

.                           (9)

Подпись:  
Рис. 1. Диалоговая форма для ввода данныхСоставим функцию Лагранжа для дискретной задачи, откуда найдем производные лагранжиана по управлениям , ,  и по фазовым переменным, чтобы выразить сопряженные функции:

 (10)

Условия трансверсальности:

.                                      (11)

Алгоритм численного решения

1. Зададим начальное управление

Подпись: Чел. А БРис. 2. Динамика подверженных заболеванию (А)и управление вакцинацией (Б) в 3-й группев зависимости от стоимости вакцинации 2. Зададим начальные значения  . По формуле (7) вычислим допустимые траектории, по формуле (6) – начальное значение интеграла.

3. Для  вычислим по формуле (10) с учетом (11) ,,.

4. Вычислим следующее приближение для управления:

.

Проверим условие (9): если оно не выполняется, делаем проекцию градиента:

если , то , если , то  ;

если , то , если , то  ;

если , то , если , то  .

5. По формуле (7) вычислим допустимые траектории , по формуле (6) – следующее значение интеграла .

6. Сравним значения  и . Если , переходим к пункту 3, присвоив  найденные в пункте 4 значения управлений. Процесс продолжается до тех пор, пока .

7. Если , возвращаемся к пункту 4 и уменьшаем шаг , например в два раза, проверяем условие (9).

Для проведения эксперимента на реальной модели были взяты статистические данные по развитию эпидемии гриппа в г. Архангельске, предоставленные Территориальным управлением по эпиднадзору. Рассмотрим четыре возрастные группы: 1-я группа – 0–2 года, 2-я группа – 3–6 лет, 3-я группа – 7–15 лет, 4-я группа – старше 15 лет. Коэффициенты смертности, средний показатель рождаемости для каждой группы вычислены по статистическим данным для г. Архангельска. Коэффициенты β найдены путем решения обратной задачи.

Программа по нахождению оптимального управления написана в среде Delphi7 на языке Object Pascal. Все параметры модели пользователь может вводить с клавиатуры, что позволяет исследовать решение задачи в зависимости от введенных параметров.

Подпись:  А БРис. 3. Динамика инфицированных (А) и управлениеизоляцией (Б) в 3-й группе в зависимостиот стоимости изоляцииПосле окончания работы программы искомое значение функции выводится в окно Результат (рис. 1), а графическое решение (фазовые переменные, программное управление, функции переключения управления) по каждой группе можно получить в окне Графики нажатием соответствующей кнопки.

Кроме того, решение дискретной задачи оптимального управления выводится в Excel, что позволяет сравнивать решения задачи при различных параметрах. Наиболее сильное влияние на решение оказывают стоимость управлений, величина временного интервала [0; T] и коэффициенты β. Зафиксировав стоимость изоляции и программы «Здоровье», выявим влияние стоимости вакцинации на решение задачи. Так как изменения во всех группах аналогичные, то для примера выясним, как меняются динамика подверженных заболеванию и управление вакцинацией в 3-й группе (рис. 2), которые сильнее всего зависят от стоимости вакцинации (d1>d2>d3).

Из рисунка видно, что с уменьшением стоимости вакцинации уменьшается остаточное на

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

Аналогично можно рассмотреть влияние стоимости изоляции (карантина) на динамику инфицированных и управление изоляцией в 3-й группе (рис. 3), которые сильнее всего зависят от стоимости изоляции (с1>с2>с3).

Программа позволяет исключить отсутствующий вид управления, чтобы проследить эффективность оставшихся видов управления и суммарные затраты на погашение эпидемии. Таким образом, было выяснено, что самым эффективным и дешевым способом управления является комплексное управления вакцинацией, карантином и программой «Здоровье».

Литература

1.   Вакцинация. Новости вакцинопрофилактики // Информационный бюллетень. М. 2003. № 3 (27).

2.   Андреева Е.А. Оптимальное управление динамическими системами. Тверь: Изд-во ТвГУ, 1999. С. 72–120.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=2342
Версия для печати
Выпуск в формате PDF (4.21Мб)
Статья опубликована в выпуске журнала № 3 за 2009 год.

Назад, к списку статей

Хотите оценить статью или опубликовать комментарий к ней - зарегистрируйтесь