Статистическая теория радиотехнических систем. Учебное пособие
.pdf
Используя обозначения
|
|
|
Ex x t 2dt , Z |
x t s t, dt, |
(3.57) |
|
|
|
|
|
|
E s2 t, dt,
алгоритм максимального правдоподобия можно записать в виде
ˆ argmax |
|
Ex |
|
|
|
2 |
|
Z 0 |
|
E 0 |
|
|
|||
|
|
|
|
|
|
||||||||||
N0 |
N0 |
|
|
||||||||||||
0 |
|
|
|
|
E 0 |
N0 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
||||||
argmax |
2 |
|
Z 0 |
|
. |
(3.58) |
|||||||||
|
|
|
|||||||||||||
0 |
|
N0 |
|
|
N0 |
|
|||||||||
Параметр называется энергетическим, если энергия сигнала является функцией информационного параметра , в остальных случаях параметр называют неэнергетическим.
Энергетическими параметрами являются амплитуда сигнала, длительность сигнала, полоса частот. Неэнергетические – это средняя частота сигнала, начальная фаза сигнала, временная задержка.
Для неэнергетического параметра алгоритм оценки максимального правдоподобия сводится в соответствии с (3.58) к максимизации значения корреляционного интеграла, т.е.
ˆ argmaxZ 0 . |
(3.59) |
0 |
|
Рис.3.4. Структурная схема устройства оценки неэнергетического параметра сигнала на фоне белого шума.
На Рис.3.4 показана структурная схема устройства оценки неэнергетического параметра сигнала на фоне белого шума.
Решающее устройство на схеме вычисляет опорное значение параметра для работы коррелятора и принимает решение о значении опорного параметра, при котором достигается максимум корреляционного интеграла.
51
Рассмотрим случай наличия мешающего параметра. Пусть мешающий параметр - начальная фаза сигнала, тогда
x t s t, , n t , |
(3.60) |
где s t, , Re S t, ej 0t , S t, - комплексная огибающая сигнала s t, , , - случайная фаза, равномерно распределенная на интервале
; .
Тогда функционал правдоподобия
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E |
x |
|
|
|
|
2 |
|
|
|
E |
||||
p x |
|
, c |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z , |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
exp |
|
Ex 2Z , E |
|
c e |
|
N0 eN0 |
|
|
|
2 .(3.61 |
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
N0 |
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
) |
|||||
|
|
Интегрируя по мешающему параметру, получим |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ex |
|
|
E |
|
|
|
1 |
|
2 |
Z , |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
p x |
|
c e N0 e |
|
N0 |
|
|
eN0 |
|
|
|
d |
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
Ex E |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
Z cos argZ |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
c e |
|
|
|
N0 |
1 |
eN0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d . |
|
(3.62) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Алгоритм оценивания в случае случайной начальной фазы имеет вид |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
ˆ argmax p x |
|
|
argmax I |
|
2 |
|
Z |
|
|
|
|
|
|
|
E |
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
exp |
|
|
, |
(3.63) |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
0 |
|
|
N |
0 |
|
|
|
|
|
|
|
|
N |
0 |
|
|
|
|||||||||||||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где Z |
X* t S t, dt, |
X t |
- комплексная огибающая наблюдаемого |
||||||||||||||||||||||||||||||||||||||||||||||||
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
сигнала. Для случая неэнергетического информационного параметра алгоритм приема сводится к виду
ˆ argmax |
|
Z 0 |
|
. |
(3.64) |
|
|
||||
0 |
|
|
|
|
|
Для вычисления комплексного корреляционного интеграла требуется квадратурный приемник.
52
Рис.3.5. Структурная схема устройства оценки неэнергетического параметра сигнала со случайной фазой на фоне белого шума.
Рассмотрим теперь способы вычисления погрешности алгоритмов оценки параметров сигналов, наблюдаемых на фоне белого гауссовского шума. В СТРТС обычно качество оценивания определяют, используя границу Рао – Крамера. Получаемые таким образом значения погрешности для оценок максимального правдоподобия являются потенциальными (асимптотическими), т.е. реальная погрешность алгоритма оценивания может быть только хуже. Для эффективных оценок погрешности, вычисленные по границе Рао – Крамера совпадают с реальной погрешностью.
Рассмотрим случай оценки векторного параметра λ 1 n T сигнала s t,λ , наблюдаемого на фоне аддитивного белого гауссовского шума. В этом случае наблюдаемый сигнал имеет вид
|
|
|
|
|
|
|
x t s t,λ n t . |
(3.65) |
||||||
Функционал правдоподобия можно записать в виде |
||||||||||||||
|
||||||||||||||
|
|
E |
x |
|
2Z λ |
E λ |
|
|||||||
|
|
|
||||||||||||
p x |
λ c exp |
|
|
|
|
|
|
|
. |
(3.66) |
||||
|
|
|
|
|
|
|||||||||
|
|
N0 |
N0 |
N0 |
|
|||||||||
|
|
|||||||||||||
Алгоритм вычисления МП оценки имеет вид |
|
|||||||||||||
λˆ argmax p x |
|
λ argmaxln p x |
|
λ . |
(3.67) |
|||||||||
|
|
|||||||||||||
|
λ |
|
|
|
|
|
λ |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||||
Вычислим значения информационной матрицы Фишера, подставляя (3.66) в (3.37) обозначая λ0 - истинное значение вектора параметров
53
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2Z λ |
E λ |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
λ |
|
М |
|
|
|
|
|
|
|
|
ln p x |
λ |
М |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(3.68) |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
ij |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
N |
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
В случае неэнергетических параметров, у которых E λ E, получим |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
|
М |
|
|
|
|
|
|
|
|
|
|
|
|
Z λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
М Z λ |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
ij |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
0 |
|
|
j |
|
|
|
|
N |
0 |
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
М x t s t,λ dt |
|
|
|
|
|
|
|
s t,λ0 s t,λ dt |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 i j |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
0 i j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
λ0 |
||||||||||||||||||||||||||||||||
|
|
|
|
2E |
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
s t,λ0 s t,λ dt |
|
|
|
q2 |
|
|
|
|
|
|
λ,λ0 |
|
λ λ0 |
(3.69) |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
N0 i j |
E |
|
|
|
|
i j |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ λ0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
где q |
|
|
|
|
2E |
|
, λ,λ0 - функция неопределенности сигнала s t,λ . |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ,λ0 |
s t,λ0 s t,λ dt, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.70) |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Часто бывает, что функция неопределенности λ,λ0 зависит только
от разности λ λ0 , в этом случае λ0 q2 0 .
В случае оценивания одного параметра погрешность оценки определяется в виде
D ˆ |
1 |
1 |
|
|
|
|
|
|
. |
(3.71) |
|
λ0 |
q2 0 |
||||
Функция неопределенности определяется только свойствами сигнала и не зависит от шума и его характеристик. Поэтому функция неопределенности используется для классификации, изучения свойств сигналов, используемых для построения РТС.
В соответствии с (3.71) нижняя граница потенциальной погрешности алгоритма оценивания определяется двумя факторами отношением сигнал/шум и степенью остроты пика функции неопределенности в нуле.
Если анализировать реальную погрешность устройств, схемы которых показаны на Рис.3.5, Рис.3.4 и алгоритмов оценки (3.59) и (3.64) соответственно следует учитывать, что помимо ошибок, связанных с влиянием шума на оценку в области истинного значения параметра, на практике возможны ошибки, вызванные появлением на входе решающего устройства нескольких локальных максимумов корреляционного интеграла, вызванных выбросами шума. При этом решающее устройство будет совершать так называемые аномальные ошибки. Т.е. решающее устройство будет принимать решения в пользу совершенно далекого от истинного значения измеряемой величины.
54
Эта ситуация возможна только при относительно низком отношении сигнал-помеха. В этом случае формулы, основанные на границе РаоКрамера, будут существенно занижать истинное значение погрешности.
Т.о. формула (3.71), как и другие формулы, основанные на границе Рао-Крамера, характеризуют погрешность МП оценок параметров сигналов на фоне аддитивных шумов только при достаточно высоких отношениях сигнал-шум или при наличии в решающем устройстве достаточно подробной информации об измеряемой величине, т.е. когда вариация измеряемого параметра в устройстве обработки не выходит за ширину главного лепестка функции неопределенности сигнала. Как правило, для реальных РТС эти условия выполняются.
3.2.1. Оценка амплитуды детерминированного сигнала
Рассмотрим пример оценки амплитуды детерминированного сигнала s t . В этом случае наблюдаемый сигнал
|
|
|
|
|
|
|
|
x t a s t n t , |
|
|
|
|
(3.72) |
|||||||||||
|
|
Энергия |
сигнала E a a2Es , |
корреляционный интеграл Z a aZ , |
||||||||||||||||||||
тогда оценка МП имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
aˆ argmax |
|
exp |
|
2aZ a |
2 |
Es |
|
|
(3.73) |
||||||||||
|
|
|
|
|
|
|
|
|
. |
|||||||||||||||
|
|
|
|
|
|
N0 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Решая |
уравнение для |
|
максимума |
|
функции |
правдоподобия |
||||||||||||||||
|
ln p x |
|
a 0, получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
a |
|
|
|
|
|
|
|
2Z |
2a |
Es |
0, |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|||||||
отсюда |
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Z |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
aˆ |
|
|
|
|
x t s t dt . |
|
|
|
(3.74) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
E |
s |
|
|
|
E |
s |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Рис.3.6. Структурная схема устройства оценки амплитуда сигнала на фоне белого шума.
Дисперсию оценки легко вычислить в виде |
|
|
|
|
|
|
|||||||
D aˆ |
|
a |
1 |
M Z M Z 2 |
1 |
|
N0Es |
|
N0 |
|
1 |
. |
(3.75) |
|
|||||||||||||
|
|
|
|
|
|
||||||||
|
|
|
Es |
Es2 |
2 |
|
2Es |
|
q2 |
|
|||
|
|
|
|
||||||||||
55
ЗАКЛЮЧЕНИЕ
Теория несмещенных оценок предполагает использование критериев несмещенности и эффективности для вычисляемых значений измеряемых параметров. Граница Рао-Крамера позволяет оценить потенциальную погрешность несмещенных оценок.
КОНТРОЛЬНЫЕ ВОПРОСЫ
1.Какие критерии для отбора оптимальных оценок используются в теории точечного оценивания?
2.Что такое граница Рао-Крамера?
3.Как с помощью информационной матрицы Фишера определить границу для погрешности одновременной оценки нескольких параметров сигнала?
4.Какими асимптотическими свойствами обладают оценки максимального правдоподобия?
5.Какие параметры сигнала называют энергетическими, а какие неэнергетическими?
6.Какие параметры сигнала называют мешающими?
7.Каков алгоритм оценки МП неэнергетического параметра сигнала, наблюдаемого на фоне БГШ.
8.Какова структурная схема алгоритма оценки МП неэнергетического параметра сигнала, наблюдаемого на фоне БГШ и имеющего неопределенную начальную фазу.
9.Как в СТРТС рассчитывают погрешность оценивания параметров?
10.Чем определяется потенциальная погрешность оценки параметров сигнала, наблюдаемого на фоне аддитивного БГШ?
11.Что такое функция неопределенности сигнала по параметру?
12.Чем определяется потенциальная погрешность оценки параметров сигнала со случайной фазой, наблюдаемого на фоне аддитивного БГШ?
13.Каков оптимальный алгоритм совместной оценки амплитуды и фазы сигнала на фоне БГШ.
14.Что такое аномальные ошибки в теории оценивания?
56
ЛЕКЦИЯ №8. Оценка комплексной амплитуды сигнала. Оценка времени запаздывания сигнала. Оценка частоты сигнала со случайной фазой. Совместная оценка времени запаздывания и частоты сигнала со случайной фазой. Расчет погрешности оценок.
3.2.2. Оценки амплитуды и фазы детерминированного сигнала
Рассмотрим пример совместной оценки амплитуды и фазы детерминированного сигнала s t . В этом случае наблюдаемый сигнал имеет вид
x t s t,a, n t , |
(3.76) |
где s t,a, Re aS t ej 0t , S t - комплексная огибающая сигнала, - начальная фаза, a - амплитуда сигнала.
Функционал правдоподобия имеет вид
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p x |
a, c exp |
|
Ex |
2aZ a2Es . |
|
|
|
|
|
|
|
(3.77) |
||||||||||||||||||
|
N0 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos arg Z , |
|||
Используя |
|
|
ранее |
полученное |
|
соотношение |
Z |
|
Z |
|
||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||
получим алгоритм оценивания в виде |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
aˆ |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
argmax c exp |
|
|
Ex 2aZ a2Es |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
ˆ |
|
|
|
|
|
a, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
Ex |
|
2a |
|
Z |
|
cos arg Z |
a |
|
|
|
Es |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
argmax exp |
|
|
|
|
2 |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.78) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
argmax 2acos arg Z a2Es |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
a, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
Z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Es |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
arg Z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
57
Рис.3.7. Структурная схема устройства совместной оценки амплитуды и фазы сигнала на фоне белого шума.
Определим погрешность совместного оценивания амплитуды и фазы сигнала. В соответствии с (3.78) запишем
|
|
|
ln p x |
|
a, ln c |
Ex |
|
|
|
|
2a |
|
|
|
Z |
|
|
cos arg Z |
a2Es |
. |
(3.79) |
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
Вычислим коэффициенты информационной матрицы Фишера |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
2Es |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
М |
|
|
|
|
|
|
|
|
|
|
ln p x |
a, |
|
|
q2 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
a a |
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
М |
|
|
|
|
|
|
|
|
|
|
|
|
|
ln p x |
a, |
|
|
|
|
|
|
М |
|
|
|
|
|
|
Z |
|
|
cos arg Z |
|
|
|||||||||||||||||||
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
2 |
М Z |
|
sin arg Z |
|
|
2 |
М Im Zej |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
2 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
М X* t S t dt ej |
|
||||||||||||||||
|
|
|
М Im |
|
|
|
|
|
|
X* t S t dt ej |
|
|
|
|
|
|
|
Im |
|
|
|
||||||||||||||||||||||||||||||||||
|
|
2 |
|
|
|
|
N0 |
2 |
|||||||||||||||||||||||||||||||||||||||||||||||
|
N0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
2a0 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2a0 |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
Im |
|
|
|
S t |
2dt ej 0 |
|
|
|
|
|
|
|
Im Es 0, |
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
|
N0 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 |
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
21 |
12 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
58
|
|
|
|
2 |
|
|
|
|
|
|
|
2a |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
22 М |
|
ln p x |
a, |
|
|
|
М |
|
|
Z |
sin arg Z |
|
|||||||||||||
|
|
N0 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|||
|
2a |
М Re Zej |
2aa0 |
Re |
1 |
|
|
S t |
2dt ej 0 |
|
|
2a0 |
Es |
a02q2, |
|||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
N |
0 |
|
|
|
N |
0 |
|
2 |
|
|
|
|
|
|
|
|
0 |
|
N |
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a a0 |
|
|
|
|
где a0, 0 - истинные значения амплитуды и фазы.
Таким образом, информационная матрица Фишера и ей обратная матрица имеют вид
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
|
|
|
|
|
|
|
q2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Φ a |
|
, |
|
0 |
|
|
|
|
Φ |
|
|
, |
|
|
q2 |
|||||||||||
0 |
0 |
|
|
2 |
2 |
, |
|
|
1 a |
0 |
0 |
|
|
1 |
|
. |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
0 |
|
a0q |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a0q |
|
|
Погрешность оценки амплитуды и фазы можно вычислить по |
(3.80) |
|||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||
формулам |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D aˆ |
|
a0, 0 |
|
, |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
q2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
D ˆ |
|
a0, 0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
(3.81) |
|||||
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
a02q2 |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
3.2.3. Оценка времени запаздывания детерминированного сигнала
Рассмотрим случай оценки времени запаздывания детерминированного сигнала s t . В этом случае наблюдаемый сигнал имеет вид
x t s t n t , |
(3.82) |
где - неизвестное время запаздывание сигнала. |
Очевидно, что время |
запаздывания является неэнергетическим параметром. В этом случае мы можем воспользоваться алгоритмом (3.59), т.е.
|
|
|
s t x t dt. |
|
|
(3.83) |
|||
ˆ argmaxZ argmax |
|
|||
Рис.3.8. Структурная схема устройства оценки временной задержки сигнала на фоне белого шума.
Определим погрешность оценки задержки, используя формулу (3.71)
59
D ˆ |
|
1 |
. |
2 |
|
||
|
0 |
||
|
q |
||
|
|
|
|
Вычислим вторую производную функции неопределенности (3.70), воспользовавшись спектральными представлениями сигналов
Ф
s(t) s( j ),
Ф |
j |
, |
s(t ) s j e |
|
1
s t s t dt.
E
Воспользуемся равенством Рэлея для преобразования Фурье, тогда
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
s* j ej s( j )d |
|
|
|
|
|
|
s( j ) |
|
2ej d . |
(3.84) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
2 E |
|
|
2 E |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Свойство дифференцирования оригинала дает возможность |
|||||||||||||||||||||||||||||||||||
определить производные функции неопределенности в виде |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
2 |
|
|
|
j |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
j |
s( j ) |
e |
d , |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
2 E |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
2 |
|
|
|
j |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
s( j ) |
|
|
e |
|
|
|
|
|
|
|
d . |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
2 E |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2d |
|
||||||||||
|
|
|
1 |
|
|
|
|
|
2 |
|
|
2 |
|
s( j ) |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
s( j ) |
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
эф , |
(3.85) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
2 E |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s( j ) |
|
d |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s t - |
|||||
где эф - |
|
эффективная полоса частот сигнала. Если сигнал |
|||||||||||||||||||||||||||||||||
узкополосный, |
как, |
|
например, |
|
радиосигнал, |
|
|
|
то |
|
|
эф 0 - несущая |
|||||||||||||||||||||||
частота радиосигнала. Если s t - низкочастотный сигнал, например видеосигнал, то эф - полоса частот видеосигнала.
Погрешность оценки времени запаздывания детерминированного сигнала имеет вид
1 |
|
D ˆ q2 эф 2 . |
(3.86) |
3.2.4. Оценка времени запаздывания сигнала со случайной начальной фазой
Рассмотрим задачу оценки времени запаздывания сигнала со случайной
60
