function [A, B]= tfar_est_Gren(Aalpha, MAR, LAR, N)% function [A, B]= tfar_est_Gren(Aalpha, MAR, LAR, N)%   This file is part of the TFPM toolbox v1.0 (c)%   michael.jachan@tuwien.ac.at and underlies the GPL.% % Estimates TFAR(MAR; LAR; N; 1/2) models from Ayy. % size(Ayy)= [6*LAR+1, 2*MAR+1], Ayy(3*LAR+1, MAR+1)= A_y[0, 0]. % Uses the Grenier estimator. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;MAR  =   2;LAR  =   1;MMA  =   0;LMA  = LAR;N    = 256;re_im= 'r';mo_no= 'n';tfpm_file_gen;%-------------alpha= 1/2;Aalpha= tfarma_ambi(Aml, Bml, N, alpha);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Ralpha= ifft([Aalpha(N/2+1:N, N/2-MAR+1:N/2+MAR+1); Aalpha(1:N/2, N/2-MAR+1:N/2+MAR+1)]);%.*[zeros(MAR, 2*MAR+1); ones(N-MAR, 2*MAR+1)];Ayy= Aalpha(N/2-3*LAR+1:N/2+3*LAR+1, N/2-MAR+1:N/2+MAR+1);% Generate R matrixR= zeros(MAR*(2*LAR+1));for n= MAR:N-1   for m= 1:MAR      for mprime= 1:MAR         for l= -LAR:LAR	    for lprime= -LAR:LAR               mm= m-mprime + MAR+1;               R( MAR*(LAR+l)+m, MAR*(LAR+lprime)+mprime )= ...               R( MAR*(LAR+l)+m, MAR*(LAR+lprime)+mprime ) + ...		   Ralpha(n-mprime+1, mm)*exp(j*2*pi*(lprime-l)*n/N) ...		                         *exp(j*2*pi*l*m/N);	    end;         end;      end;   end;end;% Generate r vectorr= zeros(MAR*(2*LAR+1), 1);for n= MAR:N-1   for m= 1:MAR      for l= -LAR:LAR         r( MAR*(LAR+l)+m )= r( MAR*(LAR+l)+m ) + Ralpha( n+1, m+MAR+1 )*exp(-j*2*pi*l*(n-m)/N);%+1!!      end;   end;end;theta= -inv(R)*r;A= [[zeros(LAR, 1); 1; zeros(LAR, 1)] flipud(rot90(reshape(theta, MAR, 2*LAR+1)))];% Reshape AM to A and Compute BB= tfar_est_b0l(Ayy, A, N);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%param_norm(A, Aml, 0)param_norm(B, Bml, 1)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%