Численные методы решения жестких систем стохастических дифференциальных уравнений

Основные принципы построения численных методов решения стохастических дифференциальных уравнений (СДУ). Определение жесткой системы СДУ. Анализ основных свойств: устойчивость, порядок сходимости и точность аппроксимации. Метод решения систем жестких СДУ.

Рубрика Математика
Вид статья
Язык русский
Дата добавления 27.11.2018
Размер файла 94,1 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Размещено на http://www.allbest.ru/

Труды ИСА РАН 2008. Т. 32 (2)

Динамика неоднородных систем, 2008

Численные методы решения жестких систем стохастических дифференциальных уравнений

Ерешко Арт.Ф., Филатова Д.В.

Вычислительный центр им.А. А. Дородницына РАН

Свентокшиская Академия в Кельцах, Польша

Аннотация

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

1. Принципы построения численных методов решения стохастических дифференциальных уравнений

Очень не многие СДУ имеют аналитические решения. В основном это - линейные СДУ с аддитивным или мультипликативным шумом или нелинейные СДУ, сводимые к линейным [1]. В связи с этим возникает практическая необходимость создания численных методов решения СДУ.

В настоящее время существует несколько подходов создания численных схем решения стохастических дифференциальных уравнений. Одной из возможностей является адаптация существующих для ОДУ схем с учетом свойств стохастических интегралов [2], другой - разработка специальных методов решения СДУ [3]. Большинство исследователей использует первый подход, поскольку теория численного решения ОДУ хорошо разработана и достаточно легко можно провести аналогии между ОДУ и СДУ.

Самым простым методом аппроксимации численного решения СДУ (с вычислительной точки зрения) является метод Эйлера, разработанный Маруямой в 1955г [4]. Эта схема удовлетворяет многим необходимым свойствам, предъявляемым к численным методам (она имеет порядок сходимости 0.5), но в тоже время обладает рядом ограничений (не всегда устойчива, ошибка аппроксимации достаточно высока и т.п.). Для устранения этих недостатков, а также повышения порядка сходимости численных схем решения СДУ были проведены и до сих пор ведутся исследования, направления которых можно представить в виде схемы (см. рис.1).

стохастическое дифференциальное уравнение жесткая система

Рис. 1. Направления создания численных методов решения СДУ

По аналогии с разработкой схем численного решения ОДУ для повышения порядка сходимости, точности аппроксимации и устойчивости можно использовать разложение в точке аппроксимации, т.е. использования производных различных порядков как переменной, так и коэффициентов дрейфа и диффузии. В литературе этот подход получил название метода Тейлора [5]. Однако, минусом схем Тейлора является то, что на каждом шаге аппроксимации требуется вычислять кратные стохастические интегралы, связанные с вышеуказанными производными. Для того, чтобы избежать вычислительные трудности можно использовать многократное деление шага аппроксимации (методы Рунге-Кутта [6]) или результаты аппроксимации предыдущих шагов (многошаговые методы [7]).

Напомним, что разработка численной схемы или алгоритма решения связана не только с математическими аспектами, а также с физической природой объекта исследования и программной реализацией (рис.1 и рис.2).

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

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

Анализ современной литературы показал, что создание численных методов решения жестких систем в большинстве случаев основано на идеях, представленных Хайрером и Ванером [8]. В своей работе они постулировали, что жесткие системы не могут быть решены явными методами, и представили подходы, основанные только на использовании неявных методов. Однако, следует отметить, что непосредственное применение этих методов всегда связано с крайне сложной процедурой определения параметров схемы, основанной на заранее выделенной области устойчивости только для рассматриваемой системы. Это обстоятельство делает предложенные подходы не приемлемыми для большинства вышеуказанных приложений, но позволяет выделить два важных математических свойства жесткости. Во-первых, все жесткие системы обладают очень широким спектром (или присутствием очень разных экспонент Ляпунова). Во-вторых, согласно теоремы единственности и существования решения, для жестких систем характерны большие значения константы Липшица.

Итак, анализ принципов создания численных схем решения систем СДУ, особенно для некорректных задач, показал необходимость поиска новых метод решения. Однако, прежде чем начать разработку собственного метода решения жесткой задачи, следует исследовать свойства уже существующих методов и проанализировать их реакцию при наличии жесткости. Поэтому, поскольку существует достаточно большая группа, основанная на разложениях Тейлора, следует пояснить такие понятия как разложение стохастического интеграла уравнения (1.15), а также методы численной аппроксимации этих интегралов.

