Дробный анализ на основе оператора Адамара Текст научной статьи по специальности «Математика»

УДК 517.3

ДРОБНЫЙ АНАЛИЗ НА ОСНОВЕ ОПЕРАТОРА АДАМАРА

В.А. Чуриков

Томский политехнический университет E-mail: vachurikov@list.ru

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

1. Введение

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

Строить дробный анализ, можно используя различные подходы, многие из которых описаны в литературе [1].

Один из наиболее простых способов построения дробного анализа предложил Адамар с помощью введённого им оператора дробного дифференцирования степенных рядов [2]. Подход Адамара, названный программой Адамара, не получил развития, хотя он представляется более естественным и интуитивно понятным, чем многие другие подходы.

Оператор Адамара можно использовать для получения, как производных, так и интегралов любого конечного вещественного порядка и действует над некоторым классом функций. Более того, простота и логичность подхода Адамара даёт основание предполагать, что в нём заложен большой потенциал для построения самодостаточных направлений дробного анализа, таких же полноценного, как и традиционный анализ. Причём, возможных направлений в дробном анализе, которые можно развивать, бесконечное множество (мощности континуума), среди которых «обычный» анализ является лишь частным, вырожденным случаем, что делает его одним из самых простых вариантов дробного анализа.

2. Оператор Адамара

Определение. Оператор й‘х порядка д, де И будем называть оператором Адамара порядка д, действующим над множеством степенных функций х", д, х, деИ, д, д=сош1; [2]

