Метод решения нелинейных дифференциальных уравнений с помощью формализации операций компьютера тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Строганов, Андрей Валентинович

• Строганов, Андрей Валентинович
• кандидат физико-математических наук
• 2012

Оглавление диссертации кандидат физико-математических наук Строганов, Андрей Валентинович

Оглавление

Введение

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

1.1 Различные подходы для получения аналитических и численных

решений дифференциальных уравнений

1.2 Численные методы решения задачи Коши, их достоинства и недостатки, разностный подход

1.3 Проблема связи получения явного и численного решений на основе разностных схем

1.4 Модельная задача для изучения свойств метода

Глава 2. Представление решения в методе компьютерной аналогии

2.1 Системы счисления и представление чисел в вычислительных

устройствах

2.2 т-система счисления

2.3 Представление решения дифференциальных уравнений в виде ряда по степеням шага аргумента

2.5 Операция переноса разрядов

2.6 Оценка отбрасываемой части ряда при приведении решения к отрезку

ряда фиксированной длины

Глава 3. Стохастические свойства решений в методе компьютерной аналогии

3.1.Случайные числа и генераторы случайных чисел. Методы оценки

качества генераторов случайных чисел

3.2. Стохастическое поведение коэффициентов в т-представлении решения

3.3 Моделирование случайного и квазислучайного поведения коэффициентов

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

3.6 Свойства самоподобия величин переноса

Глава 4. Применение метода компьютерной аналогии для решения нелинейных дифференциальных уравнений

4.1 Алгоритм построения решения методом компьютерной аналогии

4.2 Построение неполного решения для задачи (1.5)

4.3 Построение полного решения для задачи (1.5)

4.4 Система нелинейных дифференциальных уравнений

Глава 5. Задачи, не разрешимые в квадратурах или не имеющие решения в элементарных функциях. Выявление качественных свойств решения

5.1 Примеры решения задача Коши для уравнений, неразрешимых в

элементарных функциях

5.2 Задача Коши для уравнения Риккати

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

5.4 Асимптотическое решение

5.5 Осциллятор Ван дер Поля. Исследование поведения решения

Заключение

Литература

Введение диссертации (часть автореферата) на тему «Метод решения нелинейных дифференциальных уравнений с помощью формализации операций компьютера»

Введение

Многие научные, технические, экономические и другие задачи приводят к необходимости решения дифференциальных уравнений. Однако, класс дифференциальных уравнений (см. например [1-3]), допускающих аналитическое решение, достаточно узок. Поэтому обычно не удается избежать численного моделирования. Заметим, что в некоторых случаях, когда аналитическое решение уравнения существует, но требует большого объема алгебраических выкладок, компьютерные методы могут оказаться даже предпочтительнее аналитических, так что эти подходы в принципе способны взаимодополнять друг друга. И хотя применение современной компьютерной техники дает возможность разрешать многие сложные задачи, поиск аналитических (или полуаналитических) подходов не теряет актуальности. В отличие от численного, решение в явном аналитическом виде позволяет получить больше информации о задаче, увидеть решение в целом с его качественными особенностями, асимптотиками и т.д., - что важно, например, для физического описания явлений. Также существенно, что аналитические решения, даже частные, могут являться тестами, играющими большую роль при разработке численных методов.

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

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

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

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

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

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

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

Основные результаты по построению и применению этого подхода отражены в наших работах [4-15].

Опишем кратко содержание работы.

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

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

шага к [п + 1) -му коэффициенты в представлении решения в виде т-системы

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

представления решения.

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

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

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

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Строганов, Андрей Валентинович

Заключение

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

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

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

Рассмотрены простые примеры нелинейных уравнений, допускающих квадратуры, а также уравнений и систем, не имеющих известных аналитических представлений. Безусловно, метод нацелен на применение для сложных нелинейных систем уравнений. Одно из возможных перспективных

Список литературы диссертационного исследования кандидат физико-математических наук Строганов, Андрей Валентинович, 2012 год

Литература

1. Степанов В.В. Курс дифференциальных уравнений. Изд. 8-е, стереотипное. - М.: Едиториал УРСС, 2004. - 472 с.

2. Смирнов В.И. Курс высшей математики. Т.2. - М.: Наука, 1967. - 656 с.

3. Лизоркин П.И. Курс дифференциальных и интегральных уравнений с дополнительными главами анализа. - М.: Наука, 1981. - 381 с.

4. Аристов В.В. Строганов А.В. Построение решений дифференциальных уравнений с помощью метода «компьютерной аналогии» // Доклады академии наук РАН, том. 434, N2, 2010, с. 151-157.

5. Аристов В.В. Строганов А.В. Вероятностные аспекты метода «компьютерной аналогии» для решения дифференциальных уравнений // Компьютерные исследования и моделирование, 2009, Т. 1, N 1, с. 2131.

6. Аристов В.В. Строганов А.В. Перспективы метода решения дифференциальных уравнений на основе компьютерной аналогии // Прикладная информатика, 2008, 6(18), с. 91-99.

7. Aristov V.V., Stroganov A.V. A method of formalizing computer operations for solving nonlinear differential equations // Applied Mathematics and Computation, 2012, Vol. 218, p.8083-8098. doi:10.1016/ j.amc.2011.09.029.

