Научная статья на тему 'Класс кинетически моделируемых биомеханических случайных процессов'

Класс кинетически моделируемых биомеханических случайных процессов Текст научной статьи по специальности «Математика»

CC BY
111
31
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
НЕСТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ / БИОМЕХАНИЧЕСКИЕ ПРОЦЕССЫ / КИНЕТИЧЕСКОЕ МОДЕЛИРОВАНИЕ / ФУНКЦИЯ СИСТЕМЫ ДВУХ СЛУЧАЙНЫХ ВЕЛИЧИН / ПРОЗРАЧНЫЙ СТАТИСТИЧЕСКИЙ АНАЛИЗ / БИОМЕХАНИЧЕСКОЕ СОПРОВОЖДЕНИЕ РЕАЛИЗАЦИЙ / NON-STATIONARY RANDOM PROCESSES / BIOMECHANICAL PROCESSES / KINETIC SIMULATION / THE FUNCTION OF TWO RANDOM VARIABLES / CLEAR STATISTICAL ANALYSIS / BIOMECHANICAL ACCOMPANIMENT OF REALISATIONS

Аннотация научной статьи по математике, автор научной работы — Селянинов А. А.

Биомеханические процессы, связанные с биологическими системами и человеческим организмом, в силу множества факторов, определяющих их развитие со временем, являются, по сути, случайными. Просматривается класс кинетически моделируемых биомеханических случайных процессов. Особенностью класса является определенная монотонность при возрастании или убывании определяющего процесс параметра и наличие современных методов его регистрации, что позволяет произвести кинетическое моделирование реализаций случайного процесса. В результате нестационарный случайный процесс сводится к неслучайной функции системы двух случайных величин вместо системы бесконечного числа случайных величин в известном методе канонического разложения В.С. Пугачева. Для процессов, моделируемых кинетическими уравнениями первого порядка с линейно изменяющейся «константой» реакции, получены аналитические зависимости математического ожидания и дисперсии от времени. При выборе скорости биодеструкции в качестве определяющего параметра для нового восприятия процесса также оказалось возможным определение числовых характеристик в виде аналитических зависимостей от времени. Для углубленного вероятностного анализа таких нестационарных случайных процессов в работе предлагается ввести гипотезу о нормальном законе распределения определяющего процесс параметра для каждого фиксированного момента времени, для чего необходимо воспринимать текущее время в качестве параметра. Конечно, желательно подкрепить эту гипотезу хотя бы информацией об одинаковых пределах уменьшения–увеличения активности актинобактерий на момент начала процесса биодеструкции лекарственного средства. Предложенный подход позволяет значительно упростить и сделать прозрачным статистический анализ таких процессов in vitro, а также получить новый метод биомеханического сопровождения конкретных реализаций in vivo.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Селянинов А. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The biomechanical processes associated with human body and biological systems have a large body of factors determining their development in the course of time; therefore these processes are essentially random. We suggest discerning a class of biomechanical random processes that allow kinetic simulation. The main features of this class are the certain monotony under increasing or decreasing the characteristic parameter of the process, and the availability of modern methods of its registration. That allows performing a kinetic simulation of the random process realisations. As a result, the non-stationary random process is reduced to a non-random function of two random variables instead of an infinite number of random variables in the V.S. Pugachev canonical decomposition method. The time dependences of the mathematical expectation and variance were derived analytically for the processes that may be simulated by the first order kinetic equations with a linearly varying reaction factor. The determination of the numerical characteristics such as analytical functions of time has also been possible at change of biodegradation rate for the new perception of the process. We propose to introduce the hypothesis of normal distribution of determining parameter of the process for each fixed moment of time for the advanced analysis of the probable unsteady random processes; the current time is taken a parameter. Of course, this hypothesis requires at least of information about the same ranges of decreasing-increasing of the activity of actinobacteria at the start of the biodegradation process of the drug. The proposed approach allows simplifying and making clear the statistical analysis of these processes in vitro, and also getting an novel method for biomechanical accompaniment of the specific realisations in vivo.

Текст научной работы на тему «Класс кинетически моделируемых биомеханических случайных процессов»

УДК 531/534:[57+61]

КЛАСС КИНЕТИЧЕСКИ МОДЕЛИРУЕМЫХ БИОМЕХАНИЧЕСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ

А.А. Селянинов

