Компенсация помех в оптико-электронных устройствах с использованием фильтра Винера-Хопфа

УДК 519.216.3

В. А. Головков, В. А. Смирнов

КОМПЕНСАЦИЯ ПОМЕХ В ОПТИКО-ЭЛЕКТРОННЫХ УСТРОЙСТВАХ С ИСПОЛЬЗОВАНИЕМ ФИЛЬТРА ВИНЕРА—ХОПФА

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

Ключевые слова: сигнал, компенсация помехи, выборка, эффективность, фотоприемное устройство, экспериментальная установка, корреляционная функция.

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

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

Рассмотрим влияние типа выборки и ее размерности на эффективность интерполяции с использованием линейного фильтра Винера—Хопфа. Целесообразно выбирать узлы интерполяции на возможно большем временном интервале, чтобы исключить возможность совпадения узла и самого сигнала. На выходе ФПУ можно получить выборку из отдельных значений помехи и ее производных. В литературе, например [1], указывается, что если считать спектральную плотность помехи ограниченной и применить теорему Котельникова, то использование в узлах интерполяции не только значений случайного процесса, но и его первой производной позволяет в два раза увеличить интервал между узлами интерполяции. При этом необходимо знать все предыдущие отсчеты реализации случайного процесса и граничную частоту спектральной плотности случайного процесса. Различные интерполяционные представления реализаций случайных процессов изучались, например, в работе [2], где приведен список литературы, посвященный этому вопросу. В то же время целесообразно рассмотреть возможность использования простых алгоритмов, позволяющих производить интерполяцию

случайных процессов с минимальной дисперсией, используя линейную фильтрацию, базирующуюся на уравнении Винера—Хопфа [3].

Оценим эффективность алгоритма интерполяции случайного процесса с использованием в узлах интерполяции его значений и нескольких производных (ограничимся случаем двух узлов). Пусть реализация случайного процесса §(V) задана на интервале времени tе¿0 +Т, соответственно узлы интерполяции задаются в моменты времени ^ и ^ +Т . Внутри интервала Т возможно наличие сигнала. В этом временном интервале будем интерполировать значение реализации помехи § ( + т), где текущее значение т задается в интервале 0 < т < Т. Эффективность интерполяции оценим по выборке §(%), §(% +Т) либо по выборке §(^), §'00), §(% + Т), §'(tо + Т), а также по выборке §(^), §§"(%), §(*0 +Т), §'Оо +П § +Т). Винеровский вектор весовых коэффициентов [3] выборки определяется как W = Я^Р, где Я — корреляционная матрица входных сигналов, или корреляционная матрица выборки, Р — вектор-столбец КФ полезного отклика §(^ + т) и элементов выборки. Минимальное значение условной дисперсии ошибки определяется в виде о2 [§(^ + т)/§(^)...§(^ +Т)...] = о2-P*W,

где о — дисперсия случайного процесса, „*" — знак транспонирования. Пусть КФ случайного процесса Я(т) = о р(т), где р(т) — нормированная КФ. Если использована выборка

о 2р

§(%), §(^ +Т) то очевидно, что Р =

Я 2 1 Р(Т) К = о

Р(Т) 1

Если §(%), §'(tо), §(Ч +Т), §'(tо +Т), то используя [4], можно получить

, а матрица Я определяется в виде

Р*

о 2р

а матрицу Я — в виде:

1 0 р(Т) р '(Т) о -р 40) -р '(Т) -р \Т) р(Т) -р '(Т) 1 о

р'

Если использовать выборку £(%), £' (%), £''(t0), ^(t0 + Т), £'(t0 + Т), £''(t0 + Т), то вектор

Я = о2

Р*

/ / / / / / а р(т), -а р' (т),а р "(т),а р(Т-т),а р' (Т-т),а р "(Т-т) , а матрица Я будет иметь сле-

дующий вид:

Я= а2

1 0 р''

0 -р''

р"

р

р' (Т ) -р' (Т ) р(3)(Т ) 0 -р''(0) 0

р "

Для оценки точности интерполяции будем использовать величину нормированной дисперсии, обратную коэффициенту подавления помехи [5]:

К ( т)=-

[ +Т)/£(¿0), £ Щ..^ + Т), £ '(¿0 +Т)...] а2

Численные расчеты значения К (т) проведем для случайных процессов с нормированной КФ р(т) = ехр(-пт ), при таком представлении КФ энергетическая ширина спектральной плотности случайного процесса А/ = 1, этот результат можно получить из [4]. Такой процесс является дифференцируемым, причем спектральная плотность по оси частот не ограничена, хотя этот процесс относится к классу линейно-сингулярных [6].

На рис. 1 приведена К(т) для случайного процесса с КФ р(т) = ехр(-пг2); кривая 1 соответствует выборке £(¿0 +Т), 2 - £(%), £ '(¿0), £(¿0 +Т), £ % +Т), 3 - £(¿0), £ '(¿0), £ '(¿0), £^0 + Т),£' (¿0 + Т),£" (¿0 + Т) . Характер расчетных кривых на рис. 1 очевиден и не требует пояснений.

к

Рис. 1

Можно показать, что алгоритмы с использованием выборки из значений случайного процесса и их производных более эффективны, чем алгоритмы, в которых используется выборка из значений случайного процесса, расположенных вблизи узлов интерполяции типа £(Г0 -2Аt), £(^ -Аt), £(Г0), £(Т), £(Т+Аt), +2Аt), где Аt — некоторый интервал времени. Эффективность последней выборки равна эффективности выборки £(¿0), £' (¿0), £"(¿0), £(¿0 + Т), £ '(¿0 + Т), £ ''(¿0 + Т) при условии, что Аt ^ 0 . Подобные исследования были подробно описаны в работе [7] для алгоритмов прогнозирования случайных процессов, но и при интерполяции случайных процессов с распространенными КФ это положение верно.

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

