Уравнение Риккати с производной дробного переменного порядка

Исследование дифференциальных уравнений дробных порядков. Наличие в процессе эффекта памяти или нелокальности по времени. Эредитарость в уравнение Риккати. Определения производной дробного переменного порядков. Интегро-дифференциальная задача Коши.

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

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

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

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

Уравнение Риккати с производной дробного переменного порядка

Твёрдый Д.А.

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

Дифференциальные уравнения дробных порядков представляют большой интерес для исследования, так как часто находят свое применение во многих областях, таких как: математика, физика и др. [1,2]. Уравнения с производными дробного порядка принадлежат классу интегро-дифференциальных уравнений и называются по терминологии В. Вольтерра эредитарными [3]. Понятие эредитарности означает наличие в рассматриваемом процессе эффекта памяти или нелокальности по времени. Нелокальность по времени содержится в ядре интегрального оператора исходного уравнения и называется функцией памяти. В случае, когда функция памяти является степенной, то мы естественным образом переходим к уравнениям с дробными производными, которые изучаются в рамках дробного исчисления [4].

В работе объектом нашего исследования будет уравнение Риккати [5] с учетом эредитарности. Эредитарость в уравнение Риккати характеризуется производной дробного переменного порядка. Отметим, что в работе [6] авторы исследовали эредитарное уравнение Риккати, когда порядок дробной производной является константой.

Постановка задачи и метод решения. Рассмотрим следующее эредитарное уравнение [1]:

, (1)

где - функция памяти, , - время моделирования, - функция решения. Уравнение (1) является аналогом классического уравнения Рикккати [5] и оно учитывает эффект памяти.

Если функция памяти представляет собой функцию Хевисайда, то мы можем говорить, что процесс обладает полной памятью, если это функция Дирака, то память отсутствует. Поэтому мы рассмотрим функцию памяти в виде степенной функции:

, , (2)

где - гамма-функция (функция Эйлера).

Процессы с функцией памяти вида (2) называются процессами с частичной потерей памяти и заслуживают особого внимания в их изучении. Это связано с тем, что многие естественные процессы имеют степенные законы распределения, которые в большинстве случаев приводят к понятию фрактальности или фракталу [7].

Подставив функцию памяти (2) в эредитарное уравнение (1) мы приходим к следующему интегро-дифференциальному уравнению Риккати:

. (3)

В уравнении (3) введем следующее обозначение:

, (4)

которое является обобщением дробного оператора Капуто или Герасимова-Капуто.

Необходимо отметить, что существуют другие определения производной дробного переменного порядков [8]. Мы же остановимся на определении (4), поэтому уравнение (3) запишем в компактной форме:

, (5)

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

, (6)

где - const.

Поэтому постановка задачи для эредитарного уравнения Риккати в нашем случае свелась к задаче Коши (5) и (6).

Заметим, что в случае, когда - const, мы приходим к задаче Коши, рассмотренной в работе [6], если , то задача Коши (5) и (6) переходит в классическую задачу Коши для уравнения Риккати [5].

Задача Коши (5) и (6) в общем случае не имеет точного решения, поэтому мы будем использовать численные методы для ее решения. Для этого разобьём временной отрезок на равных частей, где - шаг дискретизации, и получим что , а функция решения . Аппроксимацию дробной производной (4) проведем согласно работе [9] в виде:

, (7)

Где

, .

Можно показать, что аппроксимация (7) имеет первый порядок. Интегро-дифференциальную задачу Коши (5) и (6) можно переписать в разностной постановке:

, . (8)

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

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

Пример 1. Рассмотрим сначала случай, исследованный в работе [6], когда параметр является константой. Значения управляющих параметров выберем следующими: . Построим на одном графике несколько расчетных кривых, соответствующих различным значениям (рис. 1).

Рис. 1. Семейство расчетных кривых, полученных согласно решению системы (8) при различных значениях дробного параметра

На рис. 1. приведено семейство расчетных кривых задачи Коши (5) и (6) в зависимости от значений дробного параметра : (кривая 1 - соответствует классическому решению уравнения Риккати), (кривая 2), (кривая 3), (кривая 4), (кривая 5), (кривая 6).

Можно отметить, что наличие дробного параметра , при уменьшении его значений, в исходном уравнение приводит к перестройке расчетных кривых численных решений задачи Коши (5) и (6). Это связано с тем, что наличие памяти в рассматриваемом процессе приводит к так называемым "тяжелым затягивающимся хвостам" в кривых распределении, полученных решений, например, кривая 5 из рис. 1.

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

