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

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

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

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

Добавлен: 16.10.2024

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

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

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

образований входного сигнала, например для экстраполяции, диф­ ференцирования и т. д. Для экстраполятора Ут (/) = 5 (t + Т), а для дифференциатора Кт (/) = S' (i). В общем случае линейного преобразования требуемый выходной сигнал выражается формулой

(10.45)

Г-Л

 

где Uь . . ., UN — те же случайные величины,

что и в выражении

(10.27) для входного полезного сигнала, а ф] (t),

. . ., фд, (t) — неко­

торые известные функции. Если форма сигнала известна, то N — 1. Функции фг (t) выражаются через соответствующие функции срг (/) различными формулами в зависимости от конкретной задачи. Так,

например, для экстраполятора фг (i) ~

cpr (i

-f

Т) (г =

1, . . ., N),

а для дифференциатора фг (/) == ср) (/)

(г =

1,

. . ., N).

В других

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

S

(i) = sin соД + U.2 cos со0t,

где со0—-частота сканирования. В данном случае

ср: (t) = sin сo0i, ср. (t) = cos со0t.

Требуемый выходной сигнал этой схемы представляет собой вели­ чину в одном канале и U2 в другом канале, так как эти величины пропорциональны отклонениям луча локатора от самолета в двух ка­ налах. Поэтому в данном случае фх (t) — 1, ф2 U) = 0 для одного канала и фх (/) = 0, ф2 (t) = 1 для другого канала.

Задача определения оптимальной системы для заданного преобра­ зования полезного сигнала решается совершенно так же, как задача фильтрации. Оптимальная система и в этом случае преставляет собой параллельное соединение согласованных фильтров, соответствующих различным компонентам полезного входного сигнала, с подключен­ ными к их выходам усилителями. Однако коэффициенты усиления выходных сигналов согласованных фильтров в случае преобразова­ ния сигнала будут другими, так как в выражении ошибки системы вместо функций cpr (t) будут фигурировать функции фг (t). Вследствие этого уравнения для определения оптимальных значений коэффициен­ тов усиления выходных сигналов согласованных фильтров будут отличаться от соответствующих уравнений, рассмотренных в пре­ дыдущих двух параграфах, тем, что вместо функций cpr (t) в них будут функции фг (t) .Точно так же все формулы для ошибок системы будут отличаться от соответствующих формул преды­ дущих двух параграфов, тем, что вместо функций cpr (t) в них будут функции фг (t). Таким образом, для определения коэффи­ циентов усиления выходных сигналов согласованных фильтров и для нахождения ошибок оптимальной системы преобразования си­ гнала следует в формулах (10.12), (10.13), (10.16), (10.17) заменить

271


функцию ср (t) функцией ф (i), а в формулах (10.33), (10.35), (10.36), (10.37), (10.40), (10.41), (10.42) и (10.44) заменить функции tpr (t)

функциями фг (/).

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

10.6. Оптимальные системы описываемые дифференциальными уравнениями

Оптимальные системы фильтрации сигнала известной формы и сигнала неизвестной формы, состоящего из нескольких компонент, а во многих случаях и оптимальная система преобразования сигнала описываются обыкновенными дифференциальными уравнениями. Для вывода дифференциальных уравнений оптимальной системы с бесконечной памятью в общем виде рассмотрим задачу воспроизве­ дения векторного сигнала S (t) с составляющими (/), . . ., 5„ (t), определяемого системой обыкновенных линейных дифференциаль­ ных уравнений

Si ^ i

a iiSi

(10.46)

/=i

 

( i = l ,

. . ., я)

 

