Файл: Пугачев, В. С. Основы статистической теории автоматических систем.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 — kxYx— ItipYj |
(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) Uо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* (со) + ^ |
1Ф (иог) |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