Научная статья на тему 'ИНТЕГРАЛЫ И ПРОИЗВОДНЫЕ ДРОБНОГО ПОРЯДКА НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ПРЕОБРАЗОВАНИЙ ТИПА ЛАПЛАСА С ПРИЛОЖЕНИЯМИ'

ИНТЕГРАЛЫ И ПРОИЗВОДНЫЕ ДРОБНОГО ПОРЯДКА НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ПРЕОБРАЗОВАНИЙ ТИПА ЛАПЛАСА С ПРИЛОЖЕНИЯМИ Текст научной статьи по специальности «Математика»

CC BY
113
25
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ИНТЕГРАЛ И ПРОИЗВОДНАЯ ДРОБНОГО ПОРЯДКА / ОБОБЩЕННОЕ ИНТЕГРАЛЬНОЕ ПРЕОБРАЗОВАНИЕ ЛАПЛАСА / ОПЕРАТОР ПРЕОБРАЗОВАНИЯ / СВЕРТКА ФУНКЦИЙ

Аннотация научной статьи по математике, автор научной работы — Яремко Олег Эмануилович, Яремко Наталия Николаевна

Развивается теория интегралов и производных дробного порядка. Построен аналог операционного исчисления для дифференциального оператора с кусочно-постоянными коэффициентами. Предложены различные конструкции обобщенного преобразования Лапласа. При помощи операторов преобразования установлена связь интегральных преобразований Меллина - Лапласа с обобщенным интегральным преобразованием Лапласа. Найден изоморфизм пространства оригиналов и пространства обобщенных оригиналов. Установлены формулы обращения типа Меллина - Лапласа. Доказаны теоремы о дифференцировании обобщенного оригинала и другие. Дано определение обобщенной свертки и установлена формула для ее вычисления, указана связь обобщенной и классической свертки. На основе понятия обобщенной свертки дано определение обобщенного интеграла и обобщенной производной дробного порядка. Установлены соотношения между обобщенными интегралами дробного порядка и интегралами Римана - Лиувилля дробного порядка. Для модельного уравнения теплопроводности с кусочно постоянным коэффициентом решена задача вычисления плотности теплового потока. Тепловой поток выражен в виде обобщенной производной порядка 1/2 по времени от измеренной зависимости температуры на границе.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Яремко Олег Эмануилович, Яремко Наталия Николаевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INTEGRALS AND DERIVATIVES OF FRACTIONAL ORDER BASED ON LAPLACE TYPE INTEGRAL TRANSFORMATIONS WITH APPLICATIONS

Integrals and derivatives of fractional order theory is being developed. An analogue of the operational calculus is constructed for a differential operator with piecewise constant coefficients. Various constructions of the generalized Laplace transform are proposed. The transformation operators establish a connection between the Mellin - Laplace integral transformations and the generalized Laplace integral transformation. Isomorphism between the space of originals and the space of generalized originals is found. Mellin - Laplace type inversion formulas are established. Theorems on the differentiation of the generalized original and others are proved. A definition of a generalized convolution is given and a formula for its calculation is established, a connection between the generalized and classical convolution is indicated. On the basis of the concept of generalized convolution, a definition of a generalized integral and a generalized fractional derivative is given. Relations between generalized fractional integrals and Riemann-Liouville integrals of fractional order are established. For a model equation of heat conduction with a piecewise constant coefficient, the problem of calculating the heat flow transfer is solved. The heat flow is expressed as a generalized time derivative of the order of 1/2 of the measured temperature dependence at the boundary.

Текст научной работы на тему «ИНТЕГРАЛЫ И ПРОИЗВОДНЫЕ ДРОБНОГО ПОРЯДКА НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ПРЕОБРАЗОВАНИЙ ТИПА ЛАПЛАСА С ПРИЛОЖЕНИЯМИ»

Прикладная математика & Физика, 2021, том 53, №2. С. 114-124.

УДК 517.44 DOI 10.52575/2687-0959-2021-53-2-114-124

MSC 44А10.

ИНТЕГРАЛЫ И ПРОИЗВОДНЫЕ ДРОБНОГО ПОРЯДКА НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ПРЕОБРАЗОВАНИЙ ТИПА ЛАПЛАСА С ПРИЛОЖЕНИЯМИ

О. Э. Яремко, Н. Н. Яремко

Статья представлена членом редакционной коллегии С. М. Ситником

