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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 凝聚态物理 » matlab模拟菲涅尔直边衍射
101/1返回列表
查看: 2969  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
关闭
请填处理意见
关闭