Приближенное решение сингулярного интегрального уравнения с применением рядов Чебышева

Владикавказский математический журнал 2016, Том 18, Выпуск 4, С. 15-22

УДК 517.392

ПРИБЛИЖЕННОЕ РЕШЕНИЕ СИНГУЛЯРНОГО ИНТЕГРАЛЬНОГО УРАВНЕНИЯ С ПРИМЕНЕНИЕМ РЯДОВ ЧЕБЫШЕВА

3. В. Бесаева, Ш. С. Хубежты

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

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

1. Введение

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

Наиболее ранней работой в этом направлении является работа М. А. Лаврентьева [1], в которой одна практическая задача гидродинамики сводится к сингулярному интегральному уравнению и обосновывается определенный численный метод решения таких уравнений. Про эту работу Н. И. Мусхелишвили в своей книге [2, с. 352] пишет: «Дальнейшая разработка этого и аналогичных методов приближенного решения сингулярных интегральных уравнений является, как мне кажется, одной из важнейших очередных задач теории этих уравнений». После этого рядом исследователей были разработаны различные методы численного решения сингулярных интегральных уравнений, одним из которых является «метод дискретных особенностей», предложенный С. М. Белоцерковским и обоснованный его учеником К. И. Лифановым [3-4].

Указанный метод до сих пор является одним из актуальных методов в теории приближений. Отметим, что этот метод дает приближенные значения решения в конечном числе точек. Во многих случаях требуется получить аналитическое приближение решения, годное на всем отрезке. К этому типу методов принадлежат методы, связанные с многочленами Чебышева.

В настоящей работе впервые предлагается метод с применением рядов Чебышева [5, с. 322] для приближенного решения сингулярных уравнений на отрезке интегрировния. Суть метода заключается в том, что задача решения сингулярного интегрального уравнения, после замены плотности рядом Чебышева, сводится к решению системы линейных алгебраических уравнений относительно неизвестных коэффициентов данного решения. После нахождения коэффициентов разложения С\, С2,..., Сп, приближенное решение получается в аналитическом виде, что позволяет найти значения неизвестной функции во всех точках отрезка [-1; 1].

© 2016 Бесаева 3. В., Хубежты Ш. С.

2. Ряды Чебышева

Рассмотрим функцию /, определенную на отрезке [-1,1] и принимающую действительные либо комплексные значения. Предположим, что на этом отрезке ее можно разложить в ряд по многочленам Чебышева первого рода, т. е. существуют постоянные ао, ... такие, что

/

1=00

Известно, что многочлены Чебышева Т(х) = еоз^агссовх) (1 = 0,1,...) ортогональны на отрезке [—1,1] с весом (1 — ж2)" и

1

/ (1 - ж2)"1 Tjj(x)dx =

„2W Т2,ч л,. К * = 0,

' U, А; > 0. -i ^ 2

Умножая обе части равенства (1) на (1 — ж2) ж) и интегрируя по отрезку [—1,1], получаем следующие формулы для коэффициентов ряд a (1):

1

2 Г -I

ак = - (1-х2) 2 f(x)Tk(x)dx, к = 0,1,...

-1

Отсюда с помощью подстановки ж = cos t получаем равносильную, и часто более выгодную формулу

п

2

2 [

Q>k[f] = — /(cosí) cos ktdt, к = 0,1,...

n J

о

Величину [/ ] назовем к-.м коэффициентом Чебышева функции /, а ряд

те

ак[/] Тк

к=0

— рядом Чебышева функции f. Здесь символ Y1 определяется формулой

m i

у] a,j = -ai + аг+1 + ... + ат, rri^l. j=i

3. Вычислительная схема для решения сингулярного интегрального уравнения

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

1 1

I [ + 1 [ к

п ] ъ — х п ] -1 -1

где к(х, ¿), / (х) £ Нг (а) (г ^ 10 < а ^ 1) — заданные функции (Нг (а) — класс функций, имеющих непрерывные производные до порядка г — 1, а производная по рядка г удовлетворяет условию Гёльдера с показателем а). Индексу к = 1 уравнения (2) соответствует решение вида

¥*

где является достаточно гладкой функцией на отрезке [—1,1]. Тогда уравнение (2) примет следующий вид:

1 1 Ку=- [ , 1 + - / —¡2=  к

-1 -1

Заменим неизвестную функцию ее разложением в ряд Чебышева:

те

= Е ^

к=1

где С& — неопределенные коэффициенты. Подставляя полученное разложение в уравнение (3), получаем следующее уравнение:

1 1 1 1 / \

п

_i

k=1

1

1 f k

7ГУ VT^Í5 -1

те \

]TCk Tk-i

k=1

или

те 1 1 1 Tk_im , ~ 1 1 k

^ ТгУ ^ 7гУ ^ ^ ^

к=1 _1 к=1 _1

Используя формулу обращения

1 (

7гУ л/Г^5 ¿-ж к

и квадратурную формулу Мелера

Р)

г—1 ^ '

1 1 1 1 m J 2k —1

