| 查看: 3271 | 回复: 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); 做出来的两张图没有变化。 是频率间隔和频率坐标取得不对吗？如果是，怎么改？ 求大神帮忙！谢谢！ |
回复此楼» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有8人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有11人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有7人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有8人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有9人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有10人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
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