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

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

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

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

Добавлен: 16.10.2024

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

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

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

X 1 —

Д/ + T \ 2 f t + 1

 

 

2т3k

 

/ Д< + Т \ 2 * + 2 '

+

Т + т

 

 

(Д/ + т)2

V Г + X J

+

 

r*(k+\)

 

 

L1

( - r

i f f *

3])-

<з л м >

(2/г + 3) (Л/ +

т)3

Функция

 

 

т—м

 

 

 

 

 

 

 

 

 

 

 

 

X

М *. т, г , Д0 =

0

0

0

 

 

 

S')

 

 

0

 

бф

 

 

 

X Аа£> (£) D (£') 6 (s -

g) б (s' -

 

 

s') ds ds' d ^ d f .

Г) д ^ }- S (s -

 

 

 

 

 

 

 

 

 

 

(3.195)

В результате вычисления этого интеграла получим

 

h = k2

(А - О2

 

/

М + т \2*+1

 

(2А+ 1) (Д< + х)

\

Г + т

 

 

 

 

t'J

 

 

 

х(А-

1)

 

_

/

Д<^ + Тт

\?*+2у

+

 

 

 

(Д/ + т)2

 

 

\

7 4 -т )

 

 

 

 

 

 

 

 

 

т2(* +

I)2

 

 

Д/ + т \2*+з

(3.196)

+ (2А + 3) (Д/ + т)2

 

Г + т ) ’

JJ

 

 

В вышеприведенных формулах при А/ < т необходимо прини­ мать А/ = т. Это объясняется тем, что полет не может быть управ­ ляемым в течение времени т, меньшего времени инерционности дви­ жения относительно центра массы.

Суммированием квадрата математического ожидания (3.181) и дисперсии (3.191) определяется второй начальный момент промаха. Исследование зависимости второго начального момента промаха от обобщенного коэффициента усиления показывает, что существует экстремальная точка, обеспечивающая минимум среднего квадрата промаха (с^)т1п при оптимальном коэффициенте усиления.

=

Рассмотрим числовой пример при следующих данных

[50, 61];

Gg = 0, бф =

6,28-10~6

рад2-с;

m

=

1; а

= 2 ;

п =

j

ig\ g = 9,81 м/с2;

m

= m . = 0:

ст

= 100

м;

<т. =

10

«и

пц

 

850

ч

^

=

500

м -с"1:

v

«о

Уо

м -с"1: v =

м-с

 

v

=

350

м -с"1;

D„ = 5000

 

 

 

14,28

с;

DB= 75

м,

ч

DB/vr = 0,21

с;

х =

0,2 с;

м; Г = Da/vr =

Д/ =

/ т = е т = 0 Т = 0 ц Т = 0 .

графики

второго

 

начального

момента

промаха

ракеты

 

На рис.

3.12 показаны

 

в зависимости от обобщенного коэффициента усиления к, условий стрельбы и пара­ метров ракеты.

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


Г л а в а 4

МЕТОДЫ АНАЛИЗА

 

НЕЛИНЕЙНЫХ НЕПРЕРЫВНЫХ СИСТЕМ

4.1. Задачи анализа нелинейных систем

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

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

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

Задачи анализа нелинейных динамических систем, находящихся под воздействием случайных возмущений, состоят в определении за­ конов распределения или вероятностных числовых характеристик всех или только выходных переменных, в исследовании зависимости этих характеристик от параметров системы, а также числовых ве­ личин входных сигналов, в определении точности воспроизведения или преобразования полезных сигналов. Одной из основных задач среди перечисленных является задача оценки математических ожи­ даний и корреляционных функций выходных переменных. Эта задача включает также оценки математического ожидания и . корреляцион­ ной функции ошибки воспроизведения системой полезного или задан-

Ш


ного сигнала, т. е. основную задачу точности. Общая постановка этой задачи аналогична постановке, сформулированной в гл. 2, п. 1 для линейных систем. Однако решение задачи для нелинейных замкну­ тых систем связано со значительными трудностями вычислитель­ ного характера. При этом для произвольных нелинейных систем не существует общего точного метода определения интересующих нас вероятностных характеристик. Поэтому большое значение имеют приближенные методы решения задач анализа, в том числе и задач определения точности нелинейных систем, основанные на применении обычной и статистической линеаризации функциональных зависимо­ стей. Если функциональные зависимости линеаризованы, то прибли­ женно анализ нелинейных систем может быть проведен методами ли­ нейной теории. К приближенным методам относится также метод эк­ вивалентных возмущений, основанный на аппроксимации выходных переменных любой динамической системы полиномами по случайным параметрам, если входные случайные возмущения представлены в па­ раметрическом виде через случайные параметры [23, 31].

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

4.2. Преобразование случайных функций нелинейностями

Простейшей нелинейной системой является безынерционная не­ линейность. Значение функции Y на выходе безынерционной нели­ нейности ср определяется только значением функции X на входе в тот же момент времени t:

