Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория оптико-электронных следящих систем

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
19.38 Mб
Скачать

характеризующим быстродействие ОЭСС, и минимизировать его. В качест­ ве методики, позволяющей довести вычисления до конца в общем виде, можно воспользоваться интегральной квадратической оценкой

/о = 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). Тогда и 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 * ,

 

 

 

со

 

о

 

 

 

/,Ч

,

и

2п

2п

f <t>iCos<Pid<p2,

(2.4.32)

в / ( 0 =

2

2

 

/

 

 

со

о

о

 

 

 

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