Главная страница Приемы качественного исследования [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 |