2. Разработка явного сильного метода - преобразования для решения жестких систем

Рассмотрим следующую систему

(1)

где , и . Пусть параметры заданы как , , . Для получения решения системы применим стандартную явную схему Рунге-Кутта. Результаты представлены на рис.2.

а) фазовый портрет

б)

в)

Рис.2. Численное решение системы (1)

Как можно заметить, численное решение системы (1) является нестабильным. Для улучшения численного решения используем шкалирование переменной время, вводя новую переменную [9]

, (2) где

,

Преобразуем теперь исходную систему (1), вводя дополнительное уравнение (2)

и решим новую систему еще раз для , и , используя явный метод Рунге-Кутта.

Полученные результаты свидетельствуют о значительной стабилизации результатов (см. рис.3).

Основной результат этого эксперимента можно представить следующим утверждением: пусть для исходной жесткой системы ОДУ

, , , (3)

введено вспомогательное уравнение так, что

(4)

где ,

тогда численное решение системы (4), получаемое любым явным методом, является более стабильным, чем для системы (3).

а) фазовый портрет

б)

в)

Рис. 3. Численное решение системы (1) с применением преобразования

Используем аналогичный подход для улучшения численного решения системы жестких СДУ. Рассмотрим систему, состоящую из СДУ

, , (5)

где - функция дрейфа, - функция диффузии, - приращения некоррелированных процессов Винера.

Введем новую переменную, шкалирующую переменную время

, (6)

и преобразуем систему (5) к виду

(7)

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

(8)

с параметрами , , , и .

Численное решение системы было выполнено с применением сильной явной схемы Рунге-Кутта без использования (см. рис.4. а) и с использованием преобразования (см. рис.4. б)

а) без использования преобразования

б) с использованием преобразования

Рис. 4. Фазовый портрет системы

Основной результат можно сформулировать как [10]: пусть для исходной жесткой системы СДУ

, , , (9)

введено вспомогательное уравнение так, что

(10)

где ,

тогда численное решение системы (10), получаемое любым явным методом, является более стабильным, чем для системы (9).

Литература

1. Oksendal B. Stochastic differential equations. - Berlin: Springer. - 2000. - 326p.

2. Burrage K., Burrage P., Mitsui T. Numerical solutions of stochastic differential equations - implementation and stability issues // Journal of computational and applied mathematics. - 2000. - Vol.125. - P.171 - 182.

3. Gaines J. G., Lyons T. J. Variable step size control in the numerical solution of stochastic differential equations // SIAM Journal of Applied Mathematics. - 1997. - Vol.57 (5). - P.1455 - 1484.

4. Higham D., Mao X., Stuart A. Strong convergence of Euler-type methods for nonlinear stochastic differential equations // Journal of Numerical Analysis. - 2002. - Vol.40 (3). - P.1041 - 1063.

5. Butcher J. C. Numerical methods for ordinary differential equations. - West Sussex: Wiley. - 2003. - 425 p.

6. Шалашилин В.И., Кузнецов Е.Б. Метод продолжения решения по параметру и наилучшая параметризация (в прикладной математике и механике). - М.: Эдиториал УРСС. - 1999. - 224 с.

7. Хайрер Э., Ванер Г. Решение обыкновенных дифференциальных уравнений - Жесткие и дифференциально-алгебраические задачи. - М.: Мир. - 1999. - 486с.

8. Васильев Ф.П. Численные методы решения экстремальных задач. - М.: Наука. - 1980. - 520с.

9. Filatova D., Grzywaczewski M. The problems of numerical methods for stiff stochastic systems // Abstracts of 5th Workshop of the ERCIM Working Group on Matrix Computations and Statistics: Numerical Methods for Statistics, 27 - 29.08.2004, Prague, Czech Republic. - P.7.

Размещено на Allbest.ru

...

