sym_ambi.m

用于模拟时变非平稳的ARMA过程

function Balpha= sym_ambi(Aalpha)
% function Balpha= sym_ambi(Aalpha)
%   This file is part of the TFPM toolbox v0.5 (c)
%   [e-mail] and underlies the GPL.
% 
% Computes the time-lag and frequency-lag reversed and conjugated
% ambiguity function with phase factor for alpha= 1/2!! 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   0;
LAR  =   0;
MMA  =   5;
LMA  =   3;
N    =  64;
re_im= 'r';
mo_no= 'm';
tfpm_file_gen;
%-------------
alpha= 1/2;
beta = 1/2;

%Aalpha= tfarma_ambi(Aml, Bml, N, alpha);

y= tfarma_gen(randn(N, 1), Aml, Bml, beta);
%y= [N/2:-1:1 zeros(1, N/2)].';
%y= [1 1 zeros(1, N-2)];
qYY= corr_est(y, y,   -1, alpha);
Aalpha= nm_to_ml(qYY).*tf_multiwin(N, MMA, 2*LMA, 6, min(3, N/32), 2);

Aalpha= tf_multiwin(N, MMA, 2*LMA, 6, min(3, N/32), 2);

figure(97);cesh(Aalpha)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

N= size(Aalpha, 1);
Aminmin= fft_timereverse(fft_timereverse(Aalpha.').');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(98);clf;cesh(Aminmin)
figure(99);clf;cesh(Aalpha-Aminmin)
norm(Aalpha-Aminmin)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

ml= (-N/2:N/2-1).'*(-N/2:N/2-1);
Balpha= conj(Aminmin).*exp(-2*pi*j/N*ml);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(98);clf;cesh(Balpha)
figure(99);clf;cesh(Aalpha-Balpha)
norm(Aalpha-Balpha)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%