Чувствительность характеристик автоколебаний в системе с трехпозиционным релейным элементом Текст научной статьи по специальности «Математика»

УДК 681.511.4

ЧУВСТВИТЕЛЬНОСТЬ ХАРАКТЕРИСТИК АВТОКОЛЕБАНИЙ В СИСТЕМЕ С ТРЕХПОЗИЦИОННЫМ РЕЛЕЙНЫМ ЭЛЕМЕНТОМ

Н.В. Фалдин, А.В. Моржов

Рассматриваются системы с трехпозиционным релейным элементом и нелинейным объектом управления. Разработаны методы определения функций чувствительности автоколебаний (периодической траектории, частоты, критерия устойчивости).

Ключевые слова: трехпозиционный релейный элемент, автоколебания, чувствительность, периодическая траектория, частота, критерий устойчивости.

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

При анализе и синтезе релейных систем важно располагать информацией о чувствительности системы к изменению параметров объекта управления. Зная функцию чувствительности, легко определить, как влияют отклонения параметров от номинальных значений на качественные характеристики системы. На практике такие отклонения всегда имеют место. С помощью функций чувствительности, таким образом, можно оценить работоспособность системы в реальных условиях эксплуатации. Они также позволяют выполнить синтез при задании ограничений на чувствительность системы к изменению параметров объекта управления.

В [1-4] разработаны методы исследования чувствительности систем с двухпозиционным релейным элементом. Однако аналогичные методы необходимо иметь и для систем с трехпозиционным релейным элементом. Настоящая статья посвящена разработке методов исследования чувствительности автоколебаний в системах с трехпозиционным релейным элементом и нелинейным объектом управления. Рассматривается чувствительность следующих характеристик автоколебаний: периодической траектории; критерия устойчивости.

Чувствительность периодического движения. Рассмотрим релейную автоколебательную систему с трехпозиционным релейным элементом (рис.1).

При отсутствии входного сигнала (y(t) ° 0) движение системы задается уравнениями

dx

—=f(x,au), (1)

dt

и=Ф(-ЯТ х,Х,Ъ), (2)

где х=(х\, Х2,..., хп) и f=(/1, /2,..., /п) - п-мерные векторы.

Функция Ф задается статической характеристикой трехпозицион-ного релейного элемента (рис. 2).

Рис. 1. Релейная автоколебательная система с трехпозиционным релейным элементом: у(?) - входной сигнал, ЯТ - матрица-строка

Рис. 2. Статическая характеристика трехпозиционного

релейного элемента

В равенстве (1) а - некоторый скалярный параметр. Будем предполагать, как это обычно имеет место в следящих системах, что объект управления обладает нечетной симметрией:

{(-х,а,-и)=—Г (х,а,и),

причем симметрия сохраняется при любом значении параметра а.

Фазовый годограф [5] позволяет легко определить возникающие в представленной на рис. 1 и 2 системе (у^) ° 0) автоколебания. Пусть х^ )-периодическая траектория системы (1), (2) периода 2 Т. Объект управления и функция Ф обладают нечетной симметрией. В такой системе периодическая траектория также имеет симметрию, т.е. х^+Т)=-х(?). На рис. 3 изображен вид на периоде сигнала с выхода релейного элемента и^).

477

Рис. 3. Вид на периоде сигнала с выхода релейного элемента и(1)

Малое изменение параметра а приведет к малому изменению периодической траектории. Периодическую траекторию в параметрически возмущенной системе (она также обладает симметрией) обозначим ~(г)=х(г)+дх^), здесь дх^)=(дх^ ),дх 2(1 ),...,5х п (7)).

На траекториях х(?) и ~ (?) момент ?=0 совместим с моментом переключения релейного элемента с нуля на плюс. На траектории х^) управление переключается в точках уТ и Т, а на траектории ~(?) - в моменты уТ+д(уТ) и Т+дТ, где д(уТ) и дТ - малые величины.

Для системы с номинальным значением параметра а, как следует из рис. 1 и 2, имеют место равенства

-ЯТх(уТ)=1Ь, - ЯТх(Т) = -Ь, (3)

а для параметрически возмущенной системы - равенства

- ЯТ ~ (уТ+д(уТ))=1Ь, - ЯТ (~ (Т+5Т))=-Ь. (4)

Представим уравнения (4) в виде

- ЯТ [х(уТ+д(уТ))+дх( уТ+д(уТ))]=1Ь,