Кафедра теоретической механики Пермского национального исследовательского политехнического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: luba@theormech.pstu.ac.ru

Аннотация. Биомеханические процессы, связанные с биологическими системами и человеческим организмом, в силу множества факторов, определяющих их развитие со временем, являются, по сути, случайными. Просматривается класс кинетически моделируемых биомеханических случайных процессов. Особенностью класса является определенная монотонность при возрастании или убывании определяющего процесс параметра и наличие современных методов его регистрации, что позволяет произвести кинетическое моделирование реализаций случайного процесса. В результате нестационарный случайный процесс сводится к неслучайной функции системы двух случайных величин вместо системы бесконечного числа случайных величин в известном методе канонического разложения В.С. Пугачева. Для процессов, моделируемых кинетическими уравнениями первого порядка с линейно изменяющейся «константой» реакции, получены аналитические зависимости математического ожидания и дисперсии от времени. При выборе скорости биодеструкции в качестве определяющего параметра для нового восприятия процесса также оказалось возможным определение числовых характеристик в виде аналитических зависимостей от времени. Для углубленного вероятностного анализа таких нестационарных случайных процессов в работе предлагается ввести гипотезу о нормальном законе распределения определяющего процесс параметра для каждого фиксированного момента времени, для чего необходимо воспринимать текущее время в качестве параметра. Конечно, желательно подкрепить эту гипотезу хотя бы информацией об одинаковых пределах уменьшения-увеличения активности актинобактерий на момент начала процесса биодеструкции лекарственного средства. Предложенный подход позволяет значительно упростить и сделать прозрачным статистический анализ таких процессов in vitro, а также получить новый метод биомеханического сопровождения конкретных реализаций in vivo.

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

Введение

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

В медицинской практике до сих пор основой анализа остается статистический подход к схеме: диагноз - вариант лечения - результат. Варианты лечения

© Селянинов А.А., 2012

Селянинов Александр Анатольевич, д.т.н., профессор кафедры теоретической механики, Пермь

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

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

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

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

В связи с этим предлагаются подходы к улучшению анализа определенного класса биомеханических случайных процессов:

1. Сведение анализа нестационарного случайного процесса к анализу неслучайной функции системы ограниченного числа случайных величин с применением кинетического моделирования (ожидаемый результат - существенное улучшение статистической обработки и ее прозрачность in vitro).

2. Биомеханическое сопровождение конкретной реализации случайного процесса с применением параллельного процессу кинетического моделирования при возможности мониторинга определяющего процесс параметра (ожидаемый результат -индивидуальный подход к реализации случайного процесса in vivo).

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

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

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

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

Процессы с наличием биомеханической составляющей

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

В качестве примера можно выделить ряд биомеханических процессов, которые входят в интересующий класс:

- утилизация лекарственных средств с использованием микроорганизмов (перемешивание культуральной жидкости позволяет организовать доставку питательной среды колониям микроорганизмов механическим путем) [3-5, 6, 14];

- шинирование зубного ряда (развитие пародонтита зубов нижней челюсти человека приводит к прогрессирующему расшатыванию зубов при физиологических нагрузках, операция шинирования замедляет развитие пародонтита) [7];

- устранение дефектов зубного ряда (протрузия - ретрузия передних резцов верхней челюсти человека исправляется со временем силовым давлением ортодонтических аппаратов, в частности эластопозиционерами) [18-20];

- биодоступность содержимого капсул с лекарственными средствами и биологически активными добавками (растворимость желатиновой оболочки зависит от ее влажности, величина которой в процессе сушки определяется влагопереносом от оболочки к наполнителю, а затем от наполнителя к оболочке и в окружающую среду) [2, 11];

- операция стапедопластики (скорость улучшения слуха после операции в зависимости от степени отосклероза) [13, 15];

- устранение дисплазии тазобедренного сустава (создание условий для формирования нового сустава после оперативного вмешательства посредством управления скоростью сращивания разрезанных тазобедренных костей) [16];

- послеоперационное заживление мягких тканей человека (движение инфильтрата из внутренних слоев мягких тканей не должно ограничиваться опережающим заживлением наружных слоев) [1];

- устранение дефектов зубного ряда с применением имплантатов (создание условий для снижения скорости естественной убыли кости альвеолярного отростка путем оперативного вмешательства) [10].

Особенности данного класса случайных процессов

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

