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

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

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

求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）已经有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都成夫琅禾费衍射了。