Метод Паркера – Сохацкого - Parker–Sochacki method

В математика, то Метод Паркера – Сохацкого является алгоритм для решения систем обычных дифференциальные уравнения (ODE), разработанная Г. Эдгар Паркер и Джеймс Сохацки, из Университет Джеймса Мэдисона Математический факультет. Метод производит Серия Маклорена решения систем дифференциальных уравнений с коэффициентами в алгебраической или числовой форме.

Резюме

Метод Паркера – Сохацки основан на двух простых наблюдениях:

  • Если набор ODE имеет определенную форму, то Метод Пикара можно использовать для поиска их решения в виде степенной ряд.
  • Если ОДУ не имеют требуемой формы, почти всегда можно найти расширенный набор уравнений, который действительно имеет требуемую форму, так что подмножество решения является решением исходных ОДУ.

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

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

Преимущества

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

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

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

Недостатки

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

использованная литература

  1. ^ П.Г. Warne; Д.П. Warne; J.S. Сохацкий; G.E. Паркер; Д. К. Карозерс (2006). «Явные априорные границы ошибок и адаптивный контроль ошибок для аппроксимации нелинейных дифференциальных систем начального значения» (PDF). Компьютеры и математика с приложениями. 52 (12): 1695–1710. Дои:10.1016 / j.camwa.2005.12.004. Получено 27 августа, 2017.

внешние ссылки

  • Полиномиальные ОДУ - примеры, решения, свойства (PDF), получено 27 августа, 2017. Подробное объяснение парадигмы и применения метода Паркера – Сохацкого.
  • Джозеф В. Рудмин (1998), "Применение метода Паркера – Сохацки к небесной механике", Журнал вычислительной неврологии, 27: 115–133, arXiv:1007.1677, Дои:10.1007 / s10827-008-0131-5. Демонстрация теории и использования метода Паркера – Сохацкого, включая решение классического ньютоновского Nпроблема тела с взаимным гравитационным притяжением.
  • Модифицированный метод Пикара., получено 11 ноября, 2013. Сборник статей и немного кода Matlab.