| 查看: 1888 | 回复: 7 | ||
wupeng2320金虫 (小有名气)
|
[求助]
菲涅尔衍射数值模拟 已有2人参与
|
|
最近在做菲涅尔衍射数值模拟，但在模拟上遇到了问题，希望光学加matlab的大神能不啬赐教！ 模拟中，物体是圆孔，然后菲涅尔衍射公式采用的是卷积形式：u(x,y)=exp(jkz)/(j*lamda*z)*【u0(x,y)和exp（jk/(2z)*(x^2+y^2)】。u0（x,y)是物函数，模拟中用上述圆孔做物。二维离散卷积的实现过程是通过conv2进行。但是结果出现之后不是菲涅尔，所以猜想是菲涅尔公式的实现出现了问题，还是说取样的问题？期待各位给解答下，万分感谢。 程序如下（附件是.m文件）： clear all; %做小孔。R=50，圆心（128,128） for x=1:256 for y=1:256 if ((x-128).^2+(y-128).^2<=2500) t(x,y)=1; else t(x,y)=0; end end end %显示圆孔图像 subplot(121) imshow(t,[]) %菲涅尔公式参量设置 z=200; lamda=500*10^(-3); k=2*pi/lamda; %采用的菲涅尔公式形式是：u(x,y)=exp(jkz)/(j*lamda*z)*【u0(x,y)和exp（jk/(2z)*(x^2+y^2)】的卷积 %u0（x,y)是物函数，本程序利用上述圆孔做物函数。 for x=1:256 for y=1:256 f(x,y)=exp(i*k*(x.^2+y.^2)/2/z); end end m=conv2(t,f) u=exp(i*k*z)/(i*lamda*z).*m %展示衍射图 subplot(122) imshow(u.*conj(u),[]) |
回复此楼» 猜你喜欢
一个有机合成实验室都需要哪些设备？
已经有6人回复
-
2026年国自然面上资助率
已经有20人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab模拟菲涅尔直边衍射
已经有9人回复
- 【求助】判断菲涅尔区及夫琅和费衍射区的距离，菲涅尔数
已经有5人回复
michaeldel
至尊木虫 (正式写手)
火星领导
- 应助: 64 (初中生)
- 金币: 11667.2
- 散金: 29
- 红花: 17
- 帖子: 843
- 在线: 509.2小时
- 虫号: 1819742
- 注册: 2012-05-16
- 性别: GG
- 专业: 光学
2楼2015-04-24 17:25:04
michaeldel
至尊木虫 (正式写手)
火星领导
- 应助: 64 (初中生)
- 金币: 11667.2
- 散金: 29
- 红花: 17
- 帖子: 843
- 在线: 509.2小时
- 虫号: 1819742
- 注册: 2012-05-16
- 性别: GG
- 专业: 光学
3楼2015-04-24 17:49:36
wupeng2320
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 607.4
- 散金: 143
- 红花: 6
- 帖子: 252
- 在线: 142.9小时
- 虫号: 2224647
- 注册: 2013-01-06
- 专业: 光学
4楼2015-04-24 20:31:14
wupeng2320
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 607.4
- 散金: 143
- 红花: 6
- 帖子: 252
- 在线: 142.9小时
- 虫号: 2224647
- 注册: 2013-01-06
- 专业: 光学
5楼2015-04-24 20:36:52
6楼2016-11-12 19:46:27
7楼2016-11-14 19:32:24
8楼2019-10-05 15:34:51