Научная статья на тему 'Устойчивость решений дифференциальных уравнений в частных производных параболического типа'

Устойчивость решений дифференциальных уравнений в частных производных параболического типа Текст научной статьи по специальности «Математика»

CC BY
918
102
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УСТОЙЧИВОСТЬ / СИСТЕМЫ ПАРАБОЛИЧЕСКИХ УРАВНЕНИЙ / STABILITY / SYSTEMS OF PARABOLIC EQUATIONS

Аннотация научной статьи по математике, автор научной работы — Бойков Илья Владимирович, Паксялева Оксана Геннадьевна, Романова Людмила Дмитриевна

Получены достаточные условия устойчивости систем параболических уравнений с коэффициентами, зависящими от времени и пространственных координат.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Устойчивость решений дифференциальных уравнений в частных производных параболического типа»

УДК 517.19

И. В. Бойков, О. Г. Паксялева, Л. Д. Романова

УСТОЙЧИВОСТЬ РЕШЕНИЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ ПАРАБОЛИЧЕСКОГО ТИПА

Аннотация. Получены достаточные условия устойчивости систем параболических уравнений с коэффициентами, зависящими от времени и пространственных координат.

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

Abstract. Given the sufficient conditions of stability of parabolic systems. The coefficients of parabolic systems are depended on time and space coordinates. Keywords: stability, systems of parabolic equations.

Исследованию устойчивости решений различных классов уравнений в частных производных посвящена обширная литература, подробная библиография которой приведена в [1-3]. Основным аппаратом исследования устойчивости решений линейных и нелинейных уравнений в частных производных является применение преобразования Фурье и построение обобщенных функционалов Ляпунова.

В работе [4] был предложен метод построения критериев устойчивости решений систем линейных и нелинейных дифференциальных уравнений параболического типа, основанный на исследовании спектров и логарифмических норм специальным образом построенных семейств линейных операторов. Недостатком этого метода является то обстоятельство, что он применим только к уравнениям, определенным во всем пространстве пространственных переменных.

В данной работе предложен другой метод исследования устойчивости решений систем дифференциальных уравнений параболического типа. Этот метод применим к уравнениям, заданным в конечных областях. Он является распространением метода, описанного в [5], на системы уравнений в частных производных параболического типа.

Отметим, что предлагаемые ниже критерии устойчивости могут быть распространены и на некоторые другие системы дифференциальных уравнений с распределенными параметрами.

В данной работе рассматриваются системы линейных параболических уравнений, имеющие при любых начальных значениях из рассматриваемой области решения, определенные при всех значениях t е [tg,тс).

При исследовании устойчивости решений линейных параболических уравнений ограничимся рассмотрением уравнений

dui{t,х) _ац(t,х)д^1 (t,x) + ai2(t,x)Au2(t,x) + dt

+a13(t, x)u1(t, x) + a14(t, x)u2(t, x); du2(^,x) _a21 (t,x)Au1 (t,x) + a22(t,x)Au2(t,x) +

+023^, х)щ(V, х) + 024^, х)ы2^, х) (1)

в ограниченных областях й с гладкой границей Г = дQ. Здесь

и (V, х) = (и^, х), ^2(V, х)), х = (х1,х2).

Из проводимых ниже рассуждений следует, что полученные результаты распространяются на системы размерности большей двух.

Уравнение (1) имеет тривиальное решение при нулевых начальном и граничном условиях. Будем исследовать устойчивость тривиального решения при нулевом граничном условии

и (V, х)!хег =0 (2)

и при возмущении начального условия

и (¿о, х) = и1( х). (3)

Исследование устойчивости решения системы уравнений (1) при граничном значении (2) будем проводить в пространстве функций и (V, х) = (и1 (V, х), и2 (V, х)) с нормой

Iи (t) ||2= JJ(u—(t, x) + u—(t, x))dx.

Q

Умножим первое из уравнений системы (1) на щ(ґ,х), второе - на и2 (^, х) и проинтегрируем каждое из них по рассматриваемой области интегрирования. В результате получаем систему уравнений

——JJ u2 (t, x)dx = JJ a„(t, x)ui (t, x)Aui (t, x)dx +

Q Q

+JJ ai—(t, x)ui (t, x) Au— (t, x)dx +

Q

JJ °l3(t, x)u— (t, x)dx + JJ ai4(t, x)ui(t, x)u— (t, x)dx;

+

Q Q

i d JJu—(t, x)dx = JJ a—i(t , x)u— (t, x)Aui (t, x)dx -

Q

JJ a——(t , x)u— (t, x)Au— (t, x)dx +

2 dt

Q Q

+ l

Q

+JJ a—3 (t, x)ui (t, x)u— (t, x)dx + JJ a—4 (t, x)u— (t, x)dx. (4)

Q Q

Воспользуемся интегральной теоремой Остроградского - Гаусса

J J (grad u,grad v)dx + JJ vAudx = Jv—ds, (5)

Q Q Г

где п - единичная нормаль к контуру Г.

