tfarma_est_tfywu.m

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

function AA= tfarma_est_tfywu(Ayy, MMA, N)
% function AA= tfarma_est_tfywu(Ayy, MMA, N)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 

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

Ayy= Aalpha(N/2-3*LAR+1:N/2+3*LAR+1, N/2+MMA-MAR+1:N/2+MMA+MAR+1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Dimensions
MARmax= (size(Ayy, 2)-1)/2;
LAR   = (size(Ayy, 1)-1)/6;
LL    = 2*LAR+1;

% The output
AA= zeros(LL, MARmax+1, MARmax);
BB= zeros(LL,        1, MARmax);

% Fill the Toeplitz/block-Toeplitz generator rPM and the RHS BM
rPM= [];
for mu= MARmax-1:-1:-MARmax+1%This is the right one;)
   rPM= [rPM; fliplr(Ayy(3*LAR+2-LL:3*LAR+LL, MARmax+mu+1).')];
end;

BM = [];
for m= 1:MARmax
   BM= [BM; -Ayy(2*LAR+1:4*LAR+1, MARmax+m+1)];
end;

% Solve the system
AM= invert_waxkailath(rPM, BM);

% Reshape AM to A and Compute B
for k= 1:MARmax
   Aml_e= [[zeros(LAR, 1); 1; zeros(LAR, 1)] reshape(AM(1:LL*k, 1, k), LL, k)];
   AA(:, :, k)= [Aml_e zeros(LL, MARmax-k)];
%   B0l= tfar_est_b0l(Ayy(:, MARmax+1-k:MARmax+1+k), Aml_e, N);
%   BB(:, :, k)= B0l;
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
param_norm(AA(:, :, MAR), Aml, 0)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%