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

请问高手们，下面代码中源和按一维fdtd随时间推进方式正确吗？ 代码后面还有对连接边界和角点的处理，设计到下面代码中的Eziny和Hyiny，按照葛老师 ZnXq+^Z4  
书上的方法处理的。但程序运行一百步出的图不像平面波啊。 m L,
El2  
&-B&s.,kj  
Source=zeros(1,nmax); mj^]e/s%  
M hJ;)(  
for n=1:nmax z<Y >phc  
    Source(n)=sin(2*pi*ff*n*dt); P6([[mmG  
end +ug[TV  
qcdENIy0b  
%Timestepping loop dq U.2~9  
for n=1:nmax L|
pJ\~  
     EC0M0qQ  
    for I=Ymin+1:Ymax x}*Y =Xh  
        Eziny(1,I)=Eziny(1,I)+CB*(Hyiny(1,I)-Hyiny(1,I-1)); r\f|r$i  
    end m<wEw-1.  
     oEuo@\U05v  
    for I=Ymin:Ymax-1 F)Qj<6  
        Hyiny(1,I)=Hyiny(1,I)+CQ*(Eziny(1,I+1)-Eziny(1,I)); F?!  
    end 4}>1I}!k  
mzoNXf:x  

|=U(8t  
  Eziny(1,3)=Source .. QnPgp(d<  
@[] A&)B  

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

Eziny(1,3)=Source(n); $o"PQ!z  
I54O9Aoy  
你只在(1,3)这点加入射波么？ 9GdrJ~h  
O/IW.t  
这加的是点源吧。 V;Zp3Qo!  
({cgak  
插入附件中可以加图片的。

是的，在一维Ezinx上设置Ezinx（3，1）＝source（3）设置电源，然后后面和Hyinx一起按照一维fdtd迭代，这样可以引入平面波吗？ f9D7T|J?10  
gsc*![N  
另外，对入射波要不要单独设置边界吸收条件啊？

这样不行的。 o{UwUMw5`  
你考虑一下电源的引入和平面波的引入有什么不一样的，就可以明白的。 4sZ^:h,1  
&g) `  
需要的，因为平面的特性要求了是朝一个方向传播的，不会在连接边界处发散的。  &.(ZO]  
%Si3t2W/  
而且平面波的加入应该是在一个平面，不会在一个点加入吧。

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

不太可能吧。。。 |"ck;.)  
三者是结合在一起的，彼此影响的啊。

是啊，而且Ez分量的最小值是负一，最大值是0.5，郁闷啊   ! Nr6YQH*[  
zxTm`Dh;[
 
贴不上图，更郁闷！ Su-LZ'C\  
;m@>v?zE  
上附件！

图片:planewaveTMmovie.gif

图片:plavewaveTMmesh.gif

N@Bqe{r6j  
Nbuaw[[iz  
居然贴上来了，哈哈 5/>G)&  
  $1#|<|  
+!'6:F  
还请版主、高手们多多指点

Hx Hy Ez的递推公式是什么？ OI::0KOv  
code里面怎么写的？

planewaveTM.rar (3 K) 下载次数:24 oDrfzm|[Y  
;
mjk`6p  
这是m文件，问题多多，帮我调试调试，嘿嘿！

这个会不会是柱面波？