Файл: Пугачев, В. С. Основы статистической теории автоматических систем.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