Проекционные методы решения гиперсингулярных интегральных уравнений на фракталах Текст научной статьи по специальности «Математика»

УДК 517.392

И. В. Бойков, А. И. Бойкова, П. В. Айкашев

ПРОЕКЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ ГИПЕРСИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ НА ФРАКТАЛАХ1

Аннотация.

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

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

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

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

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

I. V. Boykov, A. I. Boykova, P. V. Aykashev

PROJECTION METHODS FOR SOLVING HYPERSINGULAR INTEGRAL EQUATIONS IN FRACTALS

1 Работа выполнена при финансовой поддержке РФФИ. Грант 16-01-00594.

Abstract.

Background. Approximate methods for solving hypersingular integral equations are an actively developing area of calculus amthematics. It relates to multiple applications of hypersingular equations in aerodynamics, electrodynamics, physics, and also to one circumstance - analytical solutions of hypersingular integral equations are possible only in exceptional cases. Besides direct applications in physics and engineering, hypersingular integral equations of first kind occur in approximate solution of boundary problems of mathematical physics. Recently, the interest to studying analytical and numerical methods for solving hypersingular integral equations has significantly increased due to active application of methods of fractal geometry in radio engineering and radiolocation. It has turned out that one of main methods of fractal antennas modeling is hypersingular integral equations. The present work suggests and substantiates spline-collocation methods of zero and first orders for solving hypersingular integral equations on fractals.

Materials and methods. The study used methods of functional analysis and approximation theory. The authors considered linear one- and two-dimensional hypersingular integral equations on fractals. For determinancy in the case of a one-dimensional integral the research took the Cantor perfect set as an area of integration, and in the case of two-dimensional one - the Sierpinski carpet. The authors built projection computing schemes, substantiated on the basis of the analysis of logarithmic norms of the corresponding matrixes.

Results. The authors built three computing schemes for solving hypersingular integral equations on fractals of various types and dimensionality and obtained estimates of rapidity of convergence and error of the said schemes. The generated schemes appear to be models for building and substantiation of computing schemes on fractals of various nature.

Conclusions. The authors built and substantiated computing schemes of approximate solution of hypersingular integral equations of first and second kinds on model fractals. Model fractals are represented as the Cantor perfect set and the Sier-pinski carpet. The obtained results may be used in fractal antennas modeling.

Key words: hypersingular integral equations, projection methods, fractals, fractal antennas.

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

Данная работа посвящена приближенным методам решения гиперсингулярных интегральных уравнений на фракталах.

Работа состоит из пяти разделов и построена следующим образом. В первом разделе приведен краткий обзор приближенных методов решения

гиперсингулярных интегральных уравнений. Во втором разделе приведены определения различных видов гиперсингулярных интегралов. Третий раздел содержит вспомогательные предложения. В нем дано определение логарифмической нормы матрицы и приведены важные теоремы, связывающие нормы обратных матриц с их логарифмическими нормами. Четвертый раздел посвящен приближенному решению гиперсингулярных интегральных уравнений на фрактальных кривых. В качестве конкретной кривой взято совершенное множество Кантора. В пятом разделе исследуются приближенные методы решения гиперсингулярных интегральных уравнений на фрактальных поверхностях. В качестве конкретной поверхности взят ковер Серпинского (который иногда называют [1] пылью Серпинского).

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

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

Для приближенного решения гиперсингулярных интегральных уравнений с фиксированной сингулярностью

предложены и обоснованы [3] приближенные методы коллокации и механических квадратур.

В работе [4] предложен и обоснован сплайн-коллокационный метод со сплайнами нулевого и первого порядков для решения гиперсингулярных интегральных уравнений вида

Сплайн-коллокационный метод нулевого порядка предложен и обоснован в работе [2] для приближенного решения бигиперсингулярных интегральных уравнений вида

1. Обзор приближенных методов решения гиперсингулярных интегральных уравнений

■с) *) + f ^^ * т = f (t)

J п-F

-1

1 1

+1 t2,ТьТ2)х(тьТ2)dТ^Т2 = /p = 2,3,...

-1-1

Приближенные методы решения гиперсингулярных интегральных уравнений теории крыла в одномерном случае исследовались в работе [5], для многомерных - в работах [6, 7].

