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