param_tconv_b2.m

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

function B2= param_tconv_b2(Bml1, Bml2, N, beta)
% function B2= param_tconv_b2(Bml1, Bml2, N, beta)
%   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;
MAR  =   7;
LAR  =   5;
MMA  =   5;
LMA  =   3;
N    = 128;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
beta= -1/2;
Bml1= Aml;
Bml2= Bml;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Dimensions:
[M1, L1]= param_dim(Bml1);
[M2, L2]= param_dim(Bml2);

M= M1+M2;
L= L1+L2;

B2= zeros((M+1)*(2*L+1), (M1+1)*(2*L1+1));

% Compute B2
for m= 0:M
   for l= -L:L
      for mp= 0:M1
         for lp= -L1:L1
	    B2((m*(2*L+1))+l+L+1, mp*(2*L1+1)+lp+L1+1)= ...
		param_get(Bml2, m-mp, l-lp) * exp(j*2*pi/N*tfpm_phi(m, l, mp, lp, -beta));
	 end;
      end;
   end;
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Test B2:
norm(param_destack_ml(B2*param_stack_ml(Bml1), M, L)-param_tconv(Bml1, Bml2, N, beta))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%