
vostrikov
.pdf
10.2. Метод гармонического баланса |
|
|
|
|
|
|
|
351 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 10.2 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
7 |
|
|
1 |
|
|
… |
|
|
||
|
|
|
|
|
|
13 |
|
15 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Re |
1 |
|
– 0,083 |
0 |
|
5,983 |
|
1 |
|
|
… |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
Wл( j ) |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Im |
1 |
|
0 |
–0,135 |
0 |
|
|
0,667 |
|
… |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Wл( j ) |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Соответствующий график построен на рис. 10.11, где 4/А = 5,983. |
||||||||||||||||
|
|
Следовательно, параметры автоколебаний: |
|
|
|
|
Aa 0,669 . |
||||||||||||
|
|
a |
7 15, |
10.2.8.КОРРЕКЦИЯ АВТОКОЛЕБАНИЙ
Внекоторых случаях автоколебательный режим является желаемым режимом работы нелинейной системы, поэтому он должен иметь определенные амплитуду и частоту. Если параметры автоколебаний отличаются от требуемых, то возникает необходимость в их коррекции.
Сэтой целью можно воспользоваться следующими рекомендациями. 1. Если можно изменить значения коэффициентов линейной части,
следует попытаться подобрать их с учетом заданных параметров автоколебаний. Выбор коэффициентов линейной части осуществляется так, как описано в подразд. 10.2.5.
2.Если параметры линейной части нельзя изменить, то необходимо
еескорректировать. В этом случае на входе линейной части устанавливают дополнительное звено (корректор), которое рассчитывают любым известным в линейной теории методом синтеза.
3.При невозможности изменить линейную часть системы можно попытаться скорректировать нелинейный элемент.
10.2.9. УСЛОВИЯ ПРИМЕНИМОСТИ МЕТОДА ГАРМОНИЧЕСКОГО БАЛАНСА
Точность метода гармонического баланса зависит от точности замены нелинейного элемента эквивалентным линейным звеном, полученным в результате гармонической линеаризации. Отсюда следуют условия применимости метода гармонического баланса.

352 |
Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ |
1.Линейная часть системы должна быть низкочастотным фильтром, т. е. отфильтровывать возникающие на выходе нелинейного элемента все гармонические составляющие сигнала, кроме первой.
Для большинства систем, у которых степень полинома числителя передаточной функции меньше степени полинома ее знаменателя, это условие выполняется.
Кроме требования фильтрации, предъявляемого к линейной части, отметим случаи, когда в системе не будут возникать автоколебания.
2.При наличии однозначной статической нелинейной характеристики и передаточной функции линейной части, у которой в числителе
находится только коэффициент усиления (т. е. Wл ( p) k A( p) ), авто-
колебания в системе могут возникать только тогда, когда степень характеристического полинома n 3 .
3. В случае неоднозначной статической нелинейной характеристики и Wл ( p) k A( p) в системе может возникнуть автоколебательный режим, если n 2 .
10.3. МЕТОД МАЛОГО ПАРАМЕТРА
Этот метод обосновывает возможность пренебрежения малыми параметрами, которые присутствуют в математической модели и обычно не учитываются при проектировании реальной системы управления. Следует заметить, что подобное пренебрежение допустимо, если при этом не изменяются качественные свойства системы, а количественно они изменяются несущественно.
Рассмотрим суть метода малого параметра для системы, поведение которой описывают уравнения состояния [8]
x f (x) |
(x), x Rn . |
(10.28) |
Здесь f (x) и (x) – непрерывные дифференцируемые функции с ог-
раниченными производными, |
f (0) |
|
0 и |
(0) 0 , а нормы вектор- |
||||||
функций соизмеримы, |
|
~ |
|
(x) |
|
; – малый параметр, причем о |
||||
f (x) |
|
|
||||||||
его малости можно говорить, когда |
0,1 |
|
f (x) |
|
. |
|||||
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
Наряду с исходной системой будем рассматривать вырожденную,
которая получается из (10.26) при |
0 : |
|
x |
f (x) . |
(10.29) |