Московский государственный технологический университет (Станкин), Москва, 127055, Россия Пензенский государственный университет, Пенза, 449926, Россия

E-mail: yaremki8@gmail.com

Аннотация. Развивается теория интегралов и производных дробного порядка. Построен аналог операционного исчисления для дифференциального оператора с кусочно-постоянными коэффициентами. Предложены различные конструкции обобщенного преобразования Лапласа. При помощи операторов преобразования установлена связь интегральных преобразований Меллина - Лапласа с обобщенным интегральным преобразованием Лапласа. Найден изоморфизм пространства оригиналов и пространства обобщенных оригиналов. Установлены формулы обращения типа Меллина - Лапласа. Доказаны теоремы о дифференцировании обобщенного оригинала и другие. Дано определение обобщенной свертки и установлена формула для ее вычисления, указана связь обобщенной и классической свертки. На основе понятия обобщенной свертки дано определение обобщенного интеграла и обобщенной производной дробного порядка. Установлены соотношения между обобщенными интегралами дробного порядка и интегралами Римана - Лиувилля дробного порядка. Для модельного уравнения теплопроводности с кусочно постоянным коэффициентом решена задача вычисления плотности теплового потока. Тепловой поток выражен в виде обобщенной производной порядка 1/2 по времени от измеренной зависимости температуры на границе.

Ключевые слова: интеграл и производная дробного порядка, обобщенное интегральное преобразование Лапласа, оператор преобразования, свертка функций.

Для цитирования: Яремко О. Э., Яремко Н. Н. 2021. Интегралы и производные дробного порядка на основе интегральных преобразований типа Лапласа с приложениями. Прикладная математика & Физика, 52(2): 114-124. DOI 10.52575/2687-0959-2021-53-2-114-124.

INTEGRALS AND DERIVATIVES OF FRACTIONAL ORDER BASED ON LAPLACE TYPE INTEGRAL

TRANSFORMATIONS WITH APPLICATIONS

Oleg Iaremko, Natalia Iaremko

Article submitted by a member of the editorial board S. M. Sitnik

Moscow State Technological University (Stankin), Moscow, 127055, Russia, Penza State University, Penza, 449926, Russia E-mail: yaremki@mail.ru Received May, 17, 2020

Abstract. Integrals and derivatives of fractional order theoiy is being developed. An analogue of the operational calculus is constructed for a differential operator with piecewise constant coefficients. Various constructions of the generalized Laplace transform are proposed. The transformation operators establish a connection between the Mellin - Laplace integral transformations and the generalized Laplace integral transformation. Isomorphism between the space of originals and the space of generalized originals is found. Mellin - Laplace type inversion formulas are established. Theorems on the differentiation of the generalized original and others are proved. A definition of a generalized convolution is given and a formula for its calculation is established, a connection between the generalized and classical convolution is indicated. On the basis of the concept of generalized convolution, a definition of a generalized integral and a generalized fractional derivative is given. Relations between generalized fractional integrals and Riemann-Liouville integrals of fractional order are established. For a model equation of heat conduction with a piecewise constant coefficient, the problem of calculating the heat flow transfer is solved. The heat flow is expressed as a generalized time derivative of the order of 1/2 of the measured temperature dependence at the boundary.

Key words: fractional integral and derivative, generalized integral Laplace transform, transformation operator, convolution of functions.

For citation: Oleg Iaremko, Natalia Iaremko. 2021. Integrals and derivatives of fractional order based on Laplace type integral transformations with applications. Applied Mathematics & Physics, 52(2): 114-124. (in Russian) DOI 10.52575/2687-0959-2021-53-2-114-124.

Введение. Пусть функция f{t) задана на промежутке [0, оо) и является оригиналом, а функция .Р(р) ее изображение. Интеграл Римана - Лиувилля можно определить [15] формулой

1ат=1-Л\нР) р

С помощью теоремы о свертке последняя формула преобразуется к форме Римана - Лиувилля

I

Ш = тЬ /

где Г — гамма-функция, а — комплексное число в полуплоскости К.е а > 0. Для функции f{t) определенной на промежутке [0, оо) производная дробного порядка а определяется по правилу

