Научная статья на тему 'Расчет цифровых следящих систем комбинированного управления с использованием предельной непрерывной модели. Ч. I. построение предельной непрерывной модели'

Расчет цифровых следящих систем комбинированного управления с использованием предельной непрерывной модели. Ч. I. построение предельной непрерывной модели Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
60
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЦИФРОВАЯ СЛЕДЯЩАЯ СИСТЕМА / ПРЕДЕЛЬНАЯ НЕПРЕРЫВНАЯ МОДЕЛЬ / DIGITAL TRACKING SYSTEM / LIMITING CONTINUOUS MODEL

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Ермоленко А. И., Коршунов А. И.

В общем случае произвольного порядка непрерывной части построена предельная непрерывная модель цифровой следящей системы (ЦСС) комбинированного управления, к свойствам которой неограниченно приближаются свойства ЦСС при повышении частоты измерения рассогласования. Получено аналитическое выражение для коэффициентов дискретного фильтра, вырабатывающего добавку к задающему воздействию, компенсирующую скоростную ошибку и ошибку по ускорению ЦСС. Коэффициенты дискретного фильтра выражены через параметры непрерывной части ЦСС, что позволяет не определять ее дискретную передаточную функцию.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Ермоленко А. И., Коршунов А. И.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Calculation of digital tracking systems of combined control using a limiting continuous model. Part I. Construction of a limiting continuous model

In the general case of continuous part of an arbitrary order, a limiting continuous model of digital tracking system (DSS) of combined control is constructed, so that the DSS properties approach indefinitely to the properties of the model with increasing frequency of the mismatch measurements. A discrete filter is considered which generates an additive to the master control, compensating for the speed error and the error in the DSS acceleration by variation of parameters of the continuous part of the system. An analytical expression for the filter coefficients is derived.

Текст научной работы на тему «Расчет цифровых следящих систем комбинированного управления с использованием предельной непрерывной модели. Ч. I. построение предельной непрерывной модели»

ПРИБОРЫ И СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

УДК 62-50

DOI: 10.17586/0021-3454-2019-62-5-411-418

РАСЧЕТ ЦИФРОВЫХ СЛЕДЯЩИХ СИСТЕМ КОМБИНИРОВАННОГО УПРАВЛЕНИЯ С ИСПОЛЬЗОВАНИЕМ ПРЕДЕЛЬНОЙ НЕПРЕРЫВНОЙ МОДЕЛИ. Ч. I. ПОСТРОЕНИЕ ПРЕДЕЛЬНОЙ НЕПРЕРЫВНОЙ МОДЕЛИ

А. И. Ермоленко, А. И. Коршунов

Военно-морской политехнический институт ВУНЦВМФ „Военно-морская академия им. Н. Г. Кузнецова",

198514, Санкт-Петербург, Россия E-mail: a.i.korshunov@mail.ru

В общем случае произвольного порядка непрерывной части построена предельная непрерывная модель цифровой следящей системы (ЦСС) комбинированного управления, к свойствам которой неограниченно приближаются свойства ЦСС при повышении частоты измерения рассогласования. Получено аналитическое выражение для коэффициентов дискретного фильтра, вырабатывающего добавку к задающему воздействию, компенсирующую скоростную ошибку и ошибку по ускорению ЦСС. Коэффициенты дискретного фильтра выражены через параметры непрерывной части ЦСС, что позволяет не определять ее дискретную передаточную функцию.

Ключевые слова: цифровая следящая система, предельная непрерывная модель

Введение. Использование комбинированного управления для повышения точности систем автоматического управления предложено Г. В. Щипановым еще в середине прошлого века. В развитие теории комбинированного управления важнейший вклад внесли А. Г. Ивах-ненко, В. С. Кулебакин, А. И. Кухтенко, Б. Н. Петров, Г. М. Уланов. Особенно успешно комбинированное управление применялось при проектировании высокоточных силовых следящих приводов ракетно-артиллерийских систем [1]. Интерес к комбинированному управлению сохраняется и в наше время, см., например, [2—4]. Применение комбинированного управления особенно целесообразно, „когда задающее воздействие вычисляется в ЦВМ" [5, 6], что соответствует рассматриваемым в настоящей статье цифровым следящим системам (ЦСС).

