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

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

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

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

Добавлен: 16.10.2024

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

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

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

где Vd (со) — белый шум с интенсивностью, равной спектральной

плотности n (со) стационарной случайной последовательности

№ ( h — l).

Подставляя выражение (6.23) в формулу (6.22) и меняя порядок суммирования и интегрирования, получим

Я

т

У0 (h, е) =

J

Vd (со) Ф (ко, е) eUo/lT"dco,

(6.24)

 

Я

 

 

 

~1тм~

 

 

где введено обозначение

 

 

 

ф ((‘со,

е) =

цу (/, е) е ш,тп-

(6.24)

 

 

1=0

 

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

П

П

Dy = J

j Ф о, е) Ф (—/со", е) k Vd(со — со') е,А<(0~ш,) dcoс/со'.

(6.25)

Но корреляционная функция белого шума Vd (со) имеет вид

kVd(со — со’) = SdN б (со — со ).

Подставляя это выражение для корреляционной функции в фор­ мулу (6.25), получим

Я

Тп

 

Dy = J |Ф (/со, е) f Sn (со) dco.

(6.26)

я

 

Для вычисления интеграла в формуле (6.26) выполним некоторые

преобразования. Прежде

всего

перейдем к переменной

z = еги7п

в формуле (6.26):

 

 

 

Dy = 4 -

ф | Т

(г, е) |2 vdN (z) 4 - dz,

(6.27)

| г |=1

 

 

где учтена формула (1.60) и введено обозначение

 

y%(z) = ± S dN ( - ~ \ n z ) .

(6.28)

174


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

ций комплексного переменного формуле

1+ W

Z = -1г-1—W .

При этом окружность единичного радиуса, по которой осуществ­ ляется интегрирование по z, превращается в мнимую ось, так как

при z =^е'шГп переменная

w —

 

 

 

sC

Я )

 

а формула (6.27) принимает вид

 

 

 

 

D U =

-

/со

9NКИ

8)Г dw,

 

f

J I

(6.29)

где

 

 

 

 

 

 

(6.30)

V (w, е) = ¥

 

е) ,

pdN(ш) = vdN ( - Щ ) ■

Заменив w = tg

в

выражении

(6.29),

окончательно получим

 

 

СО

 

 

 

 

 

Dy =

J

|Y*(i£, е

) ^ ®

^

d\.

(6.31)

Функции ¥* и рм являются дробно-рациональными относительно переменной £, а их коэффициенты — действительными числами. В этом случае подынтегральное выражение в формуле (6.31) может быть представлено в следующем виде:

 

gk (iI)

 

(6.32)

 

hk № hk(—fg)

 

 

где

 

 

 

 

(г1) — ао (%)k +

ai (*’?)* 1 +

*• • +

ak\

 

gk № = b0 т - к- г +

b, (ig)2*-4+

• • • +

bk. 1.

 

В результате формула (6.26) приводится к известным уже выра­

жениям (2.46):.

 

 

gk (tj)

d£.

(6.33)

D.. = 2л/Му ' . = Т г 1 hk (i£) A* (— *‘6)

Эти интегралы вычисляют с помощью таблиц, приведенных в при­ ложении 2.

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

175


Перейдем к анализу нелинейной дискретной системы. Рассмотрим одномерную стационарную дискретно-непрерывную систему, имею­ щую дискретные части с передаточными функциями г1гд1 (г), Ч'д2 (г), непрерывную часть с передаточной функцией гРи (г) и одну безынер­ ционную нечетную нелинейность ср в непрерывной части, как изо­ бражено на рис. 6.2. Входной сигнал системы U (t) имеет полезную неслучайную составляющую ти (t) и случайную стационарную ад­ дитивную помеху N (/). Для дискретных моментов времени t = hTn, h = 0 , 1 , . . . входной сигнал

U (/г) = ти (/г) + N (/г).

Пусть идеальная операция L, характеризующая желаемый вы­ ходной сигнал уТ (/г, е), в общем случае для непрерывного выхода имеет вид выражения (6.12). Мгновенная ошибка определяется фор­ мулой (6.13).

Для анализа дискретной нелинейной системы применим метод статистической линеаризации. В рассматриваемом случае на основа­ нии формул, приведенных в п. 1.8, нечетную нелинейность предста­

вим в следующем виде:

 

Ф W = ka (tnx, Dx) + kx(mx, Dx) X \

