Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Оптимальный прием и обработка сигналов. Учебное пособие
.pdf
8.5. Потенциальная точность оценок параметров сигнала
З
Значит, потенциальная точность оценки задержки сигнала
по фазе обратно пропорциональна отношению сигнал/шум и
квадрату несущей частоты сигнала. Обычно
нозначном измерении фазы дисперсия ошибок такой оценки
меньше, чем у оценки по огибающей.
Пример 8.11. Найдем потенциальную точность оценки за-
держки
информации и по фазе, и по огибающей повышает точность
измерения задержки по сравнению с достигнутой в примерах
8.9 и 8.10. Следуя методике примеров 8.10 и 8.11, искомую
потенциальную точность выразим как
Пример 8.12. Найдем потенциальную точность оценки не-
энергетического параметра сигнала — частоты. Подставим
выражение (8.6) в выражение для
дифференцирование:
где
длительность сигнала. Видно, что потенциальная точность
оценки частоты обратно пропорциональна отношению сиг-
нал/шум и
Для сигнала
λ неравенство Рао — Крамера (8.4) обобщается на систему
неравенств
смещенной оценки
радиоимпульса по огибающей и фазе сигнала. Учет
22
Dq
DM
эфф
NtStdtNE q
00
TT
22 2
tS t dt S t dt
00
2,
NMytSt dt
0
T
22 2 2
2, 212,
0
,,
2
.
,Stλ с неслучайным векторным параметром
1
,
RJ где R — матрица корреляции ошибок не-
λ
12 .
Зэфф
T
0
с элементами
D
при
эфф
22
ln
22
— среднеквадратичная
RM
ij i i j j
, и при од-
1
и проведем
141

8. Измерение параметров сигналов
J
J
J
при ,1,ij m ;
J — информационная матрица Фишера
ij
с элементами
JM M
ij
1
Если
RJ, оценки
TT T
ln ln ln
pY pY pY
λλ λ
00 0
ij ij
при 1,im — совместно эффективны.
i
2
.
Пример 8.13. При совместно эффективной оценке двух па-
раметров (
) корреляционная матрица R ошибок имеет
2m
вид
1
DRJ H M pY
11 22 1 2 0 1
1
DRJ H M pY
22 11 1 2 0 2
2
DRRJ r DD
12 21 12 1 2
12 12
где
det
det , ln ;
J λ
det , ln ;
J λ
det , ,
2
JJJ ;
11 22 12
взаимной корреляции оценок
J
1 2 12 11 22
1
Сравним найденные выражения для
,rJJJ
и
;
2
T
T
— функция
,1 ,Hr
12 12
R с
из формулы
(8.4). В правых частях первых двух формул для
2
1
2
2
R второй
.
сомножитель дает средний квадрат ошибки эффективной
оценки одного параметра. Если ошибки оценок параметров
,0r
коррелированны (
12
), средний квадрат ошибки со-
вместно эффективных оценок больше, чем для одного параметра.
Пример 8.14. Найдем потенциальную точность совместной
оценки частоты и задержки сигнала равенства (8.6) по огибающей. Выразим элементы матрицы Фишера
142

8.5. Потенциальная точность оценок параметров сигнала
J
З
D
JM pY q
11 0
JM pY q
22 0
12 21 0
TT
t tSt St dt S t dt
где
00
222
222
JM pY qt
,, ,, ,,
ЗЗЗ З
T
ln , 2 ;
ln , 2 ;
ЗЗ
T
З
2
T
ln , 2 ,
ЗЗ
2
, от-
куда детерминант матрицы Фишера
det 4JJ J q t
J
11 22 12
2222
2
,
а элементы корреляционной матрицы ошибок совместной
оценки задержки и частоты выражаются как
DJ q t
З
det 2 ;
22
J
DJ q t
det 2 ;
J
11
DJ t q t
З
det 2 .
J
12
222
222
Для сигнала заданной формы с ростом
Тогда с ростом
дает D
(растет точность оценки частоты), и наоборот. Мерой
(падением точности оценки задержки) па-
22
2
2
2
падает [18].
точности совместной оценки частоты и задержки служит
1ЗDD
. Точность растет с увеличением отношения сиг-
нал/шум и базы
сигнала, а также при падении t (степени
корреляции между ошибками оценки задержки и частоты).
Пример 8.15. Найдем потенциальную точность оценки уг-
ловой координаты объекта фазовым методом (рис. 8.9). Пусть
плоская волна длины
A
—
A
и
1
. Расстояние между ними — d , а — угол фронта
2
падает на антенну с двумя элементами
143

