Дробное исчисление

Дробное исчисление (англ. fractional calculus) — раздел математического анализа, изучающий возможность применения операций дифференцирования и интегрирования к функциям с нецелым (дробным, иррациональным, комплексным) порядком. В отличие от классического анализа, где порядок производной или интеграла является натуральным числом (n), дробное исчисление обобщает эти понятия на произвольный действительный или комплексный показатель α. Это позволяет описывать процессы с памятью, наследственностью и нелокальными свойствами, которые не могут быть адекватно смоделированы в рамках целочисленного анализа.

История

Идея обобщения порядка производной возникла практически одновременно с созданием классического анализа. В переписке Готфрида Вильгельма Лейбница с Гийомом де Лопиталем (1695 год) Лейбниц упомянул возможность производной порядка 1/2, назвав это «парадоксом, из которого однажды извлекут полезные следствия». Однако строгая математическая теория начала формироваться лишь в XIX веке.

Ранние работы

• Нильс Хенрик Абель (1823) — при решении задачи о таутохроне (кривой, по которой тело скатывается за одинаковое время независимо от начальной точки) неявно использовал интегрирование дробного порядка.
• Жозеф Лиувилль (1832–1837) — опубликовал первые систематические работы по дробному исчислению, введя определение производной через ряды Фурье и экспоненциальные функции. Его определение, однако, было ограничено функциями, разложимыми в ряд.
• Бернхард Риман (1847, опубликовано посмертно в 1876) — разработал альтернативное определение, основанное на обобщении формулы Коши для кратного интеграла. Вместе с работой Лиувилля оно легло в основу наиболее распространённого определения — дробного интеграла Римана — Лиувилля.

Развитие в XX веке

В начале XX века теория была формализована в работах Адамара, Маршо, Летникова и Вейля. В 1960–1970-х годах интерес к дробному исчислению резко возрос в связи с развитием теории вязкоупругости, теории управления и физики сложных систем. В 1974 году в Нью-Хейвене (США) прошла первая международная конференция по дробному исчислению. В 1990-х годах появились работы, связывающие дробные производные с фрактальной геометрией и аномальной диффузией.

Основные определения

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

Дробный интеграл Римана — Лиувилля

Наиболее фундаментальное понятие — дробный интеграл порядка α > 0, который обобщает формулу Коши для n-кратного интеграла:

\[ I_a^\alpha f(x) = \frac{1}{\Gamma(\alpha)} \int_a^x (x-t)^{\alpha-1} f(t) \, dt \]

где \(\Gamma(\alpha)\) — гамма-функция Эйлера, обобщающая факториал на нецелые числа. При α = 1 формула сводится к обычному интегралу от a до x. При α = 0 — к тождественному оператору.

Дробная производная Римана — Лиувилля

Определяется как производная целого порядка от дробного интеграла:

\[ D_a^\alpha f(x) = \frac{d^n}{dx^n} I_a^{n-\alpha} f(x) \]

где n = ⌈α⌉ — наименьшее целое число, большее или равное α. Например, для α = 0.5, n = 1, и производная порядка 0.5 вычисляется как производная первого порядка от интеграла порядка 0.5.

Дробная производная Капуто

Введена Микеле Капуто в 1967 году для решения задач с начальными условиями, имеющими физический смысл. Отличается порядком операций: сначала берётся производная целого порядка, затем — дробный интеграл:

\[ ^C D_a^\alpha f(x) = I_a^{n-\alpha} \left( \frac{d^n}{dx^n} f(x) \right) \]

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

Другие определения

• Производная Грюнвальда — Летникова — основана на обобщении разностной формулы для производной на дробный порядок через биномиальные коэффициенты.
• Производная Римана — Лиувилля с переменным порядком — порядок α может зависеть от x или t.
• Дробная производная в смысле Вейля — для периодических функций, определяется через коэффициенты Фурье.
• Производная Хильфера — комбинация производных Римана — Лиувилля и Капуто.

Свойства

Дробные операторы обладают рядом свойств, аналогичных целочисленным, но с важными отличиями.

• Линейность: \(D^\alpha (af + bg) = a D^\alpha f + b D^\alpha g\).
• Полугрупповое свойство: для дробных интегралов \(I^\alpha I^\beta = I^{\alpha+\beta}\) выполняется. Для дробных производных, вообще говоря, не выполняется: \(D^\alpha D^\beta f \neq D^{\alpha+\beta} f\).
• Коммутативность: дробные интегралы коммутируют: \(I^\alpha I^\beta = I^\beta I^\alpha\). Дробные производные не коммутируют.
• Правило Лейбница: обобщение для дробного порядка имеет бесконечный ряд и не сводится к простой формуле.
• Память и нелокальность: дробная производная в точке x зависит от значений функции на всём интервале от a до x (или от -∞ до x), что отличает её от локальной целочисленной производной.

