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

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

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
1 Мб
Скачать
☆
8.3. Оценка непрерывных параметров сигнала

И
E
M
И
F N St yt St dt

   
Найдем оценки параметров радиоимпульсов, решая урав­нение (8.5). Пусть принимаемый сигнал
Пример 8.1. Оценим амплитуду A радиоимпульса
где
St ft t
нальная функция;
0t или
Пусть все параметры сигнала
A , известны. Амплитуда A — энергетический параметр. Надо
учесть зависимость
мулу (8.5) при A ное уравнение, найдем оценку амплитуды
где
EStdt
измеритель амплитуды — коррелятор входного колебания с
опорным
вание БГШ, поэтому оценка ˆA распределена нормально. Эта оценка — несмещенная, так как
Дисперсия оценки амплитуды равна дисперсии ошибки оценки амплитуды, пропорциональна интенсивности БГШ и обратно пропорциональна длительности импульса:
2,0,5, 0.

 ;
 
10
t
;
T

11
0
St E. Выражение для ˆA — линейное преобразо-
11

AM E AStntStdtA


ˆˆ ˆ
MAMA MAA N E N




T
0
0

З
И
, дифференцируя по A и решая получен-
2


,St AS t
cos

T
;
A . Подставляя выражение (8.6) в фор-
ˆ
AEytStdt
1
 
2
И
при
1
St — нормированная сиг-
1ft, если 0Иt
k

И
T
11
0
. Отсюда следует, что оптимальный
1
; 1,2,...k
1

,yt St nt
.
0,tT , (8.6)
;
,St, кроме амплитуды
,
ˆ
M
0ft, если
(8.5)

ˆ


1.
 


2
 
T

11 1
0




2

2.

01 0
131
8. Измерение параметров сигналов
З
Пример 8.2. Проведем оценку неэнергетического парамет-
ра — начальной фазы
— радиоимпульса
0
,St — при
1
условии, что остальные параметры радиоимпульса известны.
Подставив выражение (8.6) в формулу (8.5) при
ференцируя по
и решая полученное уравнение, найдем
0
, диф-
0
оценку начальной фазы:
ˆ
arctg y t f t tdt y t f t tdt


0
TT
   

00
sin cos .
ЗЗ
На рисунке 8.3 показана схема измерителя начальной фазы сигнала.
yt

T
0
T
0
Рис. 8.3. Схема измерителя начальной фазы сигнала
u
1
arctg u u
u
2
12
0
Если параметр сигнала оценить аналитически трудно, применяют поисковые методы решения уравнения правдопо­добия.
Пример 8.3. Оценим временное положение сигнала
,Stпо неэнергетическому параметру — запаздыванию
1
огибающей — с учетом того, что остальные параметры
сигнала известны. Подставив выражение (8.6) в формулу (8.5), найдем уравнение правдоподобия:
T
yt f t t dt


0
 
 
ЗЗ
cos 0.

ˆ

ЗЗ
0
Значит, оптимальный измеритель задержки (рис. 8.4) дела­ет корреляционную обработку принятого сигнала
yt с про-
132
8.3. Оценка непрерывных параметров сигнала
З
З
З
З
З
изводной огибающей полезного сигнала
,Stпо задержке.
1
Производная служит опорным сигналом. Решение уравнения
ˆ
правдоподобия — оценку
задержки — ищем методом по-
следовательных приближений.
yt

ft

 
ЗЗ
Рис. 8.4. Схема последовательного поиска
cos t
ˆ
З
оценки задержки сигнала
T
0
Схема фиксации перехода через 0
0
ˆ
В схеме, изображенной на рисунке 8.4, запоминается реа­лизация принятого сигнала. Для текущего приближения вели­чины
ˆ
выход интегратора дает корреляционный интеграл,
значение которого запоминается. По таким значениям, взятым из памяти устройства, ищется следующее приближение для
ˆ
. Операции повторяются, пока значение отсчета на выходе
интегратора не перейдет через ноль. Недостатки последова­тельного поиска — необходимость запоминать принятую реа-
лизацию и медленная сходимость. Параллельный поиск [20] основан на дискретизации области возможных изменений оцениваемого параметра. После операции поиска применяют методы оценки, указанные в пункте 8.2 (на основе многока­нальных приемников), как системы
параллельного поиска по
этому параметру.
Пример 8.4. Оценим частоту
радиоимпульса
1
,St,
считая остальные его параметры известными. Запишем урав­нение правдоподобия
T
ty t f t t dt
 
0
sin 0.

З

0
(8.7)
133
8. Измерение параметров сигналов
З
Процедура последовательного поиска по частоте реализо­вана в схеме оптимального измерителя частоты на рисунке 8.5.
yt

cos t
Обозначим уравнение (8.7) по частям и аппроксимируем производную под знаком интеграла:
TT
ytf t t dt dtytf t
 

00



