Телефон: +7 (383)-202-16-86

Статья опубликована в рамках: XI Международной научно-практической конференции «Естественные и математические науки в современном мире» (Россия, г. Новосибирск, 14 октября 2013 г.)

Наука: Математика

Секция: Дифференциальные уравнения, динамические системы и оптимальное управление

Скачать книгу(-и): Сборник статей конференции

Библиографическое описание:
Турсунов Д.А. АСИМПТОТИКА РЕШЕНИЯ СИНГУЛЯРНО ВОЗМУЩЕННОГО ОБЫКНОВЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ С ДВУМЯ ТОЧКАМИ ПОВОРОТА В КОМПЛЕКСНОЙ ПЛОСКОСТИ // Естественные и математические науки в современном мире: сб. ст. по матер. XI междунар. науч.-практ. конф. № 9-10(10). – Новосибирск: СибАК, 2013.
Проголосовать за статью
Дипломы участников
У данной статьи нет
дипломов
Статья опубликована в рамках:
 
Выходные данные сборника:


 


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


Турсунов  Дилмурат  Абдиллажанович


канд.  физ.-мат.  наук,  доцент  ОшГУ,  г.  Ош


E-maildosh2012@mail.ru


 

ASYMPTOTICS  OF  SOLUTIONS  OF  SINGULARLY  PERTURBED  ORDINARY  DIFFERENTIAL  EQUATIONS  WITH  TWO  TURNING  POINTS  IN  THE  COMPLEX  PLANE


Dilmurat  Tursunov


candidate  of  physics  and  mathematics  sciences,  associate  professor  of  Osh  State  University,  Osh


 


АННОТАЦИЯ


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


ABSTRACT


The  aim  of  this  paper  is  to  construct  a  uniform  asymptotic  solution  of  ordinary  differential  equations  with  a  simple  turning  point  in  the  complex  plane  when  a  condition  of  asymptotic  stability  is  violated.  The  method  of  stationary  phase  and  the  passes  method  are  used  to  evaluate  the  integrals.  The  uniform  asymptotic  solution  of  the  problem  was  constructed.  Asymptotic  of  the  solution  of  the  problem  essentially  depends  on  the  inhomogeneous  part  of  the  equation.


 


Ключевые  слова:  асимптотика;  точка  поворота;  условия  устойчивости;  стационарная  фаза;  точка  перевала;  малый  параметр.


Keywordsasymptotic;  turning  point;  stability  conditions;  stationary  phase;  passes  point;  small  parameter.


 


Работа  выполнена  при  поддержке  гранта  РФФИ  по  проекту  №  13-01-90903  мол_ин_нр


 


Теорию  запаздывания  потери  устойчивости  в  системах  общего  вида  построены  в  работах  А.И.  Нейштадта  [2].  Далее  появились  работы  [1],  [3]  и  др.  Численными  и  аналитическими  методами  явление  запаздывание  потери  устойчивости  изучалось  в  системах,  встречающихся  в  физике  лазеров,  химической  кинетике,  биофизике,  в  модифицированной  системе  Циглера  и  в  моделировании  безопасных  процессов  горения  с  максимальной  температуры.  Во  всех  этих  перечисленных  работах  невозмущенное  (предельное)  уравнение  имеет  тривиальное,  т.е.  гладкое  решение  в  рассматриваемой  области.  В  нашем  случае  решение  предельного  уравнения  имеет  особенность.


Рассмотрим  задачу


 


ex¢(t,e)=A(t)x(t,e)+f(t),                                  (1)


x(t0,e)=x0,                                             (2)


 

где:  A(t)  –  квадратная  матрица-функция  второго  порядка  с  элементами  аjk(t),  f(t)=colon{f1(t),f2(t)}, 


аjk(t),  fk(t)  –  аналитические  функции  в  D, 


x0=colon{x10,x20}  –  постоянный  вектор,  tÎDt=t1+it2

U1.  Пусть  матрица-функция  A(t)  имеет  комплексно-сопряженные  собственные  значения  l1,2(t)=sint±iacosta>1,  и  t0=  –p/2.

U2.  Пусть  f1(–p/2,  –a)¹0,  f1(p/2,  –a)¹0,  f2(–p/2,a)¹0,  f2(p/2,a)¹0.

В  рассматриваемом  случае  каждое  собственное  значение  имеет  по  две  простых  периодических  нулей  в  комплексной  плоскости.  Исследуемая  область  является  прямоугольником,  и  нули  собственных  значений  матрицы-функции  A(t)  находятся  на  вершинах  этого  прямоугольника.  Такой  случай  рассматривается  впервые.  Требуется  построить  асимптотику  решения  задачи  (1)—(2)  в  области  DÉ[–p/2,  p/2],  при  e®0.


Для  приведения  A(t)  к  диагональному  виду  выполняем  следующее  преобразование:  (t)A(t)B0(t)=D(t),  где 


 


