
vostrikov
.pdf
Задачи |
331 |
9.10. Решая матричное уравнение Ляпунова, проверить устойчивость системы, математическая модель которой имеет вид
x1 |
x1 |
2x2 , |
x2 |
x1 |
x2. |
9.11. Решая матричное уравнение Ляпунова, проверить устойчивость системы, математическая модель которой имеет вид
x1 |
x1 |
7x2 , |
x2 |
2x1 |
10x2. |
9.12. Вторым методом Ляпунова проверить устойчивость системы
(рис. 9.9), если W1 |
( p) |
|
1 |
|
и W2 ( p) |
5 |
|
. |
||
|
|
|
|
|
|
|||||
|
p |
2 |
3 p 1 |
|||||||
|
|
|
|
|
|
|||||
|
v |
|
|
|
|
|
|
y |
||
|
|
|
|
|
|
|
Wp() |
Wp() |
||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
2 |
|
|
Риc. 9.9. Структурная схема замкнутой системы к задаче 9.12
9.13. Вторым методом Ляпунова проверить устойчивость системы
(рис. 9.9), еслиW ( p) |
1 |
и W ( p) |
15 |
. |
|
|
|||
1 |
p |
2 |
0,5 p 1 |
|
|
|
9.14. Вторым методом Ляпунова определить устойчивость системы, уравнения которой имеют вид
x1 x2 ,
x2 5 x1 3x2 .
9.15. Вторым методом Ляпунова определить устойчивость системы, уравнения которой имеют вид
x |
x1 |
x , |
1 |
|
2 |
x |
x |
5x3 . |
2 |
1 |
2 |

332 Глава 9. УСТОЙЧИВОСТЬ НЕЛИНЕЙНЫХ СИСТЕМ
9.16. С помощью функции Ляпунова
V (x) 0,5 3x12 2x1x2 x22
проверить устойчивость системы, математическая модель которой имеет вид
x |
x |
x x |
|
x3 |
0,5x x2 |
, |
|
|||
1 |
2 |
1 |
2 |
|
1 |
|
1 |
2 |
|
|
x |
3x |
x x |
x2 x |
0,5x x2. |
||||||
2 |
2 |
|
1 |
2 |
1 |
2 |
|
|
1 |
2 |
9.17. Определить устойчивость системы, структурная схема которой представлена на рис. 9.10.
v |
|
|
|
k |
|
|
|
|
|
y |
|
1 |
y |
|||
|
|
|
|
Tp 1 |
|
|
|
|
|
|
p |
|
||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
z |
|
1 |
|
|
|
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
|
k |
(-) |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
Рис. 9.10. Структурная схема системы к задаче 9.17
9.18. Проверить условие абсолютной устойчивости системы (рис. 9.11) с однозначной нелинейной характеристикой u f ( ) , 0 f ( ) 3 ,
если
Wл ( p) |
|
|
|
6 |
|
. |
|||
|
|
|
|
|
|||||
p2 |
5 p |
1 p 1 |
|||||||
|
|
|
|
||||||
v |
u |
|
|
|
y |
||||
|
Wл(р) |
|
|||||||
|
|
НЭ |
|
|
|
|
|||
|
|
|
|
|
|||||
|
|
|
|
л ( p) |
|
|
Рис. 9.11. Структурная схема нелинейной системы к задаче 9.18

