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

книги из ГПНТБ / Пугачев, В. С. Основы статистической теории автоматических систем

.pdf
Скачиваний:
36
Добавлен:
22.10.2023
Размер:
12.95 Mб
Скачать

SnW/BliJM

 

 

SrMfil’pad

 

 

 

 

 

. .

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)

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') = — 11t' |, получаем

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 ©! 11V .—«| / —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

Соседние файлы в папке книги из ГПНТБ