太赫兹波段碲化镉介电常数的理论与实验研究

孙旺; 李粮生; 张景; 殷红成

doi:10.12000/JR17096

太赫兹波段碲化镉介电常数的理论与实验研究

DOI: 10.12000/JR17096

电磁散射重点实验室北京 100854

基金项目: 国家自然科学基金(6149069502)

详细信息

作者简介:
孙　旺(1993–)，男，硕士研究生，研究方向为太赫兹波段材料的密度泛函方法

李粮生(1981–)，男，研究员，研究方向为太赫兹技术、量子技术和统计物理

张　景(1987–)，男，博士后，研究方向为太赫兹测量技术

殷红成(1967–)，男，研究员，博士生导师，主要研究方向为电磁散射、雷达目标特性、目标识别等

通讯作者:
李粮生 liliangshengbititp@163.com

中图分类号: TN952
计量
- 文章访问数: 2390
- HTML全文浏览量: 556
- PDF下载量: 477
- 被引次数: 11
出版历程
- 收稿日期: 2017-11-03
- 修回日期: 2018-01-21
- 网络出版日期: 2018-02-01

Theoretical and Experimental Study on the Permittivity of CdTe in the Terahertz Band

Science and Technology on Electromagnetic Scattering Laboratory, Beijing 100854, China

Funds: The National Natural Science Foundation of China (6149069502)

摘要

摘要: 该文采用密度泛函理论研究了闪锌矿型碲化镉的声子色散谱、特征向量及晶格振动频率，获得了碲化镉介电常数随频率变化的理论值。通过太赫兹时域光谱系统测量碲化镉单晶的介电常数，该实验结果与局域密度近似修正、广义梯度近似修正和广义梯度近似修正的计算结果吻合较好。最后，3种近似交换关联势的计算结果之间存在一定差异性，该差异性结果表明太赫兹波段碲化镉的介电常数由电子声子耦合所主导，但是横波和纵波声子频率敏感于电子密度分布。
- 碲化镉 /
- 介电常数 /
- 反射率 /
- 横/纵光学模 /
- 时域光谱
Abstract: The phonon dispersion spectrum, eigenvector, and lattice vibration frequency of cadmium telluride with a zinc blende structure have been investigated using the density functional theory, and the permittivity of cadmium telluride crystal is numerically calculated. The permittivity of the crystal is measured using the terahertz time-domain spectroscopy system. The experimental results are consistent with the theoretical calculations on the modified local density approximation, the general gradient approximation, and the modified general gradient approximation. Finally, the differences among the three approximate exchange correlation potentials indicate that in the terahertz region, the permittivity of cadmium telluride is dominantly contributed by the coupling between electron and phonon; however, the phonon frequencies of transverse wave and longitudinal wave were sensitive to electron density distribution.
- Cadmium telluride /
- Permittivity /
- Reflectivity /
- Transverse/longitudinal optical mode /
- Time domain spectroscopy

HTML全文

1. 引言

混合源定位在无源雷达中具有广泛的应用前景^[1,2]。按照空间中的辐射源与阵列之间的相对距离，可以将辐射源分为远场源和近场源。对远场源定位，需要对波达方向(Direction Of Arrival, DOA)进行估计；对近场源定位，除了对波达方向进行估计，还需要对距离参数进行估计。阵列的近场区域小于 $2{D^2}/\lambda$ ，其中D为圆阵孔径， $\lambda$ 为波长^[3,4]。虽然远场源可以视为距离为无穷远的近场源，但是当远场源的二维DOA与近场源二维DOA相同时，如果采用近场和远场参数联合估计将无法对混合源进行识别。此外，由于远场源定位场景和近场源定位场景可以看作混合源定位场景特殊模式，因此混合源定位方法同样适用于辐射源全部为远场源定位场景或者辐射源全部为近场源定位场景。

不同于单独对远场源进行定位的方法和单独对近场源进行定位的方法，对混合源进行定位还需要对近场源和远场源进行分类和识别。文献[5]首先利用距离为无穷远时的近场源导向矢量和混合源的噪声子空间，通过一维多重信号分类(One Dimensional MUltiple SIgnal Classification, 1-D MUSIC)方法估计出远场源DOA，接着利用对称阵元的2阶统计量去除近场源的距离参数，估计出混合源DOA并通过斜投影算法对混合源进行分类，最后将近场源DOA代入后导向矢量和混合源噪声子空间，通过1-D MUSIC方法对近场源距离进行估计。文献[6]提出了基于子空间的混合源定位算法，该算法利用线性变换估计出远场源的DOA，通过正交三角分解(QR分解)得到斜投影算子避免了特征值分解，大大减少了计算复杂度，利用1-D MUSIC方法确定近场源DOA，最后将近场源DOA代入导向矢量，确定近场源的距离参数。该方法有效解决了在快拍数有限时出现饱和的情况。文献[7]利用两级对称嵌套线阵，利用距离为无穷远的近场源导向矢量和混合源的噪声子空间，通过1-D MUSIC方法估计出远场源DOA，接着利用斜投影方法去除混合源接收信号中的远场分量，计算近场源在对称阵元下的4阶累积量并进行向量化处理，估计出近场源DOA，最后利用代入近场源DOA的导向矢量和混合源噪声子空间，确定近场源的距离参数。该方法能够在阵元数相同的情况下增大阵列孔径，提高混合源参数估计精度。上述3种方法在得到远场源的波达方向的基础上，利用斜投影算法将远场源和近场源进行分离，但是在远场源参数估计的基础上利用斜投影算法得到的近场源分量会受到噪声分量的影响，引起近场源参数估计误差，导致定位精度下降。文献[8]利用协方差矩阵差分方法消除混合源中远场源和噪声分量，得到近场源差分矩阵，该方法能够有效避免噪声分量对近场源定位精度的影响，文献[9]采用降秩的方法从混合源中分离出远场源分量，文献[10]设计了3个4阶累积量矩阵消除近场源导向矢量中距离参数，进而同时估计出远场源和近场源的波达方向。但是上述方法采用线阵对混合源进行定位只能估计出一维波达方向，并且线阵在0°波达方向的分辨率最高，在60°波达方向的分辨率会降低1/2。

相较于均匀线阵只能对一维波达方向进行估计，均匀圆阵可以实现方位角和俯仰角等二维波达方向估计^[11-13]，并且在不同方位角的分辨率具有各向同性^[14,15]。在现有文献中，由于均匀圆阵的导向矩阵不具有均匀线阵导向矩阵的Vandermonde形式，不便于在数学上进行处理，因此提出的方法相对较少。均匀圆阵下对混合源进行定位需要针对相应的数学模型进行分析，文献[16]采用均匀圆阵，提出了TSMUSIC方法，该方法利用均匀圆阵下的远场源的协方差矩阵具有Hermitian和Toeplitz形式，而近场源的协方差矩阵只具有Hermitian形式，利用传播算子的正交性得到近场源的噪声子空间，进而通过2-D MUSIC方法得到近场源二维DOA的空间谱，进而根据谱峰确定近场源二维DOA；接着，利用将近场源二维DOA代入后的导向矢量和混合源的噪声子空间，通过1-D MUSIC方法确定近场源的距离；最后，利用远场源的导向矢量和混合源的噪声子空间确定远场源的二维DOA。该方法参数估计精度较高，但是计算复杂度较大。文献[17]采用均匀圆阵，利用每个阵元的频谱计算非相干的混合源在每个阵元的相位，利用均匀圆阵中心对称的特点，通过计算对角阵元的相位差来去除近场源的距离参数，得到混合源二维DOA的相位差矩阵，进而运用最小二乘法同时估计出远场源和近场源的方位角与俯仰角；接着，根据1-D MUSIC方法得到距离空间谱的收敛性，对混合源进行识别并估计出近场源的距离参数，当距离空间谱出现峰值，该混合源判断为近场源，峰值所对应的位置即为近场源的距离；若距离空间谱不收敛，则该混合源判断为远场源。基于相位差的方法能够有效减少计算复杂度，但是混合源的定位精度有所降低。

针对上述方法存在的问题，本文根据均匀圆阵中心对称的特点，提出一种基于协方差矩阵差分的混合源定位方法，该方法首先利用二维多重信号分类方法估计出远场源的方位角和俯仰角；接着利用协方差矩阵差分方法提取出近场源差分矩阵，通过改进类旋转不变估计信号参数(Estimation of Signal Parameters via Rotational Invariance Techniques like, ESPRIT-like)方法计算出近场源的方位角和俯仰角；最后，利用一维多重信号分类方法估计出近场源的距离。该方法能够消除混合源协方差矩阵中的远场源的噪声分量，从而提取出近场源差分矩阵，并且相较于基于相位差的方法，本文提出的方法能够提高混合源的定位精度。当远场源的二维DOA与近场源二维DOA相同时，如果采用近场和远场参数联合估计将无法对混合源进行识别，采用本文所提的两步法可以将远场源和近场源进行有效分离，在混合源识别上更具优势。

2. 混合源模型