8. Измерение параметров сигналов
E
волны с осью антенны. Обозначим
блюдение смеси сигнала и БГШ на
при
лезного сигнала на каждом элементе антенны
равна
подобия равно
и вектор полезного сигнала
тенциальная точность оценки угловой координаты фазовым методом равна
1, 2i . Шумы на
12
,cos 2sinSt A t d
(рис. 8.9). Энергия по-
20
Для векторного наблюдения нижняя граница Рао — Крамера
DMY
эфф
DMNStytdt
эфф
0Mntn t
A
1
Рис. 8.9. Схема фазового метода измерения
угловой координаты объекта
1ln
YNtttdt
A и 2A статистически независимы:
1
. Если
22
T
exp 2 , 0,5
00
22
2,
12 2 cos .qd
10
Плоская волна
d
T
0
,,
tStSt
S
T
02 2
0
tytyty
A и 2A :
1
, тогда
T
0
12
12
yt St nt
iii
cosSt A t
A
2
, где отношение правдо-
Sy S ,
2
T
как на-
2
2
AT .
T
. Отсюда по-
1
Значит, точность оценивания растет вместе с величиной
cosd
, где cos 1d
рения фазы.
144
требуется для однозначности изме-

8.5. Потенциальная точность оценок параметров сигнала
Пример 8.16. Найдем потенциальную точность оценки уг-
ловой координаты объекта амплитудным моноимпульсным
методом, схема которого показана на рисунке 8.10. Пусть
— направление на максимум диаграммы направленности
F
i-й антенны при 1, 2i . Диаграмма направленности
i
симметрична:
FF . Равносигнальное направление
i
определяется как
Рис. 8.10. Схема амплитудного метода измерения
Сигналы источника на входах антенн равны:
Их энергии
8.15 найдем потенциальную точность оценки угловой координаты амплитудным методом:
DMN StytdtE
эфф
St AF t
ii
2
AT N Y F Y
.
012
F
1
0
угловой координаты объекта
,cos
EATF
ii
22
2,0,5
2
0
i
12 ,
при 1, 2i .
22
. По методике примера
2
0
1
i
1
qF
2
*
Источник
F
2
2
1
0
2
ii i
1
T
ii i
0
2
2
i
1
i
1
2
1
145

8. Измерение параметров сигналов
6
p
E
где
YFF
.
iii
Если излучатель находится на равносигнальном направле-
нии
(максимум точности оценки
, достигается минимум дисперсии ошибок оценки
0
22
), равный
Чем выше скорость изменения функции
(больше
Пусть при
шагом дискретизации
nk — дискретный БГШ. Теория статистических решений
(см. п. 1.5) верна как для непрерывных сигналов, так и для
дискретных. Остаются верными и основные результаты пунктов 8.1—8.5 с заменой формы записи соответствующих величин. Так, условная ПВ
pY yi S i
При оценке непрерывного параметра λ дискретного сигнала оптимальное байесовское решение для простой функции
потерь найдем из формулы (8.1). С учетом формулы Байеса
DqF
F
), тем меньше ошибка оценки
8.6. Оценка параметров дискретного сигнала
T
0
2exp12 ,.
λλ
эфф
min
1
0,tT наблюдается выборка
K
22
TT T
YpYp pYλλλ
00 0
14
0
:
t
2
1
ap
2
.
0
F
в точке
1
T
Y отсчетов с
0
,yk Sk nkλ , 1,kK , где
i
1
.
2
0
имеем
K
22
ln , 2 0,
pyiSiE
λλ λ λ λ
ap
где
14
K
λλ
i
1
i
2
,
Si
1
— мощность отсчетов сигнала.
λλ
ˆ

8.6. Оценка параметров дискретного сигнала
E
E
Применяя метод максимального правдоподобия (см. формулу (8.2)), вместо формулы (8.1) найдем уравнение
K
22
yiSi E
i
1
,2 0
λλ λ
λλ
ˆ
.
Пример 8.17. Возьмем дискретный сигнал
при
этого сигнала
1,iK , где
Si fit it
. Мощность
10
2
AAE
, где
1
cos
З
11
K
2
Si
i
1
,SiA AS i
1
. Оценку макси-
мального правдоподобия амплитуды A найдем из уравнения
yiSiA E A A
i
1
,2 0
AA
ˆ
K
22
как
ˆ
AE yiSi
1
K
11
i
1
.
Пример 8.18. Оценку максимального правдоподобия на-
чальной фазы
где
,,Si AS i
010
,cos
Si fit it
10 0
найдем из уравнения
0
K
22
yiSi E
i
1
при 1,iK ;
;
,2 0,
З
00 0
EAE
010
ESi
10 1 0
00
2
K
i
ˆ
;
2
,
1
.
Решение этого уравнения имеет вид
K
yi f i t i t
1
ˆ
arctg
0
i
K
yi f i t i t
1
i
sin
З
cos
З
.
8.7. Оценка случайных параметров сигнала
Пусть
метр
,St — сигнал, случайный информативный пара-
которого распределен с заданной априорной ПВ
147

