Итак, будем считать, что основными действующими факторами инфекционного заболевания являются следующие величины.
1) Концентрация патогенных размножающихся антигенов V(t).
2) Концентрация антител F(t).
3) Концентрация плазматических клеток C(t).
4) Относительная характеристика пораженного органа m(t).
Переходим к построению уравнений модели. Первое уравнение будет описывать изменение числа антигенов в организме:
dV= Vdt-FVdt. (1)
Первый член в левой части этого уравнение описывает прирост антигенов dV за интервал времени dt за счет размножения .Естественно, что он пропорционален V и некоторому числу , которое будем называть коэфициентом размножения антигенов . Член FVdt описывает число антигенов ,нейтрализируемых антителами F за интервал времени dt .В самом деле, число таких вирусов, очевидно,будет пропорционально как количеству антител в организме, так и количеству антигенов; -коэфициент, связанный с вероятностью нейтрализации антигена антителами при встрече с ним. Разделив соотношение (1) на dt получим:
dV/dt=(F)V.
Второе уравнение будет описывать рост плазматических клеток.
dC=F(t-)V(t-V(t-dt-u(C-C*)dt. (2)
Первый член правой части-генерация плазмоклеток,-время,в течение которого осуществляется формирование каскада плазматических клеток,-коэфициент,учитывающий вероятность встречи антиген-антитело, возбуждение каскадной реакции и число образующихся новых клеток.Второй член во второй формуле описывает уменьшение числа плазматических клеток за счет старения, u-коэфициент,равный обратной величине их времени жизни.Разделив соотношение (2) на dt, приходим к уравнению :
dC/dt= F(t-)V(t-V(t- u(C-C*).
Для получения третьего уравнения подсчитывают баланс числа антител, реагирующих с антигеном.Исходят из соотношения:
dF=pCdt-FVdt-ufFdt. (3)
pCdt-генерация антител плазматическими клетками за интевал времени dt, p-скорость производства антител одной плазматической клеткой,FVdt-описывает уменьшение числа антител в интервале времени dt за счет связи с антигенами . ufFdt-уменьшение популяции антител за счет старения,где uf-коэфициент,обратно пропорциональный времени распада антител.Разделив (3) на dt получим:
dF/dt=pC-(uf+ V)F.
Введем в рассмотрение уравнение для относительной характеристики поражения органа- мишени.М-характеристика здорового органа.М*-соответствующая характеристика здоровой части пораженного органа Вводим в рассмотрение величину m по формуле:
m=1-M*/M
Для непораженного органа ,m равна нулю,для полностью пораженного –единице.Для этой характеристики рассмотрим уравнение(четвертое уравнение):
dm/dt=V-um
Первый член правой части характеризует степень поражения органа. V-количество антигенов, где -некоторая константа ,своя для каждого заболевания. Уменьшение этой характеристики происходит за счет восстановительной деятельности организма.
Совершенно ясно, что при сильном поражении жизненно важных органов производительность выработки антител падает. Это является роковым для организма и ведет к летальному исходу. В нашей модели фактор поражения жизненно важных органов можно учесть в уравнении (2), заменив коэффициент на произведение (m). Типичная схема для этой функции представлена на рис.1:
На этом рисунке кривая в интервале 0<=m<=m* равна 1. Это значит, что работоспособность иммунологических органов в этом интервале не зависит от тяжести болезни. Но далее их производительность быстро падает. Таким образом, приходим к следующей системе нелинейных обыкновенных дифференциальных уравнений:
dV/dt=(F)V,
dC/dt= F(t-)V(t-V(t- u(C-C*),
dF/dt=pC-(uf+ V)F,
dm/dt=V-um.
К системе уравнений присоединяют начальные данные при t=t0(V(t0),F(t0),C(t0),m(t0)).Полученную систему уравнений назовем простейшей математической моделью заболевания. Данная математическая модель может использоваться для интерпретации клинических исследований.
Литература :
1.Марчук Г.И. «Математические модели в иммунологии» М.,1994.
2.Амосов Н.М. «Искусственный разум».К.,1986.
3. «Математическое моделирование биологических процессов» М.,1989.
4.Эшби У.Р. «Введение в кибернетику».М.,1990.