均匀圆阵下的混合源模型如图1所示，以均匀圆阵的圆心为原点建立三维坐标系，均匀圆阵的半径为 $R$ ，由分布在xy平面上的M个阵元组成，混合源包含K₁个远场源和K₂个近场源，方位角参数 $\phi$ 为混合源与均匀圆阵中心的连线投影到xy平面上，相对于x坐标轴的逆时针方向旋转的角度；俯仰角参数 $\theta$ 为混合源与均匀圆阵中心的连线相对于z坐标轴旋转的角度；当混合源为近场源时，距离参数 $r$ 为近场源与均匀圆阵中心的距离。

图 1 混合源模型

Figure 1. The geometric model of mixed sources

下载: 全尺寸图片幻灯片

均匀圆阵的第 $m$ 个阵元在第 $n$ 个快拍的接收数据可以表示为

${x_m}(n) = \sum\limits_{k = 1}^{{K_1}{\rm{ + }}{K_2}} {{s_k}(n){{\rm{e}}^{{\rm{j}}{\tau _{k,m}}}}} + {w_m}(n)$

(1)

其中， ${s_k}(n)$ 表示第 $k$ 个混合源发射数据， ${w_m}(n)$ 表示零均值复白高斯过程， ${\tau _{k,m}}$ 表示第 $k$ 个混合源到达阵列中心以及到达第 $m$ 个阵元的时差，当混合源为远场源时，时差可以表示为

${\tau _{k,m}}=\frac{{2{\rm{\pi }}}}{\lambda }{\zeta _{k,m}}$

(2)

其中， $k = 1,2, \cdots ,{K_1}$ ， $\lambda$ 表示波长， $2{\rm{\pi /}}\lambda$ 表示波数， ${\zeta _{k,m}} = \cos ({\varphi _m} - {\phi _k})\sin {\theta _k}$ $,\;{\varphi _m} = 2{\rm{\pi }}(m - 1)/ M$ ，当混合源为近场源时，时差可以表示为

${\tau _{k,m}}=\frac{{2{\rm{\pi }}}}{\lambda }\left( {{r_k}{\rm{ - }}{r_{k,m}}} \right)$

(3)

其中， $k = {K_1} + 1,{K_1} + 2, \cdots ,{K_1} + {K_2}$ , ${r_{k,m}}$ 表示第 $k$ 个近场源与第 $m$ 个阵元的距离， ${r_{k,m}}$ 通过2阶泰勒级数^[11]可以展开为

$\begin{split}{r_{k,m}}& = \sqrt {\mathop {{r_k^2 }} + \mathop R\nolimits^2 - 2{r_k}R{\zeta _{k,m}}} \\ & ={r_k} - R{\zeta _{k,m}} + \frac{{{R^2}}}{{2{r_k}}}(1 - \zeta _{k,m}^2) + O\left( {\frac{{{R^2}}}{{r_k^2}}} \right)\end{split}$

(4)

将式(2)和式(4)代入式(1)，混合源模型可以简化为

$\begin{split} {x_m}(n) = &\sum\limits_{k = 1}^{{K_1}} {{s_k}(n){{\rm{e}} ^{{\rm{j}}{\eta _{k,m}}}}} \\ &+ \sum\limits_{k = {K_1} + 1}^{{K_1} + {K_2}} {{s_k}(n){{\rm{e}} ^{{\rm{j}}{\eta _{k,m}}{\rm{ - j}}{\xi _{k,m}}}}} + {w_m}(n) \end{split}$

(5)

其中， ${\eta _{k,m}} = \left( {2{\rm{\pi }}R/\lambda } \right){\zeta _{k,m}}$ , ${\xi _{k,m}} = \left( {{\rm{\pi }}{R^2}/\lambda {r_k}} \right)$ $\cdot \left( {1 - \zeta _{k,m}^2} \right)$ 。

3. 算法描述

在均匀圆阵下的近场和远场混合源数学模型的基础上，本节首先利用2-D MUSIC方法估计出远场源的方位角和俯仰角，实现远场源定位；接着利用协方差矩阵差分方法提取出近场源差分矩阵，进而通过改进的类旋转不变估计信号参数方法计算出近场源的方位角和俯仰角；最后利用1-D MUSIC方法估计出近场源的距离，实现近场源定位。

3.1 远场源二维DOA参数估计

均匀圆阵下混合源的协方差矩阵 ${\boldsymbol{E}}$ 可以由式(6)计算

${\boldsymbol{E}}=\frac{1}{N}{\boldsymbol{X}}{{\boldsymbol{X}}^{\rm{H}}}$

(6)

其中， ${\boldsymbol{X}}=\left[ {{\boldsymbol{x}}(1)\;{\boldsymbol{x}}(2)\; \cdots \;{\boldsymbol{x}}(N)} \right]$ 表示 $M \times N$ 维的混合源接收数据， ${\boldsymbol{x}}(n) = {\left[ {{x_1}(n)\;{x_2}(n)\; \cdots \;{x_M}(n)} \right]^{\rm{T}}}$ 表示第 $n$ 个快拍接收数据， ${( \cdot )^{\rm{H}}}$ 表示共轭转置，N为混合源的快拍数， ${( \cdot )^{\rm{T}}}$ 表示转置运算。

对混合源的协方差矩阵 ${\boldsymbol{E}}$ 进行特征值分解得到

${\boldsymbol{E}}={{\boldsymbol{U}}_{\rm{s}}}{{\boldsymbol{\varSigma}} _{\rm{s}}}{\boldsymbol{U}}_{\rm{s}}^{\rm{H}} + {{\boldsymbol{U}}_{\rm{n}}}{{\boldsymbol{\varSigma}} _{\rm{n}}}{\boldsymbol{U}}_{\rm{n}}^{\rm{H}}$

(7)

其中， ${{\boldsymbol{\varSigma}}_{\rm{s}}} \in {\mathbb{C}^{({K_1} + {K_2}) \times ({K_1} + {K_2})}}$ 表示 ${K_1} + {K_2}$ 个大特征值组成的对角矩阵， ${{\boldsymbol{\varSigma}}_{\rm{n}}} \in {\mathbb{C}^{(M - {K_1} - {K_2}) \times (M - {K_1} - {K_2})}}$ 表示 $M - {K_1} - {K_2}$ 个小特征值组成的对角矩阵， ${{\boldsymbol{U}}_{\rm{s}}}\in$ ${\mathbb{C}^{M \times ({K_1} + {K_2})}}$ 表示 ${K_1} + {K_2}$ 个大特征值对应的特征向量组成的混合源信号子空间， ${{\boldsymbol{U}}_{\rm{n}}} \in {\mathbb{C}^{M \times (M - {K_1} - {K_2})}}$ 表示 $M - {K_1} - {K_2}$ 个小特征值对应的特征向量组成的混合源噪声子空间。

利用3-D MUSIC方法得到混合源的空间谱函数

${\boldsymbol{V}}\left( {\phi ,\theta ,r} \right) = \frac{1}{{{\boldsymbol{a}}_{{\rm{NF}}}^{\rm{H}}\left( {\phi ,\theta ,r} \right){{\boldsymbol{U}}_{\rm{n}}}{\boldsymbol{U}}_{\rm{n}}^{\rm{H}}{{\boldsymbol{a}}_{{\rm{NF}}}}\left( {\phi ,\theta ,r} \right)}}$

(8)

其中， ${{\boldsymbol{a}} _{{\rm{NF}}}}(\phi ,\theta ,r) = {[{a_{{\rm{NF,1}}}}(\phi ,\theta ,r)\;\;{a_{{\rm{NF}},2}}(\phi ,\theta ,r)\; \cdots}$ ${a_{{\rm{NF}},M}} (\phi ,\theta ,r)]^{\rm{T}}$ 为近场源的导向矢量， ${a_{{\rm{NF}},m}} \left( {\phi ,\theta ,r} \right) =$ ${{\rm{e}}^{{\rm{j}}(2{\rm{\pi }}R/\lambda ){\zeta _m}(\phi ,\theta ) - {\rm{j}}\left( {{\rm{\pi }}{R^2}/\lambda r} \right)\left( {1 - \zeta _m^2(\phi ,\theta )} \right)}}$ ， $m =$ $1,2, \cdots ,M$ 。混合源空间谱函数 ${\boldsymbol{V}}\left( {\phi ,\theta ,r} \right)$ 会出现 $K$ 个峰值，其中 ${K_1}$ 个峰值的距离估计 $\hat r \to \infty$ ，表示 ${K_1}$ 个远场源； ${K_2}$ 个峰值的距离估计 $r \in \left[ {0.62{{(4{R^2}/\lambda )}^{0.5}},8{R^2}/\lambda } \right]$ ，表示 ${K_2}$ 个近场源。因此，远场源二维DOA的空间谱函数可以通过2-D MUSIC方法得到

${{\boldsymbol{V}}_{{\rm{FF}}}}\left( {\phi ,\theta } \right) = \frac{1}{{{\boldsymbol{a}}_{{\rm{NF}}}^{\rm{H}}\left( {\phi ,\theta ,\infty } \right){{\boldsymbol{U}}_{\rm{n}}}{\boldsymbol{U}}_{\rm{n}}^{\rm{H}}{{\boldsymbol{a}}_{{\rm{NF}}}}\left( {\phi ,\theta ,\infty } \right)}}$

(9)