Задачи |
333 |
9.19. Проверить условие абсолютной устойчивости системы (рис. 9.12)
с однозначной нелинейной характеристикой u |
|
f ( ) , |
0 f ( |
) 0,5 |
, |
||||||
если Wл ( p) |
2 |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
p3 p2 |
3 p |
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
|||
9.20. Проверить |
условие |
абсолютной |
устойчивости |
системы |
|||||||
(рис. 9.11), если Wл ( p) |
|
|
10 |
|
|
, а однозначная нели- |
|||||
|
|
|
|
|
|||||||
|
( p |
5)( p |
2)( p |
1) |
|||||||
|
|
|
|
|
|
|
|
||||
нейная характеристика |
при |
С = 5 |
имеет вид, |
показанный |
на |
||||||
рис. 9.12. |
|
|
|
|
|
|
|
|
|
|
|
|
uu |
u |
u |
|
C |
C |
|
|
|
b |
|
1 |
1 |
|
b |
|
C |
|
C |
Рис. 9.12. Статическая |
|
|
|
|
|
|
Рис. 9.13. Статическая |
|||
нелинейная характеристика |
|
|
|
|
|
нелинейная характеристика |
||||
|
|
|
|
|
|
к задачам 9.21 и 9.22 |
||||
к задаче 9.20 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
9.21. Проверить |
условие |
|
|
абсолютной |
устойчивости |
системы |
||||
(см. рис. 9.11), если Wл ( p) |
|
|
|
|
2 |
|
, а однозначная нелиней- |
|||
|
|
|
|
|
|
|||||
|
|
p |
p2 |
5 p |
4 |
|||||
|
|
|
|
|
|
|
||||
ная характеристика при C |
10, |
b |
2 |
имеет вид, изображенный на |
||||||
рис. 9.13. |
|
|
|
|
|
|
|
|
|
|
9.22. Проверить |
условие |
|
|
абсолютной |
устойчивости |
системы |
||||
(см. рис. 9.11), если Wл ( p) |
|
|
|
10( p |
1) |
|
, а однозначная нели- |
|||
|
|
|
|
|
|
|
||||
|
( p |
1)( p2 |
6 p |
1) |
||||||
|
|
|
|
|
||||||
нейная характеристика при C |
|
|
2 , b |
0,5 имеет вид, изображенный на |
||||||
рис. 9.13. |
|
|
|
|
|
|
|
|
|
|

334 Глава 9. УСТОЙЧИВОСТЬ НЕЛИНЕЙНЫХ СИСТЕМ
|
u |
|
|
|
|
|
|
|
9.23. Проверить условие абсолютной ус- |
||||||
|
|
|
|
|
|||||||||||
|
C |
|
|
|
|
|
|
|
тойчивости |
системы |
(см. |
рис. 9.11), |
если |
||
– b |
|
|
|
|
45 |
|
|
||||||||
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
Wл ( p) |
|
|
, а однозначная нели- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
45 |
|
|
|
b |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
( p |
3)( p2 p 10) |
|
|
||
|
|
|
|
C |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нейная характеристика |
C |
3 имеет вид, |
пока- |
|||
|
|
|
|
|
|
|
|
Рис. 9.14. Статиче- |
занный на рис. 9.14, где C = 3, b = 0,5. |
|
|
ская нелинейная ха- |
|
рактеристика к за- |
|
даче 9.23 |
|

Г л а в а 10
АНАЛИЗ ПРОЦЕССОВ
ВНЕЛИНЕЙНЫХ СИСТЕМАХ
Вглаве 5 мы показали, что качество работы линейной системы можно оценить различными способами без непосредственного
расчета переходного процесса. Отметим, что задача имеет смысл только для устойчивой системы, и это априори всегда предполагается.
Поскольку характер процессов нелинейной системы, как и ее устойчивость, существенно зависит от значений внешних воздействий и начальных условий, методы линейной теории неприменимы. Трудности анализа и оценивания процессов соответствуют сложности решения нелинейных дифференциальных уравнений.
Одно время казалось, что возможности второго метода Ляпунова позволят разработать регулярные процедуры оценки частных решений дифференциальных уравнений, т. е. переходных процессов автоматических систем. К сожалению, метод Ляпунова дает слишком завышенные соотношения, которые на практике не применяются.
Таким образом, к настоящему времени не разработано общей теории анализа процессов в нелинейных системах, существуют лишь методики, которые позволяют решать отдельные задачи. В этом разделе мы рассмотрим некоторые из них.
10.1. МЕТОД ФАЗОВОЙ ПЛОСКОСТИ
Метод фазовой плоскости применяется для анализа свойств систем второго порядка и основан на использовании аппарата пространства состояний. Суть метода заключается в отображении частных решений дифференциальных уравнений в совокупность фазовых траекторий.

