plot_tfpaper1.m

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

clear;tfpm;
re_im= 'i';
mo_no= 'm';

M0= 3;
L0= 2;
N0= 256;
fig1= 1;
fig2= 2;
fig3= 3;

nmin= 6;
Mmin= 1;
Lmin= 0;

nmax= 9;
Mmax= 7;
Lmax= 4;



exp= 'tfma1';
filename= sprintf('/users/mjachan/tex/prj/TF/figures/%s%02d%02d%02d%1s%1s', exp, M0, L0, log2(N0), re_im, mo_no)

figure(fig1);clf;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot MSE/VAR/BIQ over N:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MSEn= [];
VARn= [];
BIQn= [];
for n= nmin:nmax;
   N= 2^n;
   resultname= sprintf('data/%04d/%s%d%d%s%s.mat', N, exp, M0, L0, re_im, mo_no)
   load(resultname);
   esti= size(EInfo, 1);
   
   MSEnee= [];
   VARnee= [];
   BIQnee= [];
   for ee= 1:esti
      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))/sum(sum(abs(par).^2))];
%      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))];
%      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))];
%      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))];
   end;
   MSEn= [MSEn MSEnee];
   VARn= [VARn VARnee];
   BIQn= [BIQn BIQnee];
end;
MSEn= 10*log10(MSEn);
VARn= 10*log10(VARn);
BIQn= 10*log10(BIQn);

MINn= min(min([MSEn VARn BIQn]));
MINn= floor(10*MINn)/9.5;

MSEm= [];
VARm= [];
BIQm= [];
for M= Mmin:Mmax;
   resultname= sprintf('data/%04d/%s%d%d%s%s.mat', N0, exp, M, L0, re_im, mo_no)
   load(resultname);

   MSEnee= [];
   VARnee= [];
   BIQnee= [];
   for ee= 1:esti
      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))/sum(sum(abs(par).^2))];
%      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))];
%      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))];
%      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))];
   end;
   MSEm= [MSEm MSEnee];
   VARm= [VARm VARnee];
   BIQm= [BIQm BIQnee];
end;
MSEm= 10*log10(MSEm);
VARm= 10*log10(VARm);
BIQm= 10*log10(BIQm);

MINm= min(min([MSEm VARm BIQm]));
MINm= ceil(10*MINm)/9.5;

MSEl= [];
VARl= [];
BIQl= [];
for L= Lmin:Lmax;
   resultname= sprintf('data/%04d/%s%d%d%s%s.mat', N0, exp, M0, L, re_im, mo_no)
   load(resultname);

   MSEnee= [];
   VARnee= [];
   BIQnee= [];
   for ee= 1:esti
      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))/sum(sum(abs(par).^2))];
      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))/sum(sum(abs(par).^2))];
%      MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, ee)))];
%      VARnee= [VARnee; sum(sum(VAR(:, :, 1, ee)))];
%      BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, ee)))];
   end;
   MSEl= [MSEl MSEnee];
   VARl= [VARl VARnee];
   BIQl= [BIQl BIQnee];
end;
MSEl= 10*log10(MSEl);
VARl= 10*log10(VARl);
BIQl= 10*log10(BIQl);

MINl= min(min([MSEl VARl BIQl]));
MINl= ceil(10*MINl)/9.5;

MINMIN= min([MINn MINl MINm])
%MAXMAX= .55

subplot(3, 3, 1);cla;hold on
   plot(MSEn(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(MSEn(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(nmax~=nmin)
   axis([1 nmax-nmin+1 -17 2]);
   box;grid
end;
set(gca, 'XTick', 1:nmax-nmin+1)
set(gca, 'XTicklabel', 2.^(nmin:nmax))
set(gca, 'YTick', -15:5:0)
%set(gca, 'YTicklabel', ['D'; 'C'; 'B'; 'A'])
ylabel('MSE')
%xlabel('N')
%title(sprintf('M0= %d, L0= %d', M0, L0));

subplot(3, 3, 4);cla;hold on
   plot(VARn(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(VARn(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(nmax~=nmin)
   axis([1 nmax-nmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:nmax-nmin+1)
set(gca, 'XTicklabel', 2.^(nmin:nmax))
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
ylabel('VAR')
%xlabel('N')

subplot(3, 3, 7);cla;hold on
   plot(BIQn(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(BIQn(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(nmax~=nmin)
   axis([1 nmax-nmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:nmax-nmin+1)
set(gca, 'XTicklabel', 2.^(nmin:nmax))
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
ylabel('BIQ')
%xlabel('N')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot MSE/VAR/BIQ over M:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(3, 3, 2);cla;hold on
   plot(MSEm(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(MSEm(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Mmax~=Mmin)
   axis([1 Mmax-Mmin+1 -17 2]);
   box;grid
end;
set(gca, 'XTick', 1:Mmax-Mmin+1)
set(gca, 'XTicklabel', Mmin:Mmax)
set(gca, 'YTick', -15:5:0)
%set(gca, 'YTicklabel', ['D'; 'C'; 'B'; 'A'])
%ylabel('MSE')
%xlabel('M')
%title(sprintf('N0= %d, L0= %d', N0, L0));

subplot(3, 3, 5);cla;hold on
   plot(VARm(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(VARm(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Mmax~=Mmin)
   axis([1 Mmax-Mmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:Mmax-Mmin+1)
set(gca, 'XTicklabel', Mmin:Mmax)
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
%ylabel('VAR')
%xlabel('M')

subplot(3, 3, 8);cla;hold on
   plot(BIQm(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(BIQm(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Mmax~=Mmin)
   axis([1 Mmax-Mmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:Mmax-Mmin+1)
set(gca, 'XTicklabel', Mmin:Mmax)
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
%ylabel('BIQ')
%xlabel('M')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot MSE/VAR/BIQ over L:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(3, 3, 3);cla;hold on
   plot(MSEl(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(MSEl(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Lmax~=Lmin)
   axis([1 Lmax-Lmin+1 -17 2]);
   box;grid
end;
set(gca, 'XTick', 1:Lmax-Lmin+1)
set(gca, 'XTicklabel', Lmin:Lmax)
set(gca, 'YTick', -15:5:0)
%set(gca, 'YTicklabel', ['D'; 'C'; 'B'; 'A'])
%ylabel('MSE')
%xlabel('L')
%title(sprintf('N0= %d, M0= %d', N0, M0));

subplot(3, 3, 6);cla;hold on
   plot(VARl(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(VARl(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Lmax~=Lmin)
   axis([1 Lmax-Lmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:Lmax-Lmin+1)
set(gca, 'XTicklabel', Lmin:Lmax)
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
%ylabel('VAR')
%xlabel('L')

subplot(3, 3, 9);cla;hold on
   plot(BIQl(1, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'k')
   plot(BIQl(2, :), 'Linewidth', 2, 'Linestyle', '-', 'Color', 'b')
if(Lmax~=Lmin)
   axis([1 Lmax-Lmin+1 -34 4]);
   box;grid
end;
set(gca, 'XTick', 1:Lmax-Lmin+1)
set(gca, 'XTicklabel', Lmin:Lmax)
set(gca, 'YTick', -30:10:0)
%set(gca, 'YTicklabel', ['H'; 'G'; 'F'; 'E'])
%ylabel('BIQ')
%xlabel('L')

%legend(EInfo(1:esti, :))

print('-depsc2', filename);