Использование дробного интегро-дифференцирования в уравнениях электродинамики материальных сред Текст научной статьи по специальности «Физика»

УДК 534.222:629.127.4

Алифанов Роман Николаевич

ФГБОУ ВПО «Дальневосточный государственный технический рыбохозяйственный университет»

Россия, Владивосток1 Доцент, кандидат технических наук E-Mail: gidra_518@mail.ru

Карпачев Александр Афанасьевич

Военный учебно-научный центр Военно-Морского флота «Военно-морская академия имени Адмирала Флота Советского Союза Н.Г. Кузнецова»

Филиал Владивосток Профессор, доктор технических наук E-Mail: K327065@yandex.ru

Стародубцев Павел Анатольевич

Военный учебно-научный центр Военно-Морского флота «Военно-морская академия имени Адмирала Флота Советского Союза Н.Г. Кузнецова»

Филиал Владивосток Профессор, доктор технических наук E-Mail: spa1958@mail.ru

Использование дробного интегро-дифференцирования в уравнениях электродинамики материальных сред

Аннотация: В статье кратко представлен анализ основных исторических вех формирования теории динамических процессов. Он показал, что в настоящее время явно ощущается недостаточность традиционных физических моделей, потому что полное описание процессов современной обработки сигналов и полей невозможно с помощью формул классической математики, полученных на основе представления сигналов в пространствах целочисленных мер и гладких функций. Для устранения данных недостатков и в связи с созданием Б. Мандельбротом общей концепции фракталов у многих ученых возникла мысль о применении их в области радиофизики и радиолокации. А использование идей масштабной инвариантности - «скейлинга» и разделов современного функционального анализа, которые связаны с теорией множеств, теорией дробной размерности, общей топологией, геометрической теорией меры и теорией динамических систем, дополнительно открыло большие потенциальные возможности и новые перспективы в обработке многомерных сигналов и в родственных научных и технических областях, таких как нелинейная гидроакустика. Для применения формул классической математики, полученных на основе представления сигналов в пространствах целочисленных мер и гладких функций, авторами статьи предложен математический аппарат использования дробного интегро-дифференцирования в уравнениях электродинамики материальных сред. В его основу при дробном исчислении положены интегро-дифференциальные операторы Римана-Лиувилля и Капуто. При решении задачи анализа свойств электромагнитного поля в диэлектрике необходимо учитывать влияние фрактальных свойств движения зарядов в диссипативной среде на создаваемое электромагнитное поле и порядок дробного полинома, который определяет не только фрактальные свойства движения зарядов, но и, в некотором смысле, характеризует фрактальность исследуемой диссипативной среды.

Ключевые слова: Фракталы; скейлинг; дробное интегро-дифференцирование.

Идентификационный номер статьи в журнале 55TVN114

690087, г. Владивосток Ул. Луговая 52В - 501

Roman Alifanov

The Far Eastern State Technical Fishery University”

Russia, Vladivostok E-Mail: gidra_518@mail.ru

Alexander Karpachev

Military educational centre of science of Navy fleet «the Naval academy of a name of Admiral of Fleet of Soviet Union of N.G.Kuznetsova »

Branch Vladivostok Russia, Vladivostok E-Mail: K327065 @yandex.ru

Paul Starodubtsev

Military educational centre of science of Navy fleet «the Naval academy of a name of Admiral of Fleet of Soviet Union of N.G.Kuznetsova »

Branch Vladivostok Russia, Vladivostok E-Mail: spa1958@mail.ru

Use of fractional ntegrated-differentiation in the equations of electrodynamics of material environments

Abstract: The paper briefly presents an analysis of the major historical milestones formation of the theory of dynamic processe. He showed that at the moment is clearly a failure of traditional models of physical, because a complete description of the processes of modern signal processing fields, and it is impossible with the aid of classical mathematics, obtained on the basis of representation of signals in spaces integral measures and smooth functions. To address these deficiencies and in connection with the creation of the general concept by B. Mandelbrot fractals many scientists had the idea of applying them in the field of radio and radar. And the use of the ideas of scale invariance - " scaling" and sections of modern functional analysis, are associated with set theory, the theory of fractional dimensionality, general topology, geometric measure theory and the theory of dynamical systems, further opened up great potential and new perspectives in the treatment of multidimensional signals and related scientific and technical areas, such as nonlinear hydroacoustics. For the application of formulas of classical mathematics of the obtained on the basis of representation of signals in spaces of integral -represented measures and smooth functions, the authors proposed the mathematical apparatus of the use of fractional integro- differential equations of electrodynamics in material media. It is based on the fractional calculus put integrodifferential operators Riemann-Liouville and Caputo. In solving the problem analysis of the properties of the electromagnetic field in the dielectric is necessary to consider the effect of fractal properties of the motion of charges in a dissipative medium to create an electromagnetic field, and my order fractional polynomial, which determines not only the fractal properties of the motion of the charges, but also, in some sense, fractal study characterizes the dissipative environment.

