www.pudn.com > work.rar > wo.m

clear; 
cm=0.01;um=1e-6; 
N=128; 
rx=0*cm*ones(N);ry=0*cm*ones(N);rz=3*ones(N); 
ox=1*cm*ones(N);oy=1*cm*ones(N);oz=2*ones(N); 
hx=2*cm;hy=2*cm; 
dhx=hx/N;dhy=hy/N; 
w=0.63*um; 
k=2*pi/w; 
x=dhx.*(ones(N,1)*[-N/2:N/2-1]); 
y=(dhy.*(ones(N,1)*[-N/2:N/2-1]))'; 
z=0*ones(N); 
d1=sqrt((rx-x).^2+(ry-y).^2+(rz-z).^2); 
d2=sqrt((ox-x).^2+(oy-y).^2+(oz-z).^2); 
ur=exp(i*k*d1)./d1; 
uo=exp(i*k*d2)./d2; 
p=abs(ur+uo).^2; 
p=mat2gray(p); 
figure; 
imshow(p); 
p1=p.*conj(ur); 
clear;clc; 
cm=0.01;um=1e-6; 
N=128; 
hx=1*cm;hy=1*cm; 
dhx=hx/N;dhy=hy/N; 
x=dhx.*(ones(N,1)*[-N/2:N/2-1]); 
y=(dhy.*(ones(N,1)*[-N/2:N/2-1]))'; 
p1=imread('mz.jpg'); 
figure;imshow(p1); 
z=4; 
w=0.63*um; 
k=2*pi/w; 
[f1,dx1,dy1,x1,y1]=fresnell(p1,N,N,dhx,dhy,z,w); 
sit=pi/6; 
ur=exp(-j*k*y*sin(sit)); 
f1=f+ur; 
f1=abs(f1).^2; 
f1=mat2gray(f1); 
figure;imshow(f1); 
rp1=f1.*conj(ur); 
a=0.002;b=0.3; 
rp1=pure(rp1,a,b,x1,y1); 
[f2,dx2,dy2,x2,y2]=fresnell(p1,N,N,dhx,dhy,z,w); 
f2=abs(f2)^2; 
f2=mat2gray(f2); 
figure;imshow(f2);