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

Шаблон:Числення Дробове числення — розділ математичного аналізу, що вивчає різні способи задання операторів диференціювання $D$ $$Df(x)=\frac{d}{dx}f(x),$$ і інтегрування $J$ $$Jf(x)=\underset{0}{\overset{x}{\int }}f(s)ds$$ дійсного або комплексного порядку.

Історія

У прикладній математиці та математичному аналізі дробова похідна — це похідна будь-якого довільного порядку, дійсного чи комплексного. Вперше про неї згадав 1695 року Готфрід Вільгельм Лейбніц у листі, до Гійома де Лопіталя. Приблизно в той самий час Лейбніц написав одному з братів Бернулі, описуючи подібність між біноміальною теоремою та правилом Лейбніца для дробової похідної добутку двох функцій.

Дробове числення введено в одній з ранніх праць Нільса Генріка Абеля, в якій можна побачити багато його елементів: ідею дробового інтегрування та дробового диференціювання, взаємно обернений зв'язок між ними, розуміння того, що дробові диференціювання та інтегрування можна розглядати як одну й ту саму узагальнену операцію, і навіть уніфіковану нотацію для диференціювання та інтегрування довільного дійсного порядку.

Незалежно від нього, Ліувілль заклав основи предмету в статті 1832 року. Близько 1890 року самоук Олівер Гевісайд представив практичне застосування дробових диференціальних операторів до аналізу ліній електропередач. Теорія та застосування дробового числення значно розширилися протягом XIX та XX століть. Численні автори давали різні визначення дробових похідних та інтегралів.

Дробові інтеграли

Нехай $f$ — функція, визначена на $(0,+\mathrm{\infty })$. Якщо оператор $$(Jf)(x)=\underset{0}{\overset{x}{\int }}f(t)dt$$ взяти двічі від $f$, то буде $$\begin{array}{cc}\left(J^{2}f\right)(x)& =\underset{0}{\overset{x}{\int }}(Jf)(t)dt\\ & =\underset{0}{\overset{x}{\int }}(\underset{0}{\overset{t}{\int }}f(s)ds)dt.\end{array}$$

І це можна повторювати довільну кількість разів. За Шаблон:Нп $$\left(J^{n}f\right)(x)=\frac{1}{(n-1)!}\underset{0}{\overset{x}{\int }}{(x-t)}^{n-1}f(t)dt,$$ де Шаблон:Math — будь-яке натуральне число.

Використання гамма-функції замість факторіала дає такий оператор дробового інтегрування:

$$\left(J^{\alpha }f\right)(x)=\frac{1}{\mathrm{\Gamma }(\alpha )}\underset{0}{\overset{x}{\int }}{(x-t)}^{\alpha -1}f(t)dt.$$

Отриманий у такий спосіб оператор Шаблон:Math задовольняє таку умову: $$\begin{array}{cc}\left(J^{\alpha }\right)\left(J^{\beta }f\right)(x)& =\left(J^{\beta }\right)\left(J^{\alpha }f\right)(x)\\ & =\left(J^{\alpha +\beta }f\right)(x)\\ & =\frac{1}{\mathrm{\Gamma }(\alpha +\beta )}\underset{0}{\overset{x}{\int }}{(x-t)}^{\alpha +\beta -1}f(t)dt.\end{array}$$

Це відношення називають напівгруповою властивістю дробових диферінтегральних операторів.

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

Класичною формою дробового числення є Шаблон:Нп, який, по суті, є тим, що описано вище. Теорію дробового інтегрування для періодичних функцій (включаючи «граничну умову» повторення через період) дає Шаблон:Нп. Він визначений на рядах Фур'є і вимагає, щоб вільний коефіцієнт тригонометричного ряду дорівнював нулю. Інтеграл Рімана — Ліувілля існує у двох формах, верхній та нижній. На відрізку Шаблон:Closed-closed ці форми визначають як $$\begin{array}{cc}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)& =\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)\\ & =\frac{1}{\mathrm{\Gamma }(\alpha )}\underset{a}{\overset{t}{\int }}{(t-\tau )}^{\alpha -1}f(\tau )d\tau ,\\ \text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)& =\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)\\ & =\frac{1}{\mathrm{\Gamma }(\alpha )}\underset{t}{\overset{b}{\int }}{(\tau -t)}^{\alpha -1}f(\tau )d\tau .\end{array}$$ Перша форма справедлива для Шаблон:Math, а друга — для Шаблон:Math.

