登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
matlab如何处理平面波三维图
发帖
回复
4487
阅读
3
回复
[
求助
]
matlab如何处理平面波三维图
离线
wq_463
UID :20925
注册:
2008-11-06
登录:
2021-04-22
发帖:
227
等级:
仿真三级
0楼
发表于: 2009-05-15 14:54:37
根据taflove书中三维入射波的设置,
Q8DQlqHm
|fSe>uVZ
[attachment=19939]
nWMmna.5
我令θ=0,Φ=0,Ψ=0,这样平面波就可以延y轴传播,我现在想观察平面波是否满足要求,以前看到论坛上有类似的帖子
!!Aj<*%
xx=1:50;
3:8{"md@2
yy=1:50;
6IWxPt~
[xxx,yyy]=meshgrid(xx,yy);
Hd,p!_
mesh(xxx,yyy,Ez(1:50,1:50,15,1));
aj?a^}X
surf(xxx,yyy,Ez(1:50,1:50,15,1));
" t7M3i_
pause(0.1);
zxf"87se
不知道gwzhao师兄和其他的高手们能解释一下这段程序
>2]JXLq
Ez(1:50,1:50,15,1)中1该怎么理解< ..
W2z*91$
b^p"|L
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
离线
gwzhao
方恨少
UID :17098
注册:
2008-08-24
登录:
2019-01-09
发帖:
1374
等级:
荣誉管理员
1楼
发表于: 2009-05-16 14:17:50
本身这段程序没太多好说的,就是观察一个面上的变化。
|,aG%MTL
比如说要观察z = z0 通过这一点的XY平面的图形。
Z(9u<
xx=1:50; x 的范围是从1 到 50
8HZs>l
yy=1:50; y的范围也是1到50
P<Zh XN'
[xxx,yyy]=meshgrid(xx,yy); 换分成一个50 ×50 的平面
sF7^qrVQP9
mesh(xxx,yyy,Ez(1:50,1:50,15,1));
%eJGte-
surf(xxx,yyy,Ez(1:50,1:50,15,1));
qVdwfT{1J
最后的图形就是这个面上的值了。
]I L;`>Gp
_ z;q9&J)
不过这样看不够具体,你可以用颜色来表示值的大小 ,看起来比较直观点。
~`D|IWMDq
\T;\XAGr
任何一本matlab书上都有讲的。
共
条评分
逆流而上
离线
wq_463
UID :20925
注册:
2008-11-06
登录:
2021-04-22
发帖:
227
等级:
仿真三级
2楼
发表于: 2009-05-17 16:29:46
回 1楼(gwzhao) 的帖子
我想做个动画的图,能看出平面波的传播过程
6&!&\
我参考taflove的程序写了下面的一段程序,不过提示有错
n 3lE,b
我任意固定一个x=1
kEO7PK/
timestep=int2str(n);
*c)uGz'cD
tview(:,:)=squeeze(ez(1,jh_bc:upml+je,kh_bc:upml+ke));%固定x观察一个yz面上ez的变化,平面波是延y方向传播的
fS:1^A2,
rS0DSGDq
KE1ao9H8wR
subplot('position',[0.15 0.57 0.7 0.35])
suH&jE$ x
imagesc(tview');
!J=sk4T
caxis([-1.1786e+003 1.1786e+003]);
k"N(o(
colorbar;
Pub0IIs
axis image; axis xy;
vm'5s]kdh
title(['E_z(i,j,k=k_s_o_u_r_c_e), time step = ',timestep]);
Q.Aw2
xlabel('i coordinate');
J|3E- p\o
ylabel('j coordinate');
yi^b)2G
xP=/N!,#
unJ R=~E
pause(0.05)
共
条评分
离线
lhp00158
Get busy living or get busy dy
UID :68540
注册:
2010-10-30
登录:
2010-11-24
发帖:
7
等级:
旁观者
3楼
发表于: 2010-10-30 20:22:57
为什么平面波的加入人们的程序都用—inc,-low之类的,新手,见谅
共
条评分
发帖
回复