%
%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% Figure 10.30
%
% system transfer function for a multiple aperture imager
%
% units are  inverse Delta
%
M=1;
uS=4;
snr=10000;
urange=linspace(0,.3,100);
figure(1);set(gcf,'color','white');
stf1=abs(sinc(uS*urange/M).*otf(urange/M));
subplot(3,1,1);plot(urange,[abs(sinc(uS*urange/M));otf(urange/M);stf1],'-k');
title('(a)');axis([0 .3 0 1]);
line([1/(2*uS) 1/(2*uS)],[0 1],'Color','k');
M=1/2;
figure(1);set(gcf,'color','white');
stf3=abs(sinc(uS*urange/M).*otf(urange/M));
subplot(3,1,2);plot(urange,[abs(sinc(uS*urange/M));otf(urange/M);stf3],'-k');
title('(b)');axis([0 .3 0 1]);
line([1/(2*uS) 1/(2*uS)],[0 1],'Color','k');
nf=(1+snr*stf1.^2)./(1+snr*stf3.^2);
subplot(3,1,3);plot(urange,nf,'-k');
xlabel('u');title('(c)');axis([0 .3 1 100]);
line([1/(2*uS) 1/(2*uS)],[0 100],'Color','k');