param_tdir.m

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

function [TDIRnm, TDIRn0, Cnm]= param_tdir(Cml, N)
% function [TDIRnm, TDIRn0, Cnm]= param_tdir(Cml, N)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Computes root locations (time-depenent instantaneous-roots) of
% the TF parameters Cml. TDFRn0 is the normalization, so that we
% work with monic Cml. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   3;
LAR  =   2;
MMA  =   3;
LMA  =   2;
N    =  64;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
Cml= Bml;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Cnm= param_expand_nm(Cml, N);

% TDFR:
TDIRnm= [];
TDIRn0= [];
for n= 0:N-1
   cf0= Cnm(n+1, 1);
   r= roots(Cnm(n+1, :)/cf0);
   TDIRnm= [TDIRnm; r.'];
   TDIRn0= [TDIRn0; cf0];
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[TDIRnm, TDIRn0]= param_tdir(ALTI, N)
[TDIRnm, TDIRn0]= param_tdir(BLTI, N)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%