[讨论]TM波垂直入射，我这样引入源和加入平面波正确吗？谢谢！

请问高手们，下面代码中源和按一维fdtd随时间推进方式正确吗？ 代码后面还有对连接边界和角点的处理，设计到下面代码中的Eziny和Hyiny，按照葛老师 K4.GAGd  
书上的方法处理的。但程序运行一百步出的图不像平面波啊。 [UW%(N  
(hb\1wZ  
Source=zeros(1,nmax); qhIO7h  
2A,iY}R  
for n=1:nmax mh5ozv$  
    Source(n)=sin(2*pi*ff*n*dt); -C  
end 6`V2-zv$  
AhNy+p{  
%Timestepping loop 0QakFt  
for n=1:nmax D=!e6E<>@  
     +^@;J?O  
    for I=Ymin+1:Ymax 7?[{/`k~?  
        Eziny(1,I)=Eziny(1,I)+CB*(Hyiny(1,I)-Hyiny(1,I-1)); OS.oknzZZ  
    end I7Uj<a=(q  
     3lW7auH4Y{  
    for I=Ymin:Ymax-1 ?J1x'/G  
        Hyiny(1,I)=Hyiny(1,I)+CQ*(Eziny(1,I+1)-Eziny(1,I)); &OX
nZT3P  
    end Q*GJREC  
(*
XSrQ  
reU*apZ/  
  Eziny(1,3)=Source .. r6u) 6J=  
4ne95_i  

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

Eziny(1,3)=Source(n); aSy^(WN8  
SMMsXH  
你只在(1,3)这点加入射波么？ B_u1FWc  
GW[g!66^  
这加的是点源吧。 MzWVsV  
lebwGW,!  
插入附件中可以加图片的。

是的，在一维Ezinx上设置Ezinx（3，1）＝source（3）设置电源，然后后面和Hyinx一起按照一维fdtd迭代，这样可以引入平面波吗？ J0 [^hH  
A<1:vV  
另外，对入射波要不要单独设置边界吸收条件啊？

这样不行的。 s`Y8&e.Yr  
你考虑一下电源的引入和平面波的引入有什么不一样的，就可以明白的。 -$x5[6bN  
;Nd,K C0k  
需要的，因为平面的特性要求了是朝一个方向传播的，不会在连接边界处发散的。 ]]EOCGZ"  
bh[`uRC}  
而且平面波的加入应该是在一个平面，不会在一个点加入吧。

谢谢你的指点，我修改了一下程序之后，发现Hx分量像点源往外辐射波，但Hy和Ez分量都是平面波，晕了，它们三者不是被联系在一起了吗，怎么表现这么不一样啊？

不太可能吧。。。 D&{7Av  
三者是结合在一起的，彼此影响的啊。

是啊，而且Ez分量的最小值是负一，最大值是0.5，郁闷啊   ! ?o$6w(]''  
zGj0'!!-  
贴不上图，更郁闷！ w
'~f Z*  
pq#Hca[  
上附件！

图片:planewaveTMmovie.gif

图片:plavewaveTMmesh.gif

;:vbOG#aSN  
%L}9nc%~eP  
居然贴上来了，哈哈 Yb]eWLv  
  }N:0%Gk[;  
4&!`Yi_1L  
还请版主、高手们多多指点

Hx Hy Ez的递推公式是什么？ B=%cXW,  
code里面怎么写的？

planewaveTM.rar (3 K) 下载次数:24 c?A(C#~ z  
J@PwN^`  
这是m文件，问题多多，帮我调试调试，嘿嘿！

这个会不会是柱面波？