- ЯТ [(х(Т+дТ)+дх(Т+дТ)]=--Ь.

Из (5) и (3), принимая во внимание, что д(уТ), дТ и дх^) являются малыми величинами, найдем

ЯТ [х - (уТ)д(уТ)+дх - (уТ)]=0, т (6)

ЯТ [(х (Т )дТ+дх - (Т)]=0.

В равенствах (6) символом «минус» обозначены пределы слева. Далее, в них опущены величины, имеющие порядок малости выше первого. Из (6) следует

5(уТ) = - , (7)

Я7 х- (уТ)

дТ = - . (8)

х-Т

Запишем уравнение (1) для параметрически возмущенной траектории, полагая и^) фиксированной величиной

а №) + 6х(г)] = f (х + бх, а + ба, и). (9)

ш

В качестве фиксированных (см.рис.3) рассматриваются значения и(?) = А, и(1,) = -А, и(?) = 0. Будем, далее, полагать, что функция f непрерывно дифференцируема по х и а и непрерывна по и. Из (9), опуская величины, имеющие порядок малости выше первого относительно бх и ба , найдем

dбx(t) Эf(х,а,и)„ Эf(х,а,и)„ ..

-—=—4 убх+—4 уба. (10)

& Эх Эа

Равенство (10) связывает между собой отклонение траектории движения х^) с отклонением параметра а и называется уравнением в вариациях. Оно представляет собой неоднородное линейное дифференциальное уравнение с переменными во времени коэффициентами.

Траектория х(?) непрерывна в момент уТ + б(уТ). Поэтому можно записать

х - (уТ + б( уТ)) + бх - (уТ + б(уТ)) = х+(уТ + б( уТ)) + бх + (уТ + б( уТ)). (11) Учитывая непрерывность траектории х^) и опуская в (11) величины, имеющие порядок малости выше первого, получим

бх + (уТ) = бх - (уТ) + [х - (уТ) - х + (уТ )]б(уТ). (12)

Здесь и в дальнейшем символом «минус» обозначаются пределы слева, а символом «плюс» - пределы справа.

На траектории х^) управление и(^ принимает значения и=А, и=0,

и = -А. Обозначим У^,to) нормированную фундаментальную матрицу решений уравнения

dx(t) = Эf (х, а, и) (13)

Л Эх

при и = А, У2(^) при и = 0, здесь 10 - начальный момент времени.

Пусть г1^) является решением уравнения (10) при бх(0) = 0, и = А и ба = 1. Тогда

бх - (уТ)=У1(уТ ,0)бх(0)+г1(уТ )ба. (14)

В соответствии с (7), (12) и (14)

бх + (уТ) = Рбх(0) + Оба, (15)

где

(х - (уТ) - х + (уТ ))ЯТ У1( уТ ,0)

Р = У1( уТ ,0)-

Я Т х - (уТ)

О=г1(уТ )-

(х-(уТ)-х + (уТ))ЯТ гА(уТ)

х - (уТ)

Ниже в выражениях, в которых в скобках, наряду с текущим или конечным временем, указывается начальное время, оно записывается после запятой, например, У^, to).

Обозначим г 2(^ уТ) (t >уТ) решение уравнения (10) при бх(уТ)=0, и=0, ба=1. Тогда

бх - (Т) = У2 (Т, уТ )бх + (уТ) + г 2 (Т, уТ )ба.

Принимая во внимание (15), получим

бх- (Т) = М 0бх(0) + N 0ба.

(16)

Здесь

М 0 = У2(Т, уТ )Р, N 0 = У2(Т, уТ )О + г 2(Т, уТ).

Воспользуемся симметрией параметрически возмущенной траектории х^) = х^) + бх^):

- х- (Т + бТ) - бх- (Т + бТ) = х(0) + бх(0). (17)

Опуская величины, имеющие порядок малости выше первого и принимая во внимание непрерывность траектории х^), из (17) найдем

-бх- (Т) - х- (Т )бТ = бх(0).

Учитывая (16), получим равенство

Тт

-1 - М 0 +

х_ (ТМ

Я Т х - (Т)

бх(0)=

Т

N 0 -

х_ (ТN. Я Т х - (Т)

ба,

из которого следует, что

бх(0) = Ьба,

(18)

где

Ь

Т,

I - М 0 +

х_ (Т)К' М

ЯТ х - (Т)

-1

Т

N 0

х_ (Т N

Я Т х - (Т)

На периодической траектории х^) в момент t = уТ происходит переключение управления с и = А на и = 0 (рис. 3). В соответствии с уравнением (10)

бх^)=

У1 (^0)бх(0)+г1^)ба при 0 < t <уТ,

