| 查看: 2074 | 回复: 7 | ||
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初学者，请各位大神帮忙指点一二！ 谢谢了 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有121人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于matlab进行傅里叶变换和逆变换的程序，求大神指导！
已经有9人回复
- 求教matlab傅里叶变换求解
已经有5人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- matlab 傅里叶变换 频谱分析
已经有12人回复
- 如何用matlab进行快速傅里叶变换？（请教）
已经有5人回复
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
2楼2013-08-20 22:43:42
oewindbell
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 324.9
- 红花: 1
- 帖子: 120
- 在线: 213.9小时
- 虫号: 1828175
- 注册: 2012-05-21
- 性别: GG
- 专业: 光学
3楼2013-08-20 23:44:45
是什么么 
4楼2013-08-21 11:14:46
oewindbell
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 324.9
- 红花: 1
- 帖子: 120
- 在线: 213.9小时
- 虫号: 1828175
- 注册: 2012-05-21
- 性别: GG
- 专业: 光学
5楼2013-08-21 12:30:06
6楼2013-08-22 17:48:58
oewindbell
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 324.9
- 红花: 1
- 帖子: 120
- 在线: 213.9小时
- 虫号: 1828175
- 注册: 2012-05-21
- 性别: GG
- 专业: 光学
|
%%利用二维快速傅里叶变换方法计算谐振腔中的模式 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=fft2(u0); U1=U0.*H; u1=ifft2(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; 这是修改后的，结果还是有点问题 |
7楼2013-08-23 08:20:56
oewindbell
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 324.9
- 红花: 1
- 帖子: 120
- 在线: 213.9小时
- 虫号: 1828175
- 注册: 2012-05-21
- 性别: GG
- 专业: 光学
8楼2013-08-23 08:22:32