init_mcfdadjlms.m

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

% [NB,W,w,x,y,d,e,p,S,SE,yF,feF] = init_mcfdadjlms(NC,NL,Nref,Nact,Nsens,s,se,w0,xn0,d0)
%
%     Creates and initializes the variables required for 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 [size] ::
%     NC    : Required controller length in time domain
%     NL    : new samples per block
%     Nref  : number of reference signals
%     Nact  : number of actuators
%     Nsens : number of sensors
%     s     : FIR model of the secondary paths [M x Nact x Nsens]
%     se    : estimate of secondary paths [N x Nact x Nsens]
%     w0    : Initial time domain matrix of filter coef. [NC x Nref x Nact]
%     xn0   : Initial input samples block [NL x Nref]
%     d0    : Initial primary signals block [NL x Nsens]
%
% Output parameters [default] ::
%     NB   : FFT length [2.^nextpow2(NC+NL-1)]
%     W    : frequency domain filter coefficients matrix [zeros]
%     w    : time domain filter coefficients matrix [zeros]
%     x    : overlap-save input matrix [zeros]
%     y    : filter output matrix in time domain [zeros]
%     d    : primary signals block [zeros]
%     e    : error signals block [e = d + ys]
%     p    : power estimate of x(n) in freq. domain
%     S    : frequency domain secondary paths [fft(s,NB)]
%     SE   : frequency domain estimated secondary paths [fft(se,NB)]
%     yF   : output matrix buffer [zeros] 
%     feF  : filtered error matrix buffer [zeros]
%
% SEE ALSO ASPTMCFDADJLMS, ANVC_MCFDADJLMS
% Output parameters ::
%     W    : vector of filter coefficients (f-domain) [NB x Nref x Nact]
%     w    : vector of filter coefficients (t-domain) [NC x Nref x Nact]
%     se   : Estimated secondary path (f-domain) [NB x Nact x Nsens]
%     x    : vector of overlapped input samples (t-domain) [NB x Nref]
%     y    : filter output vector (t-domain) [NL x Nact]
%     d    : desired response (t-domain) [NL x Nsens]
%     e    : error vector (t-domain) [NL x Nsens]
%     p    : estimate of input vector power (f-domain) [NB x Nref]
%     yF   : output vector in freq. domain [NB x Nact]
%     feF  : filtered-e vector in freq. domain [NB x Nact]

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