2i

Положим v(t, x) = uj(t, x)ajj(t, x).

Известно, что

grad (uj (t, x)ajj(t, x) ) = ajj(t, x)grad u(t, x) + u (t, x)grad ajj(t, x).

Так как на границе Г области й uj(t, x) = 0, то из равенства (5) следует, что

Ц ajj(t, x)uj (t, x)Auj (t, x)dx = -JJ (ajj(t, x)grad uj (t, x), grad uj (t, x) )dx -й й

- JJ (uj (t, x)grad an(t, x), grad uj (t, x)) dx. й

Следовательно, если функции ajj(t, x), aj2(t, x), a2j(t, x), a22(t, x) зависят только от t, то приходим к следующей системе:

"2~~ J J uj2 (t, x)dx = -aj j (t) J J (grad uj (t, x), grad uj (t, x) )dx -й й

-aj2 (t) J J (grad uj (t, x), grad u2 (t, x)) dx + й

+JJ aj3(t, x)uj2 (t, x)dx + JJ Oj4(t, x)uj(t, x)u2 (t, x) dx;

й й

-—J Ju- (t, x)dx = _°2і(^) J J (grad uj (t, x), grad u- (t, x) )dx -й й

—a-- (t) J J (grad u- (t, x), grad u- (t, x) - dx +

й

+

J J a-з (t, x)ui (t, x)u- (t, x)dx + JJ a-4 (t, x)u- (t, x)dx. (6)

й й

Исследуем устойчивость решения системы уравнений (6) в предположении, что функции ajj (t, x), i = 1,-, j = 1,-,3,4, не зависят от x. В результате приходим к следующей системе уравнений:

—JJui- (t, x)—x = —-aj j (t) J J (grad uj (t, x))- dx — й й

—-aj-(t) JJ (grad uj (t, x), grad u- (t, x) )dx + й

+-aj3 (t) J J u- (t, x)dx + -aj4 (t) J J uj (t, x)u- (t, x)dx;

й й

j JJu2 (t, x)dx = -2«2j(t )JJ (grad uj (t, x), grad u2(t, x) )x-

й й

-2a22 (t) J J (grad u2 (t, x))2 dx + 2a23 (t) J J uj (t, x)u2 (t, x)dx + й й

+2a24 (t) J J u2 (t, x)dx. (7)

й

Сложим первое и второе уравнения системы (7). В результате приходим к уравнению

j JJ (uj2 (t, x)+ uf (t, x))dx = f (t) =

й

= -2ajj(t) J J (grad uj (t, x), grad uj (t, x))dx -й

-2(aj2 (t) + a2j(t)) J J (grad uj (t, x),grad u2 (t, x))dx + й

+2( aj4 (t) + a23 (t)) J J (uj (t, x), u2 (t, x))dx + й

+2aj3(t )JJuj2(t, x)dx + 2a24(t) JJu2 (t, x)dx -й й

-2a22 (t) J J (grad u2 (t, x), grad u2 (t, x))dx. (8)

й

Преобразуем функцию f (t) следующим образом:

f (t) = -(aj2(t) + a2j(t))) (grad uj (t, x) + grad u2 (t, x) )2 dx +

й

+ (aj2(t) + a2j(t) - 2ajj(t) ) J J (grad uj (t, x) )2 dx +

й

+ (aj2(t) + a2j(t) - 2a22 (t)) (grad u2 (t, x) )2 dx +

й

+ (aj4(t) + a23(t)) J J ((t,x) + u2 (t, x))2 dx -й

-(aj4(t) + a23 (t) - 2aj3(t)) u2 (t, x)dx -

й

-(aj4(t) + a23 (t) - 2a24 (t)) u| (t, x)dx.

й

Предыдущую функцию можно представить в следующем виде: f (t ) = - (aj2(t) + a2j(t) ))(grad u j( t, x) + grad u2(t, x) )2 dx +

й

+ (а12(^) + а2\({)- 2ац(г)) 11 (grad щ (г, х) )2 —х +

а

+ ((*) + а21 (г) - 2а22 (г)) (grad и2 (г, х) )2 )х +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а

+ (4(0 + а2з(г) )) ( (г, х) + и2(г, х) )2 2х -а

-(а^(г) + а2з(г) - 2а1з(г) )) (и 2 (г, х) + и2 (г, х) | —х +

а

+ (2а24 (г) - 2а1з(г) ))и2 (г, х) —х.

а

Предположим, что выполнены следующие условия при всех значениях

г (г > г0):