Інтеграл на додатній дійсній півосі (тобто, Шаблон:Math), виходячи з історії відкриття та використання, запропоновано назвати інтегралом Абеля — Рімана, і, в тому ж ключі, інтеграл за всією дійсною прямою названо інтегралом Ліувілля — Вейля.

Дробовий інтеграл Адамара

Дробовий інтеграл Адамара, який увів Жак Адамар, задають такою формулою: $$Error\; parsing\; sideset\; expression,\; no\; succeeding\; operator\; found𝐃f(t)=\frac{1}{\mathrm{\Gamma }(\alpha )}\underset{a}{\overset{t}{\int }}{(\log \frac{t}{\tau })}^{\alpha -1}f(\tau )\frac{d\tau }{\tau },t>a.$$

Дробовий інтеграл Атангани — Балеану

Дробовий інтеграл Атангани — Балеану для неперервної функції визначають так: $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{1-\alpha }{\mathrm{AB}(\alpha )}f(t)+\frac{\alpha }{\mathrm{AB}(\alpha )\mathrm{\Gamma }(\alpha )}\underset{a}{\overset{t}{\int }}{(t-\tau )}^{\alpha -1}f(\tau )d\tau .$$

Дробові похідні

Аналогічний процес для оператора диференціювання Шаблон:Mvar є складнішим. Можна показати, що в загальному випадку Шаблон:Mvar не є ані комутативним, ані адитивним.

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

Дробові похідні функції Гаусса — неперервна інтерполяція між функцією та її першою похідною.

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

Дробову похідну Рімана — Ліувілля обчислюють за правилом Лагранжа для диференціальних операторів. Для знаходження похідної Шаблон:Mvar-го порядку обчислюють похідну Шаблон:Mvar-го порядку від інтеграла порядку Шаблон:Math, де Шаблон:Mvar — найменше ціле число, більше за Шаблон:Mvar (тобто, Шаблон:Math). Дробові похідна та інтеграл Рімана — Ліувілля мають низку застосувань. Подібно до визначення інтеграла Рімана — Ліувілля, похідна має верхню та нижню форми: $$\begin{array}{cc}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)& =\frac{d^n}{d{t}^n}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)\\ & =\frac{d^n}{d{t}^n}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t),\\ \text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)& =\frac{d^n}{d{t}^n}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)\\ & =\frac{d^n}{d{t}^n}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t).\end{array}$$

Дробова похідна Капуто

Іншим способом обчислення дробових похідних є дробова похідна Капуто. Її ввів Мікеле Капуто у своїй статті 1967 року. На відміну від дробової похідної Рімана-Ліувілля, при розв'язуванні диференціальних рівнянь використовуючи означення Капуто не потрібно визначати початкові умови дробового порядку. Означення Капуто вводится так (тут знову Шаблон:Math): $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{1}{\mathrm{\Gamma }(n-\alpha )}\underset{0}{\overset{t}{\int }}\frac{f^{(n)}(\tau )}{{(t-\tau )}^{\alpha +1-n}}d\tau .$$