На каждой итерации перестраиваются частоты двух опор-
ных сигналов —
тем находится и запоминается корреляционный интеграл, по значениям которого из памяти устройства, находится следую­щее приближение для отсчет на выходе сумматора не перейдет через ноль.
Решим уравнение правдоподобия (8.5) методом Ньютона
[10]:
ft

ˆ


Рис. 8.5. Схема последовательного поиска

  
З З
ˆˆ
cos cos 2 0.tt
 

ˆ

cos t


0
оценки частоты сигнала
как расстройку по частоте. Интегрируем
cos

00
ˆ
cos t
ˆ
ii i
  

 ,
M ОП ДН ОП
T
0
0
T
0
0




0
. Итерации следуют до тех пор, пока
u
Схема фиксации перехода через 0
ˆ
и
cos t
 

ˆ


0
ˆ
, за-
134
8.3. Оценка непрерывных параметров сигнала
M
M


M
M
i
где i — номер итерации; 0,1,...i
i

— i -е приближение к оценке

u

ДН ОП
lnF
 
  (см. формулу (8.3)), то
i

 
F
F




i
Устройство, измеряющее рассогласование
i
текущим приближением
i

ем
, — дискриминатор. Величина
ОП
к оценке
;
— опорное значение;
ОП
ˆ
,
ОП M
ln

 
22
ln
ˆ
u
ДН ОП
1ii
. Так как

ii
M ОП
.
i


ОП
между
и опорным значени-
i

— нормиро-
ванная характеристика дискриминатора. В результате опреде­лен метод нахождения оценки параметра сигнала с помощью дискриминатора.
8.4. Следящие измерители параметров сигнала
Измерить параметр
сигнала сложнее, если его значение
изменяется со временем. Дискретизируем сигнал, принятый на отрезке
достаточно медленно — с характерным временем На каждом интервале t
эффективное значение
ний
0,T , по времени с шагом

. Пустьизменяется
t
присвоим параметрунекоторое
. Для оценки эффективных значе-
ЭФФ
применим методы, приведенные в пунктах 8.2 и 8.3.
ЭФФ
  .
t
Точность оценки можно повысить, а время вычисления оценки — уменьшить. Для этого надо: 1) формировать текущую оцен-
ку
ченных ранее; 2) перейти от измерений значения дом интервале
с учетом как текущих данных о сигнале, так и полу-
ЭФФ
ЭФФ
к оценке его изменения при переходе от од-
t
на каж-
ного интервала к другому. Эта методика реализуется в следя­щих измерителях (рис. 8.6).
135
8. Измерение параметров сигналов
6
З

M
З
ii
M ОП
i
ОП
ˆ
Фильтр­экстраполятор
yt

Дискриминатор
Рис. 8.6. Обобщенная структурная схема следящего измерителя
Пример 8.5. На рисунке 8.7 представлена схема временно-
го дискриминатора, который можно применить в следящем
измерителе задержки сигнала. СФ и амплитудный детектор Д образуют приемник. Детектор с линейной характеристикой выражает производные в конечных разностях. На входы пере-
множителя приходят сигналы
St и
СТР
St с выхода гене-
АД
ратора опорного напряжения ОП и приемника соответственно. ОП создает стробирующий импульс
где
Stft ft t
1ft, если


СТР З З
02tt


, и

0,5 ,
0ft, если 0tили
2tt
ii
yt

СФ
Д
St
АД
ˆ
S
ОП
St
СТР
Рис. 8.7. Структурная схема временного дискриминатора
Строб состоит из двух прямоугольных импульсов проти­воположной полярности. Положение середины строба на оси
i
времени определяет момент
— текущее приближение к
оценке задержки, равное напряжению на входе ОП. Если из-
.
13
8.4. Следящие измерители параметров сигнала
З
З
П
П
П
П
П
П
вестно априорное значение
i
сто
. Перемножитель образует сигнал, пропорциональный
площади перекрытия сигналов
напряжение
S интегратора оценивает величину и знак

ˆ
задержки, его можно взять вме-
St и
АД
St. Выходное
СТР
ошибки совмещения середины строба с полезным сигналом.
Пример 8.6. На рисунке 8.8 представлена обобщенная
структурная схема частотного дискриминатора с амплитуд­ным детектированием, который можно применить в следящем измерителе частоты сигнала.
ˆ
Д
Управляемый гетеродин
yt

СМ
Ф
1
Ф
2
Д
Д
1
2
ˆ
ln 2


Д
ˆ
ln 2


Д
Рис. 8.8. Структурная схема частотного дискриминатора
с амплитудным детектированием
Принятое колебание
yt сдвигается на промежуточную
частоту в смесителе СМ, после чего поступает в два канала обработки. Полосовые фильтры
частоты
ˆ

Д
ˆ
и
2
Д

Ф и
1
соответственно, где
2
Ф настроены на
2
ˆ
—
Д
текущая оценка доплеровского сдвига частоты сигнала (или ее априорное значение), а
лах обработки формируются сигналы
ˆ
ln 2

 . Каждому из них соответствует напряжение,

