https://doi.org/10.21122/2227-1031-2021-20-2-168-172 УДК 519.6517.5
Приближенное представление дилогарифмами решения одной вариационной краевой задачи для круга при граничном условии Неймана
Докт. физ.-мат. наук, проф. И. Н. Мелешко1), канд. физ.-мат. наук, доц. П. Г. Ласый1)
^Белорусский национальный технический университет (Минск, Республика Беларусь)
©Белорусский национальный технический университет, 2021 Belarusian National Technical University, 2021
Реферат. Известно, что краевые задачи для уравнений Лапласа и Пуассона эквивалентны задаче вариационного исчисления - о минимуме интеграла, для которого данное уравнение в частных производных является уравнением Эйлера - Лагранжа. Например, задача о минимуме интеграла Дирихле в единичном круге с центром в начале координат на некотором допустимом множестве функций при заданных значениях нормальной производной на окружности эквивалентна краевой задаче Неймана для уравнения Лапласа в этой области. На основе известного точного решения краевой задачи Неймана для круга с помощью специальной приближенной формулы для интеграла Дини сконструировано эффективное приближенное представление дилогарифмами решения указанной выше эквивалентной вариационной краевой задачи. Приближенная формула эффективна в том смысле, что она достаточно проста при численной реализации, устойчива, а равномерная по кругу оценка погрешности позволяет проводить вычисления с заданной точностью. Специальная квадратурная формула для интеграла Дини обладает замечательным свойством - ее коэффициенты неотрицательны. Квадратурные формулы с неотрицательными коэффициентами занимают особое место в теории приближенных вычислений определенных интегралов и ее приложениях. Естественно, что еще большую значимость это свойство приобретает, когда коэффициенты не числа, а некоторые функции. Проведенный численный анализ приближенного решения подтверждает его эффективность.
Ключевые слова: вариационная задача, краевая задача Неймана, интеграл Дини, приближенное решение, квадратурная формула
Для цитирования: Мелешко, И. Н. Приближенное представление дилогарифмами решения одной вариационной краевой задачи для круга при граничном условии Неймана / И. Н. Мелешко, П. Г. Ласый // Наука и техника. 2021. Т. 20, № 2. С. 168-172. https://doi.org/10.21122/2227-1031-2021-20-2-168-172
Approximate Dilogarithm Representation
of One Variational Boundary Value Problem Solution for Circle under the Neumann Boundary Condition
I. N. Meleshko1), P. G. Lasy1)
^Belarusian National Technical University (Minsk, Republic of Belarus)
Abstract. It is known that boundary value problems for the Laplace and Poisson equations are equivalent to the problem of the calculus of variations - the minimum of an integral for which the given partial differential equation is the Euler - Lagrange equation. For example, the problem of the minimum of the Dirichlet integral in the unit disc centered at the origin on some admissible set of functions for given values of the normal derivative on the circle is equivalent to the Neimann boundary value problem for the Laplace equation in this domain. An effective approximate dilogarithm representation of the solution of the above equivalent variational boundary value problem is constructed on the basis of the known exact solution of the Neumann Boundary value problem for a circle using a special approximate formula for the Dini integral. The approximate formula is effective in the sense that it is quite simple in numerical implementation, stable, and the error estimation, which is uniform over a circle, allows calculations with the given accuracy. A special quadrature formula for the Dini integral has
Адрес для переписки
Мелешко Иван Николаевич
Белорусский национальный технический университет
ул. Б. Хмельницкого, 9,
220013, г. Минск, Республика Беларусь
Тел.: +375 17 292-82-73
Address for correspondence
Meleshko Ivan N.
Bebrusian National Technical University 9, Б. Khmelnithskogo str., 220013, Minsk, Republic of Belarus Tel.: +375 17 292-82-73 [email protected]
Наука
итехника. Т. 20, № 2 (2021)
a remarkable property - its coefficients are non-negative. Quadrature formulas with non-negative coefficients occupy a special place in the theory of approximate calculations of definite integrals and its applications. Naturally, this property becomes even more significant when the coefficient are not number, but some functions. The performed numerical analysis of the approximate solution confirms its effectiveness.
Keywords: variational problem, boundary Neumann problem, Dini integral, approximate solution, quadrature formula
For citation: Meleshko I. N., Lasy P. G. (2021) Аpproximate Dilogarithm Representation of One Variational Boundary Value Problem Solution for Circle under the Neumann Boundary Condition. Science and Technique. 20 (2). 168-172. https://doi.org/10.21122/2227-1031-2021-20-2-168-172 (in Russian)
Введение
Вариационные методы часто применяют для решения различных задач механики, физики и математической физики [1-4]. Представляет интерес вопрос о решении вариационных задач с помощью эффективных методов решения эквивалентных краевых задач математической физики.
Рассмотрим задачу о минимуме функционала на некотором допустимом множестве 3 функций и = и(х, у), удовлетворяющих граничному условию Неймана:
Ц| gradu |2 dxdy.
Она эквивалентна краевой задаче Неймана для уравнения Лапласа в области G. В случае, когда область G - единичный круг с центром в начале координат, мы приходим к задаче о нахождении функции u = и(г, ф), удовлетворяющей уравнению Лапласа и граничному условию Неймана соответственно:
8ы 8r
Ды = 0, r <1;
= f(ф), -Л<ф<Л,
(1)
(2)
где /(ф) - непрерывная на отрезке [-л, л] функция.
Необходимым условием разрешимости краевой задачи (1), (2) является выполнение равенства
л
J f (ф^ф = 0.
Решение задачи находится в этом случае с точностью до действительной постоянной и может быть записано с помощью интеграла Дини [1, с. 598-600; 2, с. 231-232] в виде
1 л
i(r, ф) = - - J f (x)ln| t - z|dx + С,
IT
(3)
Решение краевой задачи (1), (2), а, следовательно, и соответствующей вариационной задачи будет единственным, если потребовать, чтобы множество допустимых функций 3 удовлетворяло некоторому дополнительному условию.
ОСНОВНАЯ ЧАСТЬ Приближенная формула для интеграла Дини
Сконструируем квадратурную формулу с неотрицательными коэффициентами для интеграла
- л
D(r, ф) = - - J f (x)ln| t - z|dx,
7Г J
t = e", z = гегф.
(4)
Запишем этот интеграл в виде разности двух интегралов
0(г, ф)=-Л{ /(т)1п-Цат- — } /(т)ёт. (5)
л J 11 - г | л
-л 1 1 -л
Зададим на отрезке [-л, л] систему точек
фк = кк, к = —п, п, к = 2л и аппроксимируем 2п + 1
функцию /(ф) на этом отрезке по формуле
/ (ф) - / (ф)= Е бк (ф) / (фк), (6)
k=-n
где 0k (ф) =
1, фе
0, фй
Фk
ф^
Фk -Фk +-
Подставив в представление (5) интеграла Дини вместо плотности ее приближение по (6), получим следующую квадратурную формулу:
D(r, ф) * D(r, ф) = £ Ak (r, ф) f (фk) -
k=- n
h ln2
(7)
t = en, z = reK*.
£ f fok),
k=-
Наука
итехника. Т. 20, № 2 (2021)
G
r=1
где
h
фк +2
Л(r, ф) = 1 f ln—^dx.
л J 11 - z |
(8)
Фк—
к 2
Теорема 1. Коэффициенты Ак(г, ф) квадратурной формулы (7) неотрицательны и вычисляются по формуле
a (r, ф) =
(
h ln 2 + Im
f i
L2
ze
-ц Фк-2
h ЛЛ
V v
- L2
ze
i Фк +2
(9)
^JJ
где L2(z) = E - дилогарифм Эйлера [5, 6].
к=1 k
Доказательство. Так как в единичном кру-
ге ln-
2
• > 0, то из представления (8) для ко-
^ - z|
эффициентов Ак(г, ф) следует, что все они неотрицательны для г < 1 и ф е [-л, л]. Далее запишем
Ak (r, Ф) =
h ln2 1
h
фк +2
л л
- f Reln h - ^Z 1 dx
h
Фк— к2
|D(r, ф) - D(r, ф)|< 4ю(f, h) ln 2, (10)
где ra(f, h) - модуль непрерывности функции т(ф).
Если же Дф) - непрерывно дифференцируемая на [-л, л] функция, то равномерно по всем r < 1 и ф е [-л, л]
| D(r, ф) - D (r, ф)| < 2M1h ln2, (11) где Mj = max | f '(ф) |.
фе[-л, л ]
Доказательство. Сравнивая (5) и (7), находим, что
I (f (X) - f (T))ln-2
л
D(r, ф) - D(r, ф) = 11 (f (x) - f (x))ln-
7Г J
-dx-
t - z |
ln2 л
л
I (f (X) - f (X))dx.
Оценивая эту разность по абсолютной величине и учитывая, что в единичном круге
ln-
2
|t - z |
• > 0, можно записать неравенство |D(r, ф) - D(r, ф)|<
1 л
1 f ln"
ТГ J I
л
V -л
2 , ln2 -dx +-
t - z | л
л
I dX
max | f (ф) - f (ф)|,
фе[ -л, л ]
и, значит,
фк +2 ( , чкЛ
h ln2 1 г2 ^ 1 f ^
л л
л I Re[Sft
dx =
фк -
h ln2 1
фк+
л л
1 f IEкcosк(х-ф)
л h V к=1 к
dx =
h ln2 1
- — X
л л
фк -1
f ад ^ f
=1 к2
f
sinкI ф-фк +-I-sinкI ф-фк--
1
h ln2 + Im
f f
L2
ze
-Чфк- 2
h f - L2
V
ze
i фк +2
I
h
Формула (9) доказана.
Получим оценки погрешности приближенной формулы (7).
Теорема 2. Если плотность /(ф) интеграла Дини непрерывна на отрезке [-л, л], то имеет место равномерная по г < 1 и ф е [-л, л] оценка погрешности приближенной формулы (7)
| D(r, ф) - D(r, ф)|<
<
I f n_±
V л -^i | t - z |
Поскольку
Л
-dx + 4ln2
max |f (ф)- f (ф)|.
фе[-л, л ]
(12)
ln-^ = - Reln [ 1 - - | =
11 - z | V t J
f ад 1 f z 1 ^ —k
= Re I = E—cosк(х-ф),
к=1к Vt J к=1 к
то
л 1 / I 1п-dx = V —-8т к(х-ф)Г =0
-л К - Z | {-I к2 У
и, следовательно, неравенство (12) можно переписать в виде
|5(г, ф) - 5(г, ф)|< 41п2 тах |/(ф) - /(ф)|. (13)
фе[-л ,л ]
Если функция /(ф) непрерывна на отрезке [-л, л], то
Наука
итехника. Т. 20, № 2 (2021)
h
|/(ф) - /(ф)|<ш(/, к), фе [-л, л], (14)
если же она непрерывно дифференцируема на этом промежутке, то
М,к
|f (ф) - f (ф)|<-
2
фе [-л, л]. (15)
Из (13)-(15) следуют доказываемые неравенства (10), (11).
Полученные оценки погрешности приближенной формулы (7) для интеграла Дини полезны также и тем, что в них указываются не только порядок, но и возможные значения для констант.
Приближенное представление дилогарифмами решения вариационной задачи
Приближенное решение вариационной задачи будем конструировать на основе формул (3)-(5) и (7), из которых вытекают следующие точные и приближенные равенства:
и(г, ф) = В(г, ф) + С - В(г, ф) + С; (16)
п
и (г, ф) = б(г, ф) + С = £ А (г, ф)/(фк) -
k=-n
hln2
(17)
£ f (ф:)+с,
k=-
в которых коэффициенты Ак(г, ф) определены формулами (8), (9). Формулы (9) и (17) позволяют представить приближенное решение дилогарифмами Эйлера
1
i(r, ф) = - £ Im
( f
L2
' k=—n
ze
i фk-2
h Л Л
V v
- L2
ze
Фk+2
h ЛЛЛ
(18)
f (Фk) + c.
//
Сравнив равенства (16) и (17), замечаем, что
и (г, ф) - и (г, ф) = Б(г, ф) - 3(г, ф).
Поэтому для оценки погрешности формул (17) и (18) можно воспользоваться неравенствами (10) и (11) теоремы 2.
Об устойчивости квадратурной формулы (7) для интеграла Дини
Выше отмечен тот факт, что коэффициенты Ак(г, ф) квадратурной формулы (7) неотрицательны для всех г < 1 и ф е [-л, л]. Квадратур-
Наука
итехника. Т. 20, № 2 (2021)
ные формулы с неотрицательными коэффициентами занимают особое место в теории приближенного вычисления определенных интегралов и ее приложениях [7]. Важность свойства положительности коэффициентов квад-квадратурных формул отмечается также в [8]. Естественно, что еще большую значимость это свойство приобретает, когда коэффициенты не числа, а некоторые функции, на что указывается в [9].
Проведем исследование квадратурной формулы (7) на устойчивость. Предположим, что значения плотности /(ф) интеграла Дини в узлах квадратурной формулы найдены приближенно,
т. е. вместо /(ф) имеем /(ф), так что погрешности равны вк = /(ф) - /(ф). Пусть | вк | <в,
к = -п, п. Тогда для погрешности вычисления квадратурной суммы в (7) получается неравенство
£ Ак (г, ф)/(фк) - £ Ак (г, ф)/(фк)
¡г=-п к=-п
<в£| Ак (г, ф)|,
к=-п
из которого следует, что точная верхняя грань погрешности вычисления этой квадратурной
п
суммы пропорциональна £ |Ак (г, ф)|. Оценим
к=-п
данную сумму числом, не зависящим от г и ф, учитывая, что, как мы убедились выше,
ln-
2
11 - z |
>0 и
- 1
I ln-dx = 0 и, следовательно,
J 11 - z |
П 1 n
£ |Ak (r, Ф)|=Л£
k=-n
' k=-n
h
%+-
J ln
Ф1:--
k 2
л 1
2л ln2 + i ln-dx
J 11 - z |
t - z | Л
-dx
= 2ln2.
/
Значит, при всех г < 1 и ф е [-л, л] справедливо неравенство
£ Ak (r, ф) f (Ф: ) - £ Ak (r, ф) f (Ф: )
k=-n
k=-n
< 2sln2.
Это означает, что при всех п, г и ф погрешность вычислений этой квадратурной суммы имеет тот же порядок, что и погрешность вычис-
n
ления плотности /(ф). В таких случаях говорят, что квадратурная формула численно устойчива.
Пример. Найти решение уравнения (1) при условии
du dr
= sin ф + фcos ф, - % < ф < л.
r=1
Точное решение этой задачи найдено в [10, с. 133] с помощью преобразований интеграла Шварца
и(г, ф) = Щг, ф)+С = 1т ^ - 1п(1 + г) + С, (19)
где 1п(1 + г) - ветвь логарифмической функции, принимающей на промежутке (-1, 1) действительной оси действительные значения.
Максимальная погрешность вычислений по точной формуле (19) и с помощью приближенного уравнения (18) при г = 0,1 + 0,2к, к = 0, 4;
3. Смирнов, В. И. Курс высшей математики: в 5 т. / В. И. Смирнов. М.: Наука, 1974. Т. 4. 336 с.
4. Беляев, Н. М. Методы теории теплопроводности: в 2 ч. / Н. М. Беляев, А. А. Рядно. М.: Высш. шк., 1982. Ч. 1. 327 с.
5. Бейтмен, Г. Высшие трансцендентные функции / Г. Бейтмен, А. Эрдейи. М.: Наука, 1967. 294 с.
6. Пыхтеев, Г. Н. Полилогарифмы, их свойства и методы вычисления / Г. Н. Пыхтеев, И. Н. Мелешко. Минск: Изд-во БГУ, 1976. 68 с.
7. Крылов, В. И. Приближенное вычисление интегралов / В. И. Крылов. М.: Наука, 1967. 500 с.
8. Мысовских, И. П. Интерполяционные кубатурные формулы / И. П. Мысовских. М.: Наука, 1981. 336 с.
9. Мелешко, И. Н. Квадратурные формулы с неотрицательными коэффициентами для сингулярных интегралов Коши / И. Н. Мелешко // Весщ АН БССР. Сер. ф1з.-мат. навук. 1989. № 3. С. 27-34.
10. Мелешко, И. Н. Специальные формулы для интегралов типа Коши и их приложения / И. Н. Мелешко. Минск: ВУЗ-ЮНИТИ, 1999. 197 с.
Поступила 07.10.2020 Подписана в печать 17.12.2020 Опубликована онлайн 30.03.2021
ф = — и п = 20 равна 0,002. Аналогично погрешность вычислений в тех же точках при п = 50 не превышает 0,00033, а при п = 100 она не больше чем 0,000083.
Численный эксперимент подтверждает эффективность и устойчивость квадратурной формулы.
ВЫВОДЫ
1. На основе известного точного решения краевой задачи Неймана для круга с помощью специальной квадратурной формулы для интеграла Дини сконструировано эффективное приближенное представление дилогарифмами решения эквивалентной вариационной задачи на допустимом множестве функций.
2. Получена равномерная в единичном круге оценка погрешности приближенной формулы. Проведен численный анализ квадратурной суммы на устойчивость.
ЛИТЕРАТУРА
1. Канторович, Л. В. Приближенные методы высшего анализа / Л. В. Канторович, В. И. Крылов. М.-Л.: Физ-матгиз, 1962. 708 с.
2. Лаврентьев, М. А. Методы теории функций комплексного переменного / М. А. Лаврентьев, Б. В. Шабат. М.: Наука, 1973. 736 с.
REFERENCES
1. Kantarovich L. V., Krylov V. I. (1962) Approximate Methods of Higher Analysis. Moscow-Leningrad, Fizmat-giz Publ. 708 (in Russian).
2. Lavrentiev M. A., Shabat B. V. (1973) Methods of the Theory of Functions of a Complex Variable. Moscow, Nauka Publ. 736 (in Russian).
3. Smirnov V. I. (1974) Higher Mathematics Course. Vol. 4. Moscow, Nauka Publ. 336 (in Russian).
4. Belyaev N. M., Ryadno A. A. (1982) Methods of the Theory of Heat Conduction. Part 1. Moscow, Vysshaya Shkola Publ. 327 (in Russian).
5. Bateman H., Erdelyi A. (1967) Higher Transcendental Functions. Moscow, Nauka Publ. 294 (in Russian).
6. Pykhteev G. N., Meleshko I. N. (1976) Polylogarithms, their Properties and Calculation Methods. Minsk, Publishing House of Belarusian State University. 68 (in Russian).
7. Krylov V. I. (1967) Approximate Calculation of Integral Methods. Moscow, Nauka Publ. 500 (in Russian).
8. Mysovskikh I. P. (1981) Interpolatory Cubature Formulas. Moscow, Nauka Publ. 336 (in Russian).
9. Meleshko I. N. (1989) Quadrature Formulas with Nonnegative Coefficients for Singular Cauchy Integrals. Vestsi Akademii Navuk BSSR. Ser. Fizika-Matematychnykh Na-vuk [Proceedings of the Academy of Sciences of BSSR. Physics and Mathematics Series], (3), 27-34 (in Russian).
10. Meleshko I. N. (1999) Special Formulas for Integrals of Cauchy Type and their Applications. Minsk, VUZ-UNITI Publ. 197 (in Russian).
Received: 07.10.2020 Accepted: 17.12.2020 Published online: 30.03.2021
Наука
итехника. Т. 20, № 2 (2021)