Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 406

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
Рис. 11.1.1.
/

454

гл. 11. ИССЛЕДО ВА Н И Е

ТО ЧН О СТИ

Н Е Л И Н Е Й Н Ы Х СИСТЕМ

следует, что поведение этой

системы описывается уравнениями

 

ТУ + Y

=

к (X — U),

U = аУ\

Лішеаризуя второе уравнение путем разложения нелинейной функции в ряд Тейлора в окрестности точки ту , получим следующие уравнения для мате­ матического ожидания выходной переменной системы и центрированных случайных функций:

Tmy + my = k(mx —aml), r ^° + (1 +3/cam£) У0 = *Х0.

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

л

и дисперсию

выходной

переменной

______ у . системы

У можно будет определить

V * TS+J

любым

из

методов,

изложенных

 

в § 7.2.

 

 

 

Заменяя случайные функция и / их математическими ожидания­

ми в (11.1.3), получим

ту & фо (тх, о*). (11.1.6)

Таким образом, при применении метода статистической линеари­ зации сама'форма зависимости между математическими ожиданиями

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

Для получения уравнений для центрированных случайных функций следует вычесть равенство (11.1.6) из (11.1.3). Тогда получим

У° » к, (тх, <тя) X*.

(Н.1.7)

Добавляя к таким зависимостям для всех нелинейных звеньев системы уравнения входящих в данную систему линейных систем, получим полную систему уравнений, линейных относительно центрированных случайных функций. Коэффициенты уравнений вида (11.1.7) остаются неизвестными до тех пор, пока не определе­ ны математические ожидания и дисперсии входных сигналов всех нелинейных звеньев. Поэтому методы главы 7 не могут быть непосредственно применены для определения дисперсий пере­ менных.


§ 11.1. М ЕТО ДЫ Л И Н Е А РИ ЗА Ц И И

455

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

Для решения этой задачи в общем случае можно применить метод канонических разложений (§ 7.2). Применив этот метод к линеаризованным уравнениям системы, получим для коорди­ натных функций переменных линейные уравнения с коэффи­ циентами, зависящими от неизвестных математических ожиданий и дисперсий входных сигналов нелинейных звеньев. Добавив к уравнениям, определяющим математические ожидания и коор­ динатные функции, формулы вида (7.2.20) для дисперсий входных сигналов нелинейных звеньев, получим полную бесконечную систе­ му уравнений для математических ожиданий, координатных функций и дисперсий всех переменных, описывающих поведение исследуемой системы (координат системы). Ограничиваясь конеч­ ными отрезками канонических разложений, получим конечную систему уравнений для приближенного определения математи­ ческих ожиданий, координатных функций и дисперсий координат системы. Эта система уравнений может быть решена любым числен­ ным методом на цифровой вычислительной машине или с помощью моделирующих устройств.

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

П р и м е р 11.1.2. Рассмотрим систему, состоящую из последовательно соединенных релейного элемента и апериодического звена, охваченных отри­ цательной обратной связью (рис.

11.1.2). Уравнения этой системы

 

 

 

имеют вид

 

ГА

V

7

У .

7Т + У =■Zsgn

 

- 1

Ts +1

 

 

 

 

 

 

 

 

7=

*gQV, I

(11. 1.8)

 

 

 

.X

 

 

 

 

 

Рпс. 11.1.2.

 

Произведя статистическую линеа­

 

 

 

 

 

ризацию характеристики релейно­

 

 

 

го элемента,

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

переменных и для центрированных случайных функций:

 

 

Tmy + my= (f0(mB, а0),

тв^Шх — тпу,

(11.1.9)

TYO+ li + k^mv, CT„)]yo = Ä1(m0,

оѵ)Х<>,

уо = у о _ Х о .

(11.1.10)

Коэффициенты ф0 = к0тѵ и Ä, для релейного элемента определяются форму­ лами, приведенными в приложении 5 (нелинейность 5). Первое из этих урав­ нений не может быть проинтегрировано, пока не известна величина ст„. С другой стороны, аѵ нельзя определить из второго уравнения методами главы 7, пока не известен коэффициент /сь зависящий от неизвестных тѵ и а„.


456

ГЛ . 11. И ССЛЕДО ВА Н И Е ТО ЧН О С ТИ Н Е Л И Н Е Й Н Ы Х СИСТЕМ

Применив метод канонических разложений, получим для координатных функций переменных У и V уравнения

