[发明专利]一种应用于时域有限差分法的透明激励源的实现方法

权利要求书

1.一种应用于时域有限差分法的透明激励源的实现方法，具体包括以下步骤：

步骤1、采用Yee网格，设波导结构沿z向分布，波导端口位于x-y平面，所有分布在激励源网格处的电磁场分量具有三套值，分别以下标P、N和R来区别，即E_Px、E_Py、H_Pz，E_Nx、E_Ny、H_Nz和E_Rx、E_Ry、H_Rz；P表示与正源相关，N表示与负源相关，R表示与散射场相关，下标PT表示Px或Py，对于NT、RT同理；

步骤2、

根据激励源方程f_T(i,j)来计算表示纵向磁场的辅助项f_H(i,j)，由于f_T(i,j)为已知项，所有时刻的f_H(i,j)项都可以求得

计算所有时刻的反馈电场R_NT(i,j)、R_PT(i,j)

步骤3、

更新激励源网格点处的电场

其中，表示在求解时，当涉及激励源网格上的磁场时使用正源对应的磁场值，对于同理，具体公式为

散射电场可得

然后应用标准离散麦克斯韦方程更新其余网格点的电场，同时更新完美匹配层吸收边界电场；

步骤4、

更新激励源网格点处的磁场

其中，表示在求解时，当涉及激励源网格上的电场时使用正源对应的电场值，对于同理，具体公式为

接着更新与激励源网格点相邻的横向磁场

其中，

散射磁场可得

然后按照标准离散麦克斯韦方程更新其余网格点的磁场，同时更新完美匹配层吸收边界磁场；

步骤5、重复3、4步骤，从n＝0迭代直至达到n＝nt，nt为所设置的模拟步数。

2.如权利要求1所述应用于时域有限差分法的透明激励源的实现方法，其特征在于：适用于二维时，分布在激励源网格处的三套电磁场分量分别对应为E_Py、H_Pz，E_Ny、H_Nz和E_Ry、H_Rz；电磁场分量对应公式相应退化。

3.如权利要求1所述应用于时域有限差分法的透明激励源的实现方法，其特征在于：适用于一维时，分布在激励源网格处的三套电磁场分量对应为E_Py、E_Ny和E_Ry；电磁场分量对应公式相应退化。