| 查看: 1539 | 回复: 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),[]) |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
matlab模拟菲涅尔直边衍射
已经有9人回复- 【求助】判断菲涅尔区及夫琅和费衍射区的距离，菲涅尔数
已经有5人回复
michaeldel
至尊木虫 (正式写手)
火星领导
- 应助: 64 (初中生)
- 金币: 11667.2
- 散金: 29
- 红花: 17
- 帖子: 843
- 在线: 509.2小时
- 虫号: 1819742
- 注册: 2012-05-16
- 性别: GG
- 专业: 光学
3楼2015-04-24 17:49:36
michaeldel
至尊木虫 (正式写手)
火星领导
- 应助: 64 (初中生)
- 金币: 11667.2
- 散金: 29
- 红花: 17
- 帖子: 843
- 在线: 509.2小时
- 虫号: 1819742
- 注册: 2012-05-16
- 性别: GG
- 专业: 光学
2楼2015-04-24 17:25:04
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