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

[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]

(т-Аф[1-,г(2Р-Ьг) У]

и проинтегрирована.

Состояние равновесия Ог - всегда седло. Седловая величина будет иметь значение

{Р + Qy)2 = - (3)

Существенно, что выражения Рц, + Qy для фокуса и для седла имеют одинаковое значение. Как будет показано, это определяет существенные особенности структуры разбиения фазового пространства и пространства параметров.

3. Предельные циклы. Если систему (1) записать в виде

d(p/dt = Ну, dy/dt = - Яф -Ь [iq (ф, у), (4)

где (ф, г/) = [-а -Рс(2ф) + 1Г5(ф)]г/, то значения константы ho, выделяющие кривые Сн консервативной системы, вблизи которых при малом р будут предельные циклы системы (4), определяются как корни уравнения г)(/г) = 0, где

{h)= f j" йф йг/ = - j" j [а -Ь рс (2ф) - ys (ф)] йф dy =

= - j[a-b рс(2ф)-75(ф)]г/йф. (5)

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

а) Пусть 2 ->2<Т<1. При изменении h в интервале -я(1 - Т)/2 < < -л(1 - Т)/4 кривые Си будут эллипсами, целиком расположенными в полосе (О, я/2), а интервалу (-я/4) (1 - Т)2 < < О будут соответствовать замкнутые кривые, сшитые из кусков эллипса и гиперболы, расположенных соответственно на полосах (О, л/2) и (я/2, я).

2. Состояния равновесия, Седловая величина. Система (1) имеет при Т<1 два состояния равновесия: Oi{nT/2,0) и 02(я - яТ/2,0). Состояние равновесия Oi будет неустойчивым фокусом, если для точки Oi будет

(Рф + Qv)i = - (i (а + Р - 2рГ - 7Г) = - fia > О, а < О,

устойчивым фокусом, если а > О, и центром, если а = 0. В последнем легко убедиться: при условии а = а + р - 2рТ - = О система (1) на полосе Офя/2 может быть представлена в виде



Выполняя интегрирование, ползгчим

-oFgQi), -aKh<-a/2,

-aFh), -a\/2<h<:0.

Здесь

= 4- (1 - Л% 0 (h) = /2я/ (a + h).

Fi (/i) = /2 [a Y2h + +

+ (fe + ai Имеем также

F;(fe) = /2

- + arcsm

feln >+l +

2 /2(й + а);

я , . a , , l/2ft + a -1- a

-УГ + arcsm , == + In-, -

2 /2(Л + а2) У-2Й

() = -ЙГГ?Г~

F, (- = Yn/2abi, F[ (- al2) = Y2n\ Fli- a/2) = 0.

Таким образом, функция tfi(fe) сшивается из линейной функции -oFq на интервале -а < fe < -а/2 и монотонной на интервале -а?12 < fe < О, имеюгцей с функцией -oFq одинаковый знак производной. Функция itfe) будет убывающей при а > О, возрастающей при о < О и будет тождествеппо обращаться в пуль при а = 0.

Предельных циклов в рассматриваемом случае 2 ->27<1 нет. При 0 = 0 существует внутри петли сепаратрисы континуум замкнутых кривых. При о О характер устойчивости фокуса одиозиачио определяет структуру разбиения фазового цилиндра на траектории. На рис. 205, 1, 1-4, 4 изображены структуры разбиения фазового пространства соответственно для случаев о < О, о = О и а > 0.

б) Пусть 3 -2У27<2 -У2. Фазовые траектории консервативной системы располагаются при этом не более чем на трех полосках (не выходят на полосу (-я, -я/2)). Положим также для определеппости, что Т > 0,5 (случай Т < 0,5 апалогичеи случаю Т > 0,5, по отвечает другой иоследовательпости интервалов изменения fe, соответствующих определенному числу сшиваний); для 7<0,5 будет

--{1-ТГ<:-~{Т-АТ + 2)<--2l(i f)2<0;



-2{h + a)y2a - Ь\я - 2arcsin ф£=

Имеем также

Ш = 4(2/г + b2)V2 2/2a2 - b (я-2arcsin -),

\ y2{h + a)/

ft -j- a

(- b/2) = (- b/2) = (- bV2) = 0, F (- Ь2/2) > о.

Из последнего следует, что F2{h) на интервале ~Ь/2 < h 0 есть функция, принимающая положительные значения, монотонно возрастающая и обращающаяся в нуль в единственной точке h = -Ъ12.

Выражения (7) определяют г)() как непрерывную функцию, сшитую из трех кусков с различными аналитическими представлениями на интервалах {-а, ~а/2), {-а/2, -Ь/2) и (-Ь/2,0). Функция \p{h) на интервале {-а, -Ь/2) имеет знак, противоположный знаку а, обращаясь ири а = 0 в нуль тождественно. На интервале {-Ь/2,0) функция состоит из двух слагае-

мых, одно из которых положительно, а другое имеет знак, противоположный знаку а. Выбором а значение г)() на всем интервале (-Ь/2,0) может быть сделано как положительным, так и отрицательным.

Проследим за корнями функции (h) и изменением качественной структуры разбиения фазового пространства на траектории ири возрастании о от отрицательных значений. При а<0 функция if)(fe) на всем интервале (-а, 0) будет положительной.

ДЛЯ Т > 0,5 будет

(1 Г)2 < (1 J)2 <- (Г2 4у + 2)< 0.

Выполняя интегрирование в формуле (5) для случая 0,5 < <Т<2~12, получим

- aFo(fe), - a2<fe< -а2/2,

-{h) = - «71 №). - «/2 < < - Ь/2, (7)

12{Щ - oF,{К), -b/2h0.

Здесь 6 =--я (Т - 4-Ь 2), Fo{h} и Fi (/г) имеют значение, указанное в (6), и Fzih) будет F (h) = 2 {2а - b)V2h+T + --(2fe + bf/




[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.0121