Keywords: Fractals; self similarity; fractional ntegrated-differentiation.

Identification number of article 55TVN114

Введение

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

Однако открытия второй половины ХХ-го века кардинально изменили наше представление о характере динамических процессов. В настоящее время явно ощущается недостаточность традиционных физических моделей, потому что полное описание процессов современной обработки сигналов и полей невозможно с помощью формул классической математики, полученных на основе представления сигналов в пространствах целочисленных мер и гладких функций [1-3].

В конце двадцатого века в связи с созданием Б. Мандельбротом общей концепции фракталов [1] возникла мысль о применении их в области радиофизики и радиолокации. Использование идей масштабной инвариантности - «скейлинга» и разделов современного функционального анализа, которые связаны с теорией множеств, теорией дробной размерности, общей топологией, геометрической теорией меры и теорией динамических систем, открывает большие потенциальные возможности и новые перспективы в обработке многомерных сигналов и в родственных научных и технических областях, таких как нелинейная гидроакустика.

Основная часть

В [4] представлен один из первых способов введения дробного интегро-дифференцирования в основные уравнения электродинамики материальных сред. В этой работе при дробном исчислении используются интегро-дифференциальные операторы Римана-Лиувилля и Капуто [4]. Оператор дробного интегро-дифференцирования Римана-Лиувилля порядка а е Я с началом в точке 8 функции у(у) представляется в следующем виде:

йп

Б^уО = 81^п^ — 8)—Б“~пу(0, п —1 <а < п, п е N, а> 0,

<Ап

пау(У) = 8'^п(—Г —у(в^ ^в, Г(х) — гамма — функцияЭйлера, а < 0,

Г(—а) 8 У — в\а+

Щу(У) = 0, а = 0.

Оператор Капуто (регуляризованная дробная производная) определяется

с помощью равенства

йп

80уу(У) = signn(t — 8)В8;п — у(У\ п — 1 < а < п, п е N, а > 0.

Связь между операторами Римана-Лиувилля и Капуто дается соотношением

п—1 ) ( 8)

8с“у(0 = п>(0 — Ет7^—V п—1 < а <п п еN, а > 0. (1)

Г ^ + 1 — а)

Если выполняется равенство уУ)(^) = 0, то операторы Римана-Лиувилля и Капуто тождественны. При целочисленном значении параметра а эти операторы также совпадают между собой и совпадают с обычными производными целого порядка.

В то же время в [5] вводится понятие d — оператора порядка s > 0 , действующим над

множеством степенных функций xq , для которых выполняются условия:

s, q, x є R, S, |q| — const < x> :

d~sx: xq — Г{q +1) xq—s, q £ N, q — s £ N~ (2)

r(q + 1 — s) . ()

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

В [5] получена система уравнений Максвелла в дробных производных с использованием оператора Капуто. Заменим в данных уравнениях оператор Капуто на d — оператор порядка s > 0:

rot E — —1 d st: B; div B — 0 ; (3)

т

rot H — j +1 d~st: D ; div D — p ; (4)

т

и материальные уравнения среды [6]:

B — mH ; D — ££0E,

где E - напряженность электрического поля, в/м ;

B - магнитная индукция, В^м2 {кг-с 2-А);

H - напряженность магнитного поля, а/м ;

2

м

D - электрическая индукция, Кл/

j - плотность тока, Кл/ м2 ;

р - плотность тока проводимости, Кл) м ;

Г - калибровочный коэффициент для дробной производной;

S0 - электрическая постоянная, ф/м;

S - относительная диэлектрическая проницаемость среды;

М) - магнитная постоянная, Гн/м;

