Главная страница  Приемы качественного исследования 

[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [ 91 ] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163]

три бифуркационные поверхности: две - соответствуюгцие двукратным циклам (охватывающим цилиндр) и одна - сепаратрисе, идущей из седла в седло.

Пример 2 (исследование уравнения движения самолета методом Понтрягина [33]). Введем в уравнения движения самолета (см. § 3 гл. 14) новое переменное г/ = Ур и малое р, полагая

Тогда система (1) § 3 гл. 14 может быть записана в виде dffldt = г/2 - cos ф = дН/ду,

dy/dt = - г/ зшф -f р(А; - г/2) =, дН/д + л/(г/),

(ф. г/)=г/3-г/созф = /1,

т. е. при р = О мы получаем консервативную систему.

По смыслу задачи рассматриваются лишь значения г/ > 0. Траектории этой системы изображены на рис. 149. Замкнутым

кривым, охватывающим состояние равновесия, соответствуют значения /г < О, а охватывающим цилиндр - значения h> 0.

Самому состоянию равновесия типа центра соответствует значение /1 = -2/3. В точках (-я/2, 0), (л/2, 0)- седла.

При \х¥=д в точках (-л/2, 0), (л/2, О)-по-прежнему седла, и среди сепаратрис этих седел, как и при (X = О, есть части оси у = 0. Кроме того, в этом случае, очевидно, существует еще одно состояние равновесия, координаты которого г/3, фз удовлетворяют соотношениям [у Ф 0)


Отсюда

г/2 -cos ф = О, - 8тф-Ьр(А; -созф) = 0.

2(i +iif

(6) (V)

Так как созфз= г/з> О, то, очевидно, в выражении (7) мы должны взять знак + перед корнем (при малых р будет iik- КО и подкоренное выражение больше 4р,*А;2). Таким образом, после элементарных преобразований мы получаем для третьего состояния равновесия Оз(фз, Уъ)

cos Фз = 2-(8)



Ф = arccos

то, очевидно, мы имеем, как нетрудно видеть,

{h,k) = 2 \ (А;-Зг/=) arccos 3 dt/, (И)

где ei{h) и e2(/i) -корни уравнения

3-Зг/=.3/1, -2/3/К О

(при этом 1 = ei УЗ). 18*

В силу того, что при р = 0 состояние равновесия Ог - центр, т. е. для него А > О, о = О, то, очевидно, и при достаточно малых ц для состояния Оз будет А > О, т. е. Ог - фокус.

Отметим еще, что для значений параметра, при которых сепаратриса идет из седла в седло, одновременно образуется два замкнутых контура из сепаратрис (среди которых есть части оси ф): контур, охватывающий состояние равновесия, и контур, охватывающий цилиндр. От первого из этих контуров может появиться предельный цикл, охватывающий состояние равновесия, а из второго - охватывающий цилиндр.

В согласии с методом Понтрягина для решения вопроса о циклах, охватывающих состояние равновесия, рассматриваем функцию!

(h, А;) = j j (р; + qy) йф = j j (А; - 3/) d(p dy, (9)

где двойной интеграл распространен на площадь, ограниченную кривой Я(ф, y) = h при - 2/3 < /i < 0.

Для решения вопроса о циклах, охватывающих фазовый цилиндр, рассматриваем выражение

Ah,k)== J y{k-y)dц> = 2]y{k-y)dц>, (10)

-я о

где у определяется из уравнения Я(ф, y) = h (см. гл. 12) при h>Q.

Функции if)i(/i, к) и if)2(/i, к) мы доопределим по непрерывности до значения h = Q, соответствующего сепаратрисе (консервативной системы), и доопределенную таким образом функцию будем обозначать через 1)(/г, к). Все дальнейшее посвящено изучению возможного характера функции il)(A, к) при разных к.

Мы можем записать выражения для if)i(/i, к) и if)2(, к) в одинаковом виде.

Действительно, так как из (5)



Геометрический смысл ei{h) и е2{/г) (-2/3 </г < 0) представлен на рис. 150.

Интегрируя стоягций в выражении (И) интеграл по частям, мы получим

2 {к - у) у arccos

• (fe-y)(2/ + 3)

Vdy--im-yf

,{h,k)==2[L=;MMdy.

К аналогичному виду преобразуется и выражение для "2(/Ji, к): Действительно, из (5) мы имеем

2у -L 3h

Очевидно, получаем

у Vdi- ~ (Зк - yf

,h (h 1Л 9 C - /) (2y -f- 3h) , (/г, /с) = 2 === dy

J V Чу- ~{Ш-уУ

(13)

(здесь e* - положительный корень уравнения у - Зу = 3h (1> ]3), а 62 - положительный корень уравнения г/3 + 3г/= 3/г,


Рис. 150

которое получается, если в Я(ф, y) = h подставить ф = л). Зависимость корней ei, 62, eg от h изображена на рис. 151. (Часть кривой e\{h) (начинающейся в точке А), лежащая слева от оси ординат, соответствует значениям h, - 2/3 < /г < О, т. е. замкнутым кривым, охватывающим центр; часть этой кривой, лежащая справа от оси ординат, соответствует /г > О, т. е. кривым, охватывающим цилиндр (см. рис. 151).)




[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [ 91 ] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163]

0.0128