8. Аристов В.В., Строганов А.В. Метод решения дифференциальных уравнений на основе «компьютерной аналогии». Тезисы докладов международной конференции «Тихонов и современная математика». М.: МГУ. 2006 г. —с26-27 .

9. Аристов В.В. Строганов А.В. Новый метод решения дифференциальных уравнений с формализацией математических операций компьютера // Сб. научн. трудов. «Математика. Компьютер. Образование». Ред. Г.Ю.Ризниченко. М.-Ижевск: НИЦ Регулярная и хаотич. динамика. 2006. т.2. с.295-307.

10. Аристов B.B. Строганов A.B. Полуаналитический подход для решения дифференциальных уравнений // 54-я Научно-техническая конференция МИРЭА: Сб. тр. — М.: МИРЭА, 2005. — 4.2. — Физико-математические науки. — с.35-40.

П.Аристов В.В. Строганов A.B. Развитие метода решения дифференциальных уравнений на основе «компьютерной аналогии» // 55-я Научно-техническая конференция МИРЭА: Сб. тр. — М.: МИРЭА, 2006. — 4.2. — Физико-математические науки. — с.92-95.

12. Аристов В.В. Строганов A.B. Решение дифференциальных уравнений на основе метода компьютерной аналогии // 56-я Научно-техническая конференция МИРЭА: Сб. тр. — М.: МИРЭА, 2007. — 4.2. — Физико-математические науки. — с.51-56.

13. Аристов В.В., Строганов A.B. Развитие метода «компьютерной аналогии» для решения дифференциальных уравнений. «Математика. Компьютер. Образование». Ред. Г.Ю. Ризниченко. М.-Ижевск: НИЦ Регулярная и хаотич. динамика. 2008. вып.12, т.2. с. 110-120.

14. Аристов В.В., Строганов A.B. Развитие метода решения дифференциальных уравнений с формализацией математических операций компьютера // Сб. научн. трудов. «Математика. Компьютер. Образование». Ред. Г.Ю. Ризниченко. М.-Ижевск: НИЦ Регулярная и хаотич. динамика. 2010. т.2. с.227-235.

15. Aristov V.V., Stroganov A.V. Intern. Conf. Infinite and Infinitesimal in Mathematics, Computing and Natural Sciences. Book of Abstracts. Cetraro. 2010.

16.Э.Хайрер, С.Нёрсетт, Г.Ваннер. Решение обыкновенных дифференциальных уравнений. Нежесткие задачи. М.Мир, 1990. - 512 с.

П.Годунов С.К., Рябенький B.C. Разностные схемы. М.: Наука. 1973. -440 с.

18.Самарский A.A. Теория разностных схем. М.: Наука. 1977. - 656 с.

19.Калиткин H.H. Численные методы. M.: Наука. 1978.

20.Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. М.: Наука, 1967. - 368 с.

21.Бахвалов Н.С. Численные методы. М.: Наука, 1975. - 631 с.

22.Бахвалов Н.С., Н.П.Жидков, Кобельков Г.М. Численные методы. М.: БИНОМ, 2008. - 630 с.

23.Вержбицкий В.М. Основы численных методов: Учебник для вузов/В.М. Вержбицкий. - 3-е изд., стер. - М.: Выел. Шк., 2009. - 840с.: ил.

24.Абрамович В. Числа Бернулли // Квант. — 1974. — № 6. — С. 10-14.

25.Гашков С.Б. Системы счисления и их применение — М.: МЦНМО, 2004. - 52 с.

26. Фомин C.B. Системы счисления — М.: Наука, 1987. —48 с.

27.Соболь И.М. Численные методы Монте-Карло. - М.: Наука, 1973. - 312 с.

28.Каханер Д., Моулер К., Нэш С. Численные методы и программное обеспечение. М.: Мир, 2001. - 575 с.

29.Кнут Д. Искусство программирования. Получисленные методы. - М.: ООО «И.Д. Вильяме», 2007. Т.2. - 832 с.

30.Lehmer, D. H. (1949). "Mathematical methods in large-scale computing units". Proceedings of a Second Symposium on Large-Scale Digital Calculating Machinery: 141-146. MR0044899. (journal version: Annals of the Computation Laboratory of Harvard University, Vol. 26 (1951)).

31.Бернштейн С.H. Теория вероятностей, M.: ГТТИ, 1946. - 367 с.

32.Скоробогатько В.Я. Теория ветвящихся цепных дробей и ее применение в вычислительной математике. - М.: Наука, 1983. - 312 с.

33.Султангазин У.М. Дискретные нелинейные модели уравнения Больцмана. - Алма-Ата: Наука, 1985. - 192 с.

34.Малинецкий Г.Г. Хаос. Структуры. Вычислительный эксперимент (Введение в нелинейную динамику). Изд.З, стереотипное 2002. - 256 с.

35.Aristov V.V. Direct Methods for Solving the Boltzmann Equation and Study of Nonequilibrium Flows. Dordrecht: Kluwer Academic Publishers, 2001.

36.Аристов B.B. Решение уравнения Больцмана при малых числах Кнудсена // Ж. вычисл. матем. и матем. физ. 2004, Т.44, N4. с.459-471.