В0(t)=D(t)=diag(l1(t),  l2(t)).


 


Пусть  в  области  D  выполняется  неравенство  detВ0(t)¹0. 

Задача  (1)—(2)  с  заменой  x(t,e)=B0(t)у(t,e)  принимает  вид:

 


eу'(t,e)=D(t)у(t,e)+eB(t)у(t,e)+h(t),                                   (3)


у(t0,  e)=у0,                                           (4)


 


где  .

Задачу  Коши  для  дифференциальных  уравнений  (3)-(4)  заменим  интегральным  уравнением:

 


  (5)


где  z(t,e)=y(t,e)e, 


В  работе  [4]  нами  доказана 


Теорема.  Если  для  интеграла 


 


,  (6)


 


в  некоторой  области  D  справедлива  оценка  ,  где  ,  то  для  решения  системы  интегральных  уравнений  (5)  справедлива  оценка  ||z(t,e)||£cd(e).

Значит,  достаточно  вычислить  асимптотику  интеграла  (6),  при  условиях  U1  и  U2.  Рассмотрим  собственные  значения


Re(l1(t1,0))=Re(l2(t1,0))=sint1;  Re(l1,2(t1,0))<0  –  устойчивый  интервал,  при  –p+2pk<t1<2pk;  Re(l1,2(t1,0))>0  –неустойчивый  интервал,  при  2pk<t1<p+2pk


Re(l1,2(t1,0))=0,  при  t1=pkkÎZ.


Отметим,  что  (p/2+pk,–a)  и  (p/2+pk,a),  где    kÎZ,  являются  нулями  собственных  значений  l1(t1,t2)  и  l2(t1,t2)  соответственно.


Рассмотрим  теперь  функции:  u1(t)=òl1(t)dtu2(t)=òl2(t)dt.


Если  t=t1+it2,  то  u1(t1,t2)=  –cost1(cht2+asht2)+isint1(sht2+acht2),  u2(t1,t2)=


–cost1(cht2asht2)+isint1(acht2–sht2).  Пусть  u11(t1,t2)=Re(u1(t1,t2)),  u21(t1,t2)=Re(u2(t1,t2)).  Область  D={tu11(t1,t2)£0,  u21(t1,t2)£0,  |t1|£p/2}  –  является  прямоугольником  с  вершинами  A(–p/2,–a),  B(p/2,–a),  B1(p/2,a)  и  A1(–p/2,a). 


