Diffraction of Square Aperture
Contents
Fraunhoffer pattern
func = @(x,y) rect(x/2).*rect(y/2); [out inp] = fdiffract(func,0);
subplot(1,2,2);
imshow(out);
title(sprintf('maximum value %g\n',max(max(out))));
out = fdiffract(func,0.1);
title(sprintf('maximum value %g\n',max(max(out))));
Nf = 4
out = fdiffract(func,4);
Warning: Image is too big to fit on screen; displaying at 50%
Nf = 2
out = fdiffract(func,2);
Nf = 1;
out = fdiffract(func,1);