基于matlab设计DOE元件的GS算法源代码.rar

%%IFTA FOR DOE DESIGN%% clear; %target image load N=500; target=imread('C:\DOE\test.bmp'); target=(target(:,:,1)+target(:,:,2)+target(:,:,3))/3; m=size(target,1); scale=N/m; target=imresize(target,scale); target=double(target); target=target/max(max(target)); %DOE phase estimate from random matrix DOE=2*pi*rand(N,N); s=5; %IFTA algorithm for t=1:s; DOEphase=exp(1i*DOE); %forward iteration iterf=fft2(DOEphase); intf=abs(iterf); angf=angle(iterf); A=target.*exp(1i*angf); %backward iteration iterb=ifft2(A); angb=angle(iterb); DOE=angb; error=target-intf/max(max(intf)); E=sum(sum(abs(error)))/(N*N); if E<0.05; iteration=t; break end end figure(1) colormap(gray); imagesc(target) figure(2) doef=abs((fft2(DOEphase))); colormap(gray); imagesc(doef) figure(3) colormap(gray); angf=angle(DOEphase); m11=max(max(angf)) m12=min(min(angf)) imagesc(angf)