r=512c=r;?????????????????????????%物面采樣數(shù)
Uo=zeros（cr）;??????????????????????%預設物
d=30;a=10;????????????????????????%光柵常數(shù)和縫寬
for?n=1:d:c????????????????????????%循環(huán)生成物（二維光柵）
???Uo（n:n+a:）=1;
end
for?m=1:d:r
???Uo?（:m:m+a）=1;
end
Uo=Uo（1:c1:r）;
figure
subplot（221）
imshow（Uo[]）?????????????????%顯示物分布
lamda=6328*10^（-10）;k=2*pi/lamda;????%賦值波長單位:米波矢
f=0.004;?Lo=0.001???????????????????%賦值透鏡的焦距物面的尺寸Lo單位:米
D1=0.00005????????????????????????%賦值濾波片直徑單位:米
D2=0.00005????????????????????????%賦值濾波片寬度單位:米
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面計算物光傳遞到透鏡的衍射過程（S-FFT）
xo=linspace（-Lo/2Lo/2r）;yo=linspace（-Lo/2Lo/2c）;?%賦值物面的坐標
[xoyo]=meshgrid（xoyo）;?????????????%生成物面的坐標網(wǎng)格
do=0.0041;?????????????????????????%物面到透鏡的距離do單位:米
L=r*lamda*do/Lo????????????????????%衍射光在透鏡前表面上的尺寸L單位:米
xl=linspace（-L/2L/2r）;yl=linspace（-L/2L/2c）;?%賦值透鏡前表面的坐標
[xlyl]=meshgrid（xlyl）;???????????????%生成透鏡前表面的坐標網(wǎng)格
F0=exp（j*k*do）/（j*lamda*do）*exp（j*k/2/do*（xl.^2+yl.^2））;
F=exp（j*k/2/do*（xo.^2+yo.^2））;
FU=（Lo*Lo/r/r）.*fftshift（fft2（Uo.*F））;?
U1=F0.*FU;?????????????????????????%透鏡前表面上的光場復振幅分布
I1=U1.*conj（U1）;?????????????????????%透鏡前表面上的光強分布
subplot（222）
imshow（I1[]）?colormap（pink）title（‘透鏡上的光強分布‘）
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面計算通過透鏡后的光場
U1yp=U1.*exp（-j*k.*（xl.^2+yl.^2）/2/f）;???%計算通過透鏡后的光場
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面計算通過透鏡后的光場到達后焦面的過程（S-FFT）
dlf=f
Lyp=r*lamda*dlf/L???????????????????%給出后焦面的尺寸單位:米
xf=linspace（-Lyp/2Lyp/2r）;yf=linspace（-Lyp/2Lyp/2c）;?%給出后焦面的坐標
[xfyf]=meshgrid（xfyf）;???????????????%生成后焦面的坐標網(wǎng)格
F0=exp（j*k*dlf）/（j*lamda*dlf）*exp（j*k/2/dlf*（xf.^2+yf.^2））;
F=exp（j*k/2/dlf*（xl.^2+yl.^2））;
Uf=（L*L/r/r）.*fft2（U1yp.*F）;Uf=Uf.*F0;??%?計算后焦面上的光場分布
I2=Uf.*conj（Uf）;?????????????????????%?后焦面上的光強分布
subplot（223）imshow（I2[0max（I2（:））/100]）?colormap（pink）title（‘后焦面上的光強分布‘）
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面生成三種濾波器
DD=round（D1*r/Lyp）;?????????????????%賦值濾波器直徑，單位：像素
SD=round（D2*r/Lyp/2）;????????????????%賦值濾波器寬帶，單位：像素
H1=zeros（cr）;????????????????????????%預設濾波器H1
for?n=1:c????????????????????????????%循環(huán)生成濾波器H1
???for?m=1:r
??????if?（n-c/2-1）.^2+（m-r/2-1）.^2<=（DD/2）.^2;
??????H1（nm）=1;
??????end
???end
end
figuresubplot（131）imshow（H1[]）;title（‘濾波器H1‘）
H2=zeros（cr）;????????????????????????%預設濾波器H2
H2（round（c/2）-SD:round（c/2）+SD:）=1;????%生成濾波器H2
subplot（132）imshow（H2[]）;title（‘濾波器H2‘）
H3=zeros（cr）;????????????????????????%預設濾波器H3
H3（:round（r/2）-SD:round（r/2）+SD）=1;????%生成濾波器H3
subplot（133）imshow（H3[]）;title（‘濾波器H3‘）
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%完成三個濾波，得到三個不同的光場
Uf1=H1.*Uf;
Uf2=H2.*Uf;
Uf3=H3.*Uf;
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面計算通過焦面后的光場到達像面的衍射成像過程（S-FFT）
dfi=do*f/（do-f）-f;?????????????????????%滿足物像公式
Li=r*lamda*dfi/Lyp???

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

?????文件??????27838??2016-05-06?21:20??Spatial?filtering\1.jpg

?????文件???????9489??2016-05-06?21:19??Spatial?filtering\2.jpg

?????文件??????27859??2016-05-06?21:19??Spatial?filtering\3.jpg

?????文件??????29129??2016-05-06?21:19??Spatial?filtering\4.jpg

?????文件???????4586??2016-04-19?16:35??Spatial?filtering\CX9_1.m

?????文件???????3187??2016-04-19?21:49??Spatial?filtering\CX9_1_1.m

?????目錄??????????0??2016-05-06?21:24??Spatial?filtering

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

???????????????102088????????????????????7