Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Distantsionny_kontrol_psikhofiziologicheskogo_sostoyania_s_pomoschyu_sverkhshirokopolosnoy_RLS

.pdf
Скачиваний:
25
Добавлен:
04.05.2022
Размер:
12.15 Mб
Скачать

411

Рис. 5.54. Sкоррc c ( q ), = 0, с = /2.

Рис. 5.55. Sкоррs c ( q ), = 0, с = /2

Рис. 5.56. Спектр сигнала рис. 5.51.

412

Рис.5.57. Спектр сигнала рис. 5.52

Рис. 5.58. Sкоррc c ( q ), = 0, с = 3 /4.

Рис. 5.59. Sкоррs c ( q ), = 0, с = 3 /4

Рис. 5.60. Спектр сигналов рис. 5.58 и 5.59.

Из рис. 5.47 - 5.60 можно сделать следующие основные выводы.

1. Амплитуда, форма и, следовательно, спектральный состав квадратурных компонент сигна-

лов сильно зависят от случайной начальной фазы с. Напомним, что в с включены фаза отражения и

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

личина с случайна и равномерно распределена в интервале ( , ), то при каждом сеансе радиолока-

ционного наблюдения длительностью T фаза с может принимать любое значение. Это недостаток

413

обеих квадратурных компонент сигналов. В то же время, как следует из (5.47), модуль комплексного сигнала не зависит от с. Таким образом, исключить зависимость сигнала от с можно путем форми-

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

ния поверхности тела , как видно из (5.47) (5.49), приводит только к сдвигу сигнала по оси време-

ни без изменения его формы.

2. Форма сигнала не совпадает с формой колебания поверхности тела. При принятом в (5.45)

гармоническом колебании поверхности тела спектр полезного сигнала содержит несколько гармоник.

При частоте колебания поверхности тела F = 0.25 Гц спектр полезного сигнала простирается до 1 Гц.

Можно показать, что отличие формы полезного сигнала от гармонической и, соответственно, расши-

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

длине волны 0.

3. Квадратурные компоненты полезного сигнала содержат постоянные составляющие (за ис-

ключением некоторых частных случаев). Поскольку их нельзя отделить от постоянных составляю-

щих, вызванных эхосигналами от мешающих неподвижных отражателей, их приходится исключать

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

ским потерям сигналов.

Оптимизация обработки после когерентного накопления

В соответствии с общим алгоритмом обнаружения следующая операция после когерентного накопления на М интервалах (длительностью t каждый) заключается в когерентном подавлении пас-

сивных помех в полученных отсчетах суммы собственных шумов, пассивных помех и (возможно)

полезного сигнала каждой квадратуры. Обозначив

yc

и

ys

 

 

косинусную и синусную компонен-

 

 

 

 

 

 

корр q

 

 

корр q

 

ты q-ого отсчета после когерентного накопления (5.40), получим (с заменой q 1 на q 0.5 )

 

 

A

n 1

 

[ t ( q 0.5 ) t kT

0

]2

 

 

 

yкоррc

q =

0

x( t )exp

 

п

 

 

 

 

 

2

 

 

 

 

 

 

2 k 0

 

и

 

 

 

 

 

 

 

cos 0 [ t ( q 0.5 ) t kTп 0 )] dt;

(5.50)

414

 

 

A

n 1

 

 

 

[ t ( q 0.5 ) t kT

0

]2

 

yкоррs

q =

0

x( t )exp

 

 

 

п

 

 

 

 

 

2

 

 

 

 

 

2 k 0

 

 

 

и

 

 

 

 

 

sin 0 [ t ( q 0.5 ) t kTп 0 )] dt.

(5.51)

Режекция пассивных помех от неподвижных отражателей сводится к вычитанию постоянной составляющей каждой квадратуры из всех отсчетов соответствующей квадратуры. В результату по-

лучим

 

M

 

 

M

 

 

zqc yкоррc

q (1 / M ) yкоррc

i ;

zqs yкоррs

q (1 / M ) yкоррs

i ; q = 1, 2,.. , M.

(5.52)

 

i 1

 

 

