Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Оптимальный прием и обработка сигналов. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
1 Мб
Скачать
☆
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
JJJ ;
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
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]