Анализ точности оценивания биспектра в задачах обнаружения сигналов на фоне помех Текст научной статьи по специальности «Математика»

РАДИОТЕХНИКА

УДК621.396:681.34

АНАЛИЗ ТОЧНОСТИ ОЦЕНИВАНИЯ БИСПЕКТРА В ЗАДАЧАХ ОБНАРУЖЕНИЯ СИГНАЛОВ НА ФОНЕ ПОМЕХ

ТОЦКИЙ А.В., НАУМЕНКО В.В._________

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

1. Введение

Бурное развитие и совершенствование современных методов и алгоритмов цифровой обработки сигналов наряду с прогрессом в области средств вычислительной техники позволяет специалистам ставить и решать совершенно новые задачи по обнаружению сигналов на фоне интенсивных помех, различению сигналов, распознаванию и классификации образов на основе анализа их тонкой структуры. Классический спектральный анализ, традиционно используемый при обработке процессов, основан на оценивании спектров второго порядка с помощью традиционного преобразования Фурье [1]. Для процессов, у которых спектральные компоненты статистически независимы, оценка энергетического спектра является исчерпывающей характеристикой при традиционном некогерентном спектрально-корреляционном анализе свойств таких процессов. Однако данный подход далеко не всегда приемлем в ситуациях, когда требуется решить задачи выявления частотно-фазовых зависимостей, содержащихся в наблюдении, или различения двух или нескольких процессов с использованием информативных признаков в виде различий их плотности распределения вероятности. В последнее время для решения этих задач во многих технических приложениях, таких как радиолокация [2-4], астрономия [5], распознавание образов [6], цифровая связь [7] и другие, стали широко использоваться перспективные методы и алгоритмы обработки данных с использованием кумулянтных функций третьего порядка и их преобразований Фурье - спектров третьего порядка (биспектров). Данные методы позволяют увеличить

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

Предположим, что на входе цифровой измерительной радиотехнической системы регистрируют последовательность М реализаций аддитивной смеси полезного сигнала s(i) и шума радиоканала связи nG(i) с неизвестными параметрами.

В рамках рассматриваемой модели наблюдения m-ю (m=1,2,.. ,,M) реализацию x(m)(i) представим в виде

x(m) (i) = s(i -T(m)) + nGm)(i), i = 0,2,..., I -1, (1)

где nGm)(i) - m-я реализация аддитивного стационарного гауссова шума с нулевым средним значением <nG(i)>=0 и выборочной дисперсией, величина которой ст^2 полагается априорно неизвестной; x(m) -принимающая целочисленные значения случайная задержка сигнала, максимальная величина которой не превышает I.

На основе уравнения наблюдения (1) следует решить задачу обнаружения сигнала s(i) на фоне интенсивных помех. В качестве тестовой статистики обнаружения сигнала будем использовать оценку биспектра наблюдения (1), устойчивую к воздействию гауссова шума и равную

В (p, q) = (X(m) (p)X(m) (q)X*(m) (p + q))M, (2)

где биспектр (биспектральная плотность)

Bx (p.q) = Вx (p,q)

exP[jYx(p,q)]

это комплексная

функция двух частотных переменных p=-I+1,.. ,,I-1,

РИ, 2011, № 1

3

q= -I+1,...,I-1;

Вx (p,q)

и Yx (p, q) - соответственно

амплитудный (биамплитуда) и фазовый (бифаза) биспектр рассматриваемого процесса x(m)(i);

I-1

X(m)(p) = ^ x(m)(i)exp[-j2nip/I] - преобразование

i=0

Фурье процесса (1) на конечном интервале времени его наблюдения; <.>M - процедура осреднения по ансамблю из М наблюдений.

Вероятность обнаружения сигнала по оценке биспектра (2) зависит от точности измерения последней функции. При этом интерес представляет задача определения нижней границы точности оценки биспектра -предельная минимальная погрешность оценки.

Цель данной работы заключается в исследовании точности оценивания биспектра и определения потенциального предела дисперсии оценки биспектра сигнала на фоне аддитивных гауссовых шумов.

2. Анализ точности оценивания биспектра

