param_get.m

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

function cml= param_get(Cml, m, l)
% function cml= param_get(Cml, m, l)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Delivers the (m, l)th parameter of the array Cml. Zero if (m, l)
% is outside of the range (1:M+1, -L:L). 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;%tfpm;
MAR  =   7;
LAR  =   5;
MMA  =   5;
LMA  =   3;
N    =  64;
re_im= 'i';
tfpm_file_gen;
%-------------
Cml= Aml;
m= 0;
l= 0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[M, L]= param_dim(Cml);
if(m<0 | m>M)
   cml= 0;
   return;
end;
if(l<-L | l>L)
   cml= 0;
   return;
end;
cml= Cml(L+1+l, 1+m);