- / / ¿ dt » —

где Uk_2 (ж) — многочлен Чебышева 2-го рода

sin

Uk_2(ж) =

получаем

где

к=2 к=1 j=l

Ж-,' = С08 —7Г, ] = 1,2,...,т.

2^-1, 2 т

Уравнение

1 -1

где Со — произвольная постоянная, определяющая <(Ъ). Таким образом, единственное решение <(Ъ) зависит от Со- Подставляя в (9) в место <(х) его представление (4), будем иметь

к=1

-1

х

ния Хг = — 1 + г/г, /г = г = 1,    ,п — 1, при т = п, получим систему линейных

алгебраических уравнений относительно неизвестных С1, С2,..., Сп,...

те те п

Е Скик_2(хг) + Е Ск  £ Е А;(жг,ж7)Тй_1(ж7) = /(Хг), I = 1,2,...,п - 1,

к=2 к=1 j=l

те п

к=1 j=l

В дальнейшем будем рассматривать приближенное значение функции <(Ъ) в виде

п

<(Ъ) « ^ Ск Тк-1 (ъ). к=1

Тогда вместо (11) будем иметь систему линейных алгебраических уравнений порядка п х п:

п п п

Е Скик_2{хг) + Е Сй • £ Е Кхг,Щ)Тк-г{х]) = /

к=1 5 = 1

Решение системы (12) относительно неизвестных С1, С2,..., Сп дает нам приближенное решение уравнения (2) в виде

^ — ъ2 к=

4. Обоснование вычислительной схемы и оценка погрешности

Пусть X — пространство функций ^>(i) G Hr(а). Норма в пространстве X определяется формулой

= lW)llc-.„ + *Ф '"у-у.

0<в<а

