Владикавказский математический журнал 2023, Том 25, Выпуск 1, С. 105-111
УДК 517.518.8
DOI 10.46698/m0485-4484-9134-k О МНОГОЧЛЕНАХ НАИЛУЧШЕГО ПРИБЛИЖЕНИЯ СЕГМЕНТНЫХ ФУНКЦИЙ
А. Ю. Трынин1'2
1 Саратовский национальный исследовательский государственный университет им. Н. Г. Чернышевского, Россия, 410012, Саратов, ул. Астраханская, 83; 2 Московский центр фундаментальной и прикладной математики, Московский государственный университет им. М. В. Ломоносова, Россия, 119991, Москва, Ленинские горы, 1 E-mail: tayu@rambler.ru
Аннотация. Предложен алгоритм поиска многочлена наилучшего приближения для непрерывной многозначной сегментной функции, заданной на совокупности не пересекающихся отрезков
X = (Un.UK AiО и (ULo xk таких что (LjWai Ai0 n (ULo = 0, где не пересекающиеся отрезки [aj. ,bj.] и точки xk принадлежат ограниченному отрезку [A, B] С К. Считаем, что функции f1 и f2 непрерывны на множестве X, и всюду на X значение функции f1(x) не превосходит значение функции f2(x). Оператор, ставящий в соответствие каждому x £ X отрезок [(x, f1 (x)), (x, f2(x))], будем называть сегментной функцией F(x), заданной на X. В силу непрерывности функций f1 и f2 сегментная функция F является h-полунеперывным отображением сверху. Многочлен Pm = ^Г=0 aixz наилучшего приближения в метрике Хаусдорфа на множестве X сегментной функции F с вектором коэффициентов a = (a0, a1,... , am) £ Rm+1 есть решение экстремальной задачи ^Qina^Rm+i ^Qaxx ex max(Pm(x) — f1 (x),f2(x) — Pm(x)). Методами конструктивной теории функций показано, что для любых непрерывных на X функций f1(x) < f2 (x) существует многочлен наилучшего приближения в хаусдорфовой метрике h-полунепрерывной сверху на множестве X сегментной функции F(x). Предложен алгоритм описания множества коэффициентов a многочленов наилучшего приближения сегментной функции. Получены необходимые и достаточные условия единственности многочлена наилучшего приближения сегментной функции. Приведены результаты численных экспериментов, реализованных с помощью предложенного алгоритма.
Ключевые слова: наилучшее приближение функции, аппроксимация многочленами, сегментная функция.
AMS Subject Classification: 65D15, 26E25.
Образец цитирования: Трынин А. Ю. О многочленах наилучшего приближения сегментных функций // Владикавк. мат. журн.—2023.—Т. 25, вып. 1.—C. 105-111. DOI: 10.46698/m0485-4484-9134-k.
1. Введение
В [1] подробно изложены основы теории многозначных отображений, имеющих широкий спектр прикладных возможностей, так как использование многозначных отображений позволяет в процессе реализации алгоритма сразу оценивать погрешность приближения решения. При этом подход, разработанный автором монографии [1], не предполагает использования интервальной арифметики, что снимает вопрос о построении поля интервалов или его аналогов. Вместе с тем, теория, разработанная в [1], в значительной степени справляется с многими основными задачами интервальной математики.
© 2023 Трынин А. Ю.
Пусть
(«4 \/™ \ ( «1 \/™\
UKi Ai] U U жк , U j ,bji] П IJxH = 0' (1) jl =0 ) Vfc=0 / Vji=0 / \fc=0 /
где не пересекающиеся отрезки [aj1, bj1 ] и точки принадлежат ограниченному отрезку [A, B] С R. Считаем, что функции /1 и /2 непрерывны на множестве X, и всюду на X выполняется соотношение /i(x) ^ /2(х). Оператор, ставящий в соответствие каждому х € X отрезок [(х, /1(х)), (х, /2(х))], вслед за авторами работ [2-5] будем называть на X сегментной функцией F(х). В силу непрерывности функций /1 и /2 сегментная функция F является h-полунеперывным отображением сверху (см. [1, §6]). Многочлен Рт наилучшего приближения в метрике Хаусдорфа на множестве X сегментной функции F с вектором коэффициентов а = (а0, а1,..., am) € Rm+1 есть решение экстремальной задачи
min maxmax (Рт(х) - Л(х), /2(х) - Рт(х)),
«eRm+1xex
т . (2) Рт(х) = (а, (1,х,... ,хт)) = £ а.хг.
.=0
В случае /1 = /2 задача (2) сводится к классической теореме П. Л. Чебышева для множества X.
В публикациях [2-5] приводятся необходимые и достаточные условия существования и единственности многочлена наилучшего приближения (2) для случая, когда X представляет собой либо отрезок, либо конечный набор точек. Эти работы используют методики доказательства теоремы П. Л. Чебышева и выпуклого анализа. В этом случае, даже когда X представляет собой конечный набор из n точек, количество систем, которые необходимо решить для поиска многочлена степени m, растет как Ст+2. Кроме того, многочлен, наилучшим образом приближающий сегментную функцию, одновременно представляет собой многочлен наилучшего приближения некоторой классической функции. А из результатов исследований в работах [6, 7] следует, что независимо от метода решения задачи (2) его численная реализация возможна только для небольших степеней и хороших функций, так как модули коэффициентов многочленов наилучшего приближения непрерывных функций с ростом степени многочленов могут расти быстрее геометрической прогрессии. Тем не менее в данной работе предложен другой подход к изучению этой проблемы на основе классического анализа и теории функций. Здесь приведен также ряд примеров численных экспериментов, позволяющих убедиться в конструктивности предложенного метода исследований.
2. Основной результат
Определим функцию ограниченной вариации на [A, B]: x (
Мх) = / Х(0 dC + £ 1, где х(х) = ° ^ х € X; Bn X (3)
A xl^x 1°, пРи х € [A,B] \X,
где X определено в (1), и будем рассматривать ее как меру в интеграле Стилтьеса. Сделав обозначение
ОТ)/ t\ ^ i /1(х) + /2(х) М(а, х, т, j) = 2_^ сцх--—-—.
для каждого p € N определим систему уравнений в
[ х* ¡Ща, х, т, /) аг<^ {рЩа, х, т, /)) + \ ?
А J (4)
( , ... рМ(а,х,т, /) \ , , .
х агс!? (рт(а, х,т, т)) Н—тг-т-,-^- аа(х) =0, г = 0.... ,т,
V 4 4 р2^(а,х,т,/)2 +1) ' ' ' ' '
относительно неизвестных а € Мт+1. Обозначим через Ер замыкание множества решений системы (4) при каждом р € N и
E = П U EP- (5)
k=1 p=k
Теорема 1. Для любых непрерывных на X функций fi(x) ^ f2(x) существует многочлен наилучшего приближения (2) в хаусдорфовой метрике h-полунепрерывной сверху на множестве (1) сегментной функции F (x). Он единственен тогда и только тогда, когда множество (5) одноточечное в Rm+1. Любой вектор a € E имеет координаты, являющиеся коэффициентами многочлена наилучшего приближения.
< Ввиду того, что для произвольной функции g € C[Y] справедливо соотношение
2
lim g(x) — aicAg(pg(x)) = \g{x)\,
р^те п
из которого следует (см., например, [8, гл. 1, §3]) равенство в —
{^arctg(^))} dß(x)\ = vraimax\g(x)\ = \\g\\c[Y],
^ A '
величина уклонения многочлена Pm от сегментной функции F в метрике Хаусдорфа на множестве X может быть представлена следующим образом:
( в
p(F,Pm) = lim / < х, m, /) — arctg (p&(a, x, m, /)) р^те \ J у п
A
Подынтегральная функция, как сумма (возможно нестрого) выпуклых функций с неотрицательной второй производной, в (6), как функция переменных а*, является выпуклой (возможно нестрого) на Мт+1. Поэтому при каждом р € N множество векторов коэффициентов экстремальных многочленов таких, что подынтегральная функция наименее уклоняется от нуля по норме пространства Ь2р[Х], является, в силу необходимого условия экстремума, множеством решений систем
/ Б
¿¡(¡{щ^т.П^Шыт,,»
+ Ш-Ш}2^Кх)У = 0, г = 0,...,т, (7)
или эквивалентных системам (7) систем (4). Непустота множества решений этих систем следует из необходимого условия экстремума, непрерывной дифференцируемости неотрицательной подынтегральной функции по всем переменным и ее неограниченного возрастания по любому направлению пространства Rm+2 при стремлении аргумента к бесконечности.
Так как функции fi и f2 непрерывны на множестве (1), то предельный переход в интеграле Стилтьеса обосновывается с помощью теоремы 2 из [9, гл.8, §7].
Коэффициенты экстремальных многочленов, наименее уклоняющихся от сегментной функции F(x) в метрике Хаусдорфа, в силу замкнутости множества нулей непрерывной функции составляют множество (5) точек прикосновения всевозможных пределов подпоследовательностей решений систем (4).
По теореме Кантора [8, гл. 5, § 1] множество (5) не пусто. >
3. Результаты некоторых численных экспериментов
В этом параграфе приведем результаты численных экспериментов, организованных с помощью утверждения теоремы 1. Решение систем (4) осуществлялось с помощью некоторой адаптации стандартных методов наискорейшего спуска или покоординатного спуска к решению изучаемой экстремальной задачи. На рисунке 1 изображена реализация примера, демонстрирующего, что рассматриваемый многочлен первой степени наилучшего приближения отличается от Чебышевского наилучшего приближения середин отрезков сегментной функции. На рисунке 2 график одного из многочленов двенадцатой степени наилучшим образом аппроксимирущих в метрике Хаусдорфа ту же сегментную функцию, что и на рисунке 1.
4.0 -
3.5--
3.0--
2.5--
2.0--
1.5--
1.0 -
0.5--
0.0--
0.0!
Рис. 1. n = 2, m = 1. Рис. 2. n = 2, m = 12.
Рис. 3 демонстрирует возможность предлагаемого метода обрабатывать большие массивы данных. В качестве модели сегментной функции рассматривается пара fi (xk) = sin7xk — 1, f2(xk) = sin7xk + 0.1xk + 1, где 0 ^ k ^ 1000. В качестве многочлена наилучшего приближения взята линейная функция.
Еще один очевидный пример наилучшего Хаусдорфова приближения сегментной функции fi (xk) = sin xk — 1, f2 (xk) = sin xk + 1, где 0 ^ k ^ 3, многочленом первой степени, который вычислен с помощью теоремы 1, приведен на рис. 4.
О 200 400 600 800
Рис. 3. п = 1000, т = 1.
Рис. 4. п = 3, т =1.
Рис. 5 содержит график реализации приближения сегментной функции ^(0) = [(0,1), (0,1)], ^(1) = [(1, 0), (1, 2)] многочленами наилучшего Хаусдорфова приближения первого порядка, которые вычислены с помощью теоремы 1 два раза. Это та ситуация, когда многочлен наилучшего приближения не единственен. В этом случае множество Е представляет собой отрезок [(2, — 1), (0,1)]. Адаптация метода наискорейшего спуска, о которой шла речь в начале этого параграфа, позволяет приближенно вычислять коэффициенты многочленов наилучшего приближения с границы множества Е, выбирая различные начальные итерации. Графики двух таких многочленов и получены в результате описанного численного эксперимента.
Рис. 5. Многочлены с коэффициентами с границы множества Е,
п =1, т =1.
Литература
1. Половинкин Е. С. Многозначный анализ и дифференциальные включения.—М.: Физматлит, 2014.—524 с.
2. Выгодчикова И. Ю., Дудов С. И., Сорина Е. В. Внешняя оценка сегментной функции полиномиальной полосой // Журн. вычисл. матем. и мат. физ.—2009.—Т. 49, № 7.—С. 1175-1183.
3. Выгодчикова И. Ю. О единственности решения задачи наилучшего приближения многозначного отображения алгебраическим полиномом // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика.—2006.—Т. 6, № 1-2.—С. 11-19. Б01: 10.18500/1816-9791-2006-6-1-2-11-19.
4. Выгодчикова И. Ю. Об аппроксимации многозначного отображения алгебраическим полиномом с ограничениями // Изв. вузов. Мат.—2015.—№ 2.—С. 30-34.
5. Выгодчикова И. Ю. О приближении двузначной функции алгебраическим полиномом // Изв. вузов. Мат.—2016.—№ 4.—C. 8-13.
6. Stafney J. D. A permissible restriction on the coefficients in uniform polynomial approximation to C[0,1] // Duke Math. J.—1967.—Vol. 34, № 3.—P. 393-396. DOI: 10.1215/S0012-7094-67-03443-6.
7. Хавинсон С. Я. Допустимые величины коэффициентов многочленов при равномерной аппроксимации непрерывных функций // Мат. заметки.—1969.—Т. 6, № 5.—С. 619-625.
8. Люстерник Л. А., Соболев В. И. Элементы функционального анализа.—М.: Наука, гл. ред. физ.-мат. лит-ры, 1965.—520 с.
9. Натансон И. П. Теория функций вещественной переменной.—М.: Наука, 1974.—480 с.
Статья поступила 13 января 2022 г. Трынин Александр Юрьевич
Саратовский национальный исследовательский государственный университет им. Н. Г. Чернышевского, профессор кафедры диф. уравнений и мат. экономики РОССИЯ, 410012, Саратов, ул. Астраханская, 83;
Московский центр фундаментальной и прикладной математики, Московский государственный университет им. М. В. Ломоносова, РОССИЯ, 119991, Москва, Ленинские горы, 1 E-mail: tayu@rambler.ru https://orcid.org/0000-0002-6304-4962
Vladikavkaz Mathematical Journal 2023, Volume 25, Issue 1, P. 105-111
ON THE BEST POLYNOMIALS APPROXIMATION OF SEGMENT FUNCTIONS
Trynin, A. Yu.1'2
1 Saratov State University, 83 Astrakhanskaya St., Saratov 410012, Russia;
2 Moscow Centre for Fundamental and Applied Mathematics, Lomonosov Moscow State University, 1 Leninskie Gory, Moscow 119991, Russia E-mail: tayu@rambler.ru
Abstract. An algorithm for finding the best approximation polynomial for a continuous multivalued segment function defined on a set of segments X is proposed, where X = (Un11=o[aii , bj0 ^ (Un=oXk) with (LjUK ,bj1 ]) n (Un=0 Xfc) = 0. The disjoint segments [aj ,bj1 ] and points xk belong to a bounded segment [A, B] C R. We assume that the functions f1 and /2 are continuous on the set X, and everywhere on X the value of the function f1 (x) does not exceed the value of the function /2(x). The operator assigning to each x £ X the segment [(x,f1(x)), (x,/2 (x))] will be called the segments function F(x) defined on X. Since the functions f1 and /2 are continuous, the segments function F is an upper h-semicontinuous mapping. The polynomial Pm = ^m=0 aix1 of the best approximation in the Hausdorff metric on the set X of a segment function F with a vector of coefficients a = (a0,a1,.. . ,am) £ Rm+1 is a solution to the extremal problem miHa^Rm+1 maxx ex max(Pm(x) — f1(x), /2 (x) — Pm(x)). It is shown by methods of constructive function theory that, for any functions /1(x) < /2(x) continuous on X, there exists some polynomial of best approximation in the Hausdorff metric as the segment function F(x) is upper h-semicontinuous on X. An algorithm for describing the set E of coefficients a of polynomials of the best approximation of a segment function is proposed. Necessary and sufficient conditions for the uniqueness of the polynomial of best approximation of the segment function are obtained. The results of numerical experiments carried out using the proposed algorithm are presented.
Key words: best approximation of functions, polynomial approximation, segment function.
AMS Subject Classification: 65D15, 26E25.
For citation: Trynin, A. Yu. On the Best Polynomials Approximation of Segment Functions, Vladikavkaz Math. J., 2023, vol. 25, no. 1, pp. 105-111 (in Russian). DOI: 10.46698/m0485-4484-9134-k.
References
1. Polovinkin, E. S. Mnogoznachnyy analiz i differentsial'nye vklyucheniya [Multivalued Analysis and Differential Inclusions], Moscow, Fizmatlit, 2014, 524 p. (in Russian).
2. Vygodchikova, I. Yu., Dudov, S. I. and Sorina, E. V. External Estimation of a Segment Function by a Polynomial Strip, Computational Mathematics and Mathematical Physics, 2009, vol. 49, no. 7, pp. 1119-1127. DOI: 10.1134/S0965542509070057
3. Vygodchikova, I. Yu. About the Only Solution in the Problem of the Best Plural Reflection's Approximation by Algebraic Polynomial, Izvestiya of Saratov University. New Series. Series: Mathematics. Mechanics. Informatics, 2006, vol. 6, no. 1-2, pp. 11-19 (in Russian). DOI: 10.18500/18169791-2006-6-1-2-11-19.
4. Vygodchikova, I. Yu. On Approximation of Multivalued Mapping by Algebraic Polynomial with Constraints, Russian Mathematics, 2015, vol. 59, no. 2, pp. 25-28. DOI: 10.3103/S1066369X15020048.
5. Vygodchikova, I. Yu. Approximation of Double-Valued Function by an Algebraic Polynomial, Russian Mathematics, 2016, vol. 60, no. 4, pp. 5-9. DOI: 10.3103/ S1066369X16040022.
6. Khavinson, S. Ya. Permissible Values of Coefficients of Polynomials in Uniform Approximation of Continuous Functions, Mathematical Notes of the Academy of Sciences of the USSR, 1969, vol. 6, no. 5, pp. 834-838. DOI: 10.1007/BF01101413.
7. Lyusternik, L. A. and Sobolev, V. I. Elementy funktsional'nogo analiza [Elements of Functional Analysis], Moscow, Nauka, 1965, 520 p. (in Russian).
8. Natanson, I. P. Teoriya funktsiy veshchestvennoy peremennoy [Theory of Functions of a Real Variable], Moscow, Nauka, 1974, 480 p. (in Russian).
Received January 13, 2022
Alexander Yu. Trynin Saratov State University,
83 Astrakhanskaya St., Saratov 410012, Russia, Professor;
Moscow Centre for Fundamental and Applied Mathematics, Lomonosov Moscow State University, 1 Leninskie Gory, Moscow 119991, Russia E-mail: tayu@rambler.ru https://orcid.org/0000-0002-6304-4962