Ту\-Ь [1" Ь (mo> <УѴ)] yv= ki (mB1 0 „) xv, Ѵ\ —Уѵ

(v = l, 2, ...).

 

(11.1.11)

К уравнениям (11.1.9) и (11.1.11) следует добавить

формулу

о о

 

о « = 2 Л „ К Р ,

(11.1.12)

ѵ=1

 

определяющую дисперсию входного сигнала V нелинейного звена. Ограни­ чиваясь в канонических разложениях п первыми членами, получим систему 2п + 3 совместных уравнений (11.1.9), (11.1.11) и (11.1.12), приближенно определяющую 2ге+ 3 неизвестных функций времени m„, mv, уі, ѵь . . ., уп, ѵп , ое. После нахождения этих функций дисперсия выходной переменной У определится по формуле (7.2.20).

При применении изложенного метода к автоколебательным системам, т. е. к таким системам, в которых возможны устойчивые автоколебания, математические ожидания переменных получаются в виде колеблющихся функций времени, стремящихся к некоторым периодическим функциям. Если частота колебаний математических ожиданий сравнима со средней частотой колебаний случайных составляющих входных сигналов нелинейных звеньев, то изло­ женный метод может оказаться неточным. В подобных случаях можно получить более точные результаты, применяя совместную статистическую и гармоническую линеаризацию, как было изло­ жено в § 10.7.

§ 11.2. Применение метода статистической линеаризации для исследования точности стационарных систем

Рассмотрим сначала стационарную систему с одним нелиней­ ным звеном, работающую в установившемся режиме под действием

стационарного входного сигнала

X (<). Если

нелинейное звено

не находится в цепи обратной свя­

 

 

зи и не охватывается обратной

—►ф-

4>(s)

связью (рис. 11.2.1), то дисперсия

0 ( s )

 

- J 3

Рис. 11.2.1.

Рис. 11.2.2.

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


§ 11.2. П РИ М Е Н Е Н И Е М ЕТОДА СТА ТИ СТИ ЧЕС КО Й Л И Н Е А Р И ЗА Ц И И 457

чинами (см. [54], § 3.10 или [53], § 30). После этого можно будет найти спектральную плотность выходного сигнала нелинейного звена и определить дисперсию выходной переменной системы Y по формуле (7.4.4).

Если нелинейное звено находится в цепи обратной связи (рис. 11.2.2) или охватывается обратной связью (рис. 11.2.3 и рис. 11.2.4), то точное определение дисперсии выходной перемен­ ной системы в общем случае невозможно и для ее нахождения

Рис. 11.2.3. Рис. 11.2.4.

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

Обозначая выходной сигнал нелинейного звена через U, при­ меняя статистическую линеаризацию нелинейного звена и заме­ няя в (11.1.3) функцию фо ее выражением (8.5.3): ф0 = коту7 напишем для случая нелинейного звена в цепи обратной связи уравнения (11.1.6) и (11.1.7) в виде

ти — коту, U° = kiY°,

(11.2.1)

где ко и кі являются известными функциями величин ту и оу. Формулы (11.2.1) показывают, что если считать к0 и неизвестны­ ми постоянными, то в результате статистической линеаризации нелинейное звено заменяется безынерционным усилителем с раз­ ными коэффициентами усиления полезного сигнала и помехи. Таким образом, при постоянных к0 и статистическая линеариза­ ция заменяет данную систему двумя разными линейными система­ ми для полезного сигнала и для флуктуаций. Пользуясь формулой (4.6.7), находим передаточные функции линеаризованной системы

для полезного сигнала и для

флуктуаций:

 

ФГ 00 =

ф(»)

ф '1’ (s)-

Ф(»)

( 11. 2. 2)

і + к0ф (S)

і + кіФ («)

где Ф (s) — передаточная функция линейной части системы. Считая передаточную функцию Фу0>(s) известной, применим

формулу (7.3.8) для математического ожидания выходной пере­ менной системы. Тогда, принимая во внимание, что математическое ожидание стационарной входной переменной системы X (t)

постоянно, получим

_______ Ф (0) тх

( 1 1 . 2 . 3 )

 

У 1 + ^о і т у і °у) Ф (0) ’


458

гл. 11. И ССЛЕДО ВА Н И Е ТО ЧН О СТИ Н Е Л И Н Е Й Н Ы Х СИСТЕМ

где мы теперь показываем в явной форме, что коэффициент к0 является известной функцией величин ту и ау. Уравнение (11.2.3) содержит две неизвестные величины ту т оу ѵі поэтому может быть решено только совместно с уравнением, определяющим дисперсию выходной переменной системы <з\.

Для определения дисперсии выходной переменной системы У применим формулу (7.4.4). В результате получим

2

______Ф (т)______

do,

(11.2.4)

оV

1 + ^ (ту, су) Ф (ісо)

где также показана явно зависимость к4 от ту и ау.

Решив уравнения (11.2.3) и (11.2.4) относительно неизвестных величин ту и ау, мы и найдем математическое ожидание и диспер­ сию выходной переменной рассматриваемой системы. Эти уравне­ ния могут быть решены любым приближенным методом. В частно­ сти, можно применить метод последовательных приближений или графический метод.

Для решения уравнений (11.2.3) и (11.2.4) методом последова­ тельных приближений необходимо задать исходные грубо при­ ближенные значения к0и kt и вычислить в первом приближении тпу и ау по формулам (11.2.3) и (11.2.4). После этого можно уточнить значения к0 и kt и вычислить тпу и ау во втором приближении. И так далее. Процесс заканчивается, когда два последовательных приближения совпадут в пределах принятой точности вычислений. Обычно практически бывает достаточно двух или трех прибли­ жений.

Для решения уравнений (11.2.3) и (11.2.4) графическим мето­ дом заменим уравнение (11.2.3) равноценной системой уравнений

Ч =

Л 1 +

ф (0) тх

(11.2.5)

о^ФІО)