Full Text

Введение Основным результатом работы является развитие идеи кинетического моделирования реализаций нестационарного случайного процесса при их «ламинарности» и достаточной гладкости [1, 3-7, 10, 13]. Такие процессы считаем простыми. Кинетические уравнения обладают ограниченным числом параметров, которые при появлении каждой новой реализации при одинаковых, казалось бы, начальных условиях процесса в экспериментах на повторяемость принимают новые значения, т.е., по сути, являются случайными величинами, образующими систему случайных величин. Методами статистики (чаще микростатистики) определяются числовые характеристики системы ограниченного числа случайных величин. В результате нестационарный случайный процесс представляется как обычная функция системы ограниченного числа случайных величин, что принципиально упрощает определение его числовых характеристик [11-12]. На основе такого представления можно предложить два приближенных подхода к определению числовых характеристик процесса: линеаризацию и полный анализ при известном или гипотетическом законе распределения системы случайных величин. Линеаризация позволяет определить приближенные аналитические выражения для математического ожидания и дисперсии независимо от коррелированности случайных величин в системе. Зачастую этим можно ограничиться и при известном законе распределения определяющего процесс параметра получить области реализаций процесса с определенной вероятностью. При гипотезе о нормальном законе распределения авторами был произведен такой анализ одного из случайных процессов, который показал выход части области реализаций в недопустимые значения определяющих процесс параметров [12]. Был получен некоторый усеченный вариант применения нормального закона распределения. В качестве второго приближенного подхода предлагается сразу ввести гипотезу о нормальном законе распределения имеющейся системы случайных величин (при их некоррелированности) для определения полного набора числовых характеристик случайного процесса. В работе доказана и применена лемма для вывода аналитических выражений для числовых характеристик случайного процесса. Помимо этого получены числовые характеристики скорости процесса. Отмечается тот факт, что уже в выражение для математического ожидания входят дисперсии случайных величин из системы, что означает наличие ограничений на дисперсии при использовании нормального закона или метода линеаризации. Такие ограничения были получены, также была проверена сходимость результатов обоих подходов при иллюстрации методики на одном из биотехнологических процессов утилизации лекарственных средств. 1. Кинетическая модель реализаций Обозначим через x(t) определяющий параметр, характеризующий случайный процесс, как функцию времени, через x0 - начальное значение параметра. Пусть этот параметр монотонно убывает со временем в каждой реализации и зависимость выражения ln(x/x0) от времени близка к линейной , тогда реализации случайного процесса описываются кинетическим уравнением первого порядка (1) где так называемая константа реакции K может быть функцией времени, в частности для процессов с участием живых организмов - линейной: K(t) = at + b [12]. Здесь a и b - параметры кинетического уравнения. Следует отметить, что кинетическое моделирование не является одним из способов аппроксимации экспериментальных данных, а предполагает их осознанный анализ. Интегрируя (1) с начальным условием x = x0 при t = 0, получаем выражение (2) где от реализации к реализации параметры a и b представляют собой систему двух случайных величин Aи B. Числовые характеристики системы определяются известными выражениями [3, 8]. При малом числе реализаций используются их выборочные аналоги: в качестве математических ожиданий - средние значения , ; в качестве дисперсий - выборочные дисперсии , ; в качестве средних квадратичных отклонений - стандартные отклонения , ; аналогично - корреляционная матрица . 2. Случайный процесс как функция системы ограниченного числа случайных величин Если в зависимости реализаций (2) под a и b понимаются случайные величины A и B, тогда можно записать аналогичное выражение для рассматриваемого нестационарного случайного процесса: (3) В результате нестационарный случайный процесс на основе кинетического моделирования реализаций сводится к функции системы двух случайных величин вместо функции (в виде ряда) системы бесконечного числа пар некоррелированных случайных величин с нулевыми математическими ожиданиями и одинаковыми для пар дисперсиями согласно методу канонического разложения нестационарного случайного процесса В.С. Пугачева [2, 8-9]. Следует отметить, что таким образом нестационарный случайный процесс сводится именно к элементарной функции системы двух случайных величин, а не к функции одной пары некоррелированных случайных величин с нулевыми математическими ожиданиями и одинаковыми для пары дисперсиями. В случае коррелированности величин A и B в дисперсии случайного процесса учитывается слагаемое с корреляционным моментом Kab. При стохастическом анализе интересуют числовые характеристики случайного процесса: математическое ожидание и дисперсия и среднее квадратичное отклонение σx(t), а в полном объеме корреляционные функции самого процесса и скорости его течения. При малом числе реализаций в качестве математического ожидания mx используется среднее значение , в качестве дисперсии Dx - выборочная дисперсия , в качестве среднего квадратичного отклонения - стандартное отклонение . 3. Числовые характеристики случайного процесса Для определения числовых характеристик случайного процесса можно использовать для достоверности два приближенных подхода: известный метод линеаризации [2, 9] и изначальное введение гипотезы о нормальном законе распределения исходной системы случайных величин. Для рассматриваемых в рамках данной работы простых процессов такой анализ производится аналитически. 3.1. Числовые характеристики случайного процесса на основе линеаризации Зависимость Х(t) от системы двух случайных величин A и B (3) является экспоненциальной, ее можно линеаризовать в достаточно малой окрестности средних выборки и , тогда среднее выборки (аналог математического ожидания) случайного процесса определится той же самой функциональной зависимостью (3) (время t при этом является параметром) [12]: (4) Выборочная дисперсия линеаризованной функции системы двух случайных величин принимает вид (5) откуда определяется стандартное отклонение . Пусть Kab = 0 (случайные величины A и B некоррелированы), хотя это непринципиально для метода линеаризации, тогда с применением выражения (2) из (5) получаем зависимость выборочной дисперсии от времени (6) откуда следует выражение для стандартного отклонения (7) где среднее выборки определено выражением (4). В результате на основе кинетических уравнений удалось представить нестационарный случайный процесс функцией системы двух случайных величин и определить его основные числовые характеристики (4), (6) и (7). Для данного подхода не предполагается тип закона распределения системы двух случайных величин, что не позволяет определить корреляционную функцию . Помимо числовых характеристик случайного процесса интересуют области реализаций с определенной вероятностью. Если число реализаций случайного процесса мало, то затруднительно использовать известные методы нахождения интересующего закона распределения. Исходя из желания использовать известное правило «трех сигм», можно попытаться принять гипотезу о нормальном законе распределения определяющего процесс параметра, как это сделано в работе [12]. При этом полученные методом линеаризации числовые характеристики случайного процесса подразумевают при экспоненциальной зависимости процесса от системы случайных величин A и B определенные ограничения на дисперсии и . Следует обратить внимание на то, что добавление к математическому ожиданию даже одной «сигмы» может вывести случайный процесс за пределы реально возможных реализаций. Это будет противоречить физической сущности процесса, что свидетельствует о несоответствии полученных областей реализаций заявленным вероятностям, т.е. реализуется некоторый усеченный вариант применения нормального закона распределения. 3.2. Числовые характеристики случайного процесса при гипотезе о нормальном законе распределения системы случайных величин Эта гипотеза определяет второй приближенный метод анализа рассматриваемого случайного процесса. Целью является его сравнение с методом линеаризации и определение условий сходимости для достоверности результатов. Если при линеаризации закон распределения определяющего процесс параметра может потребоваться только на этапе нахождения областей реализаций случайного процесса, то для полного анализа следует определить корреляционную функцию. Для этого изначально требуется информация о законе распределения системы случайных величин A и B и возможность определения частных законов распределения, т.е. необходима некоррелированность случайных величин в системе (Kab = 0). В связи с этим сразу введем гипотезу о нормальном законе распределения системы двух некоррелированных случайных величин и получим на основе кинетического моделирования реализаций полный набор числовых характеристик нестационарного случайного процесса. Для этого с целью удобства записи выражений используются термины математического ожидания и дисперсии, но под ними при необходимости можно понимать среднее выборки и выборочную дисперсию. Для вывода выражений числовых характеристик сформулируем следующую лемму. Лемма. Если , где с - случайная величина, распределенная по нормальному закону f(c), τ - неслучайная величина, тогда математическое ожидание случайной величины с определяется выражением (8) где mc - математическое ожидание; Dc - дисперсия случайной величины с. Замечание. Следует отметить, что в выражение (8) для математического ожидания функции случайной величины помимо математического ожидания случайной величины входит ее дисперсия. Доказательство. Производится элементарным интегрированием выражения для математического ожидания . Следствие. Если , где a, b - независимые случайные величины, распределенные по нормальным законам; x0, t - неслучайные величины, тогда числовые характеристики неслучайной функции системы двух случайных величин примут вид: (9) (10) (11) Доказательство. Рассмотрим выражение для математического ожидания неслучайной функции : , где , . Далее два раза применим лемму: при c = b, и при c = а, . В результате получим выражение (9). Рассмотрим выражение для дисперсии случайного процесса . Далее два раза применим лемму: при c = b, и при c = а, . В результате получим выражение (10). Для определения корреляционной функции рассмотрим произведение , где , . Сначала определим первое слагаемое в корреляционной функции . После чего применим два раза лемму: при c = b, и при c = а, . Добавив второе слагаемое с учетом выражения (9) получим окончательное выражение для корреляционной функции (11). Замечание 1. Корреляционная функция при является дисперсией. Сравнение с целью проверки вывода выражений (10) и (11) при дает равенство . Замечание 2. Сравним результаты линеаризации выражений (9) и (10) для математического ожидания и дисперсии полученные с изначально введенной гипотезой о нормальном распределении системы двух случайных величин, с результатами линеаризации функции системы двух случайных величин без гипотезы о распределении. Линеаризация выражений (9) и (10) состоит в разложении , по малым параметрам и , после чего линеаризованные выражения принимают вид , (12) (13) Эти разложения справедливы, если выполняется ограничение на дисперсии (14) Для сравнения приведем результаты линеаризации функции системы двух случайных величин без гипотезы о законе распределения из предыдущего раздела: (15) (16) Вывод по замечанию 2. В выражении (15) нет слагаемых с дисперсиями Da и Db по сравнению с (12), они совпадают при выполнении ограничений на дисперсии (14), а выражение (16) полностью совпадает с выражением (13). 4. Числовые характеристики скорости случайного процесса как функции системы случайных величин В качестве дальнейшего развития идеи о применении кинетического моделирования при исследовании данного случайного процесса следует отметить, что изменение скорости нестационарного процесса Х(t) со временем также является случайным процессом V(t). В качестве определяющего процесс параметра в этом случае можно выбрать скорость процесса v(t), для которой из (1) и (2) следует выражение (17) Определим числовые характеристики скорости V(t) данного процесса: среднее значение выборочную дисперсию и стандартное отклонение . Повторив процедуру линеаризации для скорости, выбранной в качестве определяющего процесс параметра, получим среднее значение реализаций (18) выборочную дисперсию (19) и стандартное отклонение для скорости процесса (20) Считая, что оба приближенных метода сошлись при выполнении ограничений на дисперсии, для определения корреляционной функции можно использовать гипотезу о нормальном распределении системы случайных величин. Тогда вместо применения леммы для определения корреляционной функции скорости процесса можно ввиду гладкости дважды продифференцировать уже известную корреляционную функцию для процесса , т.е. (21) откуда в окончательном виде (22) где Выражение для корреляционной функции скорости процесса проверили с помощью пакета программ MATLAB 7.0.1. Из выражения (22) при равенстве определяется дисперсия скорости процесса (23) Продифференцируем математическое ожидание нестационарного случайного процесса по времени (24) тогда выражение для математического ожидания скорости случайного процесса примет вид или (25) 5. Применение теории простых процессов для анализа повторяемости биологической деструкции дротаверина гидрохлорида (но-шпы) Рассмотрим применение данной методики для анализа повторяемости экспериментов по биологической деструкции дротаверина гидрохлорида, которая отвечает условиям простого нестационарного случайного процесса. Определяющими процесс параметрами являются концентрация x(t) (%) и скорость изменения концентрации v(t) (%/ сут) дротаверина гидрохлорида в культуральной жидкости как функции времени. Начальная концентрация 100 % соответствует начальной концентрации утилизируемого лекарственного средства (0,002 мг/мл). В качестве исходной информации для применения разработанной методики использовали экспериментальные данные по текущей концентрации дротаверина гидрохлорида в опытах на повторяемость процесса биологической деструкции при минимальной выборке, приведенные в работе [12] (табл. 1). Таблица 1 Изменение содержания дротаверина гидрохлорида в процессе биологической деструкции, % t, сут 0 2 5 20 30 x1 100 40,50 27,90 9,10 1,02 x2 100 7,30 3,75 2,80 0,98 x3 100 - 10,20 - 1,00 В табл. 1 приведены данные по трем реализациям одного и того же процесса биологической деструкции. При одинаковых, казалось бы, условиях эксперимента реализации отличаются друг от друга, т.е. процесс биологической деструкции носит случайный характер. Следует отметить, что в последней реализации отсутствуют данные на 2-е и 20-е сутки, что представляет определенные трудности для обычного статистического анализа. По данным табл. 1 не наблюдается осцилляций концентрации во времени, т.е. зависимости достаточно гладкие и не перекрещиваются. Изменения концентрации дротаверина гидрохлорида во времени в каждой реализации случайного процесса являются убывающими функциями, поэтому для кинетического моделирования применимы описанные ранее уравнения первого порядка (1). Параметры b и a в уравнениях для кинетических кривых (2) и (17) можно получить с применением метода наименьших квадратов. На рис. 1 приведены кинетические кривые реализаций процесса биологической деструкции для двух определяющих процесс параметров: концентрации дротаверина гидрохлорида (рис. 1, а) и скорости изменения его концентрации (рис. 1, б). а б Рис. 1. Реализации процесса биологической деструкции при экспериментах на повторяемость после кинетического моделирования: - экспериментальные точки; 1, 2, 3 - номер эксперимента Характерной особенностью данного процесса биологической деструкции является его высокая скорость в начале, постепенно убывающая со временем. Как следует из рис. 1, кинетические кривые достаточно адекватно описывают эксперименты на повторяемость изменения концентрации дротаверина гидрохлорида при биологической деструкции. Для случайного процесса (3), включающего все реализации, ускорение а и скорость процесса b приобретают смысл системы случайных величин A и B. Числовые характеристики системы двух случайных величин A и B (средние выборки, выборочные дисперсии, стандартные отклонения и корреляционную матрицу) определили по известным выражениям статистики. Результаты расчета числовых характеристик приведены в табл. 2. Следует обратить внимание на равенство нулю значения члена корреляционной матрицы, что отражает некоррелированность случайных величин A и B. Таблица 2 Числовые характеристики системы двух случайных величин А и B 1/сут2 1/сут 1/сут4 1/сут2 Kab -0,037 9 0,663 8 0,000 891 0,140 322 0 Согласно выражениям (4), (7) и (18), (20), полученным с применением метода линеаризации, на рис. 2 приведены зависимости основных числовых характеристик процесса от времени и его скорости при экспериментах на повторяемость. а б Рис. 2. Зависимости от времени средних значений и стандартных отклонений при экспериментах на повторяемость: а - процесса биологической деструкции; б - скорости процесса Максимумы и минимумы стандартных отклонений отражают поведение кинетических кривых, приведенных на рис. 1. Корреляционные функции самого процесса и его скорости, определенные выражениями (11) и (22), при значении дисперсий Da и Db из табл. 2 расходятся (стремятся к бесконечности в пределах первых 6 сут). Этого следовало ожидать из-за больших разбросов реализаций случайного процесса (табл. 1). Для корректного корреляционного анализа данного процесса биологической деструкции с применением гипотезы о нормальном законе распределения системы двух случайных величин действительно необходимо выполнение ограничений на дисперсии, приведенных в выражении (14). Путем постепенного уменьшения дисперсий Da и Db получили сходимость математических ожиданий и , дисперсий и , полученных с применением метода линеаризации и с применением гипотезы о нормальном законе распределения. Корреляционные функции для процесса и для его скорости, сошедшиеся при уменьшении дисперсий Da и Db в 280 раз, приведены на рис. 3. Штриховыми линиями на поверхностях (рис. 3) выделены зависимости от времени дисперсий и , полученные при равенстве в выражениях корреляционных функций для процесса и его скорости а б Рис. 3. Корреляционные функции нестационарного случайного процесса биологической деструкции: а - ; б - при экспериментах на повторяемость Корреляционная функция при t = 0 равна нулю, так как во всех реализациях процесс биологической деструкции идет с одной начальной концентрации. При увеличении времени до 2 сут для процесса наблюдается небольшой подъем из-за большого разброса реализаций (рис. 3, а) с дальнейшим уменьшением до нулевого значения. При корреляционная функция представляет собой зависимость дисперсии концентрации дротаверина гидрохлорида при экспериментах на повторяемость процесса от времени согласно выражению (16). Для скорости процесса корреляционная функция сначала имеет максимальное значение из-за максимальных в начале процесса разбросов реализаций скоростей биологической деструкции. При 2 сут корреляционная функция принимает нулевое значение, что соответствует практически нулевой дисперсии на данный момент времени, затем небольшой подъем (рис. 3, б) с дальнейшим уменьшением до нуля. При дисперсия скорости изменяется в соответствии с выражением (19). Корреляционный анализ показал, что по изменению концентрации процесс биологической деструкции дротаверина гидрохлорида обладает особенной чувствительностью к возмущениям на 2-е сутки, а по скорости изменения концентрации определяющее значение имеет его начало. Выводы 1. Кинетически моделируемые случайные процессы можно представить обыкновенной функцией системы случайных величин, которые являются параметрами кинетических уравнений в реализациях. 2. При гипотезе о нормальном законе распределения системы случайных величин приведена и доказана лемма, следствия из которой позволяют определить аналитические выражения для числовых характеристик случайного процесса и его скорости. 3. Показана необходимость учета ограничений на дисперсии случайных величин при использовании гипотезы о нормальном законе распределения. 4. Разработанная методика вероятностного анализа реализована для процесса биологической деструкции дротаверина гидрохлорида (но-шпы) при экспериментах на повторяемость.

About the authors

A. A Selianinov

Perm National Research Polytechnic University

M. A Osipenko

Perm National Research Polytechnic University

A. A Baranova

Perm State Pharmaceutical Academy

E. V Vikhareva

Perm State Pharmaceutical Academy

A. N Khrenkov

Perm State Pharmaceutical Academy

References

Statistics

Views

Abstract - 29

PDF (Russian) - 15

Refbacks

• There are currently no refbacks.