По правилам математической статистики оценку биспектра полагают состоятельной по вероятности, если при многократных измерениях, т. е. при проведении К экспериментов (испытаний) и осреднении выборочных оценок по К испытаниям (K ^ да) оценка биспектра стремится к биспектру.

Оценку биспектральной плотности (2), полученную в произвольном k-м (k=1,2,... ,K) испытании, представим в виде

Bxk (p. q) = (Bxmk (p. q)) = Вs(p, q) + вnk (p, q) =

\ /M

=Re [B s (p, q)] + Re[B nk (p, q)] + +j{Im[B s (p, q)] + Im[JB nk (p, q)]},

Так как проводимые далее математические выкладки подобны для Re и Im частей биспектральной оценки, то достаточно ограничиться подробным выводом формул для одной из частей, например, для вещественной компоненты, а затем получить результаты для мнимой компоненты.

Условная плотность вероятности вещественной части оценки биспектра для k-го испытания равна

P[BRe (p, q) | BRe (p, q)] = 1/ V2naRek (p, q) x

x expbrr^—r{Re[B xk (p, q)] -

2c^Rek(p,q) (5)

- < Re[B xk (p, q)] >k}2},

где

M

aRek (p, q) = 1/(M - 1)£{Re[JB nmk (p, q)] -

m = 1 M

-Re[1/M)£ В nmk (p, q)]}2

m = 1

- выборочная, т. е. полученная на ограниченной выборке объемом М реализаций дисперсия помехи в оценке биспектра наблюдения (1).

Отметим, что оценка Re < Вxk (p, q) >K^„ ^ BsRe (p, q), т. е. является несмещенной, за исключением частотных отсчетов, лежащих на осях p = 0, q = 0 и p = q.

Более удобна для анализа логарифмическая функция правдоподобия, равная

ln{P[BRe(p,q) | BRe(p,q)]} = ln{1/V2roaRek(p,q) x

x expbrr^—r{Re[]Bxk (p> q)] -

2caRek(p,q) (6)

- < Re[B xk (p, q)] >k}2}}.

где ]B s (p, q) - биспектр полезного сигнала (сигнальная составляющая); Bnk(p,q) - составляющая помехи в оценке биспектра.

Будем полагать, что в соответствии с центральной предельной теоремой Ляпунова при достаточно большом количестве независимых экспериментов К закон ^ распределения оценок Re[Bxk(p,q)] и

Im[Bxk (p, q)]} стремится к нормальному. Тогда функцию правдоподобия - многомерную (размерности К) условную плотность вероятности вектора наблюдаемой оценки Bx(p,q) при условии неизвестного биспектра В s (p, q) запишем в виде

P[Bx(p,q)|B s (p, q)] =

= P[BRe (p, q) | BRe (p, q)]P[^R” (p, q) I BIm (p, q)], (4)

где P[BRe (p, q) | BRe (p, q)] и P[Bf (p, q) | BSm (p, q)] - условные плотности вероятности векторов вещественной (Re) и мнимой (Im) частей формируемой оценки при условии неизвестных полей - сигнальных компонент В Re(p,q) и В Sm(p,q) соответственно.

В качестве условной оценки IBRR фп (p, q) вещественной части биспектра выберем такое значение переменной В Re (p, q), для которого логарифмическая функция правдоподобия, полученная в результате статистического усреднения по К наблюдениям функции (6), достигает максимума. Условную оценку В^МФП (p, q), которая удовлетворяет данному критерию качества - максимуму средней функции правдоподобия, будем считать максимально правдоподобной и обозначать МФП.

Условие максимума средней логарифмической функции правдоподобия для вещественной части оценки биспектра запишем в виде:

/ d{ln{P[BRe (p, q) | В5Кб (p, q)]}

\ дBRe (p, q)

Если существует такая несмещенная эффективная

оценка В8^фф (p, q), которая обеспечивает наименьшую среднеквадратическую ошибку по отношению к

= 0.

(7)

K

4

РИ, 2011, № 1

истинному биспектру ВRe(p,q) среди всех возможных оценок В^Мфп (p,q), т.е. при выполнении условия

(Йфф (p, q) - ]BsRe (p, q)]2) < ([В Rm фп (p, q) - B Re (p, q)]2),(B)

