[求助]FDTD的问题求助

是电流密度和电流强度的问题，如果输入考虑电流强度I，除以网格的面积以电流密度加入Maxwell方程中，这样设置的网格大小对仿真没有影响了

试试看 9e _8Z@|  
你有这方面的迭代的公式么？可以上传上来看看？

引用第10楼cem-uestc于2009-06-23 10:38发表的 : 7dq*e4z)  
是电流密度和电流强度的问题，如果输入考虑电流强度I，除以网格的面积以电流密度加入Maxwell方程中，这样设置的网格大小对仿真没有影响了 _I<LB0kgf.  

`F,*NESv  
刚刚试了试还是有影响，不知道为什么

我这样试了试，好像问题去掉了 [7sy}UH  
        double J_source=I0*(time_iter-n0)/n1*exp(-pow((time_iter-n0)/n1,2))/(dx*dy); T^1]|P  
        hx[snx][sny+1]-= dt/(u0*ur*dx)*J_source; P;DGs]PF  
        hy[snx][sny]-= dt/(u0*ur*dy)*J_source; 7 yp}  
        hx[snx][sny]+= dt/(u0*ur*dx)*J_source; n~C!PXE  
        hy[snx+1][sny]+= dt/(u0*ur*dy)*J_source; 3{FUFx  
这样可以么？ Nt/#Qu2#br  
是不是能够实现我想要的电流注入的问题。

怎么把Hx和Hy同时激励啊 Aat-938FP6  
zilM+BZ8  
dt/(u0*ur*dx)*J_source不用除dx了dt/(u0*ur)*J_source

是不是这个意思？ -FrNk>  
double J_source=I0*(time_iter-n0)/n1*exp(-pow((time_iter-n0)/n1,2))/(dx*dy); 3,[#%}1(S  
KV {J>J1  
hx[snx][sny+1]-= dt/(u0*ur)*J_source; \\{+t<?J  
 HLsG<#  
即可？ NR|t~C+  
/@`kM'1:  
其他的三个激励删除？

仿真的TM模式，一般激励Ez，产生柱面波激励出去

因为用电流密度项激励场，应加在电场的仿真中

你的公式不对吧？在常规FDTD方法中一维时： NpLZ ,|H  
Hy(n+1/2,i+1/2)=Hy(n-1/2,i+1/2)+[Ez(n,i+1)-Ez(n,i)]*dt/dx+Im(n,i+1/2)*dt /P%OXn$i/  
Ez(n+1,i)=Ez(n,i)+[Hy(n+1/2,i+1/2)-Hy(n+1/2,i-1/2)]*dt/dx+Ie(n+1/2,i)*dt 5_7y1  
在算n+1/2时刻的Hy时，n时刻的Ez是已知的，而n时刻的电流Im要么是初边界条件给定的，要么是由Im(n,i+1/2)=sigm(i+1/2)*Hy(n,i+1/2)=sigm(i+1/2)*(Hy(n+1/2,i+1/2)+Hy(n-1/2,i+1/2)算出来的。算n+1时刻的Ez时同理。由于计算Hy和Ez是交替进行的，故每个点（或者半点）只需存储一个Ez（或Hy），根本不需额外的内存消耗。 WRov7  
当采用高阶FDTD时（比如四阶时间精度、四阶空间精度的FDTD格式），需要消耗额外的内存，但由于格式精度更高，网格密度可以降低，故在只需同等精度结果时花费的代价反而有所下降（这也是发展高阶格式的目的）。

那我计算时，用Ez激励，可以通过Ez得到电压，之后在后处理时，如何得到该激励位置处的电流？因为我想得到的是该位置处的输入阻抗特性？