Файл: Методы оптимизации в статистических задачах управления..pdf

при нормальном распределении вектора х (t) оказываются свя­ занными функциональными зависимостями

Совместное интегрирование уравнений (50) и (56) с использо­ ванием зависимостей (68) может быть выполнено только с помощью ЦВМ. Причем, однократным решением этой системы определяются математическое ожидание и дисперсионная матрица вектора фазовых координат систем (32). Помимо этого, в процессе реше­ ния этих уравнений определяется матрица эквивалентных коэф­ фициентов усиления векторного нелинейного элемента Кг, которая вместе с матрицей Dx может быть исползована для решения диф­ ференциального уравнения (67). При этом однократным решением дифференциального уравнения (67) определяется сечение корре­ ляционной матрицы вектора х (t) как функция Ц при фиксирован­ ном значении t2. Общее число решений уравнения (67) для опре­

в обыкновенное дифференциальное уравнение относительно корре­ ляционной матрицы вектора фазовых координат

(69)

которое решается при начальных условиях, определяемых си­

Решение системы уравнений (70) может быть выполнено на ЦВМ методом последовательных приближений. В процессе решения этой системы уравнений определяются Кг и Dx, которые исполь­ зуются для решения уравнения (69).

Рассмотрим пример составления дифференциальных уравнений для определения элементов матрицы корреляционных функций выходных координат системы.

Пример. Нелинейная система, изображенная на рис. 3, возмущается в мо­ мент времени t = 0 стационарным белым шумом |( 0 с интенсивностью Q,

21

равной единице, И неслучайным воздействием г (t). Нелинейная система опи­ сывается дифференциальными уравнениями

//у ]

-Ф (Xx, x2)];

(71)

dx2

dt = 41 -2[*i(0-

где cp {xl7 x2) — двумерная нелинейная функция выходных координат системы, имеющая вид:

Предполагая совместный двумерный закон распределения хг (t) и х2 (t) близким к нормальному, построим систему уравнений для определения элементов матрицы корреляционных функций выходных координат системы (71) в переход­ ном режиме.

Принимая во внимание формулу (67) и используя выражения (71) и (72) при двумерном нормальном векторе' л: (t), получим уравнения относительно эле-

Рис. 3. Структурная схема нелинейной системы, содер­ жащей двумерный нелиней­ ный элемент

ментов матрицы корреляционных функций выходных координат системы в об­ ласти, где tx > t2:

dkxx Wh

22

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

3. Интегральный метод анализа точности нелинейных систем

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

Рассмотрим многомерную нелинейную систему порядка п (рис. 4). Обозначим через х (t)

вектор фазовых координат си­ стемы размерности т на входе нелинейного элемента cp (х, t) той же размерности. Введем ма­ трицу размерности \т, т] им­ пульсных переходных функций W (t, т) линейной части систе­ мы. Тогда вектор х (t) выразит­ ся следующим образом:

где В (т) — матрица переменных коэффициентов размерности [т, /]; z (т) — нормально распределенный вектор случайных воз­ действий размерности / с известными статистическими характе­ ристиками. Начальные условия вектора х (t) предполагаются нулевыми.

Принимая нормальным двумерный закон распределения век­ тора X (t), определим приближенно математическое ожидание и матрицу корреляционных функций этого вектора.

Используя формулу (73), определим математическое ожида-

О

ние тх (t) и центрированный вектор л: (t) через матрицу импульс­ ных переходных функций линейных звеньев:

23

Корреляционную функцию (78) векторного нелинейного эле­ мента ф (х, t) при двумерном нормальном законе распределения вектора х (t) можно представить в виде суммы двух составляющих:

где Ч*- (т, X) — известная функция, которая для определенного вида нелинейного элемента при двумерном нормальном распре­ делении вектора на входе выражается через математическое ожи­ дание и элементы матрицы корреляционных функций этого вектора.

24

Первое слагаемое равенства (79) представляет собой корреля­ ционную функцию векторного нелинейного преобразования, про­ порциональную корреляционной функции вектора х (t), которая получена с помощью статистической линеаризации нелинеййого элемента.

Второе слагаемое этого равенства выражает нелинейные иска­ жения вида корреляционной функции вектора х (t). Это слагаемое представляется разложением в ряд по степеням элементов матрицы корреляционной функции вектора х (t), начиная со вторых сте­ пеней. При этом используется разложение 2«-мерного нормаль­ ного закона распределения по ортогональным полиномам Чебы­ шева—Эрмита. В частном случае, когда нелинейный элемент имеет полиноминальные характеристики, функция 47 может быть определена непосредственно.

В интегральное уравнение (77) входят матрицы коэффициентов

статистической линеаризации Кі (т) и Кі (А) векторного нелиней­ ного элемента по случайной составляющей, а также функция ТР (т, А), которые при двумерном нормальном законе распределе­ ния вектора х (/) являются функциями математического ожидания и матрицы корреляционных функций этого вектора, т. е.

В силу этого интегральное уравнение (77) и уравнение (74) необходимо решать совместно с зависимостями (80).

Систему уравнений (74) и (77) можно решить методом последо­ вательных приближений или проинтегрировать методом квадра­ турных формул на цифровых вычислительных машинах. Исполь­ зуя метод прямоугольников, представим систему уравнений (74) и (77) в следующем виде [2]:

где А — шаг интегрирования по переменным т и А.

25