ml_plot.m

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

function ml_plot(Aalpha, fig)
% function 
%   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;
N    =  16;
MAR  =   3;
LAR  =   0;
MMA  =   2;
LMA  =   0;
re_im= 'r';
mo_no= 'm';
tfpm_file_gen;
%-------------
alpha= 1/2;
Aalpha= tfarma_sprf(Aml, Bml, N, alpha);
fig= 1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if(nargin==1)
   fig= 1;
end

N= size(Aalpha, 1);

figure(fig);clf
imagesc(abs(flipud(Aalpha)))
xlabel('m');
ylabel('l');
set(gca, 'XTick', [1 N/2+1 N])
set(gca, 'XTickLabel', [-N/2 0 N/2-1])
set(gca, 'YTick', [1 N/2 N])
set(gca, 'YTickLabel', [N/2-1 0 -N/2])
drawnow;