[发明专利]点源激励下导体平板与非平行介质面强散射点预估方法

权利要求书

1.一种点源激励下导体平板与非平行介质面强散射点预估方法，其特征在于，包含以下步骤：

S1、获取无限大的导体平板的格林函数；

S2、采用离散复镜像法获取半空间介质面的格林函数；

S3、对尺寸有限的导体平板进行边缘绕射修正，获取有限尺寸导体平板的格林函数；

S4、基于S2和S3，进行强散射点的位置预估；

所述的S2中，离散复镜像法是将谱域格林函数分为准静态、表面波、复镜像三个部分进行计算，具体包含以下步骤：

S21、准静态项的求取：

谱域格林函数为：

其中，k_zs表示源点所在层的z方向的波数；F(k_ρ)表示格林函数分量为k_ρ的傅里叶分量；

准静态的贡献为：

其中，当k_ρ→∞时，F₀(k_ρ)的值趋向于常数K₀，其即为准静态的贡献；

使用索末菲尔德积分恒等式，得到空间域的准静态贡献为：

其中，ρ为场点柱坐标的ρ分量；k_s为源点的波数；

S22、表面波项的求取：

提取准静态贡献后的谱域格林函数为：

提取上式中表面波极点的留数，并根据各个极点的留数得到表面波的贡献为：

其中，k_pi表示k_ρ平面上的极点位置；Resi表示极点的留数；

使用索末菲尔德积分恒等式，得到空间域的表面波贡献为：

其中，为零阶第二类汉克尔函数；

S23、复镜像项的求取：

提取准静态贡献和表面波贡献后的格林函数为：

其中，F₃(k_ρ)为提取准静态项和表面波后的格林函数分量为k_ρ的傅里叶分量；

使用索末菲尔德积分恒等式，得到空间域的复镜像贡献为：

其中，a_i为留数，b_i为极点。

2.如权利要求1所述的点源激励下导体平板与非平行介质面强散射点预估方法，其特征在于，所述的S1中，采用单一镜像法获取无限大的导体平板的格林函数。

3.如权利要求2所述的点源激励下导体平板与非平行介质面强散射点预估方法，其特征在于，所述的S22中，提取表面波极点的留数的方法，具体包含以下步骤：

S221、表面波极点定位：在k_ρ平面内框定一个矩形区域，对该矩形四条边线作数值积分；若环线积分为零，则表明该矩形区域内无极点；反之则表明矩形区域内存在极点，并再将该矩形区域分成四个子矩形，分别对四个子矩形递归的进行环线积分，直至找到所有极点；

S222、在每个极点附近取一条环线，对每条环线作数值积分，得到每个极点的留数。