SPLITTING METHOD OF THE STUDY OF MULTIPOINT BOUNDARY VALUE PROBLEMS ON THE HALF-LINE

Research article
Issue: № 7 (38), 2015
Published:
2015/08/15
PDF

Коняев Ю.А.1, Кленина Л.И.2

1доктор физико-математических наук, доцент, 2кандидат физико-математических наук, доктор педагогических наук, доцент,  Национальный исследовательский университет «Московский энергетический институт»

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

Аннотация

Изучены вопросы существования и единственности решения краевых задач для неавтономных однородных систем обыкновенных дифференциальных уравнений с T – периодической матрицей при наличии регулярных возмущений. Получены достаточные  условия устойчивости и асимптотической устойчивости решений таких задач.

Ключевые слова: системы с T - периодической матрицей, регулярные возмущения, устойчивость, асимптотика.

Konyaev Yu.A.1, Klenina  L.I.2

1PhD in Physics and mathematics, professor, 2PhD in Physics and mathematics, PhD of pedagogical Sciences, professor, National research University «Moscow power engineering Institute»

SPLITTING METHOD OF THE STUDY OF MULTIPOINT BOUNDARY VALUE PROBLEMS ON THE HALF-LINE

Abstract

Studied the existence and uniqueness of the solution of boundary value problems for non-Autonomous homogeneous systems of ordinary differential equations with T – periodic matrix in the presence of a regular perturbation.  We obtained sufficient conditions of stability and asymptotic stability of solutions of such problems.

Keywords: the systems with T - periodic matrix, the regular perturbations, the stability, the asymptotics.

  1. Существование и единственность решения краевой задачи на полуоси

Рассмотрим линейную систему

18-08-2015 10-02-12 при наличии краевых условий

18-08-2015 10-02-28

где Fj – некоторые постоянные матрицы, j=1,2,…,m.

При значении m=1 получаем известную начальную задачу, решение которой всегда существует и единственно. При 2 ≤ m n решение краевой задачи (1)–(2) не всегда существует.

При анализе многоточечных начально краевых задач вида (1)–(2) на полуоси следует выделить два момента.

  1. Существование единственного решения на любом конечном отрезке [0,t0] ⊂ R+, (t0>1)
  2. В случае существования такого решения на полуоси необходимо исследовать его на устойчивость (или асимптотическую устойчивость) при t→+∞.

Сформулируем [1] теорему об однозначной разрешимости и укажем конкретный алгоритм построения решения задачи (1)–(2) без использования традиционного аппарата функции Грина.

Теорема 1. Если для задачи (1)–(2) выполняется условие detF≠0, где 18-08-2015 10-07-37, а  Φ(t) – произвольная фундаментальная матрица соответствующей однородной (f (t) ≡ 0) системы (1), тогда краевая задача (1)–(2) однозначно разрешима на любом  [0,t0] (t0≥1)  и её решение может быть представлено в виде

18-08-2015 10-10-45

где

 18-08-2015 10-10-58

Φk(t) – квадратные матрицы (k = 1,2,…,n), содержащие столбцы фундаментальной матрицы Φ(t)  и при этом

 18-08-2015 10-12-18

Доказательство. Справедливость представления (3) следует из того, что оно является суммой общего решения Φ(t)C однородного уравнения и частного решения неоднородного уравнения (1) специального вида

18-08-2015 10-12-55

Это проверяется прямым дифференцированием последнего равенства.

В самом деле,

18-08-2015 10-13-13

При этом вектор C однозначно определяется краевыми условиями

18-08-2015 10-13-27

что и завершает доказательство теоремы 1.

Замечание 1. Теорема 1 является эффективной, т.е. решение явно выписывается формулой (3) только при известной фундаментальной матрице  Φ(t). В общем случае следует построить близкую к системе (1) другую систему =B(t)x, где разность A(t)B(t) на отрезке [0,t0] достаточно мала по некоторой норме. Тогда, если выполняется условие detF0≠0, где 18-08-2015 10-17-04, и Ψ(t)– произвольная фундаментальная матрица системы =B(t)x, то  решение задачи (1)–(2) может быть сведено к решению эквивалентного специального интегрального уравнения  x=Lx,        здесь

18-08-2015 10-18-27

