_Доклады БГУИР_
2003 июль-сентябрь Том 1, № 3
УДК 621.396.96
МЕТОДИКА АНАЛИЗА ХАРАКТЕРИСТИК СИСТЕМ РАСПОЗНАВАНИЯ НА ОСНОВЕ ПОЛИНОМИАЛЬНОЙ АППРОКСИМАЦИИ МНОГОМЕРНОЙ ПЛОТНОСТИ ВЕРОЯТНОСТИ ВЫХОДНЫХ СИГНАЛОВ
С.Н. ЯРМОЛИК, С В. ШАЛЯПИН
Военная академия Республики Беларусь ВА РБ, Минск, 220057, Беларусь
Поступила в редакцию 29 апреля 2003
Материалы статьи посвящены вопросам статистического анализа характеристик систем радиолокационного распознавания. Предложена методика аналитического расчета вероятностных показателей качества систем распознавания, основу которой составляет эффективная аппроксимация закона распределения межканальных разностей с помощью усеченного ряда по ортогональным полиномам Поллачека. В качестве весовой функции предложено использовать функцию Поллачека, адаптированную к искомому закону распределения, что позволило получить высокую точность при использовании небольшого числа членов ряда.
Ключевые слова: распознавание, плотность вероятности, ортогональные полиномы.
Введение
Согласно теории статистических решений, результат классификации объектов принимается с некоторой вероятностью ошибки, которая характеризует качество работы устройства. Работоспособность систем распознавания принято характеризовать с помощью набора вероятностей принятия правильных и ложных решений [1]:
рк/я = Р(К1\), к,^...м
где Р( А*к! А3) — условная вероятность принятия решения в пользу к-го класса, если на входе
устройства распознавания присутствует объект £-го класса.
При построении байесовских систем радиолокационного распознавания, обычно, используют следующее решающее правило [1]:
если 1к1 = ак1 + £БИТ > 0 для всех I = 1,...,М, IФ к ^ А*, (1)
где Zн — межканальные разности результатов обработки реализации портрета; Би — матрицы обработки входной гауссовской реализации портрета ^; ак1 — межканальные смещения.
Такой вариант представления алгоритма работы устройства распознавания позволяет определять условные вероятности принятия решения Рк/8 путем (М-1)-кратного интегрирования в пределах от 0 до +<» совместной плотности вероятностей М-1 межканальных разностей
+ад +ад +ад +ад
= | Pk/g(Zkl)dZkl = | Pk/g(Zk1>•••>Zk,k-1>Zk,k+1>•••>Zkм)dZk1 • • •^к,к-1^к,к+1 • • '^кЫ 0 0 0 0
(k = 1,м , g = 1,м , I = 1,м , I * k), (2)
где Pk/g( Zkl) = Pk/g(Zkl > •■■> Zk,k-l' Zk,k+l' ••■> ^м) — закон распределения случайных межканальных разностей Zkl при наличии сигнала отраженного от объекта g-го класса.
Таким образом, в основе нахождения вероятностных характеристик распознавания лежит вычисление многомерного интеграла от закона распределения случайных межканальных разностей. Очевидна сложность самого процесса интегрирования; кроме того, часто не удается аналитически выразить подынтегральную функцию.
Аппроксимация плотности вероятности выходных сигналов каналов обработки системы распознавания
Совместная плотность вероятности pk/g(Zk1/g • ..Zш/g) может быть представлена с помощью системы ортогональных полиномов:
ад ад
Pk/g( Ъ) = Фl(Zkl/g)•■ ■фm(ZШ/g)Е1••■ ЕCq1...qмQq1(Zk1/g)••■QqM(ZkM/g) , (3)
q1=0 qM=0
где фl(Zkl/ ) — весовая функция используемого полинома (/ =1,..м, l * к);
ql
Qql(Zkl/g) = Е— ортогональный полином q/-й степени; г/ — коэффициенты
p=0
ортогонального полинома; сх м — коэффициенты разложения в ряд. Коэффициенты разложения могут быть представлены в виде
q1 qM
Cq1...qM = Е^ ' ' Еrq1/s1 • •rqM/sMm{Zlc1/g ••■ZШ/g} , ,ч1=0 м=0
где m{ZSkl/g •■ •ZkM/g} математическое ожидание системы величин (Zkl/g );
показатель степени Zkl/g.
Выражение для плотности вероятности (3) может быть преобразовано как
Pk/g(Zk1/g ••■Zm/g) =
ад ад q1 qM q1 qM м
= !••• ЕЕ- ЕЕ- Е ^ • • rql/sl rql/plфl(Zkl/g)Zkl/g . (4)
q1=0 qM=0.1=0 sM=0 p1=0 pM=0 Ы
ы
Исходя из гауссовости закона распределения элементов обрабатываемого портрета , зная его функциональное преобразование, можно определить значения моментов искомого закона распределения:
Е .
тт1^--^,,} = (-*Г ^ ,,м ©(V)
м
V, =0
М / М
где 0г(V) = exp(^Аш) (2п)ыdet[I-/^\1 Хы] — многомерная характеристическая
i=i
ы
i=i
ы
функция величин Z; Хы = К3+0 Бы — определяющая матрица; V = {V п..., V к _1, V к+1,.., V М } —
41
совокупность вещественных переменных. Математические ожидания
~ £Ако£*д _ т{%ккХд } 2к =-кос.д *-имеют вид
центрированных
случайных
величин
а{ÇAkoÇ}
т{^1} = 0 , = *Г( ХХ2) , т{11,12,13} = 2&( Ш2Х3) ,
т{2~1,2:2,2:3,2:4} = ыг( х^хх^+*г( Ш2М х3х4)+*г( х{х3)^( х2х4)+*г( Ш4М х2х3),
где х) — след матрицы х .
Учитывая ортогональность полиномов, выражения для условных вероятностей принятия решений (4) примут вид
N
N q1 qM q1 qM
II ■■• II ■■• ъ
q1=0 qM=0 s1=0 sM=0p1=0 pM=0
Fk/g = I— I !■ ' ' II"' hm{Zk1/g ■ ■ZkM/g}^\rql/sl rql/pl
l=1 l*k
¡<p,(Za/g)Zp/gdZt
'kl/g
(5)
Выбор семейства ортогональных полиномов
Для качественной аппроксимации важно правильно выбрать семейство ортогональных полиномов, что позволит использовать при расчетах небольшое число членов усеченного ряда N [2].
При разложении было использовано семейство ортонормированных полиномов Полла-чека [4], представленных через степенные моменты своей весовой функции:
Qn(x) =
Ja , -А
n-1 n
h0 h1 h2 h1 h2 h3
hn-1 hn hn+1
X X
hn
hn,
2n-1 Xn
где An =
h0 h1 h1 h2
hn hn+1
n+1
h2
определитель Грама; hn = jф(x) • xndx
x dx — n-и момент весовой
функции полинома.
Весовая функция Поллачека имеет вид (рис. 1)
ф(х) = h( x,аД ) = — ( 2 kin( а )) п
( 2А-1 ) п-2а )x
|Ä(Х + ix)2,
где а, X — параметры весовой функции (0<а<п, X>0).
-mZ
а Z
1
1
h
n
Рис. 1. График функции Поллачека
Для аппроксимации используется адаптированная к параметрам аппроксимируемого закона распределения функция Поллачека с оптимально подобранными параметрами а , и X ,, которые определяются из системы уравнений:
+ад +ад м
|(к (х, а, X)(х - |к (х, а, X)хс1х)2)с1х = т(2){Zsk11/g, ... ^м[/в}, ааа Е. = 2,
—ад —ад ,=1
+ад +ад м
|(к1 (х, а, X)(х — |к (х, а, X)хСх)3)сСх = т(з){Zsk11/g, •, ZsЫ^[/g }, ааа Е. = 3.
^—ад —ад ,=1
В этом случае качественная аппроксимация плотности вероятности обеспечивается минимальным числом членов ряда. На рис. 2 в качестве примера приведены рассчитанные по предложенной методике вероятности принятия решений при обнаружении-распознавании объектов по 10-элементному частотно-коррелированному портрету от отношения сигнал/шум у
(А/ = 1 МГц).
16 18 20 П
Рис. 2. Характеристики обнаружения-распознавания объектов трех классов
Заключение
Использование ортонормированных полиномов Поллачека с весовой функцией, адаптированной к параметрам аппроксимируемого закона распределения, позволяет обеспечить качественную аппроксимацию при использовании небольшого числа членов ряда. Такой подход позволяет учесть корреляцию межканальных разностей и обеспечивает высокую точность получения вероятностных показателей качества.
Предлагаемая методика анализа может использоваться разработчиками радиотехнических систем, работающими в области обнаружения и (или) распознавания.
TECHNIQUE OF THE RECOGNITION SYSTEMS CHARACTERISTICS ANALYSIS ON THE BASIS OF MULTIVARIANT PROBABILITY DENSITY OF OUTPUT SIGNALS POLYNOMIAL APPROXIMATION
S.N. YARMOLIK, S.V. SHALIAPIN Abstract
The materials of this article are devoted to questions of the statistical analysis of recognition radar systems. The technique of analytical calculation of recognition systems probabilistic characteristic based on interchannel differences distribution law effective approximation by using truncated line on Pollachek orthogonal polynoms is developed. As weight function it is offered to use Pollachek function adapted to the required law of distribution that has allowed receiving of high accuracy by using low term of series count.
Литература
1. Охрименко А.Е. Основы радиолокации и РЭБ. Ч. 1. Основы радиолокации. М., 1983.
2. Сеге Г. Ортогональные многочлены. М., 1962.
3. PollaczekF. // Comptes Rendus de L'Acad. des Sc. Paris, Vol. 230. 1950. P. 1563-1565.