[资料共享]一维等离子体FDTD的Matlab源代码(两种方法）

如果你不考虑碰撞的话，那电磁波就没有衰减了啊。 +'V ,z  
从色散关系里面也可以看出来，碰撞对应着能量的衰减。 5H!6m_,w  
,Pj UlcO_  
不太了解你说的碰撞和等离子体频率对结果的影响差不多是从什么方面说的？对反射，入射效率还是波形？ v|IG G'r  

引用第8楼itnice于2008-10-20 11:47发表的  :
/NB;eV?  
可以先S1与Ex的关系吗？ #`qP7E w  
主要是 AGMrBd|J{  
     Sx1(i)=(1+exp(-1*u0*dt))*Sx2(i)-exp(-1*u0*dt)*Sx3(i)+(wpe^2*dt/(2*u0))*(1-exp(-1*u0*dt))*Ex(i); v5A8"&Jr  
         Ex(i)=Dx(i)/epsz-Sx1(i); } m&La4E  
         Ex(i)=Ex(i); >A.m`w  
.......

你分析一下稳定性就可以看出来为什么dt = dx/c是不行的了。 .E&-gXJ4  
`&!J6)OJ  
简单来说，如果在模拟过程中出现速度大于光速的，不稳定等现象就会出现了。 ar[*!:!  
真空中取dt = dx /c 问题不大，但等离子体模拟中dx/dt =c，等离子体中稍微有点波动，不就一下子超过光速了么？ (|_N2R!  
具体你可以推导一下看看，其实相对真空中的稳定性还是不一样的。 Kfk/pYMDq  
比如说，认为波速度在等离子中小于光速，就会得到dt = dx / c稳定的结论，这显然是不对的。 u
!DAeE  
电子速度也不能超过光速的。 %VmHw~xyF:  

引用第9楼itnice于2008-10-20 11:54发表的  : `2.2; Vk  
还有就是步长的设置为什么是dt=dx/2c而不是dt=dx/c呢。但dt=dx/c在这里好像是不稳定的

第二种方法是用的Z变换。 C ks;f
6G  
推导过程如下：

图片好像不是太清楚，可以把电子书的名字说一下吗？

好像附件不好用啊。图片也不清晰。

不是电子书，是我自己写的推导过程啊。

再试一次看看

哈哈，高人啊，都可以自己写书了 ~#E&E%s
J  
不过我还是先把第一个方法搞明白。可以说一下第一个方法用的是什么公式吗？

写错了，刚才的资料是移位算子得到的递推公式。 Vx!ZF+  
第一种方法用的是Z变化。

原来是通过H来求出D，之后再用介电常数来求出E，之后再用E来求下一步的H。不过还是不明白Sx代表的是什么