Для $\nu \in (n-1,n)$ дробова похідна Капуто має такий вигляд: $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{1}{\mathrm{\Gamma }(n-\nu )}\underset{0}{\overset{t}{\int }}(t-u{)}^{(n-\nu -1)}{f}^{(n)}(u)du,$$ і має ту перевагу, що дорівнює нулю, коли Шаблон:Mvar є константою, а її перетворення Лапласа виражається через початкові значення функції та її похідної. Крім того, похідну Капуто для Шаблон:Closed-closed визначають як $$\begin{array}{cc}\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)& =\underset{a}{\overset{b}{\int }}\varphi (\nu )[\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)]d\nu \\ & =\underset{a}{\overset{b}{\int }}[\frac{\varphi (\nu )}{\mathrm{\Gamma }(1-\nu )}\underset{0}{\overset{t}{\int }}{(t-u)}^{-\nu }{f}^{\prime }(u)du]d\nu ,\end{array}$$ де Шаблон:Mvar — вагова функція.

Дробова похідна Капуто-Фабріціо

У статті 2015 року М. Капуто та М. Фабріціо представили означення дробової похідної з несингулярним ядром для неперервно-диференційованої функції Шаблон:Mvar, заданої так: $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{1}{1-\alpha }\underset{a}{\overset{t}{\int }}{f}^{\prime }(\tau )e^{(-\alpha \frac{t-\tau }{1-\alpha })}d\tau ,$$ де Шаблон:Nowrap

Дробова похідна Атангани — Балеану

У 2016 році Атангана та Балеану запропонували диференціальні оператори на основі узагальненої функції Міттага-Лефлера Шаблон:Math. Метою було ввести дробові диференціальні оператори з несингулярним нелокальним ядром. Їхні дробові диференціальні оператори наведено нижче в сенсі Рімана — Ліувілля та Капуто відповідно для неперервно-диференційованої функції Шаблон:Math: $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{\mathrm{AB}(\alpha )}{1-\alpha }\underset{a}{\overset{t}{\int }}{f}^{\prime }(\tau )E_{\alpha }(-\alpha \frac{(t-\tau {)}^{\alpha }}{1-\alpha })d\tau .$$ Якщо функція Шаблон:Math неперервна, то похідна Атангани — Балеану в сенсі Рімана — Ліувілля має вигляд $$\text{Not implemented MEDIAWIKIEXTENSIONMATHWIKITEXVCNODESTEXNODE for sideset}f(t)=\frac{\mathrm{AB}(\alpha )}{1-\alpha }\frac{d}{dt}\underset{a}{\overset{t}{\int }}f(\tau )E_{\alpha }(-\alpha \frac{(t-\tau {)}^{\alpha }}{1-\alpha })d\tau .$$

Ядро, що використовується в дробовій похідній Атангани — Балеану, має деякі властивості кумулятивної функції розподілу. Наприклад, для всіх $\alpha \in (0,1]$ функція Шаблон:Math зростає на дійсній прямій, збігається до Шаблон:Math в Шаблон:Math, і Шаблон:Nowrap Отже, функція $x\mapsto 1-E_{\alpha }(-x^{\alpha })$ є кумулятивною функцією розподілу ймовірнісної міри на додатних дійсних числах. Таким чином, визначено розподіл, і будь-який його кратний розподіл називається Шаблон:Нп порядку Шаблон:Math. Також, усі ці розподіли ймовірностей є абсолютно неперервними. Зокрема, функція Міттага-Леффлера має окремий випадок Шаблон:Math, коли є експонентою. Таким чином, розподіл Міттага-Леффлера порядку Шаблон:Math є експоненційним розподілом.

Дробова похідна Ріса

Похідну Ріса визначають як $$ℱ\left\{\frac{{\partial }^{\alpha }u}{\partial {\left|x\right|}^{\alpha }}\right\}(k)=-{\left|k\right|}^{\alpha }ℱ\{u\}(k),$$ де $ℱ$ позначає перетворення Фур'є.

Інші типи

До класичних дробових похідних належать:

• Шаблон:Нп
• Похідна Соніна — Летнікова
• Похідна Ліувілля
• Похідна Капуто
• Похідна Адамара
• Похідна Маршо
• Похідна Ріса
• Похідна Міллера — Росса
• Шаблон:Нп
• Шаблон:Нп
• $F^{\alpha }$-похідна

