Идентификация параметров осциллирующих процессов в живой природе моделируемых дифференциальными (стр. 5 из 5)
, (9) 
, (10)где
- фиксированные весовые коэффициенты.Градиентные уравнения и соответствующие начальные условия для этих функционалов следующие:
, (11) 
, (12) 
(13) 
, 
(14)3.3 Как ускорить вычисления
Опыт реальных вычислений показывает, что минимизация функционала методом градиентных уравнений естественно делится на два этапа. На первом этапе происходит быстрое уменьшение функционала. На втором этапе это уменьшение становится все более медленным, и процесс нахождения достаточно точного приближения параметров, соответствующих локальному минимуму функционала, может потребовать неприемлемо больших затрат машинного времени.
Для того, чтобы ускорить вычисления на втором этапе, необходимо ускорить численное интегрирование исходных уравнений, уравнений в вариациях и градиентных уравнений. Исходные уравнения и уравнения в вариациях, как правило, полиномиальные и для их численного интегрирования можно использовать метод рядов Тейлора.
Градиентные уравнения не полиномиальные, и на первом из упомянутых выше этапов их естественно интегрировать методами Рунге-Кутта. На втором этапе идентифицируемые параметры изменяются медленно и правые части градиентных уравнений можно аппроксимировать полиномами по этим параметрам в окрестности некоторого их текущего значения.
Эта аппроксимация достаточно точна только на некотором промежутке изменения
, поэтому ее нужно время от времени строить заново в окрестности очередного текущего значения параметров. На соответствующих промежутках изменения 
приближенные полиномиальные градиентные уравнения можно интегрировать методом рядов Тейлора.Отметим, что построение каждой аппроксимации градиентных уравнений требует многократного численного решения исходных уравнений и уравнений в вариациях, для чего можно использовать метод рядов Тейлора.
Перейдем к формулам. Уравнения точной градиентной задачи Коши
(15) 
, 
, (16)где
, мы хотим заменить на приближенные градиентные уравнения: 
, 
, (17)где
- полином по 
, а 
- набор его коэффициентов.При этом мы хотим, чтобы величины
, (18)были достаточно малыми при
, (19)где
- некоторое фиксированное число. Коэффициенты 
поли-нома можно найти методом наименьших квадратов с функционалом: 
, (20)где
, 
, а 
- весовые коэффициенты.Отметим, что при малых
в качестве 
можно рассмотреть полином степени 3 или 4, а при больших 
и/или 
- полином степени 2.3.4 Численный эксперимент
Мы опишем здесь постановку и результаты одного из численных экспери-ментов, проведенных в полном соответствии с рассмотренной выше схемой градиентного метода. Эти результаты опубликованы в работе [4].
Обратимся к дифференциальным уравнениям для модели Лотки в п. 3.1 и в численном эксперименте будем действовать по следующей схеме:
1. Фиксируем начальные данные
, 
, 
, 
(21)и параметры
, 
, 
(22)2. При этих значениях начальных данных
и параметров 
численным интегрированием задачи Коши (1),(2) находим значения концентрации реактанта 
в моменты времени 
, 
, то есть находим 
при 
.Теперь можно имитировать «измерения» величин 
по формуле
, , (23)где
- независимые случайные величины, равномерно распределенные меж-ду 
и 
. Считаем, что 
- измерения, полученные в некотором реаль-ном эксперименте.3. Фиксируем начальное приближение:
(24)и методом градиентных уравнений находим приближенное значение точки локального минимума
.Об эффективности метода можно судить по затраченному процессорному времени и по величине относительной погрешности:
(25)Результаты этого численного эксперимента приведены на рисунках 1, 2.
4. О других методах идентификации
Ограничимся здесь ссылкой на электронную статью [5], в которой идентифицируются три неизвестных параметра в пяти кинетических уравнениях, описывающих изменение концентраций в биохимических реакциях с участи-ем различных тромбинов и их комплексов.
В этой работе рассматривается функционал МНК, использующий различные начальные данные, соответствующие измерениям для всех пяти переменных в фиксированные моменты времени
, причем все эти измерения взяты из реальных экспериментов.Для минимизации функционала используется программа VARPRO Стэнфордского университета, а численное интегрирование исходных уравнений (для вычисления функционала) проводится при помощи интегратора SDRIV1 Дэвида Кахане.
Литература
1. В. Вольтерра, «Математическая теория борьбы за существование». Москва. «Наука»,1976.
2. Э. Хайрер, С. Нёрсетт, Г. Ваннер, “Решение обыкновенных дифференциальных уравнений”, I. Нежесткие задачи. Москва. “Мир”,1990.
3. Э. Хайрер, Г. Ваннер, “Решение обыкновенных дифференциальных уравнений”, II. Жесткие и дифференциально - алгебраические задачи. Москва. “Мир”,1999.
4. L.K. Babadzanjanz, J.A. Boyle, D.R. Sarkissian, and J.Zhu, “Parameter Identification for Oscillating Chemical Reactions Modelled by Systems of ODE”, Journal of Computational Methods for Sciences and Engineering, 2002.
5. Bert W. Rust, ACMD, Robert W. Ashton, Chemical Science and Technology Laboratory, “Parameter Identifications”, 7/15/2001: http://math.nist.gov/mcsd/Reports/95/yearly/node28.html
6. R.Haberman, “Mathematical Models. Mechanical Vibrations, Population Dynamics, and Traffic Flow. Classics in Applied Mathematics, 21”, SIAM, Philadelphia, 1977.
7. A.J. Lotka, “Undamped oscillations derived from the law of mass action”, Jour. Amer. Chem. Soc. 42 (1920), 1595-1599.
8. A.J. Lotka, “Elements of Physical Biology”, Williams and Wilkins, Baltimore, 1925.