книги / Теория оптико-электронных следящих систем
..pdfхарактеризующим быстродействие ОЭСС, и минимизировать его. В качест ве методики, позволяющей довести вычисления до конца в общем виде, можно воспользоваться интегральной квадратической оценкой
/о = 7A2(t)dt,
о
описанной А.ДлКрасовским в [89]. Можно пойти и более непосредствен ным путем, если рассмотреть в качестве критерия быстродействия время окончания переходного процесса Г, удовлетворяющее условию
i Д(Г) | < 0,05 Д(0) при t > Т. |
(2.3.15) |
Однако такой критерий уже не приводит к конечным выражениям и пред полагает использование ЭВМ с процедурой поиска минимума, описанной в § 3.6, посвященном идентификации математических моделей.
Заметим здесь, что быстродействие (а следовательно, и полоса пропус кания замкнутого контура) ОЭСС ограничивается сравнительно небольшим крутящим моментом. Гораздо большими возможностями в этом смысле обладает гидравлический объемный привод. Высокое быстродействие его достигается благодаря высокому давлению рабочей жидкости (до 2000 кПа). При этом снижаются вес и габариты исполнительного элемента. Математи ческое описание гидропривода по форме не отличается от уравнений вида (2.3.1).
Индикаторный гиропривод. Обратимся теперь к наиболее распростра ненному виду исполнительных элементов —к гироскопам, работающим в режиме управления ОЭСС. Гироскопические исполнительные элементы обладают значительными преимуществами перед электрическими и гид равлическими: помимо малой инерционности и высокой точности работы они обладают принципиальной возможностью развязки следящего коор динатора по отношению к подвижному основанию и, следовательно, свой ством ’’запоминания” положения оптической оси прибора в инерциальном пространстве. Простейшей разновидностью исполнительных элементов это го типа является трехстепенный гироскоп (рис. 2.8). Конструктивно он выполняется обычно таким образом, что его ротор 1 несет на себе и один из элементов оптической части (например, рефлектор). Эллипсоид инер ции ротора сплюснут по оси х. Вместе с внутренней рамкой 2 ротор может поворачиваться вокруг оси Ог\. Внутренняя рамка может вращаться вмес те с наружной рамкой 3 вокруг оси O f, подшипники которой жестко связа ны с подвижным основанием.
Ввиду симметрии эллипсоида инерции все оси, проходящие через точку О и перпендикулярные к оси Ох> являются главными. Момент инерции относительно любой из этих осей равен /. Момент инерции относительно оси Ох обозначим через / р. При движении гироскопа точка О остается неподвижной в системе координат Oxyz, и поэтому любое движение его
можно |
рассматривать как |
вращение с мгновенной угловой скоростью |
||
со = (г, |
со'). Величину Я = / рг называют |
собственным кинетическим мо |
||
ментом ротора. |
|
|
|
|
Спроектируем |
вектор со |
на оси у, и z |
правой связки Oxyz (см. рис. |
|
2.8). Тогда 1р и lq |
будут проекциями вектора кинетического момента 0_. |
Угловые скорости р и q будут компонентами угловой скорости трехгран-
111
Рис. 2.8. Исполнительный элемент гироскопического типа с одним ротором: 1 —ро тор; 2 - внутренняя рамка; 3 - наружная рамка
ника Oxyz на оси у и z. Проекщш угловой скорости трехгранника на оси Ох обозначим г'. Она необязательно равна г, т.е. трехгранник Oxyz вращает ся вокруг оси Ох независимо от ротора /.
Согласно теореме о кинетическом моменте [18]
dQ |
= М, |
(2.3.16) |
— |
||
dt |
|
|
где М —момент внешних сил, приложенных к гироскопу. |
|
|
В |
проекциях равенства (2.3.16) на оси х, у, z кроме |
составляющих |
в Х9 0у и 6Z мы должны учесть эффект вращения вектора Я относительно осей у и z (см. рис. 2.8), а также векторов В_у и _0Z относительно оси Ох. Таким образом, получаем систему скалярных равенств
ddz
~+ H q - r ' d y = M z ,
dQy |
|
(2.3.17) |
— ~ ~Н р + г’в2 =Му, |
||
dt |
|
|
dt |
Р@у — х Мх . |
|
|
|
|
Подставляя в (2.3.17) значения 6у =Ip,6z = Iq и в х =Я, имеем |
|
|
7( ~ |
|
|
|
|
(2.3.18) |
dH |
|
|
---- - M z. |
|
|
dt |
z |
|
112
Рассмотрим более подробно физическую природу моментов, действую щих на гироскоп, управляемый по осям Of и Ог\. Заметим, что если оси Оу и Oz трехгранника Охуг совместить с осями Of и От? карданова подвеса, то г' = 0. Кроме того,
My =M ynp- M Tf 9 М2 =М упр- М Р . |
(2.3.19) |
Здесь М упр и М \ир - моменты, развиваемые датчиками под действием управляющих сигналов с каналов фазовых детекторов (см. структурные схемы на рис. 2.1 и 2.2) либо с усилителей мощности после вычитания сиг налов с площадок (см. рис. 2.5); МуР иМ \р —моменты сил сопротивления в подшипниках осей вращения рамок.
Для простоты исследования в пределах линейных дифференциальных уравнений, позволяющего сделать общие качественные выводы об устой чивости процесса сопровождения, положим 1)
|
M F = iq |
(2.3.20) |
и, считая г |
= 0 , получим уравнения движения гиропривода в виде |
|
dp |
|
|
1~^+% р+Н я=М упР, |
|
|
dt |
у |
|
I - - + £ < ? - Яр =Л/упр, |
(2.3.21) |
|
dt |
у |
|
dH
— = МХ. dt
Из этих трех уравнений лишь два первых относятся к системе управления. Третье уравнение характеризует вращение ротора во время раскрутки гиромотором либо его торможение под действием сил трения в подшип никах рамки 2 (см. рис. 2.8).
Моменты управления в режиме слежения зависят от углов рассогласо
вания. В соответствии с этим уравнения |
(2.3.21) дополняются условиями |
замыкания |
|
M lnp= - k TAy = - k T f q(r)dr, |
|
о |
(2.3.22) |
М**р =ктАг =кт/ p(f)dT. |
|
Представляется интересным вопрос о предельно допустимой доброт ности кт. Для выяснения этого значения запишем характеристическое урав нение, вытекающее из (2.3.21) и (2.3.22) :
(А+£)Х |
—(АХ + кт) |
(2.3.23) |
|
ЯХ +£Г |
(/Х+{)Х |
||
|
1) В действительности моменты сил сопротивления описываются нелинейной зави симостью. В первом приближении можно пользоваться зависимостью вида Му р =
=MyQ signp и Af^P =MZQs\gnq для более точного исследования по сравнению с (2.3.20).
8. Ю.М. Астапов |
113 |
или |
|
/ 2Х4 + 27£Х3 + (Я2 +$*)\ + 2kr H \ + k l = 0. |
(2.3.24) |
Составив определитель Гурвица четвертого порядка с коэффициентами полинома (2.3.24), получаем условие, ограничивающее добротность гиро скопического привода:
(2,3.25)
Так как величина %очень невелика, то для повышения добротности приво да ОЭСС следует увеличить отношение #//. Это достигается за счет боль шой скорости вращения ротора (30000 —60000 об/мин). В свою очередь это требует повышенной частоты питающего напряжения (500 — 1000 Гц).
Гиропривод с разгрузкой от внешних возмущений. Трехстепенный гиро скоп не способен развивать больших управляющих моментов. В то же вре мя в ряде конструкций ОЭСС вращающаяся часть может обладать значи тельной массой (например, при тяжелых телеобъективах). Внешние пере грузки со стороны подвижного основания при несбалансированности подвеса отклоняют ось оптической системы, в результате чего в управляю щем сигнале появляется ложная составляющая.
В настоящее время разработаны схемы гиропривода со специально предусмотренными мерами защиты против внешних силовых возмущений. Рассмотрим двухгироскопную схему гиропривода, которая называется гиростабилизированной платформой (рис. 2.9). Внутренняя рамка платфор мы 1 жестко связана с наружными рамками 2 и 3 гироскопов Гi и Г2. Под действием коррекционного момента, развиваемого датчиком Д1э платформа прецессирует вокруг оси у к . Момент датчика Да заставляет платформу прецессировать вокруг оси zK. При этом перемещаются тяги Ту и Т2, связанные с оптическим блоком, который подвешен в собствен ном кардановом подвесе. Таким образом осуществляется сопровождение излучающего объекта. Возмущающие внешние моменты разворачивают
Рис. 2.9. Двухроторный гироскопический исполнительный элемент: 1 —внутренняя рамка; 2, 3 - наружные рамки; 4, 5 —внутренние рамки
114
гироплатформу. При этом гироскопы стремятся сохранить неизменным направление осей вращения в инерциальном пространстве, и на датчиках П1 и П2 углов поворота внутренних рамок 4 и 5 возникают сигналы, пос
тупающие на усилители разгрузки У! |
и |
У2. |
Усиленные сигналы подаются |
на моментные датчики разгрузки |
и |
Р2. |
В результате осуществляется |
компенсация внешних возмущений.
Работа исполнительных элементов гироскопического типа основана на принципе прецессии. При этом, как мы видели в уравнениях (2.3.21), в каналах управления появляются обратные связи. Из-за различия в знаках эти связи называются антисимметричными [51, 63]. В многороторных гироскопических приводах эти связи значительно ослаблены. В дальней шем будем рассматривать общий случай гиропривода силового типа с пере даточной функцией вида
izO) |
_ цуО) |
________ К ________ |
(2.3.26) |
|
Mz (s) |
Му <хГ |
(1 +2 %Ts + T 2s2)s |
||
|
Здесь ipz (s) и ipy (s) — изображения углов поворота оптического блока ОЭСС (под действием тяг Т2 и Т^, см. рис. 2.9); Mz (s) и My (s) - изоб ражения моментов датчиков Д! и Д2 ; кт— коэффициент усиления гиро привода; £ — коэффициент демпфирования; Т — постоянная времени.
§2.4. Анализ периодических режимов в ОЭСС
снепрерывным управлением
Современные ОЭСС представляют собой сложные технические устрой ства с многочисленными связями между каналами управления. Многие элементы описываются нелинейными характеристиками. Совокупность других элементов описывается линейными дифференциальными уравне ниями высокого порядка. Основным методом анализа и конструирования ОЭСС служит математическое моделирование. В то же время для выявле ния закономерностей общего характера, определяющих динамику ОЭСС, для обзора пространства параметров элементов и возможностей выполне ния требований по качеству сопровождения необходимо иметь в распоря жении аналитические методы. Анализ проводится обычно при некоторых упрощениях математической модели ОЭСС, не вносящих существенных изменений в качественную картину процесса слежения. Разработка доста точно простой модели, допускающей математическое исследование, пред ставляет собой самостоятельную и не всегда выполнимую задачу. Поэтому часто задача исследования ограничивается созданием модели, характерис тики которой совпадают с характеристиками реальной ОЭСС лишь при не которых режимах работы.
В настоящем параграфе мы рассмотрим упрощенную модель ОЭСС с по вышенной точностью слежения и с силовым гироприводом, обладающим выраженными колебательными свойствами. Обозначим компоненты угло вого рассогласования А через Aj и Д2:
At = Acos<р, Д2 = Asinip. |
(2.4.1) |
8* |
i 15 |
Учитывая условия замыкания y z = - А! и \ру = —А2, из (2.3.26) получаем
Т2Д(13) + 2$ТА(2} + A(!1} = - к тМг, |
|
Т2Д(23) + 21ГД22) + Д(21} = —к тМу . |
(2-4.2) |
Моменты Mz (t) и M y (l) в правой части (2.4.2) |
пропорциональны напря |
жениям Mi (г) и м2(0 в каналах управления, которые находятся из выра жений (2 .1 .12 ), т.е.
Mz (t) —kj^U\(t) —АгдUm2LXm (A, ^)cos(<^ + e),
My (t) = k au2(t) = k aUmaxm(A, <p)sinQp + e),
где £д —коэффициент преобразования моментного датчика.
Формальные условия возникновения кругового предельного цикла.
Следуя [51], сделаем допущение о возможности представления функции кткя CJmахm (A, ip) в виде произведения F(А)А. Это допущение основано на предположении о гармонических колебаниях, которые могут возник нуть в замкнутом контуре ОЭСС с непрерывным управлением. Такие коле бания действительно наблюдаются при экспериментальных исследованиях и при математическом моделировании подобных систем. Поэтому в даль нейшем рассмотрим уравнения ОЭСС в виде
Т2Д(3) + 2£ ГД(.2) + Д<‘> = -F(A)Acosfc + е),
(2.4.4)
Т2 А(3) + 2^ТА(2) + Д(21} = -^(Д)Дап((р + е) |
|
|
или, учитывая обозначения (2.4.1), в виде |
|
|
Т2Д(,3) + 21ГДj2) + Д « = -F (A )(A lCose - |
A2sine), |
|
|
|
(2.4.5) |
Т 2А{3) + 2%ТА<2) + А{2) - -F (A )(A 2cose + Ajsine). |
||
При идеальных гармонических колебаниях |
|
|
Аг = AQCOSCOO^, |
А2 = A0sinco0^ |
(2.4.6) |
Здесь А0 = const и со0*~ const - амплитуда и частота этих колебаний. Заме |
||
тим, что согласно (2.4.1) и (2.4.6) |
|
|
Д= \M i + д]' |
= Д0- |
|
Поэтому в уравнениях (2.4.5) также и F(А) = F(А0) = const.
Для формальной проверки возможности существования гармонических
колебаний в ОЭСС подставим |
соотношения (2.4.6) в уравнения (2.4.5). |
При этом мы получим формальные тождества |
|
со0А0(Т2 col - l)sinco0f —2%TcooA0cosG)0t = - i 7(A0)A0cos(co0^ + e), |
|
|
(2.4.7) |
-co0A0(T 2 col - l)coso>o^ ~ |
A0sincoo/ = -F(A o)A 0sin(a;or + e). |
Так как A0 Ф 0, то (2.4.7) можно сократить на А0 и, умножая первое ра венство на sina>o^ а второе на cos сo0t и вычитая, получить
со0(Т 2со1 - 1) = i7(A0)sine. |
(2.4.8) |
Аналогично, умножая первое соотношение на cosco0t, второе на sinco0/
116
и складывая, имеем
2%Тш\ =F(A0)cose. |
(2.4.9) |
Соотношения (2 .4 .8) и (2.4.9) образуют систему, из которой можно найти частоту и>0 и амплитуду А0 предполагаемого периодического движе ния. Действительно, при со0 Ф 0, поделив (2.4.8) на (2.4.9), получаем
--------------Т 2ш20 - |
1 - tge |
|
2 |Гсо0 |
|
|
откуда находим |
|
|
w0 = ~ |
(£tge± x / F tp e + l) . |
(2.4.10) |
Амплитуду А0 после этого можно найти либо из (2.4.8), либо из (2.4.9). Можно записать и общее выражение для определения F(A0), откуда при за данной зависимости F(А) находится амплитуда А0. Например, умно жая (2.4.8) на sine, а (2.4.9) на cose и складывая, получаем
F(A0)= 2?71coocose + (T 2col - l)cj0sine. |
(2.4.11) |
Формулы (2.4.10) и (2.4.11) представляют собой искомые выражения для вычисления частоты и амплитуды автоколебаний. Остается выяснить, ка кие из значений этих параметров соответствуют реально наблюдаемым процессам.
Заметим, что устойчивое периодическое движение (автоколебания) может возникнуть лишь при неустойчивости положения равновесия, кото рому соответствует тривиальное решение системы (2.4.5) A^f) = Д2(0 = 0. Анализ устойчивости положения равновесия сводится к исследованию зна ков корней характеристического уравнения системы (2.4.5). Наиболее просто это уравнение записывается, если в (2.4.5) перейти к изображе ниям. Сохраняя символы A^s) = L{Aj(f)} и A2(s) = L {A 2(t)}9 запишем
[ T V + 2 %Ts2+s +F(0)cose]A1 (s) |
- A 2(^)F(0)sine |
|
Ai(s)F(0)sine |
[ T V + 2£7s2+ s + F(0)cose]A2(s) |
|
Отсюда,при X = s получаем искомое характеристическое уравнение |
(2.4.12) |
|
|
||
[Т2Х3 +2£ГХ2 +F(0)cose]2 +F 2(0)sin2 e= 0. |
(2.4.13) |
Непосредственное применение к уравнению шестого порядка (2.4.13) процедуры Гурвица приведет к весьма громоздкой системе неравенств. Поэтому следует использовать другой прием, позволяющий понизить порядок определителя Гурвица. Запишем (2.4.13) в форме
(Г 2Х3 +2 £ГХ2 +Х)2 +,2F(0)(r2 X3 +2 £ГХ2 + X)cose + F 2(0)=0.
(2.4.14)
Корни этого уравнения удовлетворяют условию
Г 2 Х3 + 2 £ГХ2 + Х = —F(0)cose± x/F 2(0)cos2e - F 2(0) =
= —i^Ctycose ± z‘F(0)sine. |
(2.4.15) |
117
Из (2.4.15) видно, что корни X представляют собой комплексные числа. Подставив в (2.4.15) X = а + //3 = а + /оо0 и приравнивая вещественные части, получаем
Т 2(а3 - |
Заооо) + 2£Г(а2 - |
+ а = —F(0)cose. |
(2.4.16) |
|
Очевидно, при а < |
0 положение равновесия устойчиво и автоколебаний |
|||
в системе |
не возникает. Записав определитель Гурвица |
для уравне |
||
ния (2.4.16), получаем условия устойчивости положения |
равновесия |
|||
|
< F(0) < |
2 ? |
|
(2.4.17) |
cose |
(1 -2Гы?). |
|||
|
Тcose |
|
|
Из (2.4.17) во всяком случае видно, что при F(0) < 2£/rcose положе ние равновесия неустойчиво.
Рассмотрим теперь условия устойчивости в ОЭСС с непрерывным управ лением кругового предельного цикла, описываемого равенствами (2.4.6). Для этого введем отклонения от этого предельного цикла
8i = Aj —A0 cosco0r, |
62 = А2 —А0 sin со0^ |
(2.4.18) |
|
Подставив в уравнения |
(2.4.5) вместо |
и Д2 их значения из |
(2.4.18) |
и учитывая выполнение тождеств (2.4.7), получим |
|
||
Т25 Р> + 2£Т8{? ) 5f >= —F (Д0) (5, cos е - |
52 sin е), |
|
|
T2di3) + 2£Г5£° + |
= -F (A o) (5! sin е + 62 cos е). |
(2.4.19) |
|
Эти уравнения по своей |
структуре не отличаются от уравнений |
(2.4.5), |
и условия устойчивости кругового предельного цикла сводятся к усло
виям устойчивости решения бi (г) = 62 (О — 0- Поэтому |
на основании |
|
(2.4.17) можно сразу записать эти условия в виде |
|
|
2*Гсо$ |
2% |
(2.4.20) |
-- ------ < F (A 0) < |
(1 2T ul). |
|
cose |
Г cose |
|
Здесь параметры оо0 и А0 предельного цикла должны удовлетворять урав нениям (2.4.7).
Метод малого параметра. Б ы с т р ы е и м е д л е н н ы е д в и ж е н и я . Мы рассмотрели две разновидности возможных движений в системе (2.4.5). Для более широкого исследования возможных движений в ОЭСС с непре рывным управлением воспользуемся представлениями, связанными с асим птотическими методами теории нелинейных систем (см. [3], где приводит ся метод малого параметра для систем, описываемых уравнениями второ го порядка).
При использовании асимптотических методов существенным фактором является малая величина коэффициента демпфирования %в передаточных
функциях гиропривода (2.3.26). Идея метода Ван-дер-Поля |
заключается |
в представлении уравнения второго порядка в форме |
|
х +х =д/(х, х), |
(2.4.21) |
где 1л —малый параметр; /(х , х) —некоторая, вообще говоря, нелиней ная функция своих аргументов. При ц = 0 уравнение переходит в урав нение гармонического осциллятора, решение которого x(t) = a cost +
118
+ bsint легко находится. При малых значениях параметра д решение нели нейного уравнения (2.4.21) должно мало отличаться от движения гармо нического осциллятора.
Следуя этому методу, запишем уравнения (2.4.5) в форме
д<3>+ |
F ( А) |
(Д1 cos е —Д2 sin е), |
|
||
Г2 |
|
||||
|
|
|
|
(2.4.22) |
|
|
|
|
Н Д) |
|
|
д?> + |
(Д1 sine + Д2 cose) |
|
|||
|
|
|
г2 |
|
|
или, введя обозначения |
|
|
|||
1/Г - м , |
|
|
|
||
- |
у |
д 12)- |
^ ~ ( Д 1 С05е - Д 2 зте) = дФ1 (Д1 ,Д 2 ,Д ? )), |
(2.4.23) |
|
- |
2% |
ы |
F{Д) |
т |
|
— Д2 - |
sine + Д2 cose) = дФ2(Д!, Д2, Д2 ), |
|
в форме
Д(!3) + СО2 Д(,0 = цФг (Д!, Д2, Д<2)),
Д^3) + со2Д ^ = д Ф2 (Д1 , Д2, Д^2)).
При д = 0 легко находим порождающее решение уравнений (2.4.5)
А х =Су +ai cos(oof + ipw), A2 = 02 + 02 cos(oof + <p2o)>
где Ci, c2>0i, 02, ^ю» ^20 —постоянные, определяемые начальными ус ловиями.
При д Ф 0 эти постоянные становятся медленноменяющимися функция ми времени, скоростью изменения которых можно пренебречь при их диф ференцировании. Таким образом, полагая для краткости cot + ={Pi (0> будем иметь
А/ = Ci(t) + ai (Г) cos ^ (f),
Ар* = - |
aj it) соsin iff (t), |
(2.4.25) |
Ap) = - |
Qf (t) со2 cos (r) |
(z = 1 , 2 ). |
Оценим допущенную погрешность при записи (2.4.25). Дифференцируя формально первое равенство (2.4.25), получаем точное соотношение
Ар} (г) = ср} if) + (г) cos& (f) - 0/ (Г) ^(1) (t) sin ^ (г). Приравнивая это выражение второму условию (2.4.25), находим
С,(1) (О + Я,(1) (0 cos<A-(0 - щ (0 ftp (?) sin (ft- (?) + at (?)cosing (?) = 0.
(2.4.26)
Аналогично, дифференцируя формально второе равенство (2.4.25) и при равнивая результат третьему условию (2.4.25), имеем
- 0р* (Г) соsin (г) —at (t) |
(f) соcos (t) + а%(t) со2 cos (t) = 0. (2.4.27) |
119
Наконец, дифференцируя третье равенство (2.4.25) и подставляя резуль тат в уравнения (2.4.24) с учетом приближенного выражения второй произ водной Д ^ (t ), получаем
- яр* (f) cos & (t) + я,-(г) |
(t) sin & (t) —я,-(;t) соsin (t) = |
Ф,-. (2.4.28) |
CO
Соотношения v2.4.26) — (2.4.28) позволяют оценить скорости измене ния функций Cj(t), я,-(г) и Vi(t). Действительно, из (2.4.26) и (2.4.28) получаем
(0 = 4 ф<(д 1 . |
)• |
(2.4.29) |
СО
Далее, умножая (2.4.27) на sin ^ j(r), (2.4.28) —на c o s (г) и складывая, видим, что
в/(1)(0 ■=- —2ф/ (Д1 , |
) cos (0. |
(2.4.30) |
СО
Умножая (2.4.27) на c o s ^ f ) , (2.4.28) - на sin^/(r) и вычитая, получим оценку для <p£l) (t) :
*> f40= w + 4 |
(Ai »Да. Д?4 ) *“ « (') . |
(2.4.31) |
СО
Таким образом, процессы ct (t) и я,* (г) не содержат в правых частях (2.4.29) и (2,4.30) параметра со и их следует считать медленными. Процесс ipf (t) следует считать быстрым с медленноменяющимся вторым слагаемым в (2,4.31).
Следуя далее методике Ван-дер-Поля, осредним правые части уравне ний (2.4.29) - (2.4.31) на периоде быстроменяющихся переменных (t) и кр2 ( 0 , равном 2 я/со:
/14 |
|
ц |
27Г |
27Г |
|
|
||
4 |
}(0 = |
Т~т~2 |
о |
f d * , |
|
|||
|
|
4ж |
со |
|
о |
|
|
|
|
/,Ч |
, |
и |
2п |
2п |
f <t>iCos<Pid<p2, |
(2.4.32) |
|
в / ( 0 = |
2 |
2 |
|
/ |
||||
|
|
4ж |
со |
о |
о |
|
|
|
|
1)(0 = ш + |
|
|
|
|
2тг |
|
|
^ |
|
|
|
|
/ Ф,-sin<#d<p2 . |
|
||
|
|
|
47Г2 со2 |
|
0 |
|
Подставив в (2.4.32) явные выражения функций Ф,- (2.4.23), получим укороченные уравнения типа Ван-дер-Поля [3].
Исследование решений системы (2.4.24) предполагает, что составляю щие вектора Д углового рассогласования представляются в виде разло жений
A i(0 = ci(0 + ai(0cos<p1 (r)+ |
2 |
ц к Uk, |
|
к= 1 |
|
|
оо |
(2.433) |
A 2(r)= C2( 0 + ^2(0cOS^2( 0 + |
2 |
lik Vk , |
|
к~ |
1 |
120