i 1

 

 

На рис. 5.61 5.64 для примера показаны отсчеты полезного сигнала после вычитания посто-

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

Рис. 5.61. Sкоррc c ( q ), = 0, с = 0.

Рис. 5.62. Sкоррs c ( q ), = 0, с = 0

415

Рис 5.63. Спектр сигнала на рис. 5.61(после вычитания постоянной составляющей).

Рис. 5.64. Спектр сигнала на рис. 5.62 (после вычитания постоянной составляющей).

Сравнение рис. 5.61 с рис 5.47 и рис. 5.63 с рис. 5.59 показывает, что исчезла постоянная со-

ставляющая (нулевая частота спектра), но сохранилась переменная составляющая с удвоенной часто-

той колебания поверхности тела (0.5 Гц), а также небольшая вторая гармоника этой частоты. Но так как амплитуда переменной составляющей на рис. 5.47 более чем в 2 раза меньше максимально воз-

можной (единицы), то исключение постоянной составляющей привело к значительной потере энер-

гии косинусной квадратурной компоненты.

Наоборот, в синусной компоненте на рис. 5.48 нет постоянной составляющей (см. также рис.

5.50). Поэтому ее исключение в результате операции (5.52) не изменило ни сам сигнал (рис.5.62), ни его спектр (рис. 5.64).

После исключения постоянной составляющей (режекции пассивных помех) необходимо вы-

числить модуль для каждой пары значений квадратур zqc и zqs , q = 1,…, M, а затем выполнить неко-

герентное накопление (суммирование) всех М модулей. Как следует из общей теории радиолокации,

416

при некогерентном накоплении и сильных сигналах следует вычислять модуль, а при слабых квад-

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

дует ожидать, что, в основном, придется иметь дело со слабыми сигналами (малыми значениями от-

ношения сигнал/шум). Поэтому вычисляем квадрат модуля

Uq = ( zqc )2 ( zqs )2 ,

q = 1, 2,..., M.

(5.53)

Наконец, получаем выходную статистику, подлежащую сравнению с порогом для решения за-

дачи обнаружения:

M

L = Uq .

q 1

(5.54)

На этом синтез алгоритма обнаружения, в принципе, закончен, т.е. определена структура обна-

ружителя (см. рис. 5.65).

Опорн. сигнал

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Корр.1

 

 

Когер.

 

 

Режект.

 

 

Квадра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

фильтр

 

 

тор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Некогер.

 

 

Порог

 

От

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

МШУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Корр.2

 

 

Когер.

 

 

Режект.

 

 

Квадра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

фильтр

 

 

тор

 

 

 

 

 

 

Да

Нет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Опорн. сигнал

cos

Рис. 5.65. Структура обнаружителя.

Однако осталась не определенной длительность каждого интервала когерентного накопления

t, а, следовательно, число когерентно накапливаемых импульсов n и число некогерентно накаплива-

емых отсчетов М. При частоте повторения Fп зондирующих импульсов n = Fп t, а М = Т/ t. В рас-

смотренных примерах принималось t = 0.1 с. Тогда М =200 и, например, при Fп = 15 МГц n

=1.5*106. Как будет показано далее, в некоторых случаях для удовлетворения требований электро-

417

магнитной совместимости (а также для упрощения технической реализации обработки в цифровом виде) частоту повторения целесообразно уменьшать до Fп = 2…3 МГц. В этих случаях n = (2…3)*105

импульсов, что при том же значении t приводит к снижению отношения сигнал/шум на 8.8…7.0 дБ.

Увеличение t позволяет повысить коэффициент когерентного накопления n, но приводит к сглажи-

ванию полезной переменной составляющей полезного сигнала. Кроме того, уменьшается также число некогерентно накапливаемых отсчетов М =Т/ t. Следует оптимизировать значения t и М с целью достижения наилучших характеристик обнаружения полезного сигнала.

Проведем сначала анализ полученного алгоритма, полагая, как и до сих пор, что интервал ко-

герентного накопления t достаточно мал, так что сглаживанием полезного сигнала можно прене-

