Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 203
Скачиваний: 1
где Vd (со) — белый шум с интенсивностью, равной спектральной
плотности S°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ц (г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
и для Х°, |
7° |
ТД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