应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 凝聚态物理 » matlab模拟菲涅尔直边衍射
51/1返回列表
查看: 2967  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见

» 本帖已获得的红花（最新10朵）

banuo2005

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-04 00:34:27
已阅   同方向广播   申请物理EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
3楼2012-04-04 00:37:31
已阅   申请物理EPI   回复此楼   编辑   查看我的主页
查看全部 10 个回答

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
2楼2012-04-04 00:35:55
已阅   申请物理EPI   回复此楼   编辑   查看我的主页

smy1982

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
glddoyss: 金币+5, ★有帮助, 虽然没指出我的错误，不过提供的程序也有用啦 2012-04-04 15:24:08
利用光波衍射的菲涅尔近似－菲涅尔衍射公式来计算双缝的衍射：


下面的这个程序你明白了的话，你这个问题就解决了，其他物体的衍射模拟你也就会了。
（1）实验参数设定：
    Lamda=0.6328; % um, Wavelength of the input light
DDist=200.0e+3; % um, Diffraction distance to the output plane
M=512; N=512; % Number of the sample array
dx=10.0; dy=10.0; % um, Sampling space
SWidth=4; % pixels, Width of the Slits
SDist=30; % pixels, Distance between two slits
（2）建立衍射屏函数（Slit）的抽样点阵
  DScreen=zeros(M,N); % define a M*N matrix with zero elements;
DScreen(:,(N-SDist-SWidth)/2N-SDist+SWidth)/2)=1.0; % Slit one
DScreen(:,(N+SDist-SWidth)/2N+SDist+SWidth)/2)=1.0; % Slit two
%也可用下列语句调入实现设计好的衍射屏图样（bmp格式512×512）。
DScreen=double(imread('DScreen1.bmp'));
（3）衍射屏函数的图形显示：
figure(1); % 激活1号显示窗口为当前窗口
colormap(gray);  ％设置灰度图形显示格式
imshow(Slit); ％在当前图形窗口显示衍射屏函数
  （4）建立高斯球面波函数 的抽样点阵：
    % Set coordinate of the sample array
[x,y]=meshgrid(-M*dx/2:dxM-1)*dx/2,-N*dy/2:dyN-1)*dy/2);
SphFunct=exp(i*pi*(x.^2+y.^2)/(Lamda*DDist)); % The spherical wave function
  （5） 通过计算菲涅尔衍射积分求衍射光场：
    % Simulate the Fresnel diffraction
DScreen_F=fftshift(fft2(fftshift(DScreen))); %Fourier transform of the slit function.
SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function.
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform
  （6）显示输出平面上的光强度分布：
    % Display the diffraction pattern
FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution
FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I));
imshow(FrDiffract_I/max(max(FrDiffract_I)));
  （7）保存衍射图样
赞一下(1人) 回复此楼
如果你曾歌颂黎明，那么也请你拥抱黑夜。
4楼2012-04-04 09:53:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
5楼2012-04-04 12:52:11
已阅   申请物理EPI   回复此楼   编辑   查看我的主页
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +18 大米饭！ 2026-03-15 22/1100 2026-03-21 20:20 by HH领袖
[考研] 考研化学学硕调剂，一志愿985 +5 张vvvv 2026-03-15 7/350 2026-03-21 19:23 by ColorlessPI
[考研] 材料与化工（0856）304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[基金申请] 山东省面上项目限额评审 +3 石瑞0426 2026-03-19 3/150 2026-03-21 18:59 by RRRKKK
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 材料工程（专）一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
信息提示
关闭
请填处理意见
关闭