бречь, а затем рассмотрим задачу оптимизации t и М.

5.3.5. Анализ полученного алгоритма обнаружения

Для того, чтобы построить характеристики обнаружения полученного алгоритма, следует, как обычно, найти распределения вероятностей выходной статистики L в (5.54) при отсутствии и при воздействии полезного сигнала. В отсутствие полезного сигнала на входе приемника действуют соб-

ственные шумы и пассивные помехи от неподвижных отражателей. Возможно и воздействие разного рода активных помех в рабочем диапазоне частот РЛС, однако их здесь учитывать не будем6).

Обнаружение сигналов на фоне собственных шумов в отсутствие сглаживания сигнала при когерентном накоплении.

Предположим, что пассивные помехи полностью подавлены в результате исключения посто-

янных составляющих из отсчетов квадратурных компонент (5.52).

Отсчеты квадратур шума после когерентного накопления (см. (5.50) и (5.51)) представляют со-

бой случайные гауссовские величины с нулевыми средними.

n 1

yкоррc = A0 k 0

[ t ( q 0.5 ) t kT

0

]2

 

n( t )exp

 

п

 

 

 

2

 

 

 

и

 

 

 

cos 0 [ t ( q 0.5 ) t kTп 0 ]dt ,

6)В реальном устройстве для защиты от активных помех применяется полосовой фильтр на входе приемника, рассчитанный на рабочий диапазон частот РЛС, а также полосовой фильтр после режекции пассивных помех, пропускающий только спектр полезного сигнала, например, в пределах от 0.15 Гц до 1.05 Гц в соответствии с

рис. 5.47 5.64.

418

yкоррs

sin 0

n 1

 

 

[ t ( q 0.5 ) t kTп 0 ]

2

 

 

= A0 n( t )exp

 

 

 

 

 

 

 

2

 

k 0

 

 

и

 

 

 

[ t ( q 0.5 ) t kTп 0 ]dt .

Дисперсии вычисляются, например, для косинусной составляющей, по формуле:

 

 

 

 

 

 

n 1 n 1

 

 

 

 

 

 

[ t ( q 0.5 ) t kT

0

]2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( yкоррc

)2 = A02 n( t1 )n( t2

)exp

 

1

 

п

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

k 0 p 0

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

[ t ( q 0.5 ) t pT ]2

 

 

( q 0.5 ) t kTп 0

]

exp

2

 

2

п

0

cos 0 [ t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

cos 0 [ t2 ( q 0.5 ) t pTп 0 ]dt1dt2 .

(5.55)

Черта сверху означает усреднение по множеству реализаций.

Учтем, что входные шумы белые, так что в (5.55)

n( t 1 )n( t2 ) ( N0 / 2 ) ( t 1 t2 ) ,

(5.56)

где N0 односторонняя спектральная плотность белого шума на входе приемника. Подставив

(5.56) в (5.55) и проинтегрировав, например, по t2, получим однократный интеграл от произведения экспонент и косинусов. Заметим также, что интегралы от косинуса с частотой 2 0 равны нулю. Заме-

ним переменные аналогично тому, как это сделано в (5.41), и учтем, что все слагаемые в (5.55) равны нулю при (k p) 0. В результате

( yкоррc )2 = ( yкоррc ш )2 = A02 иN0n/4 2 ,

(5.57)

где n, как и выше, коэффициент когерентного накопления. Такой же результат получим и для синусной составляющей:

( yкоррs )2 = ( yкоррc ш )2 = A02 и N0n/4 2 .

(5.58)

Видно, что благодаря стационарности шума дисперсии шума не зависят от номера q.

При воздействии сигнала дисперсия не изменяется, но появляются ненулевые средние значе-

419

ния (см. (5.48), (5.49)):

 

 

 

 

 

 

 

 

 

 

 

A

n

 

 

 

2 { R

R sin[( q 0.5 ) t ]}2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yкоррc

 

 

 

Sкоррc

 

=

0

и

 

 

Ac

exp

 

0

 

 

 

 

 

q ( c ш )

c q

 

 

 

 

 

 

2

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

и

 

 

 

 

 

4R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

sin[( q 0.5 ) t ] c ;

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.59)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

