CC BY
74
УДК 621.396.41
М.Ю. Нестеров, В.В. Мухин, А.П. Макрушин
Двухэтапный алгоритм измерения высоты
Рассмотрен алгоритм определения высоты по максимуму функционала правдоподобия. Показано, что представленный метод обладает меньшей ошибкой по сравнению с традиционным способом измерения.
Ключевые слова: радиовысотометрия, алгоритм.
В традиционном радиовысотомере (РВ) с непрерывным ЛЧМ-сигналом информация о высоте полета извлекается из сигнала биений. В первом приближении и в отсутствие вертикальной скорости частота сигнала биений /ь пропорциональна текущей высоте Н :
fb -
clm
(1)
где - девиация частоты излучаемого сигнала, с - скорость света, Тт - период модуляции.
Считаем, что полная мощность _Ротр сигнала, отраженного от однородной и изотропной шероховатой поверхности, распределена в спектральной области согласно функции P(/):
Ротр -| P( f )df ,
(2)
где
Р(/) = |8[/ - /Ь (0,ф)]ёР. (3)
Полагая, что передающая и приемная антенны высотомера совмещены, ёР - мощность сигнала, отраженного от элемента поверхности ёБ, равна
dP -
G \ 2 cos4 0 64n3H 4
ga(0^)c(0)ds,
(4)
где G - коэффициент усиления антенны, ga (0,ф) - нормированная диаграмма направленности антенны (ДНА), ст(0) - удельная ЭПР поверхности.
На рис. 1 представлена СПМ сигнала биений, полученная в натурном эксперименте в неследящем режиме при следующих условиях: высота
H «1150 м, девиация частоты излучаемого сигнала AF = 120 МГц, период модуляции Tm =1 мс, скорость полета V ~ 70 ^с, подстилающая поверхность - ровная, с покрытием типа короткой травянистой растительности, ширина диаграммы направленности антенны порядка 20°, длина волны 2,2 см. На том же рисунке изображена кривая СПМ, рассчитанная для этих же условий по описанному выше алгоритму. Диаграмма обратного рассеяния аппроксимировалась согласно [3], приложение «С»
(GRASSES - поверхность типа «луг»).
1100
1150
1200
1250
H, м
Рис. 1. Спектр сигнала биений
Для удобства частота сигнала биений приведена к соответствующей высоте. Как можно убедиться, наблюдается хорошее соответствие расчета и экспериментальных данных.
Как следует из анализа приведенного выше расчета и экспериментальных данных, огибающая спектра имеет несимметричный вид и растянута вправо из-за протяженности отражающей поверхности. Чем шире диаграмма обратного рассеяния или диаграмма антенны, тем больше кривая растягивается в область больших высот. При этом левый склон изменяется мало, стремясь к уровню 0,5 от максимума на отметке истинной высоты (см. также [2]).
Разрешающая способность РВ по высоте не может быть лучше величины
c
SH = ^—, (5)
2ДЕ
что для приведенных выше характеристик РВ составит 1,25 м. Существенно улучшить эту величину за счет увеличения полосы сигнала не представляется возможным по причине технической сложности и требований регламента радиосвязи. Однако к перспективным РВ предъявляются более жесткие требования по точности, особенно в области малых высот.
Таким образом, задача состоит в разработке метода спектрального оценивания, определяющего положение переднего фронта спектра с точностью лучше дискрета БПФ.
Следует иметь в виду, что форма спектра сигнала биений РВ на малых высотах и при относительно небольшой скорости полета достаточно характерна. Имеется резкий передний фронт и отно-
сительно затянутый, в зависимости от ширины ДНА и ДОР поверхности, задний скат. Опишем эту форму следующей аналитической моделью:
N -1
5(,х) = | А(т,х)Е(т,^т, (6)
0
где N - длина БПФ;
{-а(т-х) .
е ^ х, т> х (7)
0, т< х
- аналитическое описание затухания спектра от дальности (х - истинная высота), а
пі t-т
~21~5Г
(8)
- модель строба дальности на выходе системы обработки сигнала (после сжатия).
Неизвестными параметрами здесь являются задержка х, вообще говоря, дробная (в единицах дискрета дальности), а также показатель затухания заднего фронта спектра а . Оценку величины х предлагается проводить в 2 этапа:
1) с точностью до дискрета по дальности (высоте) определить положение переднего и заднего фронта, а также параметра а ;
2) произвести уточнение задержки х в окрестности найденного на предыдущем этапе положения переднего фронта.
Рассмотрим этот алгоритм более подробно. Предположим, что для комплексной амплитуды отдельно взятого отсчета спектра выполняется гауссовская статистика. Также предположим, что отсчеты спектра статистически независимы. На первом этапе эти предположения оправданы, поскольку параметры спектра необходимо определить лишь с той степенью точности, которая нужна для следующего шага. Исходя из этого, можно выписать логарифм функции правдоподобия для спектра, передний фронт которого соответствует отсчету под номером П1, а задний фронт - отсчету П2 :
^1 = -(( - П2 + П1)1пстП (,П2) - ( - П1)1п СТ2 (,П2) -Х1п№г , (9)
где
1
N - П2 + щ
Щ -1 2 2
IN2 +1N2
і =1 г=«2
(10)
2 1 п2-1 |у-|2
а2 =—— I ЕІ-. (11)
п2 - п1 і=щ ™і
Весовой вектор w - exp |^-а(п - пі)| n -nl...n2 -lj определяется параметром спада заднего фронта спектра
2
1 l y nl (12
а------------l^-LT-!—Г. (12)
П2 - пі ап (пі,П2 )
Следует отметить, что если спектр сигнала биений одномодален, а это действительно имеет место в отсутствие преднамеренных помех и многократных переотражений в канале «поверхность -
корпус ЛА», то нет необходимости перебирать все возможные значения n1 и n2 . Можно двигаться в
поиске от максимума спектра в сторону предполагаемого переднего фронта, при этом не рассматривая варианты спектра с эффективной шириной более максимально возможной величины для текущих условий полета. Это существенно сократит время поиска.
2 2
Теперь, определив параметры а , as и ап, и грубо оценив положение переднего фронта пі , можно записать выражение для элемента корреляционной матрицы (для одной квадратуры) при неизвестном срезе спектра x :
N-1
K(tj ,ti,x) = a2 J A2 (т,x)F(ti,t)F(tj,т)т + аП 5( ,tj). (13)
0
Отметим, что при гауссовском приближении выхода согласованной обработки сигнала это выражение сводится к табличным функциям, что позволяет вычислить его небольшими затратами (окончательная формула не приводится по причине громоздкости).
Логарифм функции правдоподобия для вектора отсчетов спектра в окрестности n1 и неизвестной задержки n1 :
F2 —ln (2tc|K| )-2 yC K-1y c-2 yS K-1y 5, (14)
где yc и ys - векторы квадратур отсчетов спектра.
Значение x, при котором достигается максимум функции правдоподобия, будет искомой оценкой текущей высоты.
Описанный алгоритм был испытан на данных, записанных в ходе натурных испытаний. Результаты испытаний подтвердили работоспособность алгоритма. Флуктуационная ошибка измерения высоты над ровной поверхностью составила ±0,5 м по уровню За (в области малых высот).
Литература
1. Rice S.O. Mathematical analysis of random noise // BSTJ. - 1945. - Vol. 24, № 1. - P. 46-156.
2. Жуковский А.П. Теоретические основы радиовысотометрии / А.П. Жуковский, Е.И. Оноприенко, В.И. Чижов. - М.: Советское радио, 1979. - 320 с.
3. Ulaby F. Handbook of Radar Scattering Statistics for Terrain / F. Ulaby, M. Dobson. - Artech House, 1989. - 357 с.
Нестеров Михаил Юрьевич
Канд. техн. наук, начальник отдела ОАО «УПКБ «Деталь», г. Каменск-Уральский Свердловской обл. Тел.: 8 (343-9) 33-92-98 Эл. почта: [email protected]
Мухин Владимир Витальевич
Начальник НИО ОАО «УПКБ «Деталь»
Тел.: 8 (343-9) 33-92-80 Эл. почта: [email protected]
Макрушин Андрей Петрович
Инженер 2-й категории ОАО «УПКБ «Деталь» Тел.: 8 (343-9) 33-92-98 Эл. почта: [email protected]
Nesterov M.Yu., Mukhin V.V., Makrushin A.P.
Flight altitude estimation via calculation of spectral moments of the beat signal
The paper presents a new method of FMCW radar altimetry. Current altitude is estimated via calculation of two-stage algorithm. The paper also provides experimental results that prove better accuracy of the method in comparison with the traditional approach.
Keywords: radar altimetry, algoritm.
