| 查看: 3313 | 回复: 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); 做出来的两张图没有变化。 是频率间隔和频率坐标取得不对吗？如果是，怎么改？ 求大神帮忙！谢谢！ |
回复此楼
Berylcczhang
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 979
- 散金: 9
- 帖子: 246
- 在线: 31.3小时
- 虫号: 4682139
- 注册: 2016-05-12
- 专业: 粒子物理学和场论
4楼2016-06-21 16:54:49
2楼2016-06-20 09:19:30
3楼2016-06-20 11:25:24
5楼2016-06-21 17:25:27
