tf_ceps2.m

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

function [Aalpha_hat, REGULARIZED, PMIN]= tf_ceps2(Aalpha, ABSMIN)
% function [Aalpha_hat, REGULARIZED, PMIN]= tf_ceps2(Aalpha, ABSMIN)
%   This file is part of the TFPM toolbox v0.5 (c)
%   [e-mail] and underlies the GPL.
% 
% Computes the complex TF cepstrum of the ambiguity function 
% A[m, l], m= -N/2:N/2-1, l= -N/2:N/2-1. ABSMIN is used to
% regularize before taking the log. Only the positive part of the
% Wigner-Ville Spectrum is used, the phase is set to zero! 

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

Psi= tf_multiwin(N, MMA, 2*LMA, 6, min(3, N/32), 2);
norm(Psi-sym_ambi(Psi))

y= tfarma_gen(randn(N, 1), Aml, Bml, beta);
Aalpha= ambi_est_cyc(y);
norm(Aalpha-sym_ambi(Aalpha))

Aalpha= Aalpha.*conj(Psi);
norm(Aalpha-sym_ambi(Aalpha))

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

Palpha= nm_to_nk(ml_to_nm(Aalpha));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(98);cesh(Palpha)
[max(max(abs(real(Palpha))))/max(max(abs(imag(Palpha))))]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Pabs= real(Palpha);

% Regularization: min(min(Palpha))->ABSMIN!
PMIN= min(min(Pabs));
mask= Pabs<ABSMIN;
Pabs= (1-mask).*Pabs + ABSMIN*mask;
REGULARIZED= sum(sum(mask));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(98);clf;
   subplot(2, 1, 1);mesh(Pang)
   subplot(2, 1, 2);mesh(Pang1)
PMIN
REGULARIZED
figure(99);clf;
   subplot(2, 1, 1);mesh(Pang)
   subplot(2, 1, 2);mesh(Pang-mod(Pang1, 2*pi))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% The real log:
Palpha_hat= log(Pabs);

Aalpha_hat= nm_to_ml(nk_to_nm(Palpha_hat));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(99);cesh(Aalpha_hat)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%