Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.pdf

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

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

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

Добавлен: 10.04.2024

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

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

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

§ 9 . 2 ]

ОП РЕДЕЛЕН И Е СИСТЕМЫ ПРИ БЕСКОНЕЧНОМ ВРЕМЕНИ

477

§ 9.2. Определение оптимальной динамической системы при бесконечном времени наблюдения

Как было указано в § 9.1, в задаче Винера на вход динамической системы поступает сумма управляющего воздействия (полезного сигнала) U (t) и возмущающего воздействия (помехи) V (t), кото­ рые являются стационарными случайными функциями времени (будем считать й=г;=0). Рассматриваемая система должна на своем выходе воспроизводить некоторый требуемый сигнал Z (t), который связан с U (t) заданным функциональным соотношением

Z(t) = NU(t),

(9.5)

где N — некоторый преобразующий оператор.

осуществимой,

Оптимальная система должна быть физически

т. е. ее реакция на единичную импульсную функцию должна обра­ титься в нуль для if 0. Следовательно, импульсная переходная функция I (t) (весовая функция) должна удовлетворять условию

l(t) — 0 при t < 0.

(9.6)

Погрешность е (t) системы определяется разностью между действительным выходным сигналом Y (t) и требуемым выходным сигналом Z (t), т. е.

s(t) = Y (t) Z (t).

(9.7)

В связи с этим условие минимума средней квадратической ошибки может быть записано в виде

М [е2 (f)] = М {[F (і) — Z (f)]2} = min.

(9.8)

Задача синтеза оптимальной системы считается решенной, если определена импульсная переходная функция 1(х) или связан­ ная с ней взаимно однозначными соотношениями передаточная функция системы L (іш):

СО

L (іш) = J I (t) e~imidt,

 

 

°

 

 

(9-9)

 

 

 

eimiL (іш) d<n.

 

 

 

 

—CO

 

 

 

Спектральная

плотность

входной

случайной функции

X (t)

в соответствии с

(1 )

определяется

соотношением

 

s x И =

S U («>) +

S V (®) +

S

« , (№) + Sna(а»),

(9.10)

где Su (ш) и (о>) — спектральные плотности случайных процес­ сов U (f )n F (t) соответственно, а Suv ( ш) — их взаимная спектраль­ ная плотность,


47 8

ВЫ БОР ОПТИМАЛЬНЫХ СХЕМ

И ПАРАМЕТРОВ ГУ

[ГЛ. 9

 

Для взаимной спектральной плотности Szx ( ш) имеем

 

 

S zx И = SzuН +

S„ (со).

(9.11)

Если спектральную плотность Sx ( ш) можно разбить на произ­ ведение двух функций С (іш) и С* со), из которых первая имеет полюсы только в верхней полуплоскости комплексного аргумента <о, а вторая — только в нижней полуплоскости, т. е. можно на­ писать

IС (іш)I2 = С (т) С* (іш) = Sx (ш),

(9.12)

то передаточную функцию оптимальной динамической системы можно представить в виде (формула Н. Винера)

ьи

-

(JU

^

ш

)

S e ~*w dt

2

Г

или

 

 

В (іш)

 

 

 

 

L (Іш) :

 

 

 

 

С (іш)

 

 

где

В (іш) = j ß (і) e~iwtdt,

о

$гх (м) eiu,tdu>.

С*(іш)

—СО

(9.13)S й

(9.14)

(9.15)

(9.16)

В том случае, когда Sx (ш) является дробно-рациональной функцией ш, для разбиения ее на множители (1 2 ) достаточно по­ ложить

 

 

С /.Л _ ,.2 Рт(“)

Рт(“)

(9.17)

 

 