Д
— расстройка по частоте. В кана-
ˆ

ln 2
Д

и
зависящее от степени перекрытия спектра сигнала, поступив­шего в канал, с АЧХ фильтра —
Ф или
1
Ф . Эти напряже-
2
137
8. Измерение параметров сигналов

A
ния одинаковы, если 0
ˆ
ln 2

 больше, чем сигнала

Д
. На выходе дискриминатора получим разность
0
ˆˆ
ln 2 ln 2
 
   .
8.5. Потенциальная точность оценок параметров сигнала
ДД
 . Напряжение сигнала
ˆ
ln 2

 , если
Д
Потенциальная точность — предел точности оценки, воз-
T
Y
можной при обработке наблюдаемой реализации
без учета
0
априорных данных. Потенциальная точность — эффективная оценка. Если параметр
сигнала скалярный и неслучайный,
потенциальная точность оценки его значения определяется
D
нижней границей
Граница
D
обратно пропорциональна информации J по
эфф
Фишеру. Значение
Рао — Крамера (см. формулу (8.4)).
эфф
D
дает дисперсию ошибок при достиже-
эфф
нии потенциальной точности оценивания. Конкретизируем эту границу. Имеем
1
22


M N St yt St dt


T
где
Y — отношение правдоподобия.
0
Пример 8.7. Найдем потенциальную точность оценки ам-

DMpY
эфф
22
2,0,5,
  
0

 
T

0
22

ln ,
MY

плитуды A радиоимпульса
(8.6) в выражение для
D
эфф
T
ln


0

T


0
,St. Подставив выражение
1
при
1
1
и проведя дифферен-
цирование, найдем
138
8.5. Потенциальная точность оценок параметров сигнала
И
D
З
З
E
З
З
З
1
1
22 2
DM AA NStdtN

A
эфф

ln 2 .




T
01 0
0
Видно, что дисперсии ошибки оценки амплитуды (см. пример 8.1) совпадает с дисперсией ошибки потенциальной оценки.
Пример 8.8. Найдем потенциальную точность оценки на-
чальной фазы
(8.6) в выражение для
радиоимпульса (8.6). Подставив выражение
0
при
D
эфф
и проведя дифферен-
0
цирование, найдем
1
22

ln 2 1 2 ,
MNEq

эфф
0



00 0
 
где q
— отношение сигнал/шум; 0qEN
; E — энергия
сигнала;
T
22 2
cos

0
0
0,5
;
И
. Видно, что потенциальная
ft ft t
 

ЗЗ
EA ft dt A

точность оценки начальной фазы определяется отношением сигнал/шум и не зависит от формы сигнала.
Пример 8.9. Найдем потенциальную точность оценки
задержки
диоимпульс находится в пределах отрезка
радиоимпульса (8.6) по огибающей. Если ра-
0,T наблюде-
ния, задержка — неэнергетический параметр. Тогда энер­гия

сигнала не зависит от
T
2
ft dt


0
ЗЗ
0
. Дифференцируя это равенство по
, откуда следует
, имеем
TT

ft ft dt ft dt
 

00
22

    
ЗЗЗ ЗЗ
Подставив выражение (8.6) в выражение для
и дифференцируя его, получим формулу Вудворда:
2
0.
D
при
эфф
139
8. Измерение параметров сигналов

N

З

E
З
З


1
22
NSt dtNE q
 
00
22

где

DM

Зэфф


2, 212,
TT


00

2,
NMytSt dt

0

2
Aft ft dt
0
T


0
,,
S t dt S t dt

ЗЗ З
ная ширина спектра сигнала

ln



T

0
T
222



0

ЗЗ
2
2
ЗЗ
22

ЗЗЗ

,ЗSt. Видно, что потенциаль-
ЗЗ
22
;  — эффектив-
ная точность оценки задержки огибающей сигнала обратно пропорциональна отношению сигнал/шум и квадрату ширины его спектра.
Пример 8.10. Найдем потенциальную точность оценки за-
держки
радиоимпульса по фазе сигнала. Пусть излучается
гармонический сигнал, и задержка при приеме оценивается по фазе принятого сигнала
где
TiT
Для однозначности измерения фазы считаем
— известное натуральное число. Энергия
k
T

0
,cos
St A t
ЗЗ
;
2T
; i — известное натуральное число.
0
0
22
,2
St dt AT

З

сигнала не зависит от
0
при
0,tT ,
, то есть
kT
, где
0З
— неэнергетический параметр. Подставим выражение для
,ЗSt в выражение для

NStSt dtNAT q
00
140
DM
Зэфф

T
2, , 12.
0

NMytSt dt

0


ЗЗЗ
D
и найдем
эфф
1
22

ln

2,
0
22 222




T

ЗЗ
22

ЗЗ
 
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]