Введение. Эредитарные процессы представляют особый интерес в связи с различными приложениями. Например, в монографии В.В. Учайкина [1] эредитарным процессам посвящена целая глава, там же приведен пример эредитарного осциллятора, который был впервые изучен итальянским математиком В. Вольтерра, а результаты исследований были приведены в его работе [2]. Эффект последействия или эредитарности характеризует зависимость текущего состояния системы от ее предыдущих состояний, необязательно от всех. Математическое описание эффекта памяти дается интегро-дифференциальным уравнением, при чем ядро этого уравнения называется функции памяти. В случае, когда функция памяти имеет степенной вид интегро-дифференциальное уравнение переходит в уравнение с производными дробных порядков, которые изучаются в рамках теории дробного исчисления [3].
Некоторые эредитарные динамические системы были описаны в монографиях [4,5]. В работе рассмотрена модель осциллятора Ван дер Поля-Дуффинга с учетом эредитарности, которая была решена с помощью конечно-разностной схемы. С помощью численного решения в зависимости от различных значений управляющих параметров были построены осциллограммы и фазовые траектории.
Постановка задачи и методика решения. Модель Ван-дер-Поля – Дуффинга (ВПД) в присутствии периодического внешнего воздействия описывается уравнением:
(1)
где 
– управляющий параметр, 
– частота, а 
- амплитуда внешнего воздействия, 
– параметр фазовой нелинейности, определяющий не изохронность колебаний. Заметим, в уравнении (1) присутствуют нелинейности характерные для осцилляторов Ван дер Поля [6] и Дуффинга [7]. Квадратичная нелинейность в коэффициенте при младшей производной, характеризует автоколебания, а кубическая нелинейность - зависимость периода колебаний от амплитуды.
Для уравнения (1) введем начальные условия:
(2)
Задача (1), (2) является классической дифференциальной задачей Коши для уравнения ВПД. Введем эредитароность в уравнение (1). Для этого рассмотрим следующее интегро-дифференциальное уравнение ВПД:
(3)
где 
и 
– функции памяти, характеризующие эредитарность.
Отметим, что если функции памяти представляют собой дельта-функции, тогда в системе отсутствует эредитарность. Если функции памяти представляют собой функции Хэвисайда тогда система обладает полной памятью. Если функции памяти являются степенными функциями, например,
(4)
где Г(t) - гамма-функция Эйлера, тогда говорят, что система - с частичной "потерей памяти".
Подставим функции памяти (4) в интегро-дифференциальное уравнение (3). В результате получим:
. (5)
Мы получили интегро-дифференциальное уравнение специального вида. Если обратиться к определению производной дробного порядка по Капуто, то мы приходим к уравнению с дробными производными следующего вида:
, (6)
где дробные дифференциальные операторы равны
определенные в смысле Герасимова-Капуто с дробными порядками 
.
Можно отметить, что в предельном случае (
) уравнение (6) переходит в классическое уравнение ВПД (1), поэтому можно считать уравнение (1) частным случаем обобщения уравнения (6). Интегро-дифференциальное уравнение ВПД (6) будем называть дробным, или фрактальным уравнением, а процесс, которые оно описывает, будем называть фрактальными, или эредитарными.
Отметим, что функции памяти в интегро-дифференциальном уравнении (5) могут быть отличными от степенных функций, что приводит к другим интегро-дифференциальным уравнениям.
Задача Коши (2), (6) в общем виде не имеет точного решения в силу того, что модельное уравнение является нелинейным, поэтому надо использовать численные методы для ее решения. В качестве численного метода возьмем метод конечно-разностных схем, так как его легко можно реализовать в любой компьютерной среде.
Будем рассматривать равномерную сетку. Для этого разобьем временной интервал 
на 
равных частей. В результате получим равномерную сетку 
, где шаг сетки 
, 
. Значения искомой функции 
, будем называть ее сеточной функцией. Аппроксимация дробных операторов задачи Коши (6) осуществляется следующим образом [5]:
,
,
, 
.
Подставим эти аппроксимации в модельное уравнение (6). Приходим к следующей конечно-разностной задаче:
(7)
Отметим, что дифференциальная задача Коши (2), (6) является жесткой при больших значений управляющего параметра , поэтому явная конечно-разностная схема (7) будет работать в случае уменьшение шага сетки 
. Оценка шага - это отдельная задача, фактически связанная с качественными свойствами схемы (7) – устойчивостью и сходимостью и в настоящей работе не рассматривалась. Будем считать, что значения параметра достаточно малы для того, чтобы в процессе вычислений не уменьшать шаг сетки , т.е. будем рассматривать лишь «мягкую» задачу Коши (2) и (6).
Результаты моделирования и их обсуждение. Рассмотрим различные примеры работы конечно-разностной схемы (8). Построим осциллограммы и фазовые траектории в зависимости от различных значений управляющих параметров нелинейной эредитарной колебательной системы ВДП.
Пример 1. (Изменение параметра 
). Значения управляющих параметров выберем следующими: 
. На рис. 1 приведены осциллограммы для примера 1, при различных значениях дробного параметра : 
Рис. 1. Осциллограммы, построенные по схеме (8) для примера 1 при различных значениях
Из рис. 1 видно, что при уменьшении значения дробного параметра , амплитуда колебаний быстрее выходит на постоянный уровень. Этот факт также подтверждают фазовые траектории (рис. 2), которые выходят на соответствующие предельные циклы.
Рис. 2. Фазовые траектории, полученные по схеме (8) для примера 1 при различных значениях
Пример 2. (Изменение параметра 
). Значения управляющих параметров выберем следующими: 
. На рис. 2 приведены осциллограммы для примера 2, при различных значениях дробного параметра : 
Рис. 3. Осциллограммы, построенные по схеме (8) для примера 2 при различных значениях
Из рис. 3 мы можем заметить, что уменьшение значений дробного параметра существенного изменения в форму осциллограмм не вносит. Колебания со временем выходят на установившийся режим с постоянной амплитудой. Поэтому можно сделать вывод, что фазовые траектории выходят практически на один и тот же предельный цикл (рис. 4).
Рис. 4. Фазовые траектории, полученные по схеме (8) для примера 2 при различных значениях
Пример 3. (Изменение параметра ). Значения управляющих параметров выберем следующими: 
.
На рис. 5 приведены осциллограммы для примера 3, при различных значениях амплитуды внешнего воздействия : 
Рис. 5. Осциллограммы, построенные по схеме (8) для примера 5 при различных значениях
Из рис. 5 видно, что отсутствие внешнего воздействия (рис.5a) приводит к росту амплитуды и начиная с некоторого момента времени ее значения выходят на постоянный уровень. Далее мы видим сложные по форме колебания, которые по-видимому, говорят о возможности много периодических решений задачи Коши (2) и (6). Это косвенно подтверждают соответствующие фазовые траектории (рис. 6).
Рис. 4. Фазовые траектории, полученные по схеме (8) для примера 5 при различных значениях
В заключении хочется отметить, что наличие различных колебательных режимов эредитарного осциллятора ВПД, требует дальнейшего его изучение. Например, интерес представляет построение карт динамических режимов и сечений Пуанкаере с целью классификации периодических решений, а также устойчивости точки покоя [8]. Другое направление исследований заключается в изучение качественных свойств конечно-разностной схемы (8) [9], а также жесткости задачи Коши (2), (6). Так же справедливо дальнейшее обобщение задачи Коши (2), (6) в случае функций 
и 
по аналогии с работами [9,10].
Автор выражает благодарность научному руководителю, к.ф.-м.н., профессору РАЕ, Р.И. Паровику за ценные советы и замечания по содержанию данной научной статьи.