应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 基础物理 » matlab程序实现菲涅尔衍射（用fft实现） 哪里不对
161/2返回列表12下一页
查看: 3303  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

你倒是跑啊

金虫 (初入文坛)

[求助] 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);


做出来的两张图没有变化。   是频率间隔和频率坐标取得不对吗？如果是，怎么改？
求大神帮忙！谢谢！
回复此楼

» 猜你喜欢
1楼 2016-06-19 16:29:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

你倒是跑啊

金虫 (初入文坛)
不要沉，求助攻！！！！谢谢啊
赞一下 回复此楼
2楼2016-06-20 09:19:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hcquan

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
你倒是跑啊: 金币+15, ★★★很有帮助, 谢谢！ 2016-06-21 17:26:40

你把距离设得太小了，z1=10，你这里的单位是um，把它设大点，程序没问题。
赞一下 回复此楼
3楼2016-06-20 11:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Berylcczhang

铜虫 (小有名气)
是傅立叶变换么？

发自小木虫IOS客户端
赞一下 回复此楼
4楼2016-06-21 16:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

你倒是跑啊

金虫 (初入文坛)
引用回帖:
4楼: Originally posted by Berylcczhang at 2016-06-21 16:54:49
是傅立叶变换么？

是的
回复此楼
5楼2016-06-21 17:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

unzchnz868
6楼2017-08-17 17:37:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

unzchnz868
7楼2017-08-17 23:11:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

unzchnz868
8楼2017-08-18 14:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

unzchnz868
9楼2017-08-18 22:43:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

unzchnz868
10楼2017-08-20 13:32:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 你倒是跑啊 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 考研调剂 +3 呼呼？~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
信息提示
关闭
请填处理意见
关闭