Обзор приближенных методов решения гиперсингулярных интегральных уравнений первого рода представлен в работе [8].

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

Для приближенного решения гиперсингулярных интегральных уравнений

аЦ)х(0 + Ъ^)| Т + т)х(^т = /^),

у (т t) у

a (t) x(t) + f d т = f (t),

Y (т-t)p

?(t) x (t) + f

h(t, t) x( t)

(т-t)p

h(t, т, x( т)) У (т-t)p

d т = f (t),

где у - единичная окружность, р = 2,3,..., предложены и обоснованы вычислительные схемы методов коллокации и механических квадратур [11].

Приближенные методы решения гиперсингулярных интегральных уравнений на предфракталах начали исследоваться в последнее время [12].

2. Вспомогательные предложения

Пусть В - банахово пространство; К - линейный оператор, действующий из В в В; I - тождественный оператор. Логарифмическая норма Л(К) оператора К определяется [13] формулой

Л (К ) = Шп(|| I + Щ --1,

Н10

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

Пусть дана вещественная матрица А = (агу}, 7,] = 1,2,.,п, в п -мерном

пространстве Яп векторов х = (х1,...,хп) с нормой

lli = Z1 xk k=1

11X12 =

Z1 xk

k=1

1/2

x 3 = max | Xk

1<k<n

x

Логарифмическая норма матрицы A равна [14]: A) = max

Л2( A) = Va A + A

J

ajj + Z i aj

v i=U

(

Л 3( A) = max

aii + Z i aij j =1,j

T

где Xтах ((A + A ) / 2) - наибольшее собственное значение матрицы

T

(A + A )/2. В данной работе используется третья норма.

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

Пусть A и B - квадратные матрицы с комплексными коэффициентами. Пусть x = (хь...,хп), у = (уь...,уп ), £ = (^1,--Лп ), П = (%•••,Пп) - п -мерные векторы с комплексными компонентами. Пусть фиксирована норма векторов и подчиненная ей операторная норма матриц; этой норме сопоставлена логарифмическая норма.

Теорема 1 [15]. Если Л(A) < 0, то матрица A неособенная и

-1 < 1/1Л (A) |.

Теорема 2 [15]. Пусть Ax = By = ц, причем Л(A) < 0, Л(B)<0.

Тогда ||* - 4 <1^+1^

IЛ(В) | |Л(А)Л(В) | В работе используется следующий класс функций.

Определение 1. Обозначим через Жг(М,[а,Ь]) множество функций, определенных на сегменте [а, Ь], имеющих непрерывные частные производные до (г — 1) -го порядка, кусочно-непрерывную функцию г -го

порядка и удовлетворяющих неравенствам | ф(к)(х)|< М, к = 0,1,.., г. 3. Определения гиперсингулярных интегралов

Определение 2 [16]. Интеграл вида Г—А(х)йХ— при целом р и

а (ь—х) р+а

0< а<1 определяет величину («конечную часть») рассматриваемого интеграла как предел при х ^ Ь суммы

X А(№ + В(х)

а (ь—х)р+а (ь—х)р+а—11'

если предположить, что А(х) имеет р производных в окрестности точки Ь .

Здесь В( х) - любая функция, на которую налагаются два условия:

а) рассматриваемый предел существует;

б) В(х) имеет по крайней мере р производных в окрестности точки

х = Ь.

Произвольный выбор В(х) никак не влияет на значение получаемого предела: условие (а) определяет значения (р -1) первых производных от В(х) в точке Ь, так что произвольный добавочный член в числителе есть

бесконечно малая величина, по меньшей мере порядка (Ь - х)р . Определение 3. Интегралом

Ь ф(Т^т < < Ь

—-, а < с < Ь,

3 (т-с)р

а

в смысле главного значения Коши Адамара называется следующий предел:

9(x)d т =

= lim ( Т-с) p v^Q

b

f + f -Ф^ + Щ + ^1(v)lnv

(т-с)р J (т-с)р v

a y ' с+v y '

где Ъ(у), ^(у) - некоторые функции, выбранные так, чтобы указанный предел существовал. При с = а отсутствует первый интеграл в правой части предела; при с = Ь отсутствует второй интеграл. Рассмотрим интеграл

Mi

ф( Т1, %2)d ^d Т2

G ~

((Т1 - t1)2 + (Т2 -12)2)p/2'

где (tl,t2) - точка области G; р (р >2) - целое число.

Пусть Я^,е) - круг с центром в точке t = (tl,t2) и с радиусом е. Регуляризацией интеграла Ьф при р > 3 называется предел

