% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% figure 8.9
%
figure(1);
set(gcf,'color','white');
colormap gray;
Hl=toeplitz([1, zeros(1,127)],[1 1 1 1 0 0 0 0,zeros(1,128-8)])/4;
[al bl cl]=svd(Hl);
subplot(2,2,1);plot(diag(bl),'-k');axis([0 128 0 1]);title('Singular value spectra');
hold on;
subplot(2,2,2);imagesc(al);axis off;axis image;title('11110000');
Hl=toeplitz([1, zeros(1,127)],[1 1 0 0 0 1 0 1,zeros(1,128-8)])/4;
[al bl cl]=svd(Hl);
subplot(2,2,1);plot(diag(bl),'-k');
subplot(2,2,3);imagesc(al);axis off;axis image;title('11000101');
Hl=toeplitz([1, zeros(1,127)],[1 1 1 0 0 1 0 0,zeros(1,128-8)])/4;
[al bl cl]=svd(Hl);
subplot(2,2,1);plot(diag(bl),'-k');
subplot(2,2,4);imagesc(al);axis off;axis image;title('11100100');