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

Баев Теория колебаниы 2015

.pdf
Скачиваний:
22
Добавлен:
12.11.2022
Размер:
2.85 Mб
Скачать

Рассмотрим теперь функцию

 

 

представляющую левую часть

характеристического уравнения

 

(2.1.82). В качестве замкнутого

4 ] ,

 

 

 

 

 

 

 

 

 

 

контура

(С), о котором говорится в принципе аргумента, возьмем

контур,

составленный

из

полуокружности

 

большого

 

радиуса,

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

плоскости корней, и мнимой

 

 

 

 

 

 

 

 

оси, как

показано на

рис. 2.1.7. Поскольку

функция

 

 

,

являясь

полиномом конечного порядка, не имеет полюсов, то

выражении

4в ]

 

 

(2.1.93) в правой части останется только величина

 

, равная числу

нулей, если таковые имеются у функции в

правой

полуплоскости

 

 

 

 

 

 

 

 

плоскости корней.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

радиус окружности

 

 

выбрать достаточно большим и

окажется, что внутри этого

контура функция

 

 

не имеет нулей, то

 

 

характеристического уравнения

это будет означать отсутствие

у

 

 

 

4 ]

 

 

 

 

 

 

корней в правой полуплоскости, т.е. означать устойчивость движения

системы.

(z)

0

Рис. 2.1.7. Контур интегрирования

Таким образом, можно сформулировать следующий критерий] устойчивости: система устойчива, если при обходе точкой замкнутого4 ]контура (С), изображенного на рис. 2.1.7, число оборотов вектора , представляющего правую часть характеристического4 уравнения, вокруг начала координат на комплексной плоскости ( равно нулю.

В соответствии с этим критерием, если положить

] ρQT,

2.1.94

111

 

то

на полуокружности

 

 

.

достаточно

 

большого радиуса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точности аппроксимировать функцию

 

 

можно с любой степенью

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ ´ 1

первым слагаемым суммы (2.1.82)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ]

 

 

 

 

]

 

 

 

 

 

 

 

 

 

 

] r ρ Q T

 

 

 

 

 

 

 

4 ]

 

2.1.95

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, следовательно, после прохождения в положительном направлении

 

полных

 

 

 

 

 

 

 

 

4 ]

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

точкой

 

 

полуокружности

 

.

 

 

аргумент

функции

 

получит

приращение

 

 

, т.е. вектор

 

 

 

в комплексной плоскости (

 

сделает

 

 

 

 

 

оборота против часовой стрелки.

 

 

 

 

 

 

 

 

 

 

%/2Тогда, согласно сформулированному выше

 

критерию,

 

система

 

 

 

 

l∞

 

 

 

 

 

4 ]

 

 

 

 

 

%/2

 

 

 

 

 

по мнимой оси от

будет устойчива, если при перемещении точки

 

(+j

 

до (

 

 

 

вектор

 

 

 

сделает

 

 

полных]оборота по часовой

стрелке.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку,

 

 

 

согласно

 

 

принятому

 

 

выше

 

ограничению,

коэффициенты

 

 

 

полинома

 

 

 

 

 

 

являются

 

 

действительными

величинами,

то

0

 

до

 

 

 

 

 

4 ], и поэтому достаточно изменять

параметр

 

от

 

4 lω 4 lω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

плоскостиω(

 

 

 

в

случае устойчивой системы должен повернуться

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

4 ]

 

 

 

 

 

 

 

против

часовой стрелки на угол

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

состоит критерий

устойчивости

 

А.В.

 

В этом,

собственно,

и

 

 

%π/4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Михайлова.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий Михайлова. Для устойчивости системы необходимо и

достаточно, чтобы при изменении

 

от 0 до

 

 

кривая, описываемая

 

концом

 

 

вектора

 

 

 

 

 

 

,

 

прошла

 

протии

 

часовой

 

стрелки

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

плоскости,

где

 

-

последовательно

 

 

квадрантов

 

комплексной

 

 

4 lω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%

 

порядок

характеристического уравнения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В качестве примера на рис. 2.1.8,а приведены типичные годографы

вектора

4 lω

для устойчивых

 

систем; на рис. 2.1.8,б изображен

 

 

 

 

 

 

 

 

 

 

 

 

 

системы

с

 

% 4,

она

неустойчива,

годограф

вектора

 

третий квадрант.

 

 

 

 

 

 

 

 

 

 

поскольку пропущен

 

 

 

 

 

 

 

 

 

 

4 lω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Число корней с положительной вещественной частью называется

порядком неустойчивости системы. Обозначим это число символом

D. Из принципа аргумента следует, что

 

 

 

 

 

2π^ %π 2Bπ,

2.1.96

где

 

угол поворота вектора

4 lω

в положительном направлении

 

изменении

 

от 0 до .

 

 

при

 

ω

 

 

 

 

 

 

 

112

 

 

n=2 jv

 

4

 

 

 

 

 

jv

