построения быстрых многоканальных RLS-алгоритмов адаптивной
фильтрации с неодинаковым числом весовых коэффициентов в кана-
лах адаптивного фильтра с использованием перестановочных матриц
рассмотрены в работе [23]. Перестановочные матрицы для случая
M
= 3
приведены в работе [8].
Вычислительные процедуры параллельных версий других бы-
стрых RLS-алгоритмов приведены ниже: Fast Transversal Filter (FTF),
Fast a Posteriori Error Sequential Technique (FAEST) и стабилизирован-
ного FAEST-алгоритма [5, 22, 24].
Параллельный FTF-алгоритм
Шаг
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 :
E
f
(
m
)
N
(
k
) =
λE
f
(
m
)
N
(
k
−
1) + e
f
(
m
)
F
(
k
)
α
f
(
m
)
H
F
(
k
)
Шаг
10 : ˜Φ
(
m
)
F
(
k
) = Φ
(
m
)
F
(
k
) I
F
−
α
f
(
m
)
H
F
(
k
)e
f
(
m
)
F
(
k
)
E
f
(
m
)
N
(
k
)
Шаг
11 : Φ
(
m
−
1)
F
(
k
) = [I
F
−
λ
−
1
˜Φ
(
m
)
F
(
k
)˜q
(
m
)
H
F
(
k
)
α
b
(
m
)
F
(
k
)]
−
1
˜Φ
(
m
)
F
(
k
)
Шаг
12 : e
b
(
m
)
F
(
k
) =
α
b
(
m
)
F
(
k
)Φ
(
m
−
1)
F
(
k
)
Шаг
13 :
E
b
(
m
)
N
(
k
) =
λE
b
(
m
)
N
(
k
−
1) + e
b
(
m
)
F
(
k
)
α
b
(
m
)
H
F
(
k
)
Шаг
14 : h
b
(
m
)
N
(
k
) = h
b
(
m
)
N
(
k
−
1) +
λ
−
1
˜T
(
m
−
1)
NF
(
k
)e
b
(
m
)
H
F
(
k
)
End for
k
ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. “Приборостроение”. 2006. № 1 39