www.pudn.com > tfarma10.rar > ambi.est.m, change:2004-12-29,size:1575b

function Qalpha= ambi.est(x, y, MAX_M_LAG, alpha) % function Qalpha= corr_est(x, y, MAX_M_LAG, alpha) % This file is part of the TFPM toolbox v1.0 (c) % michael.jachan@tuwien.ac.at and underlies the GPL. % % Computes the cyclic outer signal product of x[n] and y[n] for % alpha= +/-1/2. if(MAX_M_LAG==-1) then Qalpha is a full NxN % FFT-ready matrix. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear;tfpm; MAR = 7; LAR = 5; MMA = 0; LMA = 0; N = 16; re_im= 'r'; mo_no= 'm'; tfpm_file_gen; %------------- x= randn(N, 1); %x= (1:N).'; y= x; MAX_M_LAG= 2; alpha = 1/2; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Dimensions: N= length(x); Qalpha= corr_est_non([x(:);x(:);x(:)], [y(:);y(:);y(:)], MAX_M_LAG, alpha); if(MAX_M_LAG==-1) Qalpha= Qalpha(N+1:2*N, N+1:2*N); else Qalpha= Qalpha(N+1:2*N, :); end; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% size(Qalpha) Q0= ker_to_lag(x*y', MAX_M_LAG, alpha); norm(Q0-Qalpha) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%