и вектор C определяется по аналогии с формулой (4) как

18-08-2015 10-18-39

Оператор L, введённый в формуле (5),  при достаточной малости нормы ||A(t)B(t)||, будет сжимающим, что гарантирует однозначную разрешимость интегрального уравнения x=Lx и эквивалентной ему начально краевой  задаче (1)–(2).

  1. Анализ краевых задач на полуоси для однородных систем с T - периодической матрицей при наличии регулярных возмущений

Рассмотрим краевые задачи для однородных систем с  T – периодической матрицей A(t), представимой в виде A(t) = A0A(t)  где 18-08-2015 10-21-56 – среднее значение, а  δ – некоторый параметр.

Следуя методу расщепления [2, с.13] для произвольной квадратной n×n матрицы 18-08-2015 10-23-01, введём обозначения для её «диагональной» части 18-08-2015 10-23-15 и «бездиагональной» части  18-08-2015 10-23-25.

Теорема 2. Пусть дана однозначно разрешимая краевая задача с T периодической матрицей:

18-08-2015 10-23-50

где матричный ряд 18-08-2015 10-24-02 из T периодических матриц Ak(t) сходится абсолютно и равномерно по некоторой норме  при достаточно малых |ε|<1 и t ≥ 0.

Если матрица A0 имеет простой спектр 18-08-2015 10-25-44удовлетворяющий условиям 18-08-2015 10-25-58, тогда задача (6) при достаточно малых |ε|<1  с помощью невырожденной T периодической замены 18-08-2015 10-26-26;

18-08-2015 10-26-39

может быть приведена к краевой задаче с почти диагональной матрицей вида

18-08-2015 10-27-02

где постоянные диагональные матрицы  Λk и  T периодические матрицы Hk(t) однозначно определяются с помощью итерационного алгоритма (k = 1,2,…,N).

Доказательство. В условиях теоремы всегда существует невырожденная замена x=S0y, приводящая систему (6) к виду

18-08-2015 10-29-15

Последующее преобразование y=H(N)(t)z, невырожденное при |ε|<1, даёт нужный результат (7), если матрицы B(t), H(N)(t) и Q(t) удовлетворяют дифференциальному матричному уравнению

18-08-2015 10-32-13

Приравнивая в последнем уравнении коэффициенты при одинаковых степенях ε, получим однотипные дифференциальные матричные уравнения вида

18-08-2015 10-33-01

где   18-08-2015 10-34-11   при k=1,2,…,N.

Разбив уравнение (8) на «диагональную» и «бездиагональную» части, получим два уравнения

18-08-2015 10-34-56

Первое из этих уравнений имеет единственное периодическое решение 18-08-2015 10-35-12, если  18-08-2015 10-35-26  (k=1,2,…,N).

Второе матричное уравнение распадается на  (n2n) скалярных дифференциальных уравнений  вида: 18-08-2015 10-40-59(i, j=1,2,…,n; ij; k=1,2,…,N), имеющих в условиях теоремы 2 единственное  T – периодическое решение

18-08-2015 10-43-02

Оценка ||G(t)||≤C проверяется непосредственно. Теорема 2 доказана.

  1. Анализ устойчивости решений

Теорема 3.  Если в условиях теоремы 2 спектр 18-08-2015 10-44-13 матрицы Λ(N)(ε) в задаче (7)  удовлетворяет неравенствам 18-08-2015 10-45-07 (j=1,2,…,n; q=0,1,2,…,N),  тогда решение задачи (7)  на полуоси и эквивалентной ей задаче (6) асимптотически устойчиво.

Доказательство. Применяя дифференциальное неравенство для квадрата нормы решения задачи [2, с.14]  получаем оценку 18-08-2015 10-46-21, из которой следует асимптотическая устойчивость решения задачи (6), что и требовалось доказать.

Теорема 4. Пусть выполнены условия теоремы 2. Если спектр 18-08-2015 10-46-36 матрицы Λ(N)(ε) в задаче (7)  удовлетворяет неравенствам  18-08-2015 10-47-08 (j=1,2,…,n), то асимптотика  (ε→+0) решения краевой задачи (6) с  Т – периодической матрицей может быть представлена в виде: 18-08-2015 10-49-03, где вектор C однозначно определяется из краевых условий 18-08-2015 10-49-22, а матричные функции Λ(N)(ε) и H(N)(tj) однозначно определяются  методами теоремы 2.

