| 查看: 3424 | 回复: 8 | ||
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) 程序可以运行，但不是我要的结果。 不知哪里出问题了。  tri2.png |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : fresnel.m
- 附件 2 : test1.m
2013-09-22 11:02:24, 523 bytes
2013-09-22 11:02:26, 434 bytes
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有16人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有4人回复
-
系统查不到
已经有10人回复
-
基金不中，共勉
已经有3人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
-
2026国自然函评费到账
已经有25人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有28人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有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人回复
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
感谢参与，应助指数 +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都成夫琅禾费衍射了。  untitled.jpg |
2楼2013-09-22 16:51:58
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
ldmllion3楼2013-09-22 16:54:38
wangwei225
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 857
- 散金: 164
- 红花: 2
- 帖子: 137
- 在线: 22.9小时
- 虫号: 2661906
- 注册: 2013-09-17
- 专业: 集成电路设计与测试
4楼2013-09-22 17:04:36
ldmllion
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 178.3
- 散金: 35
- 帖子: 59
- 在线: 28.5小时
- 虫号: 2019854
- 注册: 2012-09-21
- 专业: 光学信息获取与处理
5楼2013-09-22 18:23:20
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
6楼2013-09-23 11:28:03
ldmllion
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 178.3
- 散金: 35
- 帖子: 59
- 在线: 28.5小时
- 虫号: 2019854
- 注册: 2012-09-21
- 专业: 光学信息获取与处理
7楼2013-09-23 20:41:52
smy1982
木虫 (正式写手)
- 物理EPI: 1
- 应助: 21 (小学生)
- 金币: 3128.7
- 散金: 6748
- 红花: 12
- 帖子: 853
- 在线: 450.4小时
- 虫号: 305139
- 注册: 2006-12-09
- 性别: GG
- 专业: 光学信息获取与处理
8楼2013-09-24 08:51:00
9楼2016-03-28 10:36:42
