| 查看: 1672 | 回复: 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),[]) |
回复此楼
8楼2019-10-05 15:34:51
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
20