До нових дробових похідних належать:

• Похідна Коїмбри
• Похідна Катугампола
• Похідна Гільфера
• Похідна Девідсона
• Похідна Чена
• Похідна Капуто — Фабріціо
• Похідна Атангани — Балеану

Узагальнення

Оператор Ерделі — Кобера

Оператор Ерделі — Кобера — це інтегральний оператор, який 1940 року ввели Шаблон:Нп та Шаблон:Нп, має вигляд $$\frac{x^{-\nu -\alpha +1}}{\mathrm{\Gamma }(\alpha )}\underset{0}{\overset{x}{\int }}{(t-x)}^{\alpha -1}{t}^{-\alpha -\nu }f(t)dt,$$ який узагальнює дробовий інтеграл Рімана — Ліувілля та інтеграл Вейля.

Застосування

Дробове збереження маси

Рівняння дробового збереження маси необхідне для моделювання потоку рідини, коли Шаблон:Нп недостатньо великий порівняно з Шаблон:Нп і коли потік у контрольному об'ємі є нелінійним: $$-\rho ({\mathrm{\nabla }}^{\alpha }\cdot \overrightarrow{u})=\mathrm{\Gamma }(\alpha +1)\mathrm{\Delta }{x}^{1-\alpha }\rho ({\beta }_s+\varphi {\beta }_w)\frac{\partial p}{\partial t}.$$

Електрохімічний аналіз

При вивченні окисно-відновлювальної поведінки субстрату в розчині до поверхні електрода прикладають напругу, щоб змусити електрони переходити між електродом і субстратом. Перенос електронів, що виникає в результаті, вимірюється як струм. Струм залежить від концентрації субстрату на поверхні електрода. Коли підкладка витрачається, свіжа підкладка дифундує до електрода, як описано в законах дифузії Фіка. Перетворення Лапласа другого закону Фіка дає звичайне диференціальне рівняння другого порядку (в безрозмірній формі): $$\frac{d^2}{d{x}^2}C(x,s)=sC(x,s).$$ Якщо взяти похідну від Шаблон:Math, а потім обернене перетворення Лапласа, то отримаємо таку залежність: $$\frac{d}{dx}C(x,t)=\frac{d^{\frac{1}{2}}}{d{t}^{\frac{1}{2}}}C(x,t),$$ яка пов'язує концентрацію субстрату на поверхні електрода зі струмом. Ця залежність застосовується в електрохімічній кінетиці для з'ясування механістичної поведінки. Наприклад, її використано для вивчення швидкості димеризації субстратів при електрохімічному відновленні.

Задача потоку підземних вод

У 2013—2014 роках описано деякі задачі потоку підземних вод, використовуючи поняття дробової похідної. Класичний закон Дарсі узагальнено, розглядаючи потік води як функцію похідної нецілого порядку від п'єзометричного напору. Цей узагальнений закон і закон збереження маси використали для виведення нового рівняння для потоку підземних вод.

Моделі просторово-часових дробових рівнянь дифузії

Аномальні дифузійні процеси в складних середовищах можна добре описати за допомогою моделей рівнянь дифузії дробового порядку. Часова похідна відповідає довготривалому розпаду важкого хвоста, а просторова похідна — нелокальності дифузії. Рівняння просторово-часової дробової дифузії можна записати у вигляді $$\frac{{\partial }^{\alpha }u}{\partial t^{\alpha }}=-K(-\mathrm{\Delta }{)}^{\beta }u.$$

Простим продовженням дробової похідної є дробова похідна змінного порядку, за якого Шаблон:Mvar і Шаблон:Mvar змінюються на Шаблон:Math і Шаблон:Math. Його можна застосовувати в моделюванні аномальної дифузії.

Моделі структурного згасного коливання

Дробові похідні використовують для моделювання в'язкоеластичного згасного коливання в певних типах матеріалів, таких як полімери.

