CC BY
60
№ 4 (32), 2014
Физико-математические науки. Математика
УДК 517.392
И. В. Бойков, Ю. Ф. Захарова, М. А. Семов
ПРИБЛИЖЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ ГИПЕРСИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ
Аннотация.
Актуальность и цель. Приближенные методы решения гиперсингулярных интегральных уравнений являются активно развивающимся направлением вычислительной математики, что в первую очередь связано с многочисленными приложениями гиперсингулярных интегральных уравнений к механике, аэродинамике, электродинамике, геофизике. При этом следует отметить два обстоятельства: 1) аналитическое решение гиперсингулярных интегральных уравнений возможно лишь в исключительных случаях; 2) спектр приложений гиперсингулярных интегральных уравнений постоянно расширяется. Этим обусловлена актуальность построения и обоснования численных методов решения гиперсингулярных интегральных уравнений. В настоящее время остались не разработанными методы приближенного решения нелинейных гиперсингулярных интегральных уравнений. Статья посвящена построению и обоснованию приближенного решения одного класса нелинейных гиперсингулярных интегральных уравнений методом коллокаций.
Материалы и методы. Обоснование разрешимости и сходимости метода коллокаций к приближенному решению одного класса нелинейных гиперсингулярных интегральных уравнений, определенных на замкнутых контурах, основано на применении методов функционального анализа и теории приближений.
Результаты. Предложен и обоснован метод коллокаций для приближенного решения нелинейных гиперсингулярных интегральных уравнений, определенных на замкнутых контурах. Приведены оценки быстроты сходимости и величины погрешности.
Выводы. Построена вычислительная схема, позволяющая эффективно решать прикладные задачи механики, аэродинамики, электродинамики, геофизики.
Ключевые слова: нелинейные гиперсингулярные интегральные уравнения, метод коллокаций, метод Ньютона - Канторовича.
I. V. Boykov, Yu. F. Zakharova, M. A. Semov
APPROXIMATE SOLUTION OF NONLINEAR HYPERSINGULAR INTEGRAL EQUATIONS
Abstract.
Background. Approximate methods for solving hypersingular integral equations are an actively developing section of calculus mathematics owing to numerous applications of hypersingular integral equations in mechanics, aerodynamics, electrodynamics, geophysics. Here it is necessary to mention two important circumstances:
1) analytical solution of hypersingular integral equations is possible only in exceptional cases; 2) application range of hypersingular integral equations is constantly growing. These prove the relevance of building and substantiating numerical methods for solving hypersingular integral equations. At the present time the methods of approximate solution of nonlinear hypersingular integral equations remain undeveloped. The article is devoted to building and substantiating an approximate solution of one class of nonlinear hypersingular integral equations by the method of collocations.
Physical and mathematical sciences. Mathematics
69
Известия высших учебных заведений. Поволжский регион
Materials and methods. Substantiation of solvability and convergence of the method of collocations to the approximate solution of one class of nonlinear hypersingular integral equations, determined on closed loops, is based on application of methods of functional analysis and approximation theory.
Results. The authors suggested and substantiated the method of collocations for the approximate solution of nonlinear hypersingular integral equations, determined on closed loops. The article includes the estimates of convergence rate and extent of error.
Conclusions. The authors built a computing circuit allowing to effectively calculate applied problems of mechanics, aerodynamics, electrodynamics, geophysics.
Key words: nonlinear hypersingular integral equations, method of collocations, method of Newton-Kantorovich.
Введение
Численные методы решения гиперсингулярных интегральных уравнений принадлежат к одному из наиболее активно развивающихся направлений вычислительной математики. В первую очередь это обусловлено многочисленными приложениями гиперсингулярных интегральных уравнений к задачам аэродинамики, электродинамики, ядерной физики. Решению линейных гиперсингулярных интегральных уравнений посвящено большое число работ, в которых предложены и обоснованы различные численные методы. Среди этих работ отметим монографию [1], в которой метод дискретных особенностей применяется для приближенного решения задачи бесциркулярного обтекания поверхности потоком газа; цикл работ [2-10], в которых предложены и при ряде дополнительных условий обоснованы проекционные методы решения гиперсингулярных интегральных уравнений первого и второго родов.
В то же самое время приближенные методы решения нелинейных гиперсингулярных интегральных уравнений с нелинейностями в гиперсингулярных интегралах до сих пор остаются неисследованными. Насколько авторам известно, исключением является работа [11], в которой предложен и обоснован сплайн-коллокационный метод нулевого порядка решения гиперсингулярных интегральных уравнений вида
a(t)х(t) + f h(t’TX(2T)) dт = f{t), te (-1,1).
-1 (т-t )2
В данной работе предложен и обоснован метод коллокации для решения гиперсингулярных интегральных уравнений вида
a(t)x(t) + f hU' ТX(T))dт = f (t), te Y,p = 2,3,...,
Y ( т -t)p
где y - единичная окружность с центром в начале координат в плоскости комплексной переменной z .
1. Приближенное решение линейных гиперсингулярных интегральных уравнений методом коллокации
Для обоснования приближенного решения нелинейных гиперсингулярных интегральных уравнений методом коллокации ниже используется метод Ньютона - Канторовича, основанный на линеаризации нелинейного оператора.
70
University proceedings. Volga region
№ 4 (32), 2014
Физико-математические науки. Математика
При этом будут использованы результаты, полученные в работе [12] для приближенного решения линейных гиперсингулярных интегральных уравнений методом коллокации. Приведем ниже эти результаты.
Предварительно приведем определения, используемые в работе.
Определение 1. Класс функций Гельдера Ha(M;[a,b]) (0 < a<1)
состоит из заданных на отрезке [a,b] функций f (x), удовлетворяющих во всех точках х' и х" этого отрезка неравенству
|f (х') - f (х'')|< M | х - х''\а .
В случае, когда из текста ясно, на каком множестве рассматриваются функции, вместо Ha(M;[a,b]) будем писать Ha(M). Это замечание относится и к остальным классам функций.
Определение 2. Класс Wr(M;[a,b]) состоит из функций, заданных на отрезке [a, b], непрерывных и имеющих непрерывные производные до (г - 1) -го порядка включительно и кусочно-непрерывную производную r -го
порядка, удовлетворяющую на этом отрезке неравенству | f(r)(х) |< M.
Определение 3. Класс WrHa(M;[a,b]) состоит из функций f (х), принадлежащих классу Wr(M;[a,b]) и удовлетворяющих дополнительному условию f(r)(х)е Ha(M).
Напомним определение гиперсингулярных интегралов, обобщающее понятия интеграла в смысле главного значения Коши и интеграла в смысле Адамара.
Определение 4 [13]. Интегралом
г ф(х)^т
a (т-c) p,
a < c < b,
в смысле главного значения Коши - Адамара называется следующий предел:
—V ф(т)Лт + b 9(T)dт +_|(v)
J (т-c)p J (т-c)p vp-1 ,
где ^(v) - некоторая функция, выбранная так, чтобы указанный предел существовал.
Через [X ,Y ] обозначено множество линейных ограниченных операторов, отображающих нормированное пространство X в нормированное пространство Y .
Рассмотрим гиперсингулярное интегральное уравнение
Кх = a^) + Гh(t’Т)х(pТ = f(t), tеу, (1)
Y (т-t)p
ф(т)^т _
- lim
(т-c У
Physical and mathematical sciences. Mathematics
71
Известия высших учебных заведений. Поволжский регион
где у - единичная окружность с центром в начале координат в плоскости комплексной переменной, p = 2,3,...; a(t), f (t)e Ha, 0<a<1; функция h(t,т) принадлежит классу Гельдера Ha по первой переменной и имеет производные до (p — 1) -го порядка по второй переменной, причем
д p 1h(t, т) дт p—1
e H а
0< a< 1.
Введем пространство X функций x(t), t e у, x(t) eWp 1He, в < а,
с нормой
p—1
l|x(t )|| = Z
k=0
x(k)(t)
+ sup
C( Y) ?1 ^?2; ^, r2eT
x(p—1>(t|) - x(p—1>(t2)|
111 —12 |в
Через Y обозначим пространство функций y(t), t e у, y(t) e Hp, в < а, с нормой
+
f1
sup
Ф ?2; ^1, t2eJ
I y(t1)— y(t2)l 111 —12 |P
Известно, что пространства X и Y банаховы.
Из определения гиперсингулярных интегралов и теории сингулярных операторов следует, что оператор K отображает пространство X в Y.
Будем считать, что K имеет линейный непрерывный обратный оператор K—1 e [Y, X ].
Обозначим через Xn подпространство пространства X, состоящее из функций
p—2 n+p—1 —1
xn(t )= Z akt>k ln t + Z akt>k + Z akt>k, (2)
k=0 k=p—1 k=—n+p—1
а через Yn - подпространство пространства Y, состоящее из полиномов
Уп (t)= Z о/.
k=—n
Введем узлы tk = eSk, Sk =2kn / (2n +1), k = 0,1, ...,2n , и через Pn [f] обозначим оператор проектирования непрерывных функций f e C(у) на множество интерполяционных полиномов n -го порядка, построенных на
узлах tk, k = 0,1,..,2n.
Приближенное решение уравнения (1) будем искать в виде функции
p—2 n+p—1 —1
xn(t)= Z aktk !nt + Z aktk + Z aktk,
k=0 k=p—1 k=—n+p—1
72
University proceedings. Volga region
№ 4 (32), 2014
Физико-математические науки. Математика
коэффициенты (а^} которой определяются по методу коллокации из системы линейных алгебраических уравнений, представимых в операторном виде уравнением
K х = P
L^n^n 1 n
a(t) Xn (t) + J
h(t, t) xn (x)d т (T-t)p
= Pn [ f (t)].
(3)
Справедливо следующее утверждение.
Теорема 1 [12]. Пусть оператор K е [X,Y] непрерывно обратим:
ai, f е На,h(t,т)е Наа. Тогда при n таких, что q = Cln2n /nа-в < 1, система
*
уравнений (3) имеет едиственное решение xn и справедлива оценка
II * * и,- Cln2n *
IIх -xn 11^-а-p-, где х
решение уравнения (1).
2. Метод коллокаций для нелинейных гиперсингулярных интегральных уравнений
Рассмотрим уравнение
L(х) = a(t)x(t) + Jh(t’T’X(T)}dT = f(t), tey, p = 2,3,... (4)
Y (т-t)p
Будем считать, что a (t), f (t) e На ,0< a< 1; функция h(t, T, u) удовлетворяет условию Гельдера с показателем а по первой переменной и имеет частные производные до p -го порядка включительно по двум последним переменным, удовлетворяющие условию Гельдера с показателем а .
Приближенное решение уравнения (4) будем искать в виде функции (2), коэффициенты которой определяются из системы нелинейных алгебраических уравнений
Ln ( xn ) Pn
a (t) xn (t) + J
h(t, t, x„ (T))d т (T-t)p
= p„ [ f (t)].
(5)
Исследовать приближенное решение уравнения (4) будем в пространствах X, Y и их подпространствах Xn, Yn, введенных в предыдущем разделе. Нетрудно видеть, что производные Фреше операторов L(x) и Ln (xn)
на элементе x0 соответственно равны
L'(x0)x = a(t)x(t) + fh'3(t-T;x (T))x(T) dT, JY (T-t)p
L'n (x ) xn = Pn
a t x„(t) + f h'3(t,T,x (T»x(T) dт JY (T-t)p
Physical and mathematical sciences. Mathematics
73
Известия высших учебных заведений. Поволжский регион
где т,и) означает производную по третьему аргументу функции
h(t, т, и).
Будем считать, что оператор L'(x0) имеет линейный непрерывный
обратный оператор [L'(x0)]-1 е [Y,X].
В разд. 1 было показано, что при выполнении условия
q = Cln2n / na—e < 1
оператор L'n(x0)е [Xn,Yn] имеет линейный непрерывный обратный оператор [Ln'(xo)]-1 е [Yn,Xn].
Для доказательства разрешимости системы уравнений (5) воспользуемся методом Ньютона - Канторовича.
Пусть X, Y - банаховы пространства.
Рассмотрим уравнение
Kx = 0, (6)
где K — нелинейный оператор, действующий из X в Y.
Теорема 2 [14]. Пусть X и Y - банаховы пространтства, и пусть выполнены условия:
1) || к (xo)||= no;
2) оператор K имеет произвольную Фреше в окрестности точки x0, и существует правый обратный оператор K( x0)]r1 с нормой
KX x0)]—1
= B0;
3) в сфере S :{x :|| x — x0 ||< В0П0 / (1 — q)}, q <1, выполняется условие || K(x1) — K(x2)||< q/(B0(1 + q)).
Тогда уравнение (6) имеет в S решение x , к которому сходятся приближения xn+1 = xn — [K'(xn )]r1 Kxn, n = 0,1,..., и справедлива оценка
|| x*— xn ||< qnП0B0n0/(1—q).
Теорема 3 [14]. Пусть X, Y - банаховы пространства, и пусть выполнены условия:
1) У Kx0 ||= П0;
2) оператор K имеет производную Фреше в окрестности точки x0, и
существует правый обратный оператор [K( x0)]—1 с нормой ||[ K( x0)]—1 ||= B0;
3) в сфере S :{x :|| x — x0 ||< В0П0 / (1 — q)}, q <1, выполняется условие || K'(x1) — K'(x2)||< q/%
Тогда уравнение (6) имеет в S решение x , к которому сходятся приближения xn+1 = xn — [ K( x0)]r1 Kxn, и справедлива оценка
74
University proceedings. Volga region
№ 4 (32), 2014
Физико-математические науки. Математика
II х*- х„ У< ч„ ПоАо/(1- ч).
Будем искать приближенное решение уравнения (5) модифицированным методом Ньютона - Канторовича:
хП+Х= хП - [L'„(xo)]-1(Ln (xm) - fn), m = 0,1,..., (7)
где fn = P„ [f ]-
Проверим выполнение условий теоремы 3.
Справедливость условия 2 следует при достаточно больших n из теоремы 1.
Для доказательства сходимости итерационного процесса (7) достаточно проверить выполнение условия 3 теоремы 3.
Оценим норму разности
11 L„ (х,)xn -L„ (х")xn ||с(Y)<
<|
Pn \\с(Y)
Г h3 (t, т, х'(т)) - h3 (t, т, х"(т))) xn (t) d t J (т-t) p
C (Y)
Известно, что || Pn ||с(y)< Clnn. Воспользовавшись определением гиперсингулярных интегралов и известными свойствами сингулярных интегралов, имеем
J
h3 (t, т, х'(т)) - h3 (t, т, х"(т)) хп (т) (т-t)p
d т
<
< C
J
p) (t, т, х'( t)) - h<3p) (t, т, х"(т)) xn (t)
т-t
d т
+... +
+c
h-3 (t, т, х'(т)) - h3 (t, т, х"(т)) х„р 1)(t) d т
<
< сЦх'(т) - х"(т)||X \\х„ (т)||X .
Следовательно,
|Il,„ (х) х„ - L„ ( х„ ) х„ 1C ( y ) <C ln „Iх (t) - х (т)11 x 11х„ (t)II x .
Так как (L„ (х") х„ - L„ (х„)х„) является полиномом п -го порядка, то
WL„(х)х„ - L„(х„)х„ Wy < C„e ln„ W х'(т) - х"(т) Ух У х„(т) Wx и ||L„(х") -L„(х„)|Y < Cne lnп ||х'-х" ||х .
Physical and mathematical sciences. Mathematics
75
Известия высших учебных заведений. Поволжский регион
Таким образом, если Xq достаточно хорошее начальное приближение
*
к точному решению x . то выполнены условия теоремы 3 и итерационный
*
процесс (6) сходится к решению x . Справедливо следующее утверждение.
*
Теорема 4. Пусть уравнение (4) имеет решение x (t) е Ha. Пусть a(t), f (t) е Ha, 0<a<1; функция h(t,t,u) удовлетворяет условию Гельдера
с показателем a по первой переменной и имеет частные производные до p -го порядка включительно по двум последним переменным, удовлетворяющие условию Гельдера с показателем a .
*
Пусть в шаре S(x , r) выполнены следующие условия:
1)
L( x 0)
по;
2) существует обратный оператор [L(x0)] 1 с нормой [L(x0)]'
0чП-1
= Во
(здесь х0 - начальное приближение в итерационном процессе (6)). Тогда
п, ч „in2 n ,
существует шар В(xq,r) такой, что при n таких, что q = C—-в< 1, итера-
na-p
*
ционный процесс (7) сходится к решению хп уравнения (5) и справедлива оценка
х - х„
< C
ln2n
na-e.
Список литературы
1. Вайникко, Г. М. Численные методы в гиперсингулярных интегральных уравнениях и их приложения / Г. М. Вайникко, И. К. Лифанов, Л. Н. Полтавский. -М. : Янус-К, 2001. - 508 с.
2. Бойков, И. В. Приближенное решение гиперсингулярных интегральных уравнений / И. В. Бойков, Е. Г. Романова // International Conference on Computational Mathematics. - Part first. - Novosibirsk, 2004. - P. 411-417.
3. Бойков, И. В. Коллокационный метод решения гиперсингулярных интегральных уравнений / И. В. Бойков, Е. Г. Романова // Известия высших учебных заведений. Поволжский регион. Естественные науки. - 2006. - № 5. - С. 42-50.
4. Бойков, И. В. Сплайн-коллокационный метод решения гиперсингулярных интегральных уравнений / И. В. Бойков, А. И. Бойкова // Вестник Харьковского национального университета. Сер. Математическое моделирование. Информационные технологии. Автоматизированные системы управления. - 2007. - С. 36-49.
5. Бойков, И. В. Приближенное решение гиперсингулярных интегро-дифференциальных уравнений / И. В. Бойков, Ю. Ф. Захарова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2010. -№ 1 (13). - С. 80-90.
6. Бойков, И. В. Приближенные методы решения сингулярных и гиперсингулярных интегродифференциальных уравнений / И. В. Бойков, Ю. Ф. Захарова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2012. - № 3 (23). - С. 99-113.
7. Boykov, I. V. An Approximate Solution of Hypersingular Integral Equations / I. V. Boykov, E. S. Ventsel, A. I. Boykova // Appl. Num. Math. - 2010. - № 60. -P. 607-628.
76
University proceedings. Volga region
№ 4 (32), 2014
Физико-математические науки. Математика
8. Ganesh, M. The numerical solution of a nonlinear hypersingular boundary integral equations / M. Ganesh and O. Steinbach // J. Comput. Appl. Math. - 2001. - № 131. -P. 267-280.
9. Оселедец, И. В. Приближенное обращение матриц / И. В. Оселедец, Е. Е. Тыртышников // Журнал вычислительной математики и математической физики. - 2005. - Т. 45, № 2. - C. 315-326.
10. Capobiano, M. R. Newton methods for a class of nonlinear hypersingular integral equations / M. R. Capobiano, G. Criscuolo, and P. Junghanns // Numer. Algorithms. -2010. - № 55. - P. 205-221.
11. Boykov, I. V. An approximate solution of nonlinear hypersingular integral equations / I. V. Boykov, E. S. Ventsel, V. A. Roudnev, A. I. Boykov // Applied Numerical Mathematics. - 2014. - Vol. 86, December. - P. 1-21.
12. Бойков, И. В. Приближенное решение линейных гиперсингулярных интегральных уравнений методом коллокаций / И. В. Бойков, Ю. Ф. Захарова, М. А. Семов, А. А. Есафьев // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. - № 3 (31). - С. 101-113.
13. Чикин Л . А . Особые случаи краевой задачи Римана и сингулярных интегральных уравнений / Л. А. Чикин // Ученые записки Казанского государственного университета. - 1953. - Т. 113, № 10. - С. 57-105.
14. Бойков, И. В. Приближенное решение сингулярных интегральных уравнений / И. В. Бойков. - Пенза : Изд-во ПГУ, 2004. - 316 с.
References
1. Vaynikko G. M., Lifanov I. K., Poltavskiy L. N. Chislennye metody v gipersingul-yarnykh integral’nykh uravneniyakh i ikh prilozheniya [Numerical methods in hypersingular integral equations and application thereof]. Moscow: Yanus-K, 2001, 508 p.
2. Boykov I. V., Romanova E. G. International Conference on Computational Mathematics. Part first. Novosibirsk, 2004, pp. 411-417.
3. Boykov I. V., Romanova E. G. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Estestvennye nauki [University proceedings. Volga region. Natural sciences]. 2006, no. 5, pp. 42-50.
4. Boykov I. V., Boykova A. I. Vestnik Khar’kovskogo natsional’nogo universiteta. Ser. Matematicheskoe modelirovanie. Informatsionnye tekhnologii. Avtomatizirovannye sis-temy upravleniya [Bulletin of Kharkov National University. Series: Mathematical modeling. Information technologies. Automatic control systems]. 2007, pp. 36-49.
5. Boykov I. V., Zakharova Yu. F. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2010, no. 1 (13), pp. 80-90.
6. Boykov I. V., Zakharova Yu. F. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2012, no. 3 (23), pp. 99-114.
7. Boykov I. V., Ventsel E. S., Boykova A. I. Appl. Num. Math. 2010, no. 60, pp. 607628.
8. Ganesh M. and Steinbach O. J. Comput. Appl. Math. 2001, no. 131, pp. 267-280.
9. Oseledets I. V., Tyrtyshnikov E. E. Zhurnal vychislitel’noy matematiki i matematich-eskoy fiziki [The journal of calculus mathematics and mathematical physics]. 2005, vol. 45, no. 2, pp. 315-326.
10. Capobiano M. R., Criscuolo G. and Junghanns P. Numer. Algorithms. 2010, no. 55, pp. 205-221.
11. Boykov I. V., Ventsel E. S., Roudnev V. A., Boykov A. I. Applied Numerical Mathematics. 2014, vol. 86, December, pp. 1-21.
Physical and mathematical sciences. Mathematics
77
Известия высших учебных заведений. Поволжский регион
12. 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.
13. Chikin L. A. Uchenye zapiski Kazanskogo gosudarstvennogo universiteta [Memoirs of Kazan State University]. 1953, vol. 113, no. 10, pp. 57-105.
14. Boykov I. V. Priblizhennoe reshenie singulyarnykh integral’nykh uravneniy [Approximate solution of singular integral equations]. Penza: Izd-vo PGU, 2004, 316 p.
Бойков Илья Владимирович
доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: [email protected]
Захарова Юлия Фридриховна
кандидат физико-математических наук, доцент, кафедра высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: [email protected]
Семов Михаил Александрович аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: [email protected]
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)
Zakharova Yuliya Fridrikhovna Candidate of physical and mathematical sciences, associate professor, sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)
Semov Mikhail Aleksandrovich Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)
УДК 517.392 Бойков, И. В.
Приближенное решение нелинейных гиперсингулярных интегральных уравнений / И. В. Бойков, Ю. Ф. Захарова, М. А. Семов // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. - № 4 (32). - С. 69-78.
78
University proceedings. Volga region