то уравнение правдоподобия (7) имеет единственное решение, равное ВR^ (p, q).

Уравнение правдоподобия (7) обеспечивает необходимое, но не достаточное условие максимума.

С учетом (6) уравнение (7) запишем как

3. Интерпретация результатов компьютерного моделирования

На основе теоретических результатов, полученных в разделе 2, исследуем точность оценивания биспектральной плотности в следующих условиях. Во-первых, рассмотрим влияние на точность оценки биспектра вклада отсчетов, лежащих на частотных (p=0; q=0) и диагональной осях (p=q). Во-вторых, изучим поведение точности оценки биспектральной плотности для сигналов с нулевым средним значением и сигналов с отличным от нуля средним значением. В-третьих, оценим влияние оконного сглаживания на точность оценивания биспектральной плотности.

1 h В Re (p, q) - В Re (p, q)- < Bfe (p, q) > K = 0 k£ <^Rek(p,q) . (9)

Для количественного анализа точности оценки биспектра ВRe (p, q) воспользуемся вычислением второй производной логарифмической функции правдоподобия в максимуме, взятой с обратным знаком.

Тогда погрешность оценки биспектра определим с помощью второй производной по ВRe(p,q) уравнения (9), взятой с обратным знаком. В результате получим уравнение

-1 h д2{lnP[BRe(p,q)|]BRe(p,q)]} =

K h д2В Re (p, q)

= 1 h 1 (10)

Kfcl <^^ek (p,q).

Согласно критерию Рао-Крамера дисперсия оценки вещественной части биспектра, полученная по максимуму функции правдоподобия (10), определяется условием:

var[B .'мфп (p,q)] ^ -т~к—і—’ (11)

- h------1---

Kk=1 ^Rek(p,q)

а дисперсия оценки мнимой части биспектра равна

var[]BГмфп (p. q)] ^ 1 к 1 і—. (12)

Kk=1 ^;Lk(p,q)

С учетом условий (11) и (12) предельную среднеквадратическую точность оценки комплексного биспектра определим следующим условием:

var[Bs мфп (p>q)] ^

1 к

1 h

1

1 к

1 h-

1

Kk=1 CTRek(p,q) Kk=1 CTimk(p,q)

-.(13)

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

ст^”)2 , и для различных случайных задержек т(т) исходного сигнала.

Для улучшения оценки биспектральной плотности в работах [8, 9] предложено использовать вычитание среднего значения из наблюдаемого на входе сигнала. Данный подход к улучшению биспектральной оценки теоретически позволяет исключить искажения, вызванные присутствием сигнально-зависимых

членов в оценке 13nk(p,q) (см. формулу (3)), однако это справедливо при устремлении объема выборки к бесконечности. Поэтому представляет интерес иссле-дов ание точности оценки биспектральной плотности в практических условиях ограниченного объема выборки и в условиях растекания биспектра.

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

Дисперсию ст2 (p, q) оценки биспектра определим с помощью суммы дисперсии оценки вещественной CTRe (p, q) и мнимой ст2т (p, q) частей и запишем как

ст2 (p> q) = <CTRe (p. q)+CTL (p> q) =

1 m 1 M

=M 7h{Re[Bnmk(p,q)] -Re[rMhВnmk(p,q)]}2 +

і m' . 1 M. (14)

+ M"“T ^М[Вnmk (P, q)] - Im[TT h Вnmk (P, q)]}2 .

В качестве меры точности оценивания биспектра введем понятие полного стандартного (среднеквадратичного) отклонения (TOSD), величину которого определим в виде

CTTODS ~ /hhCT2 (p,q) (15)

V p q ’

и ограничимся расчетами в пределах главной треугольной области симметрии биспектра 0 < p < q < I/2, p+q < I/2.

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

РИ, 2011, № 1

5

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

Тестовые сигналы генерировались на дискретной сетке 1=64 отсчетов со случайным сдвигом т(т) , имеющим равномерный закон распределения, и с максимальной девиацией тттХ =10 отсчетов.

