Информативными параметрами в (7) являются задержка
τ
и допле-
ровское смещение частоты
ω
д
. Случайная фаза
φ
k
−
1
является неин-
формативным параметром, поэтому для определения дискриминатора
некогерентного приемника будем использовать следующую функцию
правдоподобия, усредненную по случайной фазе
φ
к
,
f Y
M
1
|
λ
=
c
1
I
0
˜
X
(
T, τ
k
−
1
,l
, ω
l
,k
−
1
,l
)
,
(8)
где
c
1
— нормировочная константа;
λ
— измеряемый параметр;
I
0
(
x
)
—
модифицированная функция Бесселя нулевого порядка;
˜
X
(
τ
k
−
1
,l
, ω
д
,k
−
1
,l
) =
q
˜
I
2
(
τ
k
−
1
,l
, ω
д
,k
−
1
,l
) + ˜
Q
2
(
τ
k
−
1
,l
, ω
д
,k
−
1
,l
);
˜
I
(
τ
k
−
1
,l
, ω
д
,k
−
1
,l
) =
=
A
σ
2
n
M
X
l
=1
y
(
t
k
−
1
,l
) cos (
ω
0
t
k
−
1
,l
+
ω
д
,k
−
1
,l
(
l
−
1)
T
d
+
φ
k
−
1
,l
)
,
˜
Q
(
τ
k
−
1
,l
, ω
д
,k
−
1
,l
) =
=
A
σ
2
n
M
X
l
=1
y
(
t
k
−
1
,l
) sin (
ω
0
t
k
−
1
,l
+
ω
д
,k
−
1
,l
(
l
−
1)
T
d
+
φ
k
−
1
,l
)
.
Соотношение, описывающее дискриминатор задержки огибающей
сигнала в приемнике, находится по (8) с использованием аппроксима-
ции [6]
u
д
τ
(
t
k
) =
−
I
1
˜
X
(˜
τ
k
,
˜
ω
д
,k
)
I
0
˜
X
(˜
τ
k
,
˜
ω
д
,k
)
∂
˜
X
(˜
τ
k
,
˜
ω
д
,k
)
∂τ
,
где
I
1
(
z
) =
I
0
0
(
z
)
— модифицированная функция Бесселя первого по-
рядка,
˜
τ
k
,
˜
ω
д
,k
— экстраполированные на
k
-й такт оценки параметров
τ
k
, ω
д
,k
соответственно.
Алгоритм вычисления оценок задержки сигнала описывается раз-
ностными уравнениями:
ˆ
x
τ,k
= ˜
x
τ,k
+ K
τ
u
д
τ,k
−
1
;
˜
x
τ,k
= Φ
τ
ˆ
x
τ,k
−
1
,
где
ˆ
x
τ,k
— текущая оценка вектора состояния;
˜
x
τ,k
— экстраполирован-
ная оценка вектора состояния;
K
τ
— вектор-столбец коэффициентов
усиления;
Φ
τ
— переходная матрица сглаживающего фильтра в контуре
следящей системы;
u
д
τ,k
— процесс на выходе дискриминатора.
На рис. 5 приведена модель ССЗ в SystemVue. При моделировании
задавались следующие параметры входного сигнала: частота несущей
130 ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. “Приборостроение”. 2014. № 4
