function [Bnm, B, REG, PMIN]= tvma_est_cepsb(Aalpha, MMA, LMA)% function [B, REG, PMIN]= tvma_est_cepsb(Aalpha, MMA, LMA)%   This file is part of the TFPM toolbox v1.0 (c)%   michael.jachan@tuwien.ac.at and underlies the GPL.% % Estimates TVMA(MMA) models from Ralpha. size(Ralpha)= [N, N]. % Uses the complex TV cepstral recursion in IEEESPL vol9(4), apr02,% pp130. Psi is missing. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;MAR  =   0;LAR  =   0;MMA  =   3;LMA  =   2;N    = 256;re_im= 'r';mo_no= 'n';tfpm_file_gen;%-------------alpha= 1/2;beta = 1/2;%Psi= ones(N);Aalpha= tfarma_ambi(Aml, Bml, N, alpha);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Ralpha= ml_to_nm(Aalpha);% DimensionsN= size(Ralpha, 1);Bnm= [];for n= 0:N-1   [r_hat, REG, PMIN]= ti_ceps2(Ralpha(n+1, :).', 1e-4);%% Recursion in m% Initialization (m=0)%   BB= [exp(.5*r_hat(N/2+1)) zeros(1, MMA)];  BB= 1;% m<0   for m= 1:MMA      bb= 0;      for mprime= 0:m-1         bb= bb + (m-mprime)*BB(mprime+1)*...         r_hat(N/2+1+m-mprime);      end;      BB(m+1)= bb/m;   end;   Bnm= [Bnm; BB];end;Hmonic= [zeros(N, N/2) Bnm zeros(N, N/2-MMA-1)];Wmonic= abs(nm_to_nk(Hmonic)).^2;%imagesc(abs(Wmonic))B0= sqrt(real(mean( (nm_to_nk(Ralpha) ./ Wmonic).' ).'));%figure(1);clf;plot(real(B0));%figure(2);clf;plot(real(param_expand(Bml, N)))Bnm= Bnm.*(B0*ones(1, MMA+1));B= param_compress(Bnm, LMA);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%param_norm(param_expand(Bml, N), Bnm, 1)param_norm(Bml, B, 1)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%