CN105353338B - 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法 - Google Patents

宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法 Download PDF

Info

Publication number
CN105353338B
CN105353338B CN201510628069.3A CN201510628069A CN105353338B CN 105353338 B CN105353338 B CN 105353338B CN 201510628069 A CN201510628069 A CN 201510628069A CN 105353338 B CN105353338 B CN 105353338B
Authority
CN
China
Prior art keywords
mrow
msub
msup
mover
mtr
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510628069.3A
Other languages
English (en)
Other versions
CN105353338A (zh
Inventor
甄佳奇
王志芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Heilongjiang University
Original Assignee
Heilongjiang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Heilongjiang University filed Critical Heilongjiang University
Priority to CN201510628069.3A priority Critical patent/CN105353338B/zh
Publication of CN105353338A publication Critical patent/CN105353338A/zh
Application granted granted Critical
Publication of CN105353338B publication Critical patent/CN105353338B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/04—Details
    • G01S3/10—Means for reducing or compensating for quadrantal, site, or like errors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，涉及宽带信号超分辨测向中存在的阵列误差的校正方法。本发明为了解决现有的迭代最小化代价函数对阵列通道幅相不一致性误差进行校正的方法不适用于宽带信号的问题。本发明利用各个频点上的信号构建对应的优化函数，之后利用信号的空域稀疏性，通过稀疏贝叶斯学习方法分别对各个频点上的函数进行迭代优化处理，最后对所有频点上的信息进行融合估计出信号到达方向。该方法可以有效的实现阵列通道幅相不一致性误差存在时的阵列误差校正，并且利用多片数字信号处理器有效提高算法的运行速度。本发明适用于宽带信号超分辨测向中存在的阵列误差的校正领域。

Description

宽带信号超分辨测向中的阵列通道幅相不一致性误差校正 方法
技术领域
本发明涉及宽带信号超分辨测向中存在的阵列误差的校正方法。
背景技术
超分辨测向是阵列信号处理中的一个重要研究内容，在无线电监测、物联网和电子对抗等领域有着较广泛的应用。目前多数的测向方法都是以精确的掌握阵列流型为前提。而实际的测向***当中，各阵列通道的增益和长短往往不一致，导致测向估计时经常伴随着阵列通道幅相不一致性误差，这直接导致了很多的超分辨测向方法性能的恶化，甚至失效，所以有必要对其进行校正处理。
参数类的校正方法通常可以分为有源校正和自校正。有源校正可通过在空间设置方位已知的辅助信源对阵列扰动参数进行离线估计，而自校正方法通常根据某种优化函数对空间信源的方位与阵列扰动参数联合估计。较早的自校正算法只针对阵元的位置误差或阵元间幅相不一致性误差，这两种误差其实可以用相同的数学模型表示(阵元的位置误差可以看成是阵元间的相位不一致)，它们都是与方位不相关的误差。对于这类误差，A.Paulraj和T.Kallath提出了利用阵列输出协方差矩阵的特殊结构，得到幅相误差之间相互关系的线性方程组，从而可实现对均匀线阵幅相误差和信号源的到达方向估计。Benjamin Friedlander和Anthony J.Weiss利用阵列输出协方差矩阵特征分解后噪声子空间和信号子空间正交的特点，并结合多重信号分类算法，提出了一种迭代最小化代价函数对阵列幅相误差和到达方向同时估计的算法。然而它们只适用于窄带信号，对于宽带信号超分辨测向中阵列通道的幅相不一致性误差校正技术，公开发表的文献并不多见。
发明内容
本发明为了解决现有的迭代最小化代价函数对阵列通道幅相不一致性误差进行校正的方法不适用于宽带信号的问题。
宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，包括下述步骤：
步骤1：建立含有阵列通道幅相不一致性误差的阵列信号模型：
当阵列当中存在阵列通道幅相不一致性误差时，阵列输出可以表示为
X'(fi)＝A'(fi,α)S(fi)+N(fi)，i＝1,2,…,J (12)
其中，S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵；N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵，均值为0，方差为μ2(fi)；
A'(fi,α)＝[a'(fi,α1),…,a'(fi,αk),…,a'(fi,αK)]，i＝1,2,…,J (13)
为存在阵列通道幅相不一致性误差时频点fi上的阵列流型矩阵，a'(fi,αk)为存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量；
有
R'(fi)＝E{X'(fi)(X'(fi))H}，i＝1,2,…,J (14)
R'(fi)为存在阵列通道幅相不一致性误差时频点fi上的阵列接收信号协方差矩阵；
设W(fi)为阵列扰动矩阵，表示频点fi上阵列通道幅相不一致性误差；表示为
W(fi)＝[W1(fi),…,Wm(fi),…,WM(fi)]T (15)
其中
为信号频点fi上第m路通道的幅相不一致性误差，ρm(fi)、分别为频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差；
A(fi,α)＝[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵，a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量；
因此存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量可以表示为
步骤2：对含有阵列通道幅相不一致性误差的阵列信号参数进行估计：
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向，从而可得出频点fi上阵列流型矩阵的稀疏表示
其中，为频点fi上第l个稀疏信号的阵列导向矢量，相应的可获得存在阵列通道幅相不一致性误差时频点fi上阵列流型矩阵的稀疏表示
其中，为存在阵列通道幅相不一致性误差时频点fi上第l个稀疏信号的阵列导向矢量，则可得出存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的稀疏表示
其中，Λ(fi)为一个只与原信号有关的参数，与误差无关，为Λ(fi)的稀疏表示；为频点fi上阵列通道幅相不一致扰动矢量，
的协方差矩阵为
式(19)中为S(fi)的稀疏表示，
其中，为稀疏矩阵，为S(fi,kp)的稀疏表示，中只包含K个非零元素，为中的第l个元素，当且仅当时中的元素不全为零且有l＝1,2,…,L，k＝1,2,…,K；故此可以看成是S(fi)中加入了许多0元素后得到的矩阵；
设δ(fi)＝[δ1(fi),…,δl(fi),…,δL(fi)]T为中元素的方差，反映了信号的能量，即有
其中，Σ(fi)＝diag(δ(fi))，即服从均值为0，方差为δ(fi)的高斯分布；
由于可以看成是S(fi)中加入了许多0元素后得到的向量，所以δ(fi)包含了K个非零元素，并且有K<<L，根据δ(fi)，结合w(fi)和噪声方差μ2(fi)估计出从而重构出原信号，同时对误差进行校正；
根据式(19)可知，存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的概率密度为
其中，IM是M×M维的单位阵；
结合式(19)、(21)和(22)可得
采用期望最大化(Expectation Maximization,EM)方法来对w(fi)、μ2(fi)和δl(fi)进行迭代估计，得出估计值和对应的可得到以及
步骤3：利用和对阵列误差进行校正并对信号到达方向求解；
令X为一段观测时间内阵列接收到的所有频点信号的和构成的向量，由于各频点的信号具有统计独立性，因此各频点接收信号的联合概率密度为
对式(36)两端取对数有
因此令式(37)最大化即可求得信号到达方向，即信号到达方向的估计值k＝1,2,…,K，即可以通过
求得；
经过推导有
其中，Re{·}为求{·}的实部；Ω-k、分别表示从Ω和中去掉其中的第k个元素；k＝1,2,…,K；
根据的表达式可求得再根据式(16)和(15)求得W(fi)，利用W(fi)进行阵列校正可求得a'(fi,αk)和A'(fi,Ω-k)，再根据以上参数和公式(39)，能够得到经过阵列校正后的信号到达方向的估计值
本发明具有以下有益效果：
本发明提出了一种存在阵列通道幅相不一致性误差时的宽带信号超分辨测向误差校正方法，利用各个频点上的信号构建对应的优化函数，之后利用信号的空域稀疏性，通过稀疏贝叶斯学习方法分别对各个频点上的函数进行迭代优化处理，最后对所有频点上的信息进行融合估计出宽带信号到达方向。本发明可以有效的实现阵列通道幅相不一致性误差存在时的阵列误差校正，当信噪比为10dB，每个频点采样快拍数为40时，精度可达0.8°/σ。
而且本发明的方法可以用多片数字信号处理器进行处理，可以有效的提高算法的运行速度。
附图说明
图1为宽带信号超分辨测向阵列信号模型示意图；
图2为宽带信号探测***装置图；
图3为具体实施方式五的宽带信号超分辨测向装置图；
图4为具体实施方式六的宽带信号超分辨测向装置图；
图5为具体实施方式七的宽带信号超分辨测向装置图。
具体实施方式
具体实施方式一：
宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，包括下述步骤：
步骤1：建立含有阵列通道幅相不一致性误差的阵列信号模型：
当阵列当中存在阵列通道幅相不一致性误差时，阵列输出可以表示为
X'(fi)＝A'(fi,α)S(fi)+N(fi)，i＝1,2,…,J (12)
其中，S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵；N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵，均值为0，方差为μ2(fi)；
A'(fi,α)＝[a'(fi,α1),…,a'(fi,αk),…,a'(fi,αK)]，i＝1,2,…,J (13)
为存在阵列通道幅相不一致性误差时频点fi上的阵列流型矩阵，a'(fi,αk)为存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量；
有
R'(fi)＝E{X'(fi)(X'(fi))H}，i＝1,2,…,J (14)
R'(fi)为存在阵列通道幅相不一致性误差时频点fi上的阵列接收信号协方差矩阵；
设W(fi)为阵列扰动矩阵，表示频点fi上阵列通道幅相不一致性误差；表示为
W(fi)＝[W1(fi),…,Wm(fi),…,WM(fi)]T (15)
其中
为信号频点fi上第m路通道的幅相不一致性误差，ρm(fi)、分别为频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差；
A(fi,α)＝[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵，a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量；
因此存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量可以表示为
步骤2：对含有阵列通道幅相不一致性误差的阵列信号参数进行估计：
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向，从而可得出频点fi上阵列流型矩阵的稀疏表示
其中，为频点fi上第l个稀疏信号的阵列导向矢量，相应的可获得存在阵列通道幅相不一致性误差时频点fi上阵列流型矩阵的稀疏表示
其中，为存在阵列通道幅相不一致性误差时频点fi上第l个稀疏信号的阵列导向矢量，则可得出存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的稀疏表示
其中，Λ(fi)为一个只与原信号有关的参数，与误差无关，为Λ(fi)的稀疏表示；为频点fi上阵列通道幅相不一致扰动矢量，
的协方差矩阵为
式(19)中为S(fi)的稀疏表示，
其中，为稀疏矩阵，为S(fi,kp)的稀疏表示，中只包含K个非零元素，为中的第l个元素，当且仅当时中的元素不全为零且有l＝1,2,…,L，k＝1,2,…,K；故此可以看成是S(fi)中加入了许多0元素后得到的矩阵；
设δ(fi)＝[δ1(fi),…,δl(fi),…,δL(fi)]T为中元素的方差，反映了信号的能量，即有
其中，Σ(fi)＝diag(δ(fi))，即服从均值为0，方差为δ(fi)的高斯分布；
由于可以看成是S(fi)中加入了许多0元素后得到的向量，所以δ(fi)包含了K个非零元素，并且有K<<L，根据δ(fi)，结合w(fi)和噪声方差μ2(fi)估计出从而重构出原信号，同时对误差进行校正；
根据式(19)可知，存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的概率密度为
其中，IM是M×M维的单位阵；
结合式(19)、(21)和(22)可得
采用期望最大化(Expectation Maximization,EM)方法来对w(fi)、μ2(fi)和δl(fi)进行迭代估计，得出估计值和对应的可得到以及
步骤3：利用和对阵列误差进行校正并对信号到达方向求解；
令X为一段观测时间内阵列接收到的所有频点信号的和构成的向量，由于各频点的信号具有统计独立性，因此各频点接收信号的联合概率密度为
对式(36)两端取对数有
因此令式(37)最大化即可求得信号到达方向，即信号到达方向的估计值k＝1,2,…,K，即可以通过
求得；
经过推导有
其中，Re{·}为求{·}的实部；Ω-k、分别表示从Ω和中去掉其中的第k个元素；k＝1,2,…,K；
根据的表达式可求得再根据式(16)和(15)求得W(fi)，利用W(fi)进行阵列校正可求得a'(fi,αk)和A'(fi,Ω-k)，再根据以上参数和公式(39)，能够得到经过阵列校正后的信号到达方向的估计值
具体实施方式二：
本实施方式步骤1所述建立含有阵列通道幅相不一致性误差的阵列信号模型的具体步骤如下：
步骤1.1：建立理想阵列信号模型：
如图1所示，设有K个远场宽带信号sk(t)，k＝1,2,…,K，入射到M个全向阵元组成的宽带均匀直线阵列上，到达方向为α＝[α1,…,αk,…,αK]，阵元间距为d；远场宽带信号sk(t)，简称宽带信号sk(t)；
将第1个阵元作为相位参考点，在理想情况下，第m个阵元的输出表示为
其中，表示第k个宽带信号sk(t)到达第m个阵元相对于它到达相位参考点的延时，c为电磁波在真空中的传播速度，nm(t)为第m个阵元接收到的高斯白噪声；
假设宽带信号的频率范围为[fLow,fHigh]，利用离散傅里叶变换将宽带信号分成J个频点，经过窄带滤波器组将它们分开，则第i组滤波器阵列输出信号表示为
X(fi)＝A(fi,α)S(fi)+N(fi)，i＝1,2,…,J (2)
其中，fLow≤fi≤fHigh，i＝1,2,…,J；
假设在每个频点上进行了KP次采样，X(fi)的矩阵形式表示为
X(fi)＝[X(fi,1),…,X(fi,kp),…,X(fi,KP)]，i＝1,2,…,J (3)
其中，X(fi,kp)为X(fi)的第kp次数据采样矩阵，
X(fi,kp)＝[X1(fi,kp),…,Xm(fi,kp),…,XM(fi,kp)]T，i＝1,2,…,J， (4)
Xm(fi,kp)为第m个阵元在频点fi上得到的第kp次数据采样值；
A(fi,α)为理想情况下频点fi上的阵列流型矩阵，
A(fi,α)＝[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]，i＝1,2,…,J， (5)
a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量，
其中，φk是第k个信号的相位；j是复数标志；
S(fi)＝[S(fi,1),…,S(fi,kp),…,S(fi,KP)]，i＝1,2,…,J， (8)
为信号sk(t)经过傅立叶变换后的信号矢量矩阵，k＝1,2,…,K；
其中，S(fi,kp)为S(fi)的第kp次信号采样矩阵，
S(fi,kp)＝[S1(fi,kp),…Sk(fi,kp),…,SK(fi,kp)]T i＝1,2,…,J (9)
Sk(fi,kp)为第k个信号在频点fi上得到的第kp次信号采样值；
N(fi)＝[N(fi,1),…,N(fi,kp),…,N(fi,KP)]i＝1,2,…,J (10)
为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵，均值为0，方差为μ2(fi)；m＝1,2,…,M；
N(fi,kp)为N(fi)的第kp次噪声采样矩阵，
N(fi,kp)＝[N1(fi,kp),…,Nm(fi,kp),…,NM(fi,kp)]T i＝1,2,…,J (11)
Nm(fi,kp)为第m个阵元在频点fi上得到的第kp次噪声采样值；
步骤1.2：在理想阵列信号模型基础上建立含有阵列通道幅相不一致性误差的阵列信号模型：
当阵列当中存在阵列通道幅相不一致性误差时，频点fi上的阵列输出可以表示为
X'(fi)＝A'(fi,α)S(fi)+N(fi)，i＝1,2,…,J (12)
其中
A'(fi,α)＝[a'(fi,α1),…,a'(fi,αk),…,a'(fi,αK)]，i＝1,2,…,J (13)
为存在阵列通道幅相不一致性误差时频点fi上的阵列流型矩阵，a'(fi,αk)为对应的阵列导向矢量；则有
R'(fi)＝E{X'(fi)(X'(fi))H}，i＝1,2,…,J (14)
R'(fi)为存在阵列通道幅相不一致性误差时频点fi上的接收信号协方差矩阵；
设W(fi)为阵列扰动矩阵，表示频点fi上阵列通道幅相不一致性误差；表示为
W(fi)＝[W1(fi),…,Wm(fi),…,WM(fi)]T (15)
其中
为频点fi上第m路通道的幅相不一致性误差，ρm(fi)、分别为频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差；
因此存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量可以表示为
其它步骤和参数与具体实施方式一相同。
具体实施方式三：
本实施方式步骤2中所述的采用期望最大化方法来对w(fi)、μ2(fi)和δl(fi)进行迭代估计的具体步骤如下：
在期望最大化方法中的E-step步中，首先对的分布函数进行计算
其中运算符<·>表示求解条件期望；
在期望最大化方法中的M-step步中，分别求取分布函数对各未知参数的导数，即对取极值来对各未知参数求解；
分别令以上的导数为0，即可求得第p次迭代时各个未知参数的估计值
其中(p)代表迭代次数，式(28)中
为矩阵第r1行、r2列的元素，其中tr[·]表示求迹运算；
式(31)中
O(fi)＝Σ(fi)(A'(fi,Ω))H(μ2(fi)IM+A'(fi,Ω)Σ(fi)(A'(fi,Ω))H)-1X(fi)(32)
为中间变量；
Ξ(fi)
(33)
＝Σ(fi)-Σ(fi)(A'(fi,Ω))H(μ2(fi)IM+A'(fi,Ω)Σ(fi)(A'(fi,Ω))H)-1A'(fi,Ω)Σ(fi)为中间变量；
式(28)中
式(29)中
式(34)中，Ψr(fi)为中间变量，为M×M维的矩阵，只有在第±r对角线上的元素全为1，其余元素全为0，由于直接利用式(28)～(29)计算w(fi)和μ2(fi)比较复杂，因此可将式(31)～(35)代入式(28)～(29)中对等式进行化简并对w(fi)和μ2(fi)求解；
当迭代若干步后，w(fi)、μ2(fi)和δl(fi)三个量估计值的变化趋于0，此时可认为它们已经收敛，则可得出最后的估计值和对应得到以及
其它步骤和参数与具体实施方式二相同。
具体实施方式四：
本实施方式所述的ρ1(fi)＝1，W1(fi)＝1。
其它步骤和参数与具体实施方式三相同。
具体实施方式五：参照图2和图3具体说明本实施方式，
本实施方式为实现具体实施方式一至四所述方法的宽带信号探测***及实现探测的方法，
如图2所示，宽带信号探测***包括：宽带均匀直线阵列1、多通道宽带数字接收机2和宽带信号超分辨测向装置3；
如图3所示，宽带信号超分辨测向装置3包括6片数字信号处理器，即DSP，采用快速串行输入输出口，即SRIO口，组成多处理器***实现并行处理。其中，DSP3-1为主DSP，DSP3-2～DSP3-6为从DSP；宽带信号超分辨测向装置3还包括CPLD3-7、PROM3-8、FLASH3-9、SRAM3-10、JTAG3-11、电源、晶振和复位。
数字信号处理器采用Texas Instruments(TI)公司的TMS320C6678，采用6片处理器并行处理，6片DSP通过SRIO口连接，上电后PROM3-8首先将程序加载给CPLD3-7，FLASH3-9也将程序加载给这6块DSP(3-1～3-6)，之后主DSP3-1开始接收多通道宽带数字接收机2传来的J个频点的观测数据，把它们分为W组，假设J＝30，W＝6，则每片DSP可以处理U＝30/6＝5个频点的观测数据，主DSP3-1通过SRIO口将其它从DSP(3-2～3-6)负责处理的观测数据传递给它们，之后每个DSP(3-1～3-6)都按照以上理论推导的步骤进行求解，之后5片从DSP(3-2～3-6)将各自的误差估计值通过SRIO口传给主DSP3-1，主DSP3-1再利用这些结果，结合式(39)得出信号到达角度。其中SRAM3-10负责存储数据，JTAG3-11负责对DSP(3-1～3-6)进行调试，电源负责整体供电，晶振负责提供时钟，复位负责提供复位信号。
具体实施方式六：参照图2和图4具体说明本实施方式，
本实施方式为实现具体实施方式一至四所述方法的宽带信号探测***及实现探测的方法，
如图2所示，宽带信号探测***包括：宽带均匀直线阵列1、多通道宽带数字接收机2和宽带信号超分辨测向装置3；
如图4所示，宽带信号超分辨测向装置3包括6片数字信号处理器，即DSP，采用共享总线紧耦合方式组成多处理器***实现并行处理。其中，DSP3-1为主DSP，DSP3-2～DSP3-6为从DSP；宽带信号超分辨测向装置3还包括CPLD3-7、PROM3-8、FLASH3-9、SRAM3-10、JTAG3-11、电源、晶振和复位。
数字信号处理器采用Analog Device Instruments(ADI)公司的ADSP-TS201S，采用6片DSP并行处理，6片DSP通过共享总线紧耦合方式连接，上电后PROM3-8首先将程序加载给CPLD3-7对DSP(3-1～3-6)进行配置，之后FLASH3-9将程序加载给这6块DSP(3-1～3-6)，主DSP3-1开始接收多通道宽带数字接收机2传来的J个频点的观测数据，把它们分为W组，假设J＝30，W＝6，则每片DSP可以处理U＝30/6＝5个频点的观测数据，主DSP3-1通过总线将其它从DSP(3-2～3-6)负责处理的观测数据传递给它们，之后每个DSP(3-1～3-6)都按照以上理论推导的步骤进行求解，之后5片从DSP(3-2～3-6)将各自的误差估计值通过总线传给主DSP3-1，主DSP3-1再利用这些结果，结合式(39)得出信号到达角度。其中SRAM3-10负责存储数据，JTAG3-11负责对DSP(3-1～3-6)进行调试，电源负责整体供电，晶振负责提供时钟，复位负责提供复位信号。
具体实施方式七：参照图2和图5具体说明本实施方式，
本实施方式为实现具体实施方式一至四所述方法的宽带信号探测***及实现探测的方法，
如图2所示，宽带信号探测***包括：宽带均匀直线阵列1、多通道宽带数字接收机2和宽带信号超分辨测向装置3；
如图5所示，宽带信号超分辨测向装置3包括6片数字信号处理器，即DSP，采用链路口级联松耦合方式组成多处理器***实现并行处理。其中，DSP3-1为主DSP，DSP3-2～DSP3-6为从DSP；宽带信号超分辨测向装置3还包括CPLD3-7、PROM3-8、FLASH3-9、SRAM3-10、JTAG3-11、电源、晶振和复位。
数字信号处理器采用Analog Device Instruments(ADI)公司的ADSP-TS201S，采用6片处理器并行处理，6片DSP通过链路口级联松耦合方式连接，上电后PROM3-8首先将程序加载给CPLD3-7，FLASH3-9将这6片DSP的程序加载给主DSP3-1，主DSP3-1再依次将其它从DSP(3-2～3-6)的程序通过链路口一级一级传给它们，之后主DSP3-1开始接收多通道宽带数字接收机2传来的J个频点的观测数据，把它们分为W组，假设J＝30，W＝6，则每片DSP可以处理U＝30/6＝5个频点的观测数据，主DSP3-1再通过链路口将其它DSP(3-2～3-6)负责处理的观测数据一级一级逐次传递给它们，之后每个DSP(3-1～3-6)都按照以上理论推导的步骤进行求解，之后5片从DSP(3-2～3-6)将各自的误差估计值通过链路口一级一级逐次上传到主DSP3-1，主DSP3-1再利用这些结果，结合式(39)得出信号到达角度。其中SRAM3-10负责存储数据，JTAG3-11负责对DSP(3-1～3-6)进行调试，电源负责整体供电，晶振负责提供时钟，复位负责提供复位信号。

Claims (4)

1.宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，其特征在于包括下述步骤：
步骤1：建立含有阵列通道幅相不一致性误差的阵列信号模型：
当阵列当中存在阵列通道幅相不一致性误差时，阵列输出表示为
X'(fi)＝A'(fi,α)S(fi)+N(fi)，i＝1,2,…,J (12)
其中，S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵；N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵，均值为0，方差为μ2(fi)；J为频点数量；α＝[α1,…,αk,…,αK]为K个远场宽带信号入射到M个全向阵元组成的宽带均匀直线阵列上的到达方向；
A'(fi,α)＝[a'(fi,α1),…,a'(fi,αk),…,a'(fi,αK)]，i＝1,2,…,J (13)
为存在阵列通道幅相不一致性误差时频点fi上的阵列流型矩阵，a'(fi,αk)为存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量；K为远场宽带信号数量；
有
R'(fi)＝E{X'(fi)(X'(fi))H}，i＝1,2,…,J (14)
R'(fi)为存在阵列通道幅相不一致性误差时频点fi上的阵列接收信号协方差矩阵；
设W(fi)为阵列扰动矩阵，表示频点fi上阵列通道幅相不一致性误差；表示为
W(fi)＝[W1(fi),…,Wm(fi),…,WM(fi)]T (15)
其中
为信号频点fi上第m路通道的幅相不一致性误差，ρm(fi)、分别为频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差；
A(fi,α)＝[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵，a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量；
存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量表示为
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>,</mo> <msub> <mi>W</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>j&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>W</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>W</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>W</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
步骤2：对含有阵列通道幅相不一致性误差的阵列信号参数进行估计：
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向，从而得出频点fi上阵列流型矩阵的稀疏表示
<mrow> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&amp;lsqb;</mo> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>L</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow>
其中，为频点fi上第l个稀疏信号的阵列导向矢量，相应的获得存在阵列通道幅相不一致性误差时频点fi上阵列流型矩阵的稀疏表示
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&amp;lsqb;</mo> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>L</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>W</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
其中，d为阵元间距；c为电磁波在真空中的传播速度；为存在阵列通道幅相不一致性误差时频点fi上第l个稀疏信号的阵列导向矢量，则得出存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的稀疏表示
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>w</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>J</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
其中，Λ(fi)为一个只与原信号有关的参数，为Λ(fi)的稀疏表示；为频点fi上阵列通道幅相不一致扰动矢量，
的协方差矩阵为
<mrow> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>E</mi> <mo>{</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>}</mo> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>J</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
式(19)中为S(fi)的稀疏表示，KP为采样次数，kp表示第kp次采样；
其中，为稀疏矩阵，为S(fi,kp)的稀疏表示，中只包含K个非零元素，为中的第l个元素，当且仅当时中的元素不全为零且有故此看成是S(fi)中加入了许多0元素后得到的矩阵；S(fi,kp)为S(fi)的第kp次信号采样矩阵，Sk(fi,kp)为第k个信号在频点fi上得到的第kp次信号采样值；
设δ(fi)＝[δ1(fi),…,δl(fi),…,δL(fi)]T为中元素的方差，反映了信号的能量，即有
<mrow> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>~</mo> <mi>N</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mi>&amp;Sigma;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
其中，Σ(fi)＝diag(δ(fi))，即服从均值为0，方差为δ(fi)的高斯分布；
根据式(19)，存在阵列通道幅相不一致性误差时频点fi上的阵列输出信号的概率密度为
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>|</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mi>w</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>|</mo> <mi>&amp;pi;</mi> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mi>K</mi> <mi>P</mi> </mrow> </msup> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>|</mo> <mo>|</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>w</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msubsup> <mo>|</mo> <mn>2</mn> <mn>2</mn> </msubsup> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
其中，IM是M×M维的单位阵；M为全向阵元的数量；
结合式(19)、(21)和(22)得
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mi>&amp;delta;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mi>w</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>|</mo> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mi>&amp;Sigma;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mi>K</mi> <mi>P</mi> </mrow> </msup> <mo>&amp;times;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> <mi>&amp;Sigma;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
采用期望最大化方法来对w(fi)、μ2(fi)和δl(fi)进行迭代估计，得出估计值 和对应的得到以及
步骤3：利用和对阵列误差进行校正并对信号到达方向求解；
令X为一段观测时间内阵列接收到的所有频点信号的和构成的向量，由于各频点的信号具有统计独立性，因此各频点接收信号的联合概率密度为
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mover> <mi>&amp;delta;</mi> <mo>^</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mover> <mi>w</mi> <mo>^</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>|</mo> <mi>&amp;pi;</mi> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mi>J</mi> <mo>&amp;times;</mo> <mi>K</mi> <mi>P</mi> </mrow> </msup> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mo>|</mo> <mrow> <mo>(</mo> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mi>K</mi> <mi>P</mi> </mrow> </msup> <mo>&amp;times;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>36</mn> <mo>)</mo> </mrow> </mrow>
对式(36)两端取对数有
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>I</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>P</mi> <mo>(</mo> <mi>X</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mi>J</mi> <mo>&amp;times;</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mi>I</mi> <mi>n</mi> <mi>&amp;pi;</mi> <mo>-</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mi>I</mi> <mi>n</mi> <mo>|</mo> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>37</mn> <mo>)</mo> </mrow> </mrow>
令式(37)最大化即求得信号到达方向，即信号到达方向的估计值即通过
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>I</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>P</mi> <mo>(</mo> <mi>X</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mrow> <mo>&amp;part;</mo> <mi>&amp;alpha;</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>38</mn> <mo>)</mo> </mrow> </mrow>
求得；
经过推导有
<mrow> <msub> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mo>=</mo> <mi>arg</mi> <munder> <mi>max</mi> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </munder> <mo>|</mo> <mi>Re</mi> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <msub> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>)</mo> </mrow> <mo>&amp;times;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mrow> <mo>(</mo> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mrow> <mo>(</mo> <msup> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> </mrow> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;times;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>J</mi> </munderover> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <mi>A</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mover> <mi>&amp;Sigma;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;Omega;</mi> <mrow> <mo>-</mo> <mi>k</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;times;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>39</mn> <mo>)</mo> </mrow> </mrow>
其中，Re{·}为求{·}的实部；Ω-k、分别表示从Ω和中去掉其中的第k个元素；k＝1,2,…,K；
根据的表达式求得再根据式(16)和(15)求得W(fi)，利用W(fi)进行阵列校正求得a'(fi,αk)和A'(fi,Ω-k)，再根据以上参数和公式(39)，能够得到经过阵列校正后的信号到达方向的估计值
2.根据权利要求1所述的宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，其特征在于步骤1所述建立含有阵列通道幅相不一致性误差的阵列信号模型的具体步骤如下：
步骤1.1：建立理想阵列信号模型：
设有K个远场宽带信号sk(t)，k＝1,2,…,K，入射到M个全向阵元组成的宽带均匀直线阵列上，到达方向为α＝[α1,…,αk,…,αK]，阵元间距为d；远场宽带信号sk(t)，简称宽带信号sk(t)；
将第1个阵元作为相位参考点，在理想情况下，第m个阵元的输出表示为
<mrow> <msub> <mi>x</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msub> <mi>s</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>&amp;tau;</mi> <mi>m</mi> </msub> <mo>(</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>n</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中，表示第k个宽带信号sk(t)到达第m个阵元相对于它到达相位参考点的延时，c为电磁波在真空中的传播速度，nm(t)为第m个阵元接收到的高斯白噪声；
假设宽带信号的频率范围为[fLow,fHigh]，利用离散傅里叶变换将宽带信号分成J个频点，经过窄带滤波器组将它们分开，则第i组滤波器阵列输出信号表示为
X(fi)＝A(fi,α)S(fi)+N(fi)，i＝1,2,…,J (2)
其中，fLow≤fi≤fHigh，i＝1,2,…,J；
假设在每个频点上进行了KP次采样，X(fi)的矩阵形式表示为
X(fi)＝[X(fi,1),…,X(fi,kp),…,X(fi,KP)]，i＝1,2,…,J (3)
其中，X(fi,kp)为X(fi)的第kp次数据采样矩阵，
X(fi,kp)＝[X1(fi,kp),…,Xm(fi,kp),…,XM(fi,kp)]T，i＝1,2,…,J， (4)
Xm(fi,kp)为第m个阵元在频点fi上得到的第kp次数据采样值；
A(fi,α)为理想情况下频点fi上的阵列流型矩阵，
A(fi,α)＝[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]，i＝1,2,…,J， (5)
a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量，
<mrow> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>,</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>j&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>J</mi> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> <mo>=</mo> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>i</mi> </msub> <mfrac> <mi>d</mi> <mi>c</mi> </mfrac> <msub> <mi>sin&amp;alpha;</mi> <mi>k</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>J</mi> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中，φk是第k个信号的相位；j是复数标志；
S(fi)＝[S(fi,1),…,S(fi,kp),…,S(fi,KP)]，i＝1,2,…,J， (8)
为信号sk(t)经过傅立叶变换后的信号矢量矩阵，k＝1,2,…,K；
其中，S(fi,kp)为S(fi)的第kp次信号采样矩阵，
S(fi,kp)＝[S1(fi,kp),…Sk(fi,kp),…,SK(fi,kp)]T i＝1,2,…,J (9)
Sk(fi,kp)为第k个信号在频点fi上得到的第kp次信号采样值；
N(fi)＝[N(fi,1),…,N(fi,kp),…,N(fi,KP)] i＝1,2,…,J (10)
为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵，均值为0，方差为μ2(fi)；m＝1,2,…,M；
N(fi,kp)为N(fi)的第kp次噪声采样矩阵，
N(fi,kp)＝[N1(fi,kp),…,Nm(fi,kp),…,NM(fi,kp)]T i＝1,2,…,J (11)
Nm(fi,kp)为第m个阵元在频点fi上得到的第kp次噪声采样值；
步骤1.2：在理想阵列信号模型基础上建立含有阵列通道幅相不一致性误差的阵列信号模型：
当阵列当中存在阵列通道幅相不一致性误差时，频点fi上的阵列输出表示为
X'(fi)＝A'(fi,α)S(fi)+N(fi)，i＝1,2,…,J (12)
其中
A'(fi,α)＝[a'(fi,α1),…,a'(fi,αk),…,a'(fi,αK)]，i＝1,2,…,J (13)
为存在阵列通道幅相不一致性误差时频点fi上的阵列流型矩阵，a'(fi,αk)为对应的阵列导向矢量；则有
R'(fi)＝E{X'(fi)(X'(fi))H}，i＝1,2,…,J (14)
R'(fi)为存在阵列通道幅相不一致性误差时频点fi上的接收信号协方差矩阵；
设W(fi)为阵列扰动矩阵，表示频点fi上阵列通道幅相不一致性误差；表示为
W(fi)＝[W1(fi),…,Wm(fi),…,WM(fi)]T (15)
其中
为频点fi上第m路通道的幅相不一致性误差，ρm(fi)、分别为频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差；
存在阵列通道幅相不一致性误差时频点fi上第k个信号的阵列导向矢量表示为
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>a</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>,</mo> <msub> <mi>W</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>j&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>W</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>W</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> </mrow> </msup> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>W</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>a</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3.根据权利要求2所述的宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，其特征在于步骤2中所述的采用期望最大化方法来对w(fi)、μ2(fi)和δl(fi)进行迭代估计的具体步骤如下：
在期望最大化方法中的E-step步中，首先对的分布函数进行计算
其中运算符<·>表示求解条件期望；
在期望最大化方法中的M-step步中，分别求取分布函数对各未知参数的导数，即对取极值来对各未知参数求解；
<mrow> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>F</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mi>&amp;delta;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mi>w</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <mi>w</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mn>2</mn> <msup> <mi>&amp;mu;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&gt;</mo> <mi>w</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>-</mo> <mi>A</mi> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>&gt;</mo> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow> 5
<mrow> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>F</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mi>&amp;delta;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mi>w</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>M</mi> <mo>&amp;times;</mo> <mi>K</mi> <mi>P</mi> </mrow> <mrow> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>&lt;</mo> <mo>|</mo> <mo>|</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msubsup> <mo>|</mo> <mn>2</mn> <mn>2</mn> </msubsup> <mo>&gt;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>F</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>;</mo> <mi>&amp;delta;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <mi>w</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>,</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;delta;</mi> <mi>l</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>K</mi> <mi>P</mi> </mrow> <mrow> <msub> <mi>&amp;delta;</mi> <mi>l</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <msubsup> <mi>&amp;delta;</mi> <mi>l</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>&lt;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>K</mi> <mi>P</mi> </mrow> </munderover> <mo>|</mo> <msub> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mi>l</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>k</mi> <mi>p</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>&gt;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
分别令以上的导数为0，即求得第p次迭代时各个未知参数的估计值
<mrow> <msup> <mi>w</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mo>&gt;</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>-</mo> <mi>A</mi> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>&gt;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msup> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>M</mi> <mo>&amp;times;</mo> <mi>K</mi> <mi>P</mi> </mrow> </mfrac> <mo>&lt;</mo> <mo>|</mo> <mo>|</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msubsup> <mo>|</mo> <mn>2</mn> <mn>2</mn> </msubsup> <mo>&gt;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>&amp;delta;</mi> <mi>l</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>K</mi> <mi>P</mi> </mrow> </mfrac> <mo>&lt;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>K</mi> <mi>P</mi> </mrow> </munderover> <mo>|</mo> <msub> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mi>l</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>k</mi> <mi>p</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>&gt;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>
其中(p)代表迭代次数，
<mrow> <mtable> <mtr> <mtd> <mrow> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mo>&gt;</mo> <mrow> <mi>r</mi> <mn>1</mn> <mo>,</mo> <mi>r</mi> <mn>2</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>t</mi> <mi>r</mi> <mo>&amp;lsqb;</mo> <msubsup> <mi>&amp;Psi;</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>&amp;Psi;</mi> <mrow> <mi>r</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>O</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mi>O</mi> <mi>H</mi> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>+</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mi>&amp;Xi;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mi>A</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>
为矩阵第r1行、r2列的元素，其中tr[·]表示求迹运算；
式(31)中
O(fi)＝Σ(fi)(A'(fi,Ω))H(μ2(fi)IM+A'(fi,Ω)Σ(fi)(A'(fi,Ω))H)-1X(fi) (32)
为中间变量；
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>&amp;Xi;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>&amp;Sigma;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;Sigma;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <msup> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>I</mi> <mi>M</mi> </msub> <mo>+</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mi>&amp;Sigma;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mi>&amp;Sigma;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>33</mn> <mo>)</mo> </mrow> </mrow>
为中间变量；
式(28)中
<mrow> <mtable> <mtr> <mtd> <mrow> <mo>&lt;</mo> <msup> <mover> <mi>&amp;Lambda;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>-</mo> <mi>A</mi> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>&gt;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>t</mi> <mi>r</mi> <mo>&amp;lsqb;</mo> <msubsup> <mi>&amp;Psi;</mi> <mi>r</mi> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>X</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>O</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>A</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>t</mi> <mi>r</mi> <mo>&amp;lsqb;</mo> <msubsup> <mi>&amp;Psi;</mi> <mi>r</mi> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>O</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mi>O</mi> <mi>H</mi> </msup> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>+</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mi>&amp;Xi;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mi>A</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>
式(29)中
<mrow> <mtable> <mtr> <mtd> <mrow> <mo>&lt;</mo> <mo>|</mo> <mo>|</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mover> <mi>S</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msubsup> <mo>|</mo> <mn>2</mn> <mn>2</mn> </msubsup> <mo>&gt;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>|</mo> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mi>O</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msubsup> <mo>|</mo> <mi>F</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mi>K</mi> <mi>P</mi> <mo>&amp;times;</mo> <mi>t</mi> <mi>r</mi> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mi>&amp;Xi;</mi> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>A</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>&amp;Omega;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>
式(34)中，Ψr(fi)为中间变量，为M×M维的矩阵，只有在第±r对角线上的元素全为1，其余元素全为0，将式(31)～(35)代入式(28)～(29)中对等式进行化简并对w(fi)和μ2(fi)求解；
当迭代若干步后，w(fi)、μ2(fi)和δl(fi)三个量估计值的变化趋于0，则得出最后的估计值和对应得到以及
4.根据权利要求1、2或3所述的宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法，其特征在于所述的ρ1(fi)＝1，W1(fi)＝1。
CN201510628069.3A 2015-09-28 2015-09-28 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法 Expired - Fee Related CN105353338B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510628069.3A CN105353338B (zh) 2015-09-28 2015-09-28 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510628069.3A CN105353338B (zh) 2015-09-28 2015-09-28 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法

Publications (2)

Publication Number Publication Date
CN105353338A CN105353338A (zh) 2016-02-24
CN105353338B true CN105353338B (zh) 2017-10-10

Family

ID=55329339

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510628069.3A Expired - Fee Related CN105353338B (zh) 2015-09-28 2015-09-28 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法

Country Status (1)

Country Link
CN (1) CN105353338B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106371055B (zh) * 2016-10-26 2018-11-09 黑龙江大学 远近场宽带混合信号超分辨测向阵列幅相误差估计方法
CN106546948B (zh) * 2016-10-26 2019-01-11 黑龙江大学 基于幅相误差阵列的远近场宽带混合源中远场源测向方法
CN107037397A (zh) * 2017-06-21 2017-08-11 哈尔滨工业大学 一种在波达方向估计中校正多种阵列误差的方法
CN116699507A (zh) * 2023-06-07 2023-09-05 哈尔滨工程大学 基于原子范数的稀疏阵列完备模型误差自校正doa估计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0959522A1 (en) * 1998-05-19 1999-11-24 Toyota Jidosha Kabushiki Kaisha Method for determining phase correction values in radar apparatus
CN102890271A (zh) * 2012-10-25 2013-01-23 北京理工大学 一种外辐射源雷达阵列天线幅相一致性校正方法
CN203858360U (zh) * 2014-04-22 2014-10-01 南京大桥机器有限公司 一种雷达信号幅相一致性标校***
CN104931968A (zh) * 2015-06-18 2015-09-23 西安电子科技大学 基于FPGA的InSAR通道幅相误差估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5211787B2 (ja) * 2008-03-24 2013-06-12 日本電気株式会社 レーダ模擬信号発生器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0959522A1 (en) * 1998-05-19 1999-11-24 Toyota Jidosha Kabushiki Kaisha Method for determining phase correction values in radar apparatus
CN102890271A (zh) * 2012-10-25 2013-01-23 北京理工大学 一种外辐射源雷达阵列天线幅相一致性校正方法
CN203858360U (zh) * 2014-04-22 2014-10-01 南京大桥机器有限公司 一种雷达信号幅相一致性标校***
CN104931968A (zh) * 2015-06-18 2015-09-23 西安电子科技大学 基于FPGA的InSAR通道幅相误差估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Fast DOA estimation with ULA in the presence of sensor gain and phase errors;Wang, T.D，et al;《2009 International Conference on Communications, Circuits and Systems》;20091231;第395-397页 *
Two new estimation algorithms for sensor gain and phase errors based on different data models;J Jiang,et al;《IEEE Sensors Journal》;20131231;第1921-1930页 *

Also Published As

Publication number Publication date
CN105353338A (zh) 2016-02-24

Similar Documents

Publication Publication Date Title
CN104977558B (zh) 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法
CN104407335B (zh) 一种3轴交叉阵列的doa估计方法
CN107870315B (zh) 一种利用迭代相位补偿技术估计任意阵列波达方向方法
CN105353338B (zh) 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法
CN106646344A (zh) 一种利用互质阵的波达方向估计方法
CN104991236B (zh) 一种单基地mimo雷达非圆信号相干源波达方向估计方法
CN104035069B (zh) 基于部分矫正对称均匀线阵的窄带近场信号源定位方法
CN110515033B (zh) 一种基于Toeplitz矩阵恢复的欠通道测向***及方法
CN103353588B (zh) 基于天线均匀平面阵的二维波达方向角估计方法
CN110596687B (zh) 一种基于黎曼流形的单基地mimo雷达目标检测方法
CN111413668B (zh) 一种大规模阵列中基于dft增强的doa估计方法
CN103926555B (zh) 一种利用非圆信号测定天线阵列接收机幅相响应的方法
CN105242236B (zh) 宽带信号超分辨测向中的阵元位置误差校正方法
CN104793177B (zh) 基于最小二乘法的麦克风阵列测向方法
CN113835063B (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN111046591A (zh) 传感器幅相误差与目标到达角度的联合估计方法
CN104021293A (zh) 一种基于结构最小二乘法的联合到达角-频率估计方法
CN105572630B (zh) 基于多波位联合处理的单脉冲目标doa估计方法
CN105182280B (zh) 基于空域稀疏优化的宽带信号超分辨测向误差估计方法
CN105182279B (zh) 基于空域稀疏优化的宽带信号超分辨测向误差校正方法
CN104950282B (zh) 连续域内稀疏重构实现的宽带信号超分辨测向方法及装置
CN104156553A (zh) 无需信源数估计的相干信号波达方向估计方法及***
Dai et al. Enhanced second-order off-grid DOA estimation method via sparse reconstruction based on extended coprime array under impulsive noise
CN105223541B (zh) 宽带信号超分辨测向中的阵元间互耦误差校正方法
CN113589223B (zh) 基于互耦情况下嵌套阵列的测向方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171010

Termination date: 20200928