УДК 517.9
О ЧИСЛЕННОЙ ОЦЕНКЕ ПРЕДЕЛОВ МАКСИМАЛЬНЫХ СРЕДНИХ ДЛЯ ПЕРИОДИЧЕСКИХ
ФУНКЦИЙ
© 2013 А.Н. Лепилов1
Предложена схема определения погрешности численного метода вычисления предела максимального среднего для периодической функции. Рассмотрен пример вычисления предела максимального среднего.
Ключевые слова: предел максимального среднего, дифференциальное включение, периодическая функция.
Введение
Данная работа посвящена практической реализации численного метода вычисления пределов максимальных средних, предложенного в [1], и по существу является ее продолжением. Под практической реализацией имеется в виду как непосредственное вычисление значения численным методом, так и оценка погрешности этого вычисления.
1. Основные понятия
Будем говорить, что функция f содержится в классе функций F, если f : D ^ R, D = R х Rm, (t, y) ^ f (t, y), y = (ylym); T-периодическая по любой переменной t,yi,... ,ym; f G C4(D,R).
Для функции f G F рассмотрим предел максимального среднего
1 ¡-to+A
Mf = lim sup- / f (t,Y(t)) dt, (1)
7 Л Jt0
где точная верхняя грань вычисляется по всем решениям дифференциального включения
Y G G, Y(to) = yo, (2)
G = [ai, bi] х ... х [am, bm] С Rm — параллелепипед, 0 < ai <bi, i = 1, 2,... ,m.
Множество решений задачи (2) в смысле Каратеодори, определенных на промежутке [to, те), обозначим r(to,yo).
1 Лепилов Александр Николаевич (lepilov_aleksand@mail.ru), кафедра уравнений математической физики Самарского государственного университета, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
Предел максимального среднего (1) существует и не зависит от начальных условий, то есть можно считать to = 0, более того, существует и оптимальное решение задач (1), (2) [2, теорема 1].
Рассмотрим также максимальное среднее на отрезке [0, Д]
1 СЛ
Mf = sup sup - f (t,Y(t)) dt, (3)
yoEK 7er(0,yo) Д Jo
где K = K(0, T), K(yo, T) = [yoi, yoi + T] x ... x [yom, yom + T] С Rm — куб.
Задача вычисления предела максимального среднего (1) с заданной точностью может быть заменена задачей вычисления максимального среднего (3). При этом справедлива оценка предела максимального среднего (1) [2, теорема 2]
Mf — £n < M f < Mf.
Здесь £(n) = 2t0C f/Д, Д = nT, n G Z, постоянная Cf > 0 такая, что f (t,y) ^ Cf для любых (t,y) G D, n ^ to/T, to > 0 такое, что ToG содержит некоторый куб K(zo,T), zo G Rm зависит от to, to = max {T/(bi — aj)}.
Таким образом, сосредоточимся на определении максимального среднего (3), которое произведем численно.
2. Численный метод и нахождение погрешности
Зафиксируем Д = пТ. Обозначим ГЛ(0,уо) как сужение множества всех решений Г(0, уо) на отрезок [0, Д]. Решение задачи (3) существует [2], то есть существует оптимальная пара (у™ах, 7тах(Ь)), Уо"ах е К, 7тах(Ь) е ГЛ(о, у™ах), при которой достигается МЛ.
Оптимальное решение задачи (3) будем находить по принципу максимума Понтрягина [3], решая следующую задачу Коши на отрезке [0, Д]:
д
Р] = - дГ/ (*,!), Р] (Д) = 0,
г дъ (4)
аз, если рз < 0, , . . .
^ = \ ь], если р] > 0, Ъ(Д) = узЛ ] = 1,---,т-
Для численного решения задач (3), (4) введем на отрезке [0, Д] равномерную сетку Лт = {Ь0,Ь1,... } с шагом т = Д/Q, Q е Z, Д = Ь0 > Ь1 > ... > Ьд =0. Начальное условие ул в (4) в силу Т-периодичности функции / будем брать из куба К, на котором также введем равномерную сетку О^ С К с шагом Н > 0 по каждой из координат.
Приближенное решение задачи (3) будем искать методом перебора. Для каждого значения ул е О^ находим решение задачи (4) и соответствующее значение среднего [1]. Выбирая среди всех полученных средних максимальное, принимаем его за приближенное значение максимального среднего (3). Обозначим его БЛ. Пусть 5 > 0 и множество А$ = {Ь е [0, Д] : 3], 1 ^ ] ^ т, такое, что
Р (^(т < 5}.
Сформулируем условие, накладываемое на функцию /.
Условие 1. 3 целое N, V 5 > 0 и V7(Ь) е ГЛ(0, у0) 3 к : М+ ^ М+, к(5) ^ 0 при 5 ^ 0 такая, что существует конечная система интервалов (е^,^), г = 1,...,Nl,
N1
N1 ^ N, (¿1 — е^) ^ к(5), покрывающая множество .
г=1
Приведем теорему об оценке погрешности приближенного вычисления предела максимального среднего Mf. Предполагаем, что Ymax(t) в момент времени t = Д отстоит от численно найденного решения по каждой из координат по модулю не более, чем на h.
Теорема 1 [1]. Пусть функция f £ F, и выполнено условие 1. Тогда имеет место следующая оценка:
Mf — Sf\ < £i(n) + £2(т) + £з(к(6), h),
где £\(n) = 2roCf /(nT) - теоретическая погрешность; £2(т) - погрешность интегрирования (например, для метода Симпсона [4] £2(т) = т4C4/2880, C4 - постоянная, зависит от функции f); £з(к,К) = L(h^fm + Uк(6)) - погрешность, обусловленная
m
введением сеток, U = (/ (bi — ai)2)1/2, L = max max \fY. (t, y)\.
i=1
Пусть для заданных сеток Лт и методом перебора найдено численное решение задачи (4) Ys(t) = (Yi(t), ■ ■ ■ ,Ym(t)), Yв(Д) = Va), при котором достигается S^, и ps(t,Ys(t)) = (p1(t,Ys(t)), ■ ■ ■ ,psm(t,Ys(t))) и определены tj — точки переключения скорости изменения переменной Ys(t), i = 1, ■■■, Nj, j = 1, ■ ■ ■, m.
Для вычисления £з(к(6), h) из теоремы 1 необходимо определить промежутки (ci,di), i = 1, ■ ■ ■, N1, на которых \ps(t,Ys(t))\ <6, j = 1, ■ ■ ■, m. Для нахождения промежутков (ci,di), i = 1, ■■■, N1 нам нужно задать 6 и по нему их определять. Естественно считать, что 6 ^ 60, где
60 = max max\pmax(t,Ymax(t)) — ps(t,Ys (t))\, гелт J J
pmax(t,Ymax(t)) = (pmax(t,Ymax(t)),■■■,pmax(t,Ymax(t)) - решение задачи (4), отвечающее оптимальной паре (vmax,Ymax(t)), при которой достигается MA.
Пусть множество A^ = {t £ Лт : 3j, 1 ^ j ^ m, такое, что \pSj(t,Ys(t))\ ^ 6}, где 6 > 0.
На практике функцию к(6) из условия 1 достаточно определить на отрезке \60, 61], где 0 < 60 ^ 61. Поэтому сформулируем следующее условие, которое удобно проверять в процессе вычисления.
Условие 2. 3 61 > 60, V 6 £ [60,61 ], V j (j = l, ■■■,№.), для Ys(t) и заданного v > 0 должно выполняться одно из условий при t £ A^: либо p = —df (t,Ys(t))/dYj > v, либо p = —df (t,Ys(t))/dYj < —v .
Nij
tij — 26/v, tij
j=1 i=1
ства A$v.
Теорема 2. Пусть f £ F, и выполнено условие 2. Тогда в качестве сужения функции к(6) из условия 1 на отрезок [60, 61] можно взять функцию
к(6) = n(As,v), 6 £ [60,61],
в частности, для одномерной задачи (при m = l)
к(6) = 2N16/v.
Доказательство. Рассмотрим в плоскости 0tpj соответствующую j-ю составляющую ps(t,Ys(t)) (c N1j точками переключений) решения ps(t,Ys(t)). Пусть выполняется условие 2, т. е. траектория pSj(t,Ys(t)) пересекает ось времени 0t для t £ A| со скоростью, имеющей постоянный знак и по модулю больше или равной vj. Оценим для этой составляющей из условия 1 промежутки (cij,dij) и Kj(6). Для
Обзначим Ag,v = — 26/v,t\j + 26/v], ^(Ag,v) - мера Лебега множе-
этого рассмотрим г-ю точку переключения Ь], Ь] е [е4],¿.1]], решения р8(Ь,^я(Ь)). Длина отрезков [Ь],¿4]] и [е4],Ь]] оценивается сверху следующим образом:
Ь» _ е.. < 54
") 4] 4]
V] V]' " ]'
¿4] — Ь] < —, Ь] — ец < , V] < шт п{р](Ь, 78(Ь))}.
Отсюда
¿4] — е] < 25]/V4], [е],¿4]] С [Ь] — 254]/V],Ь] + 254]/V]].
при условиях |78(Д) — 7тах(Д)| < Н, т](Ь) = 7тах(Ь) (т. е. она состоит из погреш-
С 1^4] ^ 4] 4] ^г] иг]
Тогда К] (5) для ]-той координаты
25 5
12 2] < 2^ ] = К] (5)
4=1 "4]
причем 5] ^ шах{54]}, min{vг]} ^ V] > 0, 1 < г < N1].
44
Объединяя все интервалы [Ь] — 254]/V4] ,Ь] + 254]/V4] ] для всех ] координат, получаем множество А$}У.
Таким образом, если 5 ^ шах{5]}, V < шin{v]}, 1 < ] < т, то можно взять
]]
к(5) = /л(А§у). Для одномерного случая получается более простое выражение к(5) = 2N15/v. Теорема доказана.
Для использования результатов теоремы 2 необходимо численно проверять условие 2.
Оценим 51 . Напомним, решается задача от момента времени Ь = Д к Ь = 0. На первом промежутке [Ь1], Д] определяется величина 51] для решения р8(Ь) по формуле
5] = ^Лпл. |Р°таХ(Ь,7таХ(Ь)) — Р](Ь,78(Ь))| ,
,](Л) — „,тах(Д)| < Н (Ь) = А,тах(
ности выбора начального условия и погрешности выбранного метода интегрирования в задаче (4)). При прохождении р](Ь) через 51]-окрестность оси р] = 0 в плоскости 0Ьр] проверяется выполнение условия 2. Если оно выполняется, то определяем оценочный промежуток [е^], равный [Ь] — 251]/V],Ь] + 251]/V]], для которого выполняется [еу,¿1]] С [е1], ¿1]]. Далее находим, насколько максимально разойдутся численное 7](Ь) и оптимальное 7]"ах(Ь) при прохождении через отрезок [е], ¿1] ], учитывая, что на данном промежутке скорости будут максимально отличаться, и в момент Ь = е\] разница достигнет Н1 = Н + (а] — Ь])(Л{] — е\]) или Н1 = Н — (а] — Ь] )^1] — е1]).
На следующем промежутке [Ь2],е]] ищется величина 52], исходя из того, что 1^|(е1])—7]"ах(е1])| < Н1, и учитывается погрешность интегрирования. Проверяется условие 2. В случае его выполнения определяем оценочный отрезок [е2] ,¿2] ] = = [Ь2] — 252] /v2] ,Ь2] + 252] /v2] ] , содержащий [е2] , ¿2] ], и затем находим, насколько разойдутся решения задачи (2) при прохождении через [е2],¿2]], учитывая, что на нем скорости будут максимально отличаться, разница при Ь = е2] достигнет Н2, равное Н1 — (а] — Ь]^¿2] — е2]) или Н + (а] — Ь])(¿2] — е2]).
На промежутке [Ь],е'?_1]] определяем 54], исходя из того, что |т8(е4-1]) —
— 7]"аХ(е|-1])| < Н4-1, где ^-1 = 2 + (а] — Ь]^(¿|_1] — е]-]) или ^-1 = 2 — (а] —
— Ь] 1] — е4-1]), и плюс погрешность интегрирования. Проверяется условие 2.
При его выполнении находится промежуток [е],¿4]] = [Ь] — 254]/V4],Ь] + 254]/V4]]. Потом определяем, насколько могут разойтись решения задачи (2) при прохождении через [е] ,¿4] ].
Таким образом определяем все Sij до момента времени t = 0 для всех координат j = 1,..., m. Положим ¿i = max max {Sij}. Очевидно, Si ^ So.
i^j^m i^i^Nij
Теперь рассмотрим предложенный метод вычисления предела максимального среднего на конкретном примере. Пусть f(t,j(t)) = sinY(t). Требуется оценить предел максимального среднего
Ms
1 rA
lim sup— sinY(t) dt, Y G [^1,^2], y(0) = y0.
A^w y Д Jo
Для решения строим максимальное среднее (3)
MA
-^sin
1
sup sup —
yoG[0,2^] Y
Д
и задачу Коши (4) на [0, Д]
p = — cos Y(t), ш1, при p < 0, Ш2, при p ^ 0,
Yj
sin y (t) dt
o
р(Д) = 0,
Y(Д) = У A.
(5)
(6)
(7)
В среде Delphi разработана программа по вычислению значения S^n, реализующая предложенный численный метод решения задач (6) и (7). Приведем некоторые результаты вычисления для w2] = [0, 5, 2], T = 2п: h = T/104, А = 2T, т = А/(2 • 105), уд = 464 • 2п/103, £ДП = 0, 543076; h = T/103, А = 20T, т = А/(2 х х 106), уд = 440 • 2п/103, 5ДП = 0,428766; h = T/103, А = 200T, т = А/(2 • 107), уд = 433 • 2п/103, Бдп = 0, 4168015. Для сравнения, значение Ms;n, вычисленное итерационным методом [5], равно 0,4151 с точностью 10~4.
Приведем оценку погрешности £3 из теоремы 1 (для простоты только для случая А = 2T), полученной по указанной схеме. Вычисления дают следующий результат: £3(к(6), h) = 0, 019983, где (по теореме 2) 6 = 0,0032, N = 2, v = 0, 99211, к(6) = 0, 0128.
A
Литература
[1] Лепилов А.Н. Численный метод вычисления пределов максимальных средних для периодических функций // Вестник СамГУ. 2011. № 8. С. 45-49.
[2] Филатов О.П. Вычисление пределов максимальных средних для периодических функций // Вестник СамГУ. 2011. № 2. С. 75-79.
[3] Алексеев В.М., Тихомиров В.М., Фомин С.В. Оптимальное управление. М.: Физматлит, 2007. 408 с.
[4] Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989. 432 с.
[5] Кайракбаев А.К., Филатов О.П. Итерационный метод вычисления пределов максимальных средних // Журнал вычислительной математики и математической физики. 1998. Т. 38. № 10. C. 1661-1664.
Поступила в редакцию 18/XI/2013;
в окончательном варианте — 18/XI/2013.
ON NUMERICAL ESTIMATE OF LIMITS OF MAXIMAL MEAN FOR PERIODIC FUNCTIONS
© 2013 A.N. Lepilov2
The scheme of definition of an error of numerical method of calculation of limit of maximal mean for periodic function is offered. The example of calculation of limit of maximal mean is considered.
Key words: limit of maximal mean, differential inclusion, periodic function.
Paper received 18/XI/2013. Paper accepted 18/XI/2013.
2Lepilov Alexander Nikolaevich (lepilov_aleksand@mail.ru), the Dept. of Mathematics and Mechanics, Samara State University, Samara, 443011, Russian Federation.