
книги из ГПНТБ / Пугачев, В. С. Основы статистической теории автоматических систем
.pdfSnW/BliJM |
|
|
SrMfil’pad |
|
|
|
W« |
|
|
|
|
. . |
300 |
|
|
L -500м |
|
|
|
жЛ_____—------- -----------—s |
|
|||
w-J |
5 10'! |
5 |
10-1 Q ,B 0 |
|
Рис. 3.4. |
Спектральная |
плотность |
Рис. 3.5. Спектралнзация плотности |
|
тангенциальной |
составляющей ско- |
нормальной составляющей скорости |
||
|
ростн ветра |
|
ветра |
в крупномасштабном и вязком интервалах. Несмотря на это, про стота аналитических выражений этих функций служит веским осно ванием для широкого использования их в практических расчетах. Обычно значения параметров в формулах для корреляционных функ
ций и спектральных плотностей соответственно |
равны: Lr = 200 -г- |
||
-еЗОО м; ом = |
2 -ьЗ м -с-1 — для ясной погоды; |
стш = |
8ч-12 м - с '1 — |
для кучевых |
облаков и <тш = 18 ч-25 м -с-1 — для |
грозовых ус |
|
ловий [22]. |
|
|
|
3.3.Элементарные звенья
Вданном параграфе рассматривается часто встречающаяся при расчетах задача определения вероятностных характеристик выход ных сигналов элементарных звеньев, перечень и характеристики которых приведены в приложении 1. На вход элементарных звеньев действует случайный процесс типа белого шума с нулевым математи ческим ожиданием и интенсивностью G. В качестве вероятностных характеристик выходных переменных звеньев вычисляют корреля
ционную функцию и дисперсию, а в установившемся режиме —соот ветствующую корреляционной функции спектральную плот ность.
Для вычисления используется метод весовых функций. В соответ ствии с этим методом корреляционная функция
t v |
|
Ку (t, ?) = } j g (t, т) g (/', т') к х (т, т') dx dx’. |
(3.36) |
о о |
|
Так как входной сигнал есть белый шум, то его корреляционная функция является S-функцией: Кх (х, т') = G6 (т— т'). Если под ставить это выражение в формулу (3.36), то можно воспользоваться свойством 6-функции (см. п. 1.4). Причем применение этого свойства
80
возможно как по одной, так и по другой переменной. При этом возни кает двухзначное решение:
t |
(3.37) |
Ku(l,t') = G\g(l,x)g(t>,x)dx; |
О
V
Ky (t, l') = G \g (t, Y)g{t',x')dx'.
о
Неоднозначность решения можно устранить, если учесть обяза тельное требование симметрии корреляционной функции относи тельно ее аргументов Ку (t, t') — Ку {t‘, t). При учете свойства симметрии формулы (3.37) можно объединить в одно выражение
min (I, V)
Ky ( t , n = G J |
g(t, x)g(f, x)dx, |
(3.38) |
о
где функция, стоящая в верхнем пределе интеграла,
min (t, t') = |
t |
t < r , |
V |
(3.39) |
|
|
t > t ' . |
Дисперсия выходного сигнала является корреляционной функцией при равных значениях аргументов. Спектральную плотность можно вычислить как преобразование Фурье от установившегося значения корреляционной функции или по формуле
Sy (ю) = |Ф (гсо) |2 5 0, |
(3.40) |
где Ф (i(o) — частотная характеристика звена; S 0 = G/2n —• уровень спектральной плотности белого шума.
При входном сигнале в виде белого шума элементарные звенья можно рассматривать как формирующие фильтры. Поэтому спек тральные плотности выходных сигналов звеньев могут служить неко торым набором типовых характеристик, по которым можно легко решать обратную задачу построения формирующего фильтра.
Усилительное звено. Входной и выходной сигналы усилительного звена прямо пропорциональны, т. е.
Y (t) = k X (t).
Весовая функция усилительного звена g (t, т) = k8 (t — x).
Подставляя весовую функцию в формулу (3.38), получаем
|
min (t , l') |
|
к у (t, V) = k2G |
J |
6 (t — т) 6 (t' — t ) dx = k2G8 (t — t'). (3.41) |
|
о |
|
Дисперсия выходного сигнала равна бесконечности, т. е.
Dy (t) = k8 G8 ( t — i) = oo.
6 В. С. Пугачев |
81 |
Спектральная плотность в соответствии с формулой (3.40) при
| CD (гео) |2 = /г2
(3.42)
и= *2Л
Из полученных формул следует, что выходной сигнал усилитель ного звена есть случайный процесс типа белого шума с увеличенной
в /г2 раз интенсивностью входного шума.
Интегрирующее звено. Входной и выходной сигналы интегри рующего звена связаны выражением
|
|
|
Y (0 |
= |
kX (/). |
|
|
Весовая функция интегрирующего звена является единичной |
|||||||
функцией: g (/, |
т) = |
/el |
(t — т). Подставляя весовую функцию в фор |
||||
мулу (3.38) и учитывая, |
что 1 (/ — т) = |
0 при t < т , |
получаем |
||||
|
|
|
min {/, |
/') |
|
|
|
|
Ку (/, |
V) = li2G |
J |
dx = |
k~Gmin (/, /')■ |
(3-43) |
|
|
|
|
|
о |
|
|
|
Дисперсия |
выходного сигнала |
|
|
|
|||
|
|
|
Dy (t) = k2Gt. |
(3.44) |
Спектральная плотность выходного сигнала обратно пропорцио нальна квадрату частоты:
(3-45)
Из формул для корреляционной функции и дисперсии следует, что выходной сигнал интегрирующего звена представляет собой неста ционарный случайный-процесс с линейно возрастающей во времени дисперсией. Интегратор является формирующим фильтром, который служит для получения случайного процесса с корреляционной функ цией (3.43).
Дифференцирующее звено. Входной и выходной сигналы диффе
ренцирующего звена связаны соотношением |
|
||
|
Y(t) = k**JfL. |
|
|
Весовая функция |
дифференцирующего звена |
есть производная |
|
б-функции g (t, т) = |
k8 (t — т). |
Корреляционная |
функция выход |
ного сигнала |
|
|
|
min (<, (') |
|
|
|
Ку (/, /') = k2G |
| |
•dd (^ 7 Т) dx = k2G'6 (t — t'). (3.46) |
|
|
о |
|
|
Дисперсия выходного сигнала Dy — оо. Спектральная плотность выходного сигнала пропорциональна квадрату частоты:
5»<“ >= -Т5Г- |
' (3-47> |
82
Звено запаздывания. Входной и выходной сигналы связаны соотношением
Y (t) = kX (t — т3),
где т3— время запаздывания. Весовой функцией звена чистого за паздывания является 5-функция, сдвинутая на время т3: g (t, т) =
— k8 ( t — т3 — т). Корреляционная функция |
выходного |
сигнала |
min ((, Г) |
|
|
K„(t,t') = k*G J б(/ — т3 — t)S(/' |
— т3 — x)dx. |
(3.48) |
о |
|
|
Данный интеграл в области t < Д ' не равен нулю, если выпол няются условия 0 •< t — т3 < t, 0 < V — т3 <Д t, из которых следует
т3 > |
0, t > т3. Аналогично для области t' < t имеем 0 <it — т3 •< |
, |
|||
О < |
V — т3 -< V. Из |
этих |
неравенств следует, что т3 > 0 , |
? > |
т3. |
Таким образом, интеграл (3.48) равен |
|
|
|||
|
| |
1гЮ8 (t — Г) t > т3, V > т3, т3 > 0, |
(3.49) |
||
|
1 |
0 |
t < т3, V < т 3, т3 > 0. |
||
|
|
|
|||
Равенство нулю корреляционной функции в области if, |
V < |
т 3 |
определяется свойством звена чистого запаздывания. При подаче на вход звена запаздывания в нулевой момент времени белого шума выходной сигнал будет равен нулю для моментов времени, меньших времени задержки.
Дисперсия выходной переменной звена запаздывания при дей ствии на входе белого шума равна бесконечности при t, t' > т3 и нулю при t, t' < т 3.
Апериодическое звено. Входной и выходной сигналы апериоди ческого звена связаны соотношением
T Y + Y = kX (t).
Весовая функция апериодического звена g (t, т) = ke~i-i~r^ T Подставляя это выражение в формулу (3.38), получаем
W ) = ^ |
min (I, Г) |
I + Г — 2т |
|||
1 |
е |
г_ |
dx = |
||
k°-G |
- (1 + |
Г) |
|
2 min {tf-t-) |
|
~ 2Т е |
( - 1 + е ‘ |
|
■ ) - |
||
|
|
« + /*> |
_ (/' —о |
|
|
k-G |
— e |
T + e |
r |
при t < V, |
|
2T |
|
(/ + Г) |
(l - i ') |
(3.50) |
|
|
|
^ ~ + |
e_ ~ |
7 |
при t > f . |
Эти выражения можно объединить, если учесть, что выражения в показателях вторых экспонент в обоих случаях не отрицательны. Поэтому, вводя знак модуля, получаем
ь-r |
r |
Г |
_ Д±Д> |
+ |
е |
\t —t’ \ -I |
г ...(3J.51)- |
|
K g { t , n = 4 |
|
L*— |
г |
|
63.
При значении переменных t + ? > 4Г первая экспонента близка к нулю, поэтому в установившемся режиме корреляционная функция
r |
\ i - i ' |
I |
K y[t,n = w * |
т |
■ |
Дисперсия выходной переменной |
|
|
Dy(i) = ^ { 1 |
- е |
(3.52) |
В установившемся режиме дисперсия имеет постоянное значение:
Dy = k2GI2T.
Спектральная плотность выходной переменной
Sy (со) = |
k"G |
1 |
(3.53) |
2л |
1-f- (шГ)3 |
Колебательное звено. Входной и выходной сигналы колебатель ного звена связаны уравнением
Т*У + 2\ТУ + Г = kX (/).
|
Весовая функция колебательного звена при {■с |
1 |
|
|
|||
|
. |
_ Е(С-т) |
I---- ta |
|
|
|
|
|
g(t, %) = т |
t* 6 |
Г |
s in — f M * - |
т)’ |
(3,54) |
|
где |
| — коэффициент затухания; |
Т — постоянная |
времени. |
|
|||
|
Часто используют другую форму уравнения колебательного |
||||||
звена: |
|
|
|
|
|
|
|
|
Y + |
2|со„У + |
моУ = щкХ (i), |
|
|
(3.55) |
|
где |
со0 = Ц Т — собственная частота |
звена. Весовая |
функция для |
этой формы записи определяется формулой (3.54), в которой следует заменить Т на 1/ со 0. Подставляя весовую функцию (3.54) в формулу (3.38), получаем
|
I и + п min (/, /') |
Ц-х . y i - t 2 |
|
|||
Ky(t, t') = |
k2Gt |
|
|
|
(t — т) X |
|
ТЦ1- | 2) |
|
|
|
: T sin -— |
||
|
|
|
|
|
|
|
|
X sin |
/1 |
- |
£a |
(,t' — x)dx. |
|
|
|
т |
|
|
Вычисляя этот интеграл и учитывая при этом, что — (t + t') + + 2min (t, t') = — 11— t' |, получаем
Ky(t, |
|
k-G |
|
т |
t') = |
4Г(1 — £a) U-7Г\e |
|
||
|
lU + t')\ |
— _ |
|
|
— e |
T |
) cos * |
|
{t — t') -f |
84
_ £(M-n , |
r ----- |
|
|
|||
+ e |
T |
(g COS |
г ё |
( t + |
1') — |
|
— V |
1 — l 2 sin ^ * ~ |
( t + |
7 ) ) — |
|||
_ |
r |
, |
-lAFZTp |
— |
|
|
— e |
(geos-— |
|
|
|||
■ V 1 — I 2 sin |
/ Г |
If — V |
(3.56) |
|||
|
В установившемся режиме при g (f + t') > 47 получаем
_ I I / - г |
| |
г |
|
/л — 71^ ^ |
|
cos — ^ |
(f — 7) -f |
М * . 0 = щ е |
|
||
sin / 1 |
- |
6s | / - Г | |
(3.57) |
При задании уравнения колебательного звена в форме выраже ния (3.55) корреляционную функцию удобно записать в следующем виде:
Ку (t, t') = Dy cos ©i (t —t')
-|---- sin ©! 11— V .—«| / —r |
(3.58) |
где Dy — дисперсия выходного сигнала; шх = co0 ) / l — g2 — резо
нансная частота корреляционной функции; а = g/71= g(o0Дисперсия выходного сигнала колебательного звена в соответствии
с формулой (3.56)
|
|
ш |
|
Du{t) |
k-G |
т |
1 + |
m |
|
||
l V \ |
— У sin ^ |
|
21 - |
— У cos |
21 |
(3.59) |
В установившемся режиме дисперсия имеет постоянное значение:
k-G
D„ (3.60) m
Для установившегося режима спектральная плотность выходной переменной
k -G |
'1 — 21<оГ)2 |
(1 — 2g2) - f (соТ)'1■ |
(3.61) |
|
S y i ® ) — 2л |
||||
|
85
Спектральную плотность часто удобно представить в форме выра жения
5 |
у(со) |
fegQT4 |
____________1_____________ |
|
(3.62) |
||
2п |
[а2+(ш —со1)-][а3+(ш + шд)2] |
’ |
|||||
|
|
|
|||||
где а = УТ\ |
сох = ± ш 0 Y |
1 — I2— частоты, соответствующие |
|||||
максимуму спектральной |
плотности. |
осциллятор). |
|||||
Колебательное звено |
без |
затухания (линейный |
Это звено является частным случаем колебательного звена при £ = 0. Поэтому все результаты можно получить из предыдущих формул, рассматривая предел при £ —>0. Вычисляя пределы, получаем корре ляционную функцию выходной координаты:
'min (/, l') |
COS |
(t - П |
Ky m = -kw . Т |
T |
— sin min (/, |
V) COS max (/, |
V) |
T |
T |
. |
где функция max (t, t') = t при t > t' и max (t, t') Дисперсия координаты
г, ... k-G l t |
. t |
t \ |
Д / (0 — 2T V T |
®^ f |
т J ' |
Спектральная плотность координаты
Sy(tо) = |
k-G |
|
9 |
9 \ 9 у |
|
2л |
(СО” — со0) |
(3.63)
t' при t' > t.
(3.64)
(3.65)
где со о = 1IT.
3.4. Колебательное звено
Колебательное звено является математической моделью большого числа различных устройств: измерителей, маятниковых устройств, радиотехнических цепей, простейших следящих систем и т. п. По этому представляет интерес определение вероятностных характери стик выходных переменных колебательного звена при различных типах входных случайных сигналов.
Рассмотрим уравнение колебательного звена в следующей форме:
У + 2Ео)07 + шоУ = щ Х (0, |
(3.66) |
где со о, £ — параметры звена; X (t) — случайный входной сигнал. Достаточно общей моделью входного сигнала является следующее
выражение:
х (/) = 5 (0 + N (0 + |
у„6(t - g + y 0s (t - g , |
(3.67) |
где 5 (t) — полезный сигнал |
е математическим ожиданием |
ms {t) |
и корреляционной функцией |
|
|
ks (т) = Д е ~ “ 1т 1^cos cosx — -^-sincoj т |
(3.68) |
86 ■ •
В формуле (3.67) N (t) — высокочастотная помеха с нулевым математическим ожиданием и корреляционной функцией вида
|
|
kN (т) = GN8 (т ). |
|
|
|
(3.69) |
||
Уравнение (3.66) интегрируется |
при |
нулевых |
начальных усло |
|||||
виях. Начальные условия представлены |
в виде |
эквивалентного |
||||||
входного сигнала — предпоследний и последний члены в |
выраже |
|||||||
нии (3.67). |
Начальная координата |
У0 |
и |
скорость |
У0 |
являются |
||
случайными |
величинами |
с математическими ожиданиями |
тУо, тУа |
|||||
и дисперсиями DUo, DlJo. |
Слагаемые в |
формуле |
(3.67) — некорре |
лированные между собой функции.
Задачей анализа является вычисление математических ожиданий и корреляционных моментов координаты и ее производной и опреде ление моментов ошибки работы звена. Под ошибкой системы пони
мается разность |
|
Е (0 = У — YT, |
(3.70) |
где Кт = A TS Ц) — требуемый выходной сигнал. Рассмотрим слу чай, когда требуемый оператор равен единице, а требуемый выходной сигнал равен полезному сигналу, т. е. рассмотрим следящую систему.
Математическое ожидание ошибки равно следующей разности:
mE{t) = my {l)— myr (t). |
(3.71) |
|
Дисперсия ошибки |
|
|
De (0 = Dy(t) + DUy (t) - |
20m (0, |
(3.72) |
где 0дат (t) — взаимный корреляционный |
момент выходного и тре |
буемого выходного сигналов, a DUt = Ds.
Математическое ожидание выходного сигнала колебательного
звена определяется уравнением |
|
ту + 2£сo0niy + сооту = соотх(/). |
(3.73) |
Математическое ожидание входного сигнала включает математи
ческое ожидание полезного сигнала ms (t) |
= а + bt и математиче |
|
ское ожидание начальных условий. Поэтому |
|
|
тх (t) = a-{- bt + тУо6 (t — 10) + |
тУо8 (t — 10). |
(3.74) |
Для вычисления математического ожидания выходного сигнала вместо уравнения (3.73) можно воспользоваться соотношением
t |
|
ту it) = J g (t, т) тх (г) dx, |
(3.75) |
^0 |
|
где t0 — момент начала работы системы; g (t, т) — весовая функ ция системы, определяемая решением при нулевых начальных усло^- виях уравнения (3.7^), в котором правая часть представляет собой
6-функцию:
f 21». + <4 g (t, т) = 8 (f - х). (3.76)
87
Решение этого уравнения эквивалентно решению однородного уравнения при начальных условиях: g (т, т) = 0; (dg (t, т)/dt)i=x — 1. Решая однородное уравнение при указанных начальных условиях, получаем следующее выражение для весовой функции:
g (/, т) |
|
- |
е |
‘ l L — |
sin со0 У 1 — |
— т). |
(3.77) |
||||||
|
|
|
«0 I |
1 — I" |
|
|
|
|
|
|
|
||
Подставляя весовую функцию и математическое ожидание |
вход |
||||||||||||
ного сигнала (3.74) |
в формулу |
(3.75), |
получаем |
|
|||||||||
Шу(i) = ---- г. 1 |
|
t |
|
|
|
|
|
|
|
|
|
||
|
( е - |
<' - т>sin со0)/ 1 — £2 (t — т) X |
|
||||||||||
“о У 1 — V |
,) |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
to |
|
|
|
|
|
|
|
|
|
X [coo (а + Ьт) -\- m,Ja 8 (т —t0) -f Шуа8 (т — /„)] dx. |
(3.78) |
||||||||||||
Вычисляя этот интеграл, представим математическое ожидание |
|||||||||||||
выходного сигнала в виде четырех слагаемых: |
|
|
|||||||||||
где |
|
ти(0 = |
myi + |
mlj2 + |
mUt + |
my„ |
(3.79) |
||||||
|
|
■SWo (/ — /о) |
|
|
|
|
|
|
|
||||
1 |
|
|
|
|
|
|
— £2 (( — /„) + |
|
|||||
|
|
1 |
|
|
-(l sin co0 У 1 |
|
|||||||
|
|
|
i - sa |
|
|
|
|
|
|
|
|||
+ |
/ |
1 — £2 cos co0 ]/ 1 — £2 (/ — g ) |
(3.80) |
||||||||||
mUi = |
гпУое |
- |
V ~ la) |
cos (00 ]/ 1— £2 (/ — t0) ~ |
|
||||||||
— |
У 1- |
|
sin co0 |
1 — £2 (/ — /„) |
(3.81) |
||||||||
|
|
ё“ |
|
|
|
|
|
|
|
|
|||
|
m. e- |
|
a - |
10) |
|
|
|
____ |
|
(3.82) |
|||
inyt = -■ |
|
--ГТ— В— |
sin “ о У 1— £2 V — t0); |
||||||||||
|
|
ш0/1 |
- £ 2 |
|
|
|
|
|
|
|
|
||
mUi (t) |
|
|
|
|
_ |
(e - |ш0 (/ —10) Ц (i — |
|
||||||
|
|
■t0) (| sin co0 У 1 — l2 (l — 10) -f |
|
||||||||||
+ |
V |
i |
^ |
f |
cos co0 y |
i |
^ |
f |
it - |
10)) + |
|
||
+ 2E V |
y |
y |
f cos co0 y |
r |
z |
I 2(t - |
10) - |
|
|||||
- (1 - 2£2)sin co0 У Г ^ | 2 (t _ |
/„)] - |
|
|||||||||||
- 2 1 У \ - ? ) + Ы 1 |
|
|
e - | w „ |
и - t o ) |
|
||||||||
|
|
|
|
x |
|
||||||||
|
|
|
|
|
|
|
|
|
|
Vi -6* |
|
||
|
|
X (|sin w0y |
i — l2{t — to) + |
|
|
||||||||
+ |
1^1 — E2 cos co0 ] / 1 — %2{t — iQ)) |
(3.83) |
88
При |
t — t о |
имеем |
myl — ту2 = тУ4 = 0, |
туг — туд, поэтому |
|
ту (t0) = |
myQ. В установившемся режиме в формулах (3.80)—(3.83) |
||||
следует |
положить t0 = —оо. В |
этом случае |
получаем ту, = а\ |
||
тУа = b (t '■— 2 |
со0); т,л |
= mUi = |
0. Поэтому |
математическое ожи |
|
дание выходного сигнала в установившемся режиме |
|||||
|
|
% со= а + b (* — fj-)- |
(3-84) |
||
Требуемое значение математического ожидания |
|||||
|
|
|
туу = |
а + Ы. |
(3.85) |
Следовательно, математическое ожидание ошибки в соответствии |
|||||
с формулой (3.71) представляет |
собой разность выражений (3.79) |
||||
и (3.85). |
|
|
|
|
|
Для вычисления дисперсии выходного сигнала воспользуемся принципом суперпозиции и условием некоррелированности входных сигналов. В данном случае удобно применить метод уравнений мо ментов (см. п. 2.6).
Представим уравнение (3.66) в форме двух уравнений первого
порядка. Производя |
замену |
переменных Y х = Y; Y 2 = Y, полу |
||
чаем |
|
|
|
|
|
|
^1 = П; |
(3.86) |
|
Y2 = ~ |
cogFj - |
2£со0У2 + со5Х (t). |
||
|
Учитывая, что помеха и полезный сигнал некоррелированы и система является линейной, рассмотрим прохождение одной помехи. В соответствии с уравнением (2.73) получаем следующую систему уравнений относительно корреляционных моментов выходных пере менных:
®и — 2612) |
|
|
012= Э22-- 0)0011 --2^0)0012! |
j |
(3.87) |
022 = --- 4|(Оо022 — 2о)о012 ~Ь COqGдг- |
) |
|
Эту систему уравнений следует решать при начальных условиях
0ц (0) = Dy„\ 022 (О) = Dyo; 012 (0) = 012о. Решая систему уравне ний (3.87) последовательно относительно 0Ш 012, 022, получаем сле дующие уравнения:
0U + |
6^соо0ц + |
4соо (1 -)- 2|')0ц |
8§соо0п = |
2cooGat; |
|
0i2 |
6^coo0i2 |
4соо (1 ~г 2§') 012 + 8^0)0012 = 0; |
(3.88) |
||
022 -f- 6|о)о022 + |
4о)о ( 1 -f- 2^')022 -f- 8^0)0022 = |
2cOoGj\r. |
|
89