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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 基础物理 » matlab程序实现菲涅尔衍射（用fft实现） 哪里不对
161/2返回列表12下一页
查看: 3478  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了，能看出什么问题？ +3 基诺咪客 2026-08-17 3/150 2026-08-17 15:25 by xinshengwu
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 9/450 2026-08-17 14:59 by wangze12014
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +33 医学老男孩 2026-08-13 80/4000 2026-08-17 14:40 by yq199112088
[基金申请] 感觉是下周放榜了 +3 angus9576 2026-08-17 6/300 2026-08-17 14:32 by 启萌科技
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6GojgJvkDudM 2026-08-16 6/300 2026-08-17 13:58 by EYK67XJMj64E
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