Рассмотрим другой пример, когда - является функцией, в том числе и случайной.

Пример 2. Пусть функция изменяется на интервале от 0 до 1 в соответствии с нормальным законом распределения, с единичным математическим ожиданием и нулевым среднеквадратическим отклонением. Значения управляющих параметров выберем следующими: .

Рис. 2. Результаты моделирования для примера 2: а) изменение значений параметра ; b) расчетные кривые, кривая 1 - классическое решение при ; кривая численного решения системы (8); c) фазовая траектория, построенная в координатах

На рис. 2. Приведены результаты моделирования. Введения случайной функции (рис. 2a) приводит к случайной функции решения (рис. 2b), что также отражено в хаотическом режиме на фазовой плоскости.

Пример 3. Рассмотрим пример, когда

,

где значения управляющих параметров выберем следующими: , .

Рис. 3. Результаты моделирования по примеру 3: a) расчетные кривые, кривая 1 - классическое решение уравнение Риккати; кривая 2 - решение системы (8); b) фазовая траектория

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

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

Возможное продолжение исследования эредитарного уравнения Риккати связано с прикладными задачами, например в экономике [10], а также в решении обратной задачи оценки параметра по известным экспериментальным данным.

Выражаю благодарность своему научному руководителю к.ф.-м.н., профессору РАЕ Паровику Роману Ивановичу за ценные советы при планировании исследования и рекомендации по оформлению статьи.

Литература

1. Учайкин В.В. Метод дробных производных. Ульяновск: Артишок, 2008. 512 с.

2. Petras I. Fractional-Order Nonlinear Systems. Modeling, Analysis and Simulation. Beijing and Springer-Verlag Berlin Heidelberg: Springer, 2011. 218 p.

3. Volterra V. Sur les 'equations int'egro-diff'erentielles et leurs applications // Acta Mathematica. 1912. Vol. 35, no. 1. P. 295-356.

4. Нахушев А.М. Дробное исчисление и его приложения. М.: Физматлит, 2003. 272 с.

5. Riccati J. F. Animadversationes in aequationes differentiales secundi gradus. Acta Eruditorum Quae Lipside Publicantur, 1724. Supplementa 8.

6. Sweilam N. H., Khader M. M., Mahdy A. M. S. Numerical studies for solving fractional Riccati differential equation // Applications and Applied Mathematics. 2012. Т. 7. №. 2. С. 595-608.

7. Федер Е. Фракталы. М.: Мир, 1991. 254 с.

8. Паровик Р.И. Математическое моделирование линейных эредитарных осцилляторов: Петропавловск-Камчатский: КамГУ им. Витуса Беринга 2015. 178 с.

9. Паровик Р.И. Конечно-разностные схемы для фрактального осциллятора с переменными дробными порядками // Вестник КРАУНЦ. Физико-математические науки. 2015. №. 2 (11). С. 88-95.

10. Makarov D.V., Parovik R.I. Modeling of the economic cycles using the theory of fractional calculus // Journal of Internet Banking and Commerce. 2016. vol. 21. no S6.

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


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

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

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

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

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

  • Степенные ряды. Радиус сходимости. Ряды Лорана. Полюса и особые точки. Интегрирование дифференциальных уравнений при помощи степенных рядов. Общее дифференциальное уравнение Риккати. Исследование решений в окрестности полюса и существенно особой точки.

    дипломная работа [252,1 K], добавлен 15.12.2012

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

    курсовая работа [209,4 K], добавлен 04.01.2016

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

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

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

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

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

    презентация [246,0 K], добавлен 21.09.2013

  • Дифференциальные уравнения Риккати. Общее решение линейного уравнения. Нахождение всех возможных решений дифференциального уравнения Бернулли. Решение уравнений с разделяющимися переменными. Общее и особое решения дифференциального уравнения Клеро.

    курсовая работа [347,1 K], добавлен 26.01.2015

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

    презентация [77,7 K], добавлен 17.09.2013

  • Типы уравнений, допускающих понижение порядка. Линейное дифференциальное уравнение высшего порядка. Теоремы о свойствах частичных решений. Определитель Вронского и его применение. Использование формулы Эйлера. Нахождение корней алгебраического уравнения.

    презентация [103,1 K], добавлен 29.03.2016

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