Ц - относительная магнитная проницаемость среды.

Векторный потенциал A вводится стандартно:

В = rot А.

Подстановка данного выражения в уравнение (3), приводит к формулам:

1 d~st: (rot A) = -1 т т

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

і і f 1 I

rot E — — d st: (rot A) — — rot d st: A — rot — d st: A I; div (rot A) — 0

^ т )

/

rot

1

E + -d~st: A I — 0 v т )

Так как ротор от градиента любой скалярной функции р равен нулю, то выражение в скобках равно градиенту этой функции, т.е.

E + -d st: A 1 —Vp, или E — —1 d~st: A + Vp,

где р - скалярный потенциал. Далее в уравнениях (4) векторы Н и Б заменяются на

их выражения через векторы В и Е . В свою очередь, векторы В и Е заменяются на векторный и скалярный потенциалы:

rot

B

— j +1 d st: (єє0E) ; div (єє0E) — p ,

или

(

rot

rot A

mm

\

— j + - d st: т

ЄЄп

— d st: A + Vp т

(

; div

ЄЄп

— d st: A + Vp т

Выполняются следующие преобразования данных уравнений:

rot rot A — mm j — d-2s : A + :0ЄЄ

d—ss : Vp —

—p

mm j—ттр d:a+v

mmp d~ s: p

div

C 1 ^ 1 n

—d~st :A + Vp ——d~st\(div) A + divVp — — . т ) т ЄЄ г.

(5)

(6)

Далее, не ограничивая общности, принимается условие калибровочной инвариантности:

div A —

тт0ЄЄ0 J—s

d st : p I.

Использование данного условия для исключения скалярного потенциала р в уравнении (5) и, наоборот, для исключения векторного потенциала А в уравнении (6) позволяет получить следующие уравнения:

• 7-2а_. .

rot rot A — mm j-~—— d 2st: A + V (div A) ,

т

1 j-s —d t: т

f d st :pVdivVp —p , или — d-2st: p + divVp —p .

v т J єє^ т ЄЄq

Из теории поля известно, что верны следующие выражения: rot rot A — V (div) A — AA, и divVp — Ap .

С учетом первого выражения определяется уравнение для векторного потенциала: V (div A) — AA — mm j — d ~2 st: A + V (div A),

T

или

ДА ~у~~2 а 2Ч: А - - л j. (7)

(ст)

А с учетом второго выражения получают уравнение для скалярного потенциала:

Л8 л2 *^ • ™ Р

-----а !: р - —

(ст) еей

Ар-^а 2*! :р-— . (8)

где с - ]Д/Л8 - 292 792 458 м/с (скорость света в вакууме).

Уравнения (7) и (8) представляют собой уравнения с изменяющимся типом: при * -1 -гиперболический тип; при * -1/2 - параболический тип. Такие уравнения в литературе называют диффузионно-волновыми уравнениями. Их решение можно получить методом функции Грина [7].

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

