?? tfpaper_stab.m
字號:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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= []; 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); 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= .55subplot(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', 1, 'Linestyle', '--', '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 -17 2]); 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', '--', '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', '--', '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 -34 4]); 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', '--', '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', '--', '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 -34 4]); 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', '--', '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', '--', '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 -17 2]); 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', '--', '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', '--', '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 -34 4]); 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', '--', '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', '--', '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 -34 4]); 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', '--', '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', '--', '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 -17 2]); 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', '--', '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', '--', '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 -34 4]); 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', '--', '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', '--', '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 -34 4]); 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);M0= 2;L0= 2;N0= 256;nmin= 6;Mmin= 1;Lmin= 0;nmax= 9;Mmax= 4;Lmax= 4;?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -