Параллельный FTF-алгоритм
Шаг
15 :
α
F
(
k
) = d
F
(
k
)
−
h
H
N
(
k
−
1)X
NF
(
k
)
Шаг
16 : e
F
(
k
) =
α
F
(
k
)Φ
(0)
F
(
k
)
Шаг
17 : h
N
(
k
) = h
N
(
k
−
1) +
λ
−
1
˜T
(0)
NF
(
k
)e
H
F
(
k
)
Шаг
18 : ˜T
(
M
)
NF
(
k
+ 1) = ˜T
(0)
NF
(
k
)
,
Φ
(
M
)
F
(
k
+ 1) = Φ
(0)
F
(
k
)
End for
k
Параллельный FAEST-алгоритм отличается от параллельного FTF-
алгоритма тем, что в первом рекурсивно вычисляются матрицы
[Φ
(
m
)
F
(
k
)]
−
1
, а во втором — матрицы
Φ
(
m
)
F
(
k
)
. Аналогичные матри-
цы в последовательных быстрых RLS-алгоритмах отсутствуют. Они
становятся скалярными величинами, известными в теории адаптивной
фильтрации как отношения правдоподобия.
Параллельный FAEST-алгоритм
Шаг
0 :
Инициализация:
χ
N
(0) = 0
N
, . . . , χ
N
(0
−
L
+ 1) = 0
N
,
ρ
N
(0) = 0
N
, . . . , ρ
N
(0
−
L
+ 1) = 0
N
,
d
(0) = 0
, . . . , d
(0
−
L
+ 1) = 0
,
X
NF
(0) = O
NF
,
h
N
(0) = 0
N
, E
f
(
m
)
N
(0) =
δ
2
, E
b
(
m
)
N
(0) =
δ
2
λ
−
N
m
,
h
f
(
m
)
N
(0) = 0
N
,
h
b
(
m
)
N
(0) = 0
N
, m
= 1 :
M,
˜T
(
m
)
NF
(1) = 0
NF
,
Φ
(
M
)
F
(1) = S
F
For
k
= 1
,
2
, . . . , K
For
m
=
M, M
−
1
, . . . ,
1
Шаг
1 :
α
f
(
m
)
F
(
k
) = x
(
m
)
F
(
k
)
−
h
f
(
m
)
H
N
(
k
−
1)X
(
m
)
NF
(
k
)
Шаг
2 : e
f
(
m
)
F
(
k
) =
α
f
(
m
)
F
(
k
)Φ
(
m
)
F
(
k
)
Шаг
3 :
η
m
)
F
(
k
) =
α
f
(
m
)
F
(
k
)
E
f
(
m
)
N
(
k
−
1)
Шаг
4 :
˜T
(
m
)
(
N
+1)
F
(
k
) =
1
−
h
f
(
m
)
N
(
k
−
1)
η
(
m
)
F
(
k
) +
0
T
F
˜T
(
m
)
NF
(
k
)
Шаг
5 : S
(
m
)
N
+1
T
(
m
)
T
N
+1
{
˜T
(
m
)
(
N
+1)
F
(
k
)
}
=
˜Q
(
m
)
NF
(
k
)
˜q
(
m
)
F
(
k
)
Шаг
6 : h
f
(
m
)
N
(
k
) = h
f
(
m
)
N
(
k
−
1) +
λ
−
1
T
(
m
)
NF
(
k
)e
f
(
m
)
H
F
(
k
)
Шаг
7 :
α
b
(
m
)
F
(
k
) = ˜q
(
m
)
F
(
k
)
E
b
(
m
)
N
(
k
−
1)
Шаг
8 :
˜T
(
m
−
1)
NF
(
k
) = ˜Q
(
m
)
NF
(
k
) + h
b
(
m
)
N
(
k
−
1)˜q
(
m
)
F
(
k
)
Шаг
9 : [ ¯Φ
(
m
)
F
(
k
)]
−
1
= [Φ
(
m
)
F
(
k
)]
−
1
+
λ
−
1
α
f
(
m
)
H
F
(
k
)
η
(
m
)
F
(
k
)
Шаг
10 :
E
f
(
m
)
N
(
k
) =
λE
f
(
m
)
N
(
k
−
1) + e
f
(
m
)
F
(
k
)
α
f
(
m
)
H
F
(
k
)
Шаг
11 : [Φ
(
m
−
1)
F
(
k
)]
−
1
= [ ¯Φ
(
m
)
F
(
k
)]
−
1
−
λ
−
1
α
b
(
m
)
H
F
(
k
)q
(
m
)
F
(
k
)
Шаг
12 : e
b
(
m
)
F
(
k
) =
α
b
(
m
)
F
(
k
)Φ
(
m
−
1)
F
(
k
)
40 ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. “Приборостроение”. 2006. № 1