n

 

 

 

2 { R R sin[( q 0.5 ) t ]}2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yкоррs

 

 

Sкоррs

 

=

0

и

 

 

Ac

exp

 

0

 

 

 

 

 

q ( c ш )

c q

 

 

 

 

 

 

2

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

и

 

 

 

 

4R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

sin[( q 0.5 ) t ] c .

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.60)

Отношение сигнал/шум по мощности, как следует из (5.57) (5.50),

 

 

(Sкоррc

c q )2

 

A2

и

n

 

4 { R

R sin[( q 0.5 ) t ]}2

 

 

(Qc)2

=

 

 

 

 

=

c

 

 

exp

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

( yкоррc

q ш )2

 

 

 

 

 

 

 

 

 

 

N0

2

 

 

 

c

и

 

 

cos2 4Rsin[(

0

(Qs) 2

( Sкоррs

c q )2

=

 

 

 

 

 

( yкоррs

q ш )2

q 0.5 ) t ] п ;

 

A2

n

 

4 { R R sin[( q 0.5 ) t ]}2

 

=

c

и

 

 

exp

0

 

 

N0

 

 

 

 

2

2

 

 

 

 

 

2

 

 

c

и

 

sin2

 

4R

 

 

 

sin[( q 0.5 ) t ] п .

0

 

 

 

Учитывая, что для каждого колокольного импульса из последовательности эхосигналов (5.39),

длительность которого и определена по уровню exp( /4) 0.46,

A2

 

 

 

= Еи, Еи 0.71 Ри и,

и

/2 2

c

 

 

 

(5.61)

где Еи энергия импульса эхосигнала от «полезного» отражателя на входе приемника, Ри

мощность в максимуме импульса, находим

( Qс )q2

2Е n

 

4 { R

R sin[( q 0.5 ) t ]}2

 

 

и

exp

0

 

 

 

 

N0

 

2

2

 

 

 

c

и

 

 

420

 

4R

 

 

 

 

 

 

cos2

 

 

sin[( q 0.5 ) t ] п ;

 

 

 

0

 

 

 

 

 

 

 

 

 

 

(5.62)

 

 

 

 

 

 

 

 

 

( Qs )q2

 

2Е n

 

4 { R R sin[( q 0.5 ) t ]}2

 

 

=

и

 

exp

0

 

 

 

 

N0

 

 

2

2

 

 

 

 

 

c

и

 

 

sin2

 

4R

 

 

 

sin[( q 0.5 ) t ] п .

0

 

 

 

(5.63)

Видно, что благодаря когерентному накоплению сигнала и некогерентному накоплению шумов на интервалах длительностью t, отношение сигнал/шум по мощности выросло в n раз в каждом q-ом отсчете.

Следующей операцией в соответствии с синтезированным в подразделе 5.3.4 алгоритмом явля-

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

величины с нулевыми средними значениями. Однако дисперсии zqc и zqs (5.52) отличаются от (5.57) и

(5.58). Например, для косинусной компоненты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

c

 

 

 

 

c

 

c

 

 

c

 

 

 

 

c

 

 

 

c

 

 

 

 

 

 

 

zq

2

= z

2

 

yкорр qш

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

=

+ 1 / M

yкорр iш

2 yкорр qш 1 / M yкорр iш .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

i 1

 

В суммах М 1 слагаемых не зависят от yc

 

(при i q), но одно слагаемое при i = q совпа-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

корр qш

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дает с yкоррc

. Поэтому, учитывая, что

yкоррc

2

=

yкоррc

ш 2

(не зависят от q),

 

 

 

 

 

 

 

+

 

 

 

 

 

 

M

 

 

 

 

 

 

zqc 2

=

yкоррc

ш 2

 

yкоррc

ш 2

M 2

 

yкоррc

ш 2

yкоррc

ш 2

11 M .

 

 

Аналогичный результат можно получить для синусной компоненты. Учитывая (5.57) и (5.58),

имеем

zqc 2 = zqs 2 = A02 и N0n 11 M 42 .

(5.64)