Изобретение относится к гидроакустике, радиолокации и статистической радиотехнике, а именно к способам обнаружения сигналов в условиях реального флуктуирующего канала.
Реализация оптимального приема при решении задач обнаружения сигналов во многом определяется уровнем априорных знаний о принимаемом сигнале.
Если принимаемый сигнал полностью известен, то оптимальное обнаружение сигнала на фоне белого шума обеспечивает согласованная фильтрация или коррелятор, реализующие корреляционное сравнение принимаемого сигнала с копией. Однако корреляционным способам обнаружения присущ основной недостаток: в условиях реального флуктуирующего канала распространения происходит раскорреляция принимаемого сигнала и копии вследствие флуктуации параметров сигнала и как результат, значительная потеря помехоустойчивости, при этом параметры становятся случайными величинами и точно неизвестными. В случае неизвестных параметров прием ведется без учета конкретных их значений путем статистического усреднения принимаемого сигнала.
Если известны только распределения сигнала и помехи, то используют энергетические способы обнаружения, представляющие собой последовательное выполнение операций фильтрации, детектирования и интегрирования.
Недостатком таких методов является "эффект подавления малого сигнала", так как выходное отношение сигнал-помеха пропорционально входному отношению сигнал/помеха.
Для сигналов с неизвестной начальной фазой оптимальной является квадратурная обработка. Характеристики сигнал/шум свидетельствуют о проигрыше по сравнению с коррелятором, вызванный случайной начальной фазой, но сам проигрыш невелик и составляет от 1,1 до 1,34 дБ по мощности (плата за незнание фазы). Условие, при котором находится отношение правоподобия, для квадратурного приемника строится на гипотезе постоянства начальной фазы на интервале наблюдения. Данная гипотеза выполнима в случае, когда интервал постоянства начальной фазы принимаемого сигнала больше или равен времени наблюдения. Последнее проиллюстрировано на фиг.1 (интервал 1), однако известно, что при увеличении дистанции до цели для обнаружения эхо-сигналов с заданной помехоустойчивостью, необходимо увеличивать время наблюдения (интервал 2), но это приводит к флуктуации начальной фазы на интервале наблюдения за счет влияния канала распространения и, как результат, нарушению принятой гипотезы постоянства фазы к значительному снижению эффективности квадратурной обработки.
Рассмотрим процессы распpостранения сигнала в реальном флуктуирующем канале со случайными неоднородностями. Масштабы временных и пространственных изменений в канале меняются в чрезвычайно широких пределах, что приводит к флуктуациям сигнала, в большей степени этому подвержена фаза. Например, для акустического канала на частотах f 4-8 кГц среднеквадратичные значения флуктуации фазы составляют от 1,6 до 3 циклов.
При ограниченной пиковой мощности сигнала для увеличения дальности обнаружения, отраженных от цели сигналов, необходимо увеличивать длительность излучаемого импульса. Вследствие этого растет и время наблюдения. Поэтому на больших дистанциях интервал времени, в течение которого фаза может считаться случайной, но постоянной величиной, будет меньше, чем время наблюдения и поэтому применение квадратурной обработки в ее классической постановке приводит к значительному снижению помехоустойчивости оптимального приемника сигналов с неизвестной начальной фазы. Сама начальная фаза становится флуктуирующей на интервале наблюдения, причем, функция значений начальной фазы является гладкой и медленно меняющейся фиг.1 (кривая 1).
Целью изобретения является повышение помехоустойчивости квадратурного приемника в условиях, когда изменения начальной фазы на интервале наблюдения можно аппроксимировать линейным законом.
Поставленная задача достигается тем, что принимается гипотеза линейного изменения фазы на интервале наблюдения фиг.1 (интервал 2). В этом случае может быть построен оптимальный обнаружитель, где в квадратурном приемнике дополнительно реализованы два канала, в которых поводятся операции взаимокорреляционного сравнения входной реализации сигнала с квадратурами вида: в одном (V1(t)t sinω t) в другом с (V1(t)t cos ω t), затем проводится сложение и вычитание полученных корреляционных откликов с корреляционными откликами квадратурного приемника.
Проведем анализ оптимального алгоритма (правила) в случае приема сигнала с флуктуирующей неизвестной начальной фазой. Сигнал можно представить моделью
V1(t, ϕ (t)) V1(t)cos( ωo t + ϕ (t)) (1)
Флуктуирующая случайная начальная фаза сигнала считается равномерно распределенной величиной
P(ϕ(t)) / при 0≅ ϕ(t)<2π
(2)
Тогда условное отношение правдоподобия можно записать
Λ (V1, ϕ)
exp V
exp - + 2 ,
где
I1 V
V
Первый член выражения (4) много больше второго. Поэтому
I1≈ E1, E1= V
Второй интеграл выражения (3)
I2 U1(t) V1(t1ϕ(t))dt U1(t) V1(t) cos(ωot +ϕ(t))dt
U1(t)V1(t)cos ωot cos ϕ(t)dt U1(t)V1(t)sinωot sinϕ(t)dt
Как указывалось выше, функция ϕ (t) гладкая и медленно меняющаяся на интервале наблюдения, то cos ϕ (t) и sin ϕ (t) представим рядом Тейлора, ограничиваясь двумя первыми членами разложения
cosϕ (t) cos ϕo tϕoI sin ϕo
sinϕ (t) sin ϕo + tϕoI cosϕo (5)
тогда
I2= U1(t)V1(t) cosωot[cosϕo- t sin ϕo] dt
U1(t)V1(t) sinωot[ sin ϕo+ t cosϕo]dt
U1(t)V1(t) cosωot cosϕodt tU1(t)V1(t) cosωot sin ϕodt-
U1(t)V1(t) sin ωot sinϕodt tU1(t)V1(t) sinωot cosϕodt
cosϕoU1(t)V1(t)cosωot dt sin ϕotU1(t)V1(t) cosωot dt
sinϕoU1(t)V1(t)sinωot dt cos ϕotU1(t)V1(t) sinωot dt
Обозначая
X1= U1(t)V1(t)cosωot dt; X2= U1(t)V1(t) sinωot dt
Y1= U1(t)V1(t)sinωot dt; Y2= U1(t)V1(t) cosωot dt,
получаем
I2 cosϕo X1- ϕoI sinϕo Y2-sin ϕo Y1- -ϕoI cos ϕo X2 cosϕo [X1- ϕoI X2] sin ϕo -[Y1 + ϕoI Y2]
(7)
Фаза и ее производная независимы, причем, производная фазы распределена по нормальному закону с математическим ожиданием равным нулю. Функция медленно меняющаяся, гладкая и ее производная мала и, вводя обозначения в (7), получаем
X1- ϕoI X2 X
Y1 + ϕoI Y2 Y (8)
В выражениях (8) первые значения X1 и Y1 много больше, чем произведение производной фазы на Х2 и Y2 следовательно, производная не оказывает существенного влияния на результат, тогда имеем (векторная диаграмма квадратурных составляющих сигнала и шума показана на фиг.2).
I2 Xcos ϕo- Ysin ϕo Zcos( ϕo + Ψ ),
(9) где
(10)
Следовательно
Λ (U1,ϕ) exp - +
Безусловное отношение правдоподобия
Λ (U1) (U1, ϕ) ρ(ϕ)dϕ exp - + cos (ϕo+Ψ)dϕ
e dϕ (11)
Так как
excos(ϕo+Ψ) dϕ Jo(x) где Io(x) модифицированная функция Бесселя первого рода нулевого порядка, то
Λ (U1) exp - J (12) и алгоритм оптимального обнаружения сигнала при принятой гипотезе о начальной фазе запишется
Λ (U1) exp- J Λo
(13)
Ввиду монотонности Io(t) можно записать
ln J + ln Λo или Z2 ≥Zo2, где Z корень уравнения
ln J= + ln Λo
(14)
Таким образом, оптимальный обнаружитель вычисляет величину X2 + Y2= Z2 и сравнивает ее с порогом Vo2 Zo2
Следует отметить, что полученный алгоритм обнаружения сигналов реализует (как частный случай) обработку сигнала с неизвестной начальной фазой (классическую квадратурную обработку) при ϕoI 0
Останавливаясь на оценке помехоустойчивости предложенного алгоритма (приемника) нетрудно убедиться, что она незначительно отличается от помехоустойчивости квадратурной обработки в силу того, что фаза в канале флуктуирует медленно на интервале наблюдения и ее производная будет малой величиной и в (9) существенного влияния не оказывает. Следовательно, вероятности правильного обнаружения, ложной тревоги, характеристики обнаружителя оцениваются и получаются так же, как и квадратурном приемнике.
Таким образом, по сравнению с прототипом новизна предлагаемого способа обнаружения сигналов заключается в том, что на основании принятия более общей гипотезы относительно начальной фазы принимаемого сигнала выполняются следующие новые операции:
умножения квадратурных составляющих на текущее время наблюдения;
корреляционного сравнения принимаемой реализации с квадратурами вида (V1(ω)tcos (t) и (V1(t)tsin ωt);
вычитания и суммирования корреляционных откликов квадратурного приемника и вновь сформированных корреляционных откликов.
Дополнительные операции реализуют приведенный выше алгоритм и позволяют обнаруживать сигналы с неизвестной флуктуирующей начальной фазой в условиях реального канала распространения. На фиг.1 показано, что начальная фаза аппроксимирована линейной функцией, в случае, когда значение производной фазы мало на интервале наблюдения (интервал 2) помехоустойчивость оценивается также, как и в случае приема сигналов с неизвестной начальной фазой, при больших значениях производной оценка равносильна оценке помехоустойчивости при приеме сигналов с неизвестной начальной фазой и случайной амплитудой.
Устройство, реализующее предлагаемый способ, показано на фиг.3, где приняты следующие обозначения: 1 аналого-цифровой преобразователь (АЦП), 2 постоянное запоминающее устройство (ПЗУ), 3 первое взаимокорреляционное устройство (1-ВСУ), 4 четвертое взаимокорреляционное устройство (Н-ВКУ), 5 первый умножитель (1-УМ), 6 второй умножитель (2-УМ), 7 второе взаимокорреляционное устройство (В-ВКУ), 8 третье взаимокорреляционное устройство (З-ВКУ), 9 первый сумматор (1-СМ), 10 вычитатель (ВЧ), 11 первый квадратор (1-УВ), 12 второй квадратор (2/КВ), 13 второй сумматор (2-СМ), 14 пороговое устройство (ПУ), 15 устройство управления (УУ).
Устройство работает следующим образом.
На вход устройства происходит реализация U1(t) которая поступает на вход АЦП (1) с частотой дискретизации, удовлетворяющей, в соответствии с теоремой Котельникова, условии ω1 ≥ 2 ωo С выхода АЦП дискретная выборка U1(n) подается на первые входы всех ВКУ (3, 4, 7, 8). На вторые входы всех ВКУ (3, 4, 7, 8) поступает квадратура вида соответственно
1-ВКУ (3) V1(t)sin ωt со второго выхода ПЗУ (2);
2-ВКУ (7) V1(t)t cosω t с выхода первого умножителя (5);
3-ВКУ (8) V1(t)t sinω t с выхода второго умножителя (6);
4-ВКУ (4) V1(t)cosω t с первого выхода ПЗУ (6).
Корреляционные отклики с выходов взаимокорреляционных устройств (3, 7, 8, 4) поступают:
с выхода первого ВКУ/3 на первый вход первого сумматора (9),
с выхода второго ВКУ (7) на второй вход первого сумматора,
с выхода третьего ВКУ (8) на первый вход вычитателя (10),
выходы четвертого ВКУ (4) на второй вход вычитателя (10).
Сумма корреляционных откликов первого и второго ВКУ с выхода первого сумматора (9) подается на вход первого квадратора (11), с выхода которого поступает на первый вход второго сумматора (13).
Разность отливок третьего и четвертого ВКУ с выхода вычитателя (10) подается на вход второго квадратора (12), с выхода которого подается на второй вход второго сумматора (13).
На выходе второго сумматора (13) получается значение которое в пороговом устройстве (14) сравнивается с порогом Zo. Выход порогового устройства 14 является выходом устройства.
Квадратурные составляющие формируются в устройстве следующим образом.
Синфазная составляющая, хранящаяся в ПЗУ 2 со второго выхода одновременно подается на второй вход первого ВКУ 3, и на вход второго умножителя 6, на управляющий вход первого умножителя 5 поступает текущее время наблюдения с третьего выхода ПЗУ 2, с выхода второго умножителя 6 синфазная составляющая вида V1(t)tsin ωt подается на второй вход третьего ВКУ 8. Квадратурная составляющая с первого выхода ПЗУ 2 одновременно подается на второй вход четвертого ВКУ 4 и на вход первого умножителя 5, на управляющий вход первого умножителя поступает текущее время наблюдения с третьего выхода ПЗУ 2. С выхода первого умножителя 5 квадратурная составляющая вида V1(t)tcosω t подается на второй вход второго ВКУ 2. Устройство управления 15 осуществляет синхронизацию работы блоков АЦП 1, всех взаимокорреляционных устройств 3, 4, 7, 8 квадраторов 11, 12 и ПЗУ 2.
название | год | авторы | номер документа |
---|---|---|---|
СПОСОБ ОПРЕДЕЛЕНИЯ СКОРОСТИ ДВИЖЕНИЯ СУДНА ОТНОСИТЕЛЬНО ДНА | 1992 |
|
RU2037847C1 |
УСТРОЙСТВО ДЛЯ ПРИЕМА РЕЧЕВЫХ СИГНАЛОВ | 1993 |
|
RU2049426C1 |
СПОСОБ ПЕРЕДАЧИ РЕЧЕВЫХ СИГНАЛОВ | 1993 |
|
RU2049456C1 |
СПОСОБ ОБНАРУЖЕНИЯ СИГНАЛОВ | 1992 |
|
RU2249833C2 |
СПОСОБ ПОЛУЧЕНИЯ СЕЙСМИЧЕСКОГО РАЗРЕЗА | 1991 |
|
RU2022310C1 |
ФАЗОВЫЙ ИЗМЕРИТЕЛЬ ДАЛЬНОСТИ С УСТРОЙСТВОМ ИЗМЕРЕНИЯ ДИСПЕРСИИ ФАЗОВОГО СДВИГА | 1992 |
|
RU2048679C1 |
СПОСОБ ВОССТАНОВЛЕНИЯ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ СЛУЧАЙНЫХ ВОЛНОВЫХ НЕОДНОРОДНОСТЕЙ СРЕДЫ, МЕНЯЮЩИХСЯ ВО ВРЕМЕНИ | 1995 |
|
RU2099690C1 |
СПОСОБ ОПРЕДЕЛЕНИЯ СКОРОСТИ ДВИЖЕНИЯ СУДНА ОТНОСИТЕЛЬНО ДНА | 1992 |
|
RU2042152C1 |
СПОСОБ АКУСТИЧЕСКОГО ПРЕДСТАВЛЕНИЯ ПРОСТРАНСТВЕННОЙ ИНФОРМАЦИИ ДЛЯ ИНВАЛИДОВ ПО ЗРЕНИЮ | 1993 |
|
RU2085162C1 |
УСТРОЙСТВО ОБНАРУЖЕНИЯ ШУМОВЫХ ГИДРОАКУСТИЧЕСКИХ СИГНАЛОВ НА ОСНОВЕ КВАДРАТУРНОГО ПРИЕМНИКА | 2013 |
|
RU2549207C2 |
Изобретение относится к гидроакустике, радиолокации и статистической радиотехнике - к способам обнаружения сигнала в условиях реального канала распространения. Цель изобретения - повышение помехоустойчивости квадратурного приемника в условиях, когда изменение начальной фазы на интервале наблюдений можно апроксимировать линейным законом. Цель достигается тем, что принимается гипотеза литейного изменения фазы на интервале наблюдения. Сущность: реализация в квадратурном приемнике дополнительно операций умножения квадратурных составляющих на текущее время наблюдения, взаимно корреляционное сравнение принятой реализации с результатами умножения квадратурных составляющих, вычитание и суммирование корреляционных откликов квадратурного приемника и вновь сформированных коррелляционных обкликов. Схема реализации заявленного технического решения содержит аналого-цифровой преобразователь (АЦП) 1, постоянно запоминающее устройство (ПЗУ) 2, первое взаимокорреляционное устройство 3, четвертое взаимокорреляционное устройство 4, первый умножитель 5, второй умножитель 6, второе взаимокорреляционное устройство 7, четвертое взаимокорреляционное устройство 8, первый сумматор 9, вычитатель 10, первый квадратор 11, второй квадратор 12, второй сумматор 13, пороговое устройство 14, устройство управления 15. 3 ил.
СПОСОБ ОБНАРУЖЕНИЯ СИГНАЛОВ, включающий операции взаимокорреляционного сравнения принятой реализации с квадратурными составляющими сигналов, возведения в квадрат корреляционных откликов результатов сравнения принятой реализации и квадратурных составляющих и суммирования результатов возведения в квадрат, отличающийся тем, что, с целью повышения эффективности обнаружения сигналов с неизвестной начальной фазой, изменение которой на интервале наблюдения аппроксимируется линейным законом, дополнительно проводят операции умножения квадратурных составляющих на текущее время наблюдения, взаимокорреляционного сравнения принятой реализации с результатами умножения квадратурных составляющих на текущее время наблюдения, вычитания корреляционного отклика результата сравнения принятой реализации и синфазной составляющей, умноженной на текущее время, из корреляционного отклика результата сравнения принятой реализации и квадратурной составляющей и суммирования корреляционного отклика результата сравнения принятой реализации и синфазной составляющей с корреляционным откликом результата сравнения принятой реализации и квадратурной составляющей, умноженной на текущее время наблюдения, по результатам сравнения полученного значения с порогом принимают решение об обнаружении сигнала.
Лезин Ю.С | |||
Введение в теорию и технику радиотехнических систем | |||
М., Радио и связь | |||
Пневматический водоподъемный аппарат-двигатель | 1917 |
|
SU1986A1 |
САННЫЙ ВЕЛОСИПЕД С ВЕДУЩИМ КОЛЕСОМ, СНАБЖЕННЫМ ШИПАМИ | 1921 |
|
SU265A1 |
Авторы
Даты
1995-04-10—Публикация
1992-02-24—Подача