tfma_fita.m

来自「用于模拟时变非平稳的ARMA过程」· M 代码 · 共 53 行

function E= tfma_fitA(x, Mmax, Lmax, Psi, lambda)
% function E= tfma_fitA(x, Mmax, Lmax, Psi, lambda)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Searches TFMA(M, L) models for M= 1:Mmax, L= 0:Lmax. Uses CTFC
% estimator. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
N    =  64;
MAR  =   0;
LAR  =   0;
MMA  =   2;
LMA  =   1;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
beta = 1/2;
lambda= .9800;

Mmax= max([2*MMA 2]);
Lmax= max([2*LMA 2]);
%Psi= ones(N);
Psi= tf_multiwin(N, 3*Mmax, 4*Lmax, 0, 2, 1);
x= tfarma_gen(randn(N, 1), Aml, Bml, beta);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

alpha= 1/2;
beta = 1/2;
N= length(x);

Ax= fft(corr_est(x, x, -1, alpha));
Ax= [ Ax(N/2+1:N, :); Ax(1:N/2, :)].*conj(Psi);

E= 99999*ones(Mmax, Lmax+1);

for L= 0:Lmax
   [B, REG, PMIN]= tfma_est_cepsb(Ax, Mmax, L);
   for M= 1:Mmax
      [M L];
      [Bml1, lambdamax, mm, nrFDIR]= param_stabilize(B(:, 1:M+1), N, lambda);
      eprime= tfarma_inv(x, 1, Bml1/param_get(Bml1, 0, 0));
      E(M, L+1)= eprime'*eprime/N;
   end;
end;