function AA= tfarma_est_tfywu(Ayy, MMA, N)% function AA= tfarma_est_tfywu(Ayy, MMA, N)%   This file is part of the TFPM toolbox v1.0 (c)%   michael.jachan@tuwien.ac.at and underlies the GPL.% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;N    = 128;MAR  =   2;LAR  =   1;MMA  =   1;LMA  = LAR;re_im= 'r';mo_no= 'n';tfpm_file_gen;%-------------alpha= 1/2;beta= 1/2;Aalpha= tfarma_ambi(Aml, Bml, N, alpha);tf_show(Aalpha)Ayy= Aalpha(N/2-3*LAR+1:N/2+3*LAR+1, N/2+MMA-MAR+1:N/2+MMA+MAR+1);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DimensionsMARmax= (size(Ayy, 2)-1)/2;LAR   = (size(Ayy, 1)-1)/6;LL    = 2*LAR+1;% The outputAA= zeros(LL, MARmax+1, MARmax);BB= zeros(LL,        1, MARmax);% Fill the Toeplitz/block-Toeplitz generator rPM and the RHS BMrPM= [];for mu= MARmax-1:-1:-MARmax+1%This is the right one;)   rPM= [rPM; fliplr(Ayy(3*LAR+2-LL:3*LAR+LL, MARmax+mu+1).')];end;BM = [];for m= 1:MARmax   BM= [BM; -Ayy(2*LAR+1:4*LAR+1, MARmax+m+1)];end;% Solve the systemAM= invert_waxkailath(rPM, BM);% Reshape AM to A and Compute Bfor k= 1:MARmax   Aml_e= [[zeros(LAR, 1); 1; zeros(LAR, 1)] reshape(AM(1:LL*k, 1, k), LL, k)];   AA(:, :, k)= [Aml_e zeros(LL, MARmax-k)];%   B0l= tfar_est_b0l(Ayy(:, MARmax+1-k:MARmax+1+k), Aml_e, N);%   BB(:, :, k)= B0l;end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%param_norm(AA(:, :, MAR), Aml, 0)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%