clear
%程序基本參數單位m
lamta=632.8e-9;
k=2*pi/lamta;
dx=10e-5;
dy=10e-5;
xd=512;
yd=512;
z=0.05;
c1=2;
c2=2;
[mn]=meshgrid（linspace（-xd/2xd/2-1xd）linspace（-yd/2yd/2-1yd））;
dfx=1/（xd*dx）;
dfy=1/（xd*dy）;
I64=zeros（xdyd）;
I641=zeros（xdyd）;
xxd1=xd/8;
yyd1=xd/4;
r=xxd1/2;
a=xxd1;
D1=（（m-a）.^2+（n）.^2）.^（1/2）;
D2=（（m+a）.^2+（n）.^2）.^（1/2）;
i1=find（（D1<=r|D2<=r））;
I641（i1）=1;
I642=zeros（xdyd）;
B1=ones（xxd1yyd1）;
I642（[（xd/2-xxd1/2+1）:（xd/2+xxd1/2）][（yd/2-yyd1/2+1）:（yd/2+yyd1/2）]）=B1;
I64=I641+I642;
i2=find（I64>=1）;
I64（i2）=1;
figure（1）
imagesc（I64）
colormap（gray）
axis?square
axis?off


%離散光波表達式
R=zeros（xdyd）;
switch?c1
????case?1
????????R=exp（-j*k*（（dx*m）.^2+（dy*n）.^2）/2*z）;%球面波
????case?2
????????%%%用平行光照射irradiated?by?parallel?light
????????the

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2668??2003-12-18?19:33??raceway.m

-----------?---------??----------?-----??----

?????????????????2668????????????????????1