- •СОДЕРЖАНИЕ
- •ВВЕДЕНИЕ
- •1 НАЗНАЧЕНИЕ И ОСОБЕННОСТИ РАДИОТЕХНИЧЕСКИХ СИСТЕМ [1, 2, 3, 4]
- •1.1 Общая характеристика возможностей РТС
- •1.2 Основные определения и системные принципы
- •1.3 Жизненный цикл РТС
- •1.4 Классификация радиотехнических систем
- •1.5 Эффективность радиосистем
- •1.6 Основные показатели РТС
- •1.7 Повышение эффективности РТС и РЭА при проектировании
- •1.8 Характеристики сигналов и помех в РТС
- •1.9 Особенности задач оптимизации РТС
- •2.2 Оптимальное обнаружение детерминированных сигналов
- •2.2.1 Постановка задачи
- •2.2.2 Методика решения задачи обнаружения
- •2.2.3 Корреляционный оптимальный обнаружитель
- •2.2.4 Качественные показатели оптимального обнаружителя
- •2.2.5 Согласованные фильтры в оптимальных обнаружителях
- •2.2.6 Реализация согласованных фильтров
- •2.3 Обнаружение сигналов со случайными параметрами
- •2.3.1 Постановка задачи
- •2.4 Обнаружение пачек импульсов
- •2.4.2 Структурные схемы обнаружения пачек импульсов
- •2.4.3 Накопители импульсных сигналов
- •2.4.4 Характеристики обнаружения при накоплении импульсов
- •2.5 Различение сигналов
- •2.5.1 Различение детерминированных сигналов
- •2.5.2 Различение двоичных сигналов со случайными параметрами
- •2.5.3 Различение ансамблей сигналов
- •2.7 Функция неопределённости сигнала
- •3.1 Оптимальное разрешение сигналов
- •3.2 Оценка разрешающей способности
- •3.3 Диаграммы неопределенности сигналов
58
функции gopt (t) ограничена частотойFв , то по теореме Котельникова,
необходимый интервал Dt должен быть не более 1 (2Fв ) .
Рисунок 2.12 – Структурная схема согласованного фильтра на линии задержки.
Коэффициенты в блоке весовых коэффициентов выбираются в соответствии с величинами gk º gopt (kDt) , k =1,..., M . Структурная схема
согласованного фильтра на линии задержки показана на рисунке2.12. Оконечный согласованный фильтр согласован с прямоугольным импульсом, имеющим длительность Dt .
2.3Обнаружение сигналов со случайными параметрами
2.3.1 Постановка задачи
Рассмотрим задачу обнаружения сигнала, зависящего от случайных неизмеряемых параметров. Примерами такого рода сигналов являются сигнал со случайной начальной фазой, модель которого даётся выражением(1.19), сигнал со случайной начальной фазой и амплитудой(1.20), некогерентные пачки радиоимпульсов (1.21) и (1.22).
Вобщем случае при наличии в сигнале случайных параметров
наблюдаемую при обнаружении смесьr(t) можно представить в виде
r |
r |
r(t) = qs(t, b ) + n(t) , где |
b - вектор случайных параметров сигнала; |
r
плотность вероятности w(b ) предполагается неизвестной. Для определения алгоритма оптимального обнаружения необходимо вычислить отношение правдоподобия и, воспользовавшись, решающим правилом (2.20), выявить структуру оптимального обнаружителя и его качественные показатели.
|
59 |
|
|
|
2.3.2 Отношение |
правдоподобия |
при |
наличии |
случайны |
параметров сигнала |
|
|
|
|
Рассмотрим вначале |
совместную плотность вероятности дискретных |
|||
значений rk принимаемого колебания и случайных неизмеряемых параметров
r
b = b1,..., bn . Согласно теореме умножения вероятностей, можно записать:
|
r |
|
|
r |
|
|
r |
|
r |
|
|
(2.45) |
w(r, b ) = w(r)w(b | r) = w(b )w(r | b ) . |
|
|
||||||||||
Поскольку, |
по |
|
условию |
нормировки, интегрирование |
условной |
|||||||
плотности вероятности в бесконечных пределах всегда даёт единицу, после |
||||||||||||
интегрирования |
левой |
и |
право частей выражения(2.45) |
по всем |
значениям |
|||||||
r |
|
|
|
|
|
|
|
|
|
|
|
|
b = b1,..., bn получим |
|
r |
|
r |
|
|
|
|
|
|||
|
¥ |
r |
|
|
|
|
|
|
|
|
||
w(r) = ò w(b )w(r | b )d b . |
|
|
|
|
(2.46) |
|||||||
|
-¥ |
|
|
|
|
|
|
|
|
|
|
|
По условию, w(r) = wс.п (r) |
|
- плотность |
вероятности |
реализации величин |
||||||||
r1 ,..., rm при наличии суммы сигнала и помехи. Разделив обе части равенства |
||||||||||||
(2.46) на |
плотность |
вероятности реализации этих величин при условии |
||||||||||
действия одной помехи wп (r) , получим, согласно формуле (2.16), следующее |
||||||||||||
выражение для отношения правдоподобия: |
|
|
|
|
||||||||
|
¥ |
r |
|
|
r |
|
r |
|
|
|
|
|
|
|
w(r | b) |
|
|
|
|
|
|
||||
l(r) = ò w(b ) |
|
|
|
d b . |
|
|
|
|
|
|||
wп (r) |
|
|
|
|
|
|||||||
|
-¥ |
|
|
|
|
|
|
|
|
|||
Отношение плотностей вероятности под знаком интеграла можно |
||||||||||||
рассматривать |
как условное |
(частное) |
отношение |
|
правдоподобия при |
|||||||
фиксированных значениях параметров b1 ,..., bn . При этом сигнал становится |
||||||||||||
полностью |
известным, |
и |
|
для |
условного |
отношения |
правдоподобия |
|||||
справедливы выражения (2.18) и (2.19).
Чтобы получить искомое отношение правдоподобия, входящее в оптимальное решающее правило(2.20), необходимо усреднить частное отношение правдоподобия по всем значениям случайных неизмеряемых параметров с учётом их плотности вероятности:
¥ |
r |
r |
r |
l(r) = ò w(b )l(r | b )d b |
|||
-¥ |
|
|
|
или |
|
r |
r r |
¥ |
|
||
l(r(t)) = ò w(b)l(r(t) | b)d b . |
|||
-¥ |
|
|
|
Условное |
отношение |
правдоподобия в соответствии с выражением |
|
(2.19) можно представить в виде
60
|
r |
æ |
|
r |
|
|
æ 2 |
|
|
r ö |
|
|
|
|
|
|
|
|
E(b) ö |
|
|
|
|
|
|
|
|||||||
l(r(t) | b) = exp ç |
- |
|
|
÷exp ç |
|
z(r(t) | b ) ÷, |
|
|
|
|
|
|||||
N0 |
|
N0 |
|
|
|
|
|
|||||||||
где |
|
è |
|
|
ø |
è |
|
|
ø |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
T |
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
E(b ) = òs2 (t, b )dt ; |
|
|
|
|
|
|
|
|
|
(2.47) |
|
|
||||
|
0 |
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z(r(t) | b ) = òr(t)s(t, b )dt . |
|
|
|
|
|
|
|
(2.48) |
|
|
||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Окончательное |
выражение |
|
для |
отношение |
правдоподобия |
при |
||||||||||
произвольной |
плотности |
|
вероятности |
случайных |
параметров |
сигнала |
||||||||||
принимает вид |
r |
|
|
|
r |
|
|
|
|
r |
|
r |
|
|
|
|
|
¥ |
|
æ |
|
|
|
æ 2 |
ö |
|
|
|
|||||
|
|
|
E(b) ö |
|
|
|
|
|||||||||
l(r(t)) = ò w(b ) exp ç - |
|
|
÷exp ç |
|
z(r(t) | b ) ÷d b |
(2.49) |
|
|
||||||||
|
N0 |
|
|
|
||||||||||||
|
-¥ |
|
|
è |
|
ø |
|
è N0 |
|
ø |
|
|
|
|
||
Это выражение можно конкретизировать для определённых моделей |
|
|||||||||||||||
сигналов. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2.3.3 |
Отношение |
|
правдоподобия |
для |
сигнала |
со |
случайной |
|||||||||
начальной фазой |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Пусть обнаруживаемый сигнал задан в виде(1.9). Примем здесь t0 = 0 |
|
|||||||||||||||
и будем |
считать фазуj º b |
равномерно |
распределённой |
на интервале |
|
|||||||||||
[-p ,p ] с плотностью вероятности w(b ) =1
(2p ).
Для того чтобы воспользоваться общей формулой(2.49), вычислим предварительно энергию (2.47) и корреляционный интеграл(2.48) при фиксированном значении параметра b .
По формуле косинуса суммы, сигнал s(t, b ) [см. выражение (1.9)] представив в виде
s(t, b ) = S0 (t) cos(w0t +y (t)) cos b - -S0 (t) sin(w0t +y (t)) sin b .
Обозначим:
s1 (t) = S0 (t) cos(w0t +y (t)) ; s2 (t) = -S0 (t) sin(w0t +y (t))
и запишем сигнал s(t, b ) с учётом введённых обозначений |
|
s(t, b ) = s1 (t) cos b + s2 (t) sin b . |
(2.50) |
На основании этого выражения условное значение корреляционного интеграла можно записать следующим образом;
61
T
z(r(t) | b ) = òr(t)s(t, b )dt = z1 cos b + z2 sin b .
0
Здесь
T
zi = òr(t)si (t)dt = zi (r(t)), i =1, 2 |
(2.51) |
||||
|
0 |
|
|
|
|
Введём величины Z и Q , определив их соотношениями: |
|
||||
|
|
|
|
|
|
Z = |
|
z2 + z2 ; |
(2.52) |
||
|
1 |
2 |
|
|
|
cos Q = z1 |
Z ;sin Q = z2 Z . |
(2.53) |
|||
С учётом этих величин условное значение корреляционного интеграла |
|||||
приводится к виду |
|
||||
z(r(t) | b ) = Z (cos b cos Q + sin b sin Q) = |
|
||||
= Z cos(b - Q). |
(2.54) |
||||
Найдём |
условное значение энергииE(b ) сигнала. При |
малом |
|||
изменении |
|
амплитуды S0 (t) и фазы y (t) за период колебаний |
высокой |
||
частоты энергия практически не зависти от начальной фазы b , т.е.
T
E(b ) = òS02 (t) cos2 (w0t +y (t) + b )dt @
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ |
1 |
T S02 (t)dt = E . |
|
|
|
|
|
|
|
|
|
|
(2.55) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
2 |
ò0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Подставив |
полученные |
выражения |
для корреляционного интеграла |
|||||||||||||||
(2.54) и энергии (2.55) в |
формулу (2.49) и |
учитывая плотность |
вероятности |
|||||||||||||||
w(b ) , найдём выражение для безусловного отношения правдоподобия: |
||||||||||||||||||
|
|
|
|
|
æ |
|
E ö |
1 |
|
2p |
æ 2Z |
|
ö |
|
||||
l(r(t)) = exp ç |
- |
|
÷ |
|
|
|
ò exp ç |
|
cos(b - Q) ÷d b . |
|
||||||||
|
|
|
|
N0 |
|
|||||||||||||
|
|
|
|
|
è |
|
N0 ø 2p |
0 |
|
è |
|
ø |
|
|||||
Известно, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
1 |
|
2p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
exp(u cos(b - Q))d b = I0 (u) |
|
|
(2.56) |
|||||||||||
|
|
|
|
|
|
|||||||||||||
|
2p ò0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
есть модифицированная функция Бесселя нулевого порядка. Это чётная |
||||||||||||||||||
функция своего аргумента, причём |
|
I0 (0) =1 (рисунок 2.13). |
Учитывая |
|||||||||||||||
формулу |
(2.56), |
запишем |
отношение |
|
правдоподобия для сигнала со |
|||||||||||||
случайной начальной фазой в окончательном виде: |
|
|||||||||||||||||
|
|
|
|
|
æ |
|
E ö |
æ 2Z ö |
|
|
|
|
|
|||||
l(r(t)) = exp ç |
- |
|
÷ I |
0 ç |
|
÷, |
|
|
|
|
(2.57) |
|||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
è |
|
N0 ø |
è |
N0 ø |
|
|
|
|
|
||||
62
где Z определяется из соотношений (2.51) и (2.52).
2.3.4 Отношение |
правдоподобия |
для сигнала со случайными |
начальной фазой и амплитудой |
|
|
Предположим, что обнаруживаемый сигнал задан выражением(1.10), в |
||
котором t0 = 0. Пусть |
начальная фазаj º b |
равномерно распределена на |
интервале [-p ,p ], а амплитудный коэффициент A представляет случайную
величину с плотностью вероятности Рэлея при среднеквадратичном значении, равном единице.
Рисунок 2.13
Для независимых величин A и b совместная плотность вероятности определяется следующим образом:
w(b, A) = 1 2 A exp(-A2 ), A ³ 0 ; 2p
b Î(-p ,p ) .
Вычислим корреляционный интеграл z(r(t) | b, A) и энергию E(b, A) для фиксированных значений параметровb и A. По аналогии с .п2.3.3.
можно получить:
z(r(t) | b, A) = AZ cos(b - Q) ;
E(b, A) = A2 E ,
где выражения для Z , Q и E совпадают с аналогичными в п. 2.3.3. В отличие от рассмотренного случая сигнала со случайной начальной фазой энергия E(b, A) зависит от случайного амплитудного множителяA. Среднее значение энергии определяется следующим образом:
¥ |
2p |
¥ |
Eср = òdA ò E(b, A)w(b, A)d b = E ò2A3 exp(-A2 )dA = E , |
||
0 |
0 |
0 |
т.е. среднее |
значение равно энергии сигнала при частном значении |
|
A =1.
Используя общую формулу (2.49), получим выражение для отношения правдоподобия в виде
63
|
1 |
¥ |
2p |
|
|
|
æ |
|
|
|
A2 E ö |
æ 2 AZ |
|
ö |
||||||
l(r(t)) = |
|
òdA ò A exp(-A2 ) exp ç - |
|
|
|
|
|
÷exp ç |
|
cos(b - Q) ÷d b . |
||||||||||
|
|
|
|
|
|
N0 |
||||||||||||||
|
p |
0 |
0 |
|
|
|
è |
|
|
|
|
N0 ø |
è |
|
ø |
|||||
Применив формулу (2.56), а также табличный интеграл |
|
|||||||||||||||||||
¥ |
|
|
|
|
|
|
1 |
|
æ |
m |
2 |
ö |
|
|
|
|
||||
òI0 (mx) exp(-n x2 )xdx = |
exp |
|
|
|
|
|||||||||||||||
ç |
|
÷, |
|
|
|
(2.58) |
||||||||||||||
|
|
|
|
|
|
|||||||||||||||
0 |
|
|
|
|
|
|
2n |
è |
4n |
ø |
|
|
|
|
||||||
получим |
окончательное |
выражение |
отношения |
правдоподобия для |
||||||||||||||||
сигнала со случайными начальной фазой и амплитудой: |
|
|
||||||||||||||||||
l(r(t)) = |
N0 |
æ |
1 |
|
Z 2 |
|
|
ö |
|
|
|
|
|
|
|
|
||||
|
exp ç |
|
|
|
|
|
|
÷, |
|
|
|
|
|
|
(2.59) |
|||||
E + N0 |
N0 |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
è |
E + N0 ø |
|
|
|
|
|
|
|
|
||||||||
где величина Z определяется по формулам (2.51) и (2.52). |
|
|
||||
2.3.5 Структурные схемы оптимальных обнаружителей |
сигналов |
|||||
со случайными параметрами |
|
|
|
|
||
Алгоритмы оптимального обнаружения рассматриваемых сигналов |
||||||
реализуются |
на |
основании |
полученных |
выражений |
для |
отношени |
правдоподобия (2.57) |
и (2.59) и |
правила (2.20). Решающее правило можно |
||||
преобразовать, имея в виду, что сравнение с порогом остаётся в силе для любой монотонной функции от отношения правдоподобия.
Рисунок 2.14 – Структурная схема квадратурного приёмника На этом основании запишем оптимальное решающее правило:
qopt* |
ì1, |
åñëè |
Z > Z |
ï |
; |
|
(r(t)) = í |
åñëè |
|
|
(2.60) |
||
|
î0, |
Z < Zï . |
|
|||
Используя выражение (2.60), а так же |
соотношения определяющие |
|||||
величину Z , структурную схему обнаружителя можно представить в виде, изображённом на рисунке2.14. Такую схему называютквадратурным приёмником. Квадратурные каналы организуют благодаря включению
фазовращателя в цепь опорного сигнала одного |
из |
перемножителей. |
|||
Квадраторы (Кв), сумматор |
( S), вычислитель квадратного |
|
|
|
|
корня( |
) |
||||
обеспечивают формирование |
на входе ПУ в момент времениt = T |
значения |
|||
|
|
|
64 |
|
|
|
|
Z = z2 |
+ z2 |
. Наличие двух каналов исключает возможность потери |
|
1 |
2 |
|
|
полезного сигнала вследствие незнания его начальной фазы. При отсутствии приращения корреляционного интеграла, например, в первом канале, второй канал обеспечивает приращение за счёт сдвига фазы опорного сигнала на p
2 . Схема, представленная на рисунке 2.14, как и схема, изображённая на рисунке 2.1, требует знания временного положения ожидаемого сигнала. Если время запаздывания сигнала неизвестно, оптимальный обнаружитель
усложняется. |
При |
разбиении |
интервала |
неопределённости |
времени |
запаздывания |
на элементарные участки, длительность каждого из которых |
||||
определяется требуемой разрешающей способностью, можно построить многоканальный обнаружитель. Каждый канал его настраивается на сигнал с соответствующим ожидаемым запаздыванием. Решение об обнаружении сигнала принимается одновременно с оценкой времени запаздывания.
Вместо |
схемы с |
квадратурными каналами используют схему с |
||||
согласованным |
фильтром. |
Возможность |
использования |
согласованных |
||
фильтров в обнаружителях сигналов, имеющих случайные параметры, |
||||||
основывается на следующих |
|
рассуждениях. Из соотношений (2.50) и (2.51) |
||||
следует необходимость формирования квадратурных компонентов сигнала. |
||||||
Представим эти компоненты в следующем виде: |
|
|
||||
|
g |
|
g |
|
t)) ; |
|
s (t) = Re(S 0 (t) exp( jw |
t)) = Re(S* (t) exp(- jw |
|
||||
1 |
|
0 |
0 |
0 |
|
|
Рисунок 2.15 – Структурная схема обнаружителя сигнала со случайной фазой (или со случайными фазой и амплитудой)
g |
(t)exp( jw |
g |
* (t)exp(- jw t)) , |
|
s (t) = -Im(S 0 |
t)) = Im(S |
|||
2 |
0 |
0 |
0 |
|
65
|
|
· |
|
|
где |
комплексная огибающаяS0 (t) = S0 (t) exp( jy (t)). Знак «*» |
|
||
означает |
комплексно-сопряжённую |
функцию. В |
соответствии |
с |
соотношениями (2.51) и приведёнными выражениями запишем: |
|
|||
æ T z1 = Re ç ò
è 0
æ T z2 = Im ç ò
è 0
· * |
ö |
r(t) S 0 (t) exp(- jw0t)dt ÷; |
|
|
ø |
· * |
ö |
r(t) S 0 (t) exp(- jw0t)dt ÷ .
ø
Величина z представляет модуль комплексного числа:
|
|
|
T |
· * |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
z = |
z12 + z22 |
= |
òr(t) S 0 (t) exp(- jw0t)dt |
. |
(2.61) |
|||||
|
|
|
0 |
|
|
|
|
|
|
|
Используем очевидное равенство |
1 |
|
|
· * |
||||||
|
|
|
|
|
|
|
· |
|
||
r(t) = R(t) cos(w t +y |
R |
(t)) = |
|
(R(t)e jw0t + R (t)) ×e- jw0t . |
||||||
|
||||||||||
|
|
0 |
|
|
2 |
|
|
|
||
При этом подынтегральное выражение в формуле(2.61) представим в |
||||||||||
виде суммы: |
· * |
· * |
|
|
|
|
|
|||
· |
· |
|
|
|
|
|
||||
0, 5(R(t) S 0 (t) + R (t) S 0 (t) exp(- j2w0t)) ,
в которой второе слагаемое представляет быстро осциллирующую функцию с частотой2w0 . Интеграл от такой функции можно считать приближённо равным нулю. Поэтому
|
1 |
|
T · |
· * |
|
|
|
z ; |
|
òR(t) S 0 (t)dt |
. |
(2.62) |
|||
2 |
|||||||
|
|
0 |
|
|
|
||
Полученному выражению ставится в соответствие схема, содержащая
· *
согласованный фильтр (СФ) с импульсной реакцией S 0 (t) и включенный после СФ линейный детектор огибающей (ДО), который реализует операцию взятия модуля. Структурная схема обнаружителя сигнала со случайной фазой или со случайными фазой и амплитудой представлена на рисунке2.15, а, временные диаграммы в точках 1-3, поясняющие работу схемы, - на рисунке
2.15, б.
2.3.6 |
Качественные |
показатели |
обнаружения |
сигналов |
с |
||
случайными параметрами |
|
|
|
|
|
|
|
Для определения характеристик обнаружения сигналов со случайными |
|
||||||
параметрами необходимо знать условные плотности вероятности случайной |
|
||||||
величины |
Z при наличии |
и |
отсутствии сигнала |
в наблюдаемой . смеси |
|||
Плотность |
вероятности w(Z ) |
можно |
определить |
из |
выражения |
для |
|
|
|
|
|
|
|
66 |
|
|
|
|
|
|
совместной |
плотности вероятности w(Z , Q) , |
проинтегрировав |
её по всем |
|
||||||||
возможным |
значениям Q . |
Совместная плотность вероятностиw(Z , Q) |
|
|||||||||
вычисляется |
на |
|
основании совместной плотности вероятностиw(z1 , z2 ) |
|
||||||||
путём перехода к новым переменным с учётом якобиана преобразования. |
|
|||||||||||
Рассмотрим обнаружитель сигнала, содержащего случайную фразу b и |
|
|||||||||||
случайную |
амплитуду A, |
подчиняющуюся |
распределению |
Рэлея. При |
|
|||||||
известной амплитуде сигнала плотность |
вероятностиw( A) = d ( A -1) , т.е. |
|
||||||||||
принимается |
A =1. При отсутствии сигнала w( A) = d ( A) , |
т.е. принимается |
|
|||||||||
A = 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
Можно показать по аналогии с соотношениями(2.23) и (2.24), что для |
|
|||||||||||
любого фиксированного значения A средние значения случайных величин z1 |
|
|||||||||||
и z2 и их дисперсии определяются выражениями: |
|
|
|
|
||||||||
|
T |
|
|
|
|
ìAE cos b; |
|
|
|
|
|
|
< z1,2 >= Aòs(t, b)s1,2 (t)dt = |
|
|
|
|
|
|||||||
í |
|
|
|
|
|
|
||||||
|
0 |
T T |
|
|
îAE sin b; |
|
|
|
|
|
||
2 |
2 |
òò< n(t1 )n(t2 ) > s1,2 (t1 )s1,2 (t2 )dt1dt2 |
= N0 E 2 . |
|
|
|||||||
sz1,2 º s z = |
|
|
||||||||||
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
Случайные |
|
величины z1 |
и |
z2 |
некоррелированы |
и |
подчиняются |
|||||
гауссовскому |
закону |
распределения. Следовательно, эти |
величины |
|||||||||
статистически |
|
независимы. |
Некоррелированность |
|
вытекает |
из |
||||||
приближенного |
|
равенства |
|
нулю |
смешанного |
центрального |
момента |
|||||
< (z1 - < z1 >)(z2 - < z2 >) >. |
Действительно, |
этот |
момент |
сводится |
к |
|||||||
интегралу от произведения двух колебаний s1 (t) и s2 (t) , сдвинутых по фазе на p
2 , при этом интеграл практически равен нулю.
Совместная условная плотность вероятности случайных величинz1 и z2 с учётом их независимости определяется выражением
w(z1 , z2 | b, A) = w(z1 | b, A)w(z2 | b, A) =
|
1 |
æ (z - < z >)2 + (z |
2 |
- < z |
2 |
>)2 |
ö |
||||||
= |
|
exp ç |
1 |
1 |
|
|
|
|
|
÷. |
|||
2ps z2 |
|
|
|
2s z2 |
|
|
|
|
|
|
|||
|
è |
|
|
|
|
|
|
|
|
|
ø |
||
Полагая на основании выражения(2.53) |
z1 = Z cos Q, z2 = Z sin Q , |
||||||||||||
перейдём к новым переменным двухмерной плотности и вероятности |
|||||||||||||
w(Z , Q | b, A) = w(z , z |
2 |
| b, A) |
¶(z1 , z2 ) |
, |
|
||||||||
|
|
||||||||||||
|
|
|
|
1 |
|
|
|
¶(Z , Q) |
|
||||
|
|
|
|
|
|
|
|
|
|
||||
где
67
¶(z1 , z2 ) |
= |
cos Q |
sin Q |
= Z |
|
-Z sin Q |
Z cos Q |
||
¶(Z , Q) |
|
|||
есть якобиан преобразования. После преобразования получим
w(Z , Q | b, A) = |
Z |
exp(- |
1 |
- (Z 2 + A2 E2 - |
|
2ps z2 |
2sz2 |
||||
|
|
|
-2AEZ cos(b - Q))) .
Одномерная условная плотность вероятности w(Z | b, A) вычисляется путём интегрирования правой части полученного выражения по Q :
|
|
|
|
2p |
|
|
|
Z |
æ AEZ ö |
|
||
w(Z | b, A) = ò w(Z , Q | b, A)dQ = |
|
|||||||||||
|
I0 ç |
|
÷ |
´ |
||||||||
2 |
2 |
|||||||||||
|
|
|
|
0 |
|
|
|
sz |
è |
sz |
ø |
|
æ |
|
Z |
2 |
2 |
E |
2 ö |
|
|
|
|
|
|
´exp ç |
- |
|
+ A |
|
÷. |
|
|
|
|
|
||
|
|
2sz2 |
|
|
|
|
|
|
|
|||
è |
|
|
|
|
ø |
|
|
|
|
|
||
Произведя усреднение по параметрамb и A с учётом плотности вероятности w(b, A) , найдём искомую плотность вероятности величины Z :
|
1 |
2p |
¥ · |
|
w(Z ) = |
ò d b òw(Z | b, A)w( A)dA. |
|||
2p |
||||
|
0 |
0 |
||
Вид этой плотности вероятности меняется в зависимости от заданной плотности вероятности w( A) , поэтому целесообразно рассмотреть частные
случаи: |
|
|
|
|
|
|
|
|
|
|
1) при w( A) = d ( A -1) (сигнал со |
|
случайной |
начальной фазой и |
|||||||
неизвестной амплитудой) |
|
|
|
|
|
|
|
|||
|
|
Z |
æ |
|
Z 2 + E 2 ö æ ZE ö |
|
|
|||
wс.п |
(Z ) = |
|
exp ç |
- |
|
÷ I0 ç |
|
÷ |
, Z > 0 ; |
(2.63) |
sz2 |
2s z2 |
|
||||||||
|
|
è |
|
ø è s z2 |
ø |
|
|
|||
Рисунок 2.16 – Кривые условных плотностей вероятностиwп (Z ) и wс.п (Z ) при обнаружении сигнала со случайной начальной фазой.
68
Рисунок 2.17 - Семейство характеристик обнаружения.
2) при w( A) = 2 A exp(-A2 ) - рэлеевском распределении (сигнал со случайными амплитудой и начальной фазой)
wс.п |
(Z ) = |
|
2Z |
|
|
|
|
æ |
|
Z 2 |
ö |
, Z > 0 ; |
||
|
|
|
|
|
exp ç |
- |
|
÷ |
||||||
2s z2 + E 2 |
2sz2 + E2 |
|||||||||||||
|
|
|
|
è |
|
ø |
|
|||||||
3) при w( A) = d ( A) (сигнал отсутствует) |
||||||||||||||
wп (Z ) = |
Z |
exp |
æ |
- |
|
Z 2 |
ö |
, Z > 0 . |
|
|
||||
|
|
ç |
|
|
÷ |
|
|
|||||||
sz2 |
|
2sz2 |
|
|
||||||||||
|
|
|
è |
|
|
ø |
|
|
|
|
||||
Здесь при интегрировании использована табличная формула(2.58).
Кривые |
условных |
плотностей |
вероятностиw (Z ) |
и w |
(Z ) |
при |
|
|
|
п |
с.п |
|
|
обнаружении сигнала со случайной начальной фазой приведены на рисунке 2.16. Заштрихованные площади под кривыми справа от прямой, проходящей через пороговое значение Zп , соответствуют вероятностям ложной тревоги F и правильного обнаружения D .
Рисунок 2.18
Введём переменную n = Z
sz и относительный порогhп = Zп
s z .
Тогда с учётом параметра обнаружения q = 2E
N0 (отношения сигнал/шум) получим:
¥ |
|
æ |
|
q +n |
2 ö |
|
D = òn I0 ( |
|
|
||||
|
|
|||||
qn ) exp ç |
- |
|
÷dn ; |
|||
|
|
|||||
h |
è |
2 |
ø |
|||
п |
|
|
|
|
|
|
