CC BY
12
ISSN 1810-0198. Вестник Тамбовского университета. Серия: естественные и технические науки
Том 23, № 123
2018
DOI: 10.20310/1810-0198-2018-23-123-466-472 УДК 519.651
О ТОЧНОМ РЕШЕНИИ ОДНОЙ ЗАДАЧИ ОПТИМИЗАЦИИ, ПОРОЖДЕННОЙ УРАВНЕНИЕМ ЛАПЛАСА
А. Н. Мзедавее. В. И. Родионов
ФГБОУ ВО «Удмуртский государственный университет» 426034, Российская Федерация, г. Ижевск, ул. Университетская, 1 E-mail: [email protected], [email protected]
Аннотация. Определено однопараметрическое семейство конечномерных пространств, состоящих из специальных двумерных сплайнов лагранжевого типа (параметр N связан с размерностью пространства). Уравнение Лапласа порождает в каждом таком пространстве задачу минимизации функционала невязки. Доказаны существование и единственность оптимальных сплайнов. Для их коэффициентов и невязок получены точные формулы. Показано, что с ростом N минимум функционала невязки есть величина О(iV-5), а специальная последовательность, состоящая из оптимальных сплайнов, фундаментальна. Ключевые слова: интерполяция; многомерный сплайн; многочлены Чебышёва
Введение
Работа продолжает исследования [1—4]. Уравнение Лапласа u£i + u^ = 0, заданное в прямоугольнике, заменой переменных приводится к виду auti + biut = 0 (в терминах новых переменных из квадрата П = [О, I]2 ). Пусть а > 0, b > 0, а непрерывные функции ipQ, срг, д0, Q\: [0,1] ¡Е К таковы, что
Ыо) = еь(о), = й(о), = еь(1), пО) =
существуют производные ^'(0), //(0), ^(1), fg(0), eff(l), ei(0), и
а^(0) + Ы'(°) = 0, а^(0)+Ы'(1) = 0, afig(l) + Ы'(0) = 0, + Ь^'(1) = 0.
Функция u = u(£,£), / П, являющаяся решением задачи
autt + bu^ = 0, u(0. 0 = u(l,О = уиСО, u(£, 0) = u(£, 1) = в1 (£),
прсдставима в виде и = и'0' + и*1 ■ + где
uW = uW(i,O = ^(0)(l 0(1 f) + W.(l)(l t)i + vi(0)i(l f) + №(l)*f
— билинейная функция, а функции u,;1) = u: 1-)(£, £) и u1'2i = u'2'(i, £) — решения задач
autt+bu^=0, 11(0,0 =Ро(0, u(1,0=PiC0, u(f, 0) = ft,(f), u(£, 1) = &(£), (1) autt + buii = 0, и(0,£) = «юСО. u(U) = V?i(0, u(£, 0) = i0(i), u(i, 1) = l(i), (2)
соответственно. Использованы следующие обозначения:
Pi(0 = 3€a+2f) + ivf(l)(ia £), ¿ = 0,1,
qi(t)= ¿fifO0)(i3 3t2+2i) + ^'(l)(i3 t), ¿=0,1,
— это кубические полиномы, a
№(0=^(0 w(o)(i О Pi(£), ^ = o, i,
-ft(i) = ft(t) ft(0)(l £) ft(l)i 4i(f), ¿ = 0,1,
— это функции, порожденные граничными функциями исходной задачи.
В работе обсуждается специальная задача оптимизации, порожденная задачей (1), а аналогичная задача, порожденная задачей (2), симметрична, — следует лишь осуществить замены а-*Е Ь. i^t i, Pi(0He qi((D ftiGHe ^¿(ф
1. Постановка задачи построения оптимального сплайна
Задача (1) порождает задачу поиска оптимального сплайна задачи
(3)
где <г(П) — это пространство, состоящее из допустимых сплайнов (см. ниже), зависящих от коэффициентов и\, и\, г = 1,..., 3N 1 (где N — параметр, отвечающий за количество узлов разностной схемы), и определенных в квадрате П. Пусть, далее,
п=Аг 1, т = ^, к =
а точки Лд) / П таковы, что ц = ¿т, г = 0,1,..., ЗМ, з = 0,1, 2, 3.
Массив и)\_1 ® = 0,1,..., ЗЛ', у = 0,1, 2, 3, называется допустимым, если:
1) иг0 = Яо(^), ^ = для всех г = 0,1,... , ЗЛГ;
2) = ^(ЛД и^р^А,) для з = 0,1,2,3. Одномерные интерполяционные многочлены Лагранжа
С «
^к(С)=
«=0Д,2,3: афк
К а
с /М, к = 0,1, 2,3
(такие, что = для всех к, ц = 0,1, 2, 3, где — символ Кронекера) и допу-
стимый массив и* [, г = 0,1,..., ЗN, = 0,1, 2,3, порождают семейство полиномов
з з
<Л.,Ч)=П П 3,7}/Ш, к = 1.....ЛГ.
г=о ¿=о
Пусть РЧ^О^М, где я = | Зк + З, г] = | = 3£ Тогда А,) =
для всех £ = 1,...,ЛГ и = 0,1,2,3. Следовательно, полином Р*(0©является двумерным интерполяционным многочленом Лагранжа, определенным в 16 узлах полосы П* = }(*,£) / П: тзь-з ^ £ ^ Тз^, 0 ^ С ^ 1<- Значит, определена непрерывная функция П е М, такая, что = при (£,0 / П*. Другими словами,
всякий допустимый массив порождает бикубический сплайн, который мы называем аппроксимирующим. Разнообразие таких сплайнов определяется лишь наборами чисел и12, г = 1,...,ЗЛГ 1. Это означает, что аппроксимирующие сплайны образуют конечномерное пространство. Обозначим его <т(П) = ег№ (П).
2. Конечные разности аппроксимирующих сплайнов
Всякий допустимый массив и* [, г = 0,1,..., 3N, ] = 0,1, 2,3, порождает термы
х1 = и1п
и\ и{ + г4, у' = ид ЗЦ + 3^4 иг3, г = 0,1,..., ЗЛГ,
X* = хзк~г
Зк-2 хЗк-1+хЗк
Х* = у
Зк-З
X
у3*"2
у^+у3*,
10 — X ук = у3к-3
Зхзк~2+гх-
3уЫ-2 + 3уЗк-1
3^—1
X
3 к
у3\
(4)
к= 1,...,ЛГ,
и граничные элементы
V),
К * К I К
Ро = 20 +
% = 4
4 = \ Г1] ДОО(2)(Т3),_3, т) + Д&(2)(тз*-з, г) + А^2\тзк_21 т) + дЛтзк-а, т) {,
- и,к =|6>-1]Ак(3)(тзк-з,г) + д1(3)(тзь_з)г){, -г* ^Г1]^«^) д1С2)(Г3,_з,т) + Д1(2)(Т3,_2,Т) А/1^^)!, Л = | Г1] Дуз)(тз*-з,т) дЬ)(Тз,_3,г) {, „
& = 1, . . . , 71,
1,...,71, Г}' = у-е = со1 ^? = со1
е=Р1+А е=Рк0
Лг-—|— 1
Г}к = рк+х
к
9си
о* У1)
злт
X
ш
£
(5)
(б)
Применяем обозначения для конечных разностей ( I)1 СД /(£ + гт).
¿=о
3. Основные результаты Утверждение 1. Для любого и / ая (П) имеет, место равенство
За2
N
^ ^ П Л
¿=1
= 402] Х3к~3 + хзк + 22к( 60 ( 1+0 ) ] Х3к~3 + г3* + 24 + ^ (3 +50+302 ) ] Хк +§ 02] я3*"3 я3,г + 2юк ( | 0( 5+9) ] х3*"3 х3* + 2и£ {У0* + ^ ( 21 + 70 + 02 ) ] У0* {2 +1202] у3*"3 + у3* + 2г* {= £в (1+50)] у3*"3 + у3* + 2гк {Хк + % (1+70+21^2 )]Хк ( +4 д2]узк~3 узк + 2и1к{2 + в)]узк~3 узк+ 2™* {¥? + §-„ (5+76+Ы2 )]¥?(.
Утверждение 2. Коэффициенты и\, иг2, г = 1,... , 3N 1, оптимального аппроксимирующего сплайна и / сгу (П) вычислимы по формулам (4) через величины хзк, узк, к = 0,1,..., ЛГ, X*, Г*, Л'^, Ук, к = 1,..., ЛГ, удовлетворяющие системе уравнений
сзк-з + 2уХзк + хзк+3 + = О, А: = 1,... ,71, Хк = 7о ] хзк~3 + хзк + 24 {, = ¿о ] хзк~3 хзк + 2ъик{, ■зк~3 + 2хузк + узк+3 + ук=0, Хк = Ъ ] узк~3 + узк + 2гк {, Ук = б1 ] узк~3 узк + 2тк {,
к = 1,... ,п, к =
к = 1,...,ЛГ.
(7)
В системе фигурируют константы ж0, хЗР1, у0, узл (известные априори, см. (4)),
с*о =
1 1+ 02
«1 =
2 3 + 50 + 302'
3 3 + 3502
10 1 + 70 + 2102'
1 35 + 302
30 21+70+02'
. Ю 0(1+0)
7о =
А 4
1+02
2 Б + 70 + 502'
Т ^ «о + Дэ У а-'
«О Ро
71 =
X =
3 3+50 + 302' 14 0(1+50) У 1+70+2102'
«1 А
. . 14 0(5+0) оо = ^г
9 21 + 70+02'
_ . 70 0(1+0)
=тг
9 5+70+502'
(справедливо а0 > /30, «1 > у > 2, х > 2) и граничные элементы (см. (5))
'? = (1+х)]^ + ^+1{+(1 х)]^
Щ
¡4-1
т
£+1
Первая совокупность уравнений (7) имеет самостоятельный характер: ее уравнения связывают между собой лишь переменные вида х3т. Матрица совокупности имеет трех-диагоиальный вид с доминирующей главной диагональю (так как у > 2), поэтому совокупность имеет единственное решение (его легко найти методом прогонки). После этого из второй и третьей совокупностей в (7) явно вычисляются все значения и У^. Аналогично решаются четвертая (где х > 2 ), пятая и шестая совокупности в (7). Полученные значения позволяют, в конечном счете, найти величины и\, и12, 1 = 1,..., ЗN 1 (см. определения (4)), порождающие оптимальный аппроксимирующий сплайн.
Утверждение 3. Единственным решением первой и четвертой совокупностей системы (7) являются числа
1 1 Ви(у) х° + Вы(у) хЗЛ7 + [] Вы{у) ,
хзк =
узк =
ип(у)
1
к = 1,..., 71,
¿=1 п
Ды(х) у0 + Вкп(х) уш + Ц Вкг(х) ^ , к = 1, ... ,71.
¿=1
В представлении используются многочлены Чебышсва 2-го рода согласно
[5, с. 96] неравенства у > 2, х > 2 влекут [/„(у) з{= 0, Ип(х) з{= 0[ и функциональные матрицы В(ф= Вы(ф[, Ви(ф=( 1)] [4-1 ((Ф+ Е/„-Ь(№-1 (Ф{,
к,{ = 1,... ,71. Применяем символ такой, что ¿^=0 (при к < г), а иначе 5^ = 1.
Утверждение 4. Для минимума JN = min J ((J) функционала (3) имеет место точная формула
лг" 16,V J ип(у) + Е/П(х) ■
В представлении используются многочлены Чебышсва 1-го рода Tn(Q) порождающие функциональные матрицы В(ф= = 0,1,..., JV, такие, что
ШФ= ( }¿Uk ткфгп-im € ?w<№«D[.
Векторы £ и г] L- это граничные элементы (6), а для записи квадратичных форм применяется скалярное произведение )OCÜ) пространства R1+A .
Утверждение 5. Если до, Qi /С5[0,1], mo JN =0(TV-5).
Зафиксируем натуральное N0 ^ 2. Последовательность Nm = N0 2т, т = 1,2,..., порождает функциональную последовательность } йт , где um : П t 1 — это оптимальный аппроксимирующий сплайн в пространстве о"л7т(П).
Утверждение 6. Если Qo,Qi / С5[0,1], то последовательность } ит фундаментальна по норме пространства Ь2(П).
СПИСОК ЛИТЕРАТУРЫ
1. Родионов В.И. О применении специальных многомерных сплайнов произвольной степени в численном анализе // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2010. Вып. 4. С. 146-153.
2. Родионов В. И. Об одном методе построения разностных схем // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2013. Т. 18. Вып. 5. С. 2656-2659.
3. Родионов В.И., Родионова Н.В. О решении двух задач оптимизации, порожденных простейшим волновым уравнением // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2015. Т. 20. Вып. 5. С. 1406-1409.
4. Мзедавее А.Н., Родионов В.И. О решении одной задачи оптимизации, порожденной уравнением Лапласа // Современные методы прикладной математики, теории управления и компьютерных технологий: материалы Междунар. науч. конф. Воронеж, 2017. С. 252-255.
5. Суетин П.К. Классические ортогональные многочлены. М.: Наука, 1976. 328 с.
Поступила в редакцию 17 апреля 2018 г.
Прошла рецензирование 22 мая 2018 г.
Принята в печать 19 июня 2018 г.
Конфликт интересов отсутствует.
Мзедавее Асаад Насер, Удмуртский государственный университет, г. Ижевск, Российская Федерация, аспирант, кафедра информатики и математики, e-mail: [email protected]
Родионов Виталий Иванович, Удмуртский государственный университет, г. Ижевск, Российская Федерация, кандидат физико-математических наук, доцент, зав. кафедрой информатики и математики, e-mail: [email protected]
Для цитирования: Мзедавее А.Н., Родионов В.И. О точном решении одной задачи оптимизации, порожденной уравнением Лапласа // Вестник Тамбовского университета. Серия: естественные и технические науки. Тамбов, 2018. Т. 23. № 123. С. 466-472. ЕЮ1: 10.20310/1810-0198-2018-23-123-466-472
DOI: 10.20310/1810-0198-2018-23-123-466-472
ON EXACT SOLUTION OF OPTIMIZATION TASK GENERATED BY THE LAPLACE EQUATION
A. N. Mzedawee, V. I. Rodionov
Udmurt State University 1 Universitetskaya St., Izhevsk 426034, Russian Federation E-mail: [email protected], [email protected]
Abstract. A one-parameter family of finite-dimensional spaces consisting of special two-dimensional splines of Lagrangian type is defined (the parameter N is related to the dimension of the space). The Laplace equation generates in each such space the problem of minimizing the residual functional. The existence and uniqueness of optimal splines are proved. For their coefficients and residuals, exact formulas are obtained. It is shown that with increasing N, the minimum of the residual functional is O(N_5), and the special sequence consisting of optimal splines is fundamental. Keywords: interpolation; multivariate spline; Chebyshev's polynomials
REFERENCES
1. Rodionov V.I. O primenenii spetsial'nykh mnogomernykh splaynov proizvol'noy stepeni v chislennom analize [On application of special multivariate splines of any degree in the numerical analysis]. Vestnik Udmurtskogo universiteta. Matematika. Mekhanika. Komp'yuternye nauki - The Bulletin of Udmurt University. Mathematics. Mechanics. Computer Science, 2010, no. 4, pp. 146-153. (In Russian).
2. Rodionov V.I. Ob odnom metode postroeniya raznostnykh skhem [About method for constructing difference schemes]. Vestnik Tambovskogo universiteta. Seriya Estestvennye i tekhnicheskie nauki - Tambov University Reports. Series: Natural and Technical Sciences, 2013, vol. 18, no. 5, pp. 2656-2659. (In Russian).
3. Rodionov V.I., Rodionova N.V. O reshenii dvukh zadach optimizatsii, porozhdennykh pros-teyshim volnovym uravneniem [On solution of two optimization problems generated by the simplest wave equation]. Vestnik Tambovskogo universiteta. Seriya Estestvennye i tekhnicheskie nauki -Tambov University Reports. Series: Natural and Technical Sciences, 2015, vol. 20, no. 5, pp. 1406-1409. (In Russian).
4. Mzedavee A.N., Rodionov V.I. O reshenii odnoy zadachi optimizatsii, porozhdennoy urav-neniem Laplasa [On solution of optimization tasks generated by the Laplace equation]. Materialy Mezhdunarodnoy nauchnoy konferentsii «Sovremennye metody prikladnoy matematiki, teorii up-ravleniya i komp'yuternykh tekhnologiy» [Proceedings of the International Scientific Conference "Modern Methods of Applied Mathematics, Control Theory and Computer Technologies"]. Voronezh, 2017, pp. 252-255. (In Russian).
5. Suetin P.K. Klassicheskie ortogonal'nye mnogochleny [Classical Orthogonal Polynomials]. Moscow, Nauka Publ., 1976, 328 p. (In Russian).
472
А. Н. Мзеlцавее, B.H. POIHOHOB
Received 17 April 2018 Reviewed 22 May 2018 Accepted for press 19 June 2018 There is no conflict of interests.
Mzedawee Asaad Naser, Udmurt State University, Izhevsk, the Russian Federation, Post-Graduate Student, Department of Informatics and Mathematics, e-mail: [email protected]
Rodionov Vitalii Ivanovich, Udmurt State University, Izhevsk, the Russian Federation, Candidate of Physics and Mathematics, Associate Professor, Head of the Department of Informatics and Mathematics, e-mail: [email protected]
For citation: Mzedawee A.N., Rodionov V.I. O tochnom reshenii odnoi zadachi optimizatsii, porozhdennoi uravneniem Laplasa [On exact solution of optimization task generated by the Laplace equation]. Vestnik Tambovskogo universiteta. Seriya: estestvennye i tekhnicheskie nauki - Tambov University Reports. Series: Natural and Technical Sciences, 2018, vol. 23, no. 123, pp. 466-472. DOI: 10.20310/1810-0198-2018-23-123-466-472 (In Russian, Abstr. in Engl.).