Эти особенности заключаются в возможностях:

- выбора определяющего процесс параметра, для которого возможен постоянный во времени или дискретный мониторинг доступными современными экспериментальными методами анализа;

- исследования кинетики процесса и проведения параллельного процессу кинетического моделирования с некоторой аналогией с кинетикой химических процессов;

- анализа нестационарного случайного процесса как обычной функции системы ограниченного числа случайных величин с применением кинетического моделирования;

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

Уровень связанности таких постановок задач биомеханики может быть различным:

- односторонним (кинетическое уравнение как аппроксимация

экспериментальной кривой определяет граничные условия в краевой задаче или результаты решения краевой задачи влияют на скорость изменения определяющего процесс параметра при биомеханическом сопровождении случайного процесса);

- взаимным (кинетическое уравнение определяет изменение параметров процесса, а биомеханика влияет на параметры кинетического уравнения).

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

Нестационарные случайные биомеханические процессы

Рассмотрим эти процессы на примере биодеструкции дротаверина гидрохлорида клетками Я. Rhodochrous ИЭГМ 647 [12].

Обозначим через Х(^) случайный процесс изменения со временем определяющего параметра (концентрации в культуральной жидкости дротаверина гидрохлорида), для постоянного (или хотя бы дискретного) мониторинга которого имеется современное оборудование. Через х(^) обозначим реализацию случайного процесса изменения концентрации, случайность которого складывается из двух аспектов.

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

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

Второй аспект связан с погрешностью измерений современными методами начальной концентрации субстрата х0 и ее текущего значения х, не говоря уже о сложности идентификации продуктов разложения и их количественной оценки. Относительная погрешность современного метода определения концентраций реагентов с применением высокоэффективной жидкостной хроматографии находится в пределах 2%. Она ничтожно мала в сравнении с разбросом данных по повторяемости реализаций, поэтому, действительно, случайность результата биодеструкции заложена именно в биологической и биомеханической сущности процесса. При этом для каждой пробы производится 3-5 измерений и представляется усредненное значение концентрации субстрата.

В результате при экспериментальном изучении случайного процесса Х(^), исходя из требований фармакопеи [8], имеем при каждом эксперименте на повторяемость три серии значений начальной х0 и текущей х концентраций субстрата через определенные промежутки времени (табл. 1). Данные табл. 1 соответствуют одному из вариантов при поиске рационального режима биодеградации дротаверина гидрохлорида.

Таблица 1

Динамика содержания дротаверина гидрохлорида в процессе биодеградации клетками R. Rhodochrous ИЭГМ 647

t, сутки 0 5 20 30 48

х, мг/мл 0,002 0,001818 0,000712 0,000494 0,000020

0,002 0,001820 - 0,000266 0,000020

0,002 - - 0,000236 0,000230

.V, мг/мл

0,0016

0,0012

0,0008

0,0004

0 10 20 30 40 t, сутки

Рис. 1. Реализации случайного процесса биодеструкции при экспериментах

на повторяемость: n = 1, 2, 3 - номер эксперимента

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

По данным табл. 1 построены кусочно-линейные графики реализаций случайного процесса биодеструкции дротаверина гидрохлорида при экспериментах на повторяемость в виде зависимостей концентрации х от времени (рис. 1). Число реализаций связано с требованиями фармакопеи [8].

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

В предыдущих работах по кинетическому моделированию процессов биодеструкции лекарственных средств использованы кинетические уравнения первого порядка, адекватно отражающие динамику изменения начальных концентраций как исходных субстратов (парацетамола, и-аминофенола в частности), так и первичных продуктов их разложения (и-аминофенола при биодеградации парацетамола например) [3-5, 6, 14].

26 ISSN 1812-5123. Российский журнал биомеханики. 2012. Т. 16, № 4 (58): 22-35

.V, мг/мл

0,0016

0,0012

0,0008

0,0004

0 10 20 30 40 Г, сутки

Рис. 2. Реализации случайного процесса биодеструкции при экспериментах

на повторяемость после кинетического моделирования: п = 1, 2, 3 - номер эксперимента;

*

х - усредненная кривая

При биодеструкции в качестве первичного субстрата (парацетамол,

и-аминофенол, дротаверина гидрохлорид) независимо от штамма родококков,

участвующих в реакциях, процесс адекватно отражают кинетические уравнения первого порядка типа

