复杂介质球对平面波电磁散射的解析；解

引言

在过去的几十年中，研究物体的电磁散射对射频识别有着重要的意义，其中对于各向异性的介质球的电磁散射研究一直比较热。各种数值方法被提出来，如时域有限差分法、有限元法、边界元法。由于解析法可以为其他数值计算提供比较有效的数据，可以对数值计算结果的正确性进行验证，成为许多学者的追求。文献首次用T矩阵法来研究旋磁介质球电磁散射，随后该方法被其他学者用来研究单轴媒介、旋电磁介质球的散射。

本文借助T矩阵法研究复杂介质球的电磁散射。文中的复杂介质是在各项异性媒介的本征方程的基础上增加一个新的张量，在本构式子中，电场磁场之间存在耦合。

1基公式及原理

在无源均匀的媒质中电磁场满足的麦克斯韦方程为(时间因子eiωt)：

图1所示的介质球的半径为a，球心位于原点。区域0为自由空间，介电常数和磁导率分别为£0和四。区域1中的介质本征结构方程为：

介质参数的各张量的表达式分别为:

将式(2)带入到麦克斯韦方程组式(1)并对其进行化简,最终可以得出下列方程：

1.1介质球内部的电磁场展开

根据矢量球波函数的性质以及磁场所满足的麦克斯韦方程，介质球内的磁场％可以展开为:

其中k为待定参数，瓦„=旺0CmnE，E°表示入射电场的场强。将式(4)中与Blnt相关的量转换为球矢量波函数的表达式并代入式(4)最后整理得到：

上式表达的含义是：存在这样的参数k使得方程有非零解。由矩阵的知识可知：只需令式(8)的行列式为0,解出参数k记解为k,(1=1,2,3…)再用k代入式(8),求出方程不为0的解［dmnl,cmn］-1即可。可以构造一新的矢量函数V:

球内部的磁感应强度Bmt可以表示为:

αl为待定系数，由介质球体表面的边界条件决定。这样磁场就可以得到：

1.2介质球外部入射场和散射场

散射场和入射场分别定义为：ES，HS和EI，HI利用边界条件可以得到：

这样就得到了散射场中的系数amn，bmn。式中的mS=kS/

2数值结果

找到了文献中这样的数值计算结果，便可计算一个特例。令张量0I=t，退化为各向异性的介质球。从图2可以看到，其与参考文献的数据有一个较好的匹配。从而说明了本方法和程序的正确性。

图3中张量的各个参数分别为εS=2ε0，μSμt=4μ0，εSεg=0.4ε0，εSεt=4ε0，μSμg=0.4μ0，μS=2μ0，ξS=0.3/c，ξt=ξg=0，设c是光在真空的速度且x=0.5π。这样，从图中可以看到，在大约50°的地方，在H面雷达散射截面RCS达到了最小的值。图4所示是显示张量It中的参数张量参数为μS=μ0，εSεg=0.6ε0，μg=0，εSεt=2ε0，ξ=0.3/c，ξt=ξg=0，x=0.6π。图5说明球的前向散射和后向散射各参数分别为：μS=μ0，εS=2ε0，μSμt=2μ0，εSεt=2ε0，εg=0，μSμg=0.5μ0，ξ=0.3/c，ξt=ξg=0，x=0.45π。

3结语

本文给出了三个张量的复杂媒介的球散射解析解，同时给出了一些算例。所增加的张量It可以推广到更一般的形式。这是从各向异性到双各项异性迈出的一小步。接下来需要对双各向异性介质球电磁散射的解析解开展深入研究。