10.3. Метод малого параметра |
353 |
На вопрос о том, насколько процессы в системе (10.28) близки к процессам (10.29), отвечает следующее утверждение.
Теорема. Если вырожденная система (10.29) экспоненциально устойчива, то существует такое достаточно малое значение , при ко-
тором исходная система (10.26) также экспоненциально устойчива.
Для доказательства используем второй метод Ляпунова. Поскольку вырожденная система экспоненциально устойчива, для нее существует такая функция Ляпунова V (x) , что ее полная производная вдоль траек-
торий движения (10.29) удовлетворяет неравенству
V (x) |
V |
f (x) c |
|
|
|
x |
|
|
|
2 , c 0 . |
(10.30) |
|
|
|
|
|
|||||||||
T |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
x |
|
|
|
|
|
|
|
|
|
|
Используем эту же функцию V (x) |
для анализа устойчивости ис- |
||||||||||||||
ходной системы. Ее полная производная в силу (10.26) имеет вид |
|||||||||||||||
V (x) |
V |
|
|
f (x) |
|
V |
(x) . |
(10.31) |
|||||||
T |
|
T |
|||||||||||||
|
x |
|
|
|
|
|
|
|
|
x |
|
|
|||
С учетом (10.30) вместо (10.29) получим неравенство |
|
||||||||||||||
V (x) |
c |
|
|
|
x |
|
|
|
2 |
|
V |
|
(x) . |
(10.32) |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
xT |
|
xT |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Поскольку произведение |
V |
|
|
|
|
|
|
(x) также можно ограничить |
|||||||
сверху квадратом нормы, при достаточно малом значении |
вторая |
составляющая правой части (10.32) будет меньше по модулю, чем первая. При этом значение V (x) будет не просто меньше нуля, а меньше
некоторой квадратичной оценки. Следовательно, система (10.28) будет экспоненциально устойчива.
К сожалению, конкретное численное значение , начиная с которо-
го будет выполняться условие экспоненциальной устойчивости исходной системы, оценить довольно трудно. Применяемый с этой целью второй метод Ляпунова дает явно заниженные значения. В малой окрестности точки равновесия имеет смысл линеаризовать систему (10.28) и найти граничное значение параметра с помощью любого критерия
устойчивости.

354 |
Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ |
Отметим, что теорема позволяет упростить исследование устойчивости сложных нелинейных систем, так как достаточно проверить устойчивость более простой вырожденной системы.
10.4. МЕТОД РАЗДЕЛЕНИЯ ДВИЖЕНИЙ
Метод разделения движений также выражает идею малого параметра, но только по отношению к коэффициентам, определяющим инерционность отдельных элементов системы. В реальных ситуациях полное движение можно представить композицией процессов, протекающих с различными скоростями. При этом обычно пренебрегают быстрыми составляющими процесса, порожденными малыми инерционностями, что не всегда допустимо, так как может привести к ошибочным выводам относительно качественных свойств системы, в частности ее устойчивости. Метод разделения движений позволяет ответить на вопрос: «При каких условиях можно пренебречь малыми инерционностями?».
Суть метода заключается в представлении общей модели системы с малыми параметрами в виде совокупности двух подсистем, каждая из которых соответствует своему темпу процессов. В результате упрощается процедура ее анализа и синтеза: вместо сложной системы высокого порядка можно рассматривать две относительно простые подсистемы меньшего порядка.
10.4.1. ОБЩИЕ СВОЙСТВА СИСТЕМ
Математическая модель системы с разнотемповыми процессами представляет собой систему дифференциальных уравнений с малым параметром при части производных [7, 8, 10]:
x |
f (x, y), |
x Rn , |
(10.33) |
y |
|
y Rm , |
|
(x, y), |
|
где функции f ( ) и ( ) соизмеримы по норме в рабочей области про-
странства состояний.
Обсудим некоторые свойства этой системы, фазовый портрет которой изображен на рис. 10.12.

