Задача граничного управления для системы уравнений гиперболического типа Текст научной статьи по специальности «Математика»

Поскольку

2ni Jr t - z 2ni Jj t - z ^ 2n 1 j j=1

f (xt) dxtdyt t\t\P t - z

- ВД,

причем первое и второе слагаемые правой части имеют в точке 0 особенности порядка p (для первого слагаемого этот результат можно найти в [1,2], а для второго это нетрудно получить с помощью известных интегральных неравенств, см., напр., [4]), то интеграл типа Коши (3) имеет в этой точке особенность порядка не ниже p+5. Легко убедиться, что этот порядок может превосходить единицу, и в этом случае задача о скачке (1) не имеет решений в классе функций, удовлетворяющих условию (2). Тем самым доказана

Теорема 1. Существует негладкая дуга Г с началом в точке 0 и концом в точке 1 и заданная на ней функция g(t) с особенностью порядка p в начальной точке кривой, удовлетворяющая условию Гёльдера в остальных её точках, для которых задача о скачке (1) не имеет решений в классе функций, удовлетворяющих условию (2).

В связи с этим возникает потребность в достаточных условиях разрешимости этой задачи. Сейчас мы приведем одно такое условие. Опишем сначала класс кривых, к которому оно применимо.

Простую жорданову дугу Г с началом в точке 0 и концом в точке 1 мы относим к классу Z0, если существует гладкая дуга Г' с теми же началом и концом такая, что симметрическая разность ГДГ' представляет собой счетное семейство замкнутых кусочно-гладких кривых, ограничивающих попарно не пересекающиеся области с единственной точкой сгущения в начале координат.

Кроме того, нам понадобится так называемая размерность Минковского (она же box dimension) кривой Г. Её определение можно найти, например, в [5].

Теорема 2. Пусть дуга Г имеет размерность Минковского d < 2 и принадлежит классу Z0, а заданная на ней функция g(t) представима в виде g(t) = \t\-pf (t), где f e Hv(Г), p < 1 и v > d/2. Тогда задача о скачке (1) имеет решение в классе функций, удовлетворяющих условию (2).

Доказательство. Доказательство в целом аналогично доказательству разрешимости задачи о скачке на замкнутой неспрямляемой кривой в работе [3]. □

Библиографический список

1. Гахов Ф. Д. Краевые задачи. М. : Наука, 1977. 640 с. [Gakhov F. Boundary value problems. Oxford : Permagon Press, 1966.]

2. Мусхелишвили Н. И. Сингулярные интегральные уравнения. М. : Наука, 1962. 600 с. [Muskheli-shvili N. I. Singular Integral Equations / ed. by J.R. M. Radok. Leyden : Noordhoff Intern. Publish., 1977.]

3. Кац Б. А. Задача Римана на замкнутой жордановой

кривой // Изв. вузов. Математика. 1983. № 4. С. 6880. [Kats B. A. The Riemann problem on a closed Jordan curve // Soviet Math. (Izv. VUZ. Matematika). 1983. Vol. 27, № 4. P. 83-98.]

4. Векуа И. Н. Обобщенные аналитические фукции. М. : Наука, 1988. 509 с. [Vekua I. N. Generalized Analytic Functions. Oxford : Pergamon Press, 1962.]

5. Федер Е. Фракталы. М. : Мир, 1991. 282 c. [Feder J. Fractals. New York : Plenum Press, 1988.]

R

j

УДК 517.956.3

ЗАДАЧА ГРАНИЧНОГО УПРАВЛЕНИЯ ДЛЯ СИСТЕМЫ УРАВНЕНИЙ ГИПЕРБОЛИЧЕСКОГО ТИПА

Е. А. Козлова

Самарский государственный технический университет E-mail: leni2006@mail.ru

Рассмотрена задача граничного управления для системы уравнений гиперболического типа. С помощью метода Римана построены управляющие функции, переводящие объект, описываемый системой, из заданного начального состояния в финальное.