рис. 2. Установка состоит из модуля имитато-

ух

X"

17

Т

фэу

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

излучением лампы накаливания 5, имитирующим фоновое излучение. Далее излучение рассеивалось пластиной 6 и после ослабления с помощью фильтров-ослабителей 7 регистрировалось ФПУ 2.

Благодаря применению фильтров 3 и лампы накаливания 5 формировалось любое желаемое отношение сигнал/помеха. Применение фильтров 7 обеспечивало возможность изме-

нения помехи от уровня, определяемого дробовым шумом фотокатода и квантовыми флук-туациями фонового потока, до уровня „гладкой" фоновой подложки.

Для получения экспериментальных зависимостей эффективности компенсации помехи наряду с импульсным лазером применялся He—Ne-лазер ЛГН 126 непрерывного излучения с длиной волны ^ = 0,63 мкм, прерываемого модулятором со сменными секторными дисками и редуктором скорости вращения.

Фотоприемное устройство 2 выполнено на основе фотоэлектронного умножителя ФЭУ-175 (с мультищелочным катодом). Для максимально возможного подавления сетевых помех рабочее напряжение ФЭУ формировалось от источника постоянного напряжения (аккумулятора) с применением повышающих преобразователей, от него запитывалась лампа 5.

При постоянной времени ФПУ 100 мкс длительность импульса лазерного излучения была установлена равной 500 мкс. Сигнал с выхода ФПУ поступал через плату аналого-цифрового преобразователя (АЦП) в персональный компьютер, где и проводилась его обработка с целью компенсации помех и выделения принимаемого сигнала.

Частота дискретизации АЦП составляла 5 мкс. Помеха на выходе ФПУ являлась стационарной и реализации помехи в течение одной секунды при отсутствии сигнала были использованы для оценки ее корреляционной функции. Максимальное время корреляции было выбрано равным 3000 мкс. Дальнейшее увеличение времени наблюдения не приводило к изменению оценки корреляционной функции R(t), что позволило сделать вывод о высокой точности ее оценки.

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

R(t) = aexp -ат2 -b exp -Р(т-91)2 -b exp -Р(т+91)2 -c exp -в(т-02)2 -

-c exp [-р

Значения параметров a, а выбирались исходя из вида центрального лепестка оцененной КФ, b, в, с, 01, 02 — исходя из задачи интерполяции боковых лепестков КФ. Выбор такой интерполирующей функции в виде суммы гауссоид позволяет довольно просто интерполировать многомодальную КФ, так как гауссоиды одномодальны и достаточно быстро спадают (в зависимости от множителя в показателе экспоненты).

Оцененная КФ (кривая 1) и интерполирующая ее функция (2) приведены на рис. 3.

r(t), в2 ^

-2000

0

Рис. 3

Интерполяция помехи проводилась в соответствии с выражением

ад, +т) = ВДЩ,)+¿2(т)£ ' (¿0)+кз(т)£ Ч'0>+*4(тЖ* 0 +T)+ +¿5(т)£ '(¿0 + T)+¿б(т)£' (¿0 + T),

где коэффициенты к (т)при г = 1,..., 6 рассчитывались как элементы вектора W, причем в матрицу Я, выражение (3) и вектор Р значения КФ подставлялись, согласно выражению (5), для реально оцененной КФ.

На рис. 4 приведена осциллограмма импульса с аддитивной помехой на выходе ФПУ (кривая 1) и та же реализация, но при помехе, скомпенсированной (2) путем вычитания ее интерполированных значений. Оценка первой и второй производной в моменты = 0 и Т = 1000 мкс проводилась методом конечных разностей. Как видно из рис. 4, помеха была довольно хорошо компенсирована — до появления импульса сигнала она практически сведена к нулю, то же видно и

и, в 0,6

0,4 0,2

после окончания импульса.

х л

'У/

0 200 400 600 800 t, мкс

Рис. 4

Таким образом, интерполяция помехи линейным нерекурсивным фильтром Винера— Хопфа за счет оценки значений реализации и ее производных в начальной и конечной точках, а также вычитания ее из наблюдавшейся реализации позволяет компенсировать помеху, восстановив реализацию сигнала на выходе ФПУ. Применение таких методов компенсации помехи целесообразно при возможности обработки уже записанного сигнала.

СПИСОК ЛИТЕРАТУРЫ

1. Хургин Я. И., Яковлев В. П. Финитные функции в физике и технике. М.: Наука, 1971. 408 с.

2. СосулинЮ. Г. Теория обнаружения и оценивания стохастических сигналов. М.: Советское радио, 1978. 320 с.

3. Уидроу Б., Стирнз С. Адаптивная обработка сигналов. М.: Радио и связь, 1989. 440 с.

4. Тихонов В. И., Хименко В. И. Выбросы траекторий случайных процессов. М.: Наука, 1982. 303 с.

5. Ширман Я. Д., Манжос В. Н. Теория и техника обработки радиолокационной информации на фоне помех. М.: Радио и связь, 1981. 416 с.

6. Розанов Ю. А. Стационарные случайные процессы. М.: Наука, 1990. 272 с.

7. Головков В. А. Прогнозирование случайного процесса по выборке его производных // Радиотехника и электроника. 1993. № 3. С. 1049—1053.

Сведения об авторах

Владимир Алексеевич Головков — канд. техн. наук; НИИ комплексных испытаний оптико-электронных

приборов, Сосновый Бор; E-mail: golovkov_ggg@mail.ru Виктор Александрович Смирнов — НИИ комплексных испытаний оптико-электронных приборов, Сосновый Бор; старший научный сотрудник

Рекомендована институтом Поступила в редакцию

25.09.08 г.