Доказательство. Применяя дифференциальное неравенство Бернулли [2, с.15]: 18-08-2015 10-50-41, а  вектор функция  p(t) удовлетворяет вспомогательной задаче 18-08-2015 10-52-14, получаем оценку 18-08-2015 10-52-33, которая завершает доказательство теоремы 4 [2, с.15].

Замечание 2. Теорема  2 является аналогом метода усреднения для заданного класса задач.

  1. Анализ периодических систем в критических случаях

Рассмотрим систему:

18-08-2015 10-54-19

(где матричный ряд 18-08-2015 10-54-28  из непрерывных T периодических матриц Ak(t), (k=0,1,2,…), сходится абсолютно и равномерно при  |ε|<<1 и t ≥ 0).

Теорема 5.  Для однозначно разрешимой краевой задачи (9) при наличии T периодической  матрицы A0(t) произвольной жордановой структуры  при |ε|<<1 существует невырожденная  Т – периодическая замена

18-08-2015 10-56-57

приводящая задачу (9) к эквивалентной задаче с почти постоянной матрицей вида

18-08-2015 10-57-15

причем матрицы  Bk  и Hk(t) однозначно определяются с помощью итерационного алгоритма  (k=0,1,2,…,N).

Доказательство. Прямая подстановка (10) в (9) даёт нужный результат (11), если матрицы A(t), H(N)(t) и B(t)  удовлетворяют дифференциальному матричному уравнению

18-08-2015 11-00-42

Приравнивая в последнем уравнении  коэффициенты при одинаковых степенях ε, получим набор однотипных дифференциальных матричных уравнений вида 18-08-2015 11-01-22,

18-08-2015 11-01-38, которые имеют единственное  Т – периодическое решение при каждом k=1,…,N. В этом случае 18-08-2015 11-02-50.

 Оценка ||G(t,ε)||≤C0 проверяется простым вычислением.

Теорема 6.  Если в условиях теоремы  5 для однозначно разрешимой краевой задачи (9)  спектр 18-08-2015 11-04-20 матрицы B0 удовлетворяет условиям 18-08-2015 11-04-35  тогда решение задачи (9) и эквивалентной ей задачи (11) будет асимптотически устойчиво.

Доказательство.  Всегда существует невырожденная замена:

18-08-2015 11-04-52

С помощью этой замены задача (11) приводится к эквивалентной ей задаче вида

18-08-2015 11-05-33

Оценка |z(t,ε)|≤|z0|exp(-εσ1t)→0 при t→+∞, гарантирует асимптотическую устойчивость решений эквивалентных задач (12), (11) и (9), что и требовалось доказать.

Вывод. В данной статье с помощью нетрадиционного  подхода (т.е. без использования функции Грина) продемонстрирован метод расщепления для анализа краевых задач на полуоси.

Для систем с T –  периодической матрицей получены конструктивные достаточные критерии устойчивости (и асимптотической устойчивости) решений таких задач, что можно рассматривать как обобщение теоремы Ляпунова об асимптотической устойчивости по первому приближению для указанного класса задач.

Литература

  1. Коняев Ю.А. Асимптотические и аналитические методы решения некоторых классов прикладных модельных задач: Монография. – М.: изд-во Российского университета дружбы народов, – 160 с.
  2. Коняев Ю.А., Кленина Л.И. О некоторых классах многоточечных краевых задач для - периодических систем на полуоси / Современное состояние естественных и технических наук: Материалы XIX Международной научно-практической конференции (19.06.2015). – М.: Издательство «Спутник +», с. 10-16.

References

  1. Konyaev Yu. A. Asymptotic and analytical methods for solving some classes applied model problems: Monograph. – M.: publishing house of the Russian University of friendship of peoples, 2005. – 160 p.
  2. Konyaev Yu. A., Klenina I. On some classes of multipoint boundary value problems for - periodic systems on the half-line / The modern condition of natural and technical Sciences: The materials of the XIX International scientific-practical conference (19.06.2015). – M.: Publishing House "Sputnik +", 2015. p. 10-16.