dsx: xq = T(q +1 x4

(1)

r(q +1 + s)

Величину ^(q + ^ будем называть коэффи-l(q +1 + s)

циентом оператора Адамара, в котором Г(...) гамма-функция Эйлера [3].

Такие операторы обычно называют операторами дробного интегрирования или дробного дифференцирования. Операторы Адамара тоже являются

операторами дробного интегрирования или дробного дифференцирования в зависимости от значения порядка деИ.

Случаи д<0 соответствует операторам дробного дифференцирования порядка д, которые будем обозначать, как й~°х;

Для д>0, будут операторы дробного интегрирования порядка д, которые будем обозначать, как й+х или й*х.

В дальнейшем под дробными порядками дифференцирования и интегрирования будем понимать любые вещественные порядки деИ.

Определение. Операторы дифференцирования й-8х и интегрирования йдх будем называть обратными друг другу, если модули их порядков равны, |д|=|д|, а сами порядки имеют противоположные знаки, д=-д.

Из пары обратных операторов один является оператором дифференцирования, а другой оператором интегрирования.

Очевидно, что оператор нулевого порядка й °х является обратным для самого себя, или самооб-ратным.

В частности справедливо равенство й-’х-й’х=й>х=1.

Рассмотрим важные частные случаи операторов Адамара, в зависимости от д.

В задачах с дробными операторами могут участвовать операторы с различными сочетаниями порядков операторов Адамара.

Определение. Если в рассматриваемых задачах дробного анализа участвует только обратные операторы порядков д и -д, или хотя бы один из них, то будем говорить о ветви дробного анализа порядка д.

Определение. Задачи, в которых имеют место несколько различных ветвей у/>1 дробного анализа, тогда будем говорить о смешанном дробном анализе с \у ветвями.

Кратко рассмотрим некоторые важные ветви дробного анализа.

Ветвь дробного анализа с оператором нулевого порядка д=0, представляет вырожденный случай, когда имеет место только единичный оператор 1, переводящий функции сами в себя й °х: хд=1-хд=хд.

Очевидно, что порядки операторов Адамара могут иметь как целочисленные, так и нецелочисленные значения, которые, в свою очередь могут быть рациональными и иррациональными.

Определение. Если порядки дифференцирования и интегрирования в рассматриваемой ветви дробного анализа будут принадлежать только множеству натуральных чисел деИ, то такой анализ будем называть целочисленным дробным анализом, ветви которого могут принимать чётные и нечётные значения порядков.

Важный частным случаем целочисленного дробного анализа является анализ с единичным порядком интегрирования и дифференцирования, д=± 1, т. е. традиционный анализ.

В этом случае из оператора Адамара получаются формулы дифференцирования и интегрирования степенных функций традиционного анализа: й-х: хд=дхд-1, йх: хд=(д+1)-1-хд+1.

Случай, когда степень функции равна -1, а порядок оператора 1, определяется, как обычно й 'х: х^=1п х+сош1.

Другие целочисленные ветви дробного анализа является ветвями целочисленных порядков оператора Адамара д=±2, ±3, ±4,... Данные ветви анализа аналогичны случаю традиционного анализа (д=±1), но не тождественны ему.

В этом случае, используя формулу Г(д+и)=д(д+1)(д+2)...(д+и-1)Г(д), деК, ие^, получим формулу для интегрирования

ёпх: хч = ((ч + 1)(ч + 2)... (ч + п))-1 Хч +.

Формулу дифференцирования легко получить, используя формулу Г(д-и)=Г(д)(д(д-1)(д-2)...(д-и))-1, тогда получим

ё~пх: хч = ч(ч - 1)(Ч _ 2)...(ч - п) X4 .

Важными случаями дробного анализа являются ветви с рациональными порядками де р.

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

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

Интересными и наиболее простыми случаями рационального анализа представляются ветви, в которых порядки операторов обратно пропорциональны натуральным числам д=1/Я, АеМ В данном случае знаменатели Я значений порядков операторов является основным параметром данных ветвей дробного анализа.

Общими и более сложными случаями рационального анализа являются ветви, в которых порядки операторов можно представить как отношение натуральных чисел ¡=%/Я, %, Яе№

Наибольшее количество ветвей дробного анализа будет в случае иррациональных порядков оператора Адамара s.

Определение. Операторы Адамара dsx, у которых порядки s являются иррациональными числами, называются иррациональньми дробнъми операторами.

Определение. Операторы Адамара, в котором действуют иррациональные дробные операторы, будем называть иррациональным дробным анализом.

Иррациональные ветви дробного анализа представляются наиболее сложными для исследований.

Важно так же рассматривать отдельно модули порядков операторов, а именно, больше или меньше единицы |s|>1 или |s|<1.

Определение. Два оператора порядков s и 1/q называются обратными по порядку, если произведение их порядков даёт единицу sqr1=1.

Каждый оператор порядка s имеет обратный по порядку оператор со значением порядка 1/s.

Операторы традиционного анализа d ±!x (s=±1) являются обратными по порядку сами для себя.

Любая пара взаимно обратных по порядку операторов являются одновременно операторами дифференцирования или операторами интегрирования.

Определение. Для оператора dsx имеется противоположный - dsx, такой, что в сумме они дают нулевой оператор

dsx+(-dsx)=dsx-dsx=0.

3. Свойства оператора Адамара

Производная дробного порядка от константы. В

общем случае в дробном анализе производная константы не равна нулю. Константу можно представить как полином нулевой степени, или как степенную функцию, с порядком равным нулю Cx", где C вещественное число (CeR), тогда производная константы будет

d sx: Cx0 = C(-sr(-s))-1 x-.

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

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

В случае, когда порядок оператора дифференцирования и показатель степени функции равны и имеет порядок a, aeR

d-ax: xa = аГ(а) = const.

Интеграл дробного порядка от константы и интегральные полиномы

dsx: Cx0 = C(sr(s))-1 xs + Cs(x).

Здесь введены функции Cs(x)^0, которые появляются как слагаемые при интегрировании

функций и являются обобщением констант интегрирования стандартного анализа.

Определение. Функции С(х) будем называть интегральным полиномом порядка д для пары обратных операторов йх и й-х.

Основным свойством интегральных полиномов С(х) должно быть то, что при действии на них оператором дифференцирования й-х должен получаться ноль

й-х С/х)=0.

Используя это равенство, найдём степени членов интегрального полинома Сд(х) произвольного вещественного порядка д

й -х: хч = Г(ч +1) хч- = 0.

Г(ч +1 - 5)

Выполнение этого равенства возможно при соблюдении двух условий

Г(д+1)^0 и Г(д+1- з)=<».

Первое выполняется при д^-1, -2, -3, -4 ...

Для выполнения второго условия необходимо, чтобы выполнялось равенство

д+1-д=0, -1, -2, -3, -4 ...

Тогда получим д=-1+д, -2+д, -3+д, -4+д ...=-и+&

Окончательно интегральный полином Сд(х) для произвольных вещественных порядков д будет

да

С (X) = Х “пХ-п+5 =

п=1

= а1х~1+5 + а2 х~ 2+5 + а3 х~35 + а4 х~4+5 +... + апх~п+‘' +...

У полиномов дробного интегрирования С(х) бесконечное счётное множество констант интегрирования а1, а2, а3, ..., которые, в общем случае, являются вещественными (или комплексными) числами.

Для произвольных целочисленных порядков операторов дифференцирования т=1, 2, 3, ... и целочисленных показателей степеней степенных функций с показателями к>0 и к=0, 1, 2, 3, ..., полиномы интегрирования Ст(х) будут

й-тх: хк = Г(к +1) хк-т = 0.

Г(к +1 - т)

Выполнение этого равенства возможно при соблюдении двух условий.

Первое условие Г(к+1)^0 выполняется при к^-1, -2, -3, -4 ...

Для выполнения второго условия Г(к+1-т)=да необходимо, чтобы выполнялось равенство к+1-т=0, -1, -2, -3, -4 ..., которое можно переписать так к=-1+т, -2+т, -3+т, -4+т ...

В случае, когда выполняется условие к+1-т<0 или когда к<т-1, тогда значения переменной в этих точках соответствуют полюсам гамма-функции, и поэтому она обращается в бесконечность. Это значит, что для случаев целочисленных порядков показателей операторов дифференцирования и полино-

мы интегрирования будут обрываться в случае, когда k<m-1, а степени будут ограничиваться целочисленными значениями при k<m-1, m-2, m-3, ... 2, 1, 0.

Окончательно получим, что для любого оператора целочисленного порядка d -mx полиномы интегрирования Cm(x) будут иметь конечное число m слагаемых, равное порядку оператора интегрирования

m

Cm (x) =Z akxk = k=0

= a0 + a1x + a2 x2 + a3 x3 + ... + am-2 xm-2 + am-1xm-1.

Здесь даны m констант интегрирования a0, a1, a2,

... am^ am-1eR.

В случае традиционного анализа, для оператора дифференцирования d -1x равенство выполняется для случая k=0. Тогда полином интегрирования будет константой C1(x)=a0=const, a0eR.

В полиномах интегрирования с нецелым порядком необходимо задавать бесконечное счётное множество констант интегрирования, которые в различных зачах необходимо находить или задавать. Ряд Cs(x) может быть расходящимся или сходящими, а скорости сходимости и расходимости могут быть разными. Чем быстрей сходимость, то в зависимости от условий задачи, можно ограничиться несколькими первыми членами ряда Cs(x).

В случае нецелочисленных порядков s члены ряда Cs(x) могут обращаться в бесконечность в точке x=0.

Интегральным полиномом для оператора нулевого порядка (единичного или самообратного оператора) d °x является ноль.

В полиномах интегрирования Cs(x) и Cm(x) соседние слагаемые имеют показатели степеней, которые отличаются на единицу (с шагом 1).

Неопределённый интеграл дробного порядка

Из всего вышеизложенного можно ввести понятие неопределённого интеграла любого вещественного порядка от функции fx) d sx.f(x)=Fi)(x)+Cs(x).

