tfma_papera.m

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

      MBV= ['MSE'; 'BIQ'; 'VAR'];
      figure(2*(WIN_L-1)+WIN_M);clf;hold on
      for reim= 2:2
         switch(reim)
          case 1
           re_im= 'r';
          case 2
           re_im= 'i';
         end;
         for mbv= 1:3
            titlestring= sprintf('%s for TFMA(%d, %d; N), %s, %d %d', ...
                                 MBV(mbv, :), MMA, LMA, re_im, WIN_M, WIN_L);
            switch(mbv)
             case 1
              MBVNml= MSENml;
              MBVNmn= MSENmn;
             case 2
              MBVNml= BIQNml;
              MBVNmn= BIQNmn;
             case 3
              MBVNml= VARNml;
              MBVNmn= VARNmn;
             end;
            subplot(2, 3, mbv+3*(reim-1));hold on;
            subplot(1, 3, mbv);hold on;
            plot(reshape(MBVNmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
%            plot(reshape(MBVNmn(2, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
            plot(reshape(MBVNmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
      
            axis([1 nmax-nmin+1 0 .3])
%            axis([1 nmax-nmin+1 0 max(max(max(max([MBVNml MBVNmn]))))])
            set(gca, 'XTick', 1:nmax-nmin+1)
            set(gca, 'XTickLabel', 2.^(nmin:nmax))
%            title(titlestring)
	 end;
      end;
   end;
end;

s= input('next:', 's');
switch s
 case 'm'
  MMA= MMA-1;
  if(MMA<0)
     MMA= 0;
  end;
 case 'M'
  MMA= MMA+1;
  if(MMA>5)
     MMA= 5;
  end;
 case 'l'
  LMA= LMA-1;
  if(LMA<0)
     LMA= 0;
  end;
 case 'L'
  LMA= LMA+1;
  if(LMA>3)
     LMA= 3;
  end;
end;


end;



figure(1);clf;hold on
plot(reshape(MSENmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(reshape(MSENmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(BIQNmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 nmax-nmin+1 0 .4])
set(gca, 'XTick', 1:nmax-nmin+1)
set(gca, 'XTickLabel', 2.^(nmin:nmax))
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])


figure(2);clf;hold on
plot(10*log10(reshape(MSENmn(1, reim, nmin:nmax), nmax-nmin+1, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(10*log10(reshape(MSENmn(3, reim, nmin:nmax), nmax-nmin+1, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(BIQNmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, reim, nmin:nmax), nmax-nmin+1, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 nmax-nmin+1 -20 -4])
set(gca, 'XTick', 1:nmax-nmin+1)
set(gca, 'XTickLabel', 2.^(nmin:nmax))
grid;box
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])









%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   0;
LAR  =   0;
mo_no= 'n';
alpha= 1/2;
beta = 1/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

nmin= 8;
nmax= 8;

MMA= 4;
WIN_M= 1;
WIN_L= 2;
MSENml= zeros(3,   2,  4);BIQNml= MSENml;VARNml= MSENml;
MSENmn= zeros(3,   2,  4);BIQNmn= MSENmn;VARNmn= MSENmn;
N= 2^8;
for LMA= 0:3
   resultnamei= sprintf('data/%04d/tfma_ceps_papera%d%d%d%d%s.mat', ...
                         N, MMA, LMA, WIN_M, WIN_L, 'i')
   load(resultnamei)
   MSENml(:, 2, LMA+1)= MSEml;
   BIQNml(:, 2, LMA+1)= BIQml;
   VARNml(:, 2, LMA+1)= VARml;
   MSENmn(:, 2, LMA+1)= MSEmn;
   BIQNmn(:, 2, LMA+1)= BIQmn;
   VARNmn(:, 2, LMA+1)= VARmn;
end;
      
figure(3);clf;hold on
plot(reshape(MSENmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(reshape(MSENmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')

%plot(reshape(BIQNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 4 0 .2])
set(gca, 'XTick', 1:4)
set(gca, 'XTickLabel', 0:3)
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])



figure(4);clf;hold on
plot(10*log10(reshape(MSENmn(1, 2, :), 4, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(10*log10(reshape(MSENmn(3, 2, :), 4, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(BIQNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 4 -20 -4])
set(gca, 'XTick', 1:4)
set(gca, 'XTickLabel', 0:3)
grid;box
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])






tfpm;
nmin= 8;
nmax= 8;

LMA= 2;
WIN_M= 1;
WIN_L= 2;
MSENml= zeros(3,   2,  6);BIQNml= MSENml;VARNml= MSENml;
MSENmn= zeros(3,   2,  6);BIQNmn= MSENmn;VARNmn= MSENmn;
N= 2^8;
for MMA= 0:5
   resultnamei= sprintf('data/%04d/tfma_ceps_papera%d%d%d%d%s.mat', ...
                         N, MMA, LMA, WIN_M, WIN_L, 'i')
   load(resultnamei)
   MSENml(:, 2, MMA+1)= MSEml;
   BIQNml(:, 2, MMA+1)= BIQml;
   VARNml(:, 2, MMA+1)= VARml;
   MSENmn(:, 2, MMA+1)= MSEmn;
   BIQNmn(:, 2, MMA+1)= BIQmn;
   VARNmn(:, 2, MMA+1)= VARmn;
end;
      
figure(3);clf;hold on
plot(reshape(MSENmn(1, 2, :), 6, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(reshape(MSENmn(3, 2, :), 6, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')

%plot(reshape(BIQNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 4 0 .2])
set(gca, 'XTick', 1:4)
set(gca, 'XTickLabel', 0:3)
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])



figure(4);clf;hold on
plot(10*log10(reshape(MSENmn(1, 2, :), 6, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-')
plot(10*log10(reshape(MSENmn(3, 2, :), 6, 1)), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(BIQNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')
%plot(reshape(BIQNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '--')

%plot(reshape(VARNmn(1, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
%plot(reshape(VARNmn(3, 2, :), 4, 1), 'Color', 'k', 'Linewidth', 2, 'Linestyle', '-.')
axis([1 6 -20 -4])
set(gca, 'XTick', 1:6)
set(gca, 'XTickLabel', 0:5)
grid;box
set(gca, 'XTickLabel', [])
set(gca, 'YTickLabel', [])