336 Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ
Обсудим способы построения фазового портрета системы, математическая модель которой имеет вид
x |
f |
(x , x , u), |
|
||
1 |
1 |
|
1 |
2 |
(10.1) |
x |
f |
|
|
|
|
|
(x , x , u). |
|
|||
2 |
2 |
1 |
2 |
|
|
Управляющее воздействие u |
входит в правую часть (10.1) как па- |
раметр. Отметим, что если оно изменяется, то векторное поле будет управляемым. Здесь полагаем u const и его численное значение учтем в соответствующих функциях, что позволяет модель (10.1) записать в форме
x1 |
f1 |
(x1 |
, x2 ), |
(10.2) |
|
x2 |
f2 (x1, x2 ). |
||||
|
В принципе, задавая множество наборов значений x1 и x2 , можно получить поле векторов скорости (рис. 10.1) и, двигаясь вдоль них, по-
|
|
|
|
|
|
|
|
|
|
строить фазовую траекторию системы из |
|
|
xn |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
xn |
|
|
|
|
|
|
|
определенных начальных условий. |
||
|
|
|
x(0) |
|
|
|
|
|
|||
|
|
|
|
|
Таким образом, мы геометрически |
||||||
|
|
|
|
х( ) |
|
|
|
||||
|
|
|
|
|
|
x(t) |
определили решение дифференциально- |
||||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
х(ti ) |
|||||
|
|
|
|
|
|
|
|
|
|
го уравнения (10.1) для конкретных на- |
|
|
|
|
|
|
|
|
|
|
|
чальных условий. Однако в настоящее |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
x11 |
|
|||
|
|
|
|
|
|
|
|
время этот способ не находит примене- |
|||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
ния, так как наличие развитых средств |
|
|
|
|
|
|
|
|
|
|
|
вычислительной техники позволяет лег- |
|
|
Рис. 10.1. Пример построе- |
ко получить требуемую |
совокупность |
||||||||
|
решений. |
|
|||||||||
|
ния фазовой траектории |
|
|||||||||
|
Интерес представляют |
упрощенные |
|||||||||
|
|
|
|
|
|
|
|
|
|
аналитические способы построения фазового портрета системы, одним из которых является способ изоклин, когда не требуется информация о значениях вектора скорости.
Изоклиной называют линию в пространстве состояний, соединяющую все точки пространства с одинаковым наклоном фазовых траекто-
рий. Выражение для коэффициента наклона траекторий k |
имеет вид |
||||||
k |
x2 |
|
f2 |
(x1, x2 ) |
, |
(10.3) |
|
x |
|
f (x , x ) |
|||||
|
|
|
|
||||
|
1 |
1 |
1 |
2 |
|
|

10.1. Метод фазовой плоскости |
337 |
откуда чаще всего удается получить явное аналитическое описание семейства изоклин
x2 |
(x1,k) . |
Задаваясь рядом численных |
значений k из диапазона от –∞ |
до +∞, получим из (10.4) совокупность конкретных изоклин. Следует отметить, что чем больше построено изоклин, тем полнее и точнее будет фазовый портрет, технику построения которого поясним на примере.
ПРИМЕР 10.1
Построить фазовый портрет системы, поведение которой описывает уравнение
y 0,5y y 2u ,
при численном значении управляющего воздействия u 1 . Предварительно перейдем к описанию системы в переменных состоя-
ния, полагая x1 y, x2 y , |
|
x1 |
x2 , |
x2 |
x1 0,5x2 2. |
В соответствии с (10.3) определим коэффициент наклона фазовых траекторий:
k |
x2 |
|
x1 |
0,5x2 2 |
x1 |
|
|
x2 |
|
|
|
|
||
и запишем уравнение семейства изоклин |
|
x |
2 |
x1 |
. |
|
|
||
2 |
k |
0, 5 |
|
|
|
Придавая коэффициенту k определенные численные значения, получим уравнения конкретных изоклин (табл. 10.1).
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 10.1 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
0 |
|
|
1 |
|
|
–1 |
… |
|
∞ |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
arctg k |
0 |
|
|
|
|
|
|
|
|
… |
|
|
|
|
4 |
|
|
4 |
|
|
2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|||||
Изоклина |
4 2x1 |
4 |
2x1 |
4 2x1 |
… |
|
x2 0 |
|||||||
|
|
3 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|