У2 (t,уТ)бх + (уТ)+г2 (t,уТ)ба при уТ < t < Т.

(19)

Принимая во внимание полученные выше равенства, представим (19) в виде

0

0

[VI (г,0)Ь + Г1 (г)]5а при 0 < г < уТ,

5х(г) =

По определению функция чувствительности периодической траек-

[^(г, ут)(РЬ + О) + Г2(г, уТ)]5а при уТ < г < Т.

тории

Здесь

8х(г) х(г) + 8х(г) - х(г) ..

_А2= 11т —^—= р(г). (21)

йа да®0 да

р(г)=

V1(г,0)ь+г1(г) при 0 < г <уТ,

V (г, уТ )(РЬ+О)+г 2 (г, ут) при уТ < г <т. Установим чувствительность момента уТ и периода автоколебаний

2 Т.

В соответствии с равенством (7)

д(уТ) = - кТ [У1( уТ ,0)Ь + г1(уТ )]да. (22)

ЯТх-(уТ) '

Из (22) следует, что функция чувствительности

й(уТ) =-ЯТ [У1( уТ ,0)Ь + г1(уТ)]

йа ЯТ х- (уТ) '

Принимая во внимание (8) и (20), найдем

дТ = ЯТ [У2(Т,уТ)(РЬ+О)+г2(Т,уТ)]да (23)

яТ х - (Т) .

Периодическая траектория х (г) имеет период 2(Т + дТ). Чувствительность периода автоколебаний к изменению параметра а

й (2Т) =_ 2ЯТ [У2(Т, уТ )(РЬ + О) + г 2(Т, уТ)] (24)

йа яТ х- (Т) '

Чувствительность критерия устойчивости колебаний. В релейных автоколебательных системах, как правило, оценивается асимптотическая орбитальная устойчивость автоколебаний [6,7]. Для релейной системы (1), (2) устойчивость периодической траектории х(г) оценивается [5] по собственным числам матрицы

с = С 2СЬ (25)

где

с, = V] (уТ ,0) - х -(уТТкГ у1( уТ ,0), (26)

ЯТ х - (уТ)

х - (Т )ЯТ V2(T, уТ) ЯТ х - (Т)

С2 = ^(Т,уТ)-А ^ ДТ- -^ , (27)

Определим устойчивость периодической траектории х^) = х^) + бх^). Возмущенную (малым изменением начального условия) траекторию обозначим х^) = х^) + бх^). Запишем уравнение в вариациях

dбx(t) Эf (х, а + ба, и)

Обозначим

Тогда

dt Эх

^^ ч Эf (х, а, и) М(х, а, и) = —4 7 Эх

Эf (х, а + ба, и)

бх^).

= М(х + бх, а + ба, и).

Эх

Выше было установлено, что

х^)=х^)+р^ )ба

и, следовательно,

^^=М(х+рба,а+ба,и)бх. Так как ба является малой величиной, то можно записать

^ = M(„„в + »(,,„«X, (28,

где

^х,а,и)[М(х+рба,а+ба,и ] dба

5а=0

Во избежание недоразумений отметим, что ба входит в определение траектории х^), т.е. ба и бх - независимые малые величины.

В [1] установлено, что линеаризованные по а (не учитываются величины, имеющие порядок малости выше первого) нормированные фундаментальные матрицы решений (в зависимости от величины управления) уравнения (28) задаются равенствами

У ^ ,0) = У1 ^ ,0) + У1а ^ ,0)ба, (29)

У2(t, уТ) = У2& уТ) + У2а (t, уТ )ба. (3 0)

(У

В равенстве (29) У^ (^0) является матрицей [1] ,столбцы которой образованы векторами т1^),т2(t), ..., тп(t), причем каждый вектор т1 (t) является решением уравнения

— = М(х, а, А)ш + N(x, а, А)бх 0 (t) (31)

dt

при т(0) = 0. Это решение зависит от функции дх 0(г). Входящий в уравнение (31) вектор дх 0(г) представляет собой столбец матрицы Vl(г,0), причем вектору тг (г) соответствует г-ый столбец матрицы Vl (г ,0).

Далее, вектор V® (г, уТ) является матрицей, столбцы которой образованы векторами д1(г), д 2(г ),..., д п (г) (г >уТ). Каждый вектор (г) является решением уравнения

