CC BY
75
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНЦА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Метод гармонической линеаризации в задаче идентификации нелинейных динамических систем # 04, апрель 2014
DOI: 10.7463/0414.0704613
Павлов Ю. Н., Недашковский В. М., Тихомирова Е. А., Шавырин И. Б.
УДК 01.04.01
Россия, МГТУ им. Баумана pavlovrtbmstu.nir vmniu.3 ;'5:bmstu.rur elizarti^bmsturu*" shaxirin ffbrn sturu*"
Введение
Динамические системы (в основном это различного вида силовые приводы) приобретают в процессе изготовления нежелательные дефекты, которые существенно влияют на качество и точность работы объекта в целом. Наличие в кинематической схеме существенных нелинейностей типа зон нечувствительности, люфтов, сухого трения существенно ухудшают ситуацию и являются причиной появления в системе автоколебаний, что, в конечном счете, приводит к ухудшению точности и быстродействия работы объекта управления. Динамика подобных систем описывается нелинейными дифференциальными уравнениями. Проблема решения нелинейных дифференциальных уравнений при заданных коэффициентах и виде присутствующей нелинейности в настоящее время успешно решается с применением различного рода линеаризации [1-3]. Были также разработаны различные методы решения обратной задачи, т.е. определения коэффициентов нелинейных дифференциальных уравнений с использованием экспериментальных данных [4]. В данной статье предложен метод решения обратной задачи - по результатам экспериментальных данных, полученных путем подачи на систему тестовых сигналов, определяются коэффициенты дифференциального уравнения и параметры нелинейного звена.
1.Постановка задачи
Предполагается, что с реальной динамической системы, описываемой нелинейным дифференциальным уравнением второго порядка, путем подачи на вход гармонических сигналов получены амплитудно-частотная и фазо-частотная характеристики (эксперимен-
тальные данные). Необходимо создать алгоритм и программу обработки экспериментальных данных с целью определения коэффициентов дифференциального уравнения и путём статистического моделирования оценить величину погрешностей определения коэффициентов и параметров нелинейности.
2. Методика решения задачи
Предлагаемая методика решения задачи состоит из следующих этапов:
1) Выбор вида существенной нелинейности и определение коэффициента гармонической линеаризации.
2) Математическое моделирование процесса получения псевдоэкспериментальных данных.
3) Создание алгоритма, проведение вычислительного эксперимента, статистическая оценка ошибок определения коэффициентов передаточной функции динамической системы и параметров нелинейного звена.
В статье рассматривается система, описываемая дифференциальным уравнением второго порядка с нелинейным звеном типа «сухое трение». При подаче на вход системы синусоидального сигнала q sin ci в системе установятся незатухающие колебания с амплитудой A i и частотой c i, что и позволяет для решения поставленной задачи применить
метод гармонической линеаризации.
Дифференциальное уравнение второго порядка, описывающее вынужденное движение динамической системы с нелинейным звеном типа «сухое трение», имеет вид
e2x(t)+Fsigrn(t) +e0x(t) = qsinct (1)
График нелинейности типа «сухое трение» приведен на рис. 1
F i к
i i к F г
0 X
Рис.1 График нелинейности типа «сухое трение»
Неизвестные коэффициенты е0, е2, Е дифференциального уравнения (1) должны быть определены в результате решения задачи идентификации.
В статье предложен алгоритм получения оценок неизвестных коэффициентов e0, e2, F дифференциального уравнения (1) по полученному экспериментальному годографу.
В описываемом алгоритме определения неизвестных коэффициентов e0, e2, F предлагается воспользоваться методом гармонической линеаризации и нелинейность типа "сухое трение" аппроксимировать вязким трением с соответствующим коэффициентом гармонической линеаризации [5].
Вынужденное движение с использованием метода гармонической линеаризации описывается уравнением
e2x(t) + e (C)x(t) + e0x(t)—qsin ct, (2)
где коэффициент гармонической линеаризации [5]
/ \ 4F
eiС)- л( \ , (3)
A(c) - амплитуда синусоидальной составляющей выхода x(t), имеющей частоту с.
Введем обозначение
ж
Тогда
eС)-АСС (4>
Из уравнения (2) следует, что при использовании метода гармонической линеаризации частотная передаточная функция динамической системы второго порядка с сухим трением имеет вид
W (jc) --гтт1--—у .
e0 + ei [c)jc + e2(jc)
Частотную передаточную функцию W(jC) можно также записать в виде [6]
W(jc) - P(c) +jQ(c>). (5)
Здесь P(c) и QC) - вещественная и мнимая части частотной передаточной функции соответственно, которые задаются соотношениями
P(c) --, (6)
г 2\2 С
(e0 - e2c ) +
А (с)
c
О(а) =--А(о) 2 . (7)
/ 2 \ 2 С
(ео -е2а )
А (а)
Тогда квадрат значения амплитудно-частотной характеристики динамического звена для частоты а можно определить по формуле
А2(а)=Р2(а)+0(0), (8)
или с учетом (6) и (7)
Л2 (а) = ■
(e0 - e2^2)2 +
Л2 (а)
Разрешив это уравнение относительно А(а) , получим
Л2 (а) =
1 - c2
(e0 - e2^2)2
Л(а) =
л/Т—
(9) (10)
Видим, что амплитуда вынужденных колебаний имеет действительное значение лишь при условии 1 — е2 > 0.
Пример амплитудно-частотной характеристики системы приведен на рис.2
25
20
15
10
0 0,5 i 1 1,5 2 2,5 3
™(Гц)
Рис. 2. Пример амплитудно-частотной характеристики нелинейного динамической системы второго порядка
с сухим трением
Обратим внимание, что при e0 — e2 а2 = 0 или при а = сог =
V
e0
— значение Л(а) = да, e
т.е. амплитудно-частотная характеристика системы второго порядка, несмотря на наличие трения на резонансной частоте имеет разрыв.
Соотношения (9) и (10) позволяют переписать выражения (6) и (7) в виде
Р(а) = ■
1 — c2
(11)
Q(c) = —т
-л/Т—
С учетом (11) найдем выражение для значений (р(а>) фазо-частотной характеристики системы [6]:
2
c
2
2
e0 — e 2а
2
e0 — e2®
2
2
e0 — e2®
% ((а)) =
% ((а))=
0(а)_
Р(а) VI— ОС®) с
Р(а) ^Гс
, если е0 — е2 а2 >0,
, если е0 — е2а2 <0 .
С учетом (12) отсюда имеем (а) = —ат^^
л/Т—
-) + кп , к = ±1, ±2,... , если е0 — еа? >0
(13)
((а) = атс^^
-) + кп , к = ±1, ± 2,... , если е0 — е2а2 <0 .
Если учесть, что при е0 — е2а2 >0 Р(а) >0, 0(а) <0, а при е0 — е2а2 <0 Р(а) <0, ОаО) <0 , то соотношения (13) можно переписать в виде
((а) = —атС%(
л/Т—
г) , если е0 — е2 а2 >0
(14)
((а) = атсг%(
л/Т—
— п , если е0 — е2а2 <0 .
Пример фазо-частотной характеристики системы с учетом (14) приведен на рис.3
Рис. 3. Пример фазо-частотной характеристики нелинейного динамической системы второго порядка с
сухим трением
Из рис. 3 видно, что фазо-частотная характеристика имеет постоянное значение, определяемое величиной сухого трения и изменяет его скачком на резонансной частоте
ат=1
Частотная передаточная функция Ш(]а>) может быть изображена на комплексной плоскости в виде годографа [6].
2
2
2
с
2
с
2
с
2
Пример годографа системы приведен на рис.4
Рис. 4. Пример годографа нелинейной системы второго порядка с сухим трением
Видим, что годографом в этом случае являются две прямые линии: одна линия, выходящая из точки, определяемой амплитудой вынужденных колебаний при о=0 и уходящая в бесконечность, другая линия, возвращающаяся из бесконечности в ноль. Угол наклона этих прямых определяется величиной сухого трения.
3. Математическое моделирование с целью получения псевдоэкспериментальных данных.
При получении частотных характеристик реальных систем в эксперименты вмешиваются помехи, в результате которых точки годографа смещаются случайным образом. Введем обозначения для определенных с погрешностями вещественных и мнимых значений отсчетов полученных с учетом соотношений (11) для значений псевдоэкспериментального годографа Щ (]0) динамической системы для пехр значений частот
р = Р(о).., Г„хР = Р(апекр), 0 = 0(0!),..., 0пекр = 0(апекр). (15)
На рис. 5 показан пример искаженного псевдоэкспериментального годографа нелинейного звена второго порядка с сухим трением относительно годографа, приведенного на рис. 4 при наличии случайных погрешностей в псевдоэкспериментальных данных для пехр=10.
Рис. 5. Пример искаженного годографа нелинейной системы второго порядка с сухим трением при наличии случайных погрешностей в псевдоэкспериментальных данных
Решение задачи идентификации для нелинейной системы второго порядка с сухим трением будем искать в классе годографов, задаваемых моделью
1
Щм (За) =
Здесь с учетом (3) и (4)
а0 + а (а)]о + а2 (]а)7 40
(16)
а
1 (а) =
Тогда
Щм (а) =
лА(а)а 1
.40 2
а0 + 3^Т\ + а2(3а) лА(ю)
(17)
или
Щ (За) =
где а = 1, ( = 0, у = а0 - а2а>2, 8 = ■ 4
а + ЗР Г + 38
(18)
(19)
лА(а)
Отклонение А Щ г-ого отсчета Щм (Заi) годографа модели на частоте а; от г отсчета Щэ (Заi) экспериментального годографа идентифицируемой системы равна
ащ=щ (За) Щ(зщ).
С учетом (5) и (18) можно записать в виде
ащ = р + 30i -а +]Рг
(20)
Г г + 3 8г
где параметры а, Р, у1, 8 соответствуют частоте а г.
Необходимо ввести приемлемый критерий, характеризующий близость двух годографов на всей совокупности экспериментальных точек, и минимизировать его по пара
метрам модели системы а0, а2, С . В качестве критерия (меры) близости можно было бы выбрать сумму квадратов модулей расхождений а wi:
п ехр
I = ^\АЩ2 (21)
1=1
Минимизация меры I приводит к нелинейной системе уравнений для определения коэффициентов а0, а2 О модели. В работах [7, 8], применен способ, которым мы воспользуемся, и который состоит в следующем. Соотношение (20) умножим на отличный от нуля комплексный множитель (/ + }8{) :
Щ =Щ(у1+у$). (22)
Тогда с учетом (20) и (22) для а н1 и для |анг. |2 получим
щ = рп — —а +лрд а —р) (23)
|АН|2 = (Р/ — в а —а, )2 + (РА + д,г, —р )2.
В качестве меры близости годографов вместо меры I (21) примем меру J, равную
сумме квадратов модулей |анг |2 (23):
п ехр п ехр
з = £|ан, I2 = £[(р7г — вьА —а )2 + (рА + а/ —Д )2]. (24)
1=1 1=1
Мера близости J экспериментально полученного годографа системы и годографа модели системы на всей совокупности экспериментальных точек задается формулой (24), которая после некоторого преобразования принимает вид
п ехр
3 = £[(р-2 + +(р2 + е2)А2 — 2Р,га, + 20Аа, — 2РАД — 20,/,?, +а2 + Д2] . (25)
1=1
Мера J (25) является функцией параметров а0, а2, С частотной передаточной функции модели (16) или (17). Для минимизации меры J приравняем нулю частные производные от J по этим параметрам:
= 0, 3 = 0, 3 = 0. (26)
5а0 да2 до
Для рассматриваемой нелинейной системы второго порядка система уравнений (26) с учетом (19) принимает вид
пехр п ехр п ехр
,е(р2+а2)—а2 е(р2+а 2а2 =ер
,=1 ,=1 ,=1
пехр п ехр пехр - -.2 ^>2Ч 2 , тл2 ^>2Ч 4 ^ 2
ао Е (р2 + а 2)а2 — а2 Е (р2 + а 2)а4 =Е ра2, (27)
,=1 ,=1 ,=1 ■>2 .
пехр 1 п пехр Ц
с Е(р2 + а2)-^=—ПЕа
,=1
л2 (а) 4 ^ г А(а)
Первые два уравнения системы уравнений (27) приведем к матричному виду
Ш = и , (28)
где
¥11 ¥12 — ^0 0 — щ
, а = , и =
¥21 ¥22_ а 2 _ и2 _
¥ =
щ= Т(р2 + а2),
¡=1 пехр
%г=Т(р2 + а >,
¡=1 п ехр
щ = £ р,
пехр
п ехр
^12=-Е (р,2+о 2>,2
¡=1 пехр
%2=-Х(Р2 + а2)а4,
¡=1
п ехр
Х>2 .
Щ2 =
Решая матричное уравнение (28), получим
^ _ Щ22Щ1 -Щ12Щ2
т11т22 т12т21
1и2 Щ2 Щ1 т11т22 т12т21
(30)
С учетом выражения (8) третье уравнение системы уравнений (27) принимает вид
О = -
п
п ехр 2>
0
(31)
4п ехр р2 + 02 '
Покажем, что если отсчеты годографа Р, и О не содержат ошибок измерения, т.е. с учетом (11) и (15)
1 - с2
Р =-
е0 - е2®г
(32)
= -
;л/Т-
то значения параметров а0, а2, О, вычисленные по соотношениям (30) и (31) совпадают с параметрами е^ е2, ¥ дифференциального уравнения (1), описывающего вынужденное движение нелинейного динамической системы второго порядка с сухим трением.
Действительно, подставив в соотношения (30) соответствующие значения из (29) и (32), получим
^22^1 -¥12»2 = е0(1 - С ) ^^Т-^-^
,=1 -=1 (е0 - е2а, ) (е0 - е2а- )
2 2п®1>п1еЕР а2 (а2 -а2)
¥цЩ -¥2Уи1 = е2(1 - с ) ^^
,=1 2 =1 (е0 - е2а, ) (е0 - е2а2 )
п^,™ а (а2 -а) ¥1 ¥22-¥12¥21 = (1 - с ) xx"
,=1 -=1 (е0 - е2аг ) (е0 - е2а, )
а2 =
Щ Щ -Щ Щ
Щ1 1Щ2 - Щ2 Щ1
= еп
= е2 .
Щ Щ -Щ Щ
т11т22 т12т21 С учетом (32) соотношение (31) приводится к виду
1
2
2
е0 -е2аг
а0 =
G = —
n
n exp
z
Q
4n exP P2 + Q
n , 4F
-(--n exp) = F .
4nexp n
Были проведены иллюстративные вычислительные эксперименты по оценке погрешности, с которой предложенный алгоритм определяет значения параметров частотной передаточной функции нелинейной системы второго порядка с сухим трением. В экспериментах заданы конкретные значения параметров е0 = 1, е2 =1, F=0,5.
Погрешности измерения значений р, Q моделировались с помощью генератора
случайных чисел с равномерным законом плотности распределения вероятностей в различных диапазонах [-0,01, 0,01], [-0,05, 0,05], [-0,1, 0,1]. Количество отсчетов в годографе системы (т.е. количество экспериментов) nexp =10, 20, 40, 80. Для каждого конкретного количества экспериментов nexp проводилось nseriy серий этих экспериментов nseriy=25, 50, 100, 200. Диапазон частот [0, 2юг, сог =2Гц], [0, 3 сог, сог =3Гц]. В сериях экспериментов вычислялись погрешности etreQ = е0 — a, erre2 = е2 — а2, errF = F—G определения параметров e, e2, F и среднеквадратические отклонения sko для случайных величин erre0, erre2, errF, а также строились гистограммы для этих величин.
Рис. 6. Зависимость среднеквадратического отклонения sko погрешности erre0 определения коэффициента
е 0 от количества отсчетов nexp годографа системы при погрешности измерений в диапазоне [-0,10, 0,10] и
числе серий nseriy=25, 50, 100, 200
Из рис. 6 видно, что в рассматриваемом случае можно ограничиться количеством отсчетов годографа от 20 до 40.
На рис. 7 в качестве примера показана гистограмма погрешности егге0 определения
коэффициента е0.
Рис. 7. Гистограмма погрешности егге0 определения коэффициента е 0 (пехр=20, диапазон погрешности измерений [-0,05, 0,05], количество серий тепу=200, среднеквадратическое отклонение sko=0.012\)
Вычислительный эксперимент показал, что ошибки имеют нормальный закон распределения и среднеквадратическое отклонение равно sko=0,0121. Это значит, что 68% ошибок находятся в диапазоне ± зко = 0,0121, а 99% ошибок находятся в диапазоне ± Ъзко = ±0,0363. Таким образом, погрешность определения значения коэффициента е 0
=1 в иллюстративном вычислительном эксперименте не превышает диапазона погрешности измерений отсчетов годографа [-0,05, 0,05].
Аналогичные результаты были получены для погрешностей определения коэффициентов е2 и Е . Например, для погрешности егге2 определения коэффициента е2 (пехр=20, погрешность измерений в диапазоне [-0,05, 0,05], количество серий тепу=100), среднеквадратическое отклонение sko=0,0117. Для погрешности еггЕ определения коэффициента Е . и (пехр=20, погрешность измерений в диапазоне [-0,05, 0,05], количество серий тепу=100), среднеквадратическое отклонение sko=0,0l37.
На рис. 8 приведены три годографа: годограф системы, когда погрешность измерения отсчетов годографа равна нулю, годограф системы, когда погрешность измерения отсчетов годографа не равна нулю, годограф найденной модели системы
Рис. 8. Годограф системы при нулевой ошибке измерения отсчетов, годограф системы при ненулевых ошибках измерения отсчетов в диапазоне [-0,10, 0,10]; годограф найденной модели системы (количество
экспериментов пехр=20)
Выводы
В данной статье предложен алгоритм идентификации нелинейной динамической системы по экспериментально полученным отсчетам частотного годографа системы. В качестве меры близости годографов системы и модели системы выбрана сумма квадратов модулей модификации отклонений отсчетов годографа модели системы от соответствующих экспериментально полученных отсчетов системы.
В иллюстративном вычислительном эксперименте показано, что для нелинейной системы второго порядка с сухим трением можно ограничиться количеством отсчетов годографа от 20 до 40.
Иллюстративный вычислительный эксперимент показал, что при использовании предложенного алгоритма идентификации нелинейной динамической системы по частотному годографу погрешность определения значений коэффициентов частотной передаточной функции системы второго порядка с сухим трением сравнима с диапазоном погрешности измерений экспериментальных отсчетов годографа этой системы.
Список литературы
1. Полянин А.Д. , Зайцев В.Ф. , Журов А.И. Методы решения нелинейных уравнений математической физики и механики. М.: Физматлит, 2005. 256 с.
2. Тимошенко С.П. Колебания в инженерном деле: пер. с англ. М.: Наука, 1967. 444 с.
rd
[Timoshenko S. Vibration Problems in Engineering. 3 ed. D. Van Nostrand Company, Inc.Toronto New York London, 1955. 468 p.].
3. Гласко В.Б. Обратные задачи математической физики. М.: МГУ, 1984. 111 с.
4. Дейч А.И. Методы идентификации динамических объектов. М.: Энергия, 1979. 240 с.
5. Попов Е.П., Пальтов И.П. Приближенные методы исследования нелинейных автоматических систем. М.: ГИФМЛ, 1960. 790 с.
6. Основы автоматического управления / под ред. В.С. Пугачева. М.: Наука, Главная редакция физико-математической литературы, 1968. 680 с.
7. Боевкин В.И., Павлов Ю.Н. Регрессионный анализ в прикладной задаче идентификации. М.: Изд-во МГТУ им. Н.Э. Баумана, 1990. (Труды МГТУ им. Н. Э. Баумана; № 546).
8. Боевкин В.И., Недашковский В.М., Павлов Ю.Н. Идентификация линейных динамических звеньев по частотному годографу // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013 . № 9. Режим доступа: http://technomag.bmstu.ru/doc/618917.html (дата обращения 01.03.2014). DOI: 10.7463/0913.0618917
SCIENTIFIC PERIODICAL OF THH BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. N»0421200025. ISSN 1994-0408
electronic scientific and technical journal
Harmonic linearization method in the identification of nonlinear dynamical systems # 04, April 2014 DOI: 10.7463/0414.0704613
Y.N. Pavlov, V.M. Nedashkovskii, E.A.Tihomirova, I. B. Shavirin
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
pavlov?? bmstu.ru*" vmniu3 ffbm stu .rur elizarti^bmstu.rur shavirin fgjbmstu .ruc
The subject of this work is a problem of identification of nonlinear dynamic systems based on the experimental data obtained by applying the test signals to the system.
The goal is to examine the opportunity for using the obtained experimental samples of the hodograph to define the parameters of this system by the specific example of a second-order dynamical system with a type of nonlinearity "dry friction".
It was suggested to use the harmonic linearization method and to approximate the nonlinearity of "dry friction" by linear friction with the corresponding harmonic linearization coefficient.
It was assumed that the frequency transfer function of the identified system is well-known
value.
It was assumed also that there are disturbances while obtaining frequency characteristics of the real-world system. As a result of that, the points of experimentally obtained hodograph move randomly.
Searching for the solution of the identification problem is accomplished in the hodograph class, specified by the system model, which has the form of the frequency transfer function the same as the form of the frequency transfer function of the system identified.
Searching for the unknown coefficients of the frequency transfer function model of the system was carried out by minimizing the previously described and published by one of the authors proximity criterion (measure) of experimentally obtained system hodograph and the hodograph of system model for all the experimental points.
The problem solution to identify the nonlinear dynamic system identification in frequency hodograph was reduced to solving an equations system, linear relative to the unknown parameters of the frequency transfer function of the system model.
For the second-order dynamical system with nonlinearity of the type "dry friction" the simulation software was developed to provide the pseudo experimental data containing random accuracy of and determine the parameters of the system.
A computational experiment of accuracy evaluation was made with which the proposed algorithm determines the values of the system parameters.
The illustrative numerical simulation has demonstrated that the accuracy of determining the values of the coefficients transfer function is comparable with the range of measurement accuracy of experimental samples of this system hodograph.
Иллюстративный вычислительный эксперимент показал, что погрешность определения значений параметров системы сравнима с диапазоном погрешности измерений экспериментальных отсчетов годографа этой системы.
This method of the identification of nonlinear dynamic systems is not mentioned in the well-known publications.
The method of identification of nonlinear dynamical systems, which is described in the article, can be used to determine the parameters of various kinds of actuators.
Using the method of harmonic linearization and identification of dynamical systems hodo-graphs is promising for solving the problem of identification of nonlinear systems with different types of nonlinearities.
Publications with keywords: identification, harmonic linearization, nonlinear dynamical system, dry friction, frequency locus
Publications with words: identification, harmonic linearization, nonlinear dynamical system, dry friction, frequency locus
References
1. Polyanin A.D., Zaytsev V.F., Zhurov A.I. Metody resheniya nelineynykh uravneniy matematicheskoy fiziki i mekhaniki [Methods for solving nonlinear equations of mathematical physics and mechanics]. Moscow, Fizmatlit, 2005. 256 p. (in Russian).
rd
2. Timoshenko S. Vibration Problems in Engineering. 3 ed. D. Van Nostrand Company, Inc.Toronto New York London, 1955. 468 p. (Russ. ed.: Timoshenko S.P. Kolebaniya v inzhenernom dele. Moscow, Nauka, 1967. 444 p.).
3. Glasko V.B. Obratnye zadachi matematicheskoy fiziki [Inverse problems of mathematical physics]. Moscow, MGU Publ., 1984.111 p. (in Russian).
4. Deych A.I. Metody identifikatsii dinamicheskikh ob"ektov [Methods of identification of dynamic objects]. Moscow, Energiya, 1979. 240 p. (in Russian).
5. Popov E.P., Pal'tov I.P. Priblizhennye metody issledovaniya nelineynykh avtomaticheskikh system [Approximate methods of study of nonlinear automatic systems]. Moscow, GIFML Publ., 1960. 790 p. (in Russian).
6. Pugachev V.S., ed. Osnovy avtomaticheskogo upravleniya [Basics of automatic control]. Moscow, Nauka, 1968. 680 p. (in Russian).
7. Boevkin V.I., Pavlov Yu.N. Regressionnyy analiz vprikladnoy zadache identifikatsii [Regression analysis in applied problem of identification]. Moscow, Bauman MSTU Publ., 1990. (Trudy MGTU im. N. E. Baumana [Proceedings of the Bauman MSTU]; no. 546). (in Russian).
8. Boevkin V.I., Nedashkovskiy V.M., Pavlov Yu.N. [Identification of linear dynamic elements using a frequency locus]. Nauka i obrazovanie MGTU im. N.E. Baumana - Science and Education of the Bauman MSTU,, 2013, no. 9. DOI: 10.7463/0913.0618917 (in Russian).