¿ф = lim

e^Q

rf f -ф(%/2 - 5(e)-C( (

G\R(t,e) ((^1 - t1)2 + (T2 - t2)2)p/2 ep-2

где В(х), С(х) - любые функции, на которые налагаются следующие условия:

а) рассматриваемый предел существует;

б) В(х) имеет непрерывные производные до (р - 2) -го порядка в окрестности точки t;

в) функция С(х) имеет производные первого порядка в окрестности

нуля.

Достаточность этого определения доказана в [17].

4. Приближенное решение гиперсингулярных интегральных уравнений на фрактале Кантора

Классическое множество Кантора названо по имени Георга Кантора, который описал его в 1883 г. Множество Кантора строится на сегменте

Сд =[0,1]. Построение начинается с средней части (не включая концов) единичного отрезка. В результате удаляется открытый интервал (1 / 3, 2 / 3) и остается множество С = [0,1/ 3] и [2/3,1]. На следующем этапе выкидываем среднюю часть сегментов [0,1/3] и [2/3,1]. В результате получаем множество С2 =[0,1/9] и [2/9,1/3] и [2/3,7/9] и [8/9,1]. Этот процесс продолжается до бесконечности, и на каждом его этапе выкидываются средние части сегментов, составляющих множества С0, С1, С2, С3,... Предельное множество С, которое представляет собой пересечение множеств Сп, п = 0,1,.., называется классической пылью Кантора.

Нетрудно видеть, что на каждом п -м шаге построения пыли Кантора длины сегментов, составляющих множество Сп , равны (1/3)п.

Пусть п достаточно велико. Построим вычислительную схему приближенного решения гиперсингулярного интегрального уравнения:

а(')х(') + Ъ(') + [Ц(,т)х(%)с1% = /('),'е Сп, р = 2,3,... (1)

