主要关注这篇论文里关于物理光学法,几何光学法和多次散射的计算方法的内容。即文章的第六章。
赵晔. 海面与舰船目标电磁散射的建模方法研究[D].西安电子科技大学,2016.
第六章 目标和海面的复合电磁散射研究
一、研究现状分析
研究目标和粗糙面的耦合散射是一项非常具有实际应用价值的研究课题。
由于数值计算方法所剖分的面元数量过多,对算力要求较高,因此只适合频率较低情况下简单目标与粗糙面的复合散射,对于高频大范围的海船复合散射并不适用。
常用的方法是高频的近似。方法的优势在于可以高效率地处理超电大复合电磁建模问题,难点在于保证散射贡献地计算精度,并且兼顾计算效率。有如下方法:
- 传统四路径模型 将目标下方的粗糙面简化为无限大平面,然后利用镜面散射原理计算耦合散射贡献。只适合目标和微粗糙表面的复合散射计算。
- 迭代物理光学法(IPO) 能够分析目标和下方粗糙面所组成的角反射器结构的多次散射的高频近似方法。方法在物理光学电流中增加了修正电流项以近似导体目标表面因多次散射导致的实际电流。为了保证耦合散射的计算精度,IPO方法要求一定的迭代次数,并且利用消隐技术判断各个面元之间可见性也相当耗时,故处理超电大场景时依然存在计算效率低,容易发散等问题。
- 射线追踪理论(SBR) 主要根据高频电磁波的光学特性提出,能够适用于任何形状的复杂目标的电磁散射计算,但其计算的准确性依赖于射线管的数量,对于超电大目标与海面的耦合散射,依然存在计算效率低的问题。
本文结合海面面元散射场模型,通过构建海面与舰船目标复合散射的几何光学与物理光学(GO-PO)相结合的混合模型。
本文模型的优势在于GO-PO方法不仅可以用于计算复杂目标,还可以应用于目标和海面的耦合散射。方法首先将海面和舰船目标分别面元化,然后按照二次散射计算的过程,可以得到海面面元和目标上面元之间的耦合场。
二、三维目标电磁散射的GO-PO方法
1、面元对于入射波及一次反射波的可见度判断
(1)对于入射波
其中,
r
n
r_n
rn代表面元
n
n
n的坐标中心,
n
^
n
\hat{n}_n
n^n为面元法向矢量,面元
m
m
m同理。则设经过面元
m
m
m中心坐标的入射线为
r
(
t
)
=
r
m
+
k
^
i
t
r\left( t \right) =r_m+\hat{k}_it
r(t)=rm+k^it,设面元
n
n
n所在平面表示
n
^
n
(
r
−
r
n
)
=
0
\hat{n}_n\left( r-r_n \right) =0
n^n(r−rn)=0为该射线与面元
n
n
n所在平面的交点记为
P
(
x
0
,
y
0
,
z
0
)
P(x_0,y_0,z_0)
P(x0,y0,z0),满足如下公式:
{ x 0 = x m + k i ⋅ t 0 y 0 = y m + k i ⋅ t 0 z 0 = z m + k i ⋅ t 0 t 0 = n ^ n ⋅ ( r n − r m ) / ( n ^ n ⋅ k ^ i ) \left\{ \begin{array}{c} x_0=x_m+k_i\cdot t_0\\ y_0=y_m+k_i\cdot t_0\\ z_0=z_m+k_i\cdot t_0\\ \end{array} \right. \\ t_0=\hat{n}_n\cdot \left( r_n-r_m \right) /\left( \hat{n}_n\cdot \hat{k}_i \right) ⎩ ⎨ ⎧x0=xm+ki⋅t0y0=ym+ki⋅t0z0=zm+ki⋅t0t0=n^n⋅(rn−rm)/(n^n⋅k^i)
t
0
t_0
t0显然应该
t
0
<
0
t_0<0
t0<0,若交点
P
P
P在面元
n
n
n上,且入射射线和面元
m
m
m有交点,则说明有所遮挡。判断后项条件成立需要:
k
^
r
⋅
n
^
n
⩽
0
\boldsymbol{\hat{k}}_{\boldsymbol{r}}\cdot \hat{n}_n\leqslant 0
k^r⋅n^n⩽0
判断后项条件成立的方法是假设面元是如下三角型,并且赋予各顶点坐标:
则由几何知识可得,若
P
P
P在其内部,必同时满足如下三式:
{
(
r
P
N
1
×
r
N
1
N
2
)
⋅
n
^
n
>
0
(
r
P
N
2
×
r
N
2
N
3
)
⋅
n
^
n
>
0
(
r
P
N
3
×
r
N
3
N
1
)
⋅
n
^
n
>
0
\left\{ \begin{array}{c} \left( r_{PN_1}\times r_{N_1N_2} \right) \cdot \hat{n}_n>0\\ \left( r_{PN_2}\times r_{N_2N_3} \right) \cdot \hat{n}_n>0\\ \left( r_{PN_3}\times r_{N_3N_1} \right) \cdot \hat{n}_n>0\\ \end{array} \right.
⎩
⎨
⎧(rPN1×rN1N2)⋅n^n>0(rPN2×rN2N3)⋅n^n>0(rPN3×rN3N1)⋅n^n>0
只需要遍历所有的面元
n
n
n便能够判断面元
m
m
m是否受到了遮挡。算法流程图如下:
(2)对于一次反射波的可见度判断
分析方法与入射波类似。如果面元
m
m
m的一次反射波会被面元
n
n
n所遮挡,则需要满足如下条件:
(1)
k
^
r
,
m
⋅
n
^
n
⩽
0
\boldsymbol{\hat{k}}_{\boldsymbol{r},\boldsymbol{m}}\cdot \hat{n}_n\leqslant 0
k^r,m⋅n^n⩽0,即面元
m
m
m的反射射线能照射到面元
n
n
n所在的平面,同时,射线与面元
n
n
n有交点。
(2)反射波射线在到达面元
n
n
n之前不会受到其他面元的遮挡。其中条件(2)的判断方法如图所示:
在判断面元
m
m
m的反射波射线与面元
n
n
n是否有交点时,采用反向射线追踪法。即判断以面元
n
n
n的中心点
r
n
r_n
rn为起点,以
−
k
^
r
,
m
-\boldsymbol{\hat{k}}_{\boldsymbol{r},\boldsymbol{m}}
−k^r,m为方向的射线,是否与面元
n
n
n有交点。则若有交点,就保证了面元
n
n
n能被
m
m
m的反射波所照亮。将上述假设代入方程中,可以解得面元
n
n
n的反射射线与面元
m
m
m所在平面的交点
P
1
P_1
P1
{
x
1
=
x
m
−
k
r
x
,
m
⋅
t
1
y
1
=
y
m
−
k
r
y
,
m
⋅
t
1
z
1
=
z
m
−
k
r
z
,
m
⋅
t
1
t
1
=
n
^
m
⋅
(
r
n
−
r
m
)
/
(
n
^
m
⋅
k
^
r
,
m
)
\left\{ \begin{array}{c} x_1=x_m-k_{\boldsymbol{rx},\boldsymbol{m}}\cdot t_1\\ y_1=y_m-k_{\boldsymbol{ry},\boldsymbol{m}}\cdot t_1\\ z_1=z_m-k_{\boldsymbol{rz},\boldsymbol{m}}\cdot t_1\\ \end{array} \right. \\ t_1=\hat{n}_m\cdot \left( r_n-r_m \right) /\left( \hat{n}_m\cdot \hat{k}_{\boldsymbol{r},\boldsymbol{m}} \right)
⎩
⎨
⎧x1=xm−krx,m⋅t1y1=ym−kry,m⋅t1z1=zm−krz,m⋅t1t1=n^m⋅(rn−rm)/(n^m⋅k^r,m)
显然 t 1 > 0 t_1>0 t1>0,然后和入射波的方法类似可以判断点 P 1 P_1 P1是否在面元 m m m内。若满足之前所说的条件(1),下面对条件(2)进行判断。
相同的道理,对于面元 m m m和面元 l l l之间的关系,依然假设从面元 n n n上有一个方向为 k ^ m , r \hat{k}_{\boldsymbol{m},\boldsymbol{r}} k^m,r的射线,若该射线经过面元 l l l所在的平面且交点 P 2 P_2 P2在面元 l l l内部,则有所遮挡。
综上所述,得到判断面元对于一次反射波的可见度的流程图如下:
2、入射波感应的电磁流和一次散射的计算
可设入射波的电场强度矢量为:
E
i
=
p
^
E
0
exp
(
i
k
k
^
i
⋅
r
)
E^i=\hat{p}E_0\exp \left( ik\hat{k}_i\cdot r \right)
Ei=p^E0exp(ikk^i⋅r)其中
p
^
\hat{p}
p^为入射波单位极化矢量。可以将入射波的电磁场在本地坐标系中分解为水平和垂直两个分量,即:
E
h
i
=
h
^
i
′
(
p
^
⋅
h
^
i
′
)
E
0
H
h
i
=
k
^
i
×
[
h
^
i
′
(
p
^
⋅
h
^
i
′
)
E
0
]
/
η
=
−
v
^
i
′
(
p
^
⋅
h
^
i
′
)
E
0
/
η
E
v
i
=
v
^
i
′
(
p
^
⋅
v
^
i
′
)
E
0
H
v
i
=
−
h
^
i
′
(
p
^
⋅
v
^
i
′
)
E
0
/
η
E_{h}^{i}=\hat{h}_{i}^{'}\left( \hat{p}\cdot \hat{h}_{i}^{'} \right) E_0 \\ H_{h}^{i}\,\,=\,\,\hat{k}_i\times \left[ \hat{h}_{i}^{'}\left( \hat{p}\cdot \hat{h}_{i}^{'} \right) E_0 \right] /\eta =-\hat{v}_{i}^{'}\left( \hat{p}\cdot \hat{h}_{i}^{'} \right) E_0/\eta \\ E_{v}^{i}=\hat{v}_{i}^{'}\left( \hat{p}\cdot \hat{v}_{i}^{'} \right) E_0 \\ H_{v}^{i}\,\,=-\hat{h}_{i}^{'}\left( \hat{p}\cdot \hat{v}_{i}^{'} \right) E_0/\eta
Ehi=h^i′(p^⋅h^i′)E0Hhi=k^i×[h^i′(p^⋅h^i′)E0]/η=−v^i′(p^⋅h^i′)E0/ηEvi=v^i′(p^⋅v^i′)E0Hvi=−h^i′(p^⋅v^i′)E0/η
其中, h ^ i ′ = k ^ i × n ^ ∣ k ^ i × n ^ ∣ , v ^ i ′ = h ^ i ′ × k ^ i \hat{h}_{i}^{'}=\frac{\hat{k}_i\times \hat{n}}{\left| \hat{k}_i\times \hat{n} \right|},\hat{v}_{i}^{'}=\hat{h}_{i}^{'}\times \hat{k}_i h^i′=∣k^i×n^∣k^i×n^,v^i′=h^i′×k^i
因此,切向的水平极化场可以表示为:
n
^
×
E
h
=
n
^
×
(
E
h
i
+
E
h
r
)
=
n
^
×
E
h
i
(
1
+
R
h
)
n
^
×
H
h
=
n
^
×
(
H
h
i
+
H
h
r
)
=
n
^
×
(
k
^
i
⋅
E
h
i
+
k
^
r
⋅
E
h
r
R
h
)
/
η
=
−
[
(
n
^
⋅
k
^
i
)
E
h
i
+
(
n
^
⋅
k
^
r
)
E
h
r
R
h
]
/
η
=
−
(
n
^
⋅
k
^
i
)
(
1
−
R
h
)
E
h
i
/
η
\hat{n}\times E_h=\hat{n}\times \left( E_{h}^{i}+E_{h}^{r} \right) =\hat{n}\times E_{h}^{i}\left( 1+R_h \right) \\ \hat{n}\times H_h=\hat{n}\times \left( H_{h}^{i}+H_{h}^{r} \right) =\hat{n}\times \left( \hat{k}_i\cdot E_{h}^{i}+\hat{k}_r\cdot E_{h}^{r}R_h \right) /\eta \\ =-\left[ \left( \hat{n}\cdot \hat{k}_i \right) E_{h}^{i}+\left( \hat{n}\cdot \hat{k}_r \right) E_{h}^{r}R_h \right] /\eta \\ =-\left( \hat{n}\cdot \hat{k}_i \right) \left( 1-R_h \right) E_{h}^{i}/\eta
n^×Eh=n^×(Ehi+Ehr)=n^×Ehi(1+Rh)n^×Hh=n^×(Hhi+Hhr)=n^×(k^i⋅Ehi+k^r⋅EhrRh)/η=−[(n^⋅k^i)Ehi+(n^⋅k^r)EhrRh]/η=−(n^⋅k^i)(1−Rh)Ehi/η
同理可以求得切向垂直极化场:
n
^
×
H
v
=
n
^
×
H
v
i
(
1
+
R
v
)
n
^
×
E
v
=
(
n
^
⋅
k
^
i
)
⋅
(
1
−
R
h
)
H
v
i
η
\hat{n}\times H_v=\hat{n}\times H_{v}^{i}\left( 1+R_v \right) \\ \hat{n}\times E_v=\left( \hat{n}\cdot \hat{k}_i \right) \cdot \left( 1-R_h \right) H_{v}^{i}\eta
n^×Hv=n^×Hvi(1+Rv)n^×Ev=(n^⋅k^i)⋅(1−Rh)Hviη
在上述式中,
n
^
\hat{n}
n^为面元单位法相矢量,反射波的波矢量为
k
^
r
=
k
^
i
−
2
n
^
(
k
^
i
⋅
n
^
)
\hat{k}_r=\hat{k}_i-2\hat{n}\left( \hat{k}_i\cdot \hat{n} \right)
k^r=k^i−2n^(k^i⋅n^)
η
\eta
η为波阻抗,
R
h
,
R
v
R_h,R_v
Rh,Rv分别为水平和垂直极化的菲涅尔反射系数,对于理想导体,两者分别为-1,1.
有了水平和垂直极化的电磁切向场,就可以得到总的切向电磁场,根据麦克斯韦方程组就可以得到感应电流和感应磁流:
M
0
=
n
^
×
E
=
(
n
^
×
h
^
i
′
)
(
p
^
⋅
h
^
i
′
)
(
1
+
R
h
)
E
0
+
h
^
i
′
(
p
^
⋅
v
^
i
′
)
(
n
^
×
k
^
i
′
)
(
1
−
R
v
)
E
0
J
0
=
n
^
×
H
=
[
(
n
^
×
h
^
i
′
)
(
p
^
⋅
v
^
i
′
)
(
1
+
R
v
)
E
0
−
h
^
i
′
(
p
^
⋅
v
^
i
′
)
(
n
^
×
k
^
i
′
)
(
1
−
R
h
)
E
0
]
/
η
M_0\,\,=\,\,\hat{n}\times E=\left( \hat{n}\times \hat{h}_{i}^{'} \right) \left( \hat{p}\cdot \hat{h}_{i}^{'} \right) \left( 1+R_h \right) E_0+\hat{h}_{i}^{'}\left( \hat{p}\cdot \hat{v}_{i}^{'} \right) \left( \hat{n}\times \hat{k}_{i}^{'} \right) \left( 1-R_v \right) E_0 \\ J_0=\hat{n}\times H=\left[ \left( \hat{n}\times \hat{h}_{i}^{'} \right) \left( \hat{p}\cdot \hat{v}_{i}^{'} \right) \left( 1+R_v \right) E_0-\hat{h}_{i}^{'}\left( \hat{p}\cdot \hat{v}_{i}^{'} \right) \left( \hat{n}\times \hat{k}_{i}^{'} \right) \left( 1-R_h \right) E_0 \right] /\eta
M0=n^×E=(n^×h^i′)(p^⋅h^i′)(1+Rh)E0+h^i′(p^⋅v^i′)(n^×k^i′)(1−Rv)E0J0=n^×H=[(n^×h^i′)(p^⋅v^i′)(1+Rv)E0−h^i′(p^⋅v^i′)(n^×k^i′)(1−Rh)E0]/η
当面元
m
m
m能够被入射波直接照射,可以写出其远区散射场:
E
p
q
,
m
(
r
)
=
q
^
{
j
k
exp
(
j
k
R
)
4
π
R
k
^
s
×
∫
(
M
0
,
m
−
η
k
^
s
×
J
0
,
m
)
exp
[
j
k
(
k
^
i
−
k
^
s
)
⋅
r
m
′
]
d
s
′
}
=
exp
(
j
k
R
)
4
π
R
[
q
^
⋅
k
^
s
×
(
M
0
,
m
−
η
k
^
s
×
J
0
,
m
)
]
⋅
I
m
E_{pq,m}\left( r \right) =\hat{q}\left\{ \frac{jk\exp \left( jkR \right)}{4\pi R}\hat{k}_s\times \int{\left( M_{0,m}-\eta \hat{k}_s\times J_{0,m} \right) \exp \left[ jk\left( \hat{k}_i-\hat{k}_s \right) \cdot r_{m}^{'} \right] ds^{'}} \right\} \\ =\frac{\exp \left( jkR \right)}{4\pi R}\left[ \hat{q}\cdot \hat{k}_s\times \left( M_{0,m}-\eta \hat{k}_s\times J_{0,m} \right) \right] \cdot I_m
Epq,m(r)=q^{4πRjkexp(jkR)k^s×∫(M0,m−ηk^s×J0,m)exp[jk(k^i−k^s)⋅rm′]ds′}=4πRexp(jkR)[q^⋅k^s×(M0,m−ηk^s×J0,m)]⋅Im
其中, q ^ \hat{q} q^为散射波单位极化矢量,取 h ^ s \hat{h}_s h^s或 v ^ s \hat{v}_s v^s
对于相位积分项,根据Gordan方法,其具体的计算表达式为:
I
=
j
k
∫
exp
[
j
k
(
k
^
i
−
k
^
s
)
⋅
r
m
′
]
d
s
′
=
exp
(
j
k
r
0
⋅
ω
)
T
∑
m
=
1
M
(
p
^
⋅
a
m
)
exp
(
j
k
b
m
⋅
ω
)
sin
(
k
a
m
⋅
ω
/
2
)
k
a
m
⋅
ω
/
2
I=jk\int{\exp \left[ jk\left( \hat{k}_i-\hat{k}_s \right) \cdot r_{m}^{'} \right] ds^{'}} \\ =\frac{\exp \left( jkr_0\cdot \omega \right)}{T}\sum_{m=1}^M{\left( \hat{p}\cdot a_m \right) \exp \left( jkb_m\cdot \omega \right) \frac{\sin \left( ka_m\cdot \omega /2 \right)}{ka_m\cdot \omega /2}}
I=jk∫exp[jk(k^i−k^s)⋅rm′]ds′=Texp(jkr0⋅ω)m=1∑M(p^⋅am)exp(jkbm⋅ω)kam⋅ω/2sin(kam⋅ω/2)
其中, M M M为面元的边长个数; r 0 r_0 r0为面元的中心位置矢量, ω = k ^ i − k ^ s \omega=\hat{k}_i-\hat{k}_s ω=k^i−k^s, p ^ = n ^ × ω / ∣ n ^ × ω ∣ \hat{p}=\hat{n}\times \omega /\left| \hat{n}\times \omega \right| p^=n^×ω/∣n^×ω∣, a m = r m + 1 − r m , b m = ( r m + 1 + r m ) / 2 a_m = r_{m+1}-r_{m},b_m = \left(r_{m+1}+r_{m}\right)/2 am=rm+1−rm,bm=(rm+1+rm)/2, T T T为 ω \omega ω在面元上的投影长度。
综上所述,总的一次散射场可以表示为:
E
p
q
1
=
∑
m
=
1
N
E
p
q
,
m
(
r
m
)
⋅
I
v
i
s
,
m
E_{pq}^{1}=\sum_{m=1}^N{E_{pq,m}\left( r_m \right) \cdot I_{vis,m}}
Epq1=m=1∑NEpq,m(rm)⋅Ivis,m
其中,N表示目标上总的面元数, I v i s , m I_{vis,m} Ivis,m为面元 m m m对于入射波的可见度因子,能够直射则为1,否则为0
3.一次反射波感应的电磁流和二次散射的计算
显然,采取的方法是将面元
m
m
m的反射波作为入射波分析直射面元
n
n
n时的产生的感应电磁流和二次散射电磁场。当然,要想产生二次散射波,必须
I
v
i
s
,
n
m
I_{vis,nm}
Ivis,nm为1,这样可以得到总的二次散射场:
E
p
q
2
=
exp
(
j
k
R
)
4
π
R
∑
n
=
1
N
[
q
^
⋅
k
^
s
×
(
∑
m
=
1
N
I
v
i
s
,
m
⋅
M
1
,
n
m
⋅
I
v
i
s
,
n
m
−
η
k
^
s
×
∑
m
=
1
N
I
v
i
s
,
m
⋅
M
1
,
n
m
⋅
I
v
i
s
,
n
m
)
]
⋅
I
n
E_{pq}^{2}= \\ \frac{\exp \left( jkR \right)}{4\pi R}\sum_{n=1}^N{\left[ \hat{q}\cdot \hat{k}_s\times \left( \sum_{m=1}^N{I_{vis,m}\cdot M_{1,nm}\cdot I_{vis,nm}-\eta \hat{k}_s\times \sum_{m=1}^N{I_{vis,m}\cdot M_{1,nm}\cdot I_{vis,nm}}} \right) \right]}\cdot I_n
Epq2=4πRexp(jkR)n=1∑N[q^⋅k^s×(m=1∑NIvis,m⋅M1,nm⋅Ivis,nm−ηk^s×m=1∑NIvis,m⋅M1,nm⋅Ivis,nm)]⋅In
4. 目标总散射场方程
根据2,3,某个面元
m
m
m上的感应电流和感应磁流既有来自入射波的直接影响,也有其他面元的二次反射波所导致。因此目标
T
a
r
g
e
t
Target
Target的总散射场可以写出如下形式:
E
p
q
t
a
r
g
e
t
(
k
^
i
,
k
^
s
)
=
∑
n
=
1
N
exp
(
j
k
R
)
4
π
R
[
q
^
⋅
k
^
s
×
(
M
n
−
η
k
^
s
×
J
n
)
]
⋅
I
n
E_{pq}^{target}\left( \hat{k}_i,\hat{k}_s \right) =\sum_{n=1}^N{\frac{\exp \left( jkR \right)}{4\pi R}\left[ \begin{array}{c} \hat{q}\cdot \hat{k}_s\times \left( M_n-\eta \hat{k}_s\times J_n \right)\\ \end{array} \right]}\cdot I_n
Epqtarget(k^i,k^s)=n=1∑N4πRexp(jkR)[q^⋅k^s×(Mn−ηk^s×Jn)]⋅In