nonwssus_stat.m

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

function RS= nonWSSUS_stat(TAPS, W, N)
% function 
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Implements the 4D TF-correlation function of non-WSSUS mobile
% radio channels as function of the tap-TFARMA models and the
% innovations weights W[m1, l1; m2, l2]. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
N    = 256;
re_im= 'i';
mo_no= 'n';


Mch= 1;
Lch= 1;

% Tap 0:
MAR  =   2;
LAR  =   2;
MMA  =   1;
LMA  = LAR;
tfpm_file_gen;
%-------------
Aml0= Aml;
Bml0= Bml;

% Tap 1:
MAR  =   3;
LAR  =   2;
MMA  =   2;
LMA  = LAR;
tfpm_file_gen;
%-------------
Aml1= Aml;
Bml1= Bml;

TAPS= {{Aml0 Bml0} {Aml1 Bml1}}
W= zeros(Mch+1, 2*Lch+1, Mch+1, 2*Lch+1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Mch= size(TAPS, 2)-1;
Lch= (size(W, 2)-1/2)
E= randn(N*(Mch
+1), 1);