С (т- 0р С

п п

Приближенное решение уравнения (1) будем искать в виде

хп( 0= V %...,/„^„Л (0, (2)

h'-Jn

( 0 =

1, te А,- , ,

1'"'' n

0, te [0,1]\ А,1 , ,

1n

(3)

V'A

2 2 2 2 2 2 1 —i + —i2 H-----1--in; — i + —i2 H-----1--in H--

3 1 9 A 3^ 3 9 3^ 3n

^ i1, —, jn = 0,1'

Каждому сегменту А. поставим в соответствие узел

- 2 2 2 1 1 Ь ... = — Ц + — ь +-----1--\п +---.

г1 1п 3 1 9 2 3п п 2 3п

Коэффициенты а?1 ^ находятся из системы уравнений

а(Е ... )а,- ...1 + Ъ(Г ) V а ■ ■ ■ [ -—-+

г1 'п г1 1п г1 V Ь, ]2,-., 7п J (т-т ■ )р

• • ^■ А . . . ■ (т ^• • ^)

1п

+¿п Ъ Н (• • 1п , 0'1 • • ■ )% 72,. • -, ■п = /(\ • • % ), 1к =0,2, к = 1,.,п, (4) где суммирование V распространяется по всем сегментам, входящим в Сп, ¿п =3 .

Преобразуем систему (4) к следующему виду:

(sgnb(\ ■ • )

t; ■.) а ■.., + ь (ь ■. .j

l1 jn) j1 n \ j1 j

Zaf ■..j x

. j1 Jn

J1, h,--; jn

x

d t

A ■ . .j (^ tj1 ■ • 'jn ) Jl,^^■, Jn

J1 Jn

+ dn S. h (t1 ■ • % , tJ1 ■ • •Jn )ai1 ■ • %

= (sgnb( t1 ■. ,jn )) f (t1 ■. ,jn ), jl,-, jn =0,2.

(5)

Найдем условие однозначной разрешимости систем уравнений (4), (5). Представим систему (5) в матричной форме СХ = ^ и вычислим

логарифмическую норму матрицы С. Здесь X = {а,- . . },

1 п

^ = (8ВЛА^[.. п ))/(\,...,1п), 1 = 0,2.

Диагональные элементы матрицы С имеют вид У / \ / N - й Т

(sgnb(t1 ■ ■ ^ )) a (t1 ■ ■ ^ ) + b (\■ • ^ ) i

A ; i ( Т - t, ..i 1 n 1 n

+ dnh(t ■ .. , , t■

n"\\ ■ ■ -jn^ j1 ■ ■ -jn,

Нетрудно видеть, что г й т

1/2 1/3"

d т

Р -1

(2 • 3й)

n) p-1

\ ■ • % 1т-^1 • • 4п ) "1/21/3Т

Следовательно, диагональные элементы матрицы С равны

C„ = ■

b (t,

j1 ■ • i

p -1

(2 • 3n)p-1

+

+ (sgnb(^ ■ • -in )) I a ( ■ • -in ) + "3nh ( ■ • in , i ■ • \

(6)

Оценим сумму модулей вне диагональных элементов. Положим (7Ь...,¡п) = (0,...,0) и оценим

Iq■ . . Q=|b(tQ ■ . ,q)| S '

J1 ■ •' Jn

d т

A ,. .j (-tQ,...,Q)p ■'1 Jn

+ dn S

J1,-, Jn

(0 • -Q , j ■ • • Jn

где ^ означает суммирование по (у1,...,]п) ^ (0,...,0).

>1 ■ •' >п

Очень грубая, но достаточная оценка дает

'2 • 3

/

1 dT __ „ 1 ( 2 ^

0 . . . 0 <| Ъ(to . . ,о)| J ( < |Ъ(to . . .о)

2 (т- to . . .о)p p -1

3

V /

Отсюда следует, что сумма модулей любых внедиагональных

элементов .. j не превосходит

h ...г=\ъа ..■ )| У J -—-|+h<

j1 jn 1 v ji • • v ^ J (т-7 . )p 1

>1'-' j A . . . . ( \,..,jn)

1n

<|b(i . . , )|-±-2 1 n p -1

' 2 - 3" >p-1

3

V

+ H, H = max | h(t,т) |. (7)

Из выражений (6), (7) следует, что

Л (с)=^ (рЬ (2 •3")Р1 к \ • • -ъ)1 (-1+3Р-Т )+

+\•. ,и )) [а(^ • • ,и) + ^\• • *П ,\ • • ^ )) + Я

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

Теорема 3. Пусть а^), Ь(^), Ъ(),т) - ограниченные функции, | b(t) |>а > 0 при t е [0,1]. Тогда при достаточно больших п системы (4), (5) однозначно разрешимы.

Однозначная разрешимость систем уравнений (4), (5) еще не

гарантирует сходимость приближенного решения (4), (5) к точному решению

*

уравнения х ^).

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

Для простоты обозначений рассмотрим уравнение

а(0х(0 + 6(0 Г Х(Т)2 йт = /(0, t е Сп. (8)

С (т— )2

п

Мы ограничиваемся особенностью порядка р = 2, так как она наиболее часто встречается в приложениях.

Приближенное решение уравнения (8) будем искать в виде кусочно-непрерывной функции:

3n

,(t)=S S aij1-jnфц-л(t),te

C

1=1;1,-,;п

ф1 ^ (t)-

Г g¿,..,z(tхt e^,,.,

' 1' ' n 1' ' n

Q,te [Q,1]\A, -i .

1' ' n

(9)

(1Q)

Так как построение функций gl ^ . - (¿) на интервалах Д^ . -однотипно, то опишем построение функций g/ о... о(^) на интервале До . о. Функция gl о о (^) имеет вид

&1,Q,...,Q(t ) =

1, Q < t <■

,2n '

32nt - 3n +1 1

-^^ < t < JL__L

2 - 3n 32n 3n 32n

Q.JL -_!_ < t < jn;

3n 32n 3n

(11)

&2,0,...,0(t ) =

Q, 0 < t <•

2n

32nt -1 1

— < t < -1 -f

3n - 2 32n 3n 32n

1, -L-4- < t < -L.

3n 32n 3n

(12)

В качестве узлов коллокации возьмем узлы

1 12 11 V ■ = I- ■ +--■ ■ = I- ■ +----

г1,...,гп г1,...,гп 32^ г1,...,гп г1,...,гп 3п 32п

Приравнивая левые и правые части уравнения (8) в узлах ^ - ,

г1,.",гп'

порядка

получаем систему линейных алгебраических уравнений 2

n+1

a Г Сл ) ai + b Г Сл )i

xn(T)d т

Cn I т-t

l

l = 1,2, ij = 0,2, J = 0,1,...,n -1.

- f Г U) ■

(13)

Обоснование этой вычислительной схемы проводится по аналогии с обоснованием, приведенным в работе [4].

X

1

1

Теорема 4. Пусть выполнены следующие условия: 1) функции а(0, 6(0 ограничены при tе [-1,1]; 2) справедливо неравенство |а^) |> а> 0 при

* 2

t е [-1,1]; 3) уравнение (8) имеет единственное решение х (0 еЖ (М),

*

М = соп81. Тогда система уравнений (13) имеет единственное решение хп (О

и справедлива оценка

x (t)-xn(t)

< с -1

3"

5. Приближенное решение гиперсингулярных интегральных уравнений на фрактале «пыль Серпинского»

Напомним построение фрактала, который называется пылью Серпинского. Рассмотрим на плоскости OXY квадрат О = [0,1; 0,1]. Построение

фрактала заключается в последовательном выполнении следующих действий:

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

2. Множество О1 состоит из восьми квадратов со сторонами, параллельными координатным осям, и с длинами сторон равными 1/3. Пронумеруем эти квадраты начиная с квадрата, находящегося в левом нижнем углу квадрата О, двигаясь против часовой стрелки. Пронумерованные квадраты обозначим через А/ = 1,2,...,8. Каждый квадрат А/

■1 = 1,2,...,8, делится на девять равных квадратов и из каждого выбрасывается центральный квадрат. Оставшиеся в квадрате А/ /'1 = 1,2,...,8, более мелкие

квадраты пронумеруем, двигаясь против часовой стрелки начиная с левого нижнего квадрата, расположенного в Aj1 . Полученные квадраты обозначим

через А/^,/1,= 1,2,...,8. Полученное множество обозначим через О2 .

3. Описанная выше процедура проделывается с множеством О2, затем с О3 и т.д.

Замечание. Предфрактал n -го порядка обозначается через Оп.

Пересечение ^ О/ множеств О/, i = 1,2,..., называется пылью j=1

Серпинского.

Будем рассматривать гиперсингулярные интегральные уравнения на предфракталах п -го порядка. Квадраты, входящие, в предфрактал п-го

порядка, обозначим через А// / , где 1 < / < 8, j = 1,2,...,п.

12' п

Рассмотрим гиперсингулярное интегральное уравнение

a(t1,t2)x(t1,t2) + Ъ(t1,t2)J J-X(T12T2)dT1dT22 p = f (t1,t2),(t1,t2)еОп. (14)

О ((Т1 -11)2 + (Т2 -12)2)p

п

Приближенное решение уравнения (14) будем искать в виде кусочно-постоянной функции

Известия высших учебных заведений. Поволжский регион

хп ^Ъ Х2) = 2 а-1 ,-2 ,...-п ^,72,...7п (Х1, Х2 X ХЪ Х2 ,

;1,,2,...,,п

V,.../ (tl,t2) =

12 n

1, (tl,t2) e Aj1...j ,

1 n

0, (ti,t,)eQ\A: ...»

1 n

(15)

(16)

jj, i = 1,2, ...,8, J = 1,2,..., n.

Сумма в формуле (16) распространяется на индексы , ■ = 1,2,...,8, > = 1,2,..., п.

Обозначим ...-п пересечение диагоналей в квадрате Д^ - . Коэффициенты {а.^ - } определяются из решения системы

a(tj1...jn )%...»„ +

+Щ..»п ) S j,...j i i

d T,d т,

j1, j2,...,j

ji... j

j1...»n J И j1...jn

. T1 -/ I +1 T2 -1

V

f ((...in ), ij,» = 1,2,..,8, j = 1,2,...,n.

2 ^

(17)

-l

Здесь X. ■ , I = 1,2, - проекции точки Хг,./ на оси х1, I = 1,2.

■1... п 1 п

Предложенный метод справедлив для произвольных вещественных значений р. Ниже для определенности мы остановимся на случае, когда р = 3, так как сингулярность порядка р = 3 находит наибольшее применение в прикладных задачах.

Известно [7], что

Ьй

ac (( x - Xq)2+(y - Уо)2)3/2

( XQ - a)2 + (Уо - d)2 ^ (XQ - b)2 + (Уо - d)2

( X0 - a)(Уо - d) (XQ - b)(Уо - d)

( XQ- -b)2 + (уо - с)2 ^ l(xq - a)2 + (уо - с)2

+

(хо - Ь)(Уо - с) (хо - а)(Уо - с)

Диагональные элементы в системе уравнений (17) имеют вид

йхйу

л Л ,7

Д.

?(V • -n ) + b(T1 ■ -n ) i

»1 ■ • \

Г' - f , _ \2 ^

Т -t»1 ■ * -n ) + \ 2 - ^j'l ■ • -n

3/2

(18)

(19)

Из формулы (18) следует, что

dxdy

A, J. / ' 1 - * " "„ Г + (^2 - 7, . ГГ

JJ

(Л , ^2x3/2 (т1 +т2)

Г 1 А ч Лп 1 Г 1 ^ п = -25/23п, A* = 1 Г1 ^ 1 Г1

2 V 3 I 2 V 3

(20)

■ v v у v у I

Легко видеть, что при выполнении условия | Ь(^ т) |> а >0,

2

(t, т) е [0,1] , система уравнений (17) эквивалентна системе уравнений

(sgnb( j . • -п )

a (t/ .. ■ ) а/ / +

j1 . ■п) 1... п

+Ъ (\ . ■ -п ) У а. . ■ j J J-

d ^d T2

j1 . • ■ j Aj . • j [( - %. ,„ )+( - T*.. „п

2 I - \2 ^3/2

/

= sgn(b(th..4n ))f (th..4n ), ij;i = 1,2,...,8, j = 1,2,...,п. (21)

Система уравнений (21) в матричной форме имеет вид CX = F, причем структуры матрицы C, векторов X и F очевидны.

Из формул (19) и (20) следует, что диагональные элементы матрицы C

равны

(sgn b(t11 ,,,jn ))(a(tj1 .. . n) - 25/23n Ъ(tj1 .. jn)),/. : i = 0,2, j = 0,1,.,n -1.

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

Теорема 5. Пусть функции a(t,т), b(t,т) ограничены и | b(t,т) |> а >0

при (t, т) е [0,1]2. Тогда при достаточно больших значениях п системы уравнений (17) и (21) однозначно разрешимы.

Предложенные выше алгоритмы апробированы в работе [18].

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

1. Кроновер, Р. М. Фракталы и хаос в динамических системах / Р. М. Кроновер. -М. : Техносфера, 2006. - 488 с.

2. Boykov, I. V. An approximate solution of hypersingular integral equations / I. V. Boykov, E. S. Ventsel, A. I. Boykova // Applied Numerical Mathematics. - 2010. -Vol. 60 (6). - P. 607-628.

3. Бойков, И. В. Приближенное решение интегральных уравнений Фредгольма с интегралом в смысле главного значения Коши - Адамара / И. В. Бойков // Функциональный анализ и теория функций. - Вып. 7. - Казань : Изд-во КГУ, 1970. - С. 3-23.

4. Boykov, I. V. An approximate solution of nonlinear hypersingular integral equations / I. V. Boykov, E. S. Ventsel, V. A. Roudnev, A. I. Boykova // Applied Numerical Mathematics. - 2014. - Vol. 86, December. - P. 1-21.

5. Capobiano, M. R. On the numerical solution of a nonlinear integral equation of Prandtl's type / M. R. Capobiano, G. Criscuolo, P. Junghanns // Operator Theory: Advances and Applications. - 2005. - Vol. 160. -Р. 53-79.

6. Вайникко, Г. М. Численные методы в гиперсингулярных интегральных уравнениях и их приложения / Г. М. Вайникко, И. К. Лифанов, Л. Н. Полтавский. -М. : Янус-К, 2001. - 508 с.

7. Оселедец, И. В. Приближенное обращение матриц / И. В. Оселедец, Е. Е. Тыртышников // Журнал вычислительной математики и математической физики. - 2005. - Т. 45, № 2. - C. 315-326.

8. Бойков, И. В. Приближенное решение гиперсингулярных интегральных уравнений первого рода / И. В. Бойков, А. И. Бойкова, М. А. Сёмов // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. -2015. - № 3 (35). - С. 11-27.

9. Бойков, И. В. К приближенному решению сингулярных интегродифференци-альных уравнений. 1 [линейные уравнения] / И. В. Бойков, И. И. Жечев // Дифференциальные уравнения. - 1973. - Т. 9, № 8. - C. 1493-1502.

10. Бойков, И. В. К приближенному решению сингулярных интегродифференци-альных уравнений. 2 [нелинейные уравнения] / И. В. Бойков, И. И. Жечев // Дифференциальные уравнения. - 1975. - Т. 11, № 3. - C. 562-571.

11. Бойков, И . В . Приближенное решение линейных гиперсингулярных интегральных уравнений методом коллокаций / И. В. Бойков, Ю. Ф. Захарова, М. А. Сёмов, А. А. Есафьев // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. - № 3 (31). - С. 101-113.

12. Бойков, И. В. Приближенное решение гиперсингулярных интегральных уравнений на предфракталах / И. В. Бойков, А. И. Бойкова, П. В. Айкашев // Математическое и компьютерное моделирование естественно-научных и социальных проблем : сб. ст. IX Междунар. науч.-техн. конф. молодых специалистов, аспирантов и студентов. - Пенза : Изд-во ПГУ, 2015. - С. 49-58.

13. Далецкий, Ю. Л. Устойчивость решений дифференциальных уравнений в банаховом пространстве / Ю. Л. Далецкий, М. Г. Крейн. - М. : Наука, 1970. - 534 с.

14. Деккер, К. Устойчивость методов Рунге-Кутты для жестких нелинейных дифференциальных уравнений / К. Деккер, Я. Вервер. - М. : Мир, 1988.

15. Лозинский, С. М. Замечание о статье В. С. Годлевского / С. М. Лозинский // Журнал вычислительной математики и математической физики. - 1973. - № 2. -С. 457-459.

16. Адамар, Ж. Задача Коши для линейных уравнений с частными производными гиперболического типа / Ж. Адамар. - М. : Наука, 1978. - 351 с.

17. Бойков, И. В. Приближенные методы вычисления сингулярных и гиперсингулярных интегралов. Ч. 2. Гиперсингулярные интегралы / И. В. Бойков. - Пенза : Изд-во ПГУ, 2009. - 252 с.

18. Бойков, И. В. Приближенные методы решения гиперсингулярных интегральных уравнений на фракталах. Вычислительный эксперимент / И. В. Бойков, П. В. Айкашев // Математическое и компьютерное моделирование естественнонаучных и социальных проблем : сб. ст. X Междунар. науч.-техн. конф. - Пенза : Изд-во ПГУ, 2016.

References

1. Kronover R. M. Fraktaly i khaos v dinamicheskikh sistemakh [Fractals and chaos in dynamic systems]. Moscow: Tekhnosfera, 2006, 488 p.

2. Boykov I. V., Ventsel E. S., Boykova A. I. Applied Numerical Mathematics. 2010, vol. 60 (6), pp. 607-628.

3. Boykov I. V. Funktsional'nyy analiz i teoriya funktsiy [Functional analysis and functions theory]. Iss. 7. Kazan: Izd-vo KGU, 1970, pp. 3-23.

4. Boykov I. V., Ventsel E. S., Roudnev V. A., Boykova A. I. Applied Numerical Mathematics. 2014, vol. 86, December, pp. 1-21.

5. Capobiano M. R., Criscuolo G., Junghanns P. Operator Theory: Advances and Applications. 2005, vol. 160, pp. 53-79.

6. Vaynikko G. M., Lifanov I. K., Poltavskiy L. N. Chislennye metody v gipersingulyarnykh integral'nykh uravneniyakh i ikh prilozheniya [Numerical methods in hypersingular integral equations and applications thereof]. Moscow: Yanus-K, 2001, 508 p.

7. Oseledets I. V., Tyrtyshnikov E. E. Zhurnal vychisslitel'noy matematiki i matematicheskoy fiziki. [Journal of calculus mathematics and mathematical physics]. 2005, vol. 45, no. 2, pp. 315-326.

8. Boykov I. V., Boykova A. I., Semov M. A. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2015, no. 3 (35), pp. 11-27.

9. Boykov I. V., Zhechev I. I. Differentsial'nye uravneniya [Differential equations]. 1973, vol. 9, no. 8, pp. 1493-1502.

10. Boykov I. V., Zhechev I. I. Differentsial'nye uravneniya [Differential equations]. 1975, vol. 11, no. 3, pp. 562-571.

11. Boykov I. V., Zakharova Yu. F., Semov M. A., Esafev A. A. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko -matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2014, no. 3 (31), pp. 101-113.

12. Boykov I. V., Boykova A. I., Aykashev P. V. Matematicheskoe i komp'yuternoe modelirovanie estestvenno-nauchnykh i sotsial'nykh problem: sb. st. IX Mezhdunar. nauch.-tekhn. konf. molodykh spetsialistov, aspirantov i studentov [Mathematical and computer modeling of natural scientific and social problems: proceedings of IX International scientific and technical conference of young specialists, postgraduate and undergraduate students]. Penza: Izd-vo PGU, 2015, pp. 49-58.

13. Daletskiy Yu. L., Kreyn M. G. Ustoychivost' resheniy differentsial'nykh uravneniy v banakhovomprostranstve [Stability of differential equation solutions in Banach space]. Moscow: Nauka, 1970, 534 p.

14. Dekker K., Verver Ya. Ustoychivost' metodov Runge-Kutty dlya zhestkikh nelineynykh differentsial'nykh uravneniy [Stability of the Runge-Kutta methods for stiff non-linear differential equations]. Moscow: Mir, 1988.

15. Lozinskiy S. M. Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki [Journal of calculus mathematics and mathematical physics]. 1973, no. 2, pp. 457-459.

16. Adamar Zh. Zadacha Koshi dlya lineynykh uravneniy s chastnymi proizvodnymi giperbolicheskogo tipa [The Cauchy problem for linear equations with partial derivatives of hyperbolic type]. Moscow: Nauka, 1978, 351 p.

17. Boykov I. V. Priblizhennye metody vychisleniya singulyarnykh i gipersingulyarnykh integralov. Ch. 2. Gipersingulyarnye integraly [Approximate methods of singular and hypersingular integrals calculation: Part 2. Hypersingular integrals]. Penza: Izd-vo PGU, 2009, 252 p.

18. Boykov I. V., Aykashev P. V. Matematicheskoe i komp'yuternoe modelirovanie estestvenno-nauchnykh i sotsial'nykh problem: sb. st. XMezhdunar. nauch.-tekhn. konf.

[Mathematical and computer modeling of natural scientific and social problems: proceedings of X International scientific and technical conference]. Penza: Izd-vo PGU, 2016.

Бойков Илья Владимирович

доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

E-mail: boikov@pnzgu.ru

Бойкова Алла Ильинична кандидат физико-математических наук, доцент, кафедра высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 4Q)

E-mail: math@pnzgu.ru

Айкашев Павел Владимирович

аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 4Q)

E-mail: math@pnzgu.ru

Boykov Il'ya Vladimirovich Doctor of physical and mathematical sciences, professor, head of sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)

Boykova Alla Il'inichna Candidate of physical and mathematical sciences, associate professor, sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)

Aykashev Pavel Vladimirovich Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)

УДК 517.392 Бойков, И. В.

Проекционные методы решения гиперсингулярных интегральных уравнений на фракталах / И. В. Бойков, А. И. Бойкова, П. В. Айкашев // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2о16. - № 1 (37). - С. 71-86.