012 (г) + а21 (г) > 0, а12 (г) + а21(г) - 2ап (г) < 0; а12 (г) + а21(г) - 2а22 (г) < 0, «14 (г) + а2з (г) < 0; ам(0 + а2з(г) - 2а^(0>0, а^О -а^(0 < 0. (9)

Нетрудно видеть, что условия (9) непротиворечивы.

Тогда при всех значениях г (г > ¿0)

/ (0 < Ж0, (10)

где

Л(г) = -(((0 + «2з (г) - 2«1з( г)) (и( (г, х) + и| (г, х)) —х.

а

Рассмотрим уравнение

—|||и(, х) + и22 (г, х) )х = а

= - (а14(0 + а2з ( г) - 2а1з( г))) (и2 (¿, х) + и2 ( г, х)) —х. (11)

а

Очевидно, что если при всех значениях г( г > ^) выполняется условие

-(ам(0 + а2з( 0 - 2а1з( г)) < 0, (12)

то решение уравнения (11) устойчиво при любом начальном приближении.

Замечание. Если, начиная с некоторого значения T, выполняется неравенство

-(а14(0 + а2з(г)-2а1з(г))<-а, а>0, (1з)

то решение уравнения (11) устойчиво в целом.

Так как при сделанных предположениях /(г) < /1(г), то из теоремы Чаплыгина [6, с. 5] следует, что из устойчивости (асимптотической устойчивости) решения уравнения (11) следует устойчивость (асимптотическая устойчивость ) решения уравнения (8).

Функцию / (г) можно также представить в следующем виде:

f (t) = -(2(0 + a2j(t) ))(grad uj(t, x) + grad u2(t, x) )2 dx +

й

+ (aj2(t) + a2j (t) - 2ajj(t))) ) (grad uj (t, x) )2 dx +

й

+ (aj2(t) + a2j(t) - 2a22 (t)) (grad u2 (t, x))2 dx +

й

+ (aj4(t) + a23(t)) J(uj(t,x) + u2(t,x))2 dx -й

-(aj4(t) + a23(t) - 2a24 (t)) J J ((t, x) + u^ (t, x) ) -

й

- (2a24 (t) - 2aj3(t))) ) uj2 (t, x) dx.

й

Тогда при выполнении условий

aj2 (t) + a2j(t) > 0, aj2 (t) + a2j(t) - 2an (t) < 0; aj2 (t) + a2j(t) - 2a22 (t) < 0, aj4 (t) + a23 (t) < 0;

aj4(t) + a23(t) - 2a24(t) > 0, a24(t) - a23(t) > 0 ^4)

выполняется неравенство f (t) < f2 (t), где

f2(t) = -(aj4(t) + a23(t) - 2a24 (t)) (uj (t, x) + u2 (t, x)) dx.

й

Теперь, повторяя приведенные выше рассуждения, убеждаемся, что при выполнении условий ^4) система уравнений (j) устойчива, а при выполнении условий ^4) и условия

-(aj4(t) + a23(t)-2a24(t))<-a, a>0, (j5)

она асимптотически устойчива.

Таким образом, доказаны следующие утверждения.

Теорема 1. Пусть выполнены условия (9) или (14). Тогда решение системы уравнений (1) устойчиво.

Теорема 2. Пусть выполнены условия (9) и (13) или (14) и (15). Тогда

решение системы уравнений (1) асимптотически устойчиво.

Список литературы

1. Эйдельман, С. Д. Параболические системы / С. Д. Эйдельман. - М. : Наука, 1964. - 444 с.

2. Сиразетдинов, Т. К. Устойчивость систем с распределенными параметрами / Т. К. Сиразетдинов. - М. : Наука, 1987. - 232 с.

3. Шес таков, А. А. Обобщенный прямой метод Ляпунова для систем с распределенными параметрами / А. А. Шестаков. - М. : Наука, 1990. - 320 с.

4. Бойков, И. В. Об определении областей устойчивости для некоторых классов нелинейных уравнений с распределенными параметрами / И. В. Бойков // Автоматика и телемеханика. - 2001. - № 1. - С. 40-49.

5. Пу, Т. Нелинейная экономическая динамика / Т. Пу. - Ижевск : Издательский дом «Удмурдский университет», 2000. - 200 с.

6. Мартынюк, А. А. Устойчивость движения. Метод интегральных неравенств / А. А. Мартынюк, В. Лакшмикантам, С. Лила. - Киев : Наукова думка, 1989. - 272 с.

Бойков Илья Владимирович

доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет

E-mail: math@pnzgu.ru

Паксялева Оксана Геннадьевна

аспирант, Пензенский государственный университет

E-mail: math@pnzgu.ru

Романова Людмила Дмитриевна доцент, кафедра высшей и прикладной математики, Пензенский государственный университет

E-mail: math@pnzgu.ru

Boykov Ilya Vladimirovich Doctor of physico-mathematical sciences, professor, head of sub-department of higher and applied mathematics,

Penza State University

Paksyaleva Oksana Gennadyevna Postgraduate student,

Penza State University

Romanova Lyudmila Dmitrievna Associate professor, sub-department of higher and applied mathematics, Penza State University

УДК 517.19 Бойков, И. В.

Устойчивость решений дифференциальных уравнений в частных производных параболического типа / И. В. Бойков, О. Г. Паксялева, Л. Д. Романова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2009. - № 4 (12). - С. 20-26.

i Надоели баннеры? Вы всегда можете отключить рекламу.