(4

(

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=5

 

 

 

n=4

 

 

 

 

 

 

u

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

n=3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=4

 

 

 

 

 

 

 

 

 

 

а

 

 

 

4 lω

 

б

 

 

 

неустойчивой (б)

 

устойчивых (а) и

 

Рис. 2.1.8. Годографы вектора

 

 

 

 

 

 

 

 

 

систем

 

 

Поэтому

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.1.97

 

 

 

 

 

 

^ ¶2 % 2B · ;

 

здесь

 

порядок

характеристического

уравнения, а квадратные

 

означают целую часть числа (

 

)/2.

 

скобки %

 

 

 

 

 

следует еще один способ проверки

Из

критерия

 

Михайлова

 

% 2B

 

 

устойчивости системы. Выделим действительную и мнимую

составляющие комплексной функции

4 lω :

2.1.98

4 lω ¢ ω

lP ω ,

тогда, согласно критерию Михайлова, система устойчива тогда и

только тогда, когда в области ω ¤ 0 уравнения

2.1.99

¢ ω 0;

P

 

 

2.1.100

 

ω113 0

имеют только действительные, перемежающиеся корни, т.е. между двумя соседними корнями одного из этих уравнений лежит корень другого уравнения, как показано на рис. 2.1.9,а.Если же хотя бы в одном случае между двумя соседними корнями одного из уравнений (2.1.99) или (2.1.100) отсутствует корень другого уравнения, как на рис 2.1.9,б, то система неустойчива.

Критерий Михайлова не позволяет как в критерии Рауса – Гурвица получать условия устойчивости в аналитическом виде, зато он сравнительно просто справляется с системами с большим числом

степеней свободы. Хотя в критерии Михайлова параметр

должен

меняться от 0 до

 

, реально достаточно построить

соответствующие

 

чтобы

ω

графики в пределах нескольких

единиц,

 

получить

 

 

 

 

 

 

представление об устойчивости системы.

 

 

 

 

u,v ¢ ω

 

P ω

u,v

u(ω)

 

 

 

ω

v(ω)

 

ω

а

 

б

 

 

 

Рис. 2.1.9. Критерий перемежающихся корней: а – устойчивое движение, б – неустойчивое

Д-разбиения. Нередки ситуации, когда некий параметр λ системы может принимать различные значения, которые в своей совокупности образуют область (λ) его возможных значений. В этом случае, естественно, желательно знать, при каких значениях параметра λ система устойчива и покрывают ли эти значения всю область (λ) или только ее часть, или вовсе ей не принадлежат.

Для решения подобной задачи изложенные выше критерии не всегда удобны. В критерии Рауса – Гурвица могут возникнуть трудности из-за громоздкости получающихся выражений, требующих

114

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

В критерии Михайлова определение границ области устойчивости для множества значений ее параметра 4λ может. потребовать многократного построения графиков функции

В этой связи представляет несомненный интерес так называемый метод Д-разбиений, суть которого состоит в следующем.

Ограничимся ниже случаем, когда характеристическое уравнение

линейно зависит от параметра λ и его можно привести к виду

 

 

 

 

 

4 λn

0

2.1.101

или

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

λ n .

 

2.1.102

Перейдем на комплексную плоскость (λ), считая параметр λ

функцией вида

 

4 lω

 

 

 

 

 

 

 

 

λ lω

 

 

 

 

 

2.1.103

 

 

n lω ¢ ω lP ω .

Задавая

значения

 

от

 

до

 

 

, построим

в плоскости

(λ)

 

которуюω

функция (2.1.103) отобразит мнимую ось

кривую, на

 

 

 

 

 

 

плоскости корней.

Разбиение

комплексной

плоскости

(λ),

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

половину, соответствующую полуоси

 

, а затем дополнить

границу Д-разбиения зеркальным

отображением ее

половины

0 ω N ∞

 

 

 

 

 

 

относительно действительной оси.

 

 

 

 

 

 

 

 

 

 

 

 

в

Если двигаться в плоскости корней вдоль мнимой оси из

 

 

 

 

 

 

Идя

, то область устойчивости будет находиться все время слева.

 

 

по границе Д-разбиения в том же направлении, т.е. от

 

к

 

 

,

 

 

 

Тогда

при

аналогично будем штриховать и границу разбиения.

 

 

переходе границы Д-разбиения с не заштрихованной стороны на заштрихованную число корней с положительной действительной

115

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

Пусть, к примеру, Д-разбиение выглядит так, как это показано на рис. 2.1.10.

Если «разметка» областей Д-разбиения выполнена, то сразу же определяется « претендент» на область устойчивости – это область с

наименьшим порядком неустойчивости,

?. В том случае, когда

заранее известно, что существует хотя

бы одно значение параметра ,

 

^

 

 

 

области

при котором система устойчива, работа по определению

λ

устойчивости на этом заканчивается – это область с ^ ?.

 

(z)

D=r+2

jP

 

(λ)

 

 

 

 

 

 

 

 

 

0

 

 

 

 

D=r+1

¢

¢

 

 

 

 

 

¢

 

 

D=r+3

 

 

 

 

D=r

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.1.10. Пример Д-разбиения

 

 

 

 

Если же нельзя утверждать,

что

 

? = 0, то это следует

проверять. Для

 

 

 

выяснить,

удовлетворяет ли

проверки достаточно ^

 

 

 

 

разбиения с минимальным значением величины ^, и λ область является областью устойчивости. В частности, если система,

какому-либо из критериев устойчивости данное характеристическое уравнение при любом выбранном значении параметра из области Д- если да, то эта

область устойчивости: ¢ N λ N ¢ .

Д-разбиение которой изображено на рис. 2.1.10, устойчива, то искомая

В качестве примера применим метод Д-разбиений для решения следующей задачи: найти условие возбуждения RC-генератора, изображенного на рис. 2.1.11 (входным током усилителя можно пренебречь).

116

 

$

R

R

R

K

C

$ C

$2 C

Рис. 2.1.11. Схема RC-генератора

Решение задачи. С помощью известных правил электротехники

несложно построить математическую модель системы

 

 

 

 

 

¸

;

$

2

&

1

 

$

 

& 0;

2.1.104

 

$

ž

 

ž

; $

 

$

 

 

1

 

 

 

 

 

 

 

1

 

 

$

2

& 0;

2.1.105

 

ž

; $

 

$

 

& ž ; $

 

 

 

$

2

 

 

1

; $

2

$

 

&

1

 

 

 

2

& 0,

2.1.106

 

 

ž

 

 

ž ; $

 

где К – коэффициент усиления усилителя.

 

 

 

Представим токи в контурах в виде экспонент

2.1.107

 

 

 

 

 

 

 

 

$

3 Q= # 1,2,3

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

1 ž 3 3 ¸32

0;

2.1.108

3 2 ž 3 32

0;

2.1.109

3 2 ž 32 0.

2.1.110

117

Эта система

алгебраических

уравнений

будет иметь

нетривиальное решение относительно амплитуд

3 (m=1,2,3), если ее

определитель равен нулю:

 

 

 

 

 

¢2 5¢ 6¢ ¸ 1 0,

 

2.1.111

где

 

¢ ž .

 

 

 

2.1.112

 

 

 

 

 

 

 

 

Д-разбиение для характеристического уравнения (2.1.111)

изображено на рис. 2.1.12.

 

 

 

 

 

 

 

 

jv

(K) ^? 1

 

 

 

 

 

 

2

 

 

 

^

 

 

^u

 

 

 

?

 

 

?

 

 

 

 

29

 

 

 

Рис. 2.1.12. Д-разбиение RC-генератора

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

методикой, выберем

точку в рассматриваемой

области

с

 

? и в

качестве таковой возьмем точку K=

 

.

В

 

этой

точке

 

 

 

^

характеристическое уравнение принимает вид 1

 

 

 

 

2.1.113

 

¢ 5¢ 6 0.

 

 

 

 

 

 

Оба корня этого уравнения

 

и

 

 

 

отрицательны,

т.е. лежат в правой

устойчивой

полуплоскости

плоскости

корней.

¢

2

 

¢

3

 

 

 

 

 

118

 

 

 

 

 

 

 

 

Следовательно, в выбранной точке

K=

 

система устойчива, а это

значит, что

?

и что она

устойчива во всей области с

 

 

 

1

 

схема будет работать как автогенератор, если ¸ ¤ 29. Задача решена.

минимальным значением порядка неустойчивости. Таким образом,

^

0

Впринципе, аналогичный анализ устойчивости можно проводить

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

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

Любую систему можно представить в виде четырехполюсника, изображенного на рис.2.1.13.

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

uвх °

 

 

 

° uвых

°

 

°

 

 

 

 

 

 

 

 

K1

°°

°0

Рис. 2.1.13. Четырехполюсник, охваченный обратной связью

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

119

которые могут самостоятельно «подстраиваться» под окружающую их

среду.

 

 

 

 

Пусть : – передаточная функция системы

 

 

 

 

uвых = :uвх.

(2.1.114)

Аргументом

 

передаточной функции служит

параметр

показательной

функции (2.1.65), в виде которых ищутся решения

 

 

 

 

дифференциальных уравнений, описывающих движение системы. Потребуем, чтобы система была устойчива при отсутствии

обратной связи, когда ключ К переключен на клемму с нулевым

индексом и uвх = 0. В этом случае должно выполняться равенство

1

2.1.115

: 0,

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

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

У линейных систем передаточные функции являются рациональными, представляющими в общем случае частное от

деления полиномов:

 

 

 

 

 

 

 

:

 

 

,

2.1.116

 

 

 

 

 

 

 

 

 

+

 

 

+

 

 

 

 

 

 

+

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Перейдем теперь к другому режиму работы рассматриваемого

четырехполюсника, а именно: к режиму, когда включена обратная

связь и uвх=uвых. При включенной обратной связи : =1 или

4 : 1 0.

2.1.117

120