tfarma_fit.m

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

set(gca, 'XTick', 1:Lmax+1)
set(gca, 'XTickLabel', 0:Lmax)

mean= [
    sum((1:Mmax)'.*MDLM)/MM;
    sum((1:Mmax)'.*BICM)/MM;
    sum((1:Mmax)'.*AICM)/MM;
    sum((1:Mmax)'.*GICM)/MM;
    sum((1:Mmax)'.*MICM)/MM;
    sum((1:Mmax)'.*mdlM)/MM;
    sum((1:Mmax)'.*bicM)/MM;
    sum((1:Mmax)'.*aicM)/MM;
    sum((1:Mmax)'.*gicM)/MM;
    sum((1:Mmax)'.*micM)/MM;
    ];
meanMm= [meanMm mean];

mean= [
    sum((0:Lmax).*MDLL)/MM;
    sum((0:Lmax).*BICL)/MM;
    sum((0:Lmax).*AICL)/MM;
    sum((0:Lmax).*GICL)/MM;
    sum((0:Lmax).*MICL)/MM;
    sum((0:Lmax).*mdlL)/MM;
    sum((0:Lmax).*bicL)/MM;
    sum((0:Lmax).*aicL)/MM;
    sum((0:Lmax).*gicL)/MM;
    sum((0:Lmax).*micL)/MM;
    ];
meanLm= [meanLm mean];

   mse= [
       10*log10((sum(((1:Mmax)'-M).^2.*MDLM) + sum(((0:Lmax)-L).^2.*MDLL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*BICM) + sum(((0:Lmax)-L).^2.*BICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*AICM) + sum(((0:Lmax)-L).^2.*AICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*GICM) + sum(((0:Lmax)-L).^2.*GICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*MICM) + sum(((0:Lmax)-L).^2.*MICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*mdlM) + sum(((0:Lmax)-L).^2.*mdlL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*bicM) + sum(((0:Lmax)-L).^2.*bicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*aicM) + sum(((0:Lmax)-L).^2.*aicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*gicM) + sum(((0:Lmax)-L).^2.*gicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*micM) + sum(((0:Lmax)-L).^2.*micL))/MM/(M^2+L^2));
       ];
   MSEm= [MSEm mse];
   drawnow
end;

N= 256;
M= 3;
MSEl= [];
meanMl= [];
meanLl= [];
for L= 0:3
   exp= 'ord';
   resultname= sprintf('data/%04d/%s%d%d%s%s.mat', N, exp, M, L, re_im, mo_no)
   load(resultname)

   MDLM= sum(MDLhist.').';MDLL= sum(MDLhist);
   BICM= sum(BIChist.').';BICL= sum(BIChist);
   AICM= sum(AIChist.').';AICL= sum(AIChist);
   GICM= sum(GIChist.').';GICL= sum(GIChist);
   MICM= sum(MIChist.').';MICL= sum(MIChist);
   mdlM= sum(mdlhist.').';mdlL= sum(mdlhist);
   bicM= sum(bichist.').';bicL= sum(bichist);
   aicM= sum(aichist.').';aicL= sum(aichist);
   gicM= sum(gichist.').';gicL= sum(gichist);
   micM= sum(michist.').';micL= sum(michist);

figure(L+1);clf;
subplot(4, 2, 1);bar(MDLM, 1);axis([.5 Mmax+0.5 0 MM])
subplot(4, 2, 2);bar(MDLL, 1);axis([.5 Lmax+1.5 0 MM])
set(gca, 'XTick', 1:Lmax+1)
set(gca, 'XTickLabel', 0:Lmax)

subplot(4, 2, 3);bar(BICM, 1);axis([.5 Mmax+0.5 0 MM])
subplot(4, 2, 4);bar(BICL, 1);axis([.5 Lmax+1.5 0 MM])
set(gca, 'XTick', 1:Lmax+1)
set(gca, 'XTickLabel', 0:Lmax)

subplot(4, 2, 5);bar(AICM, 1);axis([.5 Mmax+0.5 0 MM])
subplot(4, 2, 6);bar(AICL, 1);axis([.5 Lmax+1.5 0 MM])
set(gca, 'XTick', 1:Lmax+1)
set(gca, 'XTickLabel', 0:Lmax)

subplot(4, 2, 7);bar(GICM, 1);axis([.5 Mmax+0.5 0 MM])
subplot(4, 2, 8);bar(GICL, 1);axis([.5 Lmax+1.5 0 MM])
set(gca, 'XTick', 1:Lmax+1)
set(gca, 'XTickLabel', 0:Lmax)

mean= [
    sum((1:Mmax)'.*MDLM)/MM;
    sum((1:Mmax)'.*BICM)/MM;
    sum((1:Mmax)'.*AICM)/MM;
    sum((1:Mmax)'.*GICM)/MM;
    sum((1:Mmax)'.*MICM)/MM;
    sum((1:Mmax)'.*mdlM)/MM;
    sum((1:Mmax)'.*bicM)/MM;
    sum((1:Mmax)'.*aicM)/MM;
    sum((1:Mmax)'.*gicM)/MM;
    sum((1:Mmax)'.*micM)/MM;
    ];
meanMl= [meanMl mean];

mean= [
    sum((0:Lmax).*MDLL)/MM;
    sum((0:Lmax).*BICL)/MM;
    sum((0:Lmax).*AICL)/MM;
    sum((0:Lmax).*GICL)/MM;
    sum((0:Lmax).*MICL)/MM;
    sum((0:Lmax).*mdlL)/MM;
    sum((0:Lmax).*bicL)/MM;
    sum((0:Lmax).*aicL)/MM;
    sum((0:Lmax).*gicL)/MM;
    sum((0:Lmax).*micL)/MM;
    ];
meanLl= [meanLl mean];

   mse= [
       10*log10((sum(((1:Mmax)'-M).^2.*MDLM) + sum(((0:Lmax)-L).^2.*MDLL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*BICM) + sum(((0:Lmax)-L).^2.*BICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*AICM) + sum(((0:Lmax)-L).^2.*AICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*GICM) + sum(((0:Lmax)-L).^2.*GICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*MICM) + sum(((0:Lmax)-L).^2.*MICL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*mdlM) + sum(((0:Lmax)-L).^2.*mdlL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*bicM) + sum(((0:Lmax)-L).^2.*bicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*aicM) + sum(((0:Lmax)-L).^2.*aicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*gicM) + sum(((0:Lmax)-L).^2.*gicL))/MM/(M^2+L^2));
       10*log10((sum(((1:Mmax)'-M).^2.*micM) + sum(((0:Lmax)-L).^2.*micL))/MM/(M^2+L^2));
       ];
   MSEl= [MSEl mse];
   drawnow
end;


test= 1;

figure(99);clf
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', 2, 'Linestyle', '--', 'Color', 'k')
   plot(MSEn(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(MSEn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEn(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(MSEn(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(MSEn(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEn(10,:), 'Linewidth', 1, 'Linestyle', '.' , 'Color', 'r')
end;
   axis([1 4 -10 -2]);box;grid
   set(gca, 'XTick', 1:4)
   set(gca, 'XTicklabel', 2.^(6:9))
   set(gca, 'YTick', -10:2:-1)
   ylabel('MSE')
   xlabel('N')

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', 2, 'Linestyle', '--', 'Color', 'k')
   plot(MSEm(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(MSEm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEm(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(MSEm(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(MSEm(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEm(10,:), 'Linewidth', 1, 'Linestyle', '.' , 'Color', 'r')
end;
   axis([1 4 -10 -2]);box;grid
   set(gca, 'XTicklabel', 2:5)
   set(gca, 'YTicklabel', ['AA'; 'BB'; 'CC'; 'DD'; 'EE'])
   ylabel('MSE')
   xlabel('M')

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', 2, 'Linestyle', '--', 'Color', 'k')
%   plot(MSEl(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(MSEl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(MSEl(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(MSEl(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(MSEl(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(MSEl(10,:), 'Linewidth', 1, 'Linestyle', '.' , 'Color', 'r')
end;
   axis([1 4 -10 -2]);box;grid
   set(gca, 'XTicklabel', 0:3)
   ylabel('MSE')
   xlabel('L')

subplot(3, 3, 4);cla;hold on
   plot(meanMn(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanMn(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanMn(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanMn(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
   plot(meanMn(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(meanMn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanMn(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanMn(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanMn(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanMn(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot([3 3 3 3], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
end;
   axis([1 4 1 5]);box;grid
   set(gca, 'XTicklabel', 2.^(6:9))
   ylabel('Mbar')
   xlabel('N')

subplot(3, 3, 5);cla;hold on
   plot(meanMm(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanMm(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanMm(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanMm(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
%   plot(meanMm(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(meanMm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanMm(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanMm(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanMm(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanMm(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot([2 3 4 5], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
end;
   axis([1 4 1 5]);box;grid
   set(gca, 'XTicklabel', 2:5)
   ylabel('Mbar')
   xlabel('M')

subplot(3, 3, 6);cla;hold on
   plot(meanMl(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanMl(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanMl(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanMl(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
%   plot(meanMl(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(meanMl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanMl(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanMl(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanMl(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanMl(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
end;
   plot([3 3 3 3], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
   axis([1 4 1 5]);box;grid
   set(gca, 'XTicklabel', 0:3)
   ylabel('Mbar')
   xlabel('L')

subplot(3, 3, 7);cla;hold on
   plot(meanLn(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanLn(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanLn(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanLn(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
   plot(meanLn(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(meanLn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanLn(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanLn(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanLn(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanLn(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
end;
   plot([2 2 2 2], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
   axis([1 4 0 4]);box;grid
   set(gca, 'XTicklabel', 2.^(6:9))
   ylabel('Lbar')
   xlabel('N')

subplot(3, 3, 8);cla;hold on
   plot(meanLm(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanLm(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanLm(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanLm(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
%   plot(meanLm(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
   plot([2 2 2 2], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
if(test)
   plot(meanLm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanLm(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanLm(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanLm(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanLm(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
end;
   axis([1 4 0 4]);box;grid
   set(gca, 'XTicklabel', 2:5)
   ylabel('Lbar')
   xlabel('M')

subplot(3, 3, 9);cla;hold on
   plot(meanLl(1, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')
   plot(meanLl(2, :), 'Linewidth', 2, 'Linestyle', ':' , 'Color', 'k')
   plot(meanLl(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')
   plot(meanLl(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')
%   plot(meanLl(5, :), 'Linewidth', 2, 'Linestyle', '.' , 'Color', 'k')
if(test)
   plot(meanLl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')
   plot(meanLl(7, :), 'Linewidth', 1, 'Linestyle', ':' , 'Color', 'r')
   plot(meanLl(8, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')
   plot(meanLl(9, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot(meanLl(10,:), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')
   plot([0 1 2 3], 'Linewidth', 1, 'Linestyle', '-', 'Color', 'k')
end;
   axis([1 4 0 4]);box;grid
   set(gca, 'XTicklabel', 0:3)
   ylabel('Lbar')
   xlabel('L')

   print -deps2 -tiff -r300 /users/mjachan/tex/prj/TF/figures/ordMSE.eps

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%