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