а-2.,. и(х,,)-{С±.Щй--(С^.Р^!), (9)

/Л£ дх /Л8 ££о

где под функцией г/(х, О понимается А или ср. Уравнение (9) - линейное, и его частное решение можно представить в виде

и(х,!) - щ■ ехр(/£х). г(!),

где г(!) - неизвестная функция, Щ - комплексная амплитуда, к - компонент

волнового вектора в направлении X. Подставляя это частное решение в (9), получаем уравнение

а-2* : г(!) + о2г(!) --°-.^М--------------------------------1-, (10)

к 880 щ ехр(/кх)

с^

где о - ,— - безразмерная частота. у1Л£

Решение уравнения (10) ищется обычным способом в виде степенного ряда [7]. Частным решением уравнения (10) является функция:

г(!) - Е2* (-о2!2*),

п

^ У

где Е2* (у) - У----------- - функция Миттага-Леффлера.

П-0 г ^ ■ п+1)

Дробное дифференцирование порядка 2* функции г (!) приводит к выражению:

ъ (_______г^2 + 2 $\п ъ (_г^2\ п

а 2* : г(!) - а-2* : Е (-о2!2*) - а 2* : У ( ° ! ) - У ( ° ) а2* :!2*п -

у Г(2sn + 1) У Г(2зп + 1)

ад / 2\п г,/о г,™ і 14 +2$п-2$ ад / ^2\1 / ^2\п-\/2$(п-1) ,-2$ ад / „2^2$\п

^ (-0 ) Г (2sn +1) • Ґ (-о ) (-О ) Ґ Ґ 2-^ (-о Ї )

л.л = А---------г.,-,.. = ^, -0 А-

„=0 Г(25п +1) Г(25п +1 - 2s) П=о Г(2s(n -1) +1) Г(-2s +1) П=о Г(2sn +1)

>-25 >—5

=-----а2Е (-а2Г25) =-----------------------------а2г(г).

Г (-25 +1) 25 (-25) Г (-25)

После подстановки результата дробного дифференцирования в уравнение (10) получаем [7]:

/~25 а>2 р(х, /) 1

25Г(-25) k ££0 и0 ехр^кх)

Отсюда определяется выражение для плотности тока проводимости:

Г2* к 2ееа

р(х, і) = щ ехр(ікх)

'0

2sГ(-2s) о

Таким образом, определяется решение

и(х, ?) = и0 • ехр(/кх) • Е2х (-а 15).

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

блуждании (5 < 1/2) затухание имеет степенную асимптотику Е25(-?2)да/^/Г^25 +1), характерную для многих фрактальных систем [7,8].

На рисунке 1 в качестве примера показаны графики функции Е2 5 (-^ ) для различных значений параметра 25 . Если параметр 25 находится в интервале от 1 до 2, то по переменной t будем иметь периодическую функцию с частотой а . Если параметр 25 находится в интервале от 0 до 1, то функция становиться монотонно убывающей. Нетрудно заметить, что параметр 25 определяет скорость убывания функции.

1,5 1

0,5 О -0,5 -1 -1,5

Рис. 1. Графики функции Миттага-Леффлера Е25 (-/25) при различных значениях параметра

25 .• 1) 25 = 2; 2) 25 = 1,9; 3) 25 = 1,8

Рассмотрим предельные случаи. Пусть 2s = 2 (гиперболический случай). Тогда функция Миттага-Леффлера преобразуется в гиперболический косинус:

/ r-^Ln

ад л,п ад \ л/ y I

ВД = = hjy) ■

П=0 Г(2n +1) П=о 2n!

Следовательно, решение u( х, t) будет иметь вид:

u(x,t) = u0 ■ exp(ikx) • E(-®2t2) = u0 ' exp(ikx) ■ ch\yj-02t21 = u0 ■ exp(ikx) • ch(iot)

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

При 2s = 1 (параболический случай) функция Миттага-Леффлера преобразуется в экспоненту:

ад п ад 1 «п

^Су) = Ъг, ^ =Ъ~; = ехР^-У)■

П=0 Г (п +1) п!

Следовательно, решение u( x, t) будет иметь вид:

u(x, t) = u0 ■ exp(ikx) ■ E(—o2t) = u0 ■ exp(ikx) ■ exp \-o2t| = u0 ■ exp(ikx - o2t)

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

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

d--st: = 0. (11)

Ц£ дх

Уравнение (11) - линейное, и его частное решение также можно представить в виде

u( х, t) = u0 ■ exp(ikx) ■ z(t),

где z(t) - неизвестная функция, Uo - комплексная амплитуда, k - компонент

волнового вектора в направлении X. Подставляя это частное решение в (11), получаем уравнение

d~2s: z(t) + cy2z(t) = 0, (12)

ckr

где о = i— - безразмерная частота. у1Ц£

Решение уравнения (12) также ищется в виде степенного ряда [8]. В качестве частного решения этого уравнения выберем дробный полином порядка 2s степени m , как определено в [4]:

т

Р2./т(О = £а '^(й+1Н, ап = сотг, (п = 0,1,2,...,т) .

п

=0

При взятии дробной производной порядка 2s от данного дробного полинома получаем:

Л: Р2„т0) = !: £а„ • <2‘п1)-1 = ±а„ • !^: (*п1)-1 = £ап Г^+ °)/2“-1

п=0 п=0

= а0

Д!) г1 + ^ Г(2<п + 1))^2.п-1 = 0 +Ц а Г(2^(п + 2))^2.(п+!)-1

Г(0) £ Г (2лп) £ ”+1 1 (2^(п +1))

Здесь Г(0) = го, поэтому первый член ряда равен нулю. Выбираем полином бесконечной степени, т.е. т = го , и потребуем выполнения условия:

