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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 基础物理 » 求大神赐教一下matlab光学仿真的一程序。
91/1返回列表
查看: 3447  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ldmllion

新虫 (小有名气)

[求助] 求大神赐教一下matlab光学仿真的一程序。

我想做一个一束光通过任意图案之后的衍射图案的matlab光学仿真。
我的思路是将光波表达式与图案相乘，然后做菲涅尔衍射。不知是不是这样，求大神指导。
主代码如下：
clear all;clc;
cm=0.01;
um=1e-6;
a=rgb2gray(imread('tri2.png'));
a=im2double(a);%图案
[M,N]=size(a);
z=1;
lambda=0.63*um;
k=2*pi/lambda
hx=0.5*cm;hy=0.5*cm;
dhx=hx/M;dhy=hy/N;
x=dhx.*(ones(N,1)*[-M/2:M/2-1]);
y=dhy.*(ones(N,1)*[-M/2:M/2-1])';
b=exp(1i.*5*atan2(y,x));%用此光照射图案a
f0=a.*b;
[f2,dx1,dy1,x1,y1]=fresnel(f0,N,N,dhx,dhy,z,lambda);%fresnel函数见附件
f2=abs(f2).^2;
f2=mat2gray(f2);
imshow(f2)

程序可以运行，但不是我要的结果。
不知哪里出问题了。求大神赐教一下matlab光学仿真的一程序。
tri2.png
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : fresnel.m
    • 2013-09-22 11:02:24, 523 bytes
  • 附件 2 : test1.m
    • 2013-09-22 11:02:26, 434 bytes

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-22 11:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smy1982

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
华丽的飘过: 金币+10, 3q 2013-10-01 21:07:31
你要什么图片啊？
clear all;clc;
cm=0.01;
um=1e-6;
a=rgb2gray(imread('tri2.png'));
a=im2double(a);
[M,N]=size(a);
z=100e-3;
lambda=0.63*um;
k=2*pi/lambda
hx=0.5*cm;hy=0.5*cm;
dhx=hx/M;dhy=hy/N;
x=dhx.*(ones(N,1)*(-M/2:M/2-1));
y=dhy.*(ones(N,1)*(-M/2:M/2-1))';
dx0=lambda*z/(M*dhx);
dy0=lambda*z/(M*dhy);

b=exp(1i.*5*atan2(y,x));
% imshow((angle(b)+pi)/(2*pi))
f0=a.*b;
[f2,dx1,dy1,x1,y1]=fresnel(f0,N,N,dhx,dhy,z,lambda);
f2=abs(f2).^2/max(max(abs(f2).^2));
f2=mat2gray(f2);
imshow(f2)

你那源程序里的z值太大了吧，1m都成夫琅禾费衍射了。
求大神赐教一下matlab光学仿真的一程序。-1
untitled.jpg
赞一下(2人) 回复此楼
如果你曾歌颂黎明，那么也请你拥抱黑夜。
2楼2013-09-22 16:51:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smy1982

木虫 (正式写手)

【答案】应助回帖

还有你那“fresnel.m ”程序用的角谱法算的，不适用于远场衍射，只能模拟近场的，太远了，就有错误了。
赞一下(1人) 回复此楼

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

ldmllion
如果你曾歌颂黎明，那么也请你拥抱黑夜。
3楼2013-09-22 16:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwei225

金虫 (小有名气)
引用回帖:
3楼: Originally posted by smy1982 at 2013-09-22 16:54:38
还有你那“fresnel.m ”程序用的角谱法算的，不适用于远场衍射，只能模拟近场的，太远了，就有错误了。

有道理
回复此楼
4楼2013-09-22 17:04:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldmllion

新虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by smy1982 at 2013-09-22 16:54:38
还有你那“fresnel.m ”程序用的角谱法算的，不适用于远场衍射，只能模拟近场的，太远了，就有错误了。

那请问我这思路有没有错，用目标光波（非平面波或球面波）与孔径光阑（任意的图案）相乘来做菲涅尔衍射？
赞一下 回复此楼
5楼2013-09-22 18:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smy1982

木虫 (正式写手)
引用回帖:
5楼: Originally posted by ldmllion at 2013-09-22 18:23:20
那请问我这思路有没有错，用目标光波（非平面波或球面波）与孔径光阑（任意的图案）相乘来做菲涅尔衍射？...

思路正确
回复此楼
如果你曾歌颂黎明，那么也请你拥抱黑夜。
6楼2013-09-23 11:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldmllion

新虫 (小有名气)
引用回帖:
6楼: Originally posted by smy1982 at 2013-09-23 11:28:03
思路正确...

我做的就是近场衍射的模拟，程序中z值取的是有点大，但我改小后，还是不是我要的结果图。文献中做的模拟图是类似于这样的（见附件）。
请问这会是哪里出问题了？
是用到菲涅尔函数出问题了吗，如果是的话，该怎么改？
谢谢了。
求大神赐教一下matlab光学仿真的一程序。-2
图片1.png


求大神赐教一下matlab光学仿真的一程序。-3
图片2.png
赞一下 回复此楼
7楼2013-09-23 20:41:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smy1982

木虫 (正式写手)
传文献看看
回复此楼
如果你曾歌颂黎明，那么也请你拥抱黑夜。
8楼2013-09-24 08:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

玉之鸣人

新虫 (初入文坛)
图片怎么来的？比方说如图的三角形
赞一下 回复此楼
9楼2016-03-28 10:36:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ldmllion 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 03:50 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 03:40 by k06BKNblrGRH
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-17 01:26 by Tql5LQhh5rLK
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 i3KBRD19ERkG 2026-09-16 3/150 2026-09-17 00:50 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-16 22:40 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 22:16 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 21:14 by MHf40yA0Uhyb
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 21:03 by MHf40yA0Uhyb
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 18:37 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
信息提示
关闭
请填处理意见
关闭