exp= 'tfarma';filename= sprintf('/users/mjachan/tex/prj/TF/figures/tfarma%02d%02d%02d%1s%1s', M0, L0, log2(N0), re_im, mo_no)figure(fig3);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);   sprintf('%s  %d', resultname, MM)      MSEnee= [];   VARnee= [];   BIQnee= [];   MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   VARnee= [VARnee; sum(sum(VAR(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   for ee= 2: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);   sprintf('%s  %d', resultname, MM)   MSEnee= [];   VARnee= [];   BIQnee= [];   MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   VARnee= [VARnee; sum(sum(VAR(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   for ee= 2: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= [];   MSEnee= [MSEnee; sum(sum(MSE(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   VARnee= [VARnee; sum(sum(VAR(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   BIQnee= [BIQnee; sum(sum(BIQ(:, :, 1, 1)))/sum(sum(abs(par(:, 1:M0)).^2))];   for ee= 2: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= .55subplot(3, 3, 1);cla;hold on   plot(MSEn(1, :), 'Linewidth', 2, 'Linestyle', 'o' , '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', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(MSEn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(MSEn(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(MSEn(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(nmax~=nmin)%   axis([1 nmax-nmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(VARn(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(VARn(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(VARn(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(VARn(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(VARn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(VARn(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(VARn(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(nmax~=nmin)%   axis([1 nmax-nmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(BIQn(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(BIQn(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(BIQn(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(BIQn(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(BIQn(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(BIQn(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(BIQn(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(nmax~=nmin)%   axis([1 nmax-nmin+1 -10 10]);   box;gridend;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', 'o' , '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', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(MSEm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(MSEm(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(MSEm(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Mmax~=Mmin)%   axis([1 Mmax-Mmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(VARm(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(VARm(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(VARm(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(VARm(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(VARm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(VARm(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(VARm(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Mmax~=Mmin)%   axis([1 Mmax-Mmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(BIQm(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(BIQm(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(BIQm(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(BIQm(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(BIQm(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(BIQm(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(BIQm(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Mmax~=Mmin)%   axis([1 Mmax-Mmin+1 -10 10]);   box;gridend;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', 'o' , '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', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(MSEl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(MSEl(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(MSEl(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Lmax~=Lmin)%   axis([1 Lmax-Lmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(VARl(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(VARl(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(VARl(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(VARl(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(VARl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(VARl(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(VARl(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Lmax~=Lmin)%   axis([1 Lmax-Lmin+1 -10 10]);   box;gridend;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', 'o' , 'Color', 'k')   plot(BIQl(2, :), 'Linewidth', 2, 'Linestyle', '-' , 'Color', 'k')   plot(BIQl(3, :), 'Linewidth', 2, 'Linestyle', '-.', 'Color', 'k')   plot(BIQl(4, :), 'Linewidth', 2, 'Linestyle', '--', 'Color', 'k')   plot(BIQl(5, :), 'Linewidth', 1, 'Linestyle', 'o' , 'Color', 'r')   plot(BIQl(6, :), 'Linewidth', 1, 'Linestyle', '-' , 'Color', 'r')   plot(BIQl(7, :), 'Linewidth', 1, 'Linestyle', '-.', 'Color', 'r')   plot(BIQl(8, :), 'Linewidth', 1, 'Linestyle', '--', 'Color', 'r')if(Lmax~=Lmin)%   axis([1 Lmax-Lmin+1 -10 10]);   box;gridend;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);