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)%   michael.jachan@tuwien.ac.at 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;