Проанализируем вначале поведение полного стандартного отклонения (SD = as (p,q) ) по графикам, приведенным на рис. 1,2 и построенным для тестового сигнала - импульса прямоугольной формы длительностью Т = 2 отсчета для фиксированного объема выборки М = 100 реализаций и различных величин среднеквадратических отклонений (СКО) шума an . Для сравнительной оценки влияния вклада величины среднего значения сигнала на рис. 2 представлены результаты, полученные в случае, когда среднее значение сигнала вычиталось перед формированием бис-пекттальной оценки.

Рис. 1. Функция SD = ст^Оз.сО ; =0.5

Рис. 2. Функция SD = аЕ(p,q) , an = 0,5; среднее значение сигнала равно нулю

Анализ результатов, представленных на рис. 1, 2, демонстрирует следующее:

- поведение функции oz (p,q) в двумерном частотном пространстве (p, q) показывает, что значительный

вклад в ее величину вносят сигнально-зависимые составляющие в оценке функции (3). Составляющая помехи в оценке биспектральной плотности сосредоточена вблизи нулевых и низких частот, так как характер изменения функции аЕ (p, q) соответствует поведению амплитудного спектра Фурье тестового сигнала |S(p)|=T|sinc(2 пpT/I)|;

- вычитание среднего значения сигнала приводит к уменьшению максимального значения функции az (p, q) . Однако, хотя вычитание среднего значения и улучшает оценку биспектральной плотности, последняя функция остается искаженной помехой при огра -ниченных на практике длительности интервала наблюдения I и объема выборки М наблюдаемых реализаций;

- случайная задержка т(т) исходного сигнала не влияет на точность оценки биспектральной плотности.

Результаты расчетов стандартных отклонений: полного (TOSD) и усеченного (TRSD), полученные для объема выборки M=100 реализаций, представлены в табл. 1 - 3.

Таблица 1

Результаты, полученные для импульса прямоугольной

формы

TOSD T °n

0,1 0,3 0,5

2 896 5,87103 1,85104

6 1,84103 1,02104 2,60 104

10 2,63103 1,26104 3,33-104

TRSD 2 760 5,04103 1,60 104

6 1,34103 7,89 103 2,13104

10 1,66103 8,62 103 2,54104

Таблица 2

Результаты, полученные для импульса треугольной

формы

TOSD T °n

0,1 0,3 0,5

4 694 5,39103 1,85104

12 1,06103 7,33-103 2,41104

20 1,30103 9,26 103 2,60 104

TRSD 4 578 4,6-103 1,60104

12 685 5,35103 1,93104

20 680 5,86 103 1,92 104

Таблица 3

Результаты, полученные для импульса гауссовой формы s(i) = exp[-a(i -1/2-1)]2

TOSD a °n

0,1 0,3 0,5

0,78 665 5,3 103 1,84 104

0,09 1,12 ■ 103 7,49 103 2,43 104

0,03 1,43103 9,44 103 2,89 104

TOSD 0,78 552 4,52 103 1,58 104

0,09 728 5,70 -103 1,95104

0,03 737 6,03103 2,08 104

6

РИ, 2011, № 1

Таблица 4

Результаты, полученные для импульса прямоугольной

формы

TOSD T CTn

0,1 0,3 0,5

2 810 5,47-103 1,98104

6 1,57103 8,61 ■ 103 2,54-104

10 1,97103 9,83103 2,84-104

TRSD 2 738 4,81 ■ 103 1/74T04

6 1,34103 7,54-103 2,20-104

10 1,66103 8,35103 2,50-104

Таблица 5

Результаты, полученные для импульса треугольной ___________________________________________формы

TOSD T CTn

0,1 0,3 0,5

4 623 4/76T03 1,83104

12 854 6,53103 2,17-104

20 861 7,04-103 2,34-104

TRSD 4 566 4,25-103 1,60104

12 710 5,60103 1,90104

20 663 5,82103 1,92-104