о„ (ші —. а

—с•

 

,•,

 

 

е

 

 

т

 

РтИ

=

11 (10

Qn(®) =

11 ("> — ѵг)";>

 

 

j= 1

 

 

/ = 1

(9.18)

 

 

 

 

 

 

 

р

 

 

п

 

к . п

=

п > -

<?; и

=

о» -

 

 

j=г

 

 

/=і

 

a=const, комплексные корни ц . и ѵг имеют положительные мни­ мые части, а т . и п1 ~ кратности соответствующих корней. Для

дробно-рациональных спектральных плотностей интеграл (13) может быть вычислен в общем виде.


§ 9 . 2 1 ОП РЕДЕЛЕНИ Е СИСТЕМЫ ПРИ БЕСКОНЕЧНОМ ВРЕМЕНИ

479

Если Z (<) является результатом применения оператора N к текущим или будущим значениям ординат функции U (t), то I65]

/у (/«>) :

1

(?„(<■>)

1 И .

(9.19)

а'1

Рт Н

 

а

Іг

 

(9.20)

 

 

 

 

г=1к=1

d lr~k

(9.21)

(lr-k)'- dJr~

 

где \ r — полюс кратности lr

Ql (<•>)

выражения p*'(-;)) SX!! (со), лежащий

в верхней полуплоскости.

Если функция Z (t) является результатом применения опера­ тора N к ординатам функции U (t) в прошлые моменты времени, то

 

 

 

 

 

 

а2

^ п Н

 

 

 

(9.22)

 

 

 

 

 

Іг

 

 

 

 

 

 

 

 

 

аг

 

 

 

 

 

 

 

 

 

*F («

 

_

ск

_

 

 

 

(9.22)

 

 

 

 

 

г—Л/.—] ѵ

 

 

 

 

 

 

 

5*г — (/; _ к) !

Д*--* г,

 

:«(<•

 

 

 

(9.24)

 

уг_к[(“

КУГp ^ (J) s » Н ]

о = Х -

 

где

.,

 

кратности

,,

 

 

QZM с

/ \

лежащий

Аг — полюс

Zr выражения

"

' охг (ш),

в нижней полуплоскости.

 

 

 

 

являющаяся мини­

 

Дисперсия ошибки оптимальной системы,

мальной для

принятых исходных данных,

определяется соотноше­

нием

 

 

 

 

 

 

 

 

 

 

 

M[s2 (7)lmin =

D[Z(0 ] - D ( F ( 0 ] =

 

 

 

 

 

 

 

 

 

 

 

с о

 

 

с о

 

 

 

 

 

 

 

 

=

J 5,(«o)Äo—

J I /у (іш)(2 ^ (ш)

(9.25)

или

 

 

 

 

 

 

 

 

 

 

(9.26)

 

 

М

[® 2

(7)]щіп=

0

 

 

2 ( 0

dt,

 

 

 

 

Кг ( )

 

2 л

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

ß (t) определяется формулой

(16).

 

 

когда управ­

 

Обобщение приведенных результатов на случай,

ляющее и возмущающее воздействия приложены к различным точ­ кам системы, дано в [47], [49]. При этом рассматривается случай, когда кроме полезного сигнала U (t) и помехи V (t) на систему действует возмущение, точка приложения которого не совпадает с точками приложения U (t) ж V (t).


480

ВЫ БОР

ОПТИМАЛЬНЫХ

СХКМ И ПАРАМЕТРОВ

ГУ

[ГЛ. 9

 

§ 9.3. Определение оптимальной динамической системы

 

 

при конечном времени наблюдения

 

 

 

Излагаемый

ниже метод оптимизации динамических систем

отличается от

рассмотренного

в предыдущем параграфе

двумя

существенными

особенностями:

во-первых, время

наблюдения

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

ляющей, содержит также составляющую в

виде полинома g (t)

от времени, т. е. математическое

ожидание

полезного сигнала

g (t) может быть задано в виде многочлена от времени

г (0 = 2

Ѵ г-

(9-27)

#= 0

 

 

В зависимости от имеющихся сведений о составляющей g (t) следует различать три случая: а) о коэффициентах кд ничего не известно; б) значения коэффициентов кд заданы; в) коэффициенты kq являются случайными и известны их вероятностные харак­ теристики.

Таким образом, в отличие от рассмотренного к § 9.2 случая Винера входная функция системы X (t) содержит неслучайное слагаемое g (t), т. е. выражается формулой

X ( t ) = g ( t ) + U ( t ) + V ( t ) ,

(9.28)

где [g (t) -\-U ([)] и F (t) по-прежнему обозначают полезный сигнал и помеху, а функции U (t) и V (t) будем считать стационарными, несвязанными случайными функциями с нулевыми математиче­ скими ожиданиями.

На выходе динамической системы получаем некоторую вели­ чину Y (t). Для импульсной переходной функции системы I (<) вместо одного условия (6 ) для рассматриваемой системы с конеч­ ной «памятью» должны выполняться следующие два условия:

I (t) — 0

при

t О 0,

(9.29)

I (t) = 0

при

t^>T.

(9.30)

Первое из них, аналогичное (6 ), является условием физической осуществимости системы; второе — отражает сделанное предполо­ жение о невосприимчивости системы к сигналам, поступающим в систему ранее чем за время Т до настоящего момента времени (конечная «память» системы).

Входной сигнал X (t) должен быть преобразован в некоторый выходной сигнал Y (t), который должен, в определенном смысле