Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Оптимальный прием и обработка сигналов. Учебное пособие
.pdf
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
ЗЗ
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]