ПІД-регулятори

Узагальнення ПІД-регуляторів для використання дробових порядків може збільшити ступінь їхньої свободи. Нове рівняння, що зв'язує керувальну змінну Шаблон:Math з виміряним значенням похибки Шаблон:Math, можна записати як $$u(t)=K_{\mathrm{p}}e(t)+K_{\mathrm{i}}{D}_t^{-\alpha }e(t)+K_{\mathrm{d}}{D}_t^{\beta }e(t),$$ де Шаблон:Mvar і Шаблон:Math — додатні дробові порядки, а Шаблон:Math, Шаблон:Math, і Шаблон:Math — невід'ємні коефіцієнти при пропорційному, інтегральному і похідному членах відповідно (іноді позначається, як Шаблон:Mvar, Шаблон:Mvar, і Шаблон:Mvar).

Рівняння акустичних хвиль для складних середовищ

Поширення акустичних хвиль у складних середовищах, таких як біологічні тканини, зазвичай передбачає згасання, що підпорядковується частотному степеневому закону. Таке явище можна описати за допомогою причинно-наслідкового хвильового рівняння, яке включає дробові похідні за часом: $${\mathrm{\nabla }}^{2}u-\frac{1}{c_0^2}\frac{{\partial }^{2}u}{\partial t^2}+{\tau }_{\sigma }^{\alpha }\frac{{\partial }^{\alpha }}{\partial t^{\alpha }}{\mathrm{\nabla }}^{2}u-\frac{{\tau }_{ϵ}^{\beta }}{c_0^2}\frac{{\partial }^{\beta +2}u}{\partial t^{\beta +2}}=0.$$

Такі моделі пов'язані із загальновизнаною гіпотезою про те, що в складних середовищах явища множинної релаксації призводять до згасання.

Дробове рівняння Шредінгера у квантовій теорії

Дробове рівняння Шредінгера має такий вигляд: $$i\hslash \frac{\partial \psi (𝐫,t)}{\partial t}=D_{\alpha }{(-{\hslash }^2\mathrm{\Delta })}^{\frac{\alpha }{2}}\psi (𝐫,t)+V(𝐫,t)\psi (𝐫,t),$$ де Шаблон:Math — хвильова функція, а Шаблон:Mvar — зведена стала Планка. Функція потенціальної енергії Шаблон:Math залежить від системи.

Шаблон:Mvar — стала з фізичною розмірністю Шаблон:Math, (при Шаблон:Math, $D_2=\frac{1}{2m}$ для частинки з масою Шаблон:Mvar). Оператор Шаблон:Math є 3-вимірною дробовою квантовою похідною Ріса, яку визначають як $$(-{\hslash }^2\mathrm{\Delta }{)}^{\frac{\alpha }{2}}\psi (𝐫,t)=\frac{1}{(2\pi \hslash {)}^3}\int d^{3}p{e}^{\frac{i}{\hslash }𝐩\cdot 𝐫}|𝐩{|}^{\alpha }\phi (𝐩,t).$$

Індекс Шаблон:Mvar у дробовому рівнянні Шредінгера є індексом Леві, Шаблон:Math.

Дробове рівняння Шредінгера змінного порядку

Як природне узагальнення дробового рівняння Шредінгера, дробове рівняння Шредінгера змінного порядку використовують для вивчення дробових квантових явищ: $$i\hslash \frac{\partial {\psi }^{\alpha (𝐫)}(𝐫,t)}{\partial t^{\alpha (𝐫)}}={(-{\hslash }^2\mathrm{\Delta })}^{\frac{\beta (t)}{2}}\psi (𝐫,t)+V(𝐫,t)\psi (𝐫,t),$$ де оператор Шаблон:Math — дробова квантова похідна Ріса змінного порядку.

Див. також

• Простір Соболєва
• Метод Кранка — Ніколсон
• Неньютонівська рідина
• Шахер Момані

Література

• Шаблон:Книга