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

3832

.pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
25.44 Mб
Скачать

ma Asin( t2 )

или sin t2 sin( t2)

ma

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

Подставив значения sin t1 и sin t2 в выражение для q(A), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

m

 

 

a 2

 

 

 

 

ma

2

 

 

 

 

 

 

q(A)

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

A

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мнимая составляющая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(A)

2

t2

 

 

cos d

 

2 m

[sin t

sin t

]

2 ma(1 m)

.

 

 

t

m

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

A

 

1

 

 

 

2

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В нормированной форме полученные выражения имеют вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

m 2

 

 

 

q(A) kн

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

A/

 

 

 

 

 

 

 

A/ a

 

 

 

 

 

 

 

a

 

 

 

 

 

 

A/ a

 

 

 

 

 

 

 

 

(4.84)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(A) kн

2(1 m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A/ a

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При m=1 из (4.84) может быть получено выражение (4.82) для одно-

значной релейной характеристики с зоной нечувствительности.

В случае безынерционных нелинейных элементов для вычисления ко-

эффициентов гармонической линеаризации необязательно получение перио-

дического сигнала на выходе /7/. Пусть нелинейный элемент имеет двух-

значную статическую характеристику (рис. 4.34) с восходящей F1 и нисхо-

дящей F2 ветвями.

В интегралах (4.73) заменим путем подстановок на переменную x:

sin x / A;

cos d dx / A;

 

 

 

Arcsin(x / A);

d

 

dx

 

,

 

 

 

A2 x2

 

 

 

 

 

где знак «+» имеет место для значений

на интервалах (0,

/2); (3 /2,

2 ),

а

знак «-» - для интервалов ( /2,

); ( ,

Рис. 4.33. Нелинейные двухзначные

3 /2). Этим интервалам (с учетом за-

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

 

 

391

мены переменных) соответствуют следующие пределы интегрирования по x: (0,А); (А,0); (0,-А); (-А,0). В результате формулы (4.73) для коэффициентов гармонической линеаризации нелинейного элемента с характеристиками,

приведенными на рис. 4.34, запишутся так:

 

 

1

 

 

 

 

A F (x)xdx

 

 

 

0

F (x)xdx

 

 

 

A

 

F (x)xdx

0

 

F (x)xdx

 

q(A)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

A

 

2

 

 

 

 

 

 

0

2

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2 x2

 

 

A2 x2

 

 

A2 x2

 

 

A2 x2

 

 

 

 

 

A

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

1

 

 

A

 

 

F (x)xdx

 

A

 

 

F (x)xdx

1

 

 

 

A

F (x) F (x)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

xdx.

 

A

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

2

 

 

 

 

 

 

 

 

 

 

2

x

2

 

 

 

A

2

x

2

 

 

A

2

2

 

 

 

 

A

 

 

 

A

 

 

A

 

 

 

 

 

 

 

 

 

 

 

A

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

A

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

0

 

 

 

 

 

 

b(A)

 

 

 

 

 

F1(x)dx F2(x)dx F2 (x)dx F1(x)dx

 

 

 

A

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

1

 

 

A

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

1

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1(x)dx F2 (x)dx

 

 

 

 

 

 

 

[F1(x) F2(x)]dx.

 

 

 

A

2

 

A

2

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

(4.85)

(4.86)

Выражения (4.85), (4.86) справедливы, если амплитуда А входного сиг-

нала охватывает весь интервал двузначности статической характеристики не-

линейного элемента. Эти формулы, по существу, означают усреднение по множеству значений x сигнала на его входе.

Обратим внимание на то, что коэффициент b(А) оказывается пропор-

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

если F1(x) F2(x) для всех x, то b(А) 0, и наоборот. Это значит, что по виду характеристики можно судить о знаке b(А) и, следовательно, о знаке фазово-

го сдвига: элемент, для которого восходящая ветвь характеристики проходит выше нисходящей (рис. 4.34, а), вносит по первой гармонической состав-

ляющей положительный фазовый сдвиг, а элемент, для которого восходящая ветвь характеристики проходит ниже нисходящей (рис. 4.34, б), вносят отри-

цательный фазовый сдвиг.

Если нелинейный элемент однозначный, то из (4.85), (4.86) получаем:

q(A)

2

A

 

F(x)

 

xdx; b(A) 0.

 

A

 

 

A2

A2 x2

392

Если нелинейный элемент имеет однозначную нечетную статическую характеристику, т. е. F(-x)=-F(x), то выражение для q(A) принимает вид

q(A)

4

A

 

F(x)

 

xdx; b(A) 0.

(4.87)

A2

0

 

 

 

A2 x2

 

Например, для идеального реле (см. (4.8)) коэффициент гармонической линеаризации, полученный в соответствии с выражением (4.87), равен:

 

4

 

A

 

 

 

 

 

4

 

 

 

 

 

 

 

A

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

m

 

 

2

 

2

 

m

 

q(A)

 

 

 

 

 

 

 

 

xdx

 

(

A

 

x

 

)

 

 

,

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

2

 

 

 

 

 

A

0

 

A

 

x

 

 

 

A

 

 

 

 

 

 

0

A

 

 

 

 

 

 

 

 

 

 

 

 

что совпадает с (4.83), но не требует предварительного анализа периодиче-

ских режимов.

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

методом гармонического баланса

Аналитическое решение уравнений нелинейных систем практически удается достаточно редко, для несложных систем невысокого порядка с про-

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

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

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

жимов их функционирования.

Полагая эквивалентный комплексный коэффициент передачи нелиней-

ного элемента не зависящим от частоты, на основании (4.79) запишем (урав-

нение Гольдфарба)

W ( j ) 1

J(A)

или (уравнение Коченбургера)

J(A) 1

W ( j )

(4.88)

(4.89)

Из выражений (4.88) и (4.89) следует, что периодический режим в сис-

теме, содержащей нелинейный элемент, может иметь место только при необ393

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

ный элемент, будет обратно соответствующему преобразованию сигнала при прохождении его через линейную часть системы в том смысле, что отноше-

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

сдвиги по фазе - обратны по направлению.

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

стики W(j ) и 1 , то согласно выражению (4.88) автоколебания возможны

J(A)

только в точках пересечения этих характеристик, поскольку только в этих точках удовлетворяется условие (4.88). Так как каждой точке годографа

W(j ) соответствует определенная частота , а каждой точке годографа 1

J(A)

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

дографов определяют сразу значения амплитуды и частоты возможных в сис-

теме периодических режимов, как это показано, например, на рис. 4.35 /10/.

Выявление всех возможных видов бифуркаций (см. также 4.6), прису-

щих той или иной исследуемой нелинейной системе, вывод критических со-

отношений параметров, определяющих границы областей с различной каче-

ственной структурой фазового пространства,

в общем случае представляют весьма слож-

ную задачу. Как уже отмечалось ранее, для систем высокого порядка существующая тео-

рия не гарантирует полного решения этой за-

дачи. Для автономных же систем второго по-

рядка такая задача принципиально может быть всегда решена. Чаще всего это связано

Рис. 4.34. Амплитудно-фазовые со значительными трудностями, так как в ка-

характеристики линейной части и нелинейного элемента

394

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

дения качественного исследования, связанного с определением бифуркаци-

онных соотношений параметров, целесообразно применять вычислительные методы с использованием ЭВМ, различные приближенные методы по ап-

проксимации функций и т. д. Однако для того чтобы какой-либо периодиче-

ский режим был не только теоретически возможен, но и реально устанавли-

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

вым, т. е. чтобы ему в фазовом пространстве системы соответствовал устой-

чивый предельный цикл. В противном случае, если найденный периодиче-

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

Суждение об устойчивости периодического режима, соответствующего точкам пересечения годографов W(j ) и 1/J((A) может быть составлено на

основании следующих соображений.

Режиму работы системы, соответствующему точке пересечения годо-

графов амплитудно-фазовых характеристик согласно - выражению (4.79), со-

ответствует амплитудно-фазовая характеристика разомкнутой эквивалентной линейной системы, находящейся на границе устойчивости. Для выяснения вопроса устойчивости режима, соответствующего, например, точке N (рис. 4.35) проведем прямую 00через точку N1, лежащую на отметке амплитуды

АN1 > АN.

 

 

 

При этом, как видно из рис. 4.35, модуль вектора

 

1

оказался

 

ON1

J (AN 1 )

меньше модуля вектора ON2 W ( j N 2 ). Поскольку аргументы этих векторов одинаковы, то аргумент вектора W( j N 2 )J(AN1) будет равен нулю. Модуль

395

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

Проведем теперь прямую 00" через точку N3, лежащую на отметке ам-

плитуды АN3 < AN1. При этом, как видно из рис. 4.35, модуль вектора

 

1

окажется больше модуля вектора ON

4 W ( j N 4 ) и модуль произ-

ON3

 

 

J (AN 3 )

 

ведения W( j N 4)J(AN3) будет меньше единицы, что свидетельствует о затуха-

нии колебаний координат системы, работающей в этом режиме.

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

ние. Значит, периодический режим, соответствующий точке N, является не-

устойчивым.

Аналогичные рассуждения относительно точки М показывают, что со-

ответствующий ей периодический режим является устойчивым.

Следовательно, в системе, характеризующейся амплитудно-фазовыми характеристиками, приведенными на рис. 4.35, установятся автоколебания частоты M и амплитуды АM.

Из приведенных рассуждений видно, что точки пересечения годогра-

фов W(j ) и 1 здесь играют роль критических точек, аналогичных точке

J(A)

(1, j0) при исследовании устойчивости при помощи критерия Найквиста, от-

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

риодического режима.

Если точка годографа 1 , соответствующая несколько увеличенной

J(A)

амплитуде (по сравнению с амплитудой в точке пересечения годографов), не охватывается годографом W(j ) , то такой периодический режим является устойчивым и наоборот.

396

Е. П. Попов предложил использовать для определения параметров ав-

токолебаний идею Крылова и Боголюбова о гармонической линеаризации /5, 85/.

Рассматривая нелинейную замкнутую систему (рис. 4.36) с отрица-

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

Рис. 4.35. Структурная схема нелинейной системы с отрицательной обратной связью

A ( p ) x B ( p ) y

а уравнение нелинейного элемента, (4.76), –

 

b(A)

 

y J (A)x q(A)

 

p x .

 

 

 

(4.90)

на основании

(4.91)

Решая эти уравнения, получаем гармоническую линеаризацию уравнения системы:

 

 

b(A)

 

(4.92)

A(p) B(p) q(A)

 

p x 0.

 

 

 

 

 

Величина находится из условия, что существует периодическое ре-

шение уравнения (4.92), т. е. что j — корень характеристического уравне-

ния

A( ) B( ) q(A)

 

b(A)

 

 

 

 

 

0.

 

 

 

 

Но если j = λ, есть корень этого уравнения, то

A( j ) B( j ) q(A) jb(A) 0.

Выделим в A(j ) и B(j ) действительную и мнимую части, т.е.

A( j ) UA ( j ) jVA ( j ), B( j ) UB ( j ) jVB ( j ). Тогда

UA ( j ) UB ( j )q(A) VB ( j )b(A) j[VA ( j ) VB ( j )q(A) UB ( j )b(A)] 0.

Приравнивая нулю порознь вещественную и мнимую части, имеем

U A ( j ) UB ( j )q(A) VB ( j )b(A) 0,

(4.93)

VA ( j ) VB ( j )q(A) UB ( j )b(A) 0.

 

 

 

При однозначных нелинейностях b(A)=0 и из (4.93) следует

397

UA ( j ) UB ( j )q(A) 0,

(4.94)

VA ( j ) VB ( j )q(A) 0.

 

 

 

Из (4.94), исключая q(A), находим уравнение относительно :

VA ( j )UB ( j ) VB ( j )UA( j ) 0.

Определив , вычисляем А из одного из уравнений (4.94) выбрав из них наиболее простое.

4.5.Статистическая линеаризация нелинейных характеристик

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

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

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

нале попытаться линеаризовать существенно нелинейные характеристики та-

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

лезных сигналов и случайных возмущений через нелинейные звенья.

Метод приближенной замены нелинейных характеристик эквивалент-

ными в вероятностном смысле линейными зависимостями называется мето-

дом статистической линеаризации.

Рассмотрим совместное прохождение плавно изменяющегося полезно-

го сигнала и помехи через нелинейное звено. На рис. 4.37, а /4/ изображена характеристика нелинейного звена.

На рис. 4.37, б кривая 1 показывает изменение полезного сигнала, кото-

рый можно считать математическим

ожиданием входного случайного

Рис. 4.36. Прохождение случайного сигнала, а кривая 2 — изменение

сигнала через нелинейный элемент полного входного сигнала

mx(t)+X0(t)с учетом помехи (флук-

398

туации) входного сигнала. На рис. 4.37, в кривая 1 показывает изменение вы-

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

вая 2 — полный выходной сигнал при наличии на входе полезного сигнала и помехи. Пунктирной кривой 3 на рис. 4.37, в показано статистическое сред-

нее значение (математическое ожидание) выходного сигнала. Из этих рисун-

ков видно, что при малом уровне флуктуации входного сигнала кривая 3 бу-

дет мало отклоняться от кривой 1. При большом уровне флуктуации входно-

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

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

нейного звена. Чем выше уровень флуктуации, тем сильнее подавляется по-

лезный сигнал. При этом величина флуктуации на выходе рассмотренного нелинейного звена также ограничивается.

Идея статистической линеаризации состоит в том, чтобы заменить ис-

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

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

смотренную физическую картину совместного прохождения полезного сиг-

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

В результате статистической линеаризации зависимость между вход-

ной и выходной переменными нелинейного звена

Y (X)

заменяется приближенной зависимостью вида

U 0 k1X0.

(4.95)

399

mu my , Du Dy.
0 k0mx,

Величина 0 представляет собой выходной полезный сигнал нелиней-

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

статистической характеристикой нелинейного звена.

Коэффициент k1 представляет собой статистический коэффициент усиления нелинейного звена по случайной составляющей.

Для нелинейных звеньев, имеющих нечетные характеристики, функция

(ро может быть представлена в форме

(4.96)

где k0 статистический коэффициент усиления нелинейного звена по ма-

тематическому ожиданию.

Величина 0 и статистические коэффициенты усиления нелинейного звена k0 и k1 определяются из основного условия сохранения величины по-

лезного сигнала и уровня флуктуации на выходе нелинейного звена. Так как полезный сигнал представляет собой математическое ожидание, а уровень флуктуации характеризуется дисперсией выходного сигнала, то величины 0, k0 и k1 можно определить из условия равенства математического ожидания mu, и дисперсии Du, случайной величины U, определяемой формулой (4.95),

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

Y нелинейного звена:

(4.97)

Так как случайная величина U представляет собой линейную функцию случайной величины X, то на основании известных свойств математических ожиданий и дисперсий математическое ожидание и дисперсия случайной ве-

личины и определяются формулами

mu 0,

Du k12Dx.

(4.98)

На основании первой из этих формул первое условие (4.97) дает

0 my.

(4.99)

400

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]