Здесь Fs)(x) - первообразная порядка s функции fx).

4. Экспоненты дробных операторов

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

Традиционная экспонента, для дробных порядков оператора Адамара уже теряет своё главное свойство - не меняться при воздействии оператора дифференцирования или интегрирования (с точностью до сложения с полиномом интегрирования). В этом легко убедиться, например, продифференцировав экспоненту exp x оператором Адамара порядка 1/2, т. е., d~1/2x: exp x^exp x.

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

Прежде, чем получить экспоненты для операторов Адамара любого порядка, обратим внимание на свойства традиционной экспоненты exp x.

Свойство членов ряда a(x) экспоненты exp x при обычном дифференцировании и интегрировании (порядки операторов ¿=±1) имеют следующие свойства:

При интегрировании каждый член ряда переводится в последующий d X: а;(х)=а;+1(х).

При дифференцировании все члены ряда, кроме первого, переводятся в предыдущий d-1x: a;+i(x)=a;(x), Ы.

Производная первого члена ряда экспоненты равна нулю d-1x: a1(x)=0, что легко получить, расписав это более подробно

d -x

1

г(1 -1 + s) _x-1-s-s

d 1 x: a1(x) =

-1- x-i = ^ = 0. Г(0) да

r(s) Г(s) Г(1 -1 + s- s)

= — x1 = 0.

да

Используя определение экспоненты любого вещественного порядка $¿0 можно получить, используя первый член разложения интегрального полинома Cs(x), т. е. функцию r-1(s)x~1+s, которую назовём стартовой функцией.

Найдём экспоненты expx для всех пар обратных операторов d±sx любого вещественного порядка $¿0

exps x = (d0x + dsx + (dsx)2 + (dsx)3 + ... +

+(dsx)n +...) x-----=|£ (dsx)

r(s)

r(s)

Здесь введён символ (dsx)n=dsx:dsx:dsx:...dsx:dsx, который обозначает последовательное действие п операторов дробного интегрирования на одну функцию.

После почленного дробного интегрирования получим ряд для экспоненты ехрх

exPs x = Z

Последнее свойство говорит о том, что первый член ряда экспоненты пропорционален первому члену полинома интегрирования. Но ввиду того, что в традиционном анализе полиномом интегрирования является константа, то тогда первый член ряда пропорционален константе. Данной константой должна быть единица, чтобы обеспечить равенство dx a1(x)=a2(x).

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

