asptmcfdadjlms.m

卡尔曼滤波器设计的一个例子

% [W,w,x,y,e,p,yF,feF] = asptmcfdadjlms(NC,W,x,xn,dn,yF,feF,S,SE,p,mu,b,c)
%
%     Performs filtering and coefficient update using the 
%     Multi Channel Frequency Domain Adjoint LMS (MCFDADJLMS) algorithm
%     for use in Active Noise & Vibration Control (ANVC) applications.
%
%     BLOCK DIAGRAM      
%                                                d(n)
%                                       ----------->-----.
%                                                        |
%            Nref                        y(n)      ys(n) V 
%      x(n) --/--->.-------------->[ w ]--/-->[ s ]---->(+)
%                  |                 ^   Nact            |
%                  |                 |                   |e(n)
%                  |                 |   fe(n)           |
%                  '-------------->[LMS]<-----[ se ]<-/--'
%                                                     Nsens
%
% Input Parameters :: 
%     NC   : controller length in time domain
%     W    : freq. domain filter coefficients matrix [NB x Nref x Nact]
%     x    : previous overlap-save input matrix [NB x Nref]
%     xn   : new input samples block [NL x Nref]
%     dn   : new primary samples block [NL x Nsens]
%     yF   : previous buffer of y(n) [NB x Nact]
%     feF  : previous buffer of fe(n) [NB x Nact]
%     S    : FIR model of the secondary paths [NB x Nact x Nsens]
%     SE   : estimated FIR model S [NB x Nact x Nsens]
%     p    : last estimated power of x(n) [NB x Nref]
%     mu   : adaptation constant
%     b    : pole of AR filter used in estimating p 
%     c    : if not zero, uses the constrained BFDAF algorithm.
%
% Output parameters ::
%     W    : updated frequency domain filter coefficients  
%     w    : updated time domain filter coefficients  
%     x    : updated overlap-save input matrix
%     y    : controller output block 
%     e    : new error block 
%     p    : updated estimate of power of x(n)
%     yF   : updated output buffer 
%     feF  : updated filtered error buffer 
%
% SEE ALSO INIT_MCFDADJLMS, ANVC_MCFDADJLMS

%       Author : John Garas PhD.
%       Version 2.1, Release October 2002.
%       Copyright (c) DSP ALGORITHMS 2000-2002.