У (0 = Ф [X (*)].

(4.1)

Нелинейное безынерционное преобразование случайного процесса, так же,как и линейное пропорциональное, не вносит дополнительных вероятностных связей. Это означает то, что, если процесс на входе безынерционного преобразования полностью характеризуется «-мер­ ным распределением, то процесс после преобразования (на выходе нелинейности) также характеризуется распределением /г-го порядка. Например, если случайный процесс на входе в нелинейность пред­ ставляет собой белый шум, то на выходе он также является белым шумом. Однако в отличие от линейного безынерционного преобразо­ вания при нелинейном преобразовании изменяется вид закона рас­ пределения.

Задача анализа применительно к заданному преобразованию (4.1) состоит в определении закона распределения функции Y (t) при заданном законе распределения функции X (t). Для опреде­ ления закона распределения функции Y(t) необходимо существова­ ние обратного функционального преобразования

X (t) = ф (Y (*)).

(4.2)

112


Если обратное преобразование (4.2) однозначно и задача функции плотности вероятности f (х, t) распределения процесса X (t), то плот­ ность вероятности к {у, t) процесса Y (t) определяется на основании

следующих рассуждений. Вероятность того,

что функция X (t) за­

ключена в интервале х 0, х 0 + dx,

равна / (х,

t) dx. Вероятность того,

что У (t) заключена в интервале у 0, у 0 +

dy,

равна к (у, t) dy. В силу

однозначной детерминированной зависимости между Y и X вероят­

ность к (у, t) dy = / (х, I) dx.

 

 

 

Разделив правую и левую части этого равенства на dy, получим

к (у, t) =

f(x, 0 - ^ - .

(4.3)

Из формулы (4.3) следует, что необходимо существование произ­

водной от обратной функции ~

^ .

Кроме того, так как плот­

ность вероятности не может быть отрицательной, в формуле (4.3) необходимо учитывать модуль производной. Окончательно формула для плотности вероятности выходной переменной принимает вид

h{y, = t) ~dy (4.4)

Если обратное преобразование неоднозначно, то в формуле (4.4) необходимо осуществить суммирование по каждой ветви неоднознач­ ной характеристики ф (у).

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

мулу

[45,

56,

68]

 

 

^ (d ll

• • -1

Уги

^1> ■■• I

k

■• • •) %nkI ^1> ■• ■>tn) JП> (4-5)

 

 

 

 

 

где xik = ф* (xit . .

т _ d (xik.......xnk)

j n

Уп)

d {‘Jv

yn) — ветвь неоднозначного преобразования; - якобиан обратного преобразования вида

dipi

дУп

дфп

дуг дуп

суммирование осуществляется по всем ветвям обратного преобразо­ вания. При этом, если определяют m-мерный закон распределения на

выходе к {уг, . . ., ут, tlt . . ., tm) и т < « ,

то в

формуле (4.5)

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

а для

определения

к (уи . . ., у,п, tlt . . ., t,n) выражение (4.5) необходимо

проинтег­

рировать по переменным хт+1, . . ., хп. Если

т >

п,

то т п

переменных у{ функционально связаны с остальными,

а функция рас­

пределения

/Н—Л

 

^ (Уи • • •> Упи ^1> • ••> ^т)

{^.Уи • • ч Уги

^ Уп+l® ( Ф/)-

 

1 = \

 

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

113


Вычисление плотностей вероятности по приведенным формулам достаточно просто для однозначных обратных преобразований. Однако на практике часто обратное преобразование многозначно пли в некоторых случаях имеет бесконечное число значений, а ее произ­ водные терпят разрыв. В этих случаях для определения одномерных функций распределения необходимы специальные приемы, а опреде­ ление многомерных функций распределения становится весьма гро­ моздким [24].

Практически валено, однако, определение одномерных моментных функций любого порядка, а также второго двумерного центрального момента — корреляционной функции на выходе нелинейности. Эта задача более простая для определения одномерных моментных функ­ ций, причем для определения одномерного момента любого порядка случайной функции на выходе произвольной нелинейности необ­ ходимо знать закон распределения входной функции:

СО

М {[К(*)]*} =

J

[ф (*)]*/(*. t)dx.

(4.6)

 

— СО

 

 

k =

1,

2.......

 

В частном случае полиноминальной характеристики ср (х) формула для k-ro момента переменной Y (t) выражает этот момент через мо­ менты переменной X (t) для того лее момента времени и может быть получена из формулы (4.6), если в нее подставить выралеение

 

 

 

N

 

 

Ф(*) = 2

 

В результате получаем

i=i

 

 

 

м

[КА] =

S

Иг • -с1кМ [Х''1+-+Ч

(4.7)

где

‘Т..... £*=1

 

 

 

 

 

 

 

 

 

со

 

 

Л4 [Хг1+-"+П] =

j х'4+-”+'*f(x, t)dx

 

 

 

 

— со

 

есть гх + • . . +

//.-начальный момент случайной функции

X (t)

для момента времени t. Из формулы (4.7) следует, что момент /е-го порядка выходной случайной функции выралеается линейно через высшие моменты входного процесса, начиная с момента /е-го порядка до момента kN-vo порядка.

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

114