где {aft} — стационарная последовательность, статистически незави симая от последовательности значений {tft}, которые упорядочены и распределены по пуассоновскому закону с параметром К. Мы исклю чили в (8.94) член с k = 0, что подробнее обсуждается ниже. Чтобы вычислить среднее значение и автокорреляцию для процесса х, выби раем некоторое значение t, а затем пронумеруем {tft} в реализации так, что t 1 — первая точка справа от t и t_x — первая точка слева от t, как показано на рис. 8.14.
На первый взгляд это может показаться странным, но длина интервала между t_x и txстатистически отлична от длины всех других
интервалов |
между импульсами последовательности. Причина этого |
в том, что |
интервал t_x, t x строится так, Чтобы включить точку t, |
в то время как для других интервалов такой необходимости нет. Плот-
Рис. 8.14. Типичная реализация случайной последовательности, ил люстрирующая нумерацию импульсов.
ность вероятности для отрезка времени х от любой точки временной оси до точки появления следующего случайного импульса задается выражением рх (а) = ^е- *-0. Поэтому, если положить tQ= t и опре делить случайную переменную = t k — tfe _ то хь имеют одина ковые распределения для всех k (включая нуль). Учитывая сказанное, можно переписать (8.94) более точно:
Х(0= |
2 a ,s ( i- t» ); |
= |
(8.95) |
k = |
— со |
Н 0 |
|
Другой корректный способ рассмотрения заключается в том, чтобы принять произвольную точку t за начало отсчета времени, тогда слу чайные величины tfc — t для k > 0 и t — t h для k < 0 соответствуют плотностям распределения (8.87) и (8.88). С учетом (8.95) можем написать
|
|
ОО |
|
ОО |
00 |
|
|
Е [х (t)] = |
2 Е |
Е ts |
= 2 |
a J s (a) KPh^ (о) da + |
|
k = |
— оо |
|
k = \ |
О |
|
|
— 1 |
_ o o |
|
|
_ o o |
|
OO |
. f \ k — l |
+ 2 |
a j S(-g)^P_ft_ ^ | ) ^ = Xa5S( 0 e - « |
2 |
+ |
k = — оо 0 |
|
|
0 |
|
k = |
1 |
+ |
Я a |
s(t)e.%t ^ |
( - U ) i—l |
|
^ |
s (t) dt = A,a q. |
|
dt = h a |
|
|
/ = |
i |
|
|
|
|
(8.96)
Для вычисления автокорреляционной функции предположим, что
{ah} |
статистически |
независимы, |
так |
что |
|
|
|
|
|
Е К |
|
а2 |
для k = j=j=О, |
|
(8.97) |
|
|
а;] |
а2 |
для k Ф /, |
|
|
Тогда имеем |
|
|
|
|
О для |
k |
или / = 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А**(^ + т, |
0 = ^[х(^ + т) х (0] = |
|
|
= |
2 2 |
|
£ [a Aa , ] £ [ s ( f - t fc+ |
T )s(* -t,)]. |
(8.98) |
|
k ~ |
— СО |
/= |
— 00 |
|
|
|
|
|
|
|
Сумма в (8.98) |
вычисляется |
путем отдельного рассмотрения членов, |
для |
которых |
k = j > |
1, |
k = j < |
|
— 1, |
k > |
j > 1, |
j > k > 1, |
j C k < — 1, k < j < 1, k > l, j < — 1, / > 1 , k < — 1. |
После |
некоторых |
преобразований |
получается |
простой |
результат: |
где |
|
|
|
^жж(х) = |
^а2 г (х) + |
ад)2, |
|
(8.99) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
со |
s(t + x)s(t)dt; |
q= |
оо |
s{t)dt. |
|
|
г(т) = |
^ |
^ |
|
|
|
|
|
— 00 |
|
|
|
|
— оо |
|
|
В заключение отметим эквивалентность этого процесса и того из процессов, описанных в § 8.2, для которого моменты появления им пульсов независимы и равномерно распределены на временной оси
[см. уравнения (8.10) и (8.11)].
Упражнение 8.12. Проделать преобразования, опущенные при выводе (8.99).
СПИСОК ЛИТЕРАТУРЫ
1.Р а г z е п Е. Stochastic process. Holden-Day, 1962.
2.B e n n e t W. R. Statistics of regenerative digital transmission. — «Bell Sys. Tech. Jour.» 1958, v. 37, p. 1501—1542.
3. |
P a p о u 1 i s A. Probability, random variables, and |
stochastic processes. |
4. |
McGraw-Hill, |
1965. |
|
plant. — «Bell Sys. |
А а г о n M. |
R. PCM Transmission in the exchange |
5. |
Tech. Jour», |
1962, 41, p. 99—141. |
techniques. — «Trans. |
L e n d e r |
A. Correlative |
digital communication |
|
IEEE». 1964, |
COM-12, p. |
128—135. |
|
6.К r e t z m e r E. R. Binary data communication by partial response trans mission. Paper № CP65-419 at the 1965 — «IEEE Communications conven tion», Boulder, Colorado.
7.J о h n s M. V. Spectral analysis of a process of randomly delayed pulses. — «Trans. IEEE», 1960, v. IT-6, № 4, p. 440—444.
8.F r a n k s L. E. A Model for the random video process. — «Bell Sys. Tech. Jour.», 1966, v. 45, p. 609—630.
9.H i l d e b r a n d F. B. Methods of applied mathematics. Prentice-Hall, 1952.
10. L e n d e r A. The duobinary technique for high-speed data transmission. —
«IEEE |
Trans |
on |
Communication and |
Electronics» 1963, v. 82, p. 214—218. |
11. L a n d i n g |
J. |
H. a n d |
B a t t i n |
R. H. Random precesses in auto |
matic |
control. |
McGraw-Hill, |
1956. |
|
9
ОПТИМАЛЬНАЯ ФИЛЬТРАЦИЯ СИГНАЛОВ
9.1. ВВЕДЕНИЕ
Разработанные в двух последних главах способы описания слу чайных процессов позволяют более полно рассмотреть вопросы про ектирования систем оптимальной фильтрации сигналов. Это обуслов лено тем, что мы теперь с помощью линейных и квадратичных функ ционалов можем описывать характеристики как детерминированных, так и случайных элементов системы. Применяя вариационные методы гл. 6, мы будем оптимизировать некоторые детерминированные эле менты системы для повышения ее эффективности. В частности, можно оптимизировать форму сигнала, применяемого в системе с АИМ, или импульсную характеристику фильтра, выделяющего сигнал из шума. В этой главе рассматривается проблема фильтрации. Мы имеем возможность исследовать широкий класс задач, требующих отыска ния оптимальных фильтров по критерию минимума среднеквадра тической ошибки. Наглядная интерпретация этих результатов осно вана на введенных ранее понятиях, связанных с пространством сиг налов.
В следующем параграфе дана общая формулировка проблемы оценки параметров сигнала при линейной фильтрации и получено условие на импульсную характеристику фильтра, обеспечивающего минимум среднего квадрата ошибки. Дальнейшие параграфы посвя щены отысканию характеристик, удовлетворяющих этому условию в частных задачах фильтрации. Некоторые из рассмотренных при меров хорошо известны и освещены в литературе по теории систем связи или систем управления. Другие примеры привлекали меньше внимания, но они поддаются тем же методам анализа.
Для упрощения и выявления главных факторов, влияющих на структуру оптимального фильтра, рассматриваемые задачи кое в чем идеализированы. В практических приложениях следует объединять черты нескольких примеров.
Нужно отметить, что во многих случаях оптимальная система об работки сигналов нелинейна. Несмотря на это, исследование линей ной фильтрации имеет большое практическое значение, так как линей ные фильтры используются чаще, чем нелинейные. Мы трактуем ус ловие линейности и условие физической реализуемости (неопережаю щий отклик) как практические ограничения, навязанные конструкто ру. Правда, задачи вначале формулируются без учета физической реализуемости. Помимо простоты, такое приближение оправдано двумя причинами. Во-первых, характеристика оптимального нереа лизуемого фильтра часто может быть хорошо аппроксимирована фи зически реализуемым фильтром с достаточно большой постоянной
времени. Во-вторых, система с нереализуемым оптимальным филь тром служит полезной основой при построении квазиоптимальных систем, удовлетворяющих условию реализуемости. Но в ряде случаев физическая реализуемость является существенным ограничением и важно оптимизацию выполнить с его учетом. Такие задачи обсуж даются в § 9.6.
9.2. МИНИМИЗАЦИЯ СРЕДНЕГО КВАДРАТА ОШИБКИ ПРИ ОЦЕНКЕ ПАРАМЕТРА
Задача оценки параметра в системах передачи сигналов иллюстри руется на рис. 9.1. Передаваемый сигнал х (t) подвергается отобра жению Т г, которое описывает влияния канала передачи. При этом х (t) преобразуется в свой образ г (t) — принимаемый сигнал, кото рый далее проходит через фильтр с импульсной характеристикой
сигнал
Рис. 9.1. К оценке параметра сигнала.
h (t, s). Импульсная характеристика фильтра должна быть выбрана так, чтобы на выходе (в момент t) получить минимум среднеквадрати ческой погрешности оценки параметра ю (i), содержащегося в сиг нале х (t). Оцениваемый параметр сигнала может быть охарактери зован отображением Т0. Иногда полезно считать, что Т0 — описание идеального канала и задача состоит в том, чтобы найти линейный фильтр, включенный последовательно с реальным каналом Ту, кото рый обеспечивает отображение, близкое к Т0.
Предполагая, что сигнал вещественный, запишем функционал среднего квадрата ошибки для момента t в виде
/ = Е [ 1у (* ) - со (/) Н = Е [| у (0 12]~ 2 Е [у (t) о (/)]+ £ [| ю (/) |2 ]. |
(9.1) |
Используя связь входа и выхода фильтра |
|
оо |
|
У (0 = § h{t, s) г (s)ds, |
(9.2) |
представим далее функционал |
/ |
как сумму квадратичного, линейного |
и постоянного функционалов |
относительно h (t, |
s), причем h (t, |
s) |
рассматривается как |
функция |
s при фиксированном |
параметре |
t: |
|
оо |
|
|
|
|
|
|
/ = |
^ kzz (s,o) h (t, |
a) h (t, s) dads— |
|
|
|
—oo |
|
|
|
|
|
|
|
oo |
|
|
|
|
|
|
—2 |
h (/, s)kza,(s, |
t)ds \-.k<o<o(t, |
t). |
(9.3) |
oo —
Как видно, эти функционалы полностью определяются через автокор реляцию и перекрестную корреляцию процессов z (t) и со (t).
Квадратичный функционал соответствует самосопряженному опе ратору, поэтому стационарные точки функционала I определяются при варьировании h (t, s) решениями уравнения
оо |
|
|
|
5 |
Kz (s>о) h (t, a) do = kza(s, t), |
(9.4) |
— oo |
|
|
|
которое получается, |
если приравнять |
нулю градиент |
функционала |
/. В большинстве важных случаев h (t, |
s) принадлежит (как функция |
s) пространству L2 ( — оо, оо ), но иногда нужно расширить функцио нальное пространство, чтобы включить в него решения уравнения
(9.4).
Принцип ортогональности
Другая трактовка проблемы оптимальной фильтрации получает ся, если функционал ошибки рассматривать на пространстве случай ных величин. Пусть Х0 — линейный оператор фильтра, удовлетворяю щего условию (9.4). Тогда (9.4) можно переписать в эквивалентной форме
|
E[{X0z{t) —со (f)} z (s)] = 0. |
(9.5) |
Это значит, |
что Х 0 нужно выбирать так, чтобы случайная |
ошибка |
у (0 — о (г1) |
была ортогональна ко всем случайным величинам {z (s); |
— оо < s < |
оо } из процесса z, соответствующего принимаемому |
сигналу. Такой результат интуитивно ясен: если бы имела место кор реляция между ошибкой и принимаемым сигналом, то при последую щей обработке можно было бы получить лучшую оценку. Используя (9.5), покажем другим способом, что (9.4) есть условие оптимальности. Пусть X — некоторый другой линейный оператор. Тогда
/ = Е [| X z (0 —(о(0 |2] = Е [| {X z (i) — X0z (0) + |
|
~\-{Х0z (*)—оа(0} |2 ] = Е [ | Х0z (0 —и (0 |2 ] + |
|
+2Е \{Х- Х0) z (t) {Хй z(t) —(о (0}J + E[\{X— X0)z(t) |2]. |
(9.6) |
В силу условия (9.5) второе слагаемое в последнем выражении равно нулю. Третье слагаемое неотрицательно, поэтому X не может