¿«[/(0] = [/«]," = [а] + 1.

Таким образом, дифференциальное исчисление дробного порядка можно построить на основе операционного исчисления. Теория дифференциального и интегрального исчисления дробного порядка рассматривается в работах [10-11], [15]-[16].

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

1 <1

=

x'(t) dt'

где функция х = x(t) имеет вид линейного сплайна с узлами to = 0, tk, k = 1,.., п

х (t) = а^ (Я (t - to) - H (t - fi)) + (a2 (t - h) + ^Ah) (Я (t - h) - H (t - t2)) +

+ (a3 (t - h) + a2At2 + aiAti) (Я (t - tt) - H (t - t2)) + ... + (a„+i (t - t„) + anAtn + ... + aiAti) H (t - tn), здесь Atk = tk ~ tfc-i, fo = 0. Оператор Dt возникает при моделировании переходных процессов.

Методы. 1.1 Обобщенное преобразование Лапласа на основе оператора дифференцирования с кусочно-постоянным множителем. Методы интегральных преобразований типа Лапласа можно использовать для построения обобщенного дифферинтеграла Римана - Лиувилля. Обобщениям операционного исчисления посвящены работы [1], [3], [4], [6]-[9], [12]-[14], [17]-[20]. Идею конструирования ядер обобщенных преобразований Лапласа на основе оператора дифференцирования с кусочно-постоянным множителем проиллюстрируем на примере классического преобразования Лапласа. Как известно, ядро в формуле обращения Меллина [20] имеет вид K{t,p) = ept, и, значит, удовлетворяет уравнению

K[(t,p) = pK(t,p),0 <t<oo

и начальному условию

К(0,р) = 1.

Будем конструировать ядро K(t,p) обобщенного преобразования Лапласа в виде решения сепаратной системы дифференциальных уравнений

K't{t,p) = akpK(t,p),tk-i <t <tk,k = 1,2, ...,n + l;f0 = 0, tn+1 = оо,

удовлетворяющих начальному условию

К(0,р) = 1

и условиям непрерывности в точках деления tk, k = 1, 2,..., п

K(tk-0,p)=K(tk + 0,p).

Вычисления приводят к формуле:

K{t,p) = eaipt, О < t < h, K(t,p) =ea2i(t-tl)+Atl,ii < t < t2,

K(t p) = e~il-a"+ll-i~i"-l+a"Ai"+-"+aiAil-1 t„ < t

где Atk = tic ~ tk-i- Пусть функция-оригинал fit) определена на промежутке [0,оо). Будем искать функцию-изображение F{p) как решение интегрального уравнения

^ /»a+ioo

f(t) = — K(t,p)F(p)dp,0<t <h.

¿Я1 J (j-ioo

» развернутом виде получаем сепаратную систему интегральных уравнении

» (7+100

А /» (7+IOO

/(t) = — / ep^F(p)dp,0<t <tb

-t /» a+ioo

f{t) = — / ep(a2{t~h)+aiAh) F (p)dp, h<t<t2,... (1)

27ri Ja-ioo

-I /• a+ioo

f(t) = - / eP(an+l(t-tn)+anAtn+...+a1At1) p ^^ ^ < f

2Я"! Ja-ioo

Выполним замену переменного в последней формуле х = x(t), где функция х = х (t) имеет вид линейного сплайна с узлами to = 0, tt, k = 1,.., п,

х (0 = att (Я (t - t0) - Н (t - ti)) + (а2 (t - ti) + aiAti) (Я (t - h) - H (t - t2)) +

+ (a3 (t - t2) + a2At2 + aiAti) (Я (t - ti) - Я (t - t2)) + ... + (a„+i (t - t„) + a„At„ + ... + aiAti) Я (t - t„), здесь At,t = tjt - tjt-i, to = 0. График сплайн-функции x = x(t) является ломаной с вершинами в точках

М0 (0,0), Mi (ti, aiAti), М2 (t2, a2At2 + aiAti),

...,Mn (t„,a„At„ + ... + aiAti).

Угловые коэффициенты звеньев ломаной равны, соответственно ai, a2,..., a„+i. Заметим, что в частном случае ai = 1,a2 = 1, ...,a„+i = 1 ломаная становится прямой. Функция t = t(x), обратная к функции х (t), есть полиномиальный сплайн первого порядка вида

t(x) = — (Я (х) - Я (х - aiAti)) + ai

х atAtt + ^ ^ ^ _ ajAti _ а2д?2)) +

+ t„ Я (x - aiAti - ... - a„At„).

a2

x - aiAti - ... - a„At„

an+i

Тогда из формулы (1) получим

__Л С (7+iOD

/(х) =f(t(x)) = — epxF (p)dp.

¿Я I J (J-ioo

По формуле для прямого преобразования Лапласа [20] получим выражение для обобщенного преобразования Лапласа

/• со

Яр) = / e~pxf(x)dx. иъ

Приведем строгие определения.

Определение 1. Функцию f (0 назовем обобщенным оригиналом, если функция f{x) = / (х)) является функцией-оригиналом.

Определение 2. Оператор ] : / —» f, действующий по правилу f(t) = f{x{t)), назовем оператором преобразования.

Определение 3. Изображением обобщенной функции-оригинала / (f) назовем выражение

L[/(f)] =l[/(f)],

где L-onepamop преобразования Лапласа.

Теорема 1. Оператор J'1, обратный к оператору преобразования, действует по правилу J_1 : /(f) —» /(f), zdef(t) = f (х-1 (f)), причем выполняются равенства

L = L -Г\1 = L-J . Теорема 2. Обобщенное преобразование Лапласа находится по формуле

/» со

L [/ (f)] = F (р) = / e-Px(f>x'(t)f(t)dt, J о

где функциях' (f) кусочно-постоянная вида

х' (f) = ai (Я (f - f0) - Я (f - fO) + a2 (Я (f ~h) -H(t - f2)) +

+a3 (Я (f - fi) - Я (f - f2)) + ... + an+iH (f - f„), Доказательство. По определению 1 имеем

/» со

L[f(t)]= e~Pxf(x)dx. J о

В приведенном интеграле выполним замену переменного х = x(t). Тогда получим

/» со

L[f(t)] = / e-Px^f(x(t))x'(t)dt. J о

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

Следствие 1. Пусть f(t) - функция-оригинал, тогда обобщенное изображение Лапласа вычисляется по формуле

F (р) = ai Г e~aiptf (t)dt + а2 [ ' e-PMt~ti)+aiAti) f (f) dt+ Jo J tt

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+a3 f 3 е-рМ'-ьНагЫг+аг&ь) f (f) Л + _

Jt2

/• CO

+a„+i / e-p(an+l(t-tn)+anAtn+...+a1At1) f

Jtn

Приведем аналог формулы Римана - Меллина.

Теорема 3. Если / (f) оригинал, то интеграл (1) сходится всюду в полуплоскости Rep > s0. При этом сходимость равномерная в любой области Rep > s > s0. Функция F (р) аполитична в полуплоскости Rep > s0. Имеет место формула обращения

/»a+ico

L~1[F(p)]=f(t)= epx<-^F(p) dp.

J (j—ioo

На основании теоремы 3 и определения функции х = x(t) получается

Следствие 2. Формула обращения Римана - Меллина для обобщенного преобразования Лапласа

л /» (j+ioo

fi(.t) = — eP^F(p)dp,0<t<tl,

¿71 J (j-ioo

л /» (J+iO0

f2 (f) = — / ep{a2{t-tl)+aihtl)F(p)dp, h<t< t2,

2тг J (J-ioo -j /»a + ioo

fn+1 (f) = — / eP(*n+l(t-tn)+anhtn+...+a1ht1)F^p-)dpJn < f

Ju-ioo

1.2 Дифференцирование обобщенного оригинала.

1 d '(t) dt

Теорема 4. Пусть Б = -гттрг4: оператор обобщенного дифференцирования, а функция Б [/(01 есть

или

обобщенный оригинал, тогда

L[D (/(f))] =-m+pF(j>).

Доказательство. Воспользуемся определением обобщенного преобразования Лапласа. Тогда, интегрируя по частям, получим

/» СО -, /» со

/ e-Px(f>x'(t)—f'(t)dt = -/(0) + р / e-Wx'WWt.

JО х \t) Jо

Свертка оригиналов.

Определение 4. Сверткой функций оригиналов/ (t), g (t) назовем функцию

f(t)*g(t)=L-1 [F(p)G(p)].

Теорема 5. Свертка вычисляется по формулам

fit) *g(t)= Г f (х"1 (х it) - х (г))) х' (г) g (г) dr J о

fit)*git)=j[fit)*git)]. Доказательство. По определению 1 запишем

/»ОО /» t

L [f it) * g (f)] = / it) / / {x-1 ix (f) - x (r))) x' (r) g (r) drdt.

J 0 J 0

Выполним замену переменного x (t) = j]

/•со /» t

L [f it) * g it)] = / e-rf / / (x"1 (p-x (r))) x' (r) g (r) drdfi. J 0 J 0

Проведем вторую замену x (г) = a

/» OO /*x_1(f7)

L [f it) * g (f)] = / e-P? / / (x"1 (/? - a)) g (x"1 (a)) J 0 J 0

По теореме о свертке для преобразования Лапласа [20] получим

L [f it) * g it)] =L[f (x"1 (a))] L [g (x"1 (a))] = L[f (t)] L [g (t)] .

Докажем, что изображения правой и левой частей во второй формуле из теоремы 5 одинаковы. По определению 4 имеем равенство

L[fit)*git)]=F(p)G(p).

С другой стороны, получим

L ■ J [fit) * git)] = L [f(t) * g(t)] =F{p)G (p). Для доказательства второй формулы в теореме 5 перепишем ее правую часть в виде

/"/((* (0 -х(т)))д(х(т))х' (т) dr. J а

Затем выполним замену переменной г = х-1 (z). Тогда

Ф)

rxw _ г1 _

/ f((x(t)-z))g(z)dz = J / f((t-z))g(z)dz J о |Уо

= /[/(0*i(t)] = fit)* git).

Теорема доказана.

Приведем формулу для обобщенной свертки в виде, удобном для приложений.

Следствие 3. Пусть функции f(t),g(t) являются обобщенными оригиналами, тогда справедлива формула для свертки

f(t) * g(t) = ai f (r)g (ait - air) dr, 0 < t < ti,

J о

rh

f(t)*g(t) = a-i / f (r)g (a2 (t - ti) + ai (ti - t)) dr +

J о

+a2 / f (r)g (a2 (t - r)) dr, ti < t < t2,

Jt i

/(f) * = ai f f (r)g (a„+i (t - t„) + a„ (t„ - t„_i) + ... + ai (tt - r)) dr + J о

r k

+a2 / / (r) (a„+i (t - tn) + an (tn - t„_i) + ... + a2 (t2 - r)) dr + ... + Jt!

+a„+1 / / (г) ^ (a„+i (t - r)) dr, t„ < t.

J t.

Доказательство. Преобразуем формулу для свертки из теоремы 5. Имеем

Ф)

f (t) * g (t) = J

/ f (r) g (t - t) dr = / f(T)g(x(t)-T)dT. J 0

Выполним замену в интеграле г = x(z), тогда получим

f (.t) * g (t) = Г / (z) £ (x(t) - x(z)) x'(z)dz. J 0

Подставим явное выражение функции г = x(z). В результате получим доказываемую формулу.

2. Обобщенный дифферинтеграл Римана - Лиувилля. 2.1 Обобщенный интеграл Римана -Лиу вилля.

Определение 5. Пусть функция f(t) является обобщенным оригиналом, а функция F(p)~ изображение функции f(t). Обобщенный интеграл Римана - Лиувилля Iaf определим по формуле

Гf{t) = Г1 [p~aF(p)] .

Из определения 5 следует полугрупповое свойство обобщенного интеграла Римана - Лиувилля

jajp - ¡а+Р-

Теорема 6. Обобщенный интеграла Римана - Лиувилля Iaf порядка а вычисляется по формуле

Iaf(t)=Jïaf(t).

Доказательство. Воспользуемся определением обобщенного интеграла Римана - Лиувилля и формулой L = LJ. Тогда получим

Г/(р) = L"1 [p-aF(p)] = J ■ L"1 [p-aF(p)] = Jïaf(t).

Следствие 4. Обобщенный интеграла Римана - Лиувилля и интеграла Римана - Лиувилля порядка а связаны формулой

Ia=JïaJ-\

Из определения обобщенного интеграла порядка а можно получить формулу типа Римана - Лиувилля. Теорема 7. Для обобщенного интеграла Римана - Лиувилля Iaf порядка а справедлива формула

Iam = ai [h f Ш Ш - air) (х(°~ДгГ 1 dT+ Jо Г(a)

^ Г2 Г ^ (*(t) - - h) + aiAtif-1 TT, ^

+a2 / /(г)-——--H (x(t) - (а2(т - h) + aiAti) dr+

Jh Г(a)

f°° (x(t) - (a„+i (r - tn) + a„At„ + ... + aiAti))"-1

+-+an+iL-rw--

•Н(х(0 - (а„+1 (г - £„) + а„А?„ + ... + а^))/ (г) с1т, 0 < г\

Доказательство. Пусть для простоты 0 < £ < тогда из определения обобщенного дробного интеграла Римана - Лиувилля имеем:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Г[/(0] =1"1

^G(p)

л /» CT+IOO ч

= e^-F(p)dp =

-, (T+ioo , / Л f!

+а2 [ 2 e-PMT~ti)+aiAti)f2(T)dT+

Jt 1

+а3 [ \-p(a3(T-h)+a2ht2+aihh)fc{r)dT+...

Jt2

+an+i £ e-p(an+1(T-tn)+anAtn+...+a1Ah)fn+i ф ^ = = at ji(T)H(t-T)--аг+

Jo ^

Г (a)

^ ^ Ы - (а2(т - t^ + агМгГ-1 w ^

+a2 / /2 (r)-——-H (ait - (а2(т - ti) + aiAti) ar+

Д Г (a)

/• 00

Jtn

(axt - (an+i (t - tn) + anAtn + ... + aiAti))"-1

Г (a)

■H(ait - (a„+i (r - f„) + anAtn + ... + aiAti))/„+i (r) dr, 0 < f < fb

Учитывая, что на промежутке 0 < f < t\ выполняется условие x(t) = a\t, доказательство закончим.

Следствие 5. Для случая линейного сплайнах = x(t) с двумя узлами t0, t\ обобщенный дробный интеграл Римана - Лиувилля имеет вид

Ia[f(t)]=ar Г /(г) J о

(ait - ait)"

о Г (a)

-dr, 0 < t < tx

I [/ М] = «1 / / (г)-—--

Л Г (а)

Г*, . (а2 - г) Г"1 , +а2у /(г)--dт, Ь < ^

Следствие 6. Для случая линейного сплайна х = с тремя узлами Н, Н обобщенный дробный интеграл Римана - Лиувилля имеет вид

1а[ГШ=аг ftf(T)^pf-^dт,0<t<t1, Jo Г (а)

I и(ю\=а! /(г)-—--йг+

Л Г (а)

Г*, . (а2 - г) Г"1 ,

+а2 / / (г)-—--йг, t1<t < Ь,

Jt1 Г (а)

Г'1,, Л («3 - +а2 - -г))""1

^ [/(01=^1 / /(г)-—--dт+

Jо Г (а)

, (а3 ~ Ь) + а2 (¿2 - г))""1 +а2 / J (г)-^т^-"г+

Л,

/tl Г (а)

•f („ ( t-

a. f ff Л , +a3 / / (г)-—--dг, f2 < f.

Л2 Г(а)

2.2 Обобщенная производная Римана - Лиувилля.

Определение 6. Для функции/(О определенной на промежутке [0, °о) производной дробного порядка а назовем функцию

Теорема 8. Оператор обобщенного дифференцирования В и модельный оператор дифференцирования Д Г) = удовлетворяют соотношению

т.е. Б = /ДГ1.

Если для функции/(^ выполняется полугрупповое свойство то выполняется полугрупповое свойство для функции/(^ Доказательство. Преобразуем данное в условии тождество

Тогда получим

. Свойство доказано.

Теорема 9. Фундаметалъное соотношение имеет вид

ша+1 = г

Доказательство. Преобразуем левую часть доказываемого тождества

В1а+1 = ]Ъг1]1а+1г1 = ]Ыа+1Г\

Аналогично преобразуем правую часть

Iя = irr1.

На основании модельного фундаметального соотношения DIa+1 = Ia доказательство закончим.

3. Уравнение теплопроводности с кусочно-постоянными коэффициентами. В качестве приложения рассмотрим задачу вычисления плотности теплового потока q{t) по измеренной зависимости температуры на границе по времени f{t). Рассмотрим смешанную краевую задачу Коши для однородного уравнения теплопроводности:

д2и

Dtu--- = 0, х > 0, t > О, (2)

dxz

и(х, 0) = 0,х > 0, (3)

u(0,t) =f(t),t > 0. (4)

Дополним условия (2)-(4) условием непрерывности решения в моменты t = tk, к = 1, 2,..., п.

lim u(t,x) = lim u(t,x)

Расмотрим также модельную смешанную краевую задачу Коши для однородного уравнения теплопроводности:

дй д2й , , с—г - —^ =0, х > 0, t > 0, (5)

dt дхz

м(х, 0) = О, х > 0, (6)

м(0, t) = f(t), t> 0. (7)

В работе [2] доказана формула для потока q(t) на границе х = 0

q(t)=kDhf(t)l, (8)

Ьг- производная порядка 1/2. Положим/[/(01 = /М- Докажем аналог формулы (8). Для этого подействуем на обе части уравнения (5) оператором /

I дй д2й\

Подставим также значение [/(?)] = [/(£)]• В итоге имеем

~ , , д2и

]ОГ1и~] -Г1—=

Учитывая, что = Би, приходим к уравнению (2). Из определения оператора ] следует, что

и(0,х) = 0. Имеем также

ллао)] =л/«],

значит,

и{1, 0)] =f{t).

Итак, мы доказали, что решение задачи (2)-(4) находится по формуле

и(1;,х) =

Отсюда следует, что

<?« =лт=1о1/21т] =Б1'2[т]

где интеграл порядка 1/2 определен в теореме 7.

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

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

Список литературы

1. Романовский П. И. 1980. Ряды Фурье. Теория поля. Аналитические и специальные функции. Преобразования Лапласа. М.: Наука, 336.

2. Яремко О. Э., Яремко Н. Н. Обобщённое двойное преобразование Лапласа и его применения для решения уравнений в частных производных. Прикладная математика и Физика. 52(4): 239-245, 2020, DOI 10.18413/2687-0959-2020-52-4-239-245.

3. Brychkov Yu. A., Prudnikov А. P., Shishov V. S. 1979, 1981. Operational calculus. Itogi Nauki i Tekhn. Ser. Mat. Anal., 16, VINITI, Moscow, 99-148; J. Soviet Math., 15:6, pp. 733-765.

4. Ermolova N. Y., Tirkkonen O. 2014. Laplace Transform of Product of Generalized Marcum Q, Bessel I, and Power Functions With Applications. IEEE Transactions on Signal Processing IEEE Trans. Signal Process. Signal Processing, IEEE Transactions on.pp. 2938-2944 Jun.

5. Jeffreys H. and Jeffreys B. 1956. Methods of Mathematical Physics, 3rd ed., Cambridge Univ. Press.

6. Ganzha E. 1.2012. On Laplace and Dini transformations for multidimensional equations with a decomposable principal symbol. Programming and Computer Software. 38: 150-155.

7. Gonzalez-Acuna, Rafael G., Gutierrez-Vega, Julio C. 2019. Transition integral transform obtained from generalization of the Fourier transform. Ain Shams Engineering Journal; 10(4): 841-845.

8. Jarad Fahd, Abdeljawad Thabet. 2018. A modified Laplace transform for certain generalized fractional operators. Results in Nonlinear Analysis, 1(2): 88-98.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

9. Koepf Wolfram, Kim Insuk, Rathie Arjun K. 2019. On a New Class of Laplace-Type Integrals Involving Generalized Hypergeometric Functions. Axioms, 8(3): 87.

10. Meerschaert M. M., Mortensen J., Wheatcraft S. W. 2006. Fractional Vector Calculus for Fractional Advection-Dispersion, Physica A, 367: 181-190.

11. Mainardi F. 2010. Fractional Calculus and Waves in Linear Viscoelasticity . Imperial College Press: London, 347.

12. Miller K. S., Ross B. 1993. An Introduction to the Fractional Calculus and Fractional Differential Equations. John Wiley and Sons, Inc.: New York, NY, USA, 366.

13. Napalkov V. V., Mullabaeva A. U. 2015. On one class of differential operators and their application. Proceedings of the Steklov Institute of Mathematics, 288: 142-155.

14. Rashidi M. M. 2009. The modified differential transforms method for solving MHD boundary-layer equations Comput. Phys. Commun., 180: 2210-2217.

15. Samko S., Kilbas A., Marichev O. 1993. Fractional Integral and Derivative. Theory and Applications.— Switzerland: Gordon and Breach.

16. Shishkina E. L., Sitnik S. M. 2020. Transmutations, Singular and Fractional. Differential Equations with Applications to Mathematical Physics. Elsevier.

17. Tsarev S. P. 2005. Generalized Laplace Transformations and Integration of Hyperbolic Systems of Linear Partial Differential Equations. In: Labahn, G. (ed.) Proc. ISSAC 2005. ACM Press, 325-331.

18. Zaikina Svetlana M. 2014. Generalized Integral Laplace Transform and Its Application to Solving Some Integral Equations. Vestnik Samarskogo Gosudarstvennogo Tehnivceskogo Universiteta. Seria: Fiziko-Matematiceskie Nauki, 1(34): 19-24.

19. Sitnik Sergei M., Yaremko Oleg, Yaremko Natalia. 2020. Transmutation Operators and Applications. Transmutation Operators Boundary Value Problems, Springer Nature Switzerland, 447-466.

References

1. Romanovsky P. 1.1980. Fourier series. Field theory. Analytical and special functions. Laplace transformations. Moscow: Nauka, 336.

2. Yaremko O. E., Yaremko N. N. Generalized double Laplace transform and its applications for solving partial differential equations. Applied Mathematics and Physics. 52(4): 239-245, 2020, DOI 10.18413/2687-09592020-52-4-239-245.

3. Brychkov Yu. A., Prudnikov A. P., Shishov V. S. 1979, 1981. Operational calculus. Itogi Nauki i Tekhn. Ser. Mat. Anal., 16, VINITI, Moscow, 99-148; J. Soviet Math., 15:6 , pp. 733-765.

4. Ermolova N. Y., Tirkkonen O. 2014. Laplace Transform of Product of Generalized Marcum Q, Bessel I, and Power Functions With Applications. IEEE Transactions on Signal Processing IEEE Trans. Signal Process. Signal Processing, IEEE Transactions on.pp. 2938-2944 Jun.

5. Jeffreys H. and Jeffreys B. 1956. Methods of Mathematical Physics, 3rd ed., Cambridge Univ. Press.

6. Ganzha E. 1.2012. On Laplace and Dini transformations for multidimensional equations with a decomposable principal symbol. Programming and Computer Software. 38: 150-155.

7. Gonzalez-Acuna, Rafael G., Gutierrez-Vega, Julio C. 2019. Transition integral transform obtained from generalization of the Fourier transform. Ain Shams Engineering Journal; 10(4): 841-845.

8. Jarad Fahd, Abdeljawad Thabet. 2018. A modified Laplace transform for certain generalized fractional operators. Results in Nonlinear Analysis, 1(2): 88-98.

9. Koepf Wolfram, Kim Insuk, Rathie Arjun K. 2019. On a New Class of Laplace-Type Integrals Involving Generalized Hypergeometric Functions. Axioms, 8(3): 87.

10. Meerschaert M. M., Mortensen J., Wheatcraft S. W. 2006. Fractional Vector Calculus for Fractional Advection-Dispersion, Physica A, 367: 181-190.

11. Mainardi F. 2010. Fractional Calculus and Waves in Linear Viscoelasticity . Imperial College Press: London, 347.

12. Miller K. S., Ross B. 1993. An Introduction to the Fractional Calculus and Fractional Differential Equations. John Wiley and Sons, Inc.: New York, NY, USA, 366.

13. Napalkov V. V., Mullabaeva A. U. 2015. On one class of differential operators and their application. Proceedings of the Steklov Institute of Mathematics, 288: 142-155.

14. Rashidi M. M. 2009. The modified differential transforms method for solving MHD boundary-layer equations Comput. Phys. Commun., 180: 2210-2217.

15. Samko S., Kilbas A., Marichev O. 1993. Fractional Integral and Derivative. Theory and Applications-Switzerland: Gordon and Breach.

16. Shishkina E. L., Sitnik S. M. 2020. Transmutations, Singular and Fractional. Differential Equations with Applications to Mathematical Physics. Elsevier.

17. Tsarev S. P. 2005. Generalized Laplace Transformations and Integration of Hyperbolic Systems of Linear Partial Differential Equations. In: Labahn, G. (ed.) Proc. ISSAC 2005. ACM Press, 325-331.

18. Zaikina Svetlana M. 2014. Generalized Integral Laplace Transform and Its Application to Solving Some Integral Equations. Vestnik Samarskogo Gosudarstvennogo Tehnivceskogo Universiteta. Seria: Fiziko-Matematiceskie Nauki, 1(34): 19-24.

19. Sitnik Sergei M., Yaremko Oleg, Yaremko Natalia. 2020. Transmutation Operators and Applications. Transmutation Operators Boundary Value Problems, Springer Nature Switzerland, 447-466.

Получена 17.05.2021

Яремко Олег Эмануилович - доктор физико-математических наук, доцент, профессор, Московский государственный технический университет «Станкин»

http://orcid.org/0000-0003-4619-0527 Вадковский пер.,1, Москва,127055, Россия E-mail: yaremki@mail.ru

Яремко Наталия Николаевна - доктор педагогических наук, доцент, профессор, Пензенский государственный университет

© http://orcid.org/0000-0003-1491-624X ул. Красная, 40, г. Пенза, 440026, Россия E-mail: yaremki@mail.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.