其中， ${{\boldsymbol{a}} _{{\rm{NF}}}}(\phi ,\theta ,\infty ) = {[{a_{{\rm{NF,1}}}}(\phi ,\theta ,\infty )\;{a_{{\rm{NF}},2}}(\phi ,\theta ,\infty )\; \cdots}$ ${\;{a_{{\rm{NF}},M}}(\phi ,\theta ,\infty )]^{\rm{T}}}$ 表示距离为无穷远情况下的近场源导向矢量，即远场源导向矢量，此时 ${a_{{\rm{NF}},m}}( {\phi ,\theta ,} {\infty } )$ $= {{\rm{e}}^{{\rm{j}}(2{\rm{\pi }}R/\lambda ){\zeta _m}(\phi ,\theta )}}$ , $m = 1,2, \cdots ,M$ 。只包含方位角和俯仰角等二维DOA。对远场源方位角和俯仰角的空间谱函数 ${{\boldsymbol{V}}_{{\rm{FF}}}}\left( {\phi ,\theta } \right)$ 进行谱峰搜索，峰值所对应的位置即为远场源的方位角估计值和俯仰角估计值，则第 $k$ 个远场源的定位结果为 $({\hat \phi _k},{\hat \theta _k},\infty )$ , $k = 1,2, \cdots ,{K_1}$ 。

3.2 混合源协方差矩阵差分

混合源协方差矩阵 ${\boldsymbol{E}}$ 可以表示为

$\begin{split} {\boldsymbol{E}}=&{{\boldsymbol{E}}_{{\rm{FF}}}}+{{\boldsymbol{E}}_{{\rm{NF}}}}+{\sigma ^2}{{\boldsymbol{I}}_M} \\ =&{{\boldsymbol{A}}_{{\rm{FF}}}}{{\boldsymbol{\varPi}} _{{\rm{FF}}}}{\boldsymbol{A}}_{{\rm{FF}}}^{\rm{H}}+{{\boldsymbol{A}}_{{\rm{NF}}}}{{\boldsymbol{\varPi}} _{{\rm{NF}}}}{\boldsymbol{A}}_{{\rm{NF}}}^{\rm{H}}+{\sigma ^2}{{\boldsymbol{I}}_M}\end{split}$

(10)

其中， ${{\boldsymbol{\varPi}} _{{\rm{FF}}}} \in {\mathbb{C}^{{K_1} \times {K_1}}}$ 表示远场源功率组成的对角矩阵， ${{\boldsymbol{\varPi}} _{{\rm{NF}}}} \in {\mathbb{C}^{{K_2} \times {K_2}}}$ 表示近场源功率组成的对角矩阵， ${\sigma ^2}$ 表示噪声功率， ${{\boldsymbol{I}}_M}$ 表示 $M \times M$ 维的单位矩阵， ${{\boldsymbol{A}}_{{\rm{FF}}}}$ 表示远场源的导向矩阵

${{\boldsymbol{A}}_{{\rm{FF}}}}=\big[{{\boldsymbol{a}}_{{\rm{FF}}}}({\phi _1},{\theta _1})\;{{\boldsymbol{a}}_{{\rm{FF}}}}({\phi _2},{\theta _2})\; \cdots \;{{\boldsymbol{a}}_{{\rm{FF}}}}({\phi _{{K_1}}},{\theta _{{K_1}}})\big]$

(11)

由于均匀圆阵具有中心对称的结构，满足 ${\varphi _{m + M/2}} = {\varphi _m} + {\rm{\pi }}$ 和 ${\zeta _{m + M/2}} = - {\zeta _m}$ ，因此，远场源的导向矢量可以表示为

$\begin{split} {{\boldsymbol{a}}_{{\rm{FF}}}}({\phi _k},{\theta _k}) = &[\underbrace {{{\rm{e}}^{{\rm{j}}{\eta _{k,1}}}}\;{{\rm{e}}^{{\rm{j}}{\eta _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{j}}{\eta _{k,M/2}}}}}_{M/2}\;\\ &\underbrace {{{\rm{e}}^{{\rm{ - j}}{\eta _{k,1}}}}\;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,M/2}}}}}_{M/2}{]^{\rm{T}}} \end{split}$

(12)

同理， ${{\boldsymbol{A}}_{{\rm{NF}}}}$ 表示近场源的导向矩阵为

$\begin{split} {{\boldsymbol{A}}_{{\rm{NF}}}}=& [{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _1},{\theta _1},{r_1})\;{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _2},{\theta _2},{r_2})\; \cdots\\ &\;{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})]\end{split}$

(13)

其中，近场源的导向矢量可以表示为

$\begin{split} {{\boldsymbol{a}}_{{\rm{NF}}}}&({\phi _k},{\theta _k},{r_k}) \\ =&[\underbrace {{{\rm{e}}^{{\rm{j}}{\eta _{k,1}}{\rm{ - j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{j}}{\eta _{k,2}}{\rm{ - j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{j}}{\eta _{k,M/2}}{\rm{ - j}}{\xi _{k,M/2}}}}}_{M/2}\;\\ &\underbrace {{{\rm{e}}^{{\rm{ - j}}{\eta _{k,1}}{\rm{ - j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,2}}{\rm{ - j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,M/2}}{\rm{ - j}}{\xi _{k,M/2}}}}}_{M/2}{]^{\rm{T}}}\end{split}$

(14)

由于远场源的导向矢量中第 $m$ 个元素和第 $m +$ $M/2$ 个元素具有共轭的结构，因此远场源的协方差矩阵 ${{\boldsymbol{E}}_{{\rm{FF}}}}$ 可以表示为

${\boldsymbol{E}}_{{\rm{FF}}}^{\rm{T}}={\boldsymbol{J}}{{\boldsymbol{E}}_{{\rm{FF}}}}{\boldsymbol{J}}$

(15)

其中，旋转变换矩阵 ${\boldsymbol{J}}$ 可以表示为

${\boldsymbol{J}} = \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{O}}_{M/2}}}&{{{\boldsymbol{I}}_{M/2}}} \\ {{{\boldsymbol{I}}_{M/2}}}&{{{\boldsymbol{O}}_{M/2}}} \end{array}} \right]$

(16)

${{\boldsymbol{O}}_{M/2}}$ 和 ${{\boldsymbol{I}}_{M/2}}$ 分别为 $M/2 \times M/2$ 维的零矩阵和单位矩阵。

由式(10)和式(15)可以推导近场源差分矩阵为

${{\boldsymbol{E}}_{{\rm{NF,diff}}}}={\boldsymbol{JEJ}}-{{\boldsymbol{E}}^{\rm{T}}}={\boldsymbol{J}}{{\boldsymbol{E}}_{{\rm{NF}}}}{\boldsymbol{J}}-{\boldsymbol{E}}_{{\rm{NF}}}^{\rm{T}}$

(17)

由于近场源协方差矩阵 ${{\boldsymbol{E}}_{{\rm{NF,diff}}}}$ 是Hermitian矩阵，满足 ${\boldsymbol{E}}_{{\rm{NF,diff}}}^{\rm{T}}={\boldsymbol{E}}_{{\rm{NF,diff}}}^{*}$ ，因此，近场源协方差矩阵 ${{\boldsymbol{E}}_{{\rm{NF,diff}}}}$ 可以重构为

$\begin{split} {{\boldsymbol{E}}_{{\rm{NF,diff}}}}=&{{\boldsymbol{JE}}_{{\rm{NF}}}}{\boldsymbol{J}}-{\boldsymbol{E}}_{{\rm{NF}}}^{*} \\ =&{{\boldsymbol{JA}}_{{\rm{NF}}}}{{\boldsymbol{\varPi}} _{{\rm{NF}}}}{\boldsymbol{A}}_{{\rm{NF}}}^{\rm{H}}{\boldsymbol{J}}-{\boldsymbol{A}}_{{\rm{NF}}}^{*}{\boldsymbol{\varPi}} _{{\rm{NF}}}^{*}{\left( {{\boldsymbol{A}}_{{\rm{NF}}}^{\rm{H}}} \right)^*} \\ = &\left[ {\begin{split} {{{\boldsymbol{JA}}_{{\rm{NF}}}}}\qquad{{\boldsymbol{A}}_{{\rm{NF}}}^{*}} \end{split}} \right]\left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{\varPi}} _{{\rm{NF}}}}}&{\rm{}}\\ {{\rm{}}}& {-{\boldsymbol{\varPi}} _{{\rm{NF}}}^{*}} \end{array}} \right] \\ & \cdot{\left[ {\begin{split} {{{\boldsymbol{JA}}_{{\rm{NF}}}}}\qquad{{\boldsymbol{A}}_{{\rm{NF}}}^{*}} \end{split}} \right]^{\rm{H}}}\\ =& {{\boldsymbol{A}}_{{\rm{NF,diff}}}}{{\boldsymbol{\varPi}} _{{\rm{NF}}}}{\boldsymbol{A}}_{{\rm{NF,diff}}}^{\rm{H}} \\ \end{split}$

(18)

其中，近场源差分矩阵的导向矩阵可以表示为

$\begin{split} {{\boldsymbol{A}}_{{\rm{NF}},}}{\rm{_{diff}}} =& [\underbrace {{\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _1},{\theta _1},{r_1})\;{\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _2},{\theta _2},{r_2})\; \cdots \;{\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})}_{{K_2}}\\ & \underbrace {{\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _1},{\theta _1},{r_1})\;{\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _2},{\theta _2},{r_2})\; \cdots \;{\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})}_{{K_2}}] \end{split}$