Ип-=М(х,а,0)д+^х,а,0)дх0(г,уТ), г>уТ,

йг

при д(уТ)=0. В равенстве (32) вектор дх 0(г, уТ) является столбцом матрицы V2(г,уТ). Вектору (г) соответствует г-ый столбец матрицы V2(г,уТ).

В соответствии с (25) - (27) матрица устойчивости периодической траектории ~х(г) задается равенством

С = С 2С1,

где

х (уТ+д(уТ ))К Т уТ+д (уТ ),0) ЯТ х (уТ+д( уТ))

С1 = уТ+д(уТ),0)^ ХТ ... 11 , (32)

С2 = ^(Т + 5Т, уТ + 5( уТ)) - ^(Т + дТ )яТ (Т + дТ, уТ + д(уТ)) (33)

х(Т + 5Т)

Рассмотрим сначала матрицу С1. Принимая во внимание равенство (29), вид функции х(г) и учитывая, что дх(г), да и д(уТ) являются малыми величинами, получим

СС1 = v1(ут ,0) + ^(уТ ,0)уда + Vх (уТ ,0)да -

- х - (уТ )ЯТ V (уТ ,0) + х- (уТ )ЯТ V (уТ ,0)уда + у ЯТ V (уТ ,0)да -ЯТ [х- (уТ) + х- (уТ)уда + уда]

х - (уТ)ЯТ V (уТ ,0)уда + х - (уТ )ЯТ V1a (уТ ,0)да ЯТ [х - (уТ) + х- (уТ)уда + уда]

у=-Я Т [V(уT ,0)Ь + г1(уТ ,0)]

ЯТ х - (Т ) '

/ ™ ^ ЭГ(х(уТ),а,А)

у=М(х(уТ),а,А)р (уТ) + 4 и ь ' .

Эа

В равенстве (34) опущены величины, имеющие порядок малости выше первого относительно да.

(34)

Здесь

Ниже равенства, в которых опущены величины, имеющие порядок малости выше первого, будем записывать, не оговаривая это особо. Представим (34) в виде

( = У1 (уТ ,0) + (У (уТ ,0)у + У а (уТ ,0))ба -

х - (уТ)ЯТ У1(уТ ,0) + (3 ба (35)

ЯТ [х (уТ) + (х (уТ)у + у)ба]

Здесь

О = х- (уТ )ЯТ У1 (уТ ,0)у + у ЯТ У1 (уТ ,0) + х - (уТ )ЯТ У&1 (уТ ,0)у +

+ х - (уТ )Я Т У а (уТ ,0). Выполним линеаризацию по ба входящей в равенство (35) дроби в окрестности точки ба=0 . В результате получим:

= ^^ + ^^аба,

где

^=У1( уТ ,0) -х - (утТ У1( уТ ,0),

1 1 ЯТ х- (уТ)

W1а = У (уТ,0)у - У а (уТ,0) - (3 6)

- ЯТх- (уТ)( - Я Т (х- (уТ)у+у )х- (уТ)ЯТУ1 (уТ,0)

[ЯТ х - (уТ )]2 '

На периодической траектории x(t) переключение управления с и=А на и=0 происходит в момент уТ, а на х ^)- в момент уТ+б(уТ) (б(уТ)=уба).

Рассмотрим матрицу С 2^):

(2 = У2(Т+бТ, уТ+б( уТ)) - ^ - (Т+бТ )ЯТ У2(Т+бТ, ут+б(уТ)). (37) 2 24 ЯТ х- (Т+бТ)

Выше было установлено

)=V2(t)+У2а (t )ба. Представим равенство (37) в виде

С2 = У2(Т+тба, уТ+■уба)+Уа (Т+цба, уТ+уба) -

[х- (Т+тба)+бх - (Т+тба)]ЯТ х (38)

ЯТ [х - (Т+тба)+бх - (Т+тба)] х[У2 (Т+тба, уТ+■уба)+Уа (Т+цба, уТ+уба)],

484

где

ЯТ [У(Т,уТ)(РЬ+О)+г 2(Т, уТ)] (39)

ЯТх-(Т) .

В равенстве (38) нормированные фундаментальные матрицы решений У2 и У® записаны следующим образом: выражение, стоящее в скобках перед запятой, задает конечное время, а после запятой - начальное.

Фундаментальная матрица решений У2 ^) формируется с точки переключения управления уТ (начальная точка У2(уТ)=I), а матрица У2^) с начальной точки уТ+б(уТ)=уТ+уба. Легко видеть (см. [8]), что У2 ^, уТ+■уба)=У2 ^ (уТ+■уба)=У2 (t )[У^1 (уТ)+У- (уТ )уба]

