Оптимальная частотно-временная фильтрация

Оптимальная частотно-временная фильтрация

Оптимальная частотно-временная фильтрация

канд. биологических наук М.П.Иванов, д-р техн. наук В.В.Кашинов

Санкт-Петербургский государственный университет

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

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

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

 (1)

и (2)

где - опорное напряжение; - импульсная переходная функция стационарного фильтра. Таким образом, ядра операторов (1,2) представлены в виде произведения k(t, t )=r (t )h(t-t ), где r - заданная известная функция, а h(t-t ) подлежит оптимизации.

При оптимизации будем использовать метод, разработанный в заметках [2-4]. В качестве критерия оптимальности выберем отношение сигнал-шум и представим значение полезного сигнала на выходе в момент максимума t0 с помощью фильтрующего свойства d - функции в виде линейного функционала

 (3)

Момент t0 заранее неизвестен и принадлежит интервалу наблюдения T, а шум на выходе частотно-временного фильтра является нестационарным. Поэтому в качестве критерия оптимальности примем функционал I отношения пиковой мощности сигнала к средней по времени и по ансамблю мощности шума Pш на выходе фильтра

 (4)

Здесь и в дальнейшем угловые скобки < T > обозначают усреднение по ансамблю.

Поскольку сомножитель h(t-t ) в ядре операторов (1) и (2) - разностный, можно использовать свойство свертки и записать

 (5)

и (6)

Таким образом, критерий оптимальности имеет вид

 (7)

где Sвых(t0) и Pш выражаются формулами (5) и (6). Здесь уже можно применить методы обобщенного вариационного исчисления [2-4]. Обобщенное уравнение Эйлера-Пуассона для функционала (7) имеет вид

 (8)

где M - коэффициент пропорциональности, не влияющий на вид коэффициента передачи фильтра.

Переходя к спектрам и обозначая соответствие функций и их преобразований Фурье;;;, получаем выражение для коэффициента передачи стационарной инерционной части оптимального частотно-временного фильтра

 (9)

где * обозначает комплексное сопряжение.

При любом выборе опорного напряжения r (t), которому соответствует спектр P(W ), не уменьшающем энергию сигнала, и любой помехи, в том числе и узкополосной, выводящей GPS из строя [1], существует коэффициент передачи K(w ) оптимального стационарного фильтра h(t ). Мощность множества пар r (t) и K(w ) может быть больше мощности континуума [3]. Даже для рассмотренного простейшего случая все ограничения для r (t) и K(w ) не определены. Из существования решений для частного случая задачи [3] следует существование множества ядер k(x,t), доставляющих функционалу (7) экстремум, причем значения этого экстремума для каждого k(x,t) из этого множества - одинаковые. Решением оптимизационной задачи будет конструктивное описание этого множества оптимальных ядер. Если r (t)=const, т.е. перемножитель отсутствует, P(w -W )=d (w -W ) и получается согласованный фильтр; если r (t)=S(t), получается коррелятор.

Таким образом, и корреляционный прием, и согласованная фильтрация являются частными предельными случаями частотно-временной фильтрации. Опорное напряжение r (t) и переходную функцию фильтра h(t-t ) следует выбирать, исходя из удобства реализации. А для осуществления оптимального приема при белом шуме применение коррелятора или согласованного фильтра обязательным не является.

Решение сформулированной задачи заведомо неоднозначное. Для описания этого множества потребуется использовать понятие функции спектральной корреляции.

Представляя знаменатель в выражении (9) в виде двойного интеграла и меняя порядок интегрирования и статистического усреднения, получаем

 (10)

Обозначим B(w 1,w 2)=; это выражение называется функция спектральной корреляции (ФСК) [5]. Если не учитывать свойства ФСК при мультипликативном воздействии на входной процесс, можно получить ошибочные результаты типа превышения потенциальной помехоустойчивости [6].

ФСК выражается через автокорреляционную функцию B(t1 ,t2 )

 (11)

Средняя мгновенная мощность B(t1,t2) нестационарного процесса может быть выражена через ФСК

 (12)

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

Средняя по времени спектральная плотность мощности нестационарного процесса может быть выражена через ФСК

 (13)

Спектральная плотность нестационарного процесса характеризует вклад составляющих в интервале частот (w + dw ) и всех коррелированных составляющих с другими частотами.

Для стационарных процессов автокорреляционная функция зависит только от разности моментов времени t = t1 vt2 , и в этом случае

 (14)

Для стационарных процессов все частотные составляющие некоррелированы.

При модуляции стационарного белого шума детерминированным опорным напряжением r (t) ФСК зависит только от разности частот

 (15)

где D w = w 1 - w 2 . Например, при стробировании стационарного белого шума периодической последовательностью импульсов средняя по времени спектральная плотность уменьшается в скважность раз; это можно наблюдать на экране анализатора спектра. Но появляется свойство, которое нельзя наблюдать на экране анализатора спектра - между спектральными составляющими появляется корреляция.

Парадокс. Предположим, что осуществляется оптимальный прием отрезка периодической последовательности импульсов на фоне белого шума. Как известно, оптимальным в данном случае является согласованный гребенчатый фильтр. Теперь включим на входе оптимального гребенчатого фильтра стробирующее устройство (перемножитель на последовательность прямоугольных импульсов единичной амплитуды) так, чтобы импульсы сигнала проходили без искажений. Спектральная плотность шума на выходе стробирующего устройства уменьшится в скважность стробов раз. Казалось бы, что отношение сигнал-шум на выходе гребенчатого фильтра должно увеличиться, но оно и так было максимально возможным, поскольку фильтр оптимальный. Разрешить парадокс помогает появление корреляции между спектральными составляющими. Ясно, что суммирование "гребенок" фильтра со сфазированными гармоническими составляющими сигнала и коррелированными составляющими шума результирующее отношение сигнал-шум не повысит.

Частотно-временная фильтрация может с успехом использоваться в спектральных дискриминаторах временных интервалов [7]. В некоторых радиоканалах, например, телеметрических каналах сверхдальней космической связи или GPS [1], отношение сигнал-шум оказывается Pс /Pш