(6.34)

После такой замены нелинейности рассматриваемая дискретная система формально может быть описана линейными разностными урав­ нениями и охарактеризована двумя структурными схемами по отно­ шению к среднему и случайному сигналам (рис. 6.3). В установив­ шемся режиме в стационарной системе, имеющей необходимый поря­ док астатизма, в дискретные моменты времени величины пгх и Dx являются постоянными. В этом случае рассматриваемая дискретная система характеризуется стационарными весовыми и передаточными функциями. Пользуясь соответствующими структурными схемами (рис. 6.3), определяем передаточные функции дискретных систем при е = 0 от входа U до двух выходов X и Y для математического ожидания и флуктуирующей центрированной составляющей. Эти передаточные функцииимеют вид соответственно: для пгх и ту

Ч'л-о (2) =

1+

Тд1 (г) Удг (3 У (г) k0 (тх, Dx) ’

(6 ‘35)

ш / . Ч _

! +

ФД1 (г) Ф„ (г) к0 (тх, Рх)

(6.36)

у0 Кг) ~

Тдг (2) Фд2(2) Фн (г) ('«*. Dx) '

 

Рис. 6.2. Нелинейная дискретно-непрерывная система

176


и для Х°,

ТД1

 

 

 

 

(*) = !+ *

 

(6.37)

 

 

Д 1 (г) '^ Д 2 (2)

(г) k1 (Wjj, ОдД.)

 

 

 

Тдг (2)

Т„ (г)

Д ’

(6.38)

 

 

 

 

 

где z =

sГ„

 

 

 

eJ

 

 

 

 

 

Для

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

д. и ту

получаем

формулы,

аналогичные выражениям (6.17):

 

 

CD

>пх = 2

^ Т Г

I1" " 1'0 (/г) 4

о .

Цо.1->

 

г=О

 

 

 

 

 

 

 

 

 

0

где

г

 

 

 

 

 

 

 

 

 

о

^Чо(е5?~п)

.0

/Уио (е Гп)

Цг.с

dsr

s= (f •И-гу =

dsr

s=0

Пусть сг — коэффициенты

ошибок,

которые

в данном

вычисляют по формулам

 

 

 

 

(6.39)

(6.40)

случае

сг = (

7 f ^Lr//

ctr ;

(г =

0, 1, ....

л)

 

 

(6.41)

Сг = {— 1)гу г Л

(г =

л + 1 , ...)•

Ф

Рис. 6.3. Структурные схемы линеаризованной

системы:

а — для математического ожидания; б — для случайной составляющей

12 В. С. Пугачев

177


Математическое ожидание ошибки тЕ выражается зависимостью, аналогичной формуле (6.18):

СО

mB = crm(r) (h) -f /ПдгЦо^-

r=0

Дисперсии Dx, Dy = DE определяются формулами, ными формуле (6.31):

00

 

2

4 ;

{ h r;, ((g)iV v№ )

1 +

I 2

— со

(6.42)

аналогич­

(6.43)

00

 

 

2

dl,

(6.44)

\ I 'IV Г PN № 1+ fe2

— со

 

 

где

 

 

Ч*х1(11) = у х1( ± ± ^ у у

 

 

- Щ г );

p A f ( ' i ) = - f ^ S V ( l r 7 l n 1 — i \ ) ’

Sjv — спектральная плотность дискретного стационарного случай­ ного процесса №.

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

Формулы (6.39) и (6.43) являются уравнениями, связывающими

величины тх и Dx. Присоединяя к ним выражения для k 0

{тх, Dx)

и Aj (rnx, Dx) и решая совместно методом последовательных

прибли­

жений или графически, определяем значения тх, Dx,

k 0, kv После

этого по

формулам (6.39), (6.40), (6.44)

вычисляют

величины ту,

т Е, Dy,

De .

дискретных

стационарных

Изложенный метод расчета точности

нелинейных систем обобщается на более сложные системы, имеющие несколько нелинейностей. При этом, как и для непрерывных систем (см. п. 2.4), объем вычислений возрастает пропорционально числу нелинейностей.

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

Рассмотренный ранее метод интегрирования уравнений применим также к дискретным нестационарным линейным системам.

Одномерная дискретная линейная нестационарная система харак­

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

 

F (А, /г) Г (/г, е) — Я (A, h) X (/г),

(6.45)

178