MATLAB模拟夫琅禾费衍射

这个是夫琅禾费衍射模拟， 包括圆孔，矩孔，圆盘，单缝， 在WinXP下MATLAB2008通过，

%本例完成夫琅和费衍射演示@wpp

clear;
clc;
close all;
%a=imread（‘FNE.bmp‘）;%调入图像
display（‘请选择物函数类型....‘）;
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）;

%计算物光传递到透镜的衍射过程

%L0=0.001; %给出输入图像的尺寸L0单位:米
L0=0.1;
x0=linspace（-L0/2L0/2r）;%给出输入图像的坐标
y0=linspace（-L0/2L0/2c）;
[x0y0]=meshgrid（x

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       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