load("result1.mat");
%ambient_dim = 100;
%sparsities = 1:100;
d = 32;%块个数
m = 2;%块内元素数
l1recovery = zeros(100,1);
for i = 1:32
    l1recovery(i) = theoretic(m,i,d);
end
m = 1:32;
n = 1:64;
figure;
contourf( m, n, result, 20, 'LineColor', 'none' );
hold on
cbar = colorbar;
myplot2 = plot( n, l1recovery, 'Color', 'w', 'Linewidth', 3 );
mylegend2 = legend(myplot2, '   Theory', 'Location', 'SouthEast' );
mylegend2.Box = 'off';
mylegend2.TextColor = [1,1,1];
mylegend2.FontSize = 12;

figure
plot(con1(:,1)/400,'*-')
xlabel("\fontname{Times New Roman} Block Sparsity \it K");
ylabel("\fontname{Times New Roman} Block Contribution Error");
savefig("contri_noise001.fig");
saveas(gca,"contri_noise001.eps");

figure
plot(con1(:,2)/400,'*-')
xlabel("\fontname{Times New Roman} Block Sparsity \it K");
ylabel("\fontname{Times New Roman} Block Contribution Error");
savefig("contri_noise01.fig");
saveas(gca,"contri_noise01.eps");

figure
plot(re1(:,1)/400,'*-')
xlabel("\fontname{Times New Roman} Block Sparsity \it K");
ylabel("\fontname{Times New Roman} Reconstruction Error");
savefig("re_noise001.fig");
saveas(gca,"re_noise001.eps");

figure
plot(re1(:,2)/400,'*-')
xlabel("\fontname{Times New Roman} Block Sparsity \it K");
ylabel("\fontname{Times New Roman} Reconstruction Error");
savefig("re_noise01.fig");
saveas(gca,"re_noise01.eps");