В случае необходимости повысить динамическую точность действующих ЦСС, например ЦСС корабельных АСУ первых поколений, использование комбинированного управления оказывается единственно возможным способом [7, 8]. Однако при низкой частоте измерения рассогласования, характерной для ЦСС АСУ первых поколений, при компенсации скоростной ошибки и ошибки по ускорению не удается получить высокое качество переходных процессов при отработке скачков задающего воздействия, хотя качество отработки ненулевых начальных условий остается высоким [8].

Повышение частоты вычисления рассогласования позволяет, как показано в работе [9], повысить качество отработки скачков задающего воздействия при компенсации скоростной ошибки и ошибки по ускорению. В упомянутой работе предложен способ расчета ЦСС

с использованием ее предельной непрерывной модели; к свойствам этой, полученной в статье модели, свойства ЦСС неограниченно приближаются при бесконечном повышении частоты вычисления рассогласования. Идея использования предельной непрерывной модели систем с периодическим высокочастотным изменением структуры [10] оказалась достаточно плодотворной также для анализа и синтеза этого сложного и пока малоизученного класса нелинейных дискретных систем. В работе [9] получена предельная непрерывная модель конкретной ЦСС, непрерывная часть которой имела третий порядок. При этом были использованы коэффициенты передаточной функции ЦСС, полученные в результате выполнения довольно трудоемкого 2- преобразования.

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

Построение предельной непрерывной модели ЦСС. Передаточная функции линеаризованной ЦСС с астатизмом первого порядка определяется выражением [11]

Ж (г) = 2

(1)

КЯ (Г>) т

где ЖНЧ (р) =-—--передаточная функции непрерывной части ЦСС, Ят (р) = 1 + 2 г^р1,

Р<2п-1( Р )

г=1

п-1

Qn-l (Р) = 1 + 2 ЧгР1 , К — коэффициент преобразования.

г=1

Разложение на простейшие дроби

ЖнЧ (Р) = К Ят (Р)

= К

где

А=1, Б=Г1-дь

А в п-1 С — + - + 2-—

Р2 Р 1=1 Р - Р1

/^0п-1(Р;), Рг

Р Р бп-1(Р) Ят (Рг) / А

ёР

корни уравне-

Рг

ния Qn-1 (Р) = 0, Яе РI < 0, г = 1, 2,..., п -1, и выполнение 2-преобразования дает

И-1 - 1 Яп-1( г)

/ С - = КТ —

Т

Ж (-) = КТ

1 + ^+1 -1'

- -1

V

Т

г=1 - - у

(- - 1)Qn-l (-)

(2)

где ёг = ехр( РгТ), Т — период вычисления рассогласования (шаг квантования по времени);

п-1 п-2 п-1

Qn_1 (-) = ^(- -= -п-1 + V а—, Яп-1 (г) = Vв-, а ь Рг — постоянные коэффициенты.

г=1 г=0 г=0

Очевидно полезное в дальнейшем равенство

Нш(- - 1)Ж(-) = (- - 1)Ж(-) | -=1 = КТ. (3)

Получение в общем виде выражений для коэффициентов числителя передаточной функции Щ-) достаточно трудоемко, сами выражения весьма громоздки и неудобны для получения предельной непрерывной модели. Для исключения этого затруднения можно использовать приближенный способ выполнения 2-преобразования, дающий точный результат при Т ^ 0. Поэтому при построении предельной непрерывной модели приближенный характер метода не является препятствием к его применению. Согласно этому методу, не требующему разложения дроби на простейшие [12], получаем

Ж (-) = --12

- -1 -_К_ Ят (Р)

- Т Р ^п-1( Р)

КЯт (Р )

- -1

Р^-1( Р)

Р=-

--1

(4)

На рис. 1 представлена структурная схема ЦСС комбинированного управления [12]. При компенсации скоростной ошибки и ошибки по ускорению передаточная функция дискретного фильтра, вводящего в закон управления первую и вторую разности задающего воздействия, имеет вид [8, 9]

