24小时热门版块排行榜

返回列表

ldmllion

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 178.3
散金: 35
帖子: 59
在线: 28.5小时
虫号: 2019854
注册: 2012-09-21
专业: 光学信息获取与处理

[求助] 求大神赐教一下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)

程序可以运行，但不是我要的结果。
不知哪里出问题了。

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

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有11人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有16人回复
初秋的晨风已经有3人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
面上没中，邀请各位路过的虫友分析一下分数已经有11人回复
Ei源刊怎么投已经有4人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）已经有6人回复
求大神帮忙设计一个工装夹具已经有6人回复
求各位大神一个问题已经有4人回复
编写matlab程序：y=aln（x）+b，求a、b及拟合的R？已经有15人回复
fluent UDF代码调用预定义程序的问题（如fortran程序）已经有15人回复
请教LSDA+U计算的问题已经有24人回复
小妹菜鸟一枚，请教大神帮忙看我的数据可以拟合成什么样的方程，多谢已经有6人回复
发光材料的基质对于激活剂是否存在能量传递啊啊啊哪位大神教教我已经有4人回复
perl程序求助，请大神们赐教~~~ 已经有9人回复
计算模拟编程用，求大神推荐笔记本已经有16人回复
小弟（高三学生）大学想学与机器人有关的专业,希望各路大神指教！.. 已经有18人回复
matlab计算面积求助！！难已经有8人回复
光学涡旋计算全息图的matlab程序已经有9人回复
ANSYS二维电场模拟电阻变化，求指教已经有4人回复
如何同时模拟发热体的辐射和传导？已经有10人回复
【求助】fluent二维模型中热流密度的设置已经有12人回复
【求助】abaqus中tie和coupling的区别？？已经有19人回复

1楼 2013-09-22 11:03:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

smy1982

木虫 (正式写手)

物理EPI: 1
应助: 21 (小学生)
金币: 3128.7
散金: 6748
红花: 12
帖子: 853
在线: 450.4小时
虫号: 305139
注册: 2006-12-09
性别: GG
专业: 光学信息获取与处理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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都成夫琅禾费衍射了。