tfpm_underspread.m

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

% 
%   This file is part of the TFPM toolbox v0.5 (c)
%   [e-mail] and underlies the GPL.
% 
% Playing round with underspread approximations... 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N    = 512;
MAR  =   1;
LAR  =   1;
MMA  =   1;
LMA  =   1;
re_im= 'r';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
beta = 1/2;

He= [];
for n= 0:N-1
   e= zeros(N, 1);
   e(n+1)= 1;
   y= real(tfarma_gen(e, Aml, Bml, beta));
   plot(y);drawnow
   He= [He y];
end;

He= ker_to_lag(He, -1, alpha);

Hu= real(tfarma_impr(Aml, Bml, N, alpha));

figure(1);clf;mesh(He(:, N/2+1:end))
view(30, 30);
axis([1 N/2 1 N min(min(He)) max(max([He Hu]))])
set(gca, 'XTick', [1 N/2])
set(gca, 'YTick', [1 N/2 N])
%set(gca, 'ZTick', [-1 0 1 2])

figure(2);clf;mesh(Hu(:, N/2+1:end))
view(30, 30)
axis([1 N/2 1 N min(min(He)) max(max([He Hu]))])
set(gca, 'XTick', [1 N/2])
set(gca, 'YTick', [1 N/2 N])
%set(gca, 'ZTick', [-1 0  2])

figure(3);clf;mesh(abs(He(:, N/2+1:end)-Hu(:, N/2+1:end)))
view(30, 30)
axis([1 N/2 1 N-1 0 max(max(abs(He(:, N/2+1:end)-Hu(:, N/2+1:end))))])
set(gca, 'XTick', [0 N/2-1])
set(gca, 'YTick', [0 N/2-1 N-1])
%set(gca, 'ZTick', [-.5 0 .5])

Le= nm_to_nk(He);
Lu= tfarma_weyl(Aml, Bml, N, alpha);

figure(4);clf;mesh(abs(Le))
view(30, 30);
axis([1 N 1 N 0 max(max(abs([Le Lu])))])
set(gca, 'XTick', [1 N/2 N])
set(gca, 'YTick', [1 N/2 N])

figure(5);clf;mesh(abs(Lu))
view(30, 30)
axis([1 N 1 N 0 max(max(abs([Le Lu])))])
set(gca, 'XTick', [1 N/2 N])
set(gca, 'YTick', [1 N/2 N])

figure(6);clf;mesh(abs(Le-Lu))
view(30, 30)
axis([1 N 1 N 0 max(max(abs(Le-Lu)))])
set(gca, 'XTick', [1 N/2 N])
set(gca, 'YTick', [1 N/2 N])