| 查看: 3180 | 回复: 15 | ||
[求助]
matlab程序实现菲涅尔衍射（用fft实现） 哪里不对 已有2人参与
|
|
%光束经菲涅尔衍射后的光强分布 clear,clc; %初始光场 c=linspace(-0.1,0.1,64); r=linspace(-0.1,0.1,64); [x,y]=meshgrid(c,r); [thero,rho]=cart2pol(x,y); l=3;w0=0.03; lambda=1.55e-6;k=2*pi/lambda; u=(sqrt(2)*rho/w0).^l.*exp(-rho.^2/w0^2).*exp(1i*l*thero)*sqrt(2/factorial(l)/pi); %光场表达式 h=abs(u); intensity=h.^2; %显示光强 figure; surf(x,y,intensity); %菲涅尔衍射（角谱衍射理论） z1=10; % 衍射距离 M=64;N=64; n=0.2; %宽度 dx=(M-1)/n;dy=(N-1)/n; %频率间隔 fx=[-31:1:32]'*ones(1,M)*dx/63; %频率坐标 fy=ones(N,1)*[-31:1:32]*dy/63; trans=exp(1i*k*z1*(ones(M,N)-lambda^2*fx.^2-lambda^2*fy.^2).^0.5); %传递函数 uu=fftshift(fft2(fftshift(u))); U=uu.*trans; UU=fftshift(ifft2(fftshift(U))); h=abs(UU); intensity=h.^2; %显示光强 figure; surf(x,y,intensity); 做出来的两张图没有变化。 是频率间隔和频率坐标取得不对吗？如果是，怎么改？ 求大神帮忙！谢谢！ |
回复此楼» 猜你喜欢
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有6人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
2楼2016-06-20 09:19:30
3楼2016-06-20 11:25:24
Berylcczhang
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 979
- 散金: 9
- 帖子: 246
- 在线: 31.3小时
- 虫号: 4682139
- 注册: 2016-05-12
- 专业: 粒子物理学和场论
4楼2016-06-21 16:54:49
5楼2016-06-21 17:25:27
 |
6楼2017-08-17 17:37:38
 |
7楼2017-08-17 23:11:07
|
|
8楼2017-08-18 14:04:51
 |
9楼2017-08-18 22:43:08
 |
10楼2017-08-20 13:32:14