8. Измерение параметров сигналов
M
ˆ
p
. Выразим смещение
ap
ˆ
параметра . Так как — случайная величина, то и
M
некоторой оценки
—
случайная величина. Значит, оценка случайного параметра —
0M
несмещенная в среднем:
метр
сигнала
,St по результатам наблюдений
. Оценим случайный пара-
T
Y . При
0
незначительных ограничениях средний квадрат ошибки любой
ˆ
оценки
не может быть меньше нижней границы, определяе-
мой неравенством Рао — Крамера [3]:
1
где
*,
для которой
DM pY
D достигает нижней границы Рао — Крамера.
MpY D
TT
ddYp Y
00
ˆ
эфф
ln ,
22
ln , ,
T
0
T
0
;
эфф
2
1
ˆ
эфф
— эффективная оценка,
Известно, что для эффективной оценки случайного параметра апостериорная ПВ должна иметь нормальное распреде-
ление
3. Так как
TT
,
pY pY p
, то
00
ap
2
ln ,
MMM
T
pY
0
222
2
ln
pY
T
0
2
ln
p
ap
.
Здесь первое слагаемое правой части учитывает информацию, найденную обработкой результатов наблюдений. Если
параметр
— детерминированный, то это слагаемое равно J
— информации по Фишеру (см. формулу (8.4)). Второе сла-
гаемое учитывает априорную информацию. Если сигнал
,Stλ имеет случайный векторный параметр λ , неравенство
Рао — Крамера обобщается на систему неравенств
R — корреляционная матрица средних квадратов ошибок
где
148
1
RJ,

8.7. Оценка случайных параметров сигнала
J
И
X
x
M
с элементами
J — матрица с элементами
ij
JM M
Если
тивны.
Пусть сигнал
мативный параметр
метр
Полагаем, что априорная ПВ
оценки параметра
отношение правдоподобия для сигнала со случайным параметром, применим методику, приведенную в пункте 2.3.
Пример 8.19. Пусть
. Отношение правдоподобия для радиоимпульса дли-
a
тельностью
(см. п. 2.3.2)
ij
RJ, то оценки
:
,, cosSt aAt t t
YaNIaNXT
RM
ij i i j j
ln , ln , ln ,
1
TT T
pY pY pY
λλ λ
00 0
ij ij
,,St
применим формулу (8.3). Чтобы выразить
, усредненное по начальной фазе, выразим как
T
exp 2 2 .
0000
ˆˆ
ˆ
при 1,im — совместно эффек-
i
имеет детерминированный инфор-
и случайный неинформативный пара-
p
ap
,
p
0
ap
2
И
при
2
при
0
задана. Для поиска
12
,
0
0,tT .
,
0
,1,ij m ;
.
,
Подставив это выражение в формулу (8.3) и дифференци-
руя его, преобразуем уравнение правдоподобия:
ˆ
2
IaXTN
10
где учтено тождество
Полученное для
При больших значениях отношения сигнал/шум, когда
ˆ
a
M
M
ˆ
IaXTN
2
00
Ix dI xdx и
10
ˆ
a уравнение можно решить численно.
M
2
T
,
И
ˆ
2
aXT N .
M
0
1x ,
149

8. Измерение параметров сигналов
И
M
И
З
З
З
X
X
найдем приближенно
ˆ
2
aXT
. Эта оценка
M
ˆ
a ампли-
туды радиоимпульса со случайной начальной фазой пропорциональна значению отсчета огибающей на выходе квадра-
tT
турного приемника (см. рис. 2.9) в момент
. Потенциаль-
ная точность оценки амплитуды определяется выражением
1
22
DM aa
a
эфф
ln ,
которое при большом отношении сигнал/шум принимает вид
DN
a
эфф
.
0
Точность оценки амплитуды сигнала со случайной начальной фазой при больших значениях отношения сигнал/шум
приближенно равна соответствующей точности для сигнала с
известной начальной фазой. При малых значениях отношения
сигнал/шум первая из этих оценок менее точна, чем вторая.
Пример 8.20. Решим задачу оценки задержки
огибаю-
щей сигнала
где
St aAt t
;
p
0
ap
0
12
;
З
,
;
0
0
.
,, cos ,
Отношение правдоподобия получено в пункте 2.3.2, где от
зависит лишь огибающая:
,,,,
Xt X t X t
З C З S З
tyA d
,cos,
C ЗЗ
tyA d
,sin.
S ЗЗ
22
t
0
t
0
Тогда из уравнения правдоподобия
,0
XT
ЗЗ
ЗЗ
ˆ
найдем
150
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]
