Файл: Шахнович, А. Р. Математические методы в исследовании биологических систем регулирования.pdf

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

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

Рис. 34. Схема управления метабо­

ц и я ф е р м е н т а Е а . П р и н а л и ч и и р а в н ы х и с х о д н ы х с у б с т р а т о в д л я с и н т е з а ф е р м е н ­

я н и е м д а ж е с л а б о г о в о з д е й с т в и я в с и с т е м е в о з н и к а е т а с и м м е т р и я , к о т о р а я н е и з б е ж н о

Наиболее тонкий подход приводит к учету клеток различных типов. Тогда рассматривается n -f- 1-камерная система (п — кле­

154

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

где ut — управление.

Однако нелинейное представление достаточно сложно для ана­ литического исследования. Поэтому при составлении математиче­ ской модели метаболического процесса представляется целесообраз­ ным использование метода квазилинеаризации (Buell, Kalaba, 1969) или более общего идентификационного подхода (Bellman, 1969).

Выше, в главе 1-1, была приведена идентификационная про­ цедура, основанная на квазилинеаризации. Эта процедура успеш­ но применяется для аналитического исследования метаболических процессов. Bellman, Jasques, Kalaba (1967) эту процедуру исполь­ зовали для определения в нелинейном дифференциальном урав­ нении кинетики химических процессов

^ = к{а-х){Ь-хУ — Кх* (ІИ-2-1)

скоростных констант к ж К при наличии экспериментальных дан­ ных.

152

с помощью четырех итерации получено значение констант с ошиб­ кой 15%..

Выше приведена линейная модель нормального метаболического процесса без внешних возмущений. Значительный интерес для практики представляет исследование процесса метаболизма при наличии внешнего возмущения — введении препарата в орга* низм.

Пусть г; — скорость введения метаболита; х — общее коли­ чество препарата в крови; у — общее количество препарата в тка­

кх — константа распада.

С учетом изложенного выше общая скорость изменения содержа­ ния препарата в крови складывается из скорости введения и ско­ рости исчезновения, т. е. дифференциальное уравнение имеет вид

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

(где к2 — константа распада в ткани). Таким образом, в линейном представлении динамика описывается уравнениями (ПІ-2-3) и (ИІ-2-4).

Одну из моделей кинетики химических процессов рассмотрели Gonsalez-Fenandes,' Atta (1968). Модель описывается двумя зависи­ мыми переменными Q и £, удовлетворяющими уравнениям в част­ ных производных

где jDr , DZ— значения коэффициентов диффузии в направлениях г и z цилиндрических координат; ср (Q) — функция концентрации кислорода Q для времени і.

153

стоянного возмущения ѵ, что сводится к решению дифференциаль­ ных уравнений.

I I I . Определение оптимального процесса метаболизма при вве­ дении возмущения с постоянной скоростью или скоростью, зави­ сящей от времени ѵ (t) для заранее сформулированного критерия качества процесса.

Требования к качеству при исследовании оптимальной дина­ мики метаболизма могут быть весьма разнообразными:

—оптимальное быстродействие (при исследовании перехода метаболита из одного состояния (количества) — в другое);

154

—оптимальное качество изменения концентрации метаболита;

—оптимальное приближение семейства случайных процессов динамики метаболизма к желаемому детерминированному процессу.

Во всех перечисленных случаях может быть записан функцио­ нал вида

т

J = ^ с(х, и, t) dt.

тересами (синтез и распад, образование и потребление продукта и т. д.), происходящее в сложных .условиях (по температуре, давлению pH и т. д.), при наличии побочных продуктов обмена (Болдуин, 1949).

Среди перечисленных выше задача I является обычной задачей анализа и имеет практическое значение. Для ее решения необхо­ димы конкретные исходные данные. Решение этой задачи (Рашевски, 1966) имеет вид: изменение количества препарата в крови описывается уравнением

а изменение количества препарата в ткани описывается уравне­ нием

Ѵ_Г_і_ ,л _ - K t \ _ 1 - e-( A '-+ g '> ' I

Эти уравнения имеют широкое применение при исследовании ме­ таболических процессов.

При решении задач класса I I I представляется целесообразным использование принципа максимума Л. С. Понтрягина (см. главу 1-2), так как здесь имеют место математическая модель, критерий качества, граничные условия. Оптимальная задача, как известно, формулируется следующим образом: определить управление, принадлежащее к классу допустимых, переводящее точку в фазо­ вом пространстве состояний, динамика которой описывается

так, чтобы критерий качества достигал экстремального значения. Рассмотрим линейную модель, которая характеризует введе­

ние, накопление и удаление лекарственных веществ.

155