со случайными начальными условиями S(- (/„) = U[ (i =

1, . . ., п).

Как известно из теории дифференциальных уравнений, общий инте­ грал системы уравнений (10.46) представляет собой линейную комби­ нацию п линейно независимых частных интегралов этой системы уравнений. Обозначим через ф1/( (Л (i = 1, . . ., я) частный интеграл системы уравнений (10.46), удовлетворяющий начальным условиям

ф(7г (*0) = 0 при i =j= k, Ф/г/г (^о) = 1- Частные интегралы такого вида, соответствующие k = 1, . . ., п, линейно независимы. Поэтому требуемый выходной сигнал со случайными начальными значениями

YTi ih) =

(*0) = Ui (t =

1,

- •

я)

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

формулой

 

YTi(t) =

Si (0 =

t

Ukxpik(t).

(10.47)

 

 

 

 

k = i

 

 

(i =

l,

 

я).

 

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

(0 = t cpiSг (t) +

Np (t) =

t

cpi t Uk%k (t) + Np (t),

1=1

i = 1

k = l

 

( p = 1, . .

.,

m)

.272


где N x (/), . . ., N,n (t) — некоррелированные один с другим белые шумы, интенсивности которых соответственно равны Glt . . ., Gm. Положив для краткости

Фрй(0 = £

cpi%k (/),

(10.48)

 

i=1

 

(Р = 1, • ■- ,

ni]

k = \ ,

n),

будем иметь

 

 

 

X p (t) =

 

 

(10.49)

(p =

1,

•■■,«*)

 

Рассматриваемая задача представляет собой задачу одновремен­ ного линейного алгебраического преобразования т сигналов в п сигналов в случае, когда функции ф,.А(t) связаны с функциями cppjfe (/) формулами (10.48). На основании изложенного в предыдущих пара­ графах оптимальная система для воспроизведения каждой состав­ ляющей требуемого выходного сигнала представляет собой парал­ лельное соединение согласованных фильтров с усилителями выход­ ных сигналов. Обозначим коэффициент усиления к-й компоненты /э-го входного сигнала в выражении i-й компоненты выходного

сигнала через Кр)}. Тогда i-я компонента выходного сигнала

/7 1 П t

Yi (/)= У У

[ фрЛ (т) Х р (т) dx.

(10.50)

 

р =1/1=1

/„

 

 

 

Подставив сюда

выражение

(10.49)

входного сигнала

X p (т)

и вводя для краткости обозначение

 

 

 

 

b$ = J( фр*(т)фрл(т) dx,

(10.51)

 

to

 

 

 

 

(k,

/ i = l ,

/г;

/9 = 1,

.... m)

 

получим

 

 

 

 

 

 

п

т

п

 

 

Yt (о = 2 u k 2 S № № +

k=l

p—1 h=1

m n

t

+ p=lh=lS S

ti1 фрлм n p (t ) dx-

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

(10.52)

273


Вычитая из этой формулы выражение (10.47) для требуемого вы­ ходного сигнала, найдем ошибку на t-м выходе системы:

 

 

 

п

 

/

т

п

 

-'ь./

+

 

 

 

 

е,(/) =2<! = 1

\ p

2= U l =2l

 

 

 

 

 

 

4-

2

2

J срр" w ^

м

dT-

 

 

(10.53)

 

 

/J = 1

/f =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t =

1,

....

п)

 

 

 

 

 

 

 

Отсюда получаем следующее выражение для среднего квадрата

ошибки на

t-м выходе

системы:

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

/ т

п

 

 

 

 

 

 

 

 

 

 

ц п = я ь

s я д а » -

 

 

 

 

 

 

 

г ,

s = l

\ р

А = 1

 

 

 

 

 

 

 

 

 

 

 

( Я )j

 

— ф/sj

+

 

 

 

 

 

 

 

 

т

 

п

 

 

 

 

 

 

 

 

 

 

 

-|-

2

Gp

2

hp/jhp/

 

,

 

 

 

(10.54)

 

 

 

р =

1

Л, (= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

( t= 1, . . . , « )

 

 

 

 

 

 

 

где yrs (г, s = 1,

. . .,

/г)

вторые начальные моменты случайных

величин t/!, . . . .

Д„.

 

 

 

 

 

 

 

смешанные

начальные

Аналогичной

формулой определяются

моменты второго

порядка составляющих

вектора

ошибки системы:

 

Yrs

т

п

 

 

 

\

!

>п

п

 

 

 

 

 

11«у — Я

X

X К Ж

- Ъ

г )

(

Я Ъ

Ф

У

+

 

Г , 5 = 1

р = 1 Л= 1

 

 

 

/ \(7 = 1 / = 1

 

 

/

 

 

 

 

4-

Я Gp

Я

bpfopibW-

 

 

 

(10.55)

 

 

 

р =

1

Л ,/ = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

(о / = 1,

. .

/I)

 

 

 

 

 

 

Дифференцируя выражение (10.54), соответствующее данному

фиксированному

значению

i,

по

параметрам

1}р)\

(р =

1, -

/п;

/г = 1, . . .,

п) и приравнивая производные

 

нулю, получим

 

после

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

уравнений для определения неизвестных

коэффициентов усиления

*1 •

 

 

 

 

Лр!1>

т

и

 

Я Yrs

=

Ьг ■ ■X

Я

G„K

r = l

р = 1 А = 1

...,

п).

(и — 1, ..., т\

s — 1,

Заметим теперь, что левые части этих уравнений не зависят от индекса и. Следовательно, и правые части не зависят от и. Это дает возможность ввести новые неизвестные:

His — Guru's

1, ... , п\ и = 1,

,, т).

(10.56)

27-;


Определив из (10.56) величины Х%] и подставив в предыдущие уравнения, получим

п

Г

т

п

 

£ Угв 1|тг— S

£ G p Хl'- ih b % ] = P / s .

Г — I

L

/3=1/1 =1

 

 

 

(s =

1,

n)

Обозначим через y~ элементы матрицы, обратной по отношению

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

(10.57)

(/"= !, . . . . п)

С помощью полученных уравнений для

формулы (10.54)

и(10.55) приводятся после преобразований к виду

Пп

% = £

= £

М ф? •

(10.58)

г/=1

ч=1

 

(г, / = 1,

....

п)

 

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

выкладок, представим все полученные формулы и уравнения в ма­ тричном виде. Вводя вектор-столбец S и матрицу А:

(10.59)

представляет собой фундаментальную матрицу решений уравнения

(10.59). Поэтому матрица 'F удовлетворяет

следующим уравнению

и начальному условию:

.

 

 

V (/„) = / ,

(10.60)

18*

275