param_weyl.m

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

function Lalpha= param_weyl(Cml, N, alpha)
% function Lalpha= param_weyl(Cml, N, alpha)
%   This file is part of the TFPM toolbox v0.5 (c)
%   [e-mail] and underlies the GPL.
% 
% Computes the generalized Weyl symbol of an underspread TFMA
% model given through Cml. 

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

%Dimensions:
[M, L]= param_dim(Cml);

% Select right alpha:
Cml= Cml.*param_phase(M, L, N, alpha);
Salpha= [ zeros(N/2-L, N); 
     [zeros(2*L+1, N/2) N*Cml zeros(2*L+1, N/2-1-M)] 
     zeros(N/2-L-1, N)
   ];
Lalpha= nm_to_nk(ml_to_nm(Salpha));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
H0= ker_to_lag(tvarma_impr(Amn, Bmn), 0, alpha);
L0= nm_to_nk(H0);
figure(97);mesh(abs(L0))
figure(98);mesh(abs(Lalpha))
figure(99);mesh(abs(L0-Lalpha))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%