Ключевые слова: граничное управление, система уравнений гиперболического типа, метод Римана.

Boundary Control Problem for the Hyperbolic System E. A. Kozlova

A boundary control problem for the hyperbolic system was considered. The control functions transferring the object described by this system from the given initial state to the final state were constructed using the Riemann method.

Keywords: boundary control, hyperbolic system, Riemann method.

© Козлова Е. А., 2013

51

ВВЕДЕНИЕ

Одной из актуальных задач современной теории управления является задача управления поведением объектов, изменение состояния которых описывается с помощью уравнений с частными производными. Цель управления состоит в том, чтобы перевести изучаемый объект из одного известного состояния в другое, влияя на некоторые его параметры. Впервые подобные задачи были сформулированы в работах А. Г. Бутковского [1] и J.-L. Lions [2].

В качестве управляющей функции может быть использована правая часть уравнения или системы уравнений (см., например, [3]). Далее будем рассматривать граничное управление, т.е. управление посредством граничных условий. В последние годы появились работы В. А. Ильина и Е. И. Моисеева, посвященные исследованию задач граничного управления для уравнения колебаний струны (в частности, [4-7]), в которых были получены в явном виде управляющие функции, переводящие струну из заданного начального состояния в заданное финальное состояние за определенное время. При этом рассматривались различные типы граничных управлений. Были исследованы задачи для уравнений колебаний неоднородной струны [8], радиально-симметричной мембраны [9], телеграфного уравнения [10].

Задачи, поставленные для уравнения колебаний струны, были обобщены на случай системы уравнений А. А. Андреевым и С. В. Лексиной в [11,12].

Настоящая работа посвящена задаче граничного управления для системы уравнений с частными производными вида

utt - Auxx + Cu = 0, (1)

где u(x, t) = (ui(x, t), u2(x, t),..., un(x,t))T — вектор-функция, A, C — постоянные квадратные матрицы порядка n, для которых верно AC = CA. Рассмотрим случай, когда матрица A имеет различные положительные собственные значения А2,..., А^ и выполняется условие А2 > ... > А^.

Пусть в прямоугольнике Q = [0,1] х [0, T] при T > 1/Ап задана система (1) и выполняются следующие начальные условия:

u(x, 0) = <°(x), ut(x, 0) = (x), 0 < x < 1, (2)

и финальные условия:

u(x,T) = <(x), ut(x, T) = ^1(x), 0 < x < 1. (3)

Найти граничные управления:

^(t) = u(0, t), v(t) = u(1, t), 0 < t < T.

Здесь <°(x), -0°(x), < (x), ^1(x), ^(t), v(t) — вектор-функции размерности n, <k(x) e C2[0,1], -k (x) e C 1[0,1], (i = 0,1); Ц,к(t),Vk(t) e C[0,T], k = î~n.

РЕШЕНИЕ ЗАДАЧИ УПРАВЛЕНИЯ

Известно [13], что существует невырожденная матрица S такая, что

S-1AS = Ja ,

где JA = diag {А2,..., А^} — жорданова нормальная форма матрицы A. Матрица перехода S составляется из собственных векторов матрицы A.

Матрицы A и C коммутируют, собственные значения A различны, поэтому C является простой и имеет тот же набор собственных векторов, что и A [13]. Таким образом, обе рассматриваемые матрицы приводятся к жордановой нормальной форме с помощью одного преобразования подобия: S-1CS = Je. Отметим, что Je = diag {с1,..., сп}.

Выполним замену w(x,t) = S-1 u(x,t), которая приведет систему (1) к виду

wtt - Jawxx + Jew = 0, (4)

начальные и финальные условия сведутся к

w(x, 0) = S-1 <°(x) = <°(x), wt(x, 0) = S-1 -°(x)= (x), 0 < x < 1, (5)

Е. Л Козлова. Задача граничного управления ррп системы уравнений гиперболического типа и

