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

24小时热门版块排行榜    

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

匿名

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

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

banuo2005

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-04 00:34:27
已阅   同方向广播   申请物理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的回帖
查看全部 10 个回答

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
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的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-17 8/400 2026-03-21 02:12 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
关闭
请填处理意见
关闭