simu_tfarma.m

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

function simu_tfarma(M, L, N, re_im, mo_no, MM)
% function simu_tfarma(M, L, N, re_im, mo_no, MM)
%   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    =  64;
M    =   2;
L    =   2;
re_im= 'i';
mo_no= 'n';
MM   =  10;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N0= N;
N= 128;
MAR  =   M;
LAR  =   L;
MMA  = M-1;
LMA  =   L;
tfpm_file_gen;
%-------------
N= N0;
alpha= 1/2;
beta = 1/2;

Bml= Bml./param_get(Bml, 0, 0);

%exp= 'TFARMA';
exp= 'tfarma';
resultname= sprintf('data/%04d/%s%d%d%s%s.mat', N, exp, M, L, re_im, mo_no)

% SNR values
SNR= [inf];

par= [Aml(:, 2:end); param_expand(Bml, N)];
ENERGY= tfarma_ambi(Aml, Bml, N, alpha);ENERGY= ENERGY(N/2+1, N/2+1);

[PAR, FLO, STA, AVG, MSE, VAR, BIQ, EInfo]= simulate(exp, N, SNR, par, ENERGY, re_im, MM, MAR, LAR, MMA, LMA);

save(resultname, 'PAR', 'par', 'FLO', 'STA', 'AVG', 'MSE', 'VAR', 'BIQ', 'MM', 'SNR', 'ENERGY', 'EInfo');