(19)

其中，第 $k$ 个导向矢量可以表示为

${\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _k},{\theta _k},{r_k}) =[\underbrace {{{\rm{e}}^{{\rm{ - j}}{\eta _{k,1}}{\rm{ - j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,2}}{\rm{ - j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,M/2}}{\rm{ - j}}{\xi _{k,M/2}}}}}_{M/2}\underbrace {{{\rm{e}}^{{\rm{j}}{\eta _{k,1}}{\rm{ - j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{j}}{\eta _{k,2}}{\rm{ - j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{j}}{\eta _{k,M/2}}{\rm{ - j}}{\xi _{k,M/2}}}}}_{M/2}{]^{\rm{T}}}$

(20)

第 $k + {K_2}$ 个导向矢量可以表示为

${\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _k},{\theta _k},{r_k}) =[\underbrace {{{\rm{e}}^{{\rm{ - j}}{\eta _{k,1}}{\rm{ + j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,2}}{\rm{ + j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{ - j}}{\eta _{k,M/2}}{\rm{ + j}}{\xi _{k,M/2}}}}}_{M/2}\underbrace {{{\rm{e}}^{{\rm{j}}{\eta _{k,1}}{\rm{ + j}}{\xi _{k,1}}}}\;{{\rm{e}}^{{\rm{j}}{\eta _{k,2}}{\rm{ + j}}{\xi _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{j}}{\eta _{k,M/2}}{\rm{ + j}}{\xi _{k,M/2}}}}}_{M/2}{]^{\rm{T}}}$

(21)

由于 ${{\boldsymbol{A}}_{{\rm{NF,diff}}}}$ 的第 $k$ 列的指数项第1部分与第 $k + {K_1}$ 列的指数项第1部分相同，各列中第 $m$ 个元素的指数项第2部分与第 $m + M/2$ 个元素的指数项第2部分相同，因此可以得到

$\begin{split} {{\boldsymbol{JA}}_{{\rm{NF,diff}}}}=[\underbrace {{\boldsymbol{D}}({\phi _1},{\theta _1}){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _1},{\theta _1},{r_1})\;\;{\boldsymbol{D}}({\phi _2},{\theta _2}){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _2},{\theta _2},{r_2})\;\; \cdots \;\;{\boldsymbol{D}}({\phi _{{K_2}}},{\theta _{{K_2}}}){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})}_{{K_2}} \\ \underbrace { {\boldsymbol{D}}({\phi _1},{\theta _1}){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _1},{\theta _1},{r_1})\;\;{\boldsymbol{D}}({\phi _2},{\theta _2}){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _2},{\theta _2},{r_2})\;\; \cdots \;\;{\boldsymbol{D}}({\phi _{{K_2}}},{\theta _{{K_2}}}){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})]}_{{K_2}} \end{split}$

(22)

其中， ${\boldsymbol{D}}({\phi _k},{\theta _k})$ 表示为

$\begin{split} {\boldsymbol{D}}({\phi _k},{\theta _k}) =& {\rm{diag}}\underbrace {[{{\rm{e}}^{{\rm{j}}2{\eta _{k,1}}}}\;{{\rm{e}}^{{\rm{j}}2{\eta _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{j}}2{\eta _{k,M/2}}}}}_{M{\rm{/2}}}\\ &\underbrace {{{\rm{e}}^{{\rm{ - j}}2{\eta _{k,1}}}}\;{{\rm{e}}^{{\rm{ - j}}2{\eta _{k,2}}}}\; \cdots \;{{\rm{e}}^{{\rm{ - j}}2{\eta _{k,M/2}}}}]}_{M{\rm{/2}}} \end{split}$

(23)

${\boldsymbol{D}}({\phi _k},{\theta _k})$ 只包含近场源的二维DOA参数，因此可以利用 ${\boldsymbol{D}}({\phi _k},{\theta _k})$ 将近场源三维位置参数估计分解为二维DOA估计和距离参数估计，能够有效降低计算复杂度。

3.3 近场源二维DOA参数估计

对仅包含近场源差分矩阵 ${{\boldsymbol{E}}_{{\rm{NF,diff}}}}$ 进行特征值分解

$\begin{split} {{\boldsymbol{E}}_{{\rm{NF,diff}}}}=&{{\boldsymbol{U}}_{{\rm{NF,diff,s}}}}{{\boldsymbol{\varSigma}}_{{\rm{NF,diff,s}}}}{\boldsymbol{U}}_{{\rm{NF,diff,s}}}^{\rm{H}} \\ &+{{\boldsymbol{U}}_{{\rm{NF,diff,n}}}}{{\boldsymbol{\varSigma}}_{{\rm{NF,diff,n}}}}{\boldsymbol{U}}_{{\rm{NF,diff,n}}}^{\rm{H}}\end{split}$

(24)

其中， ${{\boldsymbol{\varSigma}}_{{\rm{NF,diff}},{\rm{s}}}} \in {\mathbb{C}^{2{K_2} \times 2{K_2}}}$ 表示 $2{K_2}$ 个大特征值组成的对角矩阵， ${{\boldsymbol{\varSigma}}_{{\rm{NF,diff}},{\rm{n}}}} \in {\mathbb{C}^{(M - 2{K_2}) \times (M - 2{K_2})}}$ 表示 $M - 2{K_2}$ 个小特征值组成的对角矩阵， ${{\boldsymbol{U}}_{{\rm{NF,diff,s}}}} \in$ ${\mathbb{C}^{M \times 2{K_2}}}$ 表示 $2{K_2}$ 个大特征值对应的特征向量组成的近场源信号子空间， ${{\boldsymbol{U}}_{\rm{n}}} \in {\mathbb{C}^{M \times (M - 2{K_2})}}$ 表示M– $2{K_2}$ 个小特征值对应的特征向量组成的近场源噪声子空间。

存在一个满秩矩阵 ${\boldsymbol{G}}$ 使得

${{\boldsymbol{U}}_{{\rm{NF,diff,s}}}}={{\boldsymbol{A}}_{{\rm{NF,diff}}}}{\boldsymbol{G}}$

(25)

由式(22)和式(19)可以得到

$\begin{split} {{\boldsymbol{JU}}}&_{{{\rm{NF,diff,s}}}} - {\boldsymbol{\varPsi}} (\phi ,\theta ){{\boldsymbol{U}}_{{\rm{NF,diff,s}}}}\\ &={{\boldsymbol{JA}}_{{\rm{NF,diff}}}}{\boldsymbol{G}} - {\boldsymbol{\varPsi}} (\phi ,\theta ){{\boldsymbol{A}}_{{\rm{NF,diff}}}}{\boldsymbol{G}}=\left[ {{{\boldsymbol{\varOmega}} _1}\;{{\boldsymbol{\varOmega}} _2}} \right]{\boldsymbol{G}}\end{split}$

(26)

其中，

$\begin{split} {{\boldsymbol{\varOmega}} _1} = &[ {( {{\boldsymbol{D}}({\phi _1},{\theta _1}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} ){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _1},{\theta _1},{r_1})\;}\\ &\left( {{\boldsymbol{D}}({\phi _2},{\theta _2}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} \right){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _2},{\theta _2},{r_2})\; \cdots\\ &\left( {{\boldsymbol{D}}({\phi _{{K_2}}},{\theta _{{K_2}}}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} \right){\boldsymbol{J}}{{\boldsymbol{a}}_{{\rm{NF}}}}({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}})]\end{split}$

(27)

$\begin{split}{{\boldsymbol{\varOmega}} _2} =& [ {\left( {{\boldsymbol{D}}({\phi _1},{\theta _1}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} \right){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _1},{\theta _1},{r_1})\;}\\ &\left( {{\boldsymbol{D}}({\phi _2},{\theta _2}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} \right){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _2},{\theta _2},{r_2})\; \cdots \\ &\;\left( {{\boldsymbol{D}}({\phi _{{K_2}}},{\theta _{{K_2}}}) - {\boldsymbol{\varPsi}} (\phi ,\theta )} \right){\boldsymbol{a}}_{{\rm{NF}}}^*({\phi _{{K_2}}},{\theta _{{K_2}}},{r_{{K_2}}}) ]\end{split}$

(28)

当 ${\boldsymbol{D}}({\phi _k},{\theta _k}) = {\boldsymbol{\varPsi}} (\phi ,\theta )$ 时， ${{\boldsymbol{JU}}_{{\rm{NF,diff,s}}}} - {\boldsymbol{\varPsi}} (\phi ,\theta )$ ${{\boldsymbol{U}}_{{\rm{NF,diff,s}}}}=0$ 。因此，可以通过改进的ESPRIT-like方法^[15]得到近场源二维DOA的空间谱函数

${{\boldsymbol{V}}_{{\rm{NF}}}}\left( {\phi ,\theta } \right) = \frac{1}{{\det \left( {{{\boldsymbol{Q}}^{\rm{H}}}{{\boldsymbol{JU}}_{{\rm{NF,diff,s}}}} - {{\boldsymbol{Q}}^{\rm{H}}}{\boldsymbol{\varPsi}} (\phi ,\theta ){{\boldsymbol{U}}_{{\rm{NF,diff,s}}}}} \right)}}$

