Основы теории радиосистем и комплексов радиоэлектронной борьбы. Учебное пособие
.pdf
31
xi |
Z 1 |
Z 1 |
|
||
x 2 |
x1 |
x0 |
W0 |
W1 |
W 2 |
y
Рис. 3.1. Пример структуры линейного фильтра второго порядка
В фильтре, приведенном на рис. 3.1, выходная величина y формируется путем суммирования попарных произведений:
|
|
|
|
|
2 |
|
2 |
|
|
|
|
|
y W0 x2 W1x1 W2 x0 Wk x2 k |
xkW2 k , |
|
|
|
||||||
|
|
|
|
|
k 0 |
|
k 0 |
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
где W ,W ,W WT , |
|
|
|
|
T |
|
|
|
, и |
||
x , x , x xT . Тогда xT |
x x |
|
|||||||||
0 |
1 |
2 |
0 |
1 |
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
в результате мы получим, что WT x W x |
2 k |
. |
|
|
|
|
|||||
|
|
|
|
|
k |
|
|
|
|
|
|
k 0
Среднеквадратическое значение
2 S WT x S WTx T S2 2SWTx WTxxTW,
где xxT Rx соответствует корреляционной матрице сигнала, а
SxT Rxs – взаимной корреляционной матрице сигнала и шума.
Условие минимума среднего квадрата ошибки имеет вид:2 min и позволяет найти оптимальный вектор Wopt весовой обработки, минимизируя градиент:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
0 |
W , |
|
2 |
|
2R |
|
2R |
W 0 . |
||||||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
W |
|
|
|
|
|
|
|
W |
|
|
|
|
|
|
xs |
|
x |
|
|
|
Отсюда вектор |
W |
R xs |
, R |
|
R |
W , R |
1 |
|
R |
1 |
W . |
||||||||||
|
xs |
x |
R |
xs |
x |
R |
|||||||||||||||
|
|||||||||||||||||||||
|
|
|
|
|
|
R x |
|
|
|
x |
|
|
|
|
x |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
32
В результате получаем выражение Wopt Rx1Rxs , известное как уравнение Винера - Хопфа. Это же уравнение может быть записа-
|
|
|
|
|
|
|
t d Rxs t и из- |
|||||
но в непрерывном времени как |
hopt Rx |
|||||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
вестно как уравнение Фредгольма 1-го рода. |
|
|
|
|
||||||||
|
Согласно алгоритму Винера - Хопфа построим структурную |
|||||||||||
схему (см. рис. 3.2) |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
|
ˆ |
|
ˆ 1 |
|
|
|
ЛФ |
Wopt |
|||
S |
|
|
|
|
|
|||||||
|
|
R x |
|
R x |
|
|
|
|
|
|||
ˆ
R xs
Рис. 3.2 - Структурная схема оптимального фильтра
3.2. Оптимальная линейная фильтрация по критерию максимума отношения сигнал-шум
Рассмотрим линейную задачу, при которой условием синтеза является максимум мощности сигнала по отношению к шуму на выходе фильтра.
Сигнал на выходе фильтра соответствует y WTx , а величи-
на y2 эквивалентна мощности сигнала 2 |
D |
. Тогда |
|||||
|
|
|
|
x |
|
x |
|
|
|
|
|
вых |
|
вых |
|
|
|
|
|
|
|
|
|
y2 WTxxTW WT xxT W WTR |
W . |
|
|||||
|
|
|
|
x |
|
|
|
R x
Считаем, что вектор W не меняется в процессе обработки и, производя замены X S, Rx Rs , получаем значение мощности сигнала на выходе:
S2вых WTRs W , Rs Rп ш2 вх I ,
где I – единичная матрица.
Коэффициент улучшения (изменения) отношения сигнал-шум на выходе линейного фильтра
33
|
|
|
s2 |
|
|
WTR |
W |
|
|
|
|||
|
|
|
|
вых |
|
s |
|
. |
|
|
|||
|
|
2 |
|
|
WT W |
|
|
|
|||||
|
|
|
ш |
вых |
|
|
|
|
|
|
|
||
Для получения максимума этого выражения решим матричное |
|||||||||||||
уравнение: |
|
|
|
|
|
|
|
|
|
|
|
|
|
det R |
|
I |
0 |
2s |
R |
|
W . |
||||||
s |
2 |
s |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь W Wopt - |
|
вектор согласованного фильтра, который |
|||||||||||
обеспечивает максимальное отношение сигнал-шум.
В аналоговом виде частотная характеристика такого согласованного фильтра имеет вид:
Fсф j CSs j e j t .
Согласованный фильтр имеет частотную характеристику, комплексно-сопряженную со спектром сигнала. Амплитудночастотная характеристика (АЧХ) согласованного фильтра совпадает с амплитудным спектром сигнала, а фазочастотная характеристика (ФЧХ) согласованного фильтра имеет противоположный знак.
3.3. Оптимальная линейная фильтрация по критерию максимума отношения сигнал-(шум +помеха)
Рассмотрим линейную задачу выделения сигнала на фоне помехи и шума из входной смеси x(t) S(t) p(t) n(t) , при которой
условием синтеза является максимум мощности сигнала по отношению к сумме помехи и шума на выходе фильтра. Запишем, как и в предыдущем случае, исходные формулы:
|
|
T |
|
|
2 |
|
|
2 |
|
|
||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
||||||
y W |
x , |
вых |
WTR |
W , |
вых |
|
WTR |
W , |
||||
2 |
||||||||||||
|
||||||||||||
|
|
|
|
|
x |
|
2 |
s |
|
|||
|
|
|
|
|
вх |
|
|
вх S |
|
|
||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
T |
R p W |
|
вых |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
вых |
|
W |
|
|
2 |
|
|
|
|
|
|
T |
|
|
|
|||||
|
|
|
|
|
|
|
, |
вх S |
|
|
|
|
|
W |
R s W |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
2 |
|
|
|
T |
W |
|
|
|
|
|
|
|
|
|
|
. |
(3.1) |
|
|
вх |
P |
W |
|
|
|
|
2 |
|
T |
R p W |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
вых |
|
W |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вх P |
|
|
|
|
|
||
34
Числитель и знаменатель последней формулы известны из линейной алгебры как квадратичная форма, соответствующая следующему выражению:
k 1k 1
WT RW wj wl rjl . j 0l 0
Для физически реализуемых процессов эта величина всегда положительная и имеет размерность мощности. Максимум квадратичной формы находится известным методом Лагранжа и сводится к решению проблемы собственных чисел корреляционных матриц сигнала и помех. Если использовать модель нефлуктуирующего в пределах обрабатываемой выборки сигнала, то
2
WT Rs W WTS , j 0, k 1 ,
что соответствует S j S0 S1 ... Sk 1 . Тогда
Wopt R p 1 ST , R p I Wopt I 1ST ST ,
R p 2p R p 2nRn 2p R p 2nI .
Пронормировав суммарную матрицу шума и помехи, получим
R |
|
R |
|
|
2 |
I R |
|
I , |
р n |
p |
n |
p |
|||||
|
|
|
2 |
|
|
|||
|
|
|
|
|
p |
|
|
|
где – отношение шум-помеха.
При увеличении отношения шум-помеха уменьшается степень корреляции смеси помехи с шумом, что отражается на результате оптимизации, приближая параметры оптимального фильтра к параметрам согласованного фильтра [7].
3.4. Оценка параметров сигнала
Рассматриваем ситуацию, при которой процесс x(t) содержит полезный сигнал S(t, α0, α 1…), где αi – параметры сигнала, помеху и белый гауссовский шум. x(t) = S(t, α 0, α 1…) + p(t) + n(t). Пусть математическое ожидание входного процесса M x(t) 0 , а кор-
реляционная функция – R(t1,t2 ) M x(t1), x(t2 ) .
Врезультате обработки требуется найти оценку параметров α i
спредельно возможной точностью. Оценка параметров α i из-за наличия сопровождающих сигнал шумов и помех не может быть
35
вычислена с абсолютной точностью и является случайной величиной. Можно рассматривать точечные или интервальные оценки. Точечные оценки характеризуются в статистическом смысле математическим ожиданием m αi и дисперсией 2 αi. Для несмещенных значений математическое ожидание и истинное значение па-
раметра совпадают m ˆi |
i , в противном случае оценка яв- |
ляется смещенной. Смещение может быть систематическим или случайным.
Интервальные оценки дают значения неизвестного параметра с точностью до ширины интервала и могут характеризоваться вероятностью попадания истинного значения параметра в этот интервал. Точность интервальной оценки определяется шириной ин-
тервала неопределенности параметра (0… ), ( 2 , 2) в соответствии с рис. 3.3.
Рис. 3.3. Интервал неопределенности
Обычно в качестве оценки параметров принимается середина подынтервала, а в качестве меры разброса оценки используется
дисперсия 2 2 , известная из математической статистики как
12
поправка Шеппарда. |
|
|
|
|
|
|
Если рассматривать |
|
|
точечные |
|
оценки, то |
|
M ˆi i , и |
|
условием минимума дисперсии |
||||
ошиб- |
|
|
M ˆ |
|
|
min . |
ки будет выражение: D |
2 |
i |
|
|||
ˆi |
ˆi |
|
i |
|
||
Оценка с минимальной дисперсией называется эффективной оценкой. Нижняя граница дисперсии такой оценки задаётся неравенством Рао - Крамера.
36
|
|
|
|
|
|
|
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
D |
|
|
... |
|
|
ln P |
, x ,...x |
|
P |
x ,...x , dx ,...dx |
m |
, |
|
|
|||||||||||||
ˆi |
|
|
|
m |
1 |
m |
m 1 m |
1 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где Pm(α, x1,…, xm) – совместная плотность вероятности выборочного значения xi и неизвестного параметра α. Знак равенства в формуле соответствует нижней границе Рао - Крамера.
4. Дальность действия РЛС
Дальностью действия радиолокационной станции называется наибольшее расстояние между радиолокационной станцией (РЛС) и целью, на котором ее обнаружение производится с заданными вероятностями правильного обнаружения и ложной тревоги.
Дальность действия зависит от различных факторов: технических параметров станции, характеристик цели, условий распространения радиоволн, наличия и уровня различного рода помех и ряда других факторов, большинство из которых изменяется во времени случайным образом.
4.1. Дальность действия РЛС в свободном пространстве на фоне белого шума
Вначале рассмотрим дальность действия РЛС без учёта влияния земной поверхности и атмосферы на распространение радиоволн, т.е. считаем, что РЛС и цель находятся в «свободном» пространстве. Дальность действия РЛС описывается уравнением дальности, которое определяет взаимосвязь всех основных характеристик РЛС. Считаем также, что мощность излучения (Ризл.) равномерно распределена по всей площади сферы с радиусом R.
Площадь сферы – Sсф = 4 R2. Поток излучаемой мощности и мощность сигнала на входе приемника соответственно равны
Пизл Pизл Pизл2 ,
Sсф 4 R
Pпр Пизл А GAPизл .
4 R2
Последнее уравнение известно как уравнение радиосвязи в свободном пространстве, в котором G – коэффициент усиления
37
антенны, А – коэффициент направленного действия (КНД) антенны. Мощность сигнала на входе приемника обратно пропорциональна дальности Pпр 1/D2.
Для случая пассивного ответа мощность переизлученной энергии связана с эффективной поверхностью рассеивания цели (ЭПР). Если ц(м2) – площадь плоской проводящей поверхности, перпендикулярной к потоку падающей волны, то
П р Пц ц GPизл ц / 4 R2 2
4 R2
GPизл ц / 16 2 R4 П р ~ 1
R4 .
Обозначим коэффициенты направленного действия (КНД) передающей и приемной антенн G1 и G2. При использовании импульсного режима работы локатора обычно используется одна антенна на передачу и прием сигнала и, поэтому G1 = G2 =G. Излучае-
мая мощность Pр = Pпер, = 0.9, Pпр 1/R4, Pпр = ПрАпр = GPр ц
/(4 2R4), G 4 A . Мощность на входе приемника является функ-
2
цией от мощности излучения в абсолютных величинах.
Pпр = G2Pпер ц 2/(4 3R4), Pпр 1/R4. Обозначив Sa – активную площадь антенны, Pпор – пороговую мощность приемника, Sc – эффек-
тивную поверхность рассеивания цели, получим, что максимальная дальность обнаружения в свободном пространстве Rmax0 -
дальность обнаружения цели для РЛС с пассивным ответом определяется как
|
|
|
|
P G G S |
c |
|
, Rmax0 4 |
|
P S |
a1 |
S |
a 2 |
S |
c |
|
|
|
|||
Rmax0 |
4 |
|
p 1 |
2 |
|
p |
|
|
|
. |
|
|||||||||
|
|
|
Pпор |
|
|
|
2 |
Pпор |
|
(4.1) |
||||||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Для РЛС с активным ответом классическими считаются следующие выражения [3]:
38
R |
|
P G G S |
|
|
; R |
|
P S |
|
S |
|
S |
|
|
|
|
|||
p 1 2 |
c |
|
p |
a1 |
|
a 2 |
|
c |
. |
|
||||||||
P |
|
2 P |
|
|
||||||||||||||
max0 |
|
|
|
max0 |
|
|
|
|
|
|
|
|
||||||
|
|
пор |
|
|
|
|
|
|
|
|
|
|
пор |
|
|
|
|
|
|
|
|
|
|
|
|
Pр |
|
||||||||||
Отсюда следует, что R |
|
~ 4 |
Pр |
|
. Величина |
|
получи- |
|||||||||||
|
|
max0 |
|
|
Рпор |
|
|
|
|
|
|
|
|
Рпор |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ла название - энергетический потенциал РЛС.
Мощность некоррелированного шума Рш связана с пороговой мощностью на входе приемника как Рпор = kPш, где k – коэффициент различимости.
Мощность шума при отсутствии внешних помех
Рш Nш fпр kT0kш fпр .
Коэффициент kш шума зависит от технологии изготовления приёмника и находится в пределах kш = 3…5. Полоса пропускания приемного устройства обратно пропорциональна длительности импульса и . Мощность шума на входе приемника обратно пропорциональна длительности импульса Pпр.min 1/ и .
Если необходимо увеличить дальность действия РЛС Rmax0 в два раза, то для этого требуется увеличить энергетический потенциал РЛС в шестнадцать раз или мощность излучаемого сигнала увеличить на 12 дБ при сохранении чувствительности приемника.
Таким образом, в свободном пространстве дальность действия РЛС зависит от энергетического потенциала РЛС с коэффициентом пропорциональности - корень в степени ¼.
4.2. Оптимальные алгоритмы и методы обнаружения сигналов на фоне пассивных помех
Особенности пассивных помех, отличающие их от сигнала и шума, состоят в следующем [3].
1.Пассивные помехи проявляются только при работе радиопередающих устройств, при этом мощность пассивных помех пропорциональна мощности передатчика.
2.Источник пассивных помех пространственно расположен в непосредственной близости от обнаруживаемого объекта.
39
3. Мощность пассивных помех, как правило, значительно больше мощности шума на входе приемника РЛС. Рш
Рпр ,
т.е. 1.
4.Пассивную помеху от белого гауссовского шума отличает узкополосность спектра.
5.Сигналы, отражаемые от цели, отличаются от пассивных помех меньшей мощностью. Также сигналы, отражаемые от цели, имеют временные и спектральные отличия от пассивных помех.
6.Одно из основных отличий пассивных помех от сигналов целей состоит в том, что отражения от пассивных помех имеют доплеровское смещение частоты, значительно меньшее, чем смещение доплеровской частоты сигнала цели, вызванное значительной скоростью перемещения целей.
fд 2CVr f0 .
В радиолокации значение fд составляет десятки или сотни Герц и существенно зависит от несущей частоты сигнала и радиальной
скорости объекта локации Vr . Доплеровское смещение фазы принятого сигнала д = 2 fдТп. Для однозначного определения дально-
сти период повторения импульсного сигнала T |
2Rmax |
и, напри- |
||||||
|
||||||||
|
|
|
|
|
|
п |
C |
|
|
|
|
|
|
|
|
|
|
мер, |
при дальности 300 км T |
|
600 |
2 мс . Доплеровская фаза в |
||||
|
|
|
П |
|
3 105 |
|
|
|
этом |
случае |
д |
2 3,14 104 |
2 103 40 , т.е. наблюдается ее |
||||
|
|
|
|
|
|
|
|
|
существенная неоднозначность.
Так как поперечный размер элемента разрешения зависит от дальности до РЛС, то радиолокационный контраст снижается с увеличением дальности до цели.
Рассмотрим характеристики пассивных помех. Наибольшее распространение получили полуволновые дипольные противорадиолокационные отражатели (металлизированные бумажные ленты, фольга, синтетическое волокно). Параметрами отражателей являются эффективная поверхность рассеивания (ЭПР), частотная диапазонность, время развертывания, параметры спектра отражений. ЭПР одиночного диполя определяется как
40
SД 0,856 2 cos ,
где θ – угол между осью диполя и направлением излучения. Так как ориентация диполей в пространстве произвольна, то среднее
значение ЭПР одного диполя равно SД 0,172 2 [8]. Для облака из K диполей с коэффициентом разлета ς<1 ЭПР имеет вид: SД 0,172 2 K . Поперечный размер d диполя влияет на коэф-
фициент формы A. Для цилиндрических диполей длиной L коэффициент формы А=L/d. Для прямоугольных диполей А=4L/d. Максимального значения ЭПР достигает, если резонансная частота диполя совпадает с несущей частотой РЛС или, по крайней мере, кратна ей (рис. 4.1).
S/Sp
f0 |
2f0 |
3f0 |
4f0 |
5f0 |
6f0 |
f |
Рис. 4.1. Резонансная характеристика дипольного отражателя
При увеличении коэффициента формы зависимость более резкая. Обычно длина диполя выбирается несколько меньшей, чем λ/2. Для расширения диапазонноти используют диполи разной длины. Если концентрация диполей в единице объема V равна n ,
