УДК 629.78
А. В. Ксендзук, В. Ф. Фатеев, С. А. Попов
МЕТОД ОБРАБОТКИ СИГНАЛОВ В КОСМИЧЕСКИХ МНОГОПОЗИЦИОННЫХ РАДИОЛОКАЦИОННЫХ СИСТЕМАХ С СИНТЕЗОМ АПЕРТУРЫ АНТЕННЫ
Представлен метод совместной обработки сигналов в многопозиционной радиолокационной системе с синтезированием апертуры антенны, основанной на использовании малых космических аппаратов. Приведены основные соотношения и результаты моделирования обработки в многопозиционных системах.
Ключевые слова: дистанционное зондирование, многопозиционная радиолокационная синтезированная антенна, оптимальная обработка, малые космические аппараты, радиолокационное изображение.
Использование малых космических аппаратов (МКА) позволяет перейти на качественно новый уровень решения задач дистанционного зондирования за счет развертывания космической многопозиционной системы мониторинга земной поверхности [1—4]. Такая система должна состоять из приемников и передатчиков, способных осуществлять совместную обработку отраженных от объектов электромагнитных полей (рис. 1, здесь ПРД — передатчик, ПРМ — приемник). Технологическая повторяемость обусловливает относительно невысокую стоимость создания и развертывания такой многопозиционной (МП) радиолокационной синтезированной антенны (РСА).
ПРД
КА-7-
Рис. 1
Одной из основных задач при создании этой системы является разработка методов совместной обработки наблюдаемых полей, что позволит достичь наилучшего качества получаемых данных.
Методы формирования радиолокационного изображения в системах дистанционных исследований могут основываться на двух основных моделях отраженных от поверхности сигналов — стохастической и функционально-детерминированной. Эти модели базируются на представлении комплексного коэффициента отражения от объектов Ё(г), входящего в уравнение отраженного (рассеянного) сигнала в виде функционально-детерминированного или стохастического пространственного процесса [4].
Получим соотношения для оптимальной совместной обработки сигнала в многопозиционной системе с синтезированием апертуры антенны для аддитивной модели уравнения наблюдения, где первый статистический момент ) зависит от комплексного коэффициента
отражения Ё(^, г) , который можно представить в рамках функционально-детерминированных моделей. Оптимальные оценки вектора искомых функций X (электрофизических параметров поверхности как функций координат г и времени I) определяются из решения вариационного уравнения
тт 5 т ( X)
Я ти5хЬ V (*1, ¿2 ) [) - Шв (¿2, X)] ] +
0 0 TT
0 0
dti dt2 = 0,
(1)
где 5 / 5(-) — вариационная производная, К-1^, ¿2) — обратно-корреляционная функция наблюдаемого процесса ) со средним т .
Когда вектор аддитивных помех п(^) является нестационарным процессом со средним т(^) = ), он представляет собой равенство результата обработки входного сигнала в соответствии с соотношением (1)
TT
ff
0 0
R-1(ti, t2 )u(t2 )+UT (ti)R-1(ti, t2)
dt1 dt
1Ш2
где S D (t, X,) — вектор сигналов, отраженных от всей области поверхности
TT ff
0 0
^-р r-1 (ti, t2 )s d (t2, x)+sd (ti, x)r- (t,, t2)
и смещения
TT
b(r) = ff
0 0
R-' (*, t2 )ß(t2 )+ßT (/. )R-1 (/., t2 )
dtidt2;
dt^2.
(2)
(3)
(4)
Для большинства практически важных случаев справедливо предположение о стационарности п(^), по крайней мере, на интервале синтеза апертуры антенны. При этом опти-
мальный алгоритм определяется из равенства, левая часть которого имеет вид
TT
Re ff
0 0
R -'(ti, t2 )u(t2 )+UT (ti)R-i(ti, t2)
dtidt2.
(5)
Правая часть при условии узкополосности спектра сигнала многопозиционной системы, а именно при выполнении для произвольных комбинаций несущих частот условия
maxk,m I f0k - f0m I +2AF < mink,m I f0k + f0m I, (6)
может быть записана так:
1 2
(Т Т ' т т т
^ }} 58^^^ (^хща{2 ¿Г. (^х)*"^,,2)^ . (7)
V 0 0 0 0
Для определения процедуры комплексирования используем решение, учитывающее мнимые части выражений (5) и (7):
тт
И
0 0
(тт 2т
^^ к-1 (,1, ,2 )и(,2)+ит (,1 )Н-1 (,1, ,2)
5Х
5Х
&2 =
Я^т^к-1('1, ,2)8Ь ('2, ^2 £ ('1, Щ-1^,,2)
Л
(8)
V 0 0 0 0 5 )
В случае симметрии функции Н.-1^, '2) относительно транспонирования и перестановки координат ее элементов выражение (8) упрощается до
ттит ('1)Н-1('1, '2) ^¿^(г^'^ = 1 № Ь ('1, ^)Н-1('1, '2) ^'2, (9) 0 0 2 0 0 вырожденная форма которого является классическим уравнением для оптимальных оценок в РСА [6]:
тт ^ ('2,= 1 тт
2'
Я^И-Ч'ь '2) Щ^чан = 2 Я^Ь ('1)^-1('1, '2)
(10)
00 5Х — 5Х
Дальнейшее уточнение процедуры совместной обработки в МП РСА необходимо проводить с учетом (10) на основании моделей сигнала и электродинамической модели поверхности [4].
В качестве модели сигнала будем использовать вектор, координаты которого можно представить в форме
% (') = Яе} Ё [', г; X (' ,г )]8 (', г ) Ог,
Ь
где ядро 8(•) определяется параметрами излучаемых сигналов и геометрией наблюдения. Пусть модель поверхности представима в виде ряда
Ё = £[Х(',г)] £ П/[У(',г),©(',г)], (11)
/=1,..., ь
где — функция комплексирования; п — базисные функции; X (', г) — вектор параметров комплексирования; у(', г) — вектор несущественных параметров, 0(', г) — вектор электрофизических параметров поверхности.
При этом уравнение для оптимальных оценок (10) примет вид:
тт 1 тт
} } ит ('1 ('1, '2)8х ('2, У, г)О'1О'2 = - Я [Х('1, г)]х
0 0 0 0 Ь
х8Г ('1, у, г)к-1 ('1, '2 )8х ('2, У, г )О'1О'2Ог1,
(12)
где
4 г (', У, г) = £ Е П/ [', У, г, 0* (', г)]4 (', г) к=1,..., тг/=1,..., ь
— координаты модифицированного (с учетом модели поверхности) вектора опорных сигналов.
Для статической во временной области модели £(', г) выражение (12) упростится:
ТТ 1 тт
| {ит Ж—1 (¿1, /2)81 (¿2, У, Гц = -1 £[Х(г)] { 18Т (¿1, у, г)х
0 0 2,0 0 0
хя—1 , ¿2 )81 (¿2, у, Г! )^2¿Г. (13)
Правая часть уравнения (13) содержит модифицированную пространственную функцию неопределенности ^ (ПФН) для вектора функции с учетом модели (11)
тт
(г, Г1) = {{ 81 (¿1, у, г)Я-1 (¿1, ¿2 )81 (¿2, У, Г1) , (14)
0 0
которая может быть представлена в виде суммы ПФН по всем комбинациям приемников:
тт
(Г, Г1) = I ^ (Г, Г1) = I {К (¿1, Ъ гЩЧь, ^ (¿2, У, Г^
г,] =1,..., г,у=1,..., 0 0
Левая часть равенства (13) содержит в себе операции декорреляции входного векторно-
т —1
го процесса и (^) матрицей (¿1, ¿2) и последующую согласованную фильтрацию с моди-
• т *
фицированным опорным сигналом 81 (¿2,7, Г1) .
Особенностью разработанных алгоритмов является оценка не вектора X (¿, г), а связанных с ним функции комплексирования ^[Х (^, г)] . Задача определения г)] является обратной и в общем случае некорректной [6]. Однако в соответствии с общепринятой практикой [7] в качестве оптимальных оценок параметров можно принять сглаженный функцией неопределенности результат совместной обработки входного процесса.
При известных априорных значениях плотности распределения вероятности векторов оцениваемых параметров оптимальные оценки определим в рамках метода максимальной апостериорной плотности вероятности путем решения уравнения [7]
5{ р(Х)^Х 5{ I [Ж(и), X] р[и / Х(г)]^и
л
5Х () ■{1 [[ (и), Х]р(и / Х)^и = —{ р(Х)^Х^-—--. (15)
5Х(г) и Л 5Х(г)
В случае аппроксимации плотности вероятности оцениваемого параметра гауссовым распределением оптимальные оценки находятся из уравнения
5 [— 1п(С)+2 { {[Х(Г1)—тх (Г1 )]т Я—1 (Г1, Г2)х
(16)
51п{ р[и / X (г)]}
5Х (r) 2 DD
x[X(r2) - m Х (r2)]^ri^r2]=-
oX(r)
где Rx (ri, r2 ) — второй статистический момент распределения X(r) .
В качестве примера на рис. 2 приведены оценки удельной эффективной поверхности рассеяния а (7), сформированные в однопозиционной РСА (2) и МП РСА, состоящей из трех (3) и четырех (4) приемопередатчиков, а также соответствующие отношения сигнал/помеха (д). Представленные результаты получены в ходе имитационного моделирования на основе оптимальной обработки сигналов при использовании принципов, описанных в работах [5, 8].
Моделирование выполнялось при использовании данных экспериментальных исследований характеристик отраженных сигналов навигационных систем ГЛОНАСС и GPS, проведенных в МАК „Вымпел". В качестве исходных данных использовались параметры группировки навигационных спутников и экспериментально измеренный уровень отраженного от объектов навигационного сигнала относительно уровня прямого сигнала. На рис. 3 приведены
результаты определения относительном мощности по двум навигационным спутникам 13 ГЛОНАСС (7) и 26 GPS (2). Справа — созвездие наблюдаемых навигационных спутников.
Исследование оптимальных алгоритмов в зависимости от пространственной конфигурации, характера перемещения элементов МП РСА, несущих частот, вида используемых сигналов позволяет получить соответствующие значения предельных погрешностей измерений и пространственных функций неопределенности и таким образом обеспечить оптимальное планирование эксперимента дистанционного зондирования с целью получения наилучшего качества результатов.
ст, дБ
20
10
)\ i i! ч. Л1 Nil * a и
3. /ill Hi J i! 1 1 ___ 4 \ Ц4= 18,6
А í/1 1. plWV^ > í 11II а и цэ=7,4 ! : A,'! b. ¡V M h n ÍÍ I t f 2 _i2=1 [ ■ 1
1 ГГ|р 1 Mjt Ф jji Í ! n n || П i Ml fll j ал ' )|U ill I \ [ 1Л k iv № m 4\lt I J/j 1 МШV TÚ
f Ш IS N|¡1 í I'll i ^ k ■, j J i I MM ! ! У v fii гЧ * ! \ j Ц/ ;t 4 Í A tm X . Д J i f w к 1 Шш
) v. р V¡ V н? i *'i \! Í 1 И ! 1 p \ Д4 i'l v h 11VI pi i v 1 11 я V г.
600
800
1000
1200 Рис. 2
1400
1600
y, м
Р, дБ
20
15
5
0
-5
0 200 400 600 г, мс
Рис. 3
Полученные алгоритмы оптимальной обработки могут служить основой для создания различных квазиоптимальных алгоритмов, применимых при отсутствии необходимой априорной информации и/или при необходимости обеспечения измерений в масштабе реального времени. Такие алгоритмы будут обеспечивать меньшее отношение сигнал/помеха, чем оптимальный, но тем не менее существенно уменьшать уровень спекл-шума и повышать качество радиолокационных изображений по отношению к моностатическим РСА. Более того, как показывают исследования, полученные алгоритмы более устойчивы к различным помехам, а также к неточности определения исходных параметров модели поверхности.
СПИСОК ЛИТЕРАТУРЫ
1. Фатеев В. Ф., Сахно И. В. Применение навигационных КА GPS/ГЛОНАСС в составе многопозиционных РЛС обзора земной поверхности // Изв. вузов. Приборостроение. 2004. Т. 47, № 3. С. 27—31.
2. Фатеев В. Ф., Сахно И. В. Бортовая итерферометрическая многопозиционная система радиолокационного обзора Земли, использующая сигналы навигационных КА ГЛОНАСС/GPS/GALILEO // Изв. вузов. Приборостроение. 2007. Т. 50, № 6. С. 11—17.
3. Патент № 2278398 РФ. Способ получения радиолокационного изображения земной поверхности при помощи многопозиционной радиолокационной системы с синтезированной апертурой антенны / И. В. Сахно, В. Ф. Фатеев. 05.06.2004.
4. Ксендзук А. В. Энергетический потенциал ГЛОНАСС/GPS для решения задач ДЗ // Тез. докл. III науч.-практич. конф. „Применение спутниковых радионавигационных систем (GNSS) в Украине". 2002. С. 123—126.
5. Ксендзук А. В., Волосюк В. К., Зеленский А. А., и др.Цифровая обработка сигналов и изображений / Под ред. В. Ф. Кравченко. М.: Физматлит, 2007. 544 с.
6. Волосюк В. К. Физические основы дистанционного зондирования природных сред радиотехническими средствами аэрокосмического базирования. Харьков: Харьковский авиац. ин-т, 1997.
7. Фалькович С. Е., Волосюк В. К., Горбуненко О. А. Радиотехнические системы дистанционного зондирования. Харьков: Харьковский авиац. ин-т, 2002.
8. Ksendzuk A. V., Volosyuk V. K., Sologub N. S. Modeling SAR primary and secondary processing algorithms. Estimating quality of the processing techniques // 5th Europ. Conf. on Synthetic Aperture Radar EUSAR 2004. Ulm, Germany. 2004. Vol. 2. Р. 1013—1016.
Сведения об авторах
Александр Владимирович Ксендзук — д-р техн. наук; ОАО МАК „Вымпел", кафедра информационных систем, Москва; E-mail: [email protected] Вячеслав Филиппович Фатеев — д-р техн. наук, профессор; президент ОАО МАК „Вымпел", кафедра
Сергей Алексеевич Попов
информационных систем, Москва; E-mail: [email protected] канд. техн. наук; ОАО МАК „Вымпел", кафедра информационных систем, Москва; E-mail: [email protected]
Рекомендована Ученым советом ВКА им. А. Ф. Можайского
Поступила в редакцию 20.10.08 г.