%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% Figure 10.31
%
%
% system transfer function for a multiple aperture imager
%
% units are  inverse Delta
%
M=1;
uS=4;
snr=10000;
urange=linspace(0,.25,1000);
figure(10);set(gcf,'color','white');
stf1=abs(sinc(uS*urange/M).*otf(urange/M));
mse1=1./(1+stf1(1:500).^2*snr./(1+snr*stf1(501:1000)));
M=1/2;
stf3=abs(sinc(uS*urange/M).*otf(urange/M));
mse3=1./(1+stf3(1:500).^2*snr./(1+snr*stf3(501:1000)));
plot(urange(1:500),[mse1;mse3],'-k'); xlabel('u');
axis([0 0.125 0 1]);