1 1 2 (5)

ф(z) = Т1(1 - z-1) + т2(1 - z-1)2.

Рис. 1

Передаточная функция ЦСС комбинированного управления по ошибке [12] после подстановки приближенного выражения Ж (2) (см. формулу (4)) имеет вид

Ф0(z) =

0(z) _1 - ф(z)W(z)

R (z)

0BX ( z) 1 + W (z ) z2( pQn -1( p) + KRm (p))

p=-

z-1

(6)

где

R0 (z) = z2pQn-1 (p) - (T1 Tzp + t2T2p2 )KRm (p)

p=-

z -1 = (z2 - zKx{T)p +

+(z2q1 - zr1Kt1T - Kt2T2 )p2 + (z2q2 - zr2Kt1T - Kr1t2T2 )p3 +...

p=

z-1

(7)

Параметры т1, T2 дискретного фильтра ф(z) (5) должны обеспечить в установившемся

режиме безошибочную отработку равноускоренного задающего воздействия

0вх[п7]=0,5б(п7)21[п7], 0вх (z) = Z {0вх [nT ]} = 0,5sT2 z (z +1) / (z -1)3, где 8 = const — ускорение.

Определив Z-преобразование ошибки управления при отработке равноускоренного задающего воздействия, получим

0( z) = Z{0[nT ]} = Ф0 (z)

8T2 z (z + 1)

3

(8)

Если существует установившаяся ошибка при слежении с постоянным ускорением, то ее можно определить, вычислив предел

еГ2

08 = lim(z -1)0(z) = ^ lim(z -1)z(z+1) x

z^i 2 z(z -1)

(z2 - zKt1T)p + (z2q1 - zr1Kt1T - Kt2T2)p2 + (z2q2 - zr2Kt1T - Kr1t2T2)p

z2( pQn-1( p) + KRm ( p))

Очевидно, что выполнение равенств

1 - Kt1T = 0,

q1 - r1Kt1T - Kt2T 2 = 0j

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

p=-

z-1

(9)

и расположение корней уравнения

T

T

T

z -1

Q

z -1

z -1

т Л-1 ^{—) = 0 (П>

внутри центрального круга единичного радиуса на плоскости 2 является необходимым и достаточным условием равенства нулю ошибки по ускорению (08 = 0). Заметим, что расположение корней уравнения

рОп_1(р) + КЯт (р) = 0 (12)

слева от мнимой оси плоскости р — только необходимое условие расположения корней уравнения (11) внутри круга с единичным радиусом.

Решая систему линейных уравнений (10), находим

Т1 = 1/КТ, Т2 = (41 -Г1)/КТ2. (13)

Очевидно, что полученное значение Т1 обеспечивает компенсацию скоростной ошибки. Такое же значение получается при использовании точной передаточной функции ЦСС, имеющей первый порядок астатизма [12]. Значение же Т2 требует уточнения (рассматривается далее).

Для построения предельной непрерывной модели ЦСС комбинированного управления с компенсацией скоростной ошибки и ошибки по ускорению рассмотрим приближенную частотную передаточную функцию по ошибке в функции абсолютной псевдочастоты [5]

. 2 + юТ

л = —-:

Т 2

Ф 0 (jk) = Ф 0 (z),

1+ }УГ/ 2 "1-jkT/ 2

A(jk)= 1 + jk Г/ 2

f 1 + jk Г/ 2 Y

1 - jk Г/ 2 J 1 - jk Г/ 2 ^ n~1 ^ 1 - jk Г/ 2

jk

. A(jk) B(jk)

Qn -1

f jk

jk

1 - jk Г/ 21- jk Г/ 2

+(41 -1)

jk

1 - jk Г/ 2

Rm

jk

1 - jk Г/ 2

B( jk)=

f 1 + jk Г/ 2' 1 - jk Г/ 2 _

jk

Q

1 - jk Г/ 2 {1 - jk Г/ 2

jk

)

KRm

jk

1 - jk Г/ 2

Переходя к пределу при Г ^ 0, с учетом lim k = ю, где ю — круговая частота, получаем

Г ^0

Ф0 0'ю)=

_ jюQn-1 ('ю) + ('ю + (41 - r1)(7ю) )Rm ('ю)

7ю&-1(-/ю) + KRm Ою)

(14)

Нетрудно проверить, что выражение (14) представляет собой частотную передаточную функцию по ошибке непрерывной следящей системы комбинированного управления, структурная схема которой приведена на рис. 2

Рис.2

Следовательно, линеаризованная модель ЦСС комбинированного управления (см. рис. 1) при периоде вычисления рассогласования Т, стремящемся к нулю, приближается по своим свойствам к непрерывной следящей системе (см. рис. 2).

Точное значение коэффициента т2 (13) при конечном периоде Т. Точные значения т1 и т2 можно определить по точной дискретной модели ЦСС комбинированного управления из условия третьего порядка астатизма передаточной функции эквивалентной ЦСС с управлением по отклонению [13]. Полагая наличие единичного корня характеристического полинома этой передаточной функции с кратностью 3, равной порядку астатизма, получаем систему двух линейных уравнений:

0И_1(1) - Т1КТЯп_1(1) = 0, 1

20,-1 (1) + dQn-l / с2(1) - тКТЯП-! / ^(1) - Т2КТЯп-1 (1) = 0,} ()