Г (2^(п + 2) )

а = аи+1—т-1, (п = 0,1,2,3,...,) .

п п+1 Г(2^(и +1)) ( , , , , ,)

Отсюда следуют равенства:

-9 1

! st: р1со() = p2sh(t), ап = г(2ф +1)) ’ (п = 0,1,2,3,...,) .

Теперь определяем функцию:

К2. ^ ) = 1-

^2. (п+1)-1 ад ^2.п-1

„=0 Г (2&(п +1)) „=! Г (2яп)

Нетрудно убедиться в том, что для данной функции верны следующие выражения:

Л: К.) = К.) и Л: К* (0) = 0X0),

где 0- константа.

В качестве частного решения уравнения (12) выбираем функцию

^) = К2.(0) ,

при этом значение константы 0 определяется из уравнения:

02К (0) + юХ (0) = 0.

Отсюда получаем, что вЪ = -юЮ = (/ю)2, или 0 = (/ю)1^.

Таким образом, определяется решение

и(х, t) = и0 • ехр(/кх) • К2х ((/ю)12 t).

На рисунке 2 в качестве примера показаны графики функции к2. (ю, t) = (/ю)1 К2. ((/ю)11) для различных значений параметра 2. при ю = 1.

Рис. 2. Графики функции к2. (ю, () при различных значениях параметра 2.: 1) 2. = 2; 2)

2. = 1,9; 3) 2. = 1,8

Отметим, что в нашем случае функция К2. (/) при предельном значении параметра

2. = 2 (гиперболический случай) преобразуется в гиперболический синус:

) = £-

= sh(t).

„=; Г ( 2n ) п=i (2n -1)!

Следовательно, решение u( x, t) будет иметь вид:

u(x, t) = u0 • exp(ikx) • K2 ((ia)11) = u0 • exp(ikx) • sh ((ia)11)

Решение этого вида также определяет плоскую монохроматическую волну.

В параболическом случае 2s = 1 и функция K s (t) преобразуется в экспоненту:

ад > п—1 ад > п—1 ад > п

K1(t) = Y-Vt = ----=Y — = exp(t).

1() Y Г(n) Y (n — 1)! Yn! P()

Следовательно, решение u( x, t) будет иметь вид:

u(x, t) = u0 • exp(ikx) • K ((ia)1st) = u0 • exp(ikx) • exp ((ia)1st ).

В заключение отметим следующее. Значение параметра . определяет не только фрактальные свойства движения зарядов, но и, в некотором смысле, характеризует фрактальность исследуемой диссипативной среды [9,10].

ЛИТЕРАТУРА

1. Стародубцев ПА., Мироненко М.В., Григорьев В.В. Aкустическая голография и возможность ее применения в низкочастотной гидроакустической томографии океанской среды / Статья.- Владивосток, Сб. стат. пробл. и метод. разраб. и эксплуат. вооружен. и воен. техн. ВМФ.-Владивосток, МО РФ,ТОВВМУ имени С.О.Макарова,-Вып.19.- 1998 -с.131-137.

2. Стародубцев ПА., Пичугин K.A., Василенко AM. Узкополосное параболическое при-ближение для моделирования характеристик звукового поля движущейся возмущенной об-ласти / Статья.- Владивосток, Мат. 46 всерос. межвуз. научно-технической конференции. -Владивосток, ТОВМИ имени С.О. Макарова, 2003.-с.182-185.

3. Стародубцев ПА. Получение аналитического и графического отображения передаточ-ной характеристики гидроакустического канала с использованием импульсных сигналов / Статья.- Улан-Удэ, Вестник БГУ. -Вып.2.-2003.- с.107-112.

4. Потапов A.A. Фракталы, скейлинг и дробные операторы в радиотехнике и

радиоэлектроники: современное развитие / Статья.-Москва, журнал

радиоэлектроники, №1, 2010.-с. 1-98.

5. Потапов A.A., Черных ВА. Дробное исчисление A3. Летникова, теория фракталов и скейлинг / Под ред. A.A. Потапова.- М.: Физматлит, 2009.- 820 с.

6. Стародубцев ПА. Взаимодействие упругих и электромагнитных волн в проводящей морской / Статья.- Владивосток, Сб.докл. рег. научно-практ. конф. «Море-2003».-Тихоокеанская инспекция Российского Морского Регистра судоходства, 2003. -с.100 - 104.