Для сравнения в табл. 4 - 6 представлены результаты, полученные при вычитании среднего значения тестового сигнала. Длительность импульсов треугольной и гауссовой формы была выбрана для корректности сравнения в соответствии с длительностью импульса прямоугольной формы, чтобы обеспечить равенство энергии трех тестовых сигналов.

Таблица 6

Результаты, полученные для импульса гауссовой

формы

TOSD a CTn

0,1 0,3 0,5

0,78 610 5,11103 1,80104

0,09 887 6,47-103 2,24-104

0,03 938 7,25-103 2,47T04

TOSD 0,78 552 4,52-103 1,58104

0,09 728 5/70T03 2,08T04

0,03 737 6,03103 2,08T04

Сравнительный анализ результатов, приведенных в табл. 1 - 6, дает возможность сделать следующие выводы.

1. У странение отсчетов биспектра, лежащих на границах главной треугольной области, позволяет уменьшить величину TRSD по сравнению с величиной TOSD. Этот результат компьютерного моделирования соответствует теоретическому утверждению, высказанному в работе [9].

2. Для количественного анализа удобно ввести величину RTT, равную отношению TOSD к TRSD. Результаты анализа показывают, что величина RTT растет с увеличением длительности тестового сигнала и значительно превышает единицу при малых величинах СКО шума. Выигрыш, количественно выраженный в величине RTT, может достигать 90 % для сигналов треугольной и гауссовой формы для СКО шума, равного стп = 0,1 (см. табл. 2 и 3). Однако при малых отношениях сигнал-шум (ОСШ) устранение отсчетов, лежащих на границах главной треугольной области, уже не приводит к ощутимому выигрышу точности биспектралной оценки.

3. Показатель RTT ухудшается при увеличении СКО шума, так как вклад флуктуационных ошибок в оценку биспектральной плотности, когда помеха распределена на всей главной треугольной области определения биспектра, становится превалирующим.

4. В целом, удается уменьшить ошибки оценки биспектральной плотности с помощью вычитания среднего значения сигнала, что демонстрируется соответствующим сравнением данных табл. 1 и 4, 2 и 5, а также 3 и 6. Однако данный подход к улучшению биспектральной оценки результативен при малых величинах СКО шума, а с его ростом показатель выигрыша RTT незначительно отличается от единицы.

5. Показатель выигрыша RTT растет медленнее с увеличением длительности сигнала Т для оценки биспектральной плотности сигнала с отличным от нуля средним значением по сравнению с биспектральной оценкой сигнала с нулевым средним значением. Этот вывод вытекает из сравнения результатов, приведенных в табл. 1 - 3, с результатами, представленными в табл. 4 - 6. Более того, величина выигрыша RTT, соответствующая биспектральной оценке сигналов с нулевым средним значением, оказывается меньше величины RTT, рассчитанной для сигналов с отличным от нуля средним значением.

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

дисперсии е (p,q) для импульсов разной фор-

p q

мы, но одинаковой энергии.

Таблица 7

Нормированная полная суммарная дисперсия

СКО шума Прямоу- гольная Треугольная Гауссова

0,1 0,009 0,005 0,006

0,2 0,045 0,029 0,0307

0,3 0,154 0,112 0,118

0,4 0,419 0,279 0,297

0,5 1,0 0,743 0,764

РИ, 2011, № 1

7

Таблица 8

Результаты оконного сглаживания, среднее значение

сигнала не равно нулю

TOSD T CTn

0,1 0,3 0,5

2 849 5,85 103 2,01 -104

6 1,89103 9,17 103 2,71 -104

10 2,50-103 1,22 -104 3,2 -104

TRSD 2 728 5,03 103 1,74 104

6 1,39103 7,13 103 2,24 -104

10 1,58103 8,48 103 2,40 -104

Анализ данных табл. 7 позволяет сделать вывод о том, что эффект растекания биспектра проявляется наиболее ярко, во-первых, для импульса прямоугольной формы и, во-вторых, для максимальной величины СКО шума.

Для оценки эффективности уменьшения растекания биспектра с помощью оконного сглаживания рассмотрим оценки биспектральной плотности с использованием оконного сглаживания в окне Ханна вида w(i)=1/2[1-cos (2р(і-І/2- 1)/I)].