является нормированной фундаментальной матрицей решений и, следовательно,

У2 (Т+цба, уТ+уба)=У2 (Т+тба)[У-1 (уТ)+У- (уТ )уба]. (40) Далее (см. [8])

У )=-У -1(t )У (t )У -1(t).

Таким образом,

У2 (Т+тба, уТ +уба)=У2 (Т+тба)[У2"1 (уТ) - У-1 (уТ )У2 (уТ )У2-1 (уТ )уба]. Матрица У2(уТ)=I, здесь I - единичная матрица, и, следовательно У2 (Т+цба, уТ+уба)=У2 (Т+цба) [I - ^1 (уТ )Iуба]= =У2(Т+тба)Р - У2(уТ )уба].

Представим

У2(Т+тба)=У2(Т)+У2(Т )тба. Равенство (41) принимает окончательный вид

У2 (Т+тба, уТ +уба)=У2 (Т) - У2 (Т )У2 (уТ )уба+У2 (Т )|1ба=

=У2(Т)+[УУ2(Т )т - У2 (Т )У2 (уТ )у]ба. (42)

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

У2 (t, уТ)=У2 (t, уТ)+у® (t, уТ )ба, где Уа (уТ, уТ)=0. Тогда

У2а (уТ+уба, уТ )ба=[Уа (уТ, уТ)+У® (уТ )уба]ба=0. (43)

В равенстве (43), как и выше, опущено слагаемое, имеющее порядок малости относительно ба выше первого. Таким образом, матрица

У2(^ уТ+б(уТ))=У2(ОУ2~ЧуТ+б( уТ))+Уа (t, уТ+б(уТ ))ба (44)

является нормированной фундаментальной матрицей решений для параметрически возмущенной системы. В соответствии с (44) и (42)

~(Т+|да, уТ+уда)=V2(T)+05а,

где

0 = V2 (Т )|- V2 (Т )^2 ( уТ )у+V" (Т ). Упростим равенство (38):

С = v (Т)+Хда-[х-(Т)+х-(Т)||да+дх-(Т)]КТ[V2(T)+Хда] 2 2 Ят [х-(Т)+х-(Т)|да+дх-(Т)] .

В соответствии с (10)

дх - (Т)=уда,

где

ЭГ (х(Т ),а,0)

У =

М(х(Т ),а,0)р(Т)+

Эа

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

х-(Т) = п.

Здесь

П=М(х(Т ),а,0)х - (Т). Матрица С 2 принимает вид

С2 = V2(T)+.да-[х - (Т) +ТГба]КТ ^)+^да], 2 Ятх-(Т)+ЯТГда

где

Г=(п|+у). Выполним линеаризацию дроби

[х - (Т)+Гда]Я Т [V2(T)+0да] = х - (Т )ЯТ V2(T)+%да Ятх-(Т)+Я т Гда Ятх-(Т)+ЯТГда

по да в окрестности точки да=0. В результате получим, что

С=V2(T) - х ' (ТТКГ V2(T) - КТ х '(Т « - Г х - (2)КТ ^^) да+0да. (46) 2 Ятх-(Т) [Ятх-(Т)]2

В равенствах (45) и (46)

%=х - (Т )ЯТ 0+ГЯТ У2(T).

Введем обозначения:

х - (Т )Я т V2(T )

(45)

^2 = V2(T )-

Ят х - (Т)

W а = 0 - ЯТ х - (Т )% - Ят Гх - (Т )ЯТ V2(T )

2 [Ят х- (Т )]2 .

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

С2 = W2 + W2zда. (47)

Принимая во внимание (36) и (47), запишем для периодической траектории х(г) матрицу устойчивости:

С=ё2ё1 = W2 W1 + W1 + W2 W1a )да. (48)

Устойчивость определяется по собственным числам матрицы С . Если все ее собственные числа 1 у удовлетворяют неравенству

1,

<1, у =1,т,

то периодическая траектория асимпототически орбитально устойчива [5-7]. Собственные числа матрицы С являются корнями многочлена

ёе1[11 -С]. (49)

Обозначим 10, у=1,т, т £ п, собственные числа матрицы (48), если в ней положить да=0. По собственным числам 1® оценивается устойчивость периодической траектории х(г) . Представим

1 у =10у +61 у ,

где приращение д1 у порождено вариацией да параметра а. Ясно, что

