CC BY
93
МОДЕЛИРОВАНИЕ СИСТЕМ И ПРОЦЕССОВ X
УДК 621.391.826; 681.5
МАТЕМАТИЧЕСКИЕ МОДЕЛИ ЭХО-СИГНАЛОВ КОРАБЛЕЙ, НАБЛЮДАЕМЫХ ЛОКАТОРАМИ БОРТОВЫХ СИСТЕМ ОБРАБОТКИ ИНФОРМАЦИИ
А. Г. Давидчук,
ассистент Д. А. Шепета,
канд. техн. наук, доцент
Санкт-Петербургский государственный университет аэрокосмического приборостроения
Предлагается математическая модель эхо-сигналов кораблей, построенная на основе экспериментальных данных, которая позволяет учитывать флюктуации амплитуд и длительностей сигналов, а также их взаимные корреляционно-спектральные характеристики. Математическая модель основана на логарифмически нормальном многомерном законе распределения флуктуаций амплитуд и длительностей.
The mathematical model of ships’ echo - signals, constructed on the experimental data basis which allows to take into account fluctuations of amplitudes, durations of signals, and mutual correlation and spectral characteristics is offered. The mathematical model is based on the logarithmic-normal multidimensional law of distribution fluctuations of amplitudes and duration.
При построении математических моделей входных сигналов бортовых систем обработки информации применяются в основном три подхода: первый состоит в постулировании математической модели; второй основан на изучении физических явлений, обуславливающих появление сигналов; третий использует «статистический эквивалент» сигналов, построенный по экспериментальным данным. Под математической моделью в данной статье понимается статистическая модель - многомерный совместный закон распределения параметров сигналов, наблюдаемых на выходе приемного устройства бортового локатора. За основу принят третий подход, который при построении модели учитывает выводы относительно статистических характеристик наблюдаемых сигналов (приведенные в литературе по обработке экспериментальных данных).
Математическая модель информационного сигнала на входе локационного тракта
При поиске, обнаружении, идентификации, выборе, захвате на автосопровождение и сопро-
вождение корабля математическая модель информационного сигнала зависит, по крайней мере, от пяти основных факторов: характеристик самого корабля, тракта распространения электромагнитных волн, характеристик зондирующего сигнала, условий наблюдения корабля, приемного тракта бортового локатора [1, 4, 9-12]. Определяющими факторами являются характеристики корабля и вид излучаемого сигнала [9,12]. Остальные факторы можно учесть следующим образом: характеристики тракта - изменением параметров модели, условия наблюдения - введением функциональных зависимостей между условиями наблюдения и параметрами модели. Тракт обработки считаем линейным и широкополосным, что позволяет рассматривать независимо модели информационного сигнала и помех.
Под информационным сигналом, в соответствии с принятыми допущениями, будем понимать отраженный физическим объектом зондирующий сигнал, прошедший тракт распространения и приемное устройство бортового локатора. После гетеродина приемного устройства наблюдаемый ин-
формационный сигнал 5(£) с точностью до коэффициента пропорциональности можно записать в виде
N М
5(0 = л А* - *о.1 “1)) х
1=1 ]=1
хсоз(ю0£ + фси + Ф0,у), (1)
где - задержка у'-й пачки сигнала; 5(£) - норми-
рованная огибающая импульсов; со0 - промежуточная частота последнего каскада приемного устройства; М - число импульсов в наблюдаемой пачке; N - число принимаемых пачек импульсов.
Здесь принято, что за время отражения пачки условия отражения не меняются, что всегда выполняется в реальных условиях наблюдения морских целей (при импульсно-пачечном режиме). Выражение (1) представляет собой достаточно общую форму записи сигнала, позволяющую учесть фазовую, частотную, амплитудную и временную (период следования пачек импульсов) модуляцию (манипуляцию) сигнала. Частными случаями этого выражения являются выражения для одиночного импульсного сигнала, ЛЧМ-сигнала (при соответствующем предельном переходе), ФМ-сигна-ла, сигнала со случайной или детерминированной вобуляцией периода повторения.
После амплитудно-фазовой обработки сигнала у-я пачка полностью определяется характеристиками лишь одного своего импульса, а именно: его амплитудой = /г0;Аз, начальной фазой сру =
- ф0;-, задержкой пачки t0j. Задержку пачки t0j будем считать детерминированным параметром сигнала, что допустимо, если пренебречь дально-мерным шумом. Тогда каждая пачка характеризуется случайными величинами Ау, фу и огибающей £,(£).
Таким образом, при введенных допущениях информационный сигнал можно представить в виде последовательности импульсных сигналов; £-й импульс этого сигнала запишем в виде £>*(£) = = £Д£) соз(со01 + ф£), где 5;(г) - огибающая импульса; со0 - частота (может быть равной нулю); ф£ - начальная фаза (измеренная относительно сигнала, начальная фаза которого условно принята за нуль).
Данная последовательность импульсов (пачка импульсов) может быть когерентной (начальные фазы импульсов связаны друг с другом известным соотношением) и некогерентной. Рассмотрим случай некогерентной пачки импульсов, так как при наших допущениях это более общий случай (в том смысле, что для когерентной пачки достаточно задать начальную фазу лишь одного из импульсов), при этом частоту можно считать известной детерминированной величиной. Для некогерентной последовательности начальные фазы фг, £ = 1, 2, ..., ./V, представляют собой совместно независимые случайные величины, распределенные равномер-
но на интервале [0, 2л). Огибающая 5Д£) представляет собой известную функцию, зависящую от двух параметров: амплитуды сигнала^ = тах£^(£) и его длительности т£, определенной по уровню 0,5, т. е. ть = \Ь{' - ti'\, где и £/ - два крайних члена из упорядоченного ряда решения уравнения 5£(£) = = 0,5А*.
В качестве 5; обычно используют следующие функции:
- прямоугольную
в#)
0, £<-0,5тг, А, |^| — °»5т£, 0, t > 0,5т£;
(2)
треугольную
0, т- <|*|,
А +Т/~1А*> -т; ^<0, А-тГЧ*. 0<t<V,
(3)
гауссову
St(t) = Ai ехр{-41п2тг 2t2}
(4)
Если длительность зондирующих импульсов сравнима с радиолокационным размером корабля, то форма эхо-сигналов близка к гауссовской кривой [2, 5, 9].
Таким образом, наблюдаемый импульс характеризуется четырьмя параметрами: амплитудой Aj, длительностью т£, начальной фазой ф£, частотой со0. Три первых параметра случайные, а четвертый — ш0 = const. Следовательно, полезный сигнал полностью определяется случайной трехмерной величиной jj = {Ai? тг, фJ и параметром ю0. Поэтому для формирования математической модели сигнала необходимо задать параметр со0 и многомерный закон (функцию или плотность распределения) случайных величин ух, у2, ..., yN.
Для рассматриваемых некогерентных сигналов случайные многомерные величины {Ах, ij, ...,AN, Тдг} и {фх, ..., фдг} взаимно независимы, а фг - равномерно распределены на интервале [0, 2л). При любых i и N многомерная плотность распределения равна
N
/(Yi, = Tj,An, Тд.^Яфу), (5)
;=i
где
0, ф^ < 0,
(2л)-1, 0 < ф£ < 2л,
0, фу > 2л.
(6)
Для когерентной пачки выражение (6) будет содержать произведение из N - 1 дельта-функций.
Из выражения (6) следует, что для получения математической модели сигнала достаточно задать многомерную плотность распределения /(Ах, хх, А2, х2, ..., Лдг, Хдг), так как статистические характеристики флюктуаций фаз определены выражением (6). Следовательно, для математического описания полезного сигнала необходимо определить плотность распределения /(А1? т 1,А2, т2,
А\г» тту) ~ ((А-ы, х дг), N = (1, 2, ...), где Адг = {А1? А2, ...,А^}, тм= {хх, т2, тдг} (Ая и х^- векторы амплитуд и длительностей сигналов соответственно).
Плотность распределения амплитуд и длительностей локационных сигналов, отраженных от надводных объектов
В качестве одномерных функций распределения амплитуд использовались разные функции распределения: Рэлея, Рэлея-Райса, Накагами, логарифмически нормальная, хи-квадрат и др. [3, 6-8, 11, 12]. Наиболее часто в последнее время используется логарифмически нормальная плотность распределения амплитуд, так как она очень хорошо согласуется с экспериментальными данными, полученными для различных типов кораблей и различных условий их наблюдения. Эту плотность распределения и будем использовать для построения многомерной модели эхо-сигнала корабля.
Одномерная функция плотности распределения амплитуды А;- локационного импульса, распределенного по логарифмически нормальному закону, определяется по формуле
ДА) =
ехр
2 А
1п
л
2а
А
(7)
где А1 и а,- - параметры распределения, причем А1 и связаны со средним А1 и дисперсией Г)А распределения следующими соотношениями:
А; = АехР'
ч
^А.-Л2
(8)
Однако на практике вместо Г)А для характеристик флюктуаций амплитуд локационного сигна-
ла используют коэффициент вариации К А = В этом случае из выражения (8) получим
і
Б
А
А
А
А
1 + К‘
+ КІ).
(9)
Иногда вместо коэффициента вариации по амплитуде КА используют коэффициент вариации по мощности Кр . Плотность распределения мощности fp также логарифмически нормальна с параметрами:
В =
В
і+кі
= 4г’ =А/іп(1 + Я'1і)=2аА, (10)
где р. - среднее значение мощности; 1 + Кп =
2 4 ~
= (1 + КА ) ; значение Р1 определяется с учетом параметров локационной структуры по известной формуле
В=-
(4л) Д'
(И)
где ^пер, - мощность передатчика при излучении 1-го импульса; - коэффициент усиления антен-
ны по мощности (при излучении и приеме на одну антенну); - длина волны передатчика; Т1£ - коэффициент потерь, г\1 «0,25; - средняя эффек-
тивная поверхность рассеивания (ЭПР) при £-м зондировании; - расстояние между антенной и объектом.
Выражения (9)-(11) полностью определяют параметры одномерной плотности распределения амплитуд локационных сигналов.
Для аппроксимации одномерной плотности распределения длительностей импульсов эхо-сигналов кораблей предлагались различные законы распределения: хи-квадрат, нормальный, логарифмически нормальный, усеченный нормальный. Наиболее простыми являются нормальный и логарифмически нормальный. Следует отметить, что использование нормального закона некорректно в силу ненулевой вероятности появления отрицательных значений длительности. Переход от нормального закона к усеченному нормальному приводит к корректным выражениям, но настолько усложняет модель, что делает ее практически бесполезной. Все перечисленные законы (с учетом оговорки для нормального закона) достаточно хорошо согласуются с экспериментальными характеристиками, и поэтому при выборе закона (из перечисленных) следует исходить из соображений простоты модели, ее удобства для аналитических вычислений и синтеза алгоритмов моделирования.
В силу сказанного выше имеет смысл принять статистические характеристики флюктуаций хг такие же, как и для амплитуд. Физически флюктуации длительности отраженных импульсов определяются теми же причинами, что и флюктуации амплитуд.
Итак, примем, что длительность отраженного импульса распределена по логарифмически нормальному закону. Тогда функция плотности распределения определяется следующим выражением:
/Ы =
л/2тс<
-ехр
7ГСГ X;
1п2^
2ст!
(12)
где параметры распределения так же, как и в случае с амплитудой, равны
(13)
(и - средняя длительность и коэффициент вариации длительности г-го импульса соответственно).
Таким образом, одномерные плотности распределения амплитуд и длительностей эхо-сигналов кораблей полностью определены. Эти плотности распределения хорошо согласуются с экспериментальными данными, имеют достаточно простые аналитические выражения, которые удобно использовать в аналитических выкладках и при синтезе практических алгоритмов моделирования флюктуаций амплитуд и длительностей эхо-сигналов кораблей.
Двумерные плотности распределения и корреляционные функции амплитуд и длительностей
Знания одномерных функций распределения амплитуды и длительности отраженных сигналов еще недостаточно для записи двумерной функции плотности распределения этих величин. Для определения двумерных плотностей распределения амплитуд, длительностей, амплитуд и длительностей необходим очень большой объем выборочных данных, к тому же полученных при одних и тех же условиях, что практически невозможно. Поэтому на практике оценивают лишь одномерные законы распределения амплитуд и длительностей, а также ковариационные функции амплитуд, длительностей и взаимную корреляционную функцию амплитуд и длительностей. Затем подбирают такие двумерные плотности распределения, для которых маргинальные плотности распределения и ковариационные зависимости не противоречат экспериментальным данным. Эти двумерные плотности распределения и постулируются в качестве математических моделей, описывающих совместные флюктуации амплитуд и длительностей. Используем эту методику для получения соответствующих двумерных плотностей распределения.
На двумерные плотности распределения накладываем следующие требования: маргинальные плотности должны совпадать с оценками одномер-
ных плотностей для амплитуд, длительностей, амплитуд и длительностей; коэффициенты корреляции амплитуд, длительностей, амплитуд и длительностей должны совпадать с соответствующими оценками коэффициентов корреляции. Этим требованиям удовлетворяют двумерные логарифмически нормальные плотности распределения.
Для амплитуд двумерная логарифмически нормальная плотность распределения записывается в виде
Я 4,А;.) =
2коаоаА1А]^.
■ехр •<
1 — г
А,
2(!-<)
1п‘
А
6
А,
2 г
1п^- 1п
А
А;
Ч
А;
где -параметр распределения.
Коэффициент корреляции амплитуд IIА определяется следующим образом:
оооо
-Л'ХА -Я
‘
О О
А А
V
ВА ВА1
Ч/-1) _
А/ / л; /1;
е
какА)
(15)
Из этого выражения можно найти гА в виде 1п(1 + КаКаЕ^)
^1п(1 + К2^ )1п(1 + Кд_)' (16)
Двумерную плотность распределения длительностей импульсов также положим логарифмически нормальной:
/(Vх/)
2(1-С)
1п‘
1п"
2г
(17)
где гт -параметрраспределения,определяемый, по аналогии с амплитудой, через коэффициент корреляции длительностей:
г =
1п(1 + КТ Кт Я )
______________I ./ /у__________
1п(1 + К2 )1п(1 + К2 )
(18)
Двумерную совместную плотность распределения амплитуд и длительностей также будем считать логарифмически нормальной:
-ехр^ —
2(1--<А.)
1п‘
9 А; 1п -==— А;
2 г
с.
Т:Л,
т Л-1п_*-1п-=А
Ъ АУ
ст, а „
(19)
где гх А -параметрраспределения,определяемый через коэффициент корреляции амплитуд и длительностей Я
1п(1 + Кт КА_ ЯТ.А.)
^1П(1 + ^2)1П(1 + ^2,)'
(20)
Зная корреляционные функции локационных сигналов, можно определить численные значения
IIА , Яг а Ят А для различных конкретных уело-
.а I У
вий работы.
Таким образом, двумерные плотности распределения амплитуд, длительностей, амплитуд и длительностей полностью определены. Во многих практических работах на этом построение модели заканчивается, так как двумерных плотностей распределения вполне достаточно, чтобы определять корреляционные характеристики входных процессов систем обработки информации после любого линейного и нелинейного преобразования. Однако в данной работе построим многомерную совместную плотность распределения амплитуд и длительностей, которую можно использовать не только для приближенного анализа систем обработки информации, но в рамках построенной модели и для точного определения характеристик систем или, при приближенном оценивании характеристик, для определения погрешности соответствующих оценок.
Многомерная математическая модель флюктуаций амплитуд и длительностей
Для построения многомерной функции плотности распределения амплитуд и многомерной функции плотности распределения длительнос-
тей принимаемого сигнала необходимо учесть следующие основные обстоятельства: двумерные плотности распределения должны быть логарифмически нормальные; значение коэффициентов корреляции амплитуд, длительностей, амплитуд и длительностей должны быть равны заданным.
Для амплитуд и длительностей этим условиям удовлетворяют следующие функции [2,9,13,14]:
/дг (А м ) - ■
хехр*^ -
/=1
\____1 ££ , Л, аА.
{ * А, ПА;|’
/лг(тлг)
£=1
1 д(;т)
хехр<--------— УТ—-—1пз-1п^-к
1 2£>^) г=1у=1 ат;°т, Ъ т;{
где определители соответственно равны:
П(А) =
ГАт ГАЫ 2
7)(т) =
-Ь'дг
Г Г ... 1
*N1 2
(21)
(22)
(23)
(24)
~ ~ формулах и £)-ут) - алгебраические дополнения элементов гАц и гт.. в определителях и соответственно. При этом гА и гт определяются по формулам (16) и (18). 17
Для построения многомерной функции плотности распределения амплитуд и длительностей сигнала уже необходимо потребовать выполне-
ния трех условий: многомерная функция плотности распределения амплитуд должна иметь вид (21); многомерная функция плотности распределения длительностей должна иметь вид (22); взаимная корреляционная функция амплитуд и длительностей должна соответствовать заданной.
Всем этим трем условиям удовлетворяет функция
хехр-^ —
(2я)*7АглТП А^а,Ь\ 1=1
1 ££ ВЧ
2£>.
,2ЛП=1МСТА(Ч А А
+-------1п^- +-----— 1п^-1п^- +
аА% ь а, ал, Ъ Л
. , х,, х! Х1
ат ат
т( ту
(25)
у которой определитель вычисляется следующим образом:
Алг -
1 ГА 12 ” Ч* гАч гАх • Л1Х2 " ГАт Л1ТЛ/
Чх 1 •• Ч„ гАч Чт2 Чтлг
Ча • . 1 Ч^т1 Чу т ^
\а гчА ’ • Гт А 1 Гх12 • .. Г тш
V, • ^21 .. Г T2N
Г'ыА • • Гт А • Гт*г • .. 1
,(26)
а И1М - алгебраические дополнения соответствующих элементов определителя Х)2^ • При этом элементы гАт и гхА определяются по формуле (20).
В качестве многомерной плотности распределения /(Адг, т^) и будем использовать логариф-
мически нормальную плотность, определяемую выражениями (25) и (26), так как она удовлетворяет перечисленным выше требованиям, а ее частные случаи (маргинальные плотности распределения) не противоречат экспериментальным данным.
Итак, получена функция плотности распределения амплитуд и длительностей, локационного сигнала, отраженного от объекта, имеющая следующие свойства: маргинальные одномерные функции распределения амплитуд и длительностей сигналов имеют логарифмически нормальное распределение; корреляционная функция амплитуд, корреляционная функция длительностей и взаимная корреляционная функция амплитуд и длительностей могут иметь произвольный функциональный вид, при котором матрицы и |-^2Л^| положительно опреде-
лены (данные функции могут быть заданы и в виде таблиц экспериментальных данных).
Выводы
Математическая модель информационного сигнала представляет собой статистическую модель флюктуаций параметров эхо-сигнала корабля. Эта модель зависит от пяти основных факторов: характеристик объекта, тракта распространения электромагнитных волн, характеристик зондирующего сигнала, условий наблюдения объекта, приемного тракта бортового локатора. Определяющими факторами являются характеристики объекта и вид излучаемого сигнала.
Наиболее распространенными законами распределения параметров модели (амплитуды сигнала) являются: закон Рэлея, Рэлея-Райса, хи-квадрат, закон Накагами, логарифмически нормальный закон. Наиболее целесообразным является использование логарифмически нормальной модели информационного сигнала. В качестве модели флюктуаций длительности отраженных импульсов целесообразно также использовать логарифмически нормальный закон, так как он не только не противоречит экспериментальным данным, но и удобен при построении многомерных моделей.
В качестве многомерного распределения амплитуд и длительностей информационного сигнала используется многомерное логарифмически нормальное распределение, на котором и основана статистическая многомерная математическая модель эхо-сигналов корабля. Построенная модель позволяет учитывать корреляционноспектральные характеристики амплитуд и длительностей сигналов, а также совместные кор-реляционно-спектральные характеристики амплитуд и длительностей.
Литература
1. Андросов В. А., Еиатко И. В. Задачи и принципы построения стендово-имитационной среды для отработки интегрированных комплексов бортового оборудования //Радиотехника. 1996. Вып. 17.С. 120-123.
2. Веетугин А. Р., Шепета А. П. Синтез оптимальных алгоритмов классификации и выбора морских объектов в дальней зоне наблюдения // Сб. науч. трудов «Информационно-управляющие системы и сети. Структуры, моделирование, алгоритмы». СПб.: Политехника, 1999. С. 142-152.
3. Изранцев В. В., Шепета Д. А. Моделирование внешних сигналов бортовых приборных комплексов летательных аппаратов пятого поколения // йзв. вузов. Сер. Приборостроение. 2000. № 2. С.76-83.
4. Исаев С. А., Кондратенко» Г. С. Цифро-натурные и летно-модельные методы испытаний КБО //Радиотехника. 1996. Вып. 17. С. 97-100.
5. Култышев Е. И., Лемешко Н. А., Шепета А. П. Оценка средней амплитуды и средней длительности локационных импульсов, отраженных от надводного объекта// Тр. Междунар. научно-технич. конф. (Тез. докл.), май 1994. Киев: АН Украины; НПО Квант. Вып. 2. С. 5В.
6. Лемешко Н. А. Шепета Д. А. Предельные вероятности распознавания медленно флюктуирующих объектов по средней эффективной отражающей по-верхности// Тр./ СПбГААП. СПб., 1995. С.56-60,
7. Обнаружение радиосигналов / П. С. Акимов,
Ф. Ф. Евстратов, И. С. Захаров и др.; Под ред.
А. А. Колосова. М.: Радио и связь, 1989. 288с.
8. О вероятности выбора главной цели / А. А. Оводенко, А. П. Шепета и др.: Меж-вуз. сб. / ЛЭТИ. Л.у 1977. Вып. 118. С.122-124.
9. Оводенко А. А., Култышев Е. ВЦ Шепета А. П. Бортовая радиоэлектронная аппаратура/ МПЙ. М., 1989. 335 с.
10. Оценка эффективной отражающей поверхности надводной цели по коррелированной выборке / А. В. Г еллер, А. П. Шепета и др.//Тр./ ЛИАП. Л., 1974. Вып. 88. С. 96- 98.
11. Селекция и распознавание на основе локационной информации/А. Л. Горелик, Ю. Л. Барабаш, О. В. Кривошеев, С. С. Эпштейн; Под ред. А. Л. Горелика. М.: Радио и связь, 1990. 240 с.
12. Тверской Г. Н., Терентьев Г. К., Харченко И. П. Имитаторы эхо-сигналов судовых радиолокационных станций. Л.: Судостроение, 1973. 228 с.
13. Шепета Д. А. Моделирование входных сигналов бортовых информационно-измерительных систем. Ассоциация инженерных вузов. Академия космонавтики. МГАТУ им. К. Э. Циолковского// XX Гагаринские чтения: Тез. докл. молодежной научн.-технич, конф., апрель 1994. М.: МГАТУ, 1994. Ч. 4. 67 с.
14. Izrantsev V. V.» Kosenkov А. М., Shepeta D. А. Algorithms For Generating The Non-Gaussian Series With The Given Correlation Characteristics// International Symposium On Problems Of Modular Information Computer Systems And Networks. Abstracts. Moscow; - St.-Petersburg: IEEE, International Informatization Academy, Russian Academy Of Science, Moscow State University, 1997.
\ P. 58.
