%本例完成夫瑯和費(fèi)衍射演示@wpp

clear;
clc;
close?all;
%a=imread（‘FNE.bmp‘）;%調(diào)入圖像
display（‘請選擇物函數(shù)類型....‘）;
display（‘（1）	矩孔‘）;
display（‘（2）	圓孔‘）;
display（‘（3）	圓盤‘）;
display（‘（4）	單縫‘）;
yourchoose=input（‘：‘）;
switch?yourchoose
????case?1
????????a=imread（‘REC.bmp‘）;
????????tt1=‘矩孔‘;
????case?2
????????a=imread（‘ROUND.bmp‘）;
????????tt1=‘圓孔‘;
????case?3
????????a=imread（‘ROUNDpan.bmp‘）;
????????tt1=‘圓盤‘;
???case?4
????????a=imread（‘SingleSlot.bmp‘）;
????????tt1=‘單縫‘;
????otherwise
????????quit;
end

a=double（a（::1））;
%?a=double（a）;
lamda=6328*10^（-10）;
k=2*pi/lamda;%波矢
figure;
imshow（a）;title（‘Original?Image‘）;
[cr]=size（a）;

%計(jì)算物光傳遞到透鏡的衍射過程

%L0=0.001;?%給出輸入圖像的尺寸L0單位:米
L0=0.1;
x0=linspace（-L0/2L0/2r）;%給出輸入圖像的坐標(biāo)
y0=linspace（-L0/2L0/2c）;
[x0y0]=meshgrid（x

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1555??2009-04-03?17:09??Fraunhofer\fhf2.m

?????文件?????786488??2008-06-09?13:45??Fraunhofer\REC.bmp

?????文件?????786488??2008-06-09?13:41??Fraunhofer\ROUND.bmp

?????文件?????786488??2008-06-09?15:24??Fraunhofer\ROUNDpan.bmp

?????文件?????786488??2008-06-09?14:09??Fraunhofer\Side.bmp

?????文件?????786488??2008-06-09?13:44??Fraunhofer\SingleSlot.bmp

?????目錄??????????0??2009-04-03?17:05??Fraunhofer

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

??????????????3933995????????????????????7