Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

где функции

Фо = Фо(|Н„ . D i/.’

ki =

k\ (ту> D yt, D yx)\

^2 ^

^2 (^//> ^У\ r

^!/i )

имеют конкретный вид в зависимости от типа нелинейности; Dyi — дисперсия; D • — дисперсия производной функции Y х (I) в устано­

вившемся

режиме.

 

и (4.34) в уравнение (4.31)

Подставим выражения (4.32), (4.33)

и учтем,

что ср0 =

k 0my для нечетной

нелинейности. В результате

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

 

F (р) ту = Я (р) [тх — 1г0ту]\

 

(4.35)

 

 

N

sin соrt — kxYxItipYj

(4.36)

F(p)Y1 = H(p)

Х°

Из уравнения (4.35) для установившегося режима получаем

т„

Я(0)

 

 

 

(4.37)

F (0) + каН (0) т■л-;

К ~ К ( ту> Dy i ,

D'y ,).

Полигармонический

сигнал

N

можно

рассматривать

°r s'n

 

 

 

Г=1

[27 ]

 

 

как процесс со спектральной плотностью

 

 

 

 

N

 

 

 

 

 

S(co) =

 

(О— со,.) -f- 6 (со + (О,)],

 

где б (со ± сог) — дельта-функция.

 

 

 

Дисперсия Dyl и D

в установившемся режиме в рассматривае­

мом случае вычисляют по формулам линейной теории (п. 2.4), при­ нимающим вид

 

 

 

та

 

 

N

^

 

 

 

 

 

DUl =

J

I ф ( К О )

|* Sx (со) dco +

^

IФ (Ч ) I* ±

 

;

(4.38)

 

 

та

— СО

 

Г = 1

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

D-t = Г | Ф (too) I2 Sx (со) da +

V

| Ф (toor) |2 со?

?

(4.39)

где

 

 

 

 

 

 

 

 

 

 

СВ(ко) = - ■,. . .

Н

|

,---- частотная характеристика

ста-

v

'

F (ш ) +

(ко) [kx +

/г2ш]

 

r

г

 

 

тистически линеаризованной системы от входа Z до выхода Y. Выражения (4.37), (4.38), (4.39) являются уравнениями, которые

связывают величины ту, Dyi, Ь • . Для их определения необходимо эти уравнения решить совместно, учитывая формулы для /е0, k x,

К-

124


Случайный процесс Y (t) в установившемся режиме имеет вид

N

у (0 = ту + У0 (/) + £ cr IФ (ког) | sin (соrt -фг);

Г=1

ФЛ= arg® (ко,.).

Спектральная плотность этого случайного процесса

N

Sy И = I Ф ( И Г S* (со) + ^

(иог) |2 [6 (со — еог) + 6 (со + сог)].

г=1

Из последней формулы следует, что спектральная плотность со­ стоит из непрерывной и дискретной частей. Непрерывная часть может иметь экстремум иа частоте, определяемой из уравнения

dSy (со) _ а

дш ~

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

Аналогично решается задача при наличии в системе нескольких нелинейностей. Заменив нелинейности линейными безынерцион­ ными зависимостями с соответствующими статистическими коэффи­ циентами усиления и выразив при помощи формул типа выраже­ ний (4.23), (4.24) или (4.28) математические ожидания и дисперсии входных сигналов всех нелинейностей, получим необходимые урав­ нения для определения этих неизвестных. Для решения уравнений в общем случае следует применить метод последовательных прибли­ жений. В результате решения этих уравнений определяются мате­ матические ожидания и дисперсии входных сигналов нелинейностей и соответствующие статистические коэффициенты усиления. После этого по формулам типа (4.25), (4.26), (4.27), (4.29) определяют мате­ матические ожидания и дисперсии выходных переменных и ошибок системы.

4.4. Метод интегрирования уравнений системы

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

Рассмотрим существо и особенности применения метода интегри­ рования уравнений так же, как и в п. 2.5, предварительно для одно­ мерной системы с одной нелинейностью. Пусть уравнение такой

динамической нестационарной системы имеет вид

 

F (t, р) Y + R (t, р) ф (У) = Н (t, р) X ,

(4.40)

125.


где

F (t> P) =

S ar ( 0 Pr Я (*.

P) = 2 j

cr ( 0 pr;

 

r—0

г—0

 

 

m

 

 

 

н (t, p) = S

( 0 P r .

 

 

r=0

 

 

Ф (^) —• нелинейность

произвольного

вида;

X (/) — случайный

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

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

У (0) п производных

У(г) (0) до я-го порядка:

 

 

ту (0), . . . , /и]?-15

(0), Dy (0), . . . ,

£>,(!.-» (0).

Случайную функцию Y (t)

представим в виде выражения Y (t) =

= niy (/) + У0 (I). Пусть ср (У) является гладкой,дифференцируемой нелинейностью и допустима обычная ее линеаризация:

ср (У) = ср (niy) + ср' (ту) У0.

(4.41)

Подставляя выражение (4.41) в уравнение (4.40) и пользуясь принципом суперпозиции, из линеаризованного уравнения полу­

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

и

центрированной

составляющей переменной:

 

 

 

 

 

F (t,

р) ту +

R

(t,

р) ср (niy) = Н (/,

р)

тх\

(4.42)

F (t, р)

Y° -г R

(t,

р)

ср' (яд У° = Н (t,

р) Х \

(4.43)

где тх — математическое ожидание;

X й— центрированная составляющая случайной функции X (t). Уравнение для математического ожидания (4.41) является нели­ нейным и интегрируется при заданных начальных условиях для ти (0) и я — 1 ее производных. После его интегрирования определяем ту и ср' ( n i y ) . Уравнение (4.42) — линейное с переменными параметрами. Для его интегрирования и определения Dy полностью применим

метод, изложенный в и. 2.5.

В общем случае нелинейность ср (У) всегда может быть линеари­

зована статистически:

 

ф (Y) = ср0 (niy, Dy) + /г 1 (niy, Dy) У0.

(4.44)

Подставляя выражение (4.44) в уравнение (4.40) и пользуясь прин­ ципом суперпозиции, из линеаризованного уравнения получим урав­ нения для математического ожидания и случайной центрированной

составляющей

переменной

У:

 

F (t,

р) ту +

R

(t,

р) Фо (ту, Dy) = Н (t, р) тх\

(4,45)

F (t,

р)

Y° +

R

(/,

р) к, (ту, Dy) Y° = Н (t, р) Х \

(4.46)

где тх (t) — математическое ожидание; Х° (t) — центрированная со­ ставляющая случайной функции X (t).

126


Уравнение (4.45) следует интегрировать при заданных начальных условиях ту (0), . . , т ^ п_Х) (0). Для уравнения (4.46) заданы

начальные дисперсии переменной Y и ее производных У(г) до п— 1-го порядка. Это уравнение может быть проинтегрировано методом канонических разложений. Для этого представим случайную функ­ цию Х° (/) каким-либо разложением:

 

х° (0 = S

V,X, (О,

(4.47)

 

/=i

 

 

где Xj

(t) — координатные функции; V,- — случайные величины с ну­

левым

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

и корреляционной

матрицей

Ru = M \ V lV,\.

Решение уравнения (4.46) также будем изыскивать в виде разло­

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

 

(t) = 2 Vitji ((),

(4.48)

/=i

 

где tjj (t) — неизвестные координатные функции разложения. Подставив выражения (4.47) и (4.48) в уравнение (4.46) и выпол­

нив необходимые преобразования, для определения координатной функции tjj (/) получим

F (t, р) tjj + R (t, р) k x (ту,

Dy) tjj = Я (t, р) Xj.

(4.49)

( / = 1 , . . . .

ло

 

Уравнения (4.49) для координатных функций связаны между собой и с уравнением (4.46) через функции ср0 и k x от ту и DtJ. По­ этому полученные уравнения необходимо интегрировать совместно, дополнив следующими формулами для ср0, k x и Dy:

Фо = Фо (ту, Dy); k 1 = k 1 (ту, Dy)\

(4.50)

Du= Ъ RijVi ( 00/ ( 0-

t. i = i

Уравнения (4.49) необходимо интегрировать при определенных начальных условиях для координатных функций. Однако в рассмат­ риваемом случае заданы начальные дисперсии переменной Y и ее п — 1-й производной, как и в задаче, рассмотренной в п. 2.5. Поэтому для определения начальных значений координатных функций вос­ пользуемся темн же рассуждениями и выкладками. В результате начальные условия для первой координатной функции у х при раз­ ложении по действительным координатным функциям имеют вид

где D, = R lt = M [Vi\.

127


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

Re г/{г) (0) =

/ и Лг)

0.

 

у

Imr/{r) (0) =

 

(г = 0, 1, . . ., п — 1)

 

 

Начальные условия для

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

других номеров

(/ ф 1) должны быть взяты равными нулю.

 

анализа точ­

Таким образом, подставленная

основная задача

ности нелинейной нестационарной

системы сводится

к совместному

интегрированию нелинейного уравнения (4.45) для математического ожидания, N линейных уравнений (4.49) для координатных функций при определенных начальных условиях с учетом конечных соот­ ношений (4.50). Особенность применения рассмотренного метода к нелинейному исходному уравнению состоит в необходимости совмест­

ного интегрирования полученных уравнений в отличие от

линейного

исходного уравнения, рассмотренного в п.

2.5, или нелинейного, но

допускающего обычную линеаризацию. При интегрировании пере­

численных уравнений определяются ту (t)

координатные

функции

Hi (0> /= 1, . . ., N и дисперсия Dy (t).

Одновременно

с интегри­

рованием уравнений может быть вычислена корреляционная функ­ ция по формуле (2.61) или (2.62).

Если для рассматриваемой системы задана некоторая идеальная операция над полезным входным сигналом тх (t), то по формуле (2.34) могут быть рассчитаны ошибки в системе.

Изложенный метод исследования нелинейных систем применим также к многомерным системам со многими нелинейностями.

Пусть многомерная нелинейная система задана уравнениями вида

Уi =

Ф, {t, Y lt

Уя)

+

bt (t) Л'„

(4.51)

 

( t= 1..........Я)

 

 

 

где ф£— нелинейные функции произвольного вида;

статистикой,

Х[ — случайные

функции

времени

с

заданной

которые можно представить в виде какого-либо разложе­

ния по случайным параметрам:

 

 

 

 

 

 

Ni

 

(0,

 

X £ (t) = /П,. (t) +

 

 

 

(/= 1 ,

. . . . q)

 

(4-52)

где Vtj — случайные коэффициенты, имеющие равные нулю матема­ тические ожидания, заданные дисперсии и корреляцион­ ные моменты связи;

Хц (I) — известные координатные функции.

128