dx

— = -к 0) • х (1)

dt

где «константы» реакции к{() = Ш + Ь - линейные функции времени биодеградации.

Здесь параметр Ь имеет размерность [1/время] и отвечает за скорость процесса биодеструкции, а параметр а имеет размерность [1/время ] и отражает величину ускорения или замедления процесса во времени. При этом скорость процесса биодеструкции у(^ равна правой части в уравнении (1) и зависит от концентрации утилизируемого субстрата.

Результат интегрирования кинетического уравнения (1) имеет экспоненциальный вид

х = х0 • ехр(-(Ь + at /2) t), (2)

где х0 - начальная концентрация утилизируемого лекарственного средства.

На рис. 2 представлены результаты кинетического моделирования реализаций исследуемого случайного процесса биодеструкции.

Значения параметров а и Ь для уравнений (1, 2) получены по данным, приведенным в табл. 1, с применением метода наименьших квадратов. Как средняя арифметическая получена усредненная кривая х = (х1 + х2 + х3)/3, которая, в принципе, является математическим ожиданием случайного процесса биодеструкции. Конечно, необходимо количество реализаций числом значительно большим, чем три. Однако стоимость насыщения экспериментальными данными табл. 1 высока, поэтому приходится ограничиваться реальным минимумом из трех реализаций случайного процесса.

В результате для данного случайного процесса возможно определение математического ожидания mx (t), но для полной информации о процессе необходима

корреляционная функция Kx (t, t') , в частности дисперсия Dx (t), определение которых

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

Можно попытаться ввести некоторые гипотезы:

- об эргодичности случайного процесса, однако считать процесс эргодическим нельзя, так как с каждой реализацией случайного процесса изменяется математическое ожидание, а главное - изменяется дисперсия как функция времени;

- о стационарности случайного процесса, однако математическое ожидание mx (t) ^ const, хотя это не препятствует изучению случайного процесса как

стационарного при центрировании случайной функции X (t) = X (t) - mx (t), но считать

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

дисперсию Dx (t) = Dx = const можно в качестве верхней оценки со значительным ее

увеличением (особенно это неправомерно в начале процесса, когда значения определяющего параметра практически детерминированные);

- о сокращении степени случайности случайного процесса - исключить погрешность измерений концентрации первичного субстрата, оставив первый аспект о чисто биологической и биомеханической сущности процесса;

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

Исходя из этих соображений, случайный процесс явно нестационарный, и даже определение дисперсии случайного процесса Dx (t) является проблемой.

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

Случайный процесс как функция системы двух случайных величин

Как уже отмечалось, нестационарные случайные процессы чрезвычайно сложны для анализа в полном объеме. Одним из способов анализа является метод канонического разложения В.С. Пугачева [17]. Он позволяет свести нестационарный случайный процесс к функции системы бесконечного числа независимых случайных величин с определением корреляционной функции и дисперсии, зависящих от времени.

В настоящее время просматривается определенный класс случайных процессов, реализации которых подлежат кинетическому моделированию, традиционно применяемому в кинетике химических процессов [9]. В кинетической химии скорость убыли реагента описывается константой скорости реакции и порядком кинетического уравнения. Биомеханические процессы, в общем смысле биологические процессы, связаны с живой материей, которой свойственна адаптация к изменившимся условиям жизнедеятельности. Поэтому наряду с порядком кинетического уравнения (первым, нулевым, минус первым) здесь появляется ускорение (замедление) так называемой константы скорости реакции. Ряд работ позволяет считать это ускорение (замедление) константы скорости реакции постоянным в каждой реализации случайного процесса [3-5, 6, 14].

В процессах биодеструкции монотонно убывающую (вернее, невозрастающую) функцию для каждой реализации x(t), которая напоминает спадающую экспоненту (первый порядок кинетического уравнения), можно описать двумя параметрами:

а - ускорение процесса, Ь - скорость процесса в выражении так называемой константы скорости реакции из кинетики химических процессов. При постоянной величине ускорения «константа» скорости реализации случайного процесса примет вид линейной зависимости от времени к^) = at + Ь в уравнении (1), где а и Ь - неслучайные величины.

Для случайного процесса Х(^, включающего все реализации, ускорение а и скорость процесса Ь приобретают смысл случайных величин. Тогда нестационарный случайный процесс данного класса можно свести к функции системы двух случайных величин вместо функции (в виде ряда) системы бесконечного числа пар некоррелированных случайных величин с нулевыми математическими ожиданиями и одинаковыми для пар дисперсиями согласно методу канонического разложения В. С. Пугачева [17].

Следует подчеркнуть, что нестационарный случайный процесс сводится к элементарной функции системы двух случайных величин, а не к функции одной пары некоррелированных случайных величин с нулевыми математическими ожиданиями и одинаковыми для пары дисперсиями. Для случая коррелированных случайных величин а и Ь в выражение для дисперсии добавляется слагаемое с корреляционным моментом КаЬ.

В данном случайном процессе ускорение а и скорость Ь - случайные величины с математическими ожиданиями та, тЬ и соответствующими дисперсиями Ба,

По известным выражениям [17] для а и Ь в качестве случайных величин в кинетических уравнениях типа (1) при трех реализациях были определены математические ожидания та и тЬ, дисперсии Ба и БЬ, а также член корреляционной матрицы КаЬ:

математическим ожиданием понимается среднее выборки, под дисперсией -выборочная дисперсия, под средним квадратичным отклонением - стандартное отклонение.

Равенство нулю члена корреляционной матрицы КаЬ отражает некоррелированность случайных величин а и Ь, хотя в отличие от метода

В.С. Пугачева этот факт не имеет значения в предлагаемой методике определения математического ожидания тх ^) и дисперсии Бх ^) рассматриваемого класса

нестационарных случайных процессов.

Если в зависимости реализаций (2) под а и Ь понимать случайные величины, тогда можно записать аналогичное выражение для случайного процесса в виде

Зависимость Х(?) от системы двух случайных величин а и Ь экспоненциальная, поэтому линеаризуем ее в достаточно малой окрестности математических ожиданий та и тЬ , тогда математическое ожидание случайного процесса определится той же самой функциональной зависимостью [17] (время t является параметром и для любого О

Отметим тот факт, что графики зависимостей математического ожидания от времени, полученные ранее усреднением реализаций и построенные по выражению (4), достаточно близки друг к другу (рис. 2). Погрешность заложена в процедуре линеаризации.

Корреляционная функция случайного процесса Кх ^, t') в рамках данной работы мало интересует, остается вопрос о характеристике разброса реализаций случайного процесса относительно математического ожидания тх ^), т.е. о дисперсии Бх ^) .

(3)

(4)

Рис. 3. Дисперсия случайного процесса биодеструкции дротаверина гидрохлорида

при экспериментах на повторяемость

Дисперсия линеаризованной функции системы двух некоррелированных случайных величин определяется выражением [17]

В=2

2=1

дх.

V 2 ^ ша,шь

■ 22

С дх ^

дх,

дх,.

■ К, і = 1, 2; х1

а, х2

Ь, (5)

а для некоррелированных случайных величин согласно условию КаЪ = 0 выражение для дисперсии принимает вид

В =2

і=1

Удхі У

■ Ох , і = 1, 2; х1 = а, х2 = Ь .

(6)

,,тЬ

откуда определяется среднее квадратичное отклонение ах(г).

Для рассматриваемого моделирования нестационарного случайного процесса с применением кинетических уравнений первого порядка зависимость дисперсии от времени принимает вид

вх (г) = тх2(0 • г2 • (Въ + г2 • ва/4), (7)

откуда следует выражение для среднего квадратичного отклонения

Л гл / у|\0,5

°х (і) = тх (і) ■ * ■ (ВЬ + І ■ Ва/4)

(8)

На рис. 3 приведена зависимость дисперсии случайного процесса от времени.

Дисперсия определяет разброс значений реализаций случайного процесса, кривая на рис. 3 имеет максимум при г = 18 суток, что соответствует наибольшему разбросу реализаций на этот момент времени (рис. 2). Однако из-за несоответствия единиц измерения математического ожидания и дисперсии их сравнительный анализ затруднен, поэтому в качестве характеристики разброса реализаций лучше использовать среднее квадратичное отклонение а х (г).

Лл (0. мг/мл

0,0008

0,0006

0,0004

0,0002

0 10 20 30 40 сутки

Рис. 4. Зависимость среднего квадратичного отклонения случайного процесса

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

На рис. 4 приведена зависимость среднего квадратичного отклонения случайного процесса от времени. Как и следовало ожидать, максимальное значение среднее квадратичное отклонение принимает в тот же момент времени г = 18 суток, что и дисперсия.

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

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

В результате при фиксированном времени получаем сечение случайного процесса с количеством случайных величин, равным числу реализаций. Целью этого подхода является применение известного «правила трех сигм» [17].

Тогда с вероятностью 68,2% можно определить достаточно узкую полосу реализаций случайного процесса, с вероятностью 95,4% - расширенную полосу и с вероятностью 97,8% - практически возможную полосу разброса реализаций случайного процесса биодеструкции дротаверина гидрохлорида.

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

В качестве развития идеи о применении кинетического моделирования при исследовании данного случайного процесса следует отметить, что изменение скорости нестационарного процесса X(t) также является случайным процессом V(t). Поэтому помимо математического ожидания mx (t) и дисперсии Dx (t) нестационарного

случайного процесса интересует математическое ожидание mv (t) и дисперсия Dv (t)

скорости V(t) этого процесса, реализации которой определяются выражением (1)

dx

v(t) = —— = -(at + b) • x . (9)

dt

Повторив процедуру линеаризации, по выражениям, аналогичным (5) и (6), получим математическое ожидание

mv (t) = -(mb + mat) • mx (t), (10)

дисперсию

Dv (t ) = ml (t ) • {[t (mb + mat ) - 1]2 Da + 12 f0,5 t (mb + mat ) - 1] 2 Da } (1 1 )

и среднее квадратичное отклонение для скорости процесса биодеструкции

°v(t) = mx(t)-{[t(mb + mat)- 1]2 Da + t2 [0,5 t(b + mat) - ^ Da } • (12)

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

В табл. 2 приведены основные признаки такого класса и отмечены возможности решения ряда связанных с кинетическим моделированием задач.

Отметим нюансы анализа случайных процессов in vitro и in vivo.

Определение для случайного процесса биодеструкции математических ожиданий ma, mb и дисперсий Da, Db важно при наборе статистики по отработанным методически процессам (in vitro), т.е. после их окончания.

При индивидуальном прогнозировании конкретной реализации x(t) или v(t) случайного процесса биодеструкции (in vivo) с применением кинетического моделирования можно отойти от терминологии случайных процессов и случайных величин исходя из идеи параллельного процессу кинетического моделирования [14]. Ее смысл в том, что при наличии трех первых экспериментальных точек можно построить первое приближение кинетической кривой реализации случайного процесса, которое уточняется с появлением каждой новой экспериментальной точки. Причем первая из этих трех точек - начальное значение определяющего процесс параметра. Следует заметить, что при параллельном процессу кинетическом моделировании в кинетическом уравнении (1) все величины являются детерминированными, что принципиально упрощает кинетический анализ.

Следует также отметить, что для прогнозирования (in vivo) окончания процесса по предельной концентрации необходима оценка математического ожидания, дисперсии и среднего квадратичного отклонения на момент времени t = tnp, где время

окончания процесса tnp определяется через предельную концентрацию хПР полученным из уравнения (2) выражением

1

tn, =-

a

f i-------------------------Л

-b+ b2 -2a • ln ^

4 Xo ,

(13)

Таблица 2

Характеристика класса биомеханических кинетически моделируемых случайных процессов

Основные признаки класса Класс биомеханических кинетически моделируемых случайных процессов

Наличие биомеханической составляющей в случайном процессе X (ґ)

Наличие количественной характеристики определяющего параметра случайного процесса в каждой реализации х(ґ)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Измеримость современными методами диагностики определяющего параметра в реализациях х(ґ)

Нестрого возрастающая или убывающая зависимость реализаций х(ґ) от времени, достаточно гладкая для линеаризации

Достаточно малое число реализаций и количество экспериментальных точек в каждой реализации

Возможность идентификации порядка п кинетического уравнения для Жх реализации — = - к(ґ)хп, п = 1, 0, -1 Ж

Наличие статистических данных по X(t) (in vitro) Возможность периодического мониторинга x(t) (in vivo)

Основные возможности анализа Пополнение статистики для «константы реакции» к (t) -по окончании процесса Параллельное процессу кинетическое моделирование с уточнением «константы реакции» к (t) -по начальным точкам

Стохастический прогноз течения случайного процесса X(t) : mx (t), сx (t) Индивидуальный прогноз течения реализаций случайного процесса x(t)

Постановка и решение связанных с кинетическим уравнением краевых задач биомеханики

Биомеханическое сопровождение врачебных или иных вмешательств в естественное течение биомеханического процесса х(^)

Выводы

1. Биомеханические процессы, связанные с биологическими системами и человеческим организмом, в силу множества факторов, определяющих их развитие со временем, являются по сути случайными. В работе сформирован класс кинетически моделируемых биомеханических случайных процессов. Его особенностью является определенная монотонность при возрастании или убыли определяющего процесс параметра и наличие современных методов его регистрации, что позволяет произвести кинетическое моделирование реализаций случайного процесса. Отмечен ряд биомеханических процессов, причисляемых к данному классу.

2. Предложен метод сведения стохастического анализа нестационарного случайного процесса из данного класса к анализу поведения неслучайной функции системы двух случайных величин вместо системы бесконечного числа случайных величин в известном методе канонического разложения В.С. Пугачева.

3. Для нестационарных случайных процессов, моделируемых кинетическими уравнениями первого порядка с линейно изменяющейся «константой реакции», получены аналитические зависимости математического ожидания, дисперсии и среднего квадратичного отклонения от времени. Предложенный подход позволяет значительно упростить и сделать прозрачным статистический анализ таких процессов in vitro, а также получить интересный метод биомеханического сопровождения конкретных реализаций in vivo.

Благодарности

Автор выражает благодарность за помощь в обсуждении результатов исследований заведующему кафедрой теоретической механики Пермского национального исследовательского политехнического университета профессору Ю.И. Няшину, директору Института экологии и генетики микроорганизмов член-корреспонденту Российской академии наук профессору И.Б. Ившиной, заведующей кафедрой аналитической химии Пермской государственной фармацевтической академии доктору фармацевтических наук Е.В. Вихаревой.

Список литературы

1. Баранова А.А., Селянинов А.А., Вихарева Е.В. Кинетическое моделирование биомеханических процессов // Вестник Пермского национального исследовательского политехнического университета. Механика.. - 2012. - № 3. - С. 7-25.

2. Бокова С.А., Баранова А.А., Вихарева Е.В., Селянинов А.А. Разработка кинетической модели процесса сушки мягких желатиновых капсул с гидрофильным наполнителем // Вестник Пермской государственной фармацевтической академии. - 2011. - № 8. - С. 197-199.

3. Вихарева Е.В., Селянинов А.А., Данилов Ю.Л., Рудакова И.П., Нечеухина Т.А., Ившина И.Б., Няшин Ю.И. Математическая модель процесса биодеструкции парацетамола как открытой системы // Российский журнал биомеханики. - 2008. - Т. 12, № 2. - С. 41-54.

4. Вихарева Е.В., Селянинов А.А., Ившина И.Б., Няшин Ю.И. Математическое моделирование процесса биодеструкции парацетамола актинобактериями рода Rhodococcus // Российский журнал биомеханики. - 2007. - Т. 11, № 2. - С. 93-100.

5. Вихарева Е.В., Селянинов А.А. Кинетическая модель процесса биодеструкции и-аминофенола // Вопросы биологической, медицинской и фармацевтической химии. - 2010. - № 12. - С. 29-33.

6. Вихарева Е.В., Селянинов А.А., Няшин Ю.И., Ившина И.Б. Кинетическая схема процесса биодеструкции парацетамола с истекшим сроком годности // Российский журнал биомеханики. -2006. - Т. 10, № 3.- С. 72-79.

7. Гилязева Р.Ф., Селянинов А.А., Мартюшева М.В., Подгаец Р.М. Механические аспекты шинирования зубного ряда // Российский журнал биомеханики. - 2009. - Т. 13, № 3. - С. 63-72.

8. Государственная Фармакопея СССР: Вып. 2: Общие методы анализа. Лекарственное растительное сырье / ред. коллегия: М. Д. Машковский, Э.А. Бабаян, А.Н. Обоймакова, В.М.Булаев, В. А. Северцев, Б.И. Любимов, С. Д. Соколов, А.И. Тенцова. - 11-е изд., доп. - М.: Медицина, 1990. - С. 143-145.

9. Денисов Е.Т. Кинетика гомогенных химических реакций. - М.: Высшая школа, 1988. - 390 с.

10. Дубинин А.Л., Селянинов А.А. Влияние нижнечелюстного канала на состояние нижней челюсти при установке имплантата // Вестник Пермского национального исследовательского политехнического университета. Прикладная математика и механика. - 2012. - № 10. - С. 63-76.

11. Иванова Н.А., Разепина Я.А., Селянинов А.А., Вихарева Е.В., Сульдин А.В. Кинетика процесса сушки мягких желатиновых капсул с гидрофильными наполнителями // Фармация. - 2012. - № 3. -С. 40-42.

12. Каталог штаммов Региональной профилированной коллекции алканотрофных микроорганизмов / под ред. И.Б. Ившиной. - М.: Наука, 1994. - 163 с.

13. Нигматуллина С.В., Селянинов А.А. Собственные частоты имплантата стремени поршневого типа // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2012. - № 3. - С. 136-148.

14. Селянинов А.А., Вихарева Е.В. Кинетика биодеструкции лекарственных средств - производных фенола, изохинолина и карбоновых кислот // Российский журнал биомеханики. - 2010. - Т. 14, № 2. - С. 79-91.

15. Селянинов А.А., Еловиков А.М., Бородулина Т.С., Подгаец Р.М. Выбор параметров протеза стремени при стапедопластике на основе собственных частот II Российский журнал биомеханики. -2009. - Т. 13, № 4. - С. 42-53.

16. Селянинов А.А., Кислых Ф.И., Подгаец Р.М., Могилат Ю.Ю., Тузова Е.А., Хайрутдинова Ф.Ф. Биомеханические аспекты замещения дефекта нижней челюсти человека имплантатом II Российский журнал биомеханики. - 2003. - Т. 7, № 4. - С. 22-33.

17. Селянинов А.А. Статистическая механика и теория надежности. - Пермь: Изд-во Перм. гос. техн. ун-та, 2008. - 201 с.

18. Селянинов А.А., Тотьмянина А.В., Подгаец Р.М. Биомеханическое сопровождение коррекции зубного ряда с применением эластопозиционеров II Российский журнал биомеханики. - 2012. - Т. 1б, № 1. - С. 57-79.

19. Тотьмянина А.В., Селянинов А.А. Биомеханический анализ коррекции зубного ряда с применением эластопозиционеров II Вестник Пермского государственного технического университета. Прикладная математика и механика. - 2011. - № 9. - С. 183-199.

20. Тотьмянина А.В., Селянинов А.А. Кинетические уравнения перемещения зубов в эластопозиционере при протрузии центральных верхних резцов II Вестник Пермского государственного технического университета. Прикладная математика и механика. - 2012. - № 10. - С. 209-217.

CLASS OF BIOMECHANICAL RANDOM PROCESSES INVOLVING KINETIC SIMULATION

А.А. Selyaninov (Perm, Russia)

The biomechanical processes associated with human body and biological systems have a large body of factors determining their development in the course of time; therefore these processes are essentially random. We suggest discerning a class of biomechanical random processes that allow kinetic simulation. The main features of this class are the certain monotony under increasing or decreasing the characteristic parameter of the process, and the availability of modern methods of its registration. That allows performing a kinetic simulation of the random process realisations. As a result, the non-stationary random process is reduced to a non-random function of two random variables instead of an infinite number of random variables in the V.S. Pugachev canonical decomposition method. The time dependences of the mathematical expectation and variance were derived analytically for the processes that may be simulated by the first order kinetic equations with a linearly varying reaction factor. The determination of the numerical characteristics such as analytical functions of time has also been possible at change of biodegradation rate for the new perception of the process. We propose to introduce the hypothesis of normal distribution of determining parameter of the process for each fixed moment of time for the advanced analysis of the probable unsteady random processes; the current time is taken a parameter. Of course, this hypothesis requires at least of information about the same ranges of decreasing-increasing of the activity of actinobacteria at the start of the biodegradation process of the drug. The proposed approach allows simplifying and making clear the statistical analysis of these processes in vitro, and also getting an novel method for biomechanical accompaniment of the specific realisations in vivo.

Key words: non-stationary random processes, biomechanical processes, kinetic simulation, the function of two random variables, clear statistical analysis, biomechanical accompaniment of realisations.

Получено 15 августа 2012

i Надоели баннеры? Вы всегда можете отключить рекламу.