Применение

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

Физика

• Вязкоупругость: модели Зенера, Кельвина — Фойгта и Максвелла с дробными производными лучше описывают поведение полимеров и биологических тканей, чем классические модели. Дробный порядок α (0 < α < 1) интерпретируется как показатель «памяти» материала.
• Аномальная диффузия: уравнение диффузии с дробной производной по времени (субдиффузия, α < 1) или по пространству (супердиффузия, β > 2) описывает движение частиц в пористых средах, мембранах и плазме.
• Электродинамика: дробные модели диэлектрической релаксации (модель Коула — Коула) и распространения электромагнитных волн в фрактальных средах.
• Квантовая механика: дробное уравнение Шрёдингера (введено Ником Ласкиным в 2000 году) обобщает кинетическую энергию через дробный оператор Лапласа.

Техника и управление

• Теория автоматического управления: дробные ПИД-регуляторы (PI^λ D^μ) позволяют добиться большей гибкости и робастности по сравнению с классическими ПИД-регуляторами. Параметры λ и μ — дробные порядки интегрирования и дифференцирования.
• Обработка сигналов: дробное преобразование Фурье и дробная фильтрация используются для анализа нестационарных сигналов, в том числе в радиолокации и телекоммуникациях.

Биология и медицина

• Моделирование роста опухолей: дробные модели учитывают гетерогенность ткани и эффекты памяти.
• Электрофизиология: модели потенциалов действия нейронов (например, модель Ходжкина — Хаксли с дробными производными) лучше описывают адаптацию и утомление.
• Фармакокинетика: дробные модели описывают неэкспоненциальное выведение лекарств из организма.

Экономика и финансы

• Фрактальные модели рынков: дробное броуновское движение (с показателем Хёрста H ≠ 0.5) используется для моделирования временных рядов с долговременной памятью (например, цен акций).
• Опционы: дробное уравнение Блэка — Шоулза обобщает классическую модель на случай нецелого порядка.

Критика и ограничения

• Неединственность определений: отсутствие единого, общепринятого определения дробной производной затрудняет сравнение результатов разных авторов. Выбор определения часто диктуется удобством, а не физической обоснованностью.
• Интерпретация порядка: физический смысл дробного порядка α (например, 0.7) остаётся предметом дискуссий. В ряде случаев α может быть связан с фрактальной размерностью системы, но прямая связь установлена не для всех моделей.
• Вычислительная сложность: численные методы для дробных уравнений (например, метод Грюнвальда — Летникова) требуют хранения всей истории процесса, что приводит к большим затратам памяти и времени счёта. Для решения этой проблемы разрабатываются методы с «усечённой памятью» и аппроксимации рациональными функциями.
• Отсутствие физических законов: в отличие от целочисленных производных, которые непосредственно связаны с законами Ньютона и сохранения, дробные производные редко следуют из фундаментальных принципов. Чаще они являются феноменологическими подгоночными моделями.

Интересные факты

• Термин «дробное исчисление» (fractional calculus) ввёл в обиход Оливер Хевисайд в конце XIX века, хотя он использовал операторный метод, не давая строгого математического обоснования.
• Символ дробной производной часто записывают как \(_a D_x^\alpha f(x)\), где a — нижний предел интегрирования. При a = 0 (для функций, определённых на полуоси) часто используют обозначение \(_0 D_x^\alpha\).
• Гамма-функция \(\Gamma(z)\) является ключевым элементом дробного исчисления, так как она обобщает факториал на нецелые числа. При отрицательных целых z она имеет полюса, что ограничивает область применения некоторых формул.
• Существует понятие дробной производной по времени с переменным порядком, где α = α(t). Такие модели используются для описания старения материалов или изменения свойств среды во времени.

Источники

• Самко С. Г., Килбас А. А., Маричев О. И. Интегралы и производные дробного порядка и некоторые их приложения. — Минск: Наука и техника, 1987.
• Podlubny I. Fractional Differential Equations. — San Diego: Academic Press, 1999.
• Kilbas A. A., Srivastava H. M., Trujillo J. J. Theory and Applications of Fractional Differential Equations. — Amsterdam: Elsevier, 2006.
• Учайкин В. В. Метод дробных производных. — Ульяновск: Артишок, 2008.
• Mainardi F. Fractional Calculus and Waves in Linear Viscoelasticity. — London: Imperial College Press, 2010.