tfarma_fitb.m

用于模拟时变非平稳的ARMA过程

function E= tfarma_fitB(x, Mmax, Lmax, Psi, lambda)
% function E= tfarma_fitB(x, Mmax, Lmax, Psi, lambda)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Searches TFARMA(M, L) models for M= 1:Mmax, L= 0:Lmax. Uses the
% linear estimator. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
N    = 128;
MAR  =   2;
LAR  =   2;
MMA  =   1;
LMA  = LAR;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
beta = 1/2;
lambda= .9800;

Mmax= max([2*MAR 2]);
Lmax= max([2*LAR 2]);
%Psi= ones(N);
Psi= tf_multiwin(N, 6*Mmax, 4*Lmax, 0, 2, 1);
tf_show(Psi)
x= tfarma_gen(randn(N, 1), Aml, Bml, beta);
L=0;
M=2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

alpha= 1/2;
beta = 1/2;
N= length(x);

Ax= fft(corr_est(x, x, -1, alpha));
%Ax= [ Ax(N/2+1:N, :); Ax(1:N/2, :)].*conj(Psi);
Ax= [ Ax(N/2+1:N, :); Ax(1:N/2, :)];%NO KERNEL 4 ETFYW!!!!!!!!!!!

E= 99999*ones(Mmax, Lmax+1);

for L= 0:Lmax
   for M= 2:Mmax
      [M L]
      AA= tfarma_est_tfywu(Ax(N/2-3*L+1:N/2+3*L+1, N/2:N/2+2*M), M-1, N);
      [Aml1, lambdamax, mm, nrFDIR]= param_stabilize(AA(:, :, end), N, lambda);
      BB= tfarma_est_lin(x, Psi, Aml1, M-1, L, 2*M-2, 2*L);
      [Bml1, lambdamax, mm, nrFDIR]= param_stabilize(BB(:, :, end), N, lambda);
      eprime= tfarma_inv(x, Aml1, Bml1/param_get(Bml1, 0, 0));
      E(M, L+1)= eprime'*eprime/N;
   end;
end;