tfpaper_stab.m

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

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

M0= 3;
L0= 2;
N0= 256;

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

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



exp= 'tfma';
EXP= 'TFMA';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Over N
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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);
   RESultname= sprintf('data/%04d/%s%d%d%s%s.mat', N, EXP, M0, L0, re_im, mo_no);
   load(resultname);
   STAB= zeros(size(PAR));
   for mm= 1:20
      for ee= 1:esti
         
      end;
      
      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);
   sprintf('%s  %d', resultname, MM)
   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);
   sprintf('%s  %d', resultname , MM)
   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', 'k')
   plot(MSEn(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(MSEn(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEn(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEn(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(VARn(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(VARn(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(VARn(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(VARn(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(BIQn(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(BIQn(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(BIQn(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(BIQn(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(MSEm(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(MSEm(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEm(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEm(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(VARm(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(VARm(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(VARm(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(VARm(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(BIQm(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(BIQm(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(BIQm(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(BIQm(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(MSEl(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(MSEl(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEl(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEl(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(VARl(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(VARl(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(VARl(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(VARl(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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', 'k')
   plot(BIQl(3, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(BIQl(4, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(BIQl(5, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(BIQl(6, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
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);












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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot MSE/VAR/BIQ over N: