nonwssus_rs.m

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

function RS= nonWSSUS_RS(Lalpha)
% function RS= nonWSSUS_RS(Lalpha)
%   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
alpha= 1/2;
N= 256;
J= 5;
B= J+1;
offset= 200;
M= 4;
L= 3;
filename= sprintf('chpar%02d%04d%02d%02d-%04d', J, N, M, L, offset)
load(filename)
BML= B0L;
Lalpha= vtfarma_weyl(AML, BML, N, B, 1/2);
for j= 0:J
   for jp= 0:J
      mesh(abs(reshape(Lalpha(j+1, jp+1, :, :), N, N)))
      drawnow
   end;
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Dimensions
[J, J, N, N]= size(Lalpha);
J= J-1;

RS= zeros(J+1, N, J+1, N);

for tau= 0:J
   for Dtau= 0:tau
      rs= zeros(N);
      for taup= 0:J
         L1= reshape(Lalpha(tau+1     , taup+1, :, :), N, N);
         L2= reshape(Lalpha(tau-Dtau+1, taup+1, :, :), N, N);
         rs= rs + L1.*conj(L2);
      end;
      rs1= fftshift(fft(rs), 1);
      RS(tau+1, :, Dtau+1, :)= reshape(rs1.', 1, N, 1, N);
   end;
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
RRSS= fftshift(ifft(fftshift(RS, 4), [], 4), 4);

% PLOT LSF @ k= 0
k= 0;
for n= 0:N-1
   rrss= fftshift(fft(RRSS(:, :, :, n+1), N, 3), 3);
   n
   figure(1);mesh(real(   rrss(:, :, k+N/2+1)))
   figure(2);tf_show(abs(rrss(:, :, k+N/2+1)))
%   tf_show(rrss(:, :, k+N/2+1))
   drawnow
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%