Результаты компьютерного моделирования для примера тестового сигнала - импульса прямоугольной формы - представлены в табл. 8. Сравнение данных, приведенных в табл. 8 и 1, показывает, что использование окна Ханна приводит лишь к незначительному улучшению точности биспектральной оценки. Данное улучшение может достигать значения, равного приблизительно 6 % для малых величин СКО шума, и улучшение практически не наблюдается для больших величин СКО шума.

4. Выводы

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

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

Литература: 1. Бендат Дж. Применение корреляционного и спектрального анализа: Пер. с англ./Дж. Бендат,

А. Пирсол. М.: Мир, 1983. 312 с. 2. Totsky A.V. Time-frequency analysis of radar backscattered signals using phase coupled frequencies extracted from time-varying bispectrum estimates/ A. V. Totsky, I. V. Kurbatov, G.I. Khlop-ov, S.I. Khomenko, V.Ye. Morozov, J.T. Astola, K.O. Egiaz-arian// Proc. 16th International Conf. Microwaves, Radar and Wireless Communications MIKON-2006, May 22-26, 2006. Poland, Krakow. Vol. I. P. 203-206. 3. Astola J.T. Time-frequency analysis of ground surveillance Doppler radar echo signals by using short-time cross-bispectrum esti-mates/J.T. Astola, K.O. Egiazarian, G.I. Khlopov, S.I. Khomenko, I.V. Kurbatov, V.Ye. Morozov, A.L. Tepliuk, A.V. Totsky]// Proc. of International Radar Symposium IRS 2007, September 5-7, 2007. Cologne, Germany. P. 805-808. 4. Astola J. T. Application of bispectrum estimation for time-frequency analysis of ground surveillance Doppler radar echo signals /J.T. Astola, K.O. Egiazarian, G.I. Khlopov, S.I. Khomenko, I.V. Kurbatov, V.Ye. Morozov, A.V. Totsky // IEEE Transact. on Instrumentation and Measurement. 2008. Vol. 57. P. 19491957. 5. LohmannA. W. Speckle masking in astronomy: Triple correlation theory and applications / A. W. Lohmann, G. Weigelt, B. Wirnitzer// Applied Optics. 1983. Vol. 22. P. 40284037. 6. Тоцкий А.В. Распознавание летательных аппаратов по радиолокационным дальностным профилям/ А.В. Тоцкий, П.А. Молчанов, Б.Б. Поспелов// Авіаційно-космічна техніка і технологія, 5(72). 2010. С. 77-82. 7. Molchanov P.A. Application of triple correlation and bispectrum for interference immunity improvement in telecommunications systems/ P.A. Molchanov, A.V. Totsky// International Journal of Applied Mathematics and Computer Science. 2008. Vol. 18. № 3. P. 361-367. 8. Nikias C.L. Bispectral estimation: A digital signal processing framework/ C.L. Nikias, M.R. Raghuveer // Proc. IEEE. 1987. Vol. 75, No. 7, July. P. 869 -891. 9. NakamuraM. Waveform estimation from noisy signals with variable signal delay using bispectrum averaging/ M. Nakamura // IEEE Trans. on Biomedical Engineering. 1993. Vol. 40, No. 2. P. 118-127.

Поступила в редколлегию 16.01.2011

Рецензент: д-р физ.-мат. наук Сугак В.Г.

Тоцкий Александр Владимирович, д-р. техн. наук, профессор кафедры «Прием, передача и обработка сигналов» НАКУ им. Жуковского «ХАИ». Адрес: Украина, 61070, Харьков, ул. Чкалова, 17, тел. (057) 707-43-52, e-mail: totskiy@xai.edu. ua

Науменко Виктория Владимировна, аспирантка кафедры «Прием, передача и обработка сигналов» НАКУ им. Жуковского «ХАИ». Адрес: Украина, 61070, Харьков, ул. Чкалова, 17, тел. (057) 707-43-52, e-mail:

naumenko_vita@mail.ru

8

РИ, 2011, № 1