24小时热门版块排行榜    

查看: 1801  |  回复: 7
【悬赏金币】回答本帖问题,作者oewindbell将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

oewindbell

铜虫 (小有名气)

[求助] matlab 利用快速傅里叶变换计算谐振腔的模式

%%利用二维快速傅里叶变换方法计算谐振腔中的模式
clear;
clc;
lambda=1.064*10^-6;      %波长
k=2*pi/lambda;           %波数
a=0.5*10^-3;             %镜子的半径
fnum=5;                  %菲涅尔数
z=a^2/fnum/lambda;       %腔长
repet=100;               %迭代次数
M=129;N=129;
x=linspace(-a,a,M);
y=linspace(-a,a,N);
[X,Y]=meshgrid(x,y);
u0=ones(M,N);
u1=zeros(M,N);
H=exp(1j*k*z).*exp(-1j*pi*lambda*z.*(X.^2+Y.^2));
for inter=1:repet
    U0=FFT(u0);
    U1=U0.*H;
    u1=ifft(U1);
    u1=u1./max(abs(u1());
    u0=u1;
end
Ie=u1.*conj(u1);
Ie=Ie./max(Ie();
surf(X,Y,Ie);
shading interp;
axis([-a,a,-a,a,0,1]);
box on; grid off;
本人是matlab初学者,请各位大神帮忙指点一二!
谢谢了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lty0709

铜虫 (小有名气)

我能知道    是什么么
心中有梦想,眼中才有光泽!
4楼2013-08-21 11:14:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oewindbell 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +7 潇洒怡惜 2024-06-12 12/600 2024-06-18 14:03 by nano求求了
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+4 BruceChum 2024-06-15 21/1050 2024-06-18 13:03 by 冰蓝夜游神
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +8 香瓜木香 2024-06-12 14/700 2024-06-18 11:16 by 晓目崇
[基金申请] 有人中过人文社科类的博后特助吗? +3 outsider1986 2024-06-16 5/250 2024-06-18 11:10 by 袁天未然
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[硕博家园] 硕博巨婴,也许才刚刚开始 +33 SNaiL1995 2024-06-12 85/4250 2024-06-18 00:06 by 博士求上上岸
[基金申请] 化学会年会改成一年开一次了? +7 babu2015 2024-06-12 15/750 2024-06-17 20:28 by 落笔画枫
[有机交流] 苯酚和炔醇反应,mitsunobu反应 5+3 TONGMEIMEI 2024-06-12 9/450 2024-06-17 15:00 by dschong
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +10 zju2000 2024-06-12 16/800 2024-06-17 14:54 by 大天尊
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[论文投稿] 求机械类四区sci推荐 5+3 迷茫小旷 2024-06-14 4/200 2024-06-15 11:25 by bobvan
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
请填处理意见