tvma_est_cepsb.m

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

function [Bnm, B, REG, PMIN]= tvma_est_cepsb(Aalpha, MMA, LMA)
% function [B, REG, PMIN]= tvma_est_cepsb(Aalpha, MMA, LMA)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Estimates TVMA(MMA) models from Ralpha. size(Ralpha)= [N, N]. 
% Uses the complex TV cepstral recursion in IEEESPL vol9(4), apr02,
% pp130. Psi is missing. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   0;
LAR  =   0;
MMA  =   3;
LMA  =   2;
N    = 256;
re_im= 'r';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
beta = 1/2;
%Psi= ones(N);
Aalpha= tfarma_ambi(Aml, Bml, N, alpha);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ralpha= ml_to_nm(Aalpha);

% Dimensions
N= size(Ralpha, 1);

Bnm= [];
for n= 0:N-1
   [r_hat, REG, PMIN]= ti_ceps2(Ralpha(n+1, :).', 1e-4);

%% Recursion in m
% Initialization (m=0)
%   BB= [exp(.5*r_hat(N/2+1)) zeros(1, MMA)];
  BB= 1;
% m<0
   for m= 1:MMA
      bb= 0;
      for mprime= 0:m-1
         bb= bb + (m-mprime)*BB(mprime+1)*...
         r_hat(N/2+1+m-mprime);
      end;
      BB(m+1)= bb/m;
   end;
   Bnm= [Bnm; BB];
end;

Hmonic= [zeros(N, N/2) Bnm zeros(N, N/2-MMA-1)];
Wmonic= abs(nm_to_nk(Hmonic)).^2;
%imagesc(abs(Wmonic))


B0= sqrt(real(mean( (nm_to_nk(Ralpha) ./ Wmonic).' ).'));
%figure(1);clf;plot(real(B0));
%figure(2);clf;plot(real(param_expand(Bml, N)))
Bnm= Bnm.*(B0*ones(1, MMA+1));
B= param_compress(Bnm, LMA);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
param_norm(param_expand(Bml, N), Bnm, 1)
param_norm(Bml, B, 1)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%