10.2. Метод гармонического баланса |
339 |
а все линейные звенья объединены в одно с передаточной функцией Wл ( p) . Структурная схема рассматриваемых систем изображена на рис. 10.3.
v |
|
u |
y |
|
|
НЭ |
W (рp)) |
|
|
||
|
|
|
лл |
Рис. 10.3. Структурная схема исследуемой системы
Метод гармонического баланса пригоден для исследования автоколебательных систем практически любого порядка, но требует обеспечения условий хорошей фильтрации возникающих на выходе нелинейного звена гармонических составляющих сигнала (выше первой гармоники).
В основе расчетных соотношений метода гармонического баланса лежит способ гармонической линеаризации нелинейного элемента [6, 19], который мы далее и рассмотрим.
10.2.2. МЕТОД ГАРМОНИЧЕСКОЙ ЛИНЕАРИЗАЦИИ
Под линеаризацией понимают приближенную замену нелинейной функции линейной таким образом, чтобы по како- му-то выбранному показателю обе эти функции
совпадали. В способе гармонической линеаризации нели-
нейный элемент (рис. 10.4) заменяется квазилинейным звеном, параметры которого определяются при синусоидальном входном сигнале
Asin( t)
из условия равенства амплитуд первых гармоник на выходе нелинейного элемента и эквивалентного ему линейного звена.
Рассмотрим процедуру линеаризации для нелинейного элемента, уравнение которого примем в виде
u f ( , ) . |
(10.5) |

340 |
Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ |
При поступлении на его вход гармонического сигнала (10.4) на выходе звена в установившемся режиме также будет периодический, но несинусоидальный сигнал
u |
f Asin( |
t), A cos( |
t) f ( A, |
t) . |
(10.6) |
Разложим его в ряд Фурье [40] и получим |
|
|
|||
u |
u0 |
bk sin(k t) |
ck cos(k |
t) , |
(10.7) |
|
k 1 |
|
|
|
|
где будем полагать |
u0 0 , |
что справедливо для симметричной нели- |
нейной характеристики (10.6).
С учетом (10.6) коэффициенты ряда Фурье (10.7) определяются известными соотношениями
1 |
|
2 |
|
|
|
||
bk |
|
|
f ( A, |
t)sin(k |
t) d ( |
t), |
|
|
|
||||||
|
|
|
|
0 |
|
|
|
1 |
|
2 |
|
|
|
||
ck |
|
|
f ( A, |
t)cos(k |
t) d ( |
t). |
|
|
|||||||
|
|
|
|
0 |
|
|
|
Используем только первые члены ряда разложения в (10.7), пренебрегая высшими гармониками, и получим
|
u |
b1 sin( |
t) |
c1 cos( t) . |
(10.8) |
||||
Учтем, что |
Asin( t) , а |
A |
cos( |
t) , следовательно, |
|
||||
|
|
sin( |
t) |
|
|
|
, |
|
|
|
|
|
A |
(10.9) |
|||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||
|
|
cos( |
t) |
|
|
|
. |
|
|
|
|
|
|
A |
|
После подстановки (10.9) в (10.8) получим выражение для выходного сигнала нелинейного звена
u |
b1 |
|
c1 |
, |
|
A |
|
A |
|
||
|
|
|