登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
matlab如何处理平面波三维图
发帖
回复
4486
阅读
3
回复
[
求助
]
matlab如何处理平面波三维图
离线
wq_463
UID :20925
注册:
2008-11-06
登录:
2021-04-22
发帖:
227
等级:
仿真三级
0楼
发表于: 2009-05-15 14:54:37
根据taflove书中三维入射波的设置,
]ikomCg
8z2Rry w
[attachment=19939]
/KCPpERk{
我令θ=0,Φ=0,Ψ=0,这样平面波就可以延y轴传播,我现在想观察平面波是否满足要求,以前看到论坛上有类似的帖子
]]0,|My7
xx=1:50;
FS"eM"z
yy=1:50;
X[](Kj^`<
[xxx,yyy]=meshgrid(xx,yy);
qQ0C ?
mesh(xxx,yyy,Ez(1:50,1:50,15,1));
t*G/]
surf(xxx,yyy,Ez(1:50,1:50,15,1));
x%N\5 V1
pause(0.1);
P[ r];e
不知道gwzhao师兄和其他的高手们能解释一下这段程序
80pid[F
Ez(1:50,1:50,15,1)中1该怎么理解< ..
?F7o!B
y "w|g~x]c
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
离线
gwzhao
方恨少
UID :17098
注册:
2008-08-24
登录:
2019-01-09
发帖:
1374
等级:
荣誉管理员
1楼
发表于: 2009-05-16 14:17:50
本身这段程序没太多好说的,就是观察一个面上的变化。
/ $ :j
比如说要观察z = z0 通过这一点的XY平面的图形。
l>gI&1)%
xx=1:50; x 的范围是从1 到 50
j(:I7%3&(*
yy=1:50; y的范围也是1到50
7uorQfR?
[xxx,yyy]=meshgrid(xx,yy); 换分成一个50 ×50 的平面
f3j{V N
mesh(xxx,yyy,Ez(1:50,1:50,15,1));
>(w2GD?
surf(xxx,yyy,Ez(1:50,1:50,15,1));
\0l>q ,
最后的图形就是这个面上的值了。
0P^L }VVX
oI:o"T77sA
不过这样看不够具体,你可以用颜色来表示值的大小 ,看起来比较直观点。
D\w h;r
&6 s) X
任何一本matlab书上都有讲的。
共
条评分
逆流而上
离线
wq_463
UID :20925
注册:
2008-11-06
登录:
2021-04-22
发帖:
227
等级:
仿真三级
2楼
发表于: 2009-05-17 16:29:46
回 1楼(gwzhao) 的帖子
我想做个动画的图,能看出平面波的传播过程
WW+l' 6.
我参考taflove的程序写了下面的一段程序,不过提示有错
k#8Ti"0
我任意固定一个x=1
ES~^M840f
timestep=int2str(n);
6S6f\gAM
tview(:,:)=squeeze(ez(1,jh_bc:upml+je,kh_bc:upml+ke));%固定x观察一个yz面上ez的变化,平面波是延y方向传播的
H>9CW<8
alh >"9~!
p4IZ
subplot('position',[0.15 0.57 0.7 0.35])
Ol"3a|
imagesc(tview');
)"f N!9,F
caxis([-1.1786e+003 1.1786e+003]);
Uv YF[@
colorbar;
v_e9}yI
axis image; axis xy;
W$U0[^1
title(['E_z(i,j,k=k_s_o_u_r_c_e), time step = ',timestep]);
&l$Q^g
xlabel('i coordinate');
Om{l>24i.\
ylabel('j coordinate');
O }9KJU
b'MSkEiQG
*k"|i*{
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之类的,新手,见谅
共
条评分
发帖
回复