61 у ®0 при да®0. Собственные числа 10 являются корнями многочлена (49) при да=0.

Пусть 10 - простое (не кратное) собственное число матрицы (25)

(матрицы (48) при да=0). Положим в многочлене (49) 1=1° +д1 у . Поскольку д1 у и да являются зависимыми малыми величинами, то в многочлене ёе1;[(10 +д1 у )1 - С ] можно приравнять к нулю произведения дад1, а также (д1 у)п при п > 2 как имеющие порядок малости выше первого. Собственное число 10 является корнем многочлена (49) при да=0 и поэтому

после указанных сокращений получим многочлен первого порядка относительно д1 у, т.е. равенство вида

/у (10 )д1 у = Гу (10у )да, (50)

где /у (10) и Гу (10) - некоторые многочлены.

487

Из (50) следует, что

Р (I0/)

б1 / = ] ба,

] // (Ц)

т.е. функция чувствительности собственного числа

л 7 а°7+б1 7-1°7 р7(А°7)

-Л. = 11т —7-7-7 = 7 7 . (51)

dа ба®° ба // (1 °)

В [1] показано, что равенства (50) и (51) справедливы и в случае кратного собственного числа . Но тогда (51) приводит к неопределенности вида «ноль», деленный на «ноль», и требуется ее раскрытие, например, по правилу Лопиталя.

Получение функций чувствительности рассмотрено при изменении одного скалярного параметра объекта управления. При изменении в объекте нескольких параметров функции чувствительности по каждому параметру определяются изложенным выше способом. Отклонения характеристик системы, обусловленные изменением нескольких параметров, находятся путем суммирования отклонений, порожденных вариацией каждого из параметров.

В работе для получения функций чувствительности используются точные методы исследования автоколебаний, это гарантирует справедливость самих функций чувствительности.

Работа выполнена при финансовой поддержке РФФИ (грант №1408-00662).

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

1. Моржов А.В., Фалдин Н.В. Функции чувствительности характеристик автоколебаний в релейных системах с нелинейным объектом управления // Известия РАН. Теория и системы управления. 2014. №6. С. 14 - 24.

2. Фалдин Н.В., Моржов А.В. Чувствительность ошибки слежения к изменению параметров объекта управления в релейной автоколебательной системе // Мехатроника, автоматизация, управление. 2015. Т. 16. № 2. С. 81 - 88.

3. Фалдин Н.В., Моржов А.В. Чувствительность вынужденных периодических движений релейной системы к изменению параметров объекта управления // Мехатроника, автоматизация, управление. 2015. Т. 16. № 11. С. 721 - 730.

4. Моржов А.В., Фалдин Н.В. Чувствительность ошибки слежения в релейной системе, работающей в режиме вынужденных колебаний // Известия РАН. Теория и системы управления. 2016. №3. С. 84 - 96.

488

5. Фалдин Н.В. Релейные системы автоматического управления /под. ред. К.А. Пупкова, Н.Д. Егупова // Математические модели, динамические характеристики и анализ систем автоматического управления. М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. С. 573 - 636.

6. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М.: Физматгиз, 1959. 914 с.

7. Ким Д.П. Теория автоматического управления: учеб. пособие для вузов. Т. 2. Многомерные, нелинейные, оптимальные и адаптивные системы. М.: Физматлит, 2004. 464 с.

8. Демидович Б.П. Лекции по математической теории устойчивости. М.: Наука, 1967. 472 с.

Фалдин Николай Васильевич, д-р техн. наук, проф., [email protected], Россия, Тула, Тульский государственный университет,

Моржов Александр Владимирович, канд. техн. наук, доц., morzhovamail.ru, Россия, Тула, Тульский государственный университет

SENSITIVITY OF CHARACTERISTICS OF SELF-OSCILLA TIONS IN SYSTEMS WITH

THREE-POSITION RELAY BLOCK

N.V. Faldin, A.V. Morzhov

Systems with three-position relay element and nonlinear plant is offered. Methods of determining of self-oscillations sensitivity functions (periodic trajectory, frequency, criterion of stability) are designed.

Key words: three-position relay element, self-oscillations, sensitivity, periodic trajectory, frequency, criterion of stability.

Nikolay Vasilyevich Faldin, doctor of technical sciences, professor, nvfal-dinayandex. ru, Russia, Tula, Tula State University,

Alexander Vladimirovich Morzhov, candidate of technical sciences, docent, morzhova mail.ru, Russia, Tula, Tula State University