п-1 п-2

где 0п-1(2X Кп-1(2) — полиномы Ж(г) (см. формулу (2)), 0п-1(1) = П (1 - 4) =1 + Е ai,

1=1 1=0

Кп-1(1) = Е Р.-, ¿0п-\/¿2(1) = п - 1 + , ёЯп-1/¿2(1) = £ ¡в.

.=0 .=1 .=1

Очевидно, что основная трудность определения точных значений т1 и т2 описанным в работе [9] методом заключается в необходимости определения точных значений коэффициентов а. и р.. Другой подход к вычислению т1 и т2 состоит в определении скоростной ошибки и ошибки по ускорению. Как показано в [12], в общем случае астатизма ЦСС первого порядка Т1 = 1/ КТ . В частном случае, согласно [8], значение т2 выражено через К, Т и параметры непрерывной части (постоянные времени).

Для определения ошибки по ускорению в общем случае при скомпенсированной скоростной ошибке требуется вычислить предел

ва = ВШ(, -1) ь£ 2 (2 +1)1 - т1(1 - 2 -)П-) = ет2,;т 1 - (1 --1)У М' КТ =

8 _ ^ ' О 1\3 1 . Т7-7-/ \

2-1 2 (2 - 1)3 1 + Ж (2) 2 - 1)[2 - 1 + (2 - 1)Ж(г)]

= е^ ^ 2 - (2 - т, .0)/КТ = £ ,Ш |1 - 4 [(.. - 1)Ж (2)/КТ ]} =

КТ 2-1 (2 - 1) КТ 2-1 [

ЕТ2 Г С Г - - 1 п-1

81 -НшЬ-■

КТ 2-1 сЬ

1 - а-а (2 -1)+1X е..

Т Т-=1 . 2 - ^

— ^ + ^ |. (16)

КТ I Т 1

При вычислении предела использовано соотношение (3) и правило Лопиталя. Дифференцирование функции (2 - 1)Ж (2) / КТ выполнено с учетом выражения (2).

Для компенсации ошибки по ускорению к 9вх [пТ] необходимо прибавить ее вторую

разность: А^ [пТ] = 0вх [пТ] - 20вх [(п - 1)Т] + 0вх[(п - 2)Т] = 8Т—[п2 - 2(п -1)2 + (п - 2)2 ] = вТ2,

умноженную на коэффициент т2, обеспечивающий равенство добавки и ошибки по ускорению 08:

Т28Т2 =08 =fT(l+qT¿)=KT--(Г + ?,-„). (17)

Из уравнения (17) следует выражение

Т2 = ^ ■ С8)

не содержащее коэффициентов передаточной функции Ж (2) (2).

Таким образом, в случае компенсации обеих ошибок (скоростной и по ускорению)

Ф(,) = -Т<1-,-уТКр0-2-')2. (19)