(29)

其中， $\det \left( \cdot \right)$ 表示行列式的值， ${\boldsymbol{Q}}$ 为 $M \times 2{K_2}$ 维随机满秩矩阵， ${\boldsymbol{J}}$ 为3.2节确定的 $M \times M$ 维旋转变换矩阵，对角矩阵 ${\boldsymbol{\varPsi}} (\phi ,\theta )$ 为

$\begin{split} {\boldsymbol{\varPsi}} (\phi ,\theta ) = {\rm{diag([}}\underbrace {{{\rm{e}} ^{{\rm{j}}2{\eta _1}(\phi ,\theta )}}\;{{\rm{e}} ^{{\rm{j}}2{\eta _2}(\phi ,\theta )}}\; \cdots \;{{\rm{e}} ^{{\rm{j}}2{\eta _{M/2}}(\phi ,\theta )}}}_{M/2}\;\\ \underbrace {{{\rm{e}} ^{{\rm{ - j}}2{\eta _1}(\phi ,\theta )}}\;{{\rm{e}} ^{{\rm{ - j}}2{\eta _2}(\phi ,\theta )}}\; \cdots \;{{\rm{e}} ^{{\rm{ - j}}2{\eta _{M/2}}(\phi ,\theta )}}}_{M/2}{\rm{])}}\end{split}$

(30)

${\rm{diag}}\left( \cdot \right)$ 表示将向量转化为对角矩阵运算。对近场源方位角和俯仰角的空间谱函数 ${{\boldsymbol{V}}_{{\rm{NF}}}}\left( {\phi ,\theta } \right)$ 进行谱峰搜索，峰值所对应的位置即为近场源的方位角估计值 ${\hat \phi _k}$ 和俯仰角估计值 ${\hat \theta _k}$ , $k = 1,2, \cdots ,{K_2}$ 。

3.4 近场源距离参数估计

在估计出近场源的方位角和俯仰角的基础上，通过1-D MUSIC方法得到第 $k$ 个近场源的距离空间谱函数

$\begin{split} &{{\boldsymbol{G}}_{{\rm{NF}}}}({\hat \phi _k},{\hat \theta _k},{r_k}) =\\ &\quad\frac{1}{{{{\boldsymbol{a}}_{{\rm{NF}}}}\left({{\hat \phi }_k},{{\hat \theta }_k},{r_k}\right){{\boldsymbol{U}}_{\rm{n}}}{\boldsymbol{U}}_{\rm{n}}^{\rm{H}}{{\boldsymbol{a}}_{{\rm{NF}}}}\left({{\hat \phi }_k},{{\hat \theta }_k},{r_k}\right)}} \end{split}$

(31)

其中， ${{\boldsymbol{U}}_{\rm{n}}}$ 为混合源噪声子空间， ${{\boldsymbol{a}}_{{\rm{NF}}}}({\hat \phi _k},{\hat \theta _k},{r_k})$ 表示第 $k$ 个近场源的方位角估计值 ${\hat \phi _k}$ 和俯仰角估计值 ${\hat \theta _k}$ 代入近场源导向矢量。对第 $k$ 个近场源的距离空间谱函数 ${{\boldsymbol{G}}_{{\rm{NF}}}}({\hat \phi _k},{\hat \theta _k},{r_k})$ 进行谱峰搜索，峰值位置即为第 $k$ 个近场源相对于阵列中心的距离估计值 ${\hat r_k}$ ，则第 $k$ 个近场源的定位结果为 $({\hat \phi _k},{\hat \theta _k},{\hat r_k})$ , $k =$ $1,2, \cdots ,{K_2}$ 。

4. 算法分析

4.1 计算复杂度分析

对于均匀圆阵下的近场源定位，传统的3-D MUSIC方法通过三维谱峰搜索估计出近场源方位角、俯仰角和距离，本文首先通过改进的ESPRIT-like方法得到近场源的二维DOA空间谱，利用二维谱峰搜索估计出近场源的方位角和俯仰角，进而利用1-D MUSIC方法得到近场源的距离空间谱，利用一维谱峰搜索估计出近场源的距离。本文提出的方法需要 ${M^2}N$ 次乘法和 ${M^2}(N - 1)$ 次加法计算混合源协方差矩阵，需要 $2{M^2}N$ 次乘法和 ${M^2}$ 次减法计算近场源差分矩阵，需要 $O({M^2})$ 量级的计算复杂度分别对混合源协方差矩阵以及近场源差分矩阵进行特征值分解，需要 $O({M^2})$ 量级的计算复杂度得到远场源二维DOA空间谱，需要 $4{K_2}{M^2} + 8K_2^2M$ 次乘法、 $(4{K_2}M + 8K_2^2)(M - 1)$ 次加法和 $4K_2^2$ 次减法得到近场源二维DOA空间谱，需要 $O(M)$ 量级的计算复杂度得到近场源距离空间谱。

4.2 混合源能量差异敏感性分析

由于本文所提的两步法通过协方差矩阵差分去除混合源的远场源和噪声分量，进而提取出近场源差分矩阵，因此对近场源进行估计时不受远场源和噪声能量的影响。本文第1步在对远场源进行估计时，当远场源的能量远远小于近场源的能量，此时通过第1步的方法对远场源进行定位的性能急剧恶化。为了解决两个源能量差异敏感的问题，可以首先对噪声能量进行估计，去除混合源协方差矩阵中的噪声分量；利用第2步的近场源估计结果，通过斜投影算法去除无噪混合源协方差矩阵中近场源分量，进而得到远场源分量。具体步骤如下：

混合源的协方差矩阵可以分解为

${\boldsymbol{R}} = \mathop {\left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{R}}_{11}}}& {{{\boldsymbol{R}}_{12}}}\\ {{{\boldsymbol{R}}_{21}}}& {{{\boldsymbol{R}}_{22}}} \end{array}} \right]}\limits^{K\qquad M - K} {\rm{ }}_{\scriptstyle\hfill\atop \scriptstyle M - K\hfill}^K$

(32)

噪声的能量可以通过式(33)计算

$\hat {\boldsymbol{\sigma}} _{\rm{n}}^2=\frac{{{\rm{tr}}\left\{ {{{\boldsymbol{R}}_{{\rm{22}}}}{\boldsymbol Z}} \right\}}}{{{\rm{tr}}\left\{ {\boldsymbol Z} \right\}}}$

(33)

其中， ${\boldsymbol Z}={{\boldsymbol{I}}_{M - K}}-{{\boldsymbol{R}}_{{\rm{21}}}}{\boldsymbol{R}}_{{\rm{21}}}^\dagger$ , ${\boldsymbol{R}}_{{\rm{21}}}^\dagger ={\left( {{\boldsymbol{R}}_{{\rm{21}}}^{\rm{H}}{{\boldsymbol{R}}_{{\rm{21}}}}} \right)^{ - 1}} {\boldsymbol{R}}_{{\rm{21}}}^{\rm{H}}$ ，(·)^†表示伪逆运算，tr(·)表示矩阵的迹， ${{\boldsymbol{I}}_{M - K}}$ 表示(M–K)维单位矩阵。因此，无噪的混合源协方差矩阵可以表示为

$\bar {\boldsymbol{R}}={{\boldsymbol{R}}_{{\rm{NF}}}} + {{\boldsymbol{R}}_{{\rm{FF}}}} = {\boldsymbol{R}} - \hat {\boldsymbol\sigma} _{\rm{n}}^{\rm{2}}{{\boldsymbol{I}}_M}$

(34)

