[已解决]关于平面波源的加入问题

引用第9楼wq_463于2009-03-24 13:28发表的 : }!;s.[y  
fortran写的啊.... [表情]
勉强着看吧,谢谢版主了

|(pRaiJ  
%ZNI:Uh  
O(∩_∩)O哈哈~，标准的MATLAB代码，不是Fortran

   ~~oh my god！

看的不是很明白，有些参数不知道怎么出来的

gwzhao写的那个程序是根据葛德彪书第二版的例子，把fortran改编为matlab程序的，所以最好结合葛德彪讲平面波源加入的那些章节看

引用第13楼herosword于2009-03-26 10:59发表的 : C_ \q?>  
gwzhao写的那个程序是根据葛德彪书第二版的例子，把fortran改编为matlab程序的，所以最好结合葛德彪讲平面波源加入的那些章节看

.c#y%S  
)~V4+*<  
谢谢 herosword 师兄

漂亮的很 'SYo_!  
谢谢提供 r|*&GHo L  
受教不浅

,olwwv_8G  
里面的fdtd迭代公式系数是怎么一回事啊？请大侠赐教！ gX/NtO%  
谢谢！

FE2(i,j)=1./((Media(1,i)+Media(1,j))-pi*0.5*(Media(3,i)+Media(3,j))/WL); U/}AiCdj@  
       #-l!`\@  
FH2(i,j)=1./((Media(2,i)+Media(2,j))-pi*0.5*(Media(4,i)+Media(4,j))/WL);  "Nk`RsW  
这里为什么不是： AaJ,=eQ  
FE2(i,j)=1./3.0e8/((Media(1,i)+Media(1,j))-pi*0.5*(Media(3,i)+Media(3,j))/WL); pE9aT5 L  
       .Y=Z!Q  
FH2(i,j)=1./3.0e8/((Media(2,i)+Media(2,j))-pi*0.5*(Media(4,i)+Media(4,j))/WL); 8A`p  

OMIGA=pi/WL; 75QXkJu  
你认为这个程序中的这个有没有问题啊？ X^?|Sz<^E  
首先我认为它的量纲就不对。

应该是对的吧。 a|6x!p2X  
OMIGA 应该是角度 ， pi 也是角度 ， WL是没有量纲的值。 HIWmh4o/.  
WL是将一个波长分成WL份。相当于数值波长了。 波长 = 步长 × WL 。 YKX>@)Dxv  
foP>w4pB