CC BY
30
Владикавказский математический журнал 2011, Том 13, Выпуск 2, С. 56-62
УДК 517.956
О КВАДРАТУРНЫХ ФОРМУЛАХ ДЛЯ СИНГУЛЯРНЫХ ИНТЕГРАЛОВ С ВЕСОВЫМИ ФУНКЦИЯМИ
Ш. С. Хубежты, Л. Ю. Плиева, З. В. Бесаева
Рассматриваются сингулярные интегралы с весовыми функциями. Строятся квадратурные формулы для сингулярных интегралов. Доказываются новые формулы обращения, аналогичные формулам для сингулярных интегралов с Чебышевскими многочленами. Указываются применения построенных квадратурных формул к численному решению сингулярных интегральных уравнений.
Ключевые слова: сингулярный интеграл, формулы обращения, квадратурные формулы. Рассмотрим сингулярный интеграл с весовой функцией
1
в{(р,х) = - ж €(-1,1),
п 7 £ — х
(1)
-1
где — плотность, удовлетворяющая условию Гельдера [4], а р(£) = (1 — £)а(1 + (а, в > —1) — весовая функция. Как известно (см. [1, 2, 5]), многочлены Якоби
р (а,в) =
(_1)« Я
{ > ■(! — х)~а{1 + ж)-/3
2пп!
(1 — х)а+п(1 + х)в+п
(2)
являются ортогональными с весом р(х) = (1 — х)а(1 + х)в на отрезке [—1,1]. Кроме сингулярного интеграла (1) рассмотрим несобственный интеграл
1
/(/) = ^ J р(х)/(х)(1х.
(3)
1
Для таких интегралов построены квадратурные формулы (см. [1,2]) наивысшей алгебраической степени точности вида
1 Л
(4)
1
- / рШ(х) = ^ АкПхк) + КП(Л П11 к=1
где в роли узлов хк (к = 1, 2,..., п) используются нули многочлена Якоби, коэффициенты Лк (к = 1, 2,... ,п) вычисляются по формуле
Ак = ~1 р(%)
р (а,в)
(х)
(х — хк )РП(а,в) (хк)
Ях.
2011 Хубежты Ш. С., Плиева Л. Ю., Бесаева З. В.
1
В классических работах [3, 4] для р(х) = (1—х)а(1+х)в рассматриваются случаи (а = —тр /3 = (а = — /3 = —(а = 0, /3 = 0), (а = /3 = — Они встречаются чаще всего на практике. Для последнего случая а = ¡3 = — \ [1] получена квадратурная формула
1
_1 k=1
где Rn(f) = 22n(12ra)!/(2") (V), -1 < V < 1; = cos J^j (A; = 1, 2,... , n) — корни ортогонального многочлена;
sm arccos ж
W 22ra(n!)2
Для случая а = — Ь, /3 = Ь, т. е. для интеграла
1
^/ут^/(ж)сгж' (6) -1
таких формул нет, но рекомендована подстановка [2] х = — после чего интеграл сводится к виду
-¡\ Г—^Я-х)йх, (7)
П У V 1 + х -1
который можно вычислять по формуле (5).
Далее мы решаем задачу: построить для интеграла (6) квадратурные формулы наивысшей алгебраической степени точности.
Сперва построим ортогональные многочлены на отрезке [—1,1] с весовой функцией
р(х) = • Этот многочлен обозначим через Рп °'5,0'5') (ж), тогда для любого многочлена Q(x) степени меньше п должно быть выполнено условие ортогональности 1 -- 1
-1 -1 Рассмотрим многочлен
5(х) = (1 + х)РПГ°'5;°'5) (х). Степень его равна (п +1) и он ортогонален на [—1,1] с весом 2 всякому многочлену С^(х) степени меньше п. Если его разложить по многочленам Чебышева первого рода Т&(х) (к = 1, 2,..., п), то в разложении, ввиду указанной ортогональности, должны обратиться в нуль все коэффициенты при многочленах (х) до степени п—1 включительно и разложение должно иметь форму 5(х) = СпТп(х)+Сп+1Тп+1 (х). Кроме того, так как 5(х) должно нацело делиться на 1+х, при х = —1 должно быть 5(-1) = СгаТга(—1) + Сп+1Тп+ (—1) = (—1)п(Сп — Сп+1) = 0. Поэтому Сп+1 = Сп и
р(-0,5;0,5)^ _ £< Тп+\{х) + Тп{ж)
п
<Т I г
п
ж+ 1
Из свойств старших коэффициентов ортогональных многочленов имеем
_ (2п)!
^П -
22n(n!)2
2
1
Если положить x = cos $ и воспользоваться тем, что Tk(x) = cos(k arccos x) = cos k$, то для изучаемого полинома Якоби получаем
р(-0,5;0)5)Ы = JMLCOS^?
1 J 22ra(n!)2 cos § ' 1 J
„„„ 1?
^ V"".
Корни его есть
2k-1 ,, . 0 ,
Хь = cos-7Г (k = 1,2,..., n).
2n +1 v ' ' ' ;
Используя многочлены (8), получим следующую квадратурную формулу наивысшей
алгебраической степени точности [10]
где
_1 k=1
Введем обозначения
Sn(x) = , сп(х) =-, (10)
Sin ^ COS ^
где $ = arccos ж. Ниже приводятся конкретные многочлены Со(ж) = 1; Ci(x) = 2ж - 1; C2 (ж) = 4ж2 - 2ж - 1; С3(ж) = 8ж3 - 4ж2 - 4ж + 1; C4 (ж) = 16ж4 - 8ж2 - 12ж2 +4ж + 1; Сб(ж) = 32ж5 - 16ж4 - 32ж3 + 12ж2 + 6ж - 1;
50 (ж) = 1;
51 (ж) = 2ж + 1;
52 (ж) = 4ж2 + 2ж - 1;
53 (ж) = 8ж3 + 4ж2 - 4ж - 1;
54 (ж) = 16ж4 + 8ж2 - 12ж2 - 4ж + 1;
55 (ж) = 32ж5 + 16ж4 - 32ж3 - 12ж2 + 6ж + 1;
В классических работах (см. [3, 4, 7]) широко используются следующие формулы обращения:
1
i f—L=Mldt = Un-1(x), и = 1,2,..., 7Г J \/1 — t2 t — X
-1 (11)
- [ Vl^^^dt = -Tn+1(x), n = 0,1,..., п J t — x
-1
где Tn(x) = cos(n arccos ж), Un(x) = sm(("-+|^r^cosa:) _ ортогональные многочлены Чебы-шева I-го и II-го рода [5].
Нами доказана следующая
Теорема. Для многочленов 5п(х) и Сп(х) справедливы следующие формулы обращения
i
i Í + = (12)
п J VI — 11 — x -i
i
- Í \ ]J,S+n}f) dt =-Cn(x) (-1<ж<1, n = 0,1,2,...).
п j y i i tt x -i
< Применим упомянутое выше представление
_ Tn(x) + Tra+i(ж)
Сга(ж) - ГТж '
получаем
i
1 Л/1 + í 1 Гп(*) + Tn+1 (t) dt = 1 Г 1 rra(¿) + Тга+1 (¿) ^
7Г J yl-tt-X l+t IT J yjl -12 í-ж
ii
Отсюда, используя формулы (11), имеем
i f ¡i + tcn(t) dt = Un_l{x) + Un{x)j
п J Vi — t t — x -i
из которой после некоторого упрощения с учетом x = cos получаем
1 } fl + t Cn(t) _ sinnfl + sin(n + l)fl _ sin (g^ arceos _
— / A /--^^ — -, :- — -T~¡-7- — O (X ).
IT J\l-tt-x VI- COS21? sin (i arceos x)
Аналогично из представления
"" 1 — Ж '
c учетом (11) получаем
i ,- i
i f LZ^lд = I í iTn(t)-T^t) = _
TTjyi + tt-x тгУ VT^í5 t-x ii
sinn#-sin(n + l)# COS^^-1? sin # eos |
= —Cn(x). >
Далее, для сингулярных интегралов с весовыми функциями р(х) = и р(%) =
^построены следующие квадратурные формулы:
On -I- 1
nj Vi — tt — x 2n + 1^ x — xk /1o4
-i k=i (13)
2k — 1
X COS 2 + ^^ 7г[Сп(ж) +/Xfc(a;)](/?(a;fc) + Rn((p;x),
i
где
цк{х) = У] Aa Cn^ - АкС'п{хк)
' гр _ гр
<Aj <AJ @
о=1
2k - 1 4 2 2k - 1
хк = cos-7Г, Ак =-cos —--7г;
2n + 1 2n + 1 2(2n + 1) '
1 / fT^trt) 2 " (-i)yi-^
9n 4- 1
nj V1 + tt - X 2n + 1 ^ x - Xfc -1 k=1 (14)
X sin ^ 1 [ - ¿'„(ж) + Дй(ж)] + /?„(</?; ж),
где
Цк(х) = jTAa-?^-AkS'n( хк),
а=1
2kn . 4 2 kn ж& = cos ---, Ак = --- sin
2п + 1' 2п + 1 2п + 1
Дадим оценку погрешности полученных квадратурных формул.
Теорема. Если плотность имеет производные до г-го порядка включительно (г > 1) и (¿) удовлетворяет условию Гёльдера Н(а) (0 < а ^ 1), то для квадратурных формул (13) и (14) справедлива оценка
\3(<р,х)-3п(<р,х)\^0' 1
n
r+a— 1
равномерная для всех -1 ^ x ^ 1, где x) — соответствующая квадратурная сумма. < Пусть ^>(t) е Hr(а). Тогда
)=++• • •+ы
t
+ ¥>(r) Ы + ^ J(t- uy-1 (u) - [xk)) du.
Xk
Отсюда для остатка квадратурных формул (13) и (14) получаем
п Хк+1 *
к=0 Хк Хк
= -1, хп+1 = 1 и I</?М(и) - </?(г)Ы < - ж^|а. А это при р(Ь) = 11=1
где Жо = —I, жга+1 = L и \cpv > [и) — cpv >{хк)\ ^ A\u — хк\ . А это при р[Ъ) = или
P(t) = I Дает оценку
1
nr
< О { —г ) -
Формулы (5), (9), (13) и (14) могут быть успешно применены к численному решению задач плоской теории упругости. Например, задача о трещинах сводится к интегральному уравнению вида [6]:
1 1
1 [ + - [ K(x,t)<p(t)dt = f(x), же (-1,1), (15)
J t — x п J
-1 -1
где неизвестная функция ^>(t) ищется в виде
V® = V® = УТТ^0^'
или
В первом случае сингулярное интегральное уравнения (15) примет вид 1 _ 1
п J VI — tt — x п J VI — t 11
1 [ ¡l + tip0{t) dt+ 1 Г =/(ж), же (-1,1). (16)
Используя квадратурные формулы (9) и (13) для (16), получим дискретное уравнение 1 / 2к — 1
cos (2(2п + 1)^) [с«(ж)
k=1 k (17)
^TlEcos2 (^^к)к(х,хк)Мхк) = Пх).
Придавая параметру ж последовательно значения Ж1,Жг,... ,ж„, ж & = cos 7Г, получим систему линейных алгебраических уравнений порядка n х n относительно неизвестных (x1), ^o(x2), . . . , (x„).
Кроме решения уравнения (15) обычно требуется еще вычислить компоненты напряжений и смещений. Они выражаются с помощью интегралов типа Коши
1
1
1
Следовательно нужно вычислить интегралы типа Коши с весовыми функциями = Р&) = Для указанных интегралов от ортогональных многочленов
5п(х) и Сп(х) получены следующие формулы [9]
1
1 г ii + tcn{t) i(z-V^T)n
^>-dt = «1*- У*'-1Г(! + , _ ^¡JTT),
m J VI - tt-z V-г2 - 1
1
1
i r [T^tsn(t) i{z-jznT)n
™ di = + г + ^¡JTT),
vri у V 1 + i i - -г V-г2 - 1
1
где под л/-г2 — 1 подразумевается фиксированная ветвь, однозначная в плоскости с разрезом вдоль отрезка [—1,1] действительной оси и принимающая на (1, то) действительные значения.
Литература
1. Крылов В. И. Приближенные вычисления интегралов.—М.: Наука, 1967.—500 с.
2. Крылов В. И., Шульгин Л. Т. Справочная книга по численному интегрированию.—М.: Наука, 1966.—370 с.
3. Иванов В. В. Теория приближенных методов и ее применение к численному решению сингулярных интегральных уравнений.—Киев: Наукова думка, 1968.—288 с.
4. Лифанов И. К. Метод сингулярных интегральных уравнений и численный эксперимент.—М.: ТОО Янус, 1995.—520 с.
5. Пашковский С. Вычислительные применения многочленов и рядов Чебышева.—М.: Наука, 1983.— 384 с.
6. Панасюк В. В., Саврук М. П., Дацыщин А. П. Распределение напряжений около трещин в пластинах и оболочках.—Киев: Наукова думка, 1976.—444 с.
7. Ахиезер Н. И. О некоторых формулах обращения сингулярных интегралов // Изв. АН СССР. Математика.—1945.—Т. 9.—С. 275-290.
8. Сеге Г. Ортогональные многочлены.—М.: Физматгиз, 1962.—500 с.
9. Пыхтеев Г. Н. Точные методы вычисления интегралов типа Коши.—Новосибирск: Наука, 1980.— 120 с.
10. Плиева Л. Ю., Бесаева З. В. Об одной квадратурной формуле и некоторе ее применение // Тр. междунар. симпозиума (МДОЗМФ-2009).—Харьков-Херсон, 2009.—С. 141-144.
Статья поступила 14 февраля 2011 г. Хубежты Шалва Соломонович
Южный математический институт ВНЦ РАН и РСО-А, заведующий лаб. математического моделирования РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22 E-mail: [email protected]
Плиева Любовь Юрьевна
Южный математический институт ВНЦ РАН и РСО-А, младший научный сотрудник РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22 E-mail: [email protected]
Бесаева Зарина Вячеславовна
Юго-Осетинский государственный университет им. А. А. Тибилова, преподаватель
ЮЖНАЯ ОСЕТИЯ, 383570, Цхинвал, Московская, 8 E-mail: [email protected]
QUADRATURE FORMULAS FOR SINGULAR INTEGRALS WITH WEIGHT FUNCTIONS
Khubezhty Sh. S., Plieva L. Yu., Besaeva Z. V.
We study singular integrals with weight functions and construct quadrature formulas for such integrals. Inversion formulas similar to the formulas for singular integrals with Chebyshev polynomials are proved. Application of obtained quadrature formulas to the numerical solution of singular integral equations are pointed out.
Key words: singular integral, inversion formulas, quadrature formulas.