由于斜投影矩阵 ${{\boldsymbol{E}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}$ 满足

${{\boldsymbol{E}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}{{\boldsymbol{A}}_{{\rm{NF}}}}={{\boldsymbol{A}}_{{\rm{NF}}}},\;{{\boldsymbol{E}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}{{\boldsymbol{A}}_{{\rm{FF}}}}={{\boldsymbol{O}}_{M \times (K - {K_1})}}$

(35)

并可以通过式(36)计算

${\hat {\boldsymbol{E}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}={\hat {\boldsymbol{A}}_{{\rm{NF}}}}{\left( {\hat {\boldsymbol{A}}_{{\rm{NF}}}^{\bf{H}}{{\bar {\boldsymbol{R}}}^\dagger }{{\hat {\boldsymbol{A}}}_{{\rm{NF}}}}} \right)^{ - 1}}\hat {\boldsymbol{A}}_{{\rm{NF}}}^{\bf{H}}{\bar {\boldsymbol{R}}^\dagger }$

(36)

其中， ${{\boldsymbol{O}}_{M \times (K - {K_1})}}$ 表示 $M \times (K - {K_1})$ 维零矩阵， ${\hat {\boldsymbol{A}}_{{\rm{NF}}}}$ 表示将近场源参数估计结果代入后的导向矩阵， ${\bar {\boldsymbol{R}}^\dagger }={\left( {{{\bar {\boldsymbol{R}}}^{\rm{H}}}\bar {\boldsymbol{R}}} \right)^{ - 1}}{\bar {\boldsymbol{R}}^{\rm{H}}}$ 。因此，远场源的接收矩阵可以表示为

$\begin{split} {{\boldsymbol{X}}_{{\rm{FF}}}}={{\boldsymbol{A}}_{{\rm{FF}}}}{{\boldsymbol{S}}_{{\rm{FF}}}}=&\left( {{{\boldsymbol{I}}_M}-{{\hat {\boldsymbol{E}}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}} \right) \\ &\cdot \left( {{{\boldsymbol{A}}_{{\rm{NF}}}}{{\boldsymbol{S}}_{{\rm{NF}}}} + {{\boldsymbol{A}}_{{\rm{FF}}}}{{\boldsymbol{S}}_{{\rm{FF}}}}} \right)\end{split}$

(37)

进而远场源的协方差矩阵可以通过式(38)计算

$\begin{split} {\hat {\boldsymbol{R}}_{{\rm{FF}}}}={{\boldsymbol{X}}_{{\rm{FF}}}}{\boldsymbol{X}}_{{\rm{FF}}}^{\rm{H}}=&\left( {{{\boldsymbol{I}}_M}-{{\hat {\boldsymbol{E}}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}} \right) \\ &\cdot \bar {\boldsymbol{R}}{\left( {{{\boldsymbol{I}}_M}-{{\hat {\boldsymbol{E}}}_{{{\boldsymbol{A}}_{{\rm{NF}}}}\left| {{{\boldsymbol{A}}_{{\rm{FF}}}}} \right.}}} \right)^{\rm{H}}}\end{split}$

(38)

此时，可以减少两个源的能量差异对远场源二维DOA估计的影响。

4.3 算法总结

基于协方差矩阵差分的混合源定位方法流程图如图2所示，对于远场源定位，本文利用均匀圆阵下的混合源接收数据得到混合源协方差矩阵，利用远场源导向矢量以及经过特征值分解得到的混合源噪声子空间，通过2-D MUSIC方法估计出远场源的方位角和俯仰角；对于近场源定位，本文利用混合源协方差矩阵和旋转矩阵得到近场源差分矩阵，利用随机满秩矩阵、对角矩阵以及经过特征值分解的近场源信号子空间，通过改进的ESPRIT-like方法计算出近场源的方位角和俯仰角，利用近场源的导向矢量和混合源的噪声子空间，通过1-D MUSIC方法估计出近场源的距离。本文根据均匀圆阵下的远场源在中心对称阵元的导向参数具有共轭性质，利用协方差矩阵差分方法去除混合源的远场源分量，得到近场源差分矩阵。由于协方差矩阵差分方法没有利用第1步的远场源二维DOA参数估计结果，因此，第2步近场源参数估计性能与远场源DOA估计性能无关。或者说，第1步和第2步的顺序也可以进行互换，可以先通过混合源协方差矩阵差分方法得到近场源差分矩阵，估计出近场源位置参数，再利用2-D MUSIC方法估计远场源DOA参数。另外，在多个混合源能量相同的场景下，由于第1步利用远场源导向矢量与噪声子空间的正交性，因此不会将近场源判断为远场源。当近场源的距离临近 $2{D^2}/\lambda$ ，由于噪声和模型误差的影响，第1步可能会将近场源的方向判断为远场源。当远场源的二维DOA与近场源二维DOA相同时，如果采用近场和远场参数联合估计将无法对混合源进行识别，采用本文所提的两步法可以将远场源和近场源进行有效分离，在混合源识别上更具优势。

图 2 基于协方差矩阵差分的混合源定位方法流程图

Figure 2. Flow chart of the mixed source localization based on the covariance differencing method

下载: 全尺寸图片幻灯片

5. 仿真实验

本节通过MATLAB仿真验证提出方法对混合源的定位效果，仿真中均匀圆阵的阵元个数为8，半径为0.5 m，以均匀圆阵的中心为原点建立三维坐标系，混合源包含一个远场源和一个近场源，远场源的位置为(50°, 30°, $\infty$ )，近场源的位置为(55°, 35°, 10 m)。

5.1 混合源定位有效性

这一小节用于验证本文提出方法得到的空间谱对远场和近场混合源的定位效果，在信噪比为20 dB，快拍数为600时，图3为通过本文提出方法得到的混合源定位的空间谱，图3(a)为通过2-D MUSIC方法得到的远场源二维DOA空间谱，方位角估计范围为0.1°～360°，间隔为0.1°，俯仰角估计范围为0°～90°，间隔0.1°，可以看出二维空间谱中出现的一个峰值对应的位置为(50.0°, 30.2°)，即为远场源方位角和俯仰角估计。图3(b)为通过改进的ESPRIT-like方法得到的近场源二维DOA空间谱，同样地，方位角估计范围为0.1°～360°，间隔为0.1°，俯仰角估计范围为0°～90°，间隔0.1°，可以看出二维空间谱中出现的一个峰值对应的位置为(55.0°, 35.1°)，即为近场源方位角和俯仰角估计。图3(c)为通过1-D MUSIC方法得到的近场源一维距离空间谱，近场源距离估计范围为0.1～30.0 m，间隔为0.1 m，图中红色的线表示近场源真实距离，可以看出距离空间谱中出现的一个峰值对应的位置为9.9 m，即为近场源距离估计。综上可知，本文所提协方差矩阵差分方法能够从混合源中提取出近场源分量，并通过改进的ESPRIT-like方法近场源三维位置参数中分离出二维DOA，最终实现混合源的有效定位。

图 3 混合源空间谱

Figure 3. Spatial spectrum of mixed sources

下载: 全尺寸图片幻灯片

5.2 不同信噪比下混合源定位精度

本节验证所提方法在不同信噪比下对混合源的定位性能，并与TSMUSIC方法^[16]、相位差方法^[17]以及克拉美罗界(Cramer-Rao Lower Bound, CRLB)进行对比，其中TSMUSIC方法是利用2-D MUSIC方法分别估计出远场源和近场源的二维DOA，利用1-D MUSIC方法分别估计出近场源的距离；相位差方法是利用均匀圆阵下对角阵元的相位差，通过最小二乘法同时估计出远场源和近场源的二维DOA，根据1-D MUSIC方法对混合源进行识别并估计出近场源的距离参数；CRLB是参数任意无偏估计量的下限。在仿真中利用均方根误差(Root Mean Square Error，RMSE)对算法的性能进行分析，RMSE可由式(39)计算

${\rm{RMSE}}\left( \beta \right) = \sqrt {\frac{1}{L}\sum\limits_{l = 1}^L {{{\left( {{{\hat \beta }_l} - \beta } \right)}^2}} }$

(39)

其中，L表示独立重复仿真次数， ${\hat \beta _l}$ 表示第 $l$ 次的混合源参数估计结果， $l = 1,2, ··· ,L$ , $\beta$ 表示混合源的真实位置参数。当快拍数为600时，图4(a)—图4(e)分别为通过600次独立重复仿真得到在不同信噪比下的远场源方位角RMSE、远场源俯仰角RMSE、近场源方位角RMSE、近场源俯仰角RMSE以及近场源距离RMSE，从图4可以看出混合源的RMSE随着信噪比的增大而降低，表明混合源参数估计精度不断提高。另外，还可以看出本文所提基于协方差矩阵差分方法的混合源定位精度比基于相位差方法的定位精度高。由于本文方法与TSMUSIC方法都采用2-D MUSIC方法对远场源进行定位，因此两种方法的远场源定位精度相当，本文所提改进ESPRIT-like方法对近场源二维DOA参数估计精度略低于TSMUSIC方法参数估计精度。此外，由于近场源距离估计是在近场源二维DOA估计基础上进行的，相较于近场源二维DOA估计误差，近场源距离估计误差相对较大，从图4中也可以看出，在信噪比为20 dB时，近场源的二维DOA估计误差接近0.01°，而近场源的距离误差接近于0.1 m。

图 4 不同信噪比下混合源RMSE

Figure 4. RMSE of mixed sources versus SNR

下载: 全尺寸图片幻灯片

5.3 不同快拍数下混合源定位精度

本节进一步验证在不同快拍数下混合源的定位性能，本文所提方法同样与TSMUSIC方法、相位差方法以及克拉美罗界进行了对比。当信噪比为20 dB时，图5(a)—图5(e)分别为通过600次独立重复仿真得到在不同快拍数下的远场源方位角RMSE、远场源俯仰角RMSE、近场源方位角RMSE、近场源俯仰角RMSE以及近场源距离RMSE，从图5中可以看出混合源RMSE随着快拍数的增大而降低，表明混合源参数估计精度逐渐提高。另外，从图中可以看出当快拍数小于600时，混合源RMSE快速下降，当快拍数超过600时，混合源RMSE逐渐平缓。此外，从图5可以看出本文提出方法的混合源定位精度比基于相位差方法的定位精度高，对于远场源参数估计精度，本文提出的方法与TSMUSIC方法相当，对于近场源参数估计精度，本文提出的方法略低于TSMUSIC方法。

图 5 不同快拍数下的混合源RMSE

Figure 5. RMSE of mixed sources versus snapshots

下载: 全尺寸图片幻灯片

5.4 混合源能量差异性对比

由于本文所提的两步法通过协方差矩阵差分去除混合源的远场源和噪声分量，进而提取出近场源差分矩阵，因此对近场源进行估计时不受远场源和噪声能量的影响。图6(a)—图6(c)给出了当近场源和远场源能量比(Near-field and Far-field source Power Ratio, NFPR)分别为1.0, 0.2和0.1时的近场源空间谱，可以看出近场源空间谱出现一个峰值，峰值所对应的位置即为近场源二维DOA。仿真结果表明，通过本文提出的算法对近场源进行估计时不受远场源和噪声能量的影响。

图 6 近场源空间谱

Figure 6. Spatial spectrum of near-field source

下载: 全尺寸图片幻灯片

为了解决远场源对能量差异敏感的问题，通过对信噪比和近场源进行估计，去除混合源协方差矩阵中的噪声分量和近场源分量，得到式(38)远场源协方差矩阵，进而估计远场源二维DOA。图7(a)—图7(c)给出了当近场源和远场源能量比分别为1, 5和10时的近场源空间谱，可以看出远场源空间谱出现了一个峰值，峰值所对应的位置即为远场源二维DOA。仿真结果表明，通过本文提出的算法对近场源进行估计时不受远场源和噪声能量的影响。仿真结果表明，所提算法能够有效地解决远场源能量差异敏感的问题。

图 7 远场源空间谱

Figure 7. Spatial spectrum of far-field source

下载: 全尺寸图片幻灯片

6. 结论

本文所提均匀圆阵下利用混合源协方差矩阵差分的方法能够消除混合源中的近场源和噪声分量，提取出近场源差分矩阵，实现远场源和近场源的有效分离，利用2-D MUSIC方法估计出远场源的方位角和俯仰角，通过改进的ESPRIT-like方法提取并估计出近场源的方位角和俯仰角，运用1-D MUSIC方法估计出近场源的距离参数。仿真结果表明，相较于基于相位差的混合源定位算法，本文所提方法能够提高混合源的定位精度。由于协方差矩阵差分方法没有利用第1步的远场源二维DOA参数估计结果，因此，第2步近场源参数估计性能与远场源DOA估计性能无关。或者说，第1步和第2步的顺序也可以进行互换，可以先通过混合源协方差矩阵差分方法得到近场源差分矩阵，估计出近场源位置参数，再利用2-D MUSIC方法估计远场源DOA参数。另外，在多个混合源能量相同的场景下，由于第1步利用远场源导向矢量与噪声子空间的正交性，因此不会将近场源判断为远场源。当近场源的距离临近 $2{D^2}/\lambda$ ，由于噪声和模型误差的影响，第1步可能会将近场源的方向判断为远场源。当远场源的二维DOA与近场源二维DOA相同时，如果采用近场和远场参数联合估计将无法对混合源进行识别，采用本文所提的两步法可以将远场源和近场源进行有效分离，在混合源识别上更具优势。

图 1 CdTe晶体结构：红球为Te原子，蓝球为Cd原子，绿色箭头表示晶体的主轴方向

Figure 1. The crystal structure of CdTe: The red spheres indicate Te atoms. The blue spheres indicate Cd atoms. The green arrow indicates the direction of crystal major axis

下载: 全尺寸图片幻灯片

图 2 总能量与截断能的依赖关系。黑色箭头表示选取GGA的截断能位置 $E_{\rm{c}}^{{\rm{PBE/PW91}}}{\rm{ = }}680 \ {{\rm{eV}}}$ 。红色箭头表示选取LDA的截断能位置 $E_{\rm{c}}^{{\rm{LDA}}}{\rm{ = }}1225\ {{\rm{eV}}}$

Figure 2. The total energy vs the cut-off energy. The black arrow indicates the value of GGA cut-off energy $E_{\rm{c}}^{{\rm{PBE/PW91}}}{\rm{ = }}680 \ {{\rm{eV}}}$ . The red arrow indicates the value of LDA cut-off energy $E_{\rm{c}}^{{\rm{LDA}}}{\rm{ = }}1225\ {\rm eV}$

下载: 全尺寸图片幻灯片

图 3 总能量随k-point数量的依赖关系。黑色箭头表示合适的k-point选择

Figure 3. The total energy as a function of k-point. Black arrows indicate the selected k-point

下载: 全尺寸图片幻灯片

图 4 LDA, GGA-PBE和GGA-PW91方法计算的总能量与尝试晶格常数a的依赖关系。黑色箭头表示LDA计算总能量最小值的位置 $a_{{\rm{LDA}}}^{{\rm{CA {^{\scriptsize{-}}} PZ}}} = 0.64424 \ {\rm{nm}}$

Figure 4. The total energy (calculated by LDA, GGA-PBE, and GGA-PW91 methods) vs the test lattice constant. The black arrow indicates the minimum of total energy, where $a_{{\rm{LDA}}}^{{\rm{CA {^{\scriptsize{-}}} PZ}}} = 0.64424 \ {\rm{nm}}$

下载: 全尺寸图片幻灯片

图 5 CdTe的声子色散谱，使用3种交换关联势为LDA, PBE, PW91

Figure 5. Phonon dispersions of CdTe with three exchange correlations LDA, PBE, PW91

下载: 全尺寸图片幻灯片

图 6 声子波矢方向与原子振动方向平行(a), (d)与垂直(b), (c)示意图。(a)和(b)中的声子波矢方向为CdTe主轴方向。(c)和(d)中的声子波矢方向为y轴方向。黑色箭头表示原子振动方向。绿色箭头表示声子传播方向

Figure 6. Schematic diagram: the phonon wave vectors are parallel (a), (d) and perpendicular (b), (c) to directions of atomic vibration. (a)-(b) the phonon wave vectors are parallel to the direction of CdTe principal axis. (c)-(d) The phonon wave vectors are parallel to the y-axis direction. Black arrows indicate the direction of atomic vibration. Green arrows indicate the direction of phonon propagation

下载: 全尺寸图片幻灯片

图 7 全空气太赫兹系统实验装置示意图。β-BBO为I型硼酸钡，PM为离轴抛物面镜，HV为高电压调制器，PMT为光电倍增管

Figure 7. Schematic diagram of full-air terahertz experimental system. β-BBO is type I barium borate, PM is an off-axis parabolic mirror, HV is a high-voltage modulator, and PMT is a photomultiplier tube

下载: 全尺寸图片幻灯片

图 8 理论计算与实验测量结果

Figure 8. The results of theoretical calculations and experimental measurements

下载: 全尺寸图片幻灯片

表 1 CdTe单晶材料： ${\text{ω} _{\text{LO}}}$ , ${\text{ω} _{\text{TO}}}$ , $\text{ε}$ (∞)和 $\text{γ}$ 理论和实验值

Table 1. CdTe single crystal: the theoretical and experimental values of ${\text{ω} _{{\text{LO}}}}$ , ${\text{ω} _{{\text{TO}}}}$ , $\text{ε}$ (∞) and $\text{γ}$

方法	光学声子横波模式 ${\omega _{{\rm{TO}}}}$ (THz)	光学声子纵波模式 ${\omega _{{\rm{LO}}}}$ (THz)	高频极限介电常数 $\varepsilon (\infty )$	碰撞频率 $\gamma$ (THz)	晶格常数 ${a_0}$ (nm)
LDA	4.402	5.111	7.55792	0.13	0.64424
GGA-PBE	4.518	5.301	6.62839	0.12	0.66594
GGA-PW91	4.533	5.309	6.72645	0.12	0.66541
实验值	4.25	5.01	6.93	0.125	0.64827^[33]

下载: 导出CSV

参考文献(33)

[1]	Lloyd-Hughes J and Jeon T I. A review of the terahertz conductivity of bulk and Nano-materials[J]. Journal of Infrared,Millimeter,and Terahertz Waves, 2012, 33(9): 871–925. DOI: 10.1007/s10762-012-9905-y.
[2]	程伟, 王迎新, 赵自然. 光电导太赫兹源新进展[J]. 激光与红外, 2011, 41(6): 597–604. DOI: 10.3969/j.issn.1001-5078.2011.06.001. Cheng Wei, Wang Ying-xin, and Zhao Zi-ran. New research progress of photoconductive terahertz source[J]. Laser&Infrared, 2011, 41(6): 597–604. DOI: 10.3969/j.issn.1001-5078.2011.06.001.
[3]	Withayachumnankul W, Png G M, Yin X X, et al. T-Ray sensing and imaging[J]. Proceedings of the IEEE, 2007, 95(8): 1528–1558. DOI: 10.1109/JPROC.2007.900325.
[4]	Scalari G, Maissen C, Turčinková D, et al. Ultrastrong coupling of the cyclotron transition of a 2D electron gas to a THz metamaterial[J]. Science, 2012, 335(6074): 1323–1326. DOI: 10.1126/science.1216022.
[5]	Stantchev R I, Sun B Q, Hornett S M, et al. Noninvasive, near-field terahertz imaging of hidden objects using a single-pixel detector[J]. Science Advances, 2016, 2(6): e1600190. DOI: 10.1126/sciadv.1600190.
[6]	Dengler R J, Cooper K B, Chattopadhyay G, et al.. 600 GHz imaging radar with 2 cm range resolution[C]. The IEEE/MTT-S International Microwave Symposium, Honolulu, HI, USA, 2007: 1371–1374.
[7]	Cooper K B, Dengler R J, Llombart N, et al. THz imaging radar for standoff personnel screening[J]. IEEE Transactions on Terahertz Science and Technology, 2011, 1(1): 169–182. DOI: 10.1109/TTHZ.2011.2159556.
[8]	Cooper K B, Dengler R J, Llombart N, et al. Penetrating 3-D imaging at 4- and 25-m range using a submillimeter-wave radar[J]. IEEE Transactions on Microwave Theory and Techniques, 2008, 56(12): 2771–2778. DOI: 10.1109/TMTT.2008.2007081.
[9]	Cooper K B, Dengler R J, Chattopadhyay G, et al. A high-resolution imaging radar at 580 GHz[J]. IEEE Microwave and Wireless Components Letters, 2008, 18(1): 64–66. DOI: 10.1109/LMWC.2007.912049.
[10]	Llombart N, Cooper K B, Dengler R J, et al. Confocal ellipsoidal reflector system for a mechanically scanned active terahertz imager[J]. IEEE Transactions on Antennas and Propagation, 2010, 58(6): 1834–1841. DOI: 10.1109/TAP.2010.2046860.
[11]	Sheen D M, Hall T E, Severtsen R H, et al.. Standoff concealed weapon detection using a 350-GHz radar imaging system[C]. SPIE Passive Millimeter-Wave Imaging Technology XIII, Orlando, Florida, United States, 2010, 7670: 767008.
[12]	Sheen D M, McMakin D L, Hall T E, et al.. Active millimeter-wave standoff and portal imaging techniques for personnel screening[C]. 2009 IEEE Conference on Technologies for Homeland Security, Boston, MA, USA, 2009: 440–447.
[13]	Robertson D A, Marsh P N, Bolton D R, et al.. 340-GHz 3D radar imaging test bed with 10-Hz frame rate[C]. SPIE Passive and Active Millimeter-Wave Imaging XV, Baltimore, Maryland, 2012, 8362: 836206.
[14]	Wanke M C, Mangan M A, and Foltynowicz R J. Atmospheric propagation of THz radiation[R]. Albuquerque, NM, USA: Sandia National Laboratories, 2005: 2005–6389.
[15]	戴宁, 葛进, 胡淑红, 等. 太赫兹探测技术在遥感应用中的研究进展[J]. 中国电子科学研究院学报, 2009, 4(3): 232–237. DOI: 10.3969/j.issn.1673-5692.2009.03.002. Dai Ning, Ge Jin, Hu Shu-hong, et al. The development of the terahertz detection techniques in the applications of remote sensing[J]. Journal of CAEIT, 2009, 4(3): 232–237. DOI: 10.3969/j.issn.1673-5692.2009.03.002.
[16]	刘其军, 刘正堂, 冯丽萍, 等. 闪锌矿型CdTe电子结构和光学性质的第一性原理[J]. 中国科学院研究生院学报, 2009, 26(5): 615–620. Liu Qi-jun, Liu Zheng-tang, Feng Li-ping, et al. First-principle calculations of electronic structure and optical properties of Zinc blende CdTe[J]. Journal of the Graduate School of the Chinese Academy of Sciences, 2009, 26(5): 615–620.
[17]	孙立忠, 陈效双, 郭旭光, 等. CdTe和HgTe能带结构的第一性原理计算[J]. 红外与毫米波学报, 2004, 23(4): 271–275. DOI: 10.3321/j.issn:1001-9014.2004.04.007. Sun Li-zhong, Chen Xiao-shuang, Guo Xu-guang, et al. First principles calculation of the band structure of CdTe and HgTe[J]. Journal of Infrared and Millimeter Waves, 2004, 23(4): 271–275. DOI: 10.3321/j.issn:1001-9014.2004.04.007.
[18]	王琰, 侯延冰, 唐爱伟, 等. 不同稳定剂对水溶性CdTe纳米晶光学性质的影响[J]. 发光学报, 2008, 29(1): 171–175. Wang Yan, Hou Yan-bing, Tang Ai-wei, et al. Influence of different stabilizers on optical properties of water-soluble CdTe nanocrystals[J]. Chinese Journal of Luminescence, 2008, 29(1): 171–175.
[19]	Merad A E, Kanoun M B, Merad G, et al. Full-potential investigation of the electronic and optical properties of stressed CdTe and ZnTe[J]. Materials Chemistry and Physics, 2005, 92(2/3): 333–339. DOI: 10.1016/j.matchemphys.2004.10.031.
[20]	吴文智, 郑植仁, 金钦汉, 等. 水溶性CdTe量子点的三阶光学非线性极化特性[J]. 物理学报, 2008, 57(2): 1177–1182. DOI: 10.7498/aps.57.1177. Wu Wen-zhi, Zheng Zhi-ren, Jin Qin-han, et al. The property of third-order optical nonlinear susceptibility of water soluble CdTe quantum dots[J]. Acta Physica Sinica, 2008, 57(2): 1177–1182. DOI: 10.7498/aps.57.1177.
[21]	Vodopyanov L K. Optical studies of II–VI alloy lattice dynamics[J]. Journal of Alloys and Compounds, 2004, 371(1/2): 72–76. DOI: 10.1016/j.jallcom.2003.05.007.
[22]	Deligoz E, Colakoglu K, and Ciftci Y. Elastic, electronic, and lattice dynamical properties of CdS, CdSe, and CdTe[J]. Physica B:Condensed Matter, 2006, 373(1): 124–130. DOI: 10.1016/j.physb.2005.11.099.
[23]	Polit J, Sheregii E M, Cebulski J, et al. Phonon and vibrational spectra of hydrogenated CdTe[J]. Journal of Applied Physics, 2006, 100(1): 013521. DOI: 10.1063/1.2211368.
[24]	Robouch B V, Zajdel P, Kisiel A, et al. Analysis of the phonon line profile of hydrogenated CdTe[J]. Journal of Physics:Condensed Matter, 2008, 20(32): 325217. DOI: 10.1088/0953-8984/20/32/325217.
[25]	Stergiou V C, Kontos A G, and Raptis Y S. Anharmonic effects and Faust-Henry coefficient of CdTe in the vicinity of the energy gap[J]. Physical Review B, 2008, 77(23): 235201. DOI: 10.1103/PhysRevB.77.235201.
[26]	Tan J J, Ji G F, Chen X R, et al. The high-pressure phase transitions and vibrational properties of zinc-blende XTe (X=Zn, Cd, Hg): Performance of local-density-approximation density functional theory[J]. Computational Materials Science, 2010, 48(4): 796–801. DOI: 10.1016/j.commatsci.2010.03.037.
[27]	Ceperley D M and Alder B J. Ground state of the electron gas by a stochastic method[J].Physical Review Letters, 1980, 45(7): 566–569. DOI: 10.1103/PhysRevLett.45.566.
[28]	Perdew J P and Zunger A. Self-interaction correction to density-functional approximations for many-electron systems[J]. Physical Review B, 1981, 23(10): 5048–5079. DOI: 10.1103/PhysRevB.23.5048.
[29]	Perdew J P and Wang Y. Accurate and simple analytic representation of the electron-gas correlation energy[J]. Physical Review B, 1992, 45(23): 13244–13249. DOI: 10.1103/PhysRevB.45.13244.
[30]	Perdew J P, Burke K, and Ernzerhof M. Generalized gradient approximation made simple[J]. Physical Review Letters, 1996, 77(18): 3865–3868. DOI: 10.1103/PhysRevLett.77.3865.
[31]	黄昆, 韩汝琦. 固体物理[M]. 北京: 高等教育出版社, 1988: 104–111. Huang Kun and Han Ru-qi. Solid State Physics[M]. Beijing: Higher Education Press, 1988: 104–111.
[32]	Monkhorst H J and Pack J D. Special points for Brillonin-zone integrations[J]. Physical Revie B, 1976, 13(12): 5188–5192. DOI: 10.1103/PhysRevB.13.5188.
[33]	Rabadanov M K, Verin I A, Ivanov Y M, et al. Refinement of the atomic structure of CdTe single crystals[J]. Crystallography Reports, 2001, 46(4): 636–641. DOI: 10.1134/1.1387130.

施引文献

期刊类型引用(5)

1.	马助兴，张立硕，徐红元，郑焕坤，赵智龙，康哲. 面向变电站智能巡检5G大规模MIMO融合定位方法. 全球能源互联网. 2023(03): 308-315 . 百度学术
2.	任笑莹，王英民，张立琛，王奇，廉杰. 均匀双圆环阵虚拟成阵方法研究. 水下无人系统学报. 2023(05): 735-745 . 百度学术
3.	苏晓龙，户盼鹤，刘天鹏，彭勃，程耘，刘振. 基于深度展开ISTA网络的混合源定位方法. 信号处理. 2022(10): 2082-2091 . 百度学术
4.	母采凤，李森，吕梦然. 脉冲噪声环境下基于矩阵差分的远近场混合源定位. 信号处理. 2022(11): 2342-2349 . 百度学术
5.	张国鑫，易伟，孔令讲. 基于1比特量化的大规模MIMO雷达系统直接定位算法. 雷达学报. 2021(06): 970-981 . 本站查看