Заключение. Получена предельная непрерывная модель ЦСС комбинированного управления с компенсацией скоростной ошибки и ошибки по ускорению при любом порядке ее непрерывной части.

Близость свойств ЦСС и ее предельной непрерывной модели позволяет в общем случае использовать для ее расчета любые методы синтеза непрерывных систем, достаточно хорошо разработанные, а удовлетворение ЦСС заданным требованиям обеспечивать посредством выбора высокой частоты вычисления рассогласования.

список литературы

1. Яворский В. Н., Бессонов А. А., Потапов А. М. Проектирование инвариантных следящих приводов. М.: Высш. школа, 1963. 428 с.

2. Лукьянова Г. В., Никифоров В. О. Алгоритм компенсации внешних детерминированных возмущений: операторный метод синтеза // Научно-технический вестник СПбГУ ИТМО. 2003. № 10. С. 5—10.

3. Бобцов А. А., Лукьянова Г. В., Никифоров В. О. Алгоритм компенсации внешнего гармонического возмущения неизвестной частоты для систем активной виброзащиты // Изв. вузов. Приборостроение. 2007. Т. 50, № 11. С. 39—43.

4. Никифоров В. О., Лукьянова Г. В. Следящая система комбинированного управления // Научно-технический вестник СПбГУ ИТМО. 2011. Т. 76, № 6. С. 39—43.

5. Бесекерский В. А. Цифровые автоматические системы. М.: Наука, 1976. 575 с.

6. Бесекерский В. А., Попов Е. П. Теория систем автоматического регулирования. М.: Наука, 1972. 768 с.

7. Ермоленко А. И. Повышение точности цифровых следящих систем при вводе данных в торпедное оружие на маневре уклонения // Материалы XV молодежной науч.-техн. конф. „Взгляд в будущее — 2017". СПб: ЦКБ МТ „Рубин", 2017. С. 239—249.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

8. Ермоленко А. И., Коршунов А. И. Повышение динамической точности цифровых следящих систем АСУТП методом комбинированного управления. Ч. I. Низкий темп вычисления рассогласования // Изв. вузов. Приборостроение. 2018. Т. 61, № 4. С. 309—215.

9. Ермоленко А. И., Коршунов А. И. Повышение динамической точности цифровых следящих систем АСУТП методом комбинированного управления. Ч. II. Высокий темп вычисления рассогласования // Изв. вузов. Приборостроение. 2018. Т. 61, № 12. С. 1036—1044.

10. Коршунов А. И. Предельная непрерывная модель системы с высокочастотным периодическим изменением структуры // Изв. вузов. Приборостроение. 2009. Т. 52, № 9. С. 42—48.

11. Коршунов А. И. Основы теории управления. Руководство к курсовому проектированию. Петродворец: ВМУРЭ им. А. С. Попова, 1998. С. 136.

12. Коршунов А. И. Основы теории управления. Часть II. Основы теории и систем автоматического управления. Петродворец: ВМПИ, 2017. 167 с.

13. Бесекерский В. А., Федоров С. М. Применение эквивалентной передаточной функции при расчете следящих систем комбинированного управления // Тр. I Междунар. конгресса по автоматическому управлению. М.: Изд-во АН СССР, 1961. Т I. С. 154—165.

Сведения об авторах курсант; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", факультет САУ д-р техн. наук, профессор; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", кафедра радиоэлектроники; E-mail: a.i.korshunov@mail.ru

Поступила в редакцию 17.09.18 г.

Ссылка для цитирования: Ермоленко А. И., Коршунов А. И. Расчет цифровых следящих систем комбинированного управления с использованием предельной непрерывной модели. Ч. I. Построение предельной непрерывной модели // Изв. вузов. Приборостроение. 2019. Т. 62, № 5. С. 411—418.

CALCULATION OF DIGITAL TRACKING SYSTEMS OF COMBINED CONTROL USING A LIMITING CONTINUOUS MODEL.

PART I. CONSTRUCTION OF A LIMITING CONTINUOUS MODEL

A. I. Ermolenko, A. I. Korshunov