Через Хп обозначим пространство функций многочленов (рп(Ь) = ^ а^¿к с нормой

к=0

С[_1;1] + sup ---—-.

0<в<а

Будем считать, что оператор К действует го пространства X в X и имеет ограниченный обратный оператор К-1. Обозначим через Рп оператор, проектирующий пространство X на пространство Хп по формуле

Известно [9, с. 342], что для константы Лебега справедлива оценка ||Рп|| ^ С 1пп. Приближенное решение уравнения (3) будем искать в виде функции

п

= Ск Тк_1 (¿). к=1

Ск

представленной в операторной форме уравнением

p

1 n

l

l

11 УМ dt+l-i 1 -1 -1

Pn [k(x,%n(i)] di

= Pn[f

где Рп[^(£)] — оператор проектирования на множестве интерполяционных полиномов п

мулами (6) и (7) уравнение (13) можно переписать в виде

J+=pn[f\x)\. (14) -1 - -1 -

Применим метод коллокации к уравнению (3). В операторной форме он имеет вид

K-n^n —

1 Vn{t)

't-x

di + P

vr J -1

k(x, i) ipn(i) di

= PXf

ii

n

1

1

1

1

1

1

Оценим норму разности

К(рп - Кп(рп

VT J VT^ -1

k(x,t) ^„(t) dt - P

vr J VT^ -1

k(x, t) ^>n(t) dt

+

W

-1

Px

(k(x,t) - (x,t))(t) dt

-1

(k(x,t) - (x,t))(t) dt

= Ii + /2,

где k„(x,t) — полином наилучшего равномерного приближения степени n — 1 по переменной x к функции k(x,t). Нетрудно видеть, что

h ^Сп^.ЖхШ^пЬ Ь^СпЦРг^-ШХШФА.

И, следовательно,

К(рп - Kra¥>J < СпРЦРпЦЕ^Лк^Щф.

где Еп(к(х,Ъ)) = Еп(к(х,Ъ)).

Из последнего неравенства и общей теории приближенных методов для обратимых операторов [10, с. 517] следует, что при п таких, что д = Сп13\\К~1\\\\Рп\\Еп_1(к(х,1)) < 1,

существует обратный оператор Кп1. Оценим норму разности \\Кп — Кп\\. Очевидно

K-n^n Кn^n I

PX

1

1

VT J

(k(x,t) - Pn[k(x,t)]) ^n(t) dt

-1

^ C^WPnWElikix^t^WipnW,

где En = max En[k(x,t)]. Пусть существует такое По, что при n ^ По выполняется неравенство

n\\^n(k(x, t)) +

Cn'3 (||Pra||<(fc(a:,i)) + \\Рп\\Е1(к{х,1))) < 1. Из теоремы Банаха [10, с. 211] следует, что

- р* || < Спе (\\Pn\\X(k(x,t)) + \\Pn\\El{k{x,t))) ,

где и ^n — решения уравнений

Теорема. Пусть оператор К имеет обратный, функции k(x, t), f (t) G Hr (a) и n такой,

что

Cnß\\K-l\\ (Wn(k{x,t)) +Etn(k(x,t))^j Inn < 1. Тогда система (13) имеет единственное решение и

f - ^nII < O

n

ln n

r+a—в I

1

1

1

1

1

1

1

1

1

1

1

1

1

5. Примеры

Приведем несколько простых примеров, иллюстрирующих изложенный выше метод. Пример 1. Рассмотрим уравнение

1 1 - / , 1 • ей + - / ; 1 • (ж + г) М) (И = х -1 -1

с дополнительным условием

1

- I , 1 •<р(г)<и = 1, -1

т. е. Со = 1. Точным решением этого уравнения является функция = 1. В этом

случае решая систему (12) при п = 3 получаем: С1 = 1 С2 = 1, 324547  10_8, С3 = —3, 973643  10_8 и приближенным решением будет функция

^з (Ь) = С1То(0 + С2Т1(*) + СзТ2(Ь) = 1 + (1, 324547  10_8 )Ь — (3, 973643  10_8 )(2Ь2 — 1) » 1.

Пример 2. Рассмотрим уравнение

1 1

1 [ 1 ¥>(*), 1 / 1 , / \ 7 3

- / ,  ей + - / ,  ж + £ (¿) ей = -

-1 -1

с дополнительным условием

1

11

- / ¥>(*) <Й = О,

ТГУ ^Г"^ ' -1

т. е. С0 = 0. Точное решение = ¿.При п = 3 получаем: С1 =0 С2 = 1, С3 = 0, т. е.

^з (Ь) = ЭД  0 + Т1 (Ь)  1 + Т2 (Ь)  0 = Ь.

Пример 3. Рассмотрим уравнение

1 1

1 [ 1 ¥>(*), 1 / 1 , / \ 7 3

- / , • ей + - / , • ж + ^ т м = -х

тгУ ^г^ ¿-ж тгУ уТ^Т5 ^ ; ^ 2 -1 -1

с дополнительным условием

1

1 1 1 -1

т. е. Со = т?. Точное решение </?(£) = ¿2. При п = 3 получаем: С1 = 0,4999999, С2 = 6, 698393  10_^ Сз = 0, 5000001, т. е.

¥>з(*) = С1  1 + С2Ь + Сз (2Ь2 — 1) = 0,4999999 + (6, 698393  10_8) Ь + 0, 5000001(2Ь2 — 1) » Ь2.

Вычисления проводились на языке QBasic с использованием метода Гаусса для решения систем линейных алгебраических уравнений. Очевидно, что погрешность для всех |e(i)| ^ 0, 2  10-6. Эти численные результаты показывают, что выше изложенная вычислительная схема реализуется хорошо.

Литература

1. Лаврентьев М. А. О построении потока, обтекающего дугу заданной формы // Тр. ЦАГИ.—1932.— Вып. 118.—С. 3-56.

2. Мусхелишвили Н. И. Сингулярные интегральные уравнения.—М.: Наука, 1968.—511 с.

3. Велоцерковский С. М., Лифанов И. К. Численные методы в сингулярных интегральных уравнениях.—М.: Наука, 1985.—252 с.

4. Лифанов И. К. Метод сингулярных интегральных уравнений и численный эксперимент.—М.: Янус, 1995.-520 с.

5. Пашковский С. Вычислительные применения многочленов и рядов Чебышева.—М.: Наука, 1983.— 384 с.

6. Крылов В. И. Приближенное вычисление интегралов.—М.: Наука, 1967.—500 с.

7. Хубежты Ш. С. Квадратурные формулы для сингулярных интегралов и некоторые их применения.—Владикавказ: ЮМ II ВНЦ РАН, 2011.—232 с.

8. Суетин П. К. Классические ортогональные многочлены.—М.: Наука, 1979.—406 с.

9. Натансон И. 17. Конструктивная теория функции.—M.-JL: ГИФМЛ, 1949.—688 с.

10. Канторович Л. В., Акилов Г. 17. Функциональный У НУЛЮ. —М.: Наука, 1984.—750 с.

Статья поступила 30 марта 2015 г. Бесаева Зарина Вячеславовна

Юго-Осетинский государственный университет им. А. А. Тибилова, преподаватель кафедры математики ЮЖНАЯ ОСЕТИЯ, 100001, г. Цхинвал, Московская, 8 E-mail: besaeva85@mail.ru

Хубежты Шалва Соломонович

Северо-Осетинский государственный университет им. К. Л. Хетагурова,

профессор кафедры математического анализа

РОССИЯ, 362025, Владикавказ, ул. Ватутина, 40;

Южный математический институт — филиал ВНЦ РАН,

ведущий научный сотрудник отдела математического моделирования

РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22

E-mail: shalva57@rambler. ru

APPROXIMATE SOLUTION OF A SINGULAR INTEGRAL EQUATION USING CHEBYSHEV SERIES

Besaeva Z. V., Khubezhty S. S.

A new method of approximate solution of singular integral equations with application of Chebyshev series is offered. Decomposition coefficients are determined by means of the solution of systems of simple algebraic equations. A justification of the constructed computational scheme is given and error estimate is deduced. The method is illustrated by test examples.

Key words: singular integral, Chebyshev series, approximation of integral, error estimation, orthogonal polynomial, coefficient of expansion.