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

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

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

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

Добавлен: 16.10.2024

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

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

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

где

F (Л, h) = Ц аг (/г) А/,;

Г=1

m

Я (А, /г) = 2 6Г(/г) А*;

Г=1

А/, — разностный оператор г -г о порядка.

Для применения метода необходимо представить случайное воз­ мущение X (h) разложением по случайным параметрам подобно вы­ ражению (2.51):

N

 

X (Л) = тк(h) + Ц v,x, (/г).

(6.46)

/=1

 

Кроме того, должны быть заданы

начальные

математические ожи­

дания и дисперсии для переменной

Y (0) и /г —

1 разностей.

 

Как и в п. 2.5, искомую переменную Y отыскиваем в форме

К Сh, в) = ту (/г, е) + £ У,у, (А, е).

(6.47)

 

l=i

 

 

Применяя принцип суперпозиции и используя выкладки, ана­ логичные приведенным в п. 2.5 для определения ту (Л, е) и коорди­ натных функций ijj (h, е). получаем уравнения

F (A, h) ту (h, е)

=

Я (A,

/i) m* (Я);

(6.48)

К (А, Л) у/ (h, е) =

Я (А, Л) х,- {К).

(6.49)

(/ = 1, •

Я)

 

 

 

Уравнения (6.48) и (6.49) аналогичны исходному уравнению (6.45).

Следовательно, для определения

функций

ту (h, е)

и y-t (h,

е),

входящих в формулу (6.47), в этом случае также необходимо N +

1

раз решить исходную систему уравнений при определенных началь­ ных условиях.

Уравнение (6.48) следует интегрировать при заданных начальных условиях для переменной ту (0) и ее п — 1 разности. При интегри­ ровании уравнения (6.49) начальные условия могут быть учтены

аналогично тому, как это сделано в п.

2.5.

В результате

уравнение

(6.49) для определения функции

у х (/г,

е)

необходимо интегрировать

при ненулевых начальных условиях

 

 

 

Re А'У1(0) = У

;

im Агу, (0) = 0,

(6.50)

= 0 , 1, . . ., п 1)

 

где Dj = М [К?] — дисперсия первого случайного коэффициента разложения (6.46); Ядг^(0)— дисперсия г-ой разности функции Y

в

начальный момент. При интегрировании

остальных уравнений

(/

ф 1) начальные условия следует полагать

равными нулю.

12*

179.


После решения разностных уравнений (6.48) и (6.49) определяют т.и (/г, е) и г/у- (/г, е) для любого момента времени t — (h + в) Тп,

h = 0, 1, . . ., О < в <; 1.

Корреляционную функцию переменной Y (/г, в) вычисляют по следующей формуле:

N N

Ку ( К И,

1и, в2) =

£

£

KiVyj(hi,

&i)yv (h2, е2),

(6.51)

 

 

/=1 V=1

 

 

 

 

где K/v = М [ K / V v

1. При /г2 =

/г15

е2 = B j

и з

формулы (6.51) по­

лучаем дисперсию переменной

Y (А, е).

 

X (/г) представлена

В частном случае, если случайная функция

каноническим разложением,

из выражения

(6.51) получаем

 

Ку (hi, еь 1ц, е2)

=

N

 

_

 

(6.52)

У, Djtjj (Ль ех) г/„ (/г2, е2) ,

где Dj = М [К/].

Если для рассматриваемой дискретной системы задана некоторая идеальная линейная операция преобразования полезного сигнала

tnx (li) вида

 

Ут(h, е) = L (А, /г, в) пгх (/г),

(6.53)

то ошибка в системе определяется формулой (2.62). Для математи­ ческого ожидания и дисперсии ошибки дискретной системы спра­ ведливы формулы, рассмотренные в п. 2.5 и записанные для дискрет­ ного аргумента.

Так решается основная задача оценки точности дискретной одно­ мерной линейной системы при использовании метода интегрирова­ ния уравнений состояния.

Рассмотрим нелинейные системы. Уравнение одномерной дискрет­ ной нелинейной системы, аналогичное выражению (4.40), имеет вид

F (А, К) Y (Л, в) + R (А, К) <p IY (/г, в) ] = Я (A, h) X (/г), (6.54)

где ф — нелинейная характеристика.

Для уравнения (2.39) должны быть заданы начальные условия такие же, как для уравнения (6.45).

Случайное возмущение необходимо представить в виде выраже­

ния (2.33),

а решение для

Y (/г, в) отыскивается в форме

 

 

N

Y (h,

в) = my (h, в) +

К0 (/г, в); К0 (/г, в) = £ V,y, (А, в). (6.55)

 

 

/=1

