УДК 621.384.32
Е. Г. Лебедько, Нгуен Ву Тунг
ПОГРЕШНОСТЬ ИЗМЕРЕНИЯ РАДИАЛЬНОЙ ПРОТЯЖЕННОСТИ ПОВЕРХНОСТИ ПРИ ЕЕ МОНОИМПУЛЬСНОМ ОБЛУЧЕНИИ
Представлен анализ основных случайных составляющих погрешности определения радиальной протяженности облучаемой поверхности, к которым относятся: шумовая погрешность, вызванная сдвигом отсчетных точек под действием шумов; погрешность, обусловленная уровнем отсчета сигнала при априорной неопределенности формы облучаемой поверхности и ее отражающих свойств; аномальная погрешность, вызванная принятием шумового выброса за сигнальный.
Ключевые слова: локация, дисперсия оценок, импульсная характеристика, радиальная протяженность.
При исследовании планет и астероидов Солнечной системы используются беспилотные аппараты и автоматические транспортные средства, которые могут спускаться на платформах с борта космических аппаратов. При их посадке, во избежание опрокидывания, возникает необходимость определения структуры рельефа посадочной поверхности. Особенно важно при этом знать перепад высот посадочной поверхности. Оперативный анализ ее рельефа должен осуществляться с достаточно больших высот; погрешности изменения перепада высот составляют десятые доли метра. Определенные перспективы для работы таких навигационных систем открывает трансформация пространственных признаков во временные в отраженном сигнале при моноимпульсной локации исследуемой поверхности. В этом случае параметры отраженного сигнала содержат информацию о пространственной структуре облученной поверхности и ее радиальной протяженности при однократном зондировании. При этом радиальная протяженность облучаемой поверхности характеризуется длительностью импульсной характеристики поверхности на нулевом уровне отсчета.
Под импульсной характеристикой поверхности понимается отраженный сигнал от поверхности, облучаемой импульсом в виде дельта-функции. Импульсная характеристика поверхности определяется зависимостью [1 ]
§э () = Иг (у) с^2 у5 ^*-^ Ж,
где 5 — площадь облучаемой поверхности; у — угол между направлением облучения и нормалью к элементу поверхности ds; г (у) — коэффициент яркости элемента поверхности; 5 () — дельта-функция; с — скорость света.
Такой важный информационный параметр, как длительность импульсной характеристики поверхности изменяется вследствие того, что диаграммы излучения и приема отличаются от идеальных, зондирующий импульс не является дельта-функцией, а в условиях априорной неопределенности появления облучаемой поверхности не удается оптимизировать приемно-усилительный тракт. Перечисленные факторы определяют систематическую погрешность измерений, которую можно в определенной степени скомпенсировать в измерительной схеме приемной системы.
Наряду с этими систематическими погрешностями возникают случайные составляющие погрешностей, к которым относятся:
1) шумовая погрешность, вызванная сдвигом отсчетных точек под действием шумов;
2) погрешность, обусловленная уровнем отсчета сигнала при априорной неопределенности формы облучаемой поверхности и ее отражающих свойств;
3) аномальная погрешность, обусловленная принятием шумового выброса за сигнальный;
4) методическая погрешность схемы измерения временного интервала;
5) погрешность, возникающая под действием дестабилизирующих факторов, к которой, в частности, относится погрешность, вызванная условиями распространения излучения (турбулентностью атмосферы, если она имеется).
Настоящая статья посвящена анализу 1, 2 и 3-й случайных составляющих погрешностей при работе навигационной системы в отсутствие атмосферы.
1. Так как искомый информационный параметр содержится в длительности принимаемого сигнала, то дисперсия шумовой погрешности, наряду с дисперсиями оценки фиксации временного положения сигнала по фронту и спаду, в силу корреляции шума в этих точках определяется зависимостью
О = Офр +°2п -В2 (т3 ) ,
2 2
где Офр и осп — дисперсии оценки фиксации (запаздывания) при пересечении заданного уровня фронтом и спадом сигнала на выходе приемно-усилительного тракта; В2 (т) — корреляционная функция оценок при длительности Т3 сигнала на уровне отсчета (фиксации).
Будем считать, что сигнал на выходе приемно-усилительного тракта симметричен, тогда
°2 = 2офр -В2 (т3 ).
В условиях приема на фоне белого гауссова шума с энергетическим спектром О дисперсия оценки фиксации временного положения сигнала по фронту определяется зависимостью
G J[g (t)]2 dt
2
°фр =
]
J *2(t-T0)g(t0-t)dt
где g (t) — импульсная характеристика приемно-усилительного тракта, который вследствие априорной неопределенности структуры облучаемой поверхности настроен на прием сигнала с параметрами зондирующего импульса; s2 (t - r0j) — производная интегральной свертки
зондирующего импульса и системной импульсной характеристики поверхности, которая определяется зависимостью
g* (t) = JJ Vi (вx - A*, вy - Ay ) у2 (вx - A*, вy - Ay) r (у) cos2 y5 (t-2z/c) ds,
здесь Vi (x, y) и у2 (x, У) — диаграммы полей излучения и приема соответственно.
В силу того что форма оптического импульса заданной энергии незначительно влияет на обнаружение и оценку параметров принимаемых сигналов, для описания диаграмм полей излучения и приема, а также аппроксимации зондирующего импульса будем при расчетах использовать гауссовы функции. В этом случае получим следующие выражения:
пг-тЗ ( Л \
2 2GT20 to
°фр = 2 2 2 ( ) exp ~ п t0а (y) 1т20.
2
B (Тз ) = -
>/2 Gt20 exp ( /2 > ,т ¿0 пТ2 ^ Т20 ) exp ( Т2 ^ Т30 п 2 1 2т20 )
a2 (у)п2¿о2
здесь т20=
2 2 Т2 +т2
(2)
где То и Т3 — длительности зондирующего импульса и системном им-
пульсной характеристики на уровне 0,5 от максимального значения; Т30 = ¿сп -¿фр, где и ¿сп — математические ожидания оценок фиксации сигнала по фронту и спаду соответ-
¿фр
ственно; а (у) — значение принимаемого сигнала.
Определенный интерес представляет зависимость дисперсии оценки и корреляционной
функции от уровня фиксации Пр. В этом случае ¿о =
2ln p
п
Т20, и соотношения (1), (2) при-
нимают следующий вид:
2
афр =
Gt20 Р Р
2 2 2 аш т20 Р
B (Р) =
V2 •пa (у)lnp •пa (y)lnp
2 2
Gt20 °шт20
(3)
(4)
л/2 • па2 (у)р 21п р л/2 • па 2 (у)р 21п р'
где аШ = С/ Т20 — дисперсия шумов на выходе приемно-усилительного тракта.
Таким образом, с учетом выражений (3) и (4) дисперсия шумовой погрешности определяется зависимостью
а2 =
_2 2 аш т20
(
V2 •пa (y)lnp
2 p2 - 4 1 p )
(5)
22
На рис. 1 приведен график зависимостей а , афр и В2 от уровня фиксации сигнала.
Нормировка проводилась при условии, что
_2 2 аш т20
V2 •п а1 (у)
= 1.
2 2 „2 а , афр , В ■
20 — —1 — —
10
1,2
3
Рис. 1
Как видно из приведенных графиков, ве-
22
личины а и афр имеют минимальное значение при 1/р = 1/у[ё. Минимальное значение дисперсии оценки фиксации сигнала будет при максимальном значении производной сигнала (для гауссовой функции этот уровень составляет 1Д/ё ).
В соответствии с формулой (5) и рис. 1 можно сделать заключение о том, что основной вклад в шумовую погрешность вносит дисперсия оценки фиксации по фронту.
2. Истинное значение временной протяженности облучаемой поверхности (Т1)
пропорционально длительности импульсной характеристики, отсчитанной на нулевом уровне. При заданном уровне отсчета значения длительностей импульсных характеристик для
2
о
различных поверхностей при одинаковой их радиальной протяженности различны и зависят от формы поверхности.
Этот вопрос чрезвычайно важен, так как даже в теоретически идеальном случае, когда поверхность облучается равномерным круглым полем с временной функцией в виде 5-функции и отраженный сигнал является импульсной характеристикой поверхности, вследствие воздействия шума уровень отсчета необходимо выбирать, исходя из заданных вероятностных характеристик принятия решения. А так как уровень отсчета не нулевой, и форма облучаемой поверхности неизвестна, то обусловленная этим факторами погрешность является случайной и скомпенсировать ее не удается.
Следует отметить, что, например, при выборе уровня отсчета, равным 0,5 от максимального значения импульсной характеристики, относительная погрешность определения временной протяженности поверхности может достигать 50 % даже для элементарных поверхностей, рассмотренных в работе [2]. Анализ этой погрешности произведем на примере двух поверхностей — наклонной плоскости и конуса, импульсные характеристики которых наиболее известны и длительность их описывается различными функциями от уровня отсчета (что и вызывает погрешность). При этом временная протяженность поверхностей одинакова.
На рис. 2 приведены диаграммы изменения относительной погрешности А определения временной протяженности поверхности в зависимости от уровня отсчета при двух значениях длительности зондирующего импульса: а — т0 =11/2, б — т0 = 2^/3 . На графиках кривая 1 соответствует уровню отсчета, равному 0,5, а кривые 2 и 3 — уровням 0,35 и 0,25 соответственно.
0,4
0,2
0
Рис. 2
Диаграммы построены для трех позиций: первая а — для импульсных характеристик поверхностей, когда их первичные информационные параметры не искажены; вторая Ь — для отраженных сигналов, когда искажение информационных параметров обусловлено длительностью и формой зондирующего импульса, а также полями излучения и приема; третья с — для выходных сигналов, в которых происходит дополнительное искажение первичных информационных параметров в приемно-усилительном тракте.
Как видно из представленных диаграмм, относительная погрешность определения радиальной протяженности облучаемой поверхности сокращается как по мере роста искажений информационных параметров, так и с уменьшением уровня отсчета. Следует отметить, что минимизация этой погрешности требует энергетических затрат и соответствующей технической реализации системы.
3. В условиях, когда отношение порог/шум недостаточно велико, принятие шумового выброса за сигнальной может привести к результатам, отличным от истинных значений оценок длительности принимаемого сигнала. Ошибки, обусловленные принятием шумового
1 1 1 1 1 1 1 1
2 1 1 1 1 1 1 1 1
з 1 1 1 1 к 1 X. 1
1 | I
а Ь с
выброса за сигнальный и выводящие оценку параметра за пределы протяженности сигнальной функции, называются аномальными.
Дисперсия аномальной погрешности определяется зависимостью
2
_2 = р2 __ ° а Г д.тр° т '
где Рл тр — вероятность ложной тревоги, ат — среднее квадратическое значение длительности шумового выброса на пороговом уровне С .
Так как вероятность ложной тревоги в задачах анализа рельефа поверхности может со-
*5... 7
ставлять величину 10 ... 10 , а следовательно, отношение порога к шуму До = С/ аЕ то для С > 0 большинство положительных выбросов имеют малую длительность. В этом случае поведение дифференцируемого случайного процесса при С > 0 аппроксимируют параболой [3, 4], и плотность вероятности длительности выбросов случайного процесса определяется соотношением [5]
W (т, С )=4 (-¿0>2 Т exp
1
2 2
-ттНО )н2 т
с
; т > 0, до =->> 1
(6)
1 »
где -bJ =-— | ®2G(c)dш — мощность производной случайного процесса; аШ = B2 (о);
2па
О (о) — энергетический спектр шумов.
С учетом выражения (6) дисперсия длительности шумового выброса
ат =
|[т-т(С)]2 4(-¿0
Нотexp
~(-bJ )н0 Т2
d т ,
где т(С) — средняя длительность выброса случайного процесса на уровне С :
(C )="^[1"Ф(До)] exp
( 2 ^
До 2
v
1 До
здесь Ф(но )=— f exp 2п J
( ^ ^ v"2 )
dx — интеграл вероятностей.
Произведем анализ аномальной погрешности при условиях оптимальной фильтрации принимаемого сигнала гауссовой формы на фоне белого нормального шума. Тогда
а T=i
8т2О
пн2
л2^2 (erfc(но ))2 exp (н2)
-\/уerfc(но )exp
( 2^ До
2
v )
+1
(7)
где erfc(x) = 1 - erf (x), erf (x) — функция Крампа.
Результаты расчетов аТ по формуле (7) сведены в таблице.
ш
Но аТ, с, при т20=2-10-8 с аТ, с, при т20=5-10-8 с аТ, с, при т20=10 7 с
3 1,06-Ю-8 2,649-1O-8 5,299-Ш-8
4 O,789-1O-8 1,995-Ю-8 3,989-Ю-8
5 O,638-1O-8 1,596-iO-8 3Д92-10-8
7 O,456-1O-8 1,14-1O-8 2,28-1O-8
Как следует из приведенных рассуждений, не всегда можно пренебрегать аномальной погрешностью: необходимо ее учитывать при оценке временной протяженности облучаемой поверхности и искать пути ее снижения посредством повышения порогового уровня и временного стробирования.
Анализ случайных составляющих погрешности измерения радиальной протяженности облучаемой поверхности показал, что наиболее опасной следует считать погрешность, обусловленную уровнем отсчета при априорной неопределенности формы облучаемой поверхности и ее отражающих свойств. При этом даже при равновероятном появлении пространственной структуры облучаемой поверхности плотность вероятностей этой погрешности остается неизвестной и приходится рассматривать максимальное ее значение.
СПИСОК ЛИТЕРАТУРЫ
1. Лебедько Е. Г., Порфирьев Л. Ф., Хайтун Ф. И. Теория и расчет импульсных и цифровых оптико-электронных систем. Л.: Машиностроение, 1984.
2. Лебедько Е. Г., Нгуен Ву Тунг. Определение параметров отраженных сигналов в моноимпульсной лазерной системе анализа рельефа поверхности // Изв. вузов. Приборостроение. 2007. Т. 50, № 5. С. 63—69.
3. Кузнецов П. И., Стратонович Р. Л., Тихонов В. И. О длительности выбросов случайной функции // ЖТФ. 1954. Т. 24, вып. 1
4. Стратонович Р. Л. Избранные вопросы теории флуктуаций в радиотехнике М.: Сов. радио, 1961.
5. Тихонов В. И. Выбросы случайных процессов. М.: Наука, 1970.
Сведения об авторах
Евгений Георгиевич Лебедько — д-р техн. наук, профессор; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра оптико-электронных приборов и систем; E-mail: eleb@rambler.ru
Нгуен Ву Тунг — аспирант; Санкт-Петербургский государственный университет инфор-
мационных технологий, механики и оптики, кафедра оптико-электронных приборов и систем; E-mail: eleb@rambler.ru
Рекомендована кафедрой Поступила в редакцию
оптико-электронных приборов и систем 26.03.08 г.