4488阅读
3回复

[求助]matlab如何处理平面波三维图

UID ：20925

注册：2008-11-06
登录：2021-04-22
发帖：227
等级：仿真三级

0楼发表于: 2009-05-15 14:54:37

根据taflove书中三维入射波的设置,

[attachment=19939]
我令θ=0，Φ=0，Ψ=0，这样平面波就可以延y轴传播,我现在想观察平面波是否满足要求,以前看到论坛上有类似的帖子
xx=1:50;
yy=1:50;
  [xxx,yyy]=meshgrid(xx,yy);
  mesh(xxx,yyy,Ez(1:50,1:50,15,1));
  surf(xxx,yyy,Ez(1:50,1:50,15,1));
  pause(0.1);
不知道gwzhao师兄和其他的高手们能解释一下这段程序
Ez(1:50,1:50,15,1)中1该怎么理解< ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

1楼发表于: 2009-05-16 14:17:50

本身这段程序没太多好说的，就是观察一个面上的变化。
比如说要观察z = z0 通过这一点的XY平面的图形。
xx=1:50;  x 的范围是从1 到 50
yy=1:50;  y的范围也是1到50
  [xxx,yyy]=meshgrid(xx,yy); 换分成一个50 ×５０　的平面
  mesh(xxx,yyy,Ez(1:50,1:50,15,1));
  surf(xxx,yyy,Ez(1:50,1:50,15,1));
最后的图形就是这个面上的值了。

不过这样看不够具体，你可以用颜色来表示值的大小，看起来比较直观点。

任何一本matlab书上都有讲的。

逆流而上

离线wq_463

UID ：20925

注册：2008-11-06
登录：2021-04-22
发帖：227
等级：仿真三级

2楼发表于: 2009-05-17 16:29:46

回 1楼(gwzhao) 的帖子

我想做个动画的图，能看出平面波的传播过程
我参考taflove的程序写了下面的一段程序，不过提示有错
我任意固定一个x＝1
timestep=int2str(n);
    tview(:,:)=squeeze(ez(1,jh_bc:upml+je,kh_bc:upml+ke));%固定x观察一个yz面上ez的变化，平面波是延y方向传播的


    subplot('position',[0.15 0.57 0.7 0.35])
    imagesc(tview');
    caxis([-1.1786e+003  1.1786e+003]);
    colorbar;
    axis image; axis xy;
    title(['E_z(i,j,k=k_s_o_u_r_c_e), time step = ',timestep]);
    xlabel('i coordinate');
    ylabel('j coordinate');


    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之类的，新手，见谅

发帖回复