%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% Figure 4.15
%
% high pass impulse response plots
%
fnumber=2;
ratioP=1/.2;
apertureRange=-1:.01:1;
[x,y]=meshgrid(apertureRange,apertureRange);
figure(1);set(gcf,'color','white');
colormap(1-gray(256));
field=sqrt(x.^2+y.^2);
subplot(2,2,1);imagesc(apertureRange,apertureRange,P(field)-P(ratioP*field));
title('Transfer Function');axis square;
subplot(2,2,3);plot(apertureRange,P(apertureRange)-P(ratioP*apertureRange),'-k');
axis([-1 1 0 1.5]);
xlabel('1/wavelength');
fourierRange=-40:.1:40;
[x,y]=meshgrid(fourierRange,fourierRange);
field=sqrt(x.^2+y.^2);
subplot(2,2,2);imagesc(fourierRange,fourierRange,jinc(field./fnumber)-jinc(field./(ratioP*fnumber)));
title('impulse response');axis square;
subplot(2,2,4);plot(fourierRange,jinc(fourierRange./fnumber)-jinc(fourierRange./(ratioP*fnumber)),'-k');
xlabel('wavelengths');
axis([-40 40 -1.5 .5]);