Naval Polytechnic Institute "N. G. Kuznetsov Naval Academy", 198514, St, Petersburg, Russia E-mail: a.i.korshunov@mail.ru

In the general case of continuous part of an arbitrary order, a limiting continuous model of digital tracking system (DSS) of combined control is constructed, so that the DSS properties approach indefinitely to the properties of the model with increasing frequency of the mismatch measurements. A discrete filter is considered which generates an additive to the master control, compensating for the speed error and the error in the DSS acceleration by variation of parameters of the continuous part of the system. An analytical expression for the filter coefficients is derived.

Keywords: digital tracking system, limiting continuous model

REFERENCES

1. Yavorskiy V.N., Bessonov F.F., Potapov A.M. Proektirovanie invariantnykh sledyashchikh privodov (Design of Invariant Tracking Drives), Moscow, 1963, 428 p. (in Russ.)

2. Luk'yanova G.V., Nikiforov V.O. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2003, no. 10, pp. 5-10. (in Russ.)

3. Bobtsov A.A., Luk'yanova G.V., Nikiforov V.O. Journal of Instrument Engineering, 2007, no. 11(50), pp. 39-43. (in Russ.)

4. Nikiforov V.O., Luk'yanova G.V. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2011, no. 6(76), pp. 39-43. (in Russ.)

5. Besekerskiy V.A. Tsifrovye avtomaticheskie sistemy (Digital Automatic Systems), 1976, 576 p. (in Russ.)

6. Besekerskiy V.A., Popov E.P. Teoriya sistem avtomaticheskogo regulirovaniya (The Theory of Automatic Control Systems), Moscow, 1972, 768 p. (in Russ.)

7. Ermolenko A.I. Materialy XV molodezhnoy nauchno-tekhnicheskoy konferentsii "Vzglyad v budushchee-2017" (Proceedings of the XV Scientific and Technical Youth Conference "Looking to the Future-2017"), St. Petersburg, 2017, рр. 239-249. (in Russ.)

8. Ermolenko A.I., Korshunov A.I. Journal of Instrument Engineering, 2018, no. 4(61), pp. 309-315. (in Russ.)

9. Ermolenko A.I., Korshunov A.I. Journal of Instrument Engineering, 2018, no. 12(61), pp. 1036-1044. (in Russ.)

10. Korshunov A.I. Journal of Instrument Engineering, 2009, no. 9(52), pp. 42-48. (in Russ.)

11. Korshunov A.I. Osnovy teorii upravleniya. Rukovodstvo k kursovomu proektirovaniyu (Bases of the Theory of Management. The Management to Course Design), Petrodvorets, 1998, 136 р. (in Russ.)

12. Korshunov A.I. Osnovy teorii upravleniya. Osnovy teorii i sistem avtomaticheskogo upravleniya (Bases of the Theory of Management. Bases of the Theory and Systems of Automatic Control), Petrodvorets, 2017, 167 р. (in Russ.)

13. Besekerskiy V.A., Fedorov S.M. Trudy I mezhdunarodnogo kongressa Mezhdunarodnoy federatsii po avtomaticheskomu upravleniyu (Proceedings of the I International ingress of the International Federation on Automatic Control), Moscow, 1961, vol. I, pp. 154-165. (in Russ.)

Артем Игоревич Ермоленко Анатолий Иванович Коршунов

Data on authors

Artem I. Ermolenko — Cadet, Naval Polytechnic Institute "N. G. Kuznetsov Naval Academy",

Department of Automatic Control Systems Anatoly I. Korshunov — Dr. Sci., Professor; Naval Polytechnic Institute "N. G. Kuznetsov Naval

Academy", Department of Radio Electronics; E-mail: a.i.korshunov@mail.ru

For citation: Ermolenko A. I., Korshunov A. I. Calculation of digital tracking systems of combined control using a limiting continuous model. Part I. Construction of a limiting continuous model. Journal of Instrument Engineering. 2019. Vol. 62, N 5. P. 411—418 (in Russian).

DOI: 10.17586/0021-3454-2019-62-5-411-418

i Надоели баннеры? Вы всегда можете отключить рекламу.