Для применения излагаемого метода нелинейность ф необходимо линеаризовать любым способом (см. п. 1.8). Рассмотрим более общий случай статистической линеаризации нелинейности для дискретной переменной:

Ф [Y (h, в)) = ф0 + k iY 0 (h, в).

(6.56)

Подставляя выражение (6.56) в уравнение (6.54) и пользуясь принципом суперпозиции, с учетом формулы (6.55) получим уравне­

180



ния для математического ожидания и координатных функций пере­ менной:

F (А, /г) 1пу (/г, е) + R (А, /г) ср0 = Я (А, /г) mA.(/г); (6.57)

Т7 (А, /г) г/у (/г, е) + 7? (A, h) kxlJj (/г, е) = Я (А, /г) Xj (Л). (6.58)

(/ = 1, • ■ Я)

Уравнения (6.58) для координатных функций любого номера оди­ наковы, но связаны между собой и с уравнением (6.57), что обуслов­ лено зависимостью ср0 и 1гх от ту и Dy. Поэтому полученные раз­ ностные уравнения необходимо интегрировать совместно, дополнив формулами

Фо =

Фо (h, е),

Dy (/г,

в)];

k x =

[m„ (/г, е),

(Л,

в)];

 

n

_

(6.59)

А ,=

£ Яг/г/£(Л, s) г/, (/г, в).

 

/=1

 

 

При интегрировании разностных уравнений (6.57) начальные условия должны быть заданы в виде математических ожиданий пере­ менной У и ее п — 1 разностей для начального момента времени. Начальные условия для уравнений (6.58) задаются и выбираются в соответствии с формулами (6.50).

Изложенный метод исследования точности .нестационарных диск­ ретных систем обобщается на многомерные системы и нелинейные системы со многими нелинейностями аналогично тому, как это сде­ лано в п. 2.5 и 4.4 для непрерывных систем.

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

Кдискретным линейным системам можно применить метод со­ ставления и решения уравнений корреляционных моментов. Для этого их уравнения должны быть приведены к форме выражения (1.18). Рассмотрим для упрощения обозначений дискретный выход

(в = °)

АУ, (А) = £ aki 0h) Уг (/г) + Ьк (/г) Vk (Л),

(6.60)

/= 1

(k = 1, • ■., п)

где Vk (/г) — дискретный белый шум, имеющий отличное от нуля математическое ожидание. Применяя операцию математического ожидания к уравнению (6.60), получим систему разностных уравне­ ний для определения математических ожиданий всех переменных

Атук (A) = £

aki (h) У; (А) + bk (К) V,. (/г).

(6.61)

t=1

 

(k

= 1, . . ., /г)

 

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

181


Путем почленного вычитания уравнений (6.61) из уравнений (6.60) получим систему разностных уравнений для центрированных случайных составляющих

AYl (A) = £ aki (/г) У? (/г) + bk (А) Vl (А).

(6.62)

1= 1

(/г = 1, . . ., п)

Пользуясь уравнениями (6.62), составим разностные уравнения для определения корреляционных моментов переменных. Как и для непрерывных систем, обозначим через

0// (А) = М [У? (Л) У/ (Л)]

корреляционный момент переменных У? (/г.) и У9 (/г) и запишем формулу первой разности для этого момента

А0,у (Л) = М [ДУ° (/г) У/ (/г) + АУ? (/г) У? (Л) +

+ АУ? (Л) АУ/ (Л)]•

(6.63)

(', / = К ■• »)

Используя уравнения (6.62) и подставляя выражения для А У? (/г) в формулу (6.63), получим

Д0,-у (А) = S [а/Г (А) в{/ (/г) + а/г (/г) 0Г/ (/г)] +

/•=1

+ bt (А) М [У? (/г) У? (/г)] + Ь,- (А) М [У? (/г) У? (А)] +

+ Ъ, (Л) М [У? (А) У/ (А)] + £

(А) а/р(А) 0гр (А) +

г, р=1

 

+bi (h ) £ air(h)M\Y°r (h)V0i(h)} +

Г= 1

+ 6, (Л) £ аip(li) М [Ур (А) У? (А)].

(6.64)

p = i

 

Значения выходных переменных физически возможной дискретной многомерной системы выражаются через входные переменные фор­ мулой

y°t(h) = t £ Ян (А, г) У? (г),

(6.65)

;= 1 г=о

где |f.; (/г, /-) — весовые коэффициенты дискретной линейной системы, соответствующие i-му выходу и I-му входу.

Вычислим корреляционные моменты М [ У9 (А) У9(А)], исполь­ зуя формулу (6.65):

М [У? (А) У/ (А)] = f

(А, г) М [У? (А) У/ (А)],

(6.66)

;= 1 г—0

(i, / = 1, . . . . /г).

182