ш(х,Т) = Б-1ф1(х) = ф1 (х), Шг(х,Т) = Б-1ф1(х) = ф\х), 0 < х < I, (6)

соответственно, а граничные управления примут вид

ф(Ь) = ш(0, ь) = б-1^(г), ф(ь) = ш(1, г) = Б-1и(г), о < г < т.

Из (4)-(6) для к-й компоненты вектор-функции ш(х,Ь) получаем задачу управления в области Q, состоящую из уравнения

(шк)и - Х2к (шк)хх + окШк = 0, (7)

начальных условий:

Шк(х, 0) = фк(х), (шШк)г(х, 0)= фк(х), 0 < х < I, (8)

и финальных условий:

Шк(х, Т) = фк(х), (Шк)г(х, Т) = фк(х), 0 < х < I. (9)

Уравнение (7) является уравнением гиперболического типа второго порядка и имеет два семейства характеристик х + ХкЬ = С1 и х — ХкЬ = С2 [14].

Характеристики х — \кЬ = I — \кТ и х + ХкЬ = ХкТ разбивают прямоугольник Q на четыре области:

Ак = \I — Хк(Т — Ь) < х < Хк(Т — Ь), 0 < х < I, 0 < Ь < Т —

I 2Хк >

А2к = {Т + ^ < Ь < Т — х, 0 < х < ¡/2\, Хк Хк

Ак = {Хк(Т — Ь) < х < I — Хк(Т — Ь),Т — — < Ь < Т},

А4к = {Т — х < Ь < Т + 1/2 < х < А.

Хк Хк

Для построения решения задачи (7)-(9) расширим промежуток задания функций фк(х), Ф0(х) до отрезка [I — ХкТ, ХкТ] и положим Ак = {I — Хк(Т — Ь) < х < Хк(Т — Ь), 0 < Ь < Т — 1/(2Хк)}. Теперь в области Ак функция Шк(х,Ь), как решение задачи Коши, определена начальными условиями (8):

к к х+\к г

Шк(хЛ)= фк(х + Хк+ фк(х - Хкг) _ кР1 ^ ((х — г)2 — Хк(2)) фк (г) ¿г+

х—\к г

х+\к г

+¿7 0 ((х — г)2 — Хк Ь2)) Фк (г) (Ю)

а в области А к — финальными условиями (9):

= ф к (х + Хк (Т — Ь)) + ф к (х — Хк (Т — Ь)) ш к (х,Ь) — 2

х+\к (т-г)

Ск(Т — Ь) [ кР (2; ((х — г)2 — М

4Хк 3

х—\к (т-г) х+\к (т-г)

щ ((х — г )2 — Хк (Т — Ь)2)) ф к (г) йг—

^ / оР1 ((х — г)2 — Хк(Т — Ь)2)) (г) йг. (11)

х—\к (т-г)

Здесь оР1 (а; и) — вырожденная гипергеометрическая функция [15]. Решения построены методом Римана [14].

Поскольку точка (1/2, Т — 1/(2Хк)) принадлежит одновременно Ак и Ак, продолжения 'фк(х), ф0(х) необходимо выбирать таким образом, чтобы решения двух задач Коши в данной точке совпадали.

Введем константы аО° = / — АОТ, ак° = 0, „О1 = АОТ, ак1 = / и функции

2х — а к

гк(1)= й(2» — + й(а°) — I 0^ ц(г — „к,){г — 2» + „о)) й(г) ^

2 а;> 2

2ж —аЗ

/ 0^1 (г — „О)(г — 2» + „О)) (г) ¿г,

кк+т

» + т + А Ск(» — т — /) [ /, /Л , /г + /

^ (», т) = 'О ( ^^^ + ^ 4а2 " у 0Ъ ^2; (г — » — т)(т — » + /^ г^¿г,

0

Обозначим 6о(»,() = 0^ (^4^ (» — Ак(Т — '))(» + Ао(Т — () — .

,2 „ д 4 2

В областях АО и АО получаем две задачи с данными на характеристиках, решения которых имеют следующий вид:

(», () = Vю(», —Ао(Т — ()) + V!0(», Ак(Т — () — /) — 6ог00(//2) для Д2 (12)

и

^(», () = ^(», Ак(Т — () — /) + ^ (», —Ак(Т — ()) — 6кг11(1/2) для АО. (13)

Тогда управляющие функции Д 2(() = (0,() и /к(() = (/,() определяются из (10), (11)

Лк ^

0ГМ 2; 7Т2 — А к'

Дк(() = /(А ^ ^^ — «1 I 0^ (2; 4Ак (г2 — Ак'2^ й(2) +

— Лк t

Лк t

/ 0Ъ (1; 4Ак (г2 — Ак(2)) </0(г) ¿г, (14)

—Лк t

0 0 1+Лfct

.!(')= + А2') + — Ак') — ^ I 0 у, (2;4Ак ((, — 2)2 — Ак'2 ^ 0О(2) ^

1 — Лк t

1+Лк t

+2Ао / 01; 4АО ((/ — г)2 — Ак(2)) <(г) ¿г (15)

при 0 < ' < Т — //А о и из (12), (13)

До(') = 'Г () + г!0 (М^) — 6.(0,')гГ (2) +

—Лк(Т—1)

+ Со(А °(Г- () — в ! 0ъ (2; 4Ао(г + А к(Т — ())(/ — А к(Т — ())) ^ (¿г—

0

Лк (к—1

— Ск(А^;4лТо — О) ] 0(2;4Ао(г — Ао(Т — () + ,)(Г — г»' (¿г, (16)

0

/О(') = + г!^ — 62 —

Лк(Т—1)

—Ск(Ак(Т — () — г^ / _С1_(г — А (Т —'))(, (Т — () — /Л г0^. ¿г—

4А2

(2; 4Лo(г — А к(Т — ())(Ак(Т — () — /^ ¿г-

0

Е. Л Козлова. Задача граничного управления ррп системы уравнений гиперболического типа

ck(Хк(T - t))

l-Xk (T-t)

4Ak

0F1 (2;^ (z + X к (T - t) - 1)(t - T)) rf( ] dz

2

(17)

при T - l/\k < t < T.

Пусть матрица перехода S = (sj)nj=i • Тогда

Uk (x, t) = (ski... skn )(wi(x, t)... Wn(x, t))T, Mk(t) = (ski . . . Skn)(Ai (t) . . . Дn(t))T, Vk(t) = (ski . . . Skn)(Vi (t) . . . ¿>n(t))J

(18)

Формулы (14)-(18) позволяют определить все компоненты искомых вектор-функций и V(£). В рассмотренном случае Т > 1/Хп задача (1)-(3) не имеет единственного решения, т.е. функции и V(£) определяются неоднозначно и зависят от выбора продолжений начальных условий (2) на отрезок [I — Т, ХкТ].

Для значений п = 1, Х1 = 1, Т = I полученное решение согласуется с результатами В. А. Ильина и Е. И. Моисеева.

Библиографический список

1. Бутковский А. Г. Теория оптимального управления системами с распределенными параметрами. М. : Наука, 1965. 476 с. [Butkovskiy А. G. Distributed Control Systems. New York : American Elsevier Pub. Co., 1969. 446 p.]

2. Lions J.-L. Controle optimal de systemes gouvernes par des equalions aux derivees partielles. Paris : Dunod Gauthier-Villars, 1968. 426 p.

3. Агошков В. И. Методы оптимального управления и сопряженных уравнений в задачах математической физики. М. : Ин-т выч. мат. РАН, 2003. 256 с. [Agoshkov V. I. Optimal control methods and the method of adjoint equations in problems of mathematical physics. Moscow : Inst. Vych. Mat. RAN, 2003. 256 p.]

4. Ильин В. А. Граничное управление процессом колебаний на двух концах в терминах обобщенного решения волнового уравнения с конечной энергией // Дифференциальные уравнения. 2000. Т. 36, № 11. С. 1513-1528. [Il'in V. A. Two-endpoint boundary control of vibrations described by a finite-energy generalized solution of the wave equation // Differ. Equ. 2000. Vol. 36, № 11. P. 1659-1675.]

5. Ильин В. А. Граничное управление процессом колебаний струны на двух концах при условии существования конечной энергии // Докл. АН. 2001. Т. 376, № 3. С. 295-299. [Il'in V. A. Boundary control of the string oscillation process at two ends under conditions of the existence of a finite energy // Dokl. Math. 2001. Vol. 63, № 1. P. 38-41.]

6. Ильин В. А., Моисеев Е. И. Оптимизация за произвольный достаточно большой промежуток времени граничного управления колебаниями струны упругой силой // Дифференциальные уравнения. 2006. Т. 42, № 12. С. 1699-1711. [Il'in V. A., Moiseev E. I. Optimization of the boundary control of string vibrations by an elastic force on an arbitrary sufficiently large time interval // Differ. Equ. 2006. Vol. 42, № 12. P. 17751786.]

7. Ильин В. А., Моисеев Е. И. Оптимизация граничных управлений смещением на двух концах струны

за произвольный достаточно большой промежуток времени // Докл. АН. 2007. Т. 417, № 2. С. 160-166. [Il'in V. A., Moiseev E. I. Optimization of boundary controls by displacements at two ends of a string during an arbitrary sufficiently large time interval // Dokl. Math.

2007. Vol. 76, № 3. P. 828-834.]

8. Боровских А. В. Формулы граничного управления неоднородной струной. I // Дифференциальные уравнения. 2007. Т. 43, № 1. С. 64-89. [Borovskikh A. V. Formulas of boundary control of an inhomogeneous string. I // Differ. Equ. 2007. Vol. 43, № 1. P. 69-95.]

9. Ильин В. А., Моисеев Е. И. Граничное управление радиально-симметричными колебаниями круглой мембраны // Докл. АН 2003. Т. 393, № 6. С. 730-734. [Il'in V. A., Moiseev E. I. A boundary control of radially symmetric oscillations of a round membrane // Dokl. Math. 2003. Vol. 68, № 3. P. 421-425.]

10. Ильин В. А., Моисеев Е. И. Граничное управление на двух концах процессом, описываемым телеграфным уравнением // Докл. АН. 2004. Т. 394, № 2. С. 154158. [Il'in V. A., Moiseev E. I. A boundary control at two ends by a process described by the telegraph equation // Dokl. Math. 2004. Vol. 69, № 1. P. 33-37.]

11. Андреев А. А., Лексина С. В. Задача граничного управления для системы волновых уравнений // Вестн. Самар. гос. техн. ун-та. Сер. Физ.-мат. науки.

2008. № 1(16) С. 5-10. [Andreev A. A., Leksina S. V. The boundary control problem for the system of wave equations // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki. 2008. № 1(16). P. 5-10.]

12. Андреев А. А., Лексина С. В. Задача граничного управления в условиях первой краевой задачи для системы гиперболического типа второго порядка // Дифференциальные уравнения 2011. Т. 47, № 6. С. 843-849. [Andreev A. A., Leksina S. V. Boundary control problem for the first boundary value problem for a second-order system of hyperbolic type // Differ. Equ. 2011. Vol. 47, № 6. P. 848-854.]

13. Lancaster P. Theory of matrices. New York : Acad. Press, 1969. 316 p.

0

14. Бицадзе А. В. Некоторые классы уравнений в частных производных. М. : Наука, 1981. 448 с. [Bitsad-ze A. V. Some Classes of Partial Differential Equations. New York : Gordon and Breach Science Publishers, 1988.]

15. Erdélyi A., Magnus W., Oberhettinger F., Trico-mi F. G. Higher transcendental functions. Vol. I / ed. H. Bateman. New York; Toronto; London : McGraw-Hill Book Co, Inc., 1953. 302 p.

УДК 517.518.82

ПРИБЛИЖЕНИЕ ФУНКЦИИ БОЛЬЦАНО МНОГОЧЛЕНАМИ БЕРНШТЕЙНА

И. А. Козлова

Калужский государственный университет им. К.Э.Циолковского E-mail: irena1983.83@mail.ru

В настоящей работе рассматривается функция Больцано f (x), которая является непрерывной и недифференцируемой. Данная функция определяется как предел последовательности ломаных и для ее построения используются вспомогательные функции, представляющие собой ломаные. В работе получена оценка модуля непрерывности функции Больцано. Из полученной оценки следует, что данная функция принадлежит классу Липшица порядка 1/2 с константой M = 6, т.е. f (x) <G 6 Lip 1/2. Для функции Больцано при a = 1 и h = 1 построена последовательность многочленов Бернштейна и получена оценка погрешности приближения функции Больцано многочленами Бернштейна.

Ключевые слова: функция Больцано, модуль непрерывности, многочлены Бернштейна, оценка погрешности приближения.

Approximation of Boltsano Function by Means of Bernstein Polynomials

I. A. Kozlova

In given work is considered Boltsano function f (x), which can be represented in rows. Boltsano function is continuous and not differentiable. It is received the estimation of the module of continuity of Boltsano function. From the estimation of the module of continuity follows that function f (x) belongs to the Lipschitz class Lip 1/2 with the constant 6, i.e. f(x) e 6Lipl/2. For the Boltsano function for a = l and h = l it is presented the sequence of Bernstein polynomials and it is proved the estimation of the error of approximation for Boltsano function by means of Bernstein polynomials.

Key words: Boltsano function, module of continuity, Bernstein polynomials, estimation the error of approximation.

1. МОДУЛЬ НЕПРЕРЫВНОСТИ ФУНКЦИИ БОЛЬЦАНО

Функция Больцано f (x) строится следующим образом [1]. Определяются вспомогательные функции /o(x),/i(x), f2(x),..., fn(x),.... Графиком функции f0(x) является отрезок OA4, где 0(0;0), А4(а; h). Заменим отрезок OA4 ломаной OA1A2A3A4 так, что точки A1, A2, A3 имеют координаты

Ai (а/4, -h/2), A2 (а/2, 0), A3 (3а/4, h/2).

Функцию, имеющую график OA1A2A3A4, обозначим через fi(x). По функции f1 (x) строим функцию f2(x), аналогично заменив каждый из отрезков OA1, A1A2, A2A3, A3A4 соответствующими ломаными (рис. 1).

Повторим эту операцию n раз, придем к функции fn(x). Определим функцию Больцано f(x)

вначале в точках вида

x = ка/4п, 0 < k ^ 4n, (1)

n = 0,1, 2, 3,..., k — целое, полагая f (x) = fn (x). Нам остается еще определить f (x) для значений x, отличных от точек вида (1). Это можно сделать, полагая f (x) = lim f (t), где t — точки вида (1).

t^x

Теперь f (x) определена во всем отрезке [0,а] и является непрерывной функцией в этом промежутке.

В дальнейшем будем рассматривать функцию Больцано при h =1 и а = 1.

От геометрического способа задания функции

Рис. 1. Схема построения вспомогательных Больцано перейдем к аналитическому. Для эт°г° функций Больцано fc(ж) (пунктирная линия), используем разложение аргумента x в 4-ичную fo(ж) (отрезок OÄ4) и fi(x) (ломаная OÄiÄ4) дробь: b/c = (0, а1 а2 ... ап .. .)4. Разряды этой дроби

© Козлова И. А., 2013