Параллельные вычисления в RLS-алгоритмах адаптивной фильтрации - page 9

...
χ( m)
N (k − L ) = xT
N1 (k − L − 1), xT
N2 (k − L − 1), . . . ,
xT
Nm ( k − L − 1), xT
Nm+1 ( k − L) , . . . , xT
NM (k − L )
T
,
...
χ( M )
N ( k − L) = xT
N1( k − L − 1), xT
N2(k − L − 1), . . . ,
xT
Nm ( k − L − 1), . . . , xT
NM (k − L − 1)
T
;
ρ(0)
N (k) = ρN ( k),
ρ(1)
N (k) = [pT
N1 (k − 1), pT
N2 (k), . . . , pT
Nm (k), . . . ,
pT
NM (k)] T ,
...
ρ( m)
N (k) = [pT
N1( k − 1), pT
N2 (k − 1), . . . , pT
Nm (k − 1),
pT
Nm +1 (k) , . . . , pT
NM (k)] T ,
...
ρ( m)
N (k) = [pT
N1( k − 1), pT
N2 (k − 1), . . . ,
pT
Nm (k − 1), . . . , pT
NM (k − 1)]T ;
ρ(0)
N (k − L) = ρN (k − L ),
ρ(1)
N (k − L) = [pT
N1( k − L − 1), pT
N2 (k − L ), . . . ,
pT
Nm ( k − L) , . . . , pT
NM (k − L)] T ,
...
ρ(m)
N ( k − L ) = [pT
N1 (k − L − 1), pT
N2 (k − L − 1), . . . ,
pT
Nm ( k − L − 1), pT
Nm +1 (k − L ), . . . , pT
NM (k − L )]T ,
...
ρ(m)
N ( k − L ) = [pT
N1 (k − L − 1), pT
N2 (k − L − 1), . . . ,
pT
Nm ( k − L − 1), . . . , pT
NM (k − L − 1)]T .
В быстром RLS-алгоритме матрицы S(m)
N +1 и T (m)
N+1 являются
перестановочными. Умножение векторов на эти матрицы не требует
арифметических операций. Некоторые детали аналогичного способа
38 ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. “Приборостроение”. 2006. № 1