exp x = (d0 x + d1 x + d2 x + d3 x +... + dnx + ...)x0 =

= (1 + d x + d x + d x +... + d x + ...)1 =

^ m ^ , x x2 x3 x4 xn

> dnx 1:1 = 1 + —+— + — + — +... + — +... =

J 1! 2! 3! 4! n!

да vn v-» Л

= Z nr

Распространим аналогичные свойства на экспоненты любого вещественного порядка s, которые в дальнейшем будем обозначать как expx.

В степенном ряде экспоненты expx для пары обратных дробных операторов Адамара d ±sx, между соседними членами имеют место соотношение при дробном интегрировании dsx: a;(x)=a;+1(x). Первым членом ряда является ненулевая функция a1(x)=r-1(s)x-1+s, которая пропорциональна первому члену ряда полинома интегрирования Cs(x), при действии на который оператора дифференцирования порядка s d-sx, должно выполняться равенство, справедливое для традиционной экспоненты, а именно

x-1+s x-1+ 2s-x-1+ 3s-----x- + + +... + -

r(s) Г(2 s) Г (3 s)

Г( ns)

- +...

Определение. Экспоненту ехрх будем называть дробной экспонентой.

Дробная экспонента, по сути, представляет бесконечное множество экспонент любого конечного вещественного порядка 5. Мощность этого множества равна мощности множества всех возможных порядков оператора Адамара, или мощность континуума.

Определение. Для конкретных вещественных порядков 5 дробную экспоненту будем называть частными экспонентами порядка 5.

С помощью дробной экспоненты можно получить экспоненту любого вещественного порядка для любой пары обратных операторов, подставив вместо 5 конкретное значение модуля их порядков 5| обратных операторов.

Теорема. Для каждой пары обратных операторов Адамара d±sx имеется своя частная экспонента ехрх, причём единственная и, отличная от экспонент других пар обратных операторов, или ехрх=ехр,х, 5=# и ехрх^ехрх, 5^.

Интеграл порядка 5 от экспоненты ехрх будет

d~sx:expx=expx+ С/х).

Дифференцируя правую часть оператором d~sx, получим d~sx:(expx+C/x))=expx, и в частности

d~sx:expx=expx.

Рассмотрим свойства дробной экспоненты.

Теорема. Ряд дробной экспонента expx с порядками 5>1 является сходящимся, с радиусом сходимости Я равным бесконечности (Я=<»).

Теорема. Ряд дробной экспоненты expx с порядками s<1 имеют особую точку x=0, в которой ряд расходится, а в остальных точках являются сходящимися рядом, с радиусом сходимости R, равным бесконечности (R=<»).

Заметим, что все степенные ряды с положительными целочисленными степенями в точке x=0 всегда сходятся [4].

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

cosax = 2(expa (ix) + expa (-ix))

sin« x = ^(expa (ix) - exPa (-ix)),

2i

sin x cos„ x

tga x=—— ,ctga x =—

cos„x

sin„x 1

sec „X = -

cos „X

sin„ X

Здесь даны обозначения тригонометрических функций порядка 5: со8^ sinsx, tgx, ctgx, secsx, cosecsx соответственно косинус, синус, тангенс, котангенс, секанс, косеканс.

Гиперболические функции можно ввести так

С^Х = 2(ехРа(х) + exPа (-х))

ShaX = 2(еХРа (Х) - еХРа (-Х))>

^ах сКх

th„x = ——, сш„ х = -

sch„x =

ch„x ch„x’

csch„ x =

sh„x 1

sh„x

Здесь гиперболические функции порядка 5: chsx, shsx, cthsx, schsx, cschsx, соответственно косинус,

синус, тангенс, котангенс, секанс, косеканс.

Определение. Функции cos.x, sinx, tgsx, ctgsx, secsx, cosecsx будем называть дробными тригонометрическими функциями, а chx, shsx, thsx, schsx, cschsx -

дробными гиперболическими функциями порядка 5.

При таком определении дробные тригонометрические и гиперболические функции в частном случае, когда порядок их 5=1, дают соответствующие традиционные тригонометрические и гиперболические функциям.

СПИСОК ЛИТЕРАТУРЫ

1. Сашко С.Г., Килбас А.А., Маричев О.И. Интегралы и производные дробного порядка. - Минск: Наука и техника, 1987. -687 с.

2. Hadamar J. Essai sur l’etude des fonctions donnees par leur development de Taylor // J. math. pures et appl. - 1892. - V. 8. - Ser. 4. -P. 101-186.

3. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. Т. 1. - М.: Наука, 1973. - 295 с.

4. Кудявцев Л.Д. Курс математического анализа. Т 1. - М.: Высшая школа, 1981. - 687 с.

Поступила 07.11.2007г.