7. Чуриков ВА. Дополнительные главы анализа. Дробное интегрирование и дробное дифференцирование на основе d-оператора: учебное пособие / Томск: Изд-во Томского политехнического университета, 2010. - 118 с.

8. Стародубцев ПА., Пичугин K.A. Оптимизация процесса обработки гидролокационной информации методами нелинейной фильтрации/ Статья. -Москва, Наукоемкие технологии. - № 5.-Т.5,-2004.-с.45-50.

9. Потапов A.A. Фракталы, скейлинг и дробные операторы как основа новых методов об-работки информации и конструирования фрактальных радиосистем // Технология и кон-струирование в электронной аппаратуре. 2008. № 5(77). с. 3

- 19.

10. Стародубцев ПА., Карпачев A.A., Халаев Н.Л. Новые физические объяснения

решений уравнений электродинамики материальных сред/ Статья. -

Владивосток, Проблемы и методы разработки и эксплуатации вооружения и военной техники ВМФ: сб. научных трудов. - Владивосток: ВУНЦ ВМФ «ВMA» (филиал, г. Владивосток). - Вып.82, ч.2,-2013.-с.7-12.

Рецензент: Карасев Владимир Владимирович, Профессор кафедры Судовождение, к.т.н., «Дальневосточный государственный технический рыбохозяйственный университет», ФГБОУ ВПО «Дальрыбвтуз».

REFERENCES

1.

2.

3.

4.

5.

6.

7.

8.

9.

10.

Starodubtcev P.A., Mironenko M.V., Grigoriev V.V. Acoustic holography and the possibility of its application in low-frequency sonar imaging of the ocean environment / Article . - Vladivostok, Sat stat . probl . method . razrab . and operation . forearmed. and military . tehn. VMF. Vladivostok, Russian Defense Ministry , TOVVMU named Makarov , - Vyp.19 . - 1998 p.131 -137 .

Starodubtcev P.A., Pichugin K.A., Vasilenko A.M. Narrowband parabolic approximation for modeling the behavior of the sound field of a moving perturbed the field / Article . - Vladivostok, mate . 46 All-Russia. Intercollege . scientific and technical conference . Vladivostok , TOVVMU named after SO Makarova , 2003. -p.182 -185 .

Starodubtcev P.A. Preparation of analytical and graphical display peredatoch -tion characteristics of sonar channel using pulsed signals / Article . - Ulan-Ude, BGU Herald . Vyp.2. - 2003 . - P.107 -112 .

Potapov A.A. Fractals, scaling and fractional operators in the radio and electronics : the modern development / Art. - Moscow, electronics magazine , № 1, 2010.-s. 1-98 .

Potapov A.A., Chernykh V.A. Fractional Calculus A.V. Letnikova , the theory of fractals and scaling / Ed. A.A. Potapov. - Fizmatlit, Moscow, 2009 . - 820 .

Starodubtcev P.A. Interaction of elastic and electromagnetic waves in conducting marine / Article . - Vladivostok, Sb.dokl . reg . scientific and practical . conf. "Sea -2003 ." Pacific Inspectorate of the Russian Maritime Register of Shipping , 2003 . , p.100 - 104.

Tchurikov V.A. Additional chapters of analysis. Fractional integration and fractional differentiation on the basis of d- operator : a tutorial / Tomsk: Tomsk Polytechnic University in 2010 . - 118 .

Starodubtcev P.A., Pichugin K.A. Optimization of processing sonar information methods for nonlinear filtering / Article . - Moscow, High Tech . - № 5.- T.5 , - 2004.

- p.45- 50 .

Potapov A.A. Fractals, scaling and fractional operators as a basis for new methods of information processing and design of fractal radio / / Technology and concentration struirovanie in electronic equipment . 2008 . Number 5 (77). s. 3 - 19.

Starodubtcev P.A. Karpachev A.A., Halaew N.L. Physical explanation of solutions of equations of electrodynamics of material media / article . - Vladivostok, Problems and methods of development and operation of weapons and military equipment of the Navy : Sat scientific papers. - Vladvostok VUNTS Navy " WMA " (Branch , Vladivostok) . - Vyp.82 , p.2 - 2013. - p.7- 12.