Подобные документы

  • Определение дифференциальных уравнений в частных производных параболического типа. Приведение уравнения второго порядка к каноническому виду. Принцип построения разностных схем. Конечно-разностный метод решения задач. Двусторонний метод аппроксимации.

    дипломная работа [603,8 K], добавлен 24.01.2013

  • Анализ методов решения систем дифференциальных уравнений, которыми можно описать поведение материальных точек в силовом поле, законы химической кинетики, уравнения электрических цепей. Этапы решения задачи Коши для системы дифференциальных уравнений.

    курсовая работа [791,0 K], добавлен 12.06.2010

  • Определение и анализ многошаговых методов, основы их построения, устойчивость и сходимость. Постановка задачи Коши для обыкновенных дифференциальных уравнений. Метод Адамса, значение квадратурных коэффициентов. Применение методов прогноза и коррекции.

    контрольная работа [320,8 K], добавлен 13.03.2013

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

    контрольная работа [138,9 K], добавлен 08.01.2016

  • Методы оценки погрешности интерполирования. Интерполирование алгебраическими многочленами. Построение алгебраических многочленов наилучшего среднеквадратичного приближения. Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений.

    лабораторная работа [265,6 K], добавлен 14.08.2010

  • Изучение численных методов приближенного решения нелинейных систем уравнений. Составление на базе вычислительных схем алгоритмов; программ на алгоритмическом языке Фортран - IV. Приобретение практических навыков отладки и решения задач с помощью ЭВМ.

    методичка [150,8 K], добавлен 27.11.2009

  • Решение эллиптических и параболических дифференциальных уравнений в частных производных. Суть метода Кранка-Николсона и теории разностных схем для теплопроводности. Построение численных методов с помощью вариационных принципов, описание Matlab и Mathcad.

    курсовая работа [1,4 M], добавлен 13.03.2011

  • Обобщенные решения линейных дифференциальных уравнений. Основные примеры построения фундаментальных решений линейных дифференциальных операторов с постоянными коэффициентами, метод преобразования Фурье. Преимущества использования методов спуска.

    курсовая работа [1,1 M], добавлен 10.04.2014

  • Дифференциальные уравнения как математический инструмент моделирования и анализа разнообразных явлений и процессов в науке и технике. Описание математических методов решения систем дифференциальных уравнений. Методы расчета токов на участках цепи.

    курсовая работа [337,3 K], добавлен 19.09.2011

  • Численные методы решения систем линейных уравнений: Гаусса, простой итерации, Зейделя. Методы аппроксимации и интерполяции функций: неопределенных коэффициентов, наименьших квадратов. Решения нелинейных уравнений и вычисление определенных интегралов.

    курсовая работа [322,7 K], добавлен 27.04.2011

  • Метод Эйлера: сущность и основное содержание, принципы и направления практического применения, определение погрешности. Примеры решения задачи в Excel. Метод разложения решения в степенной ряд. Понятие и погрешность, решение с помощью метода Пикара.

    контрольная работа [129,0 K], добавлен 13.03.2012

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

    дипломная работа [395,4 K], добавлен 10.06.2010

  • Общая постановка задачи решения обыкновенных дифференциальных уравнений, особенности использования метода Адамса в данном процессе. Решение системы обыкновенных дифференциальных уравнений методом Адамса и точным методом, сравнение полученных результатов.

    курсовая работа [673,6 K], добавлен 27.04.2011

  • Решение дифференциальных уравнений в частных производных. Метод минимальных невязок, минимальных поправок, скорейшего спуска, сопряженных градиентов. Алгоритмы и блок-схемы решения. Руководство пользователя программы. Решение системы с матрицей.

    курсовая работа [380,3 K], добавлен 21.01.2014

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

    курсовая работа [1,7 M], добавлен 03.11.2011

  • Рассмотрение теории дифференциальных уравнений. Выделение классов уравнений с систем, решения которых не имеют подвижных критических особых точек. Установление достаточности найденных условий путем сравнения с классическими системами типа Пенлеве.

    курсовая работа [137,0 K], добавлен 01.06.2015

  • Основные понятия и теоремы систем линейных уравнений, характеристика методов их решения. Критерий совместности общей системы. Структура общих решений однородной и неоднородной систем. Матричный метод решения и обобщение. Методы Крамера и Гаусса.

    курсовая работа [154,5 K], добавлен 13.11.2012

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

    курсовая работа [820,5 K], добавлен 22.08.2010

  • Дифференциальное уравнение первого порядка, разрешенное относительно производной. Применение рекуррентного соотношения. Техника применения метода Эйлера для численного решения уравнения первого порядка. Численные методы, пригодные для решения задачи Коши.

    реферат [183,1 K], добавлен 24.08.2015

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

    дипломная работа [1,2 M], добавлен 16.12.2008

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.