10.4. Метод разделения движений |
355 |
y |
(x, y) = 0 |
x |
Рис. 10.12. Пример фазового портрета системы с разнотемповыми процессами
Как видим, процессы имеют две фазы движения: из произвольных
начальных условий к поверхности |
(x, y) |
0 |
и вдоль нее. |
Выясним, |
|
как соотносятся при этом векторы скорости y |
и x [8]. |
|
|||
1. Во |
всем пространстве, |
кроме |
окрестности поверхности |
||
(x, y) |
0 , вектор скорости ориентирован почти параллельно коорди- |
||||
нате y . |
Это означает, что скорость изменения переменных |
y много |
выше, чем x . Они связаны следующим приближенным соотношением:
y |
|
|
|
1 |
|
|
|
x |
|
|
|
x, y : |
0 . |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. Если изображающая точка системы движется вдоль поверхности (x, y) 0 , то скорости изменения переменных будут соизмеримы. В
самом деле, поскольку в этом случае справедливо условие
t
x(t), y(t) 0 ,
полная производная функции ( ) по времени будет равна нулю, т. е.
(x, y) |
|
x |
|
|
y 0 . |
T |
y |
T |
|||
|
x |
|
|
Векторы y и x связаны конечным соотношением, что и означает их
соизмеримость.
Таким образом, при движении из произвольных начальных состояний вначале процесс развивается быстро в силу большей скорости

356 Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ
изменения переменной y . При подходе к (x, y) 0 модуль вектора y уменьшается, а на поверхности y и x становятся соизмеримыми, и
вдоль нее изображающая точка движется с «нормальной» скоростью. Интуитивно понятно, что пренебречь быстрыми движениями можно
в случае, когда процессы сходятся к поверхности (x, y) 0 за сущест-
венно меньшее время, чем полная длительность процессов системы. Конкретные рекомендации на этот счет и дает метод разделения движений.
10.4.2. ВЫДЕЛЕНИЕ ОТДЕЛЬНЫХ СОСТАВЛЯЮЩИХ ДВИЖЕНИЯ
Для исследования свойств необходимо научиться выделять отдельные составляющие движения: быструю и медленную.
Получим сначала подсистему медленных движений. С этой целью будем рассматривать асимптотику системы (10.33), полагая в ней
0 ,
|
x f (x, y0 ), |
(10.34) |
|
|
0 |
(x, y0 ), |
|
|
|
||
где через |
y0 обозначены значения y , соответствующие поверх- |
||
ности ( ) |
0 . |
|
|
|
Если теперь из последнего уравнения (10.34) найти y0 как функцию |
y0 |
(x) и подставить в первое, то получим описание независимой |
подсистемы медленных движений в виде дифференциального уравнения состояния n-го порядка
x f (x, (x)) f 0(x) . |
(10.35) |
Для выделения подсистемы быстрых движений удобно ввести масштаб времени
1t ,
что позволяет «растянуть» процессы и исследовать быструю составляющую.

10.4. Метод разделения движений |
|
|
357 |
||||
Поскольку t |
, описание исходной системы (10.33) в новом вре- |
||||||
мени принимает вид |
|
|
|
||||
|
|
|
dx |
|
f (x, y), |
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
||||
|
|
|
|
|
|||
|
|
|
dy |
|
(x, y) |
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
||
|
|
|
|
|
|
||
или окончательно |
|
|
|
||||
|
|
|
dx |
|
f (x, y), |
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
(10.36) |
||
|
|
|
|
|
|||
|
|
|
dy |
|
|
||
|
|
|
|
(x, y). |
|
||
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
||
|
|
|
|
|
|
||
Рассмотрим |
теперь асимптотику системы |
(10.34), полагая в ней |
0 . При этом модель быстрых процессов в новом времени τ пред-
ставляет собой следующее дифференциальное уравнение состояния m -го порядка:
x |
const, |
|
dy |
(10.37) |
|
|
(x, y). |
|
d |
||
|
Возвращаясь к обычному времени t , получим описание подсистемы быстрых движений в виде
x const ,
(10.38)
y (x, y).
Отметим, что системы (10.35) и (10.38) независимы и описывают две составляющие общего процесса исходной системы, причем сумма порядков подсистем медленных и быстрых движений равна порядку полной системы (10.31). Степень этого разделения зависит от численного значения .


