傅里叶叠层成像matlab程序

clc;clear all;
%% simulate the forward imaging process of Fourier ptychography
%% simulate the high resoluiton complex object
objAmp = double（imread（‘cameraman.tif‘））;
phase = double（imread（‘westconcordorthophoto.png‘））;
phase = pi*imresize（phase[256 256]）./max（max（phase））;
obj = objAmp.*exp（1i.* phase）;
imshow（abs（obj）[]）;title（‘Input complex object‘）;
F = @（x） fftshift（fft2（x））;
invF = @（x） ifft2（ifftshift（x））;

%% setup the parameters for the coherent imaging system
wLength = 0.63e-6;
k0 = 2*pi/wLength;
spsize = 2.75e-6; % sampling pixel size of the CCD
psize = spsize / 4; % final pixel size of the reconstruction
NA = 0.08;

%% create the wave vectors for the LED illumiantion 
arraysize = 15;
kx = k0.* repmat（linspace（0.3-0.3arraysize）1arraysize）;
ky = k0.* re