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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 凝聚态物理 » matlab模拟菲涅尔直边衍射
101/1返回列表
查看: 2935  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

匿名

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

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

banuo2005

» 猜你喜欢

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

匿名

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

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
3楼2012-04-04 00:37:31
已阅   申请物理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   回复此楼   编辑   查看我的主页

smy1982

木虫 (正式写手)
引用回帖:
5楼: Originally posted by glddoyss at 2012-04-04 12:52:11:
这个不很懂啊。学matlab不久。。  能帮忙看看我函数定义哪里出了问题吗？
M文件这个
function CS=fresnelC(a)
CS=quad(@fun1,0,a);
function fun=fun1(u)
fun=cos(0.5.*pi.*u.^2)+i.*sin(0.5.*pi.*u.^2);
...

你这个程序不是用run来运行的，你的是个函数。得用函数命令。
在matlab命令行里直接输入比如：CS=fresnelC(200e-3)；
但是你这个程序里面，“l(n)=abs(F(x1(n))+sqrt(0.5)*exp(i.*pi./4));%表示的是菲涅尔积分，f(n)=F(x1(n));%得到菲涅尔积分的列阵
” 这个F（）是什么函数啊？没有给出来啊？
赞一下 回复此楼
如果你曾歌颂黎明，那么也请你拥抱黑夜。
6楼2012-04-04 20:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhfeng

木虫 (小有名气)
glddoyss: 回帖置顶 2012-04-06 19:28:32
对于quad函数，积分限应该是常数，a必须赋值。可以看一下帮助文件
赞一下(1人) 回复此楼
7楼2012-04-04 23:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

banuo2005
送鲜花一朵
8楼2012-04-05 20:47:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
9楼2012-04-06 19:28:49
已阅   申请物理EPI   回复此楼   编辑   查看我的主页

卡洛比

新虫 (初入文坛)
??? Error using ==> times
Matrix dimensions must agree.

Error in ==> Fresnelyanshe at 25
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform傅里叶逆变换（球面函数及狭缝函数的乘积）

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> Fresnelyanshe at 25
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform傅里叶逆变换（球面函数及狭缝函数的乘积）
同学 你这个程序运行下去出现了如上错误，是不是你球面函数和狭缝函数的取样尺寸不同导致此上错误呢？
赞一下 回复此楼
10楼2012-04-25 10:36:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 glddoyss 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] [复试调剂]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 6/300 2026-03-16 19:38 by z1z2z3879
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