10.4. Метод разделения движений |
359 |
а затем в нормальном времени: |
|
x1 |
const , |
x2 |
const , |
x3 |
x3 3u const . |
Отметим, что вместо нелинейной системы третьего порядка получили две независимые линейные подсистемы второго и первого порядков.
10.4.3. ОСНОВНЫЕ ТЕОРЕМЫ МЕТОДА РАЗДЕЛЕНИЯ ДВИЖЕНИЙ
Полученные в результате рассмотренной процедуры разделения независимые подсистемы быстрых и медленных движений можно исследовать отдельно. Вывод относительно свойств исходной системы позволяют сформулировать теоремы, которые мы далее приведем без доказательства [7, 10] и в упрощенном виде.
1. Теорема 1. Если подсистема быстрых движений (10.38) экспоненциально устойчива, то для любого как угодно малого момента вре-
мени t0 0 найдется такое достаточно малое значение параметра , что траектория движения исходной системы (10.33) будет находиться к моменту времени t0 в какой угодно малой окрестности по-
верхности (x, y) 0 .
Доказательство теоремы основано на использовании второго метода Ляпунова. Поскольку подсистема быстрых движений (10.38) экспоненциально устойчива, существующая для нее функция Ляпунова используется при оценке устойчивости полной системы (10.31) [8].
Теорема 1 технически означает, что при экспоненциальной устойчивости подсистемы быстрых движений поведение исходной системы будет близко к поведению медленных движений, т. е. процессы системы (10.33) можно с достаточной точностью оценивать по (10.35).
2. Теорема 2. Если подсистемы быстрых и медленных движений порознь экспоненциально устойчивы, то существует такое достаточно малое значение , что исходная система (10.33) также экспо-
ненциально устойчива [10].
Таким образом, эти теоремы устанавливают связь между свойствами полной системы и двумя упрощенными и дают обоснование замене исследования точного решения системы с малым параметром анализом приближенного, а именно решения подсистемы медленных движений.

360 Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ
ПРИМЕР 10.7
Используя метод разделения движений, определить устойчивость системы
x |
x |
x2 |
2x |
, |
|
1 |
1 |
|
1 |
2 |
|
x |
2x2 |
4x |
5u . |
||
2 |
|
1 |
|
2 |
|
Полагая = 0 , получим |
|
|
|
|
|
x |
x |
|
x2 |
2x |
, |
1 |
1 |
1 |
2 |
||
0 |
2x2 |
4x |
5u , |
||
|
1 |
|
2 |
|
|
откуда следует |
|
|
|
|
|
x |
0, 5x2 |
1, 25u . |
|||
2 |
|
1 |
|
|
|
В результате описание подсистемы медленных движений будет иметь вид x1 x1 2,5u.
Получили линейное уравнение первого порядка, которое соответствует устойчивой подсистеме.
Уравнение подсистемы быстрых движений следующее:
x1 const ,
x2 4x2 5u const.
Как видим, это линейное уравнение первого порядка, причем подсистема быстрых движений устойчива.
Согласно приведенным теоремам исходная система будет экспоненциально устойчива, а процессы в ней имеют монотонный характер в соответствии с описанием подсистемы медленных движений.
10.4.4.УСЛОВИЕ РАЗДЕЛИМОСТИ ДВИЖЕНИЙ
Кнастоящему времени метод разделения движений получил широкое применение, поскольку позволяет обоснованно представить сложную нелинейную систему в виде совокупности нескольких более простых подсистем. Это существенно упрощает задачу анализа и, что гораздо важнее, синтеза.
Наличие разнотемповых процессов в физической системе обычно обусловлено использованием в ней больших коэффициентов или