Решение  предельной  (вырожденной)  системы:  (t)=  –A-1(t)f(t)  в  четырех  точках  A(–p/2,–a),  B(p/2,–a),  B1(p/2,a)  и  A1(–p/2,a)  имеет  особенность,  при  f(t)¹0.  Кроме  того  это  решение  не  удовлетворяет  начальному  условию  (t0x0.


Таким  образом,  мы  здесь  тоже  сталкиваемся  с  двумя  сингулярностью:  первое  это  наличие  пограничного  слоя,  а  второе  —  соответствующее  невозмущенное  (предельное)  уравнение  имеет  особенность  в  точках  A(–p/2,–a),  B(p/2,–a),  B1(p/2,a)  и  A1(–p/2,a).  Каждое  уравнение  системы  (1)  имеет  по  две  точки  поворота,  в  целом  система  имеет  четыре  точек  поворота.  Поэтому  рассматриваемую  задачу  можно  называть  бисингулярной.


Перейдем  к  оценке  интегралов


 


,


 

где:  L  —  пути  интегрирования,  соединяющие  точки  (–p/2,0)  и  tÎD.  Пути  интегрирования  L  симметричны  относительно  действительной  оси.  Область  D  разобьём  на  подобласти:


 


Н00={t:  0£t1+p/2£d,  d£t2+a£2a},  Н10={t:0£t1+p/2£degt2+a=deg},


Н01={t:–p/2+d£t1u11(t1,t2)£(elne)/2,  d£|(t1–p/2)+i(t2+a)|,  t2£a},


Н11={t:–((1–2g)elne)/2£t2+a£deg,  t1+p/2=deg},


Н12={t:  |t1|£p/2–deg,  u11(t1,t2)=((1–2g)elne)/2},


Н13={t:–((1–2g)elne)/2£t2+a£deg,  t1=p/2–deg},


Н14={t:–  deg£t1–p/2£((1–2g)elne)/2,  t2+a=deg},


Н15={tu11(t1,t2)=((1–2g)elne)/2,  deg£t2+a£2a},


Н20={t:0£t1+p/2£d,  0£t2+a£d},  Н21={t:  |t1|£p/2–d,  (elne)/2£u11(t1,t2)  £0},


Н22={t:–d£t1–p/2£0,  (elne)/2£u11(t1,t2)£0},  0<d<<1  —  достаточное  малое  число,  0£g<1/2,  Н0=Н00ÈН01Н1=Н2=D=Н0ÈН1ÈН2.


Лемма  1.  Если  tÎH00,  то  для  интеграла  J(t,e)  справедлива  оценка


 


|J(t,e)|£ce.                 (7)


 


Доказательство.  Путь  интегрирования  L  состоит  из  отрезка  прямой  t2=(t1t0)t2/(t1t0),  соединяющей  точки  (t0,0)  и  tÎН00.  Так  как  в  области  Н00  имеет  место  условие:  Re(l1(t1,t2))=sint1(cht2+asht2)£–c<0,  то  имеем: 


 


.


 


Справедливость  этой  оценки  также  можно  доказать  и  интегрируя  по  частям  J(t,e).  ¨


Лемма  2.  Если  tÎH10,  то  для  интеграла  J(t,e)  справедлива  оценка


 


|J(t,e)|£ce1–g,  где  0£g<1/2.                                     (8)


 


Доказательство.  Путь  интегрирования  L  тот  же  что  в  лемме  1.  Интегрируя  J(t,e)  по  частям  имеем:  |J(t,e)|£e|h1(t)/l1(t)|+O(e).


Учитывая,  что  l1(t)=O((t+p/2+ia))  и  сeg£|l1(t)|£Ceg,  то  следует  оценка: 


 


|eh1(t)/l1(t)|=O(e1–g).¨


 


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


 


,


 

где:  L=L1ÈL2ÈL3J(t,e)=j1(t,e)+j2(t,e)+j3(t,e), 


L1  —  отрезок  линии  t1=–p/2,  –a£t2£0; 


L2  —  отрезок  линии  t2=  –a,  –p/2£t1£t1


L3  —  отрезок  линии  t1=t1,  –a£t1£t2;


Вычисляя  эти  интегралы,  получим  оценки:


 


;  (9)


;  (10)


 


1.  Пусть  –p/2£t1£–p/2+de1/2,  0£t2+a£de1/2,  тогда  |j3(t,e)|£ce1/2.        (11)


2.  Пусть  |t1|£p/2–de1/2,  тогда  |j3(t,e)|£ce.              (12)


3.  Пусть  p/2–de1/2£t1£p/2,  0£t2+a£de1/2,  тогда  |j3(t,e)|£ce1/2.                     (13)


4.  Пусть  p/2–de1/2£t1£p/2,  de1/2£t2+a£2a,  тогда


 


.                    (14)


 


Из  этих  оценок  следует  справедливость  следующих  лемм.


Лемма  3.  Если  tÎH20ÈH21ÈH22,  то  для  интеграла  J(t,e)  справедлива  оценка:  |J(t,e)|£ce1/2.


Доказательство  следует  из  (9)—(11),  (13),  (14).


Лемма  4.  Если  tÎH12ÈH13ÈH14ÈH15,  то  для  J(t,e)  справедлива  оценка:


|J(t,e)|£ce1–g.


Доказательство  следует  из  (9),  (10),  (12),  (14).


Лемма  5.  Если  tÎH01,  то  для  J(t,e)  справедлива  оценка:  |J(t,e)|£ce.


Доказательство  следует  из  (9),  (10),  (12),  (14).


Следовательно,  справедлива 


Теорема.  Пусть  выполняются  условия  U1U2.  Тогда  задача  (1)-(2)  имеет  единственное  решение  и  для  него  справедлива  асимптотическая  оценка:


 


,


 


где 


 


 


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


1.Алыбаев  К.С.  Метод  линий  уровня  исследования  сингулярно  возмущенных  уравнений  при  нарушении  условия  устойчивости:  Дис.  …  д-ра  физ.-мат.  наук:  01.01.02.  Жалалабат,  2001.  —  203  с.


2.Нейштадт  А.И.  Асимптотическое  исследование  потери  устойчивости  равновесия  при  медленном  прохождении  пары  собственных  чисел  через  мнимую  ось  //  Успехи  мат.  наук.  —  1985.  —  Т.  40,  —  Вып.  5.  —  С.  300—301.


3.Турсунов  Д.А.  Асимптотика  решений  сингулярно  возмущенных  уравнений  в  случае  смены  устойчивости,  когда  собственные  значения  имеют  n-кратный  полюс:  Дис.  …  канд.  физ.-мат.  наук:  01.01.02.  Ош,  2005.  —  110  с.


4.Турсунов  Д.А.  Асимптотика  решения  сингулярно  возмущенной  задачи  с  периодической  точкой  поворота  //  интернет  журнал  ВАК  КР.  [Электронный  ресурс]  —  Режим  доступа.  —  URL:  http://195.38.189.154:81/jurnal/  (дата  обращения  27.08.2013).

Проголосовать за статью
Дипломы участников
У данной статьи нет
дипломов

Оставить комментарий

Уважаемые коллеги, издательство СибАК с 30 марта по 5 апреля работает в обычном режиме