
- •3. Математичне моделювання в імунології
- •3.1. Основні компоненти імунної відповіді
- •3.2. Форми інфекційного захворювання
- •3.3. Проста модель інфекційного захворювання
- •3.4. Існування та невід’ємність розв’язку
- •3.5. Стаціонарні розв’язки
- •Розв’язавши систему, одержимо
- •3.6. Стійкість стану здорового організму
- •3.9. Вплив температурної реакції організму на динаміку захворювання
- •Список літератури
- •3.1. Основні компоненти імунної відповіді............................................. 57
3.4. Існування та невід’ємність розв’язку
Розглянемо модель імунної відповіді
,
, (3.8)
,
,
де t 0, 0 m 1, функція (m) визначена згідно (3.6).
Теорема 3.1. [2, c. 55-56] Нехай параметри моделі (3.8) і початкові умови (3.7) – невід’ємні числа.
Тоді для всіх t > 0 існує диференційовний, за винятком окремих точок, розв’язок задачі (3.8), (3.7).
Зауважимо, що внаслідок розриву початкової функції V(t) при t = 0, компонента розв’язку V(t) тільки неперервна для t = , але з ростом t її гладкість, як і всього розв’язку системи (3.8), зростає.
Теорема 3.2. Якщо виконуються умови теореми 3.1, то розв’язок початкової задачі (3.8), (3.7) невід’ємний для всіх t 0.
Доведення.
Оскільки
,
то
і при
для деякого
.
Із
рівняння (3.1) маємо:
,
Звідки
,
(3.9)
для всіх
,
якщо
V0
0.
Зокрема, V(t)
0, якщо V0
= 0 і V(t)
> 0, якщо V0
> 0.
Розглянемо тепер (3.4) як лінійне неоднорідне рівняння із неод-норідністю V(t). Розв’язок такого рівняння записується у вигляді [9]:
. (3.10)
Оскільки m0 0, 0 і V(t) 0 для t 0, то m(t) 0 для t 0. Зокрема, m(t) > 0, якщо m0 > 0 або V0 > 0.
Для t [0, ), враховуючи умови (3.7), рівняння (3.2) запишемо у вигляді
.
Розв’язок цього лінійного рівняння
C(t)
=
,
оскільки
.
Із
неперервності розв’язку F(t)
випливає, що F(t)
> 0 на деякому інтервалі
,
якщо F0
> 0. Якщо на
існує точка t1,
в якій F(t1)
= 0 і t1
–
перша
така точка. Тоді F
(t1)
0. Але на підставі (3.3)
що суперечить припущенню і тому F(t) 0 для t [0, ].
Якщо F 0 = 0, то F (t) > 0 для t > 0. Тому F(t) > 0 для t (0, t2). Далі доведення аналогічне.
Нехай t [, 2). Оскільки (m(t)) 0, V(t) 0 для t 0 і F(t – ) 0, t [, 2), то
= (m(t))F(t
– )V(t
– )
– c(C
– C*)
–c(C
– C*).
Звідси для t [, 2] одержимо
C(t)
0. (3.11)
Тепер можна показати, що F(t) 0 для [, 2]. Ці міркування можна продовжити для t [2, 3], [3, 4], ..., що завершує доведення теореми.
Наслідок 1. Якщо V0 > 0, то V(t) > 0 і m(t) > 0 для всіх t 0.
Нерівності випливають з (3.9) і (3.10).
Наслідок 2. Якщо в умовах теореми 3.2 C0 C*, то C(t) C* для всіх t 0.
Справді,
C(t)
= (C0
– C*)
+ C
*
C
*, t
[0, ].
Для t > доведення одержується із нерівності (3.11).
Біологічний зміст одержаної нерівності полягає в тому, що якщо концентрація плазмоклітин в початковий момент перевищує C *, то такий стан імунної системи збережеться і надалі.
Наслідок 3. Якщо C0 = C*, то C(t) = C* для t [0, ] і C(t) > C* для t > .
Наслідок
4.
Якщо V0
> 0, 0 <
C0
< C*,
то
таке, що 0
< C(t)
<
C*,
t
(0, t2);
C(t2)
= C*
і C(t)
> C*
для t
> t2.
Справді,
(t)
= –c(C0
– C*)
> 0,
тому функція (3.10) зростає на (0,
1),
0 <
< .
Оскільки C(1)
= C
* + (C0
– C*)
< C*,
то t2
> .
Якби
t2
= +,
то
= 0,
причому C(t)
< C*.
Звідси
випливає, що
(t)
+0, F(t)
+0,
C(t).
Оскільки С()
> C0,
то T
>
таке, що
t
T
> 0.5C0
> 0. Тому
F(T)
> F()
+ 0.5C0(T
– ),
що суперечить F(t)
+0,
коли
t
.
Нехай t3 > t2 – перша точка, в якій C(t3) = C*, C(t3) 0. Тоді
> 0,
що суперечить припущенню.