CN105182280A - 基于空域稀疏优化的宽带信号超分辨测向误差估计方法 - Google Patents
基于空域稀疏优化的宽带信号超分辨测向误差估计方法 Download PDFInfo
- Publication number
- CN105182280A CN105182280A CN201510628112.6A CN201510628112A CN105182280A CN 105182280 A CN105182280 A CN 105182280A CN 201510628112 A CN201510628112 A CN 201510628112A CN 105182280 A CN105182280 A CN 105182280A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- msup
- prime
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000005457 optimization Methods 0.000 title claims abstract description 13
- 230000008878 coupling Effects 0.000 claims abstract description 36
- 238000010168 coupling process Methods 0.000 claims abstract description 36
- 238000005859 coupling reaction Methods 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims description 101
- 239000013598 vector Substances 0.000 claims description 48
- 238000006880 cross-coupling reaction Methods 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 12
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 2
- 241000764238 Isis Species 0.000 claims 1
- 230000007850 degeneration Effects 0.000 claims 1
- 238000006467 substitution reaction Methods 0.000 claims 1
- 235000019800 disodium phosphate Nutrition 0.000 description 39
- 238000001514 detection method Methods 0.000 description 10
- 239000013078 crystal Substances 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000009795 derivation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
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/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
基于空域稀疏优化的宽带信号超分辨测向误差估计方法,涉及宽带信号超分辨测向中存在的阵列误差的估计方法。本发明为了解决现有的阵列误差的估计方法不适用于宽带信号的问题。本发明针对宽带信号超分辨测向时同时存在阵元间互耦、阵列通道幅相不一致性以及阵元位置误差时,利用各个频点上的信号构建对应的优化函数,之后利用信号的空域稀疏性,分别对各个频点上的函数进行迭代优化处理,最后对所有频点上的信息进行融合估计出阵元间互耦、阵列通道幅相不一致性以及阵元位置误差。该方法可以准确的估计出阵元间互耦、阵列通道幅相不一致性以及阵元位置误差。本发明适用于宽带信号超分辨测向中存在的阵列误差的估计。
Description
技术领域
本发明涉及宽带信号超分辨测向中存在的阵列误差的估计方法
背景技术
超分辨测向是阵列信号处理中的一个重要研究内容,广泛的应用于雷达、物联网和声呐等领域。常规的超分辨测向方法都需要准确的掌握阵列流型。而实际的测向系统当中往往存在高频振荡、放大器以及长度不一的通道等,并且有时会伴随着阵元位置的扰动,导致测向估计时经常伴随着阵元间互耦、阵列通道幅相不一致性以及阵元位置误差,这直接导致了很多的超分辨测向方法的性能恶化,甚至失效,所以有必要对各种误差进行校正,而校正前需要首先对这些误差进行估计。
FriedlanderB和WeissAJ基于子空间原理,提出了一种信源方位、阵元间互耦、阵元增益和相位扰动交替迭代估计的阵列误差估计与校正技术。Song采用盲估计方法,对声矢量传感器的幅相不一致性、方向不一致性以及阵元位置误差进行了迭代估计,得到了较好的效果。Wang等国外学者利用均匀线阵互耦误差矩阵的带状Toeplitz结构,提出了一种阵元间互耦误差估计方法。Yang针对ULA和均匀圆阵,基于凸松弛原理研究了互耦条件下误差估计问题。Liao通过将互耦误差转换成与角度相关的复阵列增益,提出了一种互耦误差与信号到达方向联合估计算法。然而以上方法都是针对于窄带信号超分辨测向误差的估计,对于宽带信号测向误差的估计,尤其是有关多种误差同时存在时的估计技术,公开发表的文献并不多见。
发明内容
本发明为了解决现有的阵列误差的估计方法不适用于宽带信号的问题。
基于空域稀疏优化的宽带信号超分辨测向误差估计方法,包括下述步骤:
步骤1:建立同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号模型:
当阵列中同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时,阵列输出可以表示为
X”'(fi)=A”'(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)A(fi,α)S(fi)+Λ(3)(fi)w(3)(fi)+N(fi),i=1,2,…,J(12)
=W(1)(fi)W(3)(fi,α)·A(fi,α)S(fi)+Λ(2)(fi)w(2)(fi)+N(fi)
=W(2)(fi)W(3)(fi,α)·A(fi,α)S(fi)+Λ(1)(fi)w(1)(fi)+N(fi)
其中,“·”表示Hadamard积,即两个矩阵对应位置的元素相乘,最后得到与原矩阵相同维数的矩阵,A″′(fi,α)为同时存在以上三种误差时频点fi上的阵列流型矩阵
A”'(fi,α)=[a”'(fi,α1),…,a”'(fi,αk),…,a”'(fi,αK)](13)
A(fi,α)=[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵,a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量;
A”'(fi,α)与理想情况下频点fi上的阵列流型矩阵的关系为
A”'(fi,α)=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)(14)
相应的同时存在以上三种误差时频点fi上第k个信号的阵列导向矢量为
a”'(fi,αk)=W(1)(fi)W(2)(fi)W(3)(fi,αk)·a(fi,αk)(15)
式(12)中,S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵;N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵,均值为0,方差为μ2(fi);则同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时频点fi上的接收信号协方差矩阵
R″′(fi)=E{X″′(fi)(X″′(fi))H},i=1,2,…,J(16)
W(1)(fi)为假设阵列当中只存在阵元间互耦误差时的阵列扰动矩阵,
其中,cq(fi)表示间距为q、信号频率为fi时的阵元间的互耦系数,q=1,2,…,Q;
定义阵元间互耦扰动矢量为w(1)(fi)=[c1(fi),…,cQ(fi)]T;Λ(1)(fi)为仅存在阵元间互耦误差时一个只与原信号有关的参数,与误差无关;
W(2)(fi)为假设阵列当中只存在阵列通道幅相不一致性误差时的幅相阵列扰动矩阵,
W(2)(fi)=diag([W1(fi),…,Wm(fi),…,WM(fi)]T)(18)
其中
为信号频点fi上第m路通道的幅相不一致性误差,ρm(fi)、分别为在信号频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差;
定义频点fi上阵列通道幅相不一致扰动矢量为Λ(2)(fi)为仅存在阵列通道幅相不一致性误差时频点fi上一个只与原信号有关的参数,与误差无关;
W(3)(fi,α)=[W(3)(fi,α1),…,W(3)(fi,αk),…,W(3)(fi,αK)]为假设阵列当中只存在阵元位置误差时频点fi上的阵元位置误差扰动矩阵,其中
为只存在阵元位置误差时频点fi上、第k个信号的阵元位置误差扰动矩阵,其中
为第k个信号到达第m个阵元时,由阵元位置误差引入的信源传播时延误差,Δdm为第m个阵元的真实位置与测量位置间的偏差,与信号频率fi无关;
定义频点fi上的阵元位置误差扰动矢量为w(3)(fi)=[Δd2,…,ΔdM]T;Λ(3)(fi)为仅存在阵元位置误差时一个只与原信号有关的参数,与误差无关;
步骤2:对同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号参数进行估计:
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向,可得出频点fi上阵列流型矩阵的稀疏表示
其中, 为频点fi上第l个稀疏信号的阵列导向矢量,同时得出阵元位置误差扰动矩阵的稀疏表示 其中 为频点fi上、第l个稀疏信号的阵元位置误差扰动矩阵,为第l个稀疏信号到达第m个阵元时,由阵元位置误差扰动引入的信源传播时延误差,
相应的可获得同时存在以上三种误差时频点fi上阵列流型矩阵的稀疏表示
其中, 为同时存在以上三种误差时频点fi上、第l个稀疏信号对应的阵列导向矢量,则可得出同时存在以上三种误差时频点fi上的阵列输出信号的稀疏表示
其协方差矩阵为
式(22)中 为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)估计出从而重构出原信号,同时对误差进行估计;
根据式(22)可知,同时存在以上三种误差时频点fi的阵列输出信号的概率密度为
IM为M×M维的单位阵;结合式(22)、(24)和(25)可得
采用期望最大化(ExpectationMaximization,EM)方法来对w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)进行迭代估计,得出估计值和对应的可得到 以及
本发明具有以下有益效果:
本发明提出了一种基于空域稀疏优化的宽带信号阵列误差估计方法,当阵元间互耦、阵列通道幅相不一致性以及阵元位置误差同时存在时,利用各个频点上的信号构建对应的优化函数,之后利用信号的空域稀疏性,分别对各个频点上的函数进行迭代优化处理,最后对所有频点上的信息进行融合估计出阵元间互耦、阵列通道幅相不一致性以及阵元位置误差,并且利用多片数字信号处理器对该方法进行实现,有效提高了该方法的运行速度。
当信噪比为10dB,每个频点采样快拍数为40时,阵元间互耦误差估计精度:实部0.2/σ,虚部0.2/σ;阵列通道幅相不一致性误差估计精度:实部0.15/σ,虚部0.15/σ;阵元位置误差估计精度:0.05d/σ(d为阵元间距)。
附图说明
图1为宽带信号超分辨测向阵列信号模型示意图;
图2为宽带信号探测系统装置图;
图3为具体实施方式五的宽带信号超分辨测向装置图;
图4为具体实施方式六的宽带信号超分辨测向装置图;
图5为具体实施方式七的宽带信号超分辨测向装置图。
具体实施方式
具体实施方式一:
基于空域稀疏优化的宽带信号超分辨测向误差估计方法,包括下述步骤:
步骤1:建立同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号模型:
当阵列中同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时,阵列输出可以表示为
X”'(fi)=A”'(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)A(fi,α)S(fi)+Λ(3)(fi)w(3)(fi)+N(fi),i=1,2,…,J(12)
=W(1)(fi)W(3)(fi,α)·A(fi,α)S(fi)+Λ(2)(fi)w(2)(fi)+N(fi)
=W(2)(fi)W(3)(fi,α)·A(fi,α)S(fi)+Λ(1)(fi)w(1)(fi)+N(fi)
其中,“·”表示Hadamard积,即两个矩阵对应位置的元素相乘,最后得到与原矩阵相同维数的矩阵,A″′(fi,α)为同时存在以上三种误差时频点fi上的阵列流型矩阵
A”'(fi,α)=[a”'(fi,α1),…,a”'(fi,αk),…,a”'(fi,αK)](13)
A(fi,α)=[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵,a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量;
A”'(fi,α)与理想情况下频点fi上的阵列流型矩阵的关系为
A”'(fi,α)=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)(14)
相应的同时存在以上三种误差时频点fi上第k个信号的阵列导向矢量为
a”'(fi,αk)=W(1)(fi)W(2)(fi)W(3)(fi,αk)·a(fi,αk)(15)
式(12)中,S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵;N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵,均值为0,方差为μ2(fi);则同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时频点fi上的接收信号协方差矩阵
R″′(fi)=E{X″′(fi)(X″′(fi))H},i=1,2,…,J(16)
W(1)(fi)为假设阵列当中只存在阵元间互耦误差时的阵列扰动矩阵,
其中,cq(fi)表示间距为q、信号频率为fi时的阵元间的互耦系数,q=1,2,…,Q;
定义阵元间互耦扰动矢量为w(1)(fi)=[c1(fi),…,cQ(fi)]T;Λ(1)(fi)为仅存在阵元间互耦误差时一个只与原信号有关的参数,与误差无关;
W(2)(fi)为假设阵列当中只存在阵列通道幅相不一致性误差时的幅相阵列扰动矩阵,
W(2)(fi)=diag([W1(fi),…,Wm(fi),…,WM(fi)]T)(18)
其中
为信号频点fi上第m路通道的幅相不一致性误差,ρm(fi)、分别为在信号频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差;
定义频点fi上阵列通道幅相不一致扰动矢量为Λ(2)(fi)为仅存在阵列通道幅相不一致性误差时频点fi上一个只与原信号有关的参数,与误差无关;
W(3)(fi,α)=[W(3)(fi,α1),…,W(3)(fi,αk),…,W(3)(fi,αK)]为假设阵列当中只存在阵元位置误差时频点fi上的阵元位置误差扰动矩阵,其中
为只存在阵元位置误差时频点fi上、第k个信号的阵元位置误差扰动矩阵,其中
为第k个信号到达第m个阵元时,由阵元位置误差引入的信源传播时延误差,Δdm为第m个阵元的真实位置与测量位置间的偏差,与信号频率fi无关;
定义频点fi上的阵元位置误差扰动矢量为w(3)(fi)=[Δd2,…,ΔdM]T;Λ(3)(fi)为仅存在阵元位置误差时一个只与原信号有关的参数,与误差无关;
步骤2:对同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号参数进行估计:
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向,可得出频点fi上阵列流型矩阵的稀疏表示
其中, 为频点fi上第l个稀疏信号的阵列导向矢量,同时得出阵元位置误差扰动矩阵的稀疏表示 其中 为频点fi上、第l个稀疏信号的阵元位置误差扰动矩阵,为第l个稀疏信号到达第m个阵元时,由阵元位置误差扰动引入的信源传播时延误差,
相应的可获得同时存在以上三种误差时频点fi上阵列流型矩阵的稀疏表示
其中, 为同时存在以上三种误差时频点fi上、第l个稀疏信号对应的阵列导向矢量,则可得出同时存在以上三种误差时频点fi上的阵列输出信号的稀疏表示
其协方差矩阵为
式(22)中 为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)估计出从而重构出原信号,同时对误差进行估计;
根据式(22)可知,同时存在以上三种误差时频点fi的阵列输出信号的概率密度为
IM为M×M维的单位阵;结合式(22)、(24)和(25)可得
采用期望最大化(ExpectationMaximization,EM)方法来对w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)进行迭代估计,得出估计值和对应的可得到 以及
具体实施方式二:
本实施方式步骤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个信号的阵列导向矢量,
a(fi,αk)=[1,exp(-jφk),…,exp(-j(M-1)φk)]T,i=1,2,…,J,(6)
其中,φ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)]Ti=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)=[N1(fi,kp),…,Nm(fi,kp),…,NM(fi,kp)]Ti=1,2,…,J(11)
为频点fi上的第kp次噪声采样矩阵,其中Nm(fi,kp)为第m个阵元在频点fi上得到的第kp次噪声采样值;
步骤1.2:在理想阵列信号模型基础上建立同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号模型:
当阵列中同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时,频点fi上的阵列输出可以表示为
X”'(fi)=A”'(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)S(fi)+N(fi)
=W(1)(fi)W(2)(fi)A(fi,α)S(fi)+Λ(3)(fi)w(3)(fi)+N(fi),i=1,2,…,J(12)
=W(1)(fi)W(3)(f,α)·A(fi,α)S(fi)+Λ(2)(fi)w(2)(fi)+N(fi)
=W(2)(fi)W(3)(f,α)·A(fi,α)S(fi)+Λ(1)(fi)w(1)(fi)+N(fi)
其中,“·”表示Hadamard积,即两个矩阵对应位置的元素相乘,最后得到与原矩阵相同维数的矩阵,A″′(fi,α)为同时存在以上三种误差时频点fi上的阵列流型矩阵
A”'(fi,α)=[a”'(fi,α1),…,a”'(fi,αk),…,a”'(fi,αK)](13)
它与理想情况下频点fi上的阵列流型矩阵的关系为
A”'(fi,α)=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)(14)
相应的同时存在以上三种误差时频点fi上第k个信号的阵列导向矢量为
a”'(fi,αk)=W(1)(fi)W(2)(fi)W(3)(fi,αk)·a(fi,αk)(15)
式(12)中S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵,N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵,均值为0,方差为μ2(fi);相应的协方差矩阵R″′(fi)为
R″′(fi)=E{X″′(fi)(X″′(fi))H},i=1,2,…,J(16)
W(1)(fi)为假设阵列当中只存在阵元间互耦误差时的阵列扰动矩阵,
其中,cq(fi)表示间距为q、信号频率为fi时的阵元间的互耦系数,q=1,2,…,Q;
定义阵元间互耦扰动矢量为w(1)(fi)=[c1(fi),…,cQ(fi)]T;Λ(1)(fi)为仅存在阵元间互耦误差时一个只与原信号有关的参数,与误差无关;
W(2)(fi)为假设阵列当中只存在阵列通道幅相不一致性误差时的幅相阵列扰动矩阵,
W(2)(fi)=diag([W1(fi),…,Wm(fi),…,WM(fi)]T)(18)
其中
为信号频点fi上第m路通道的幅相不一致性误差,ρm(fi)、分别为在信号频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差,;
定义阵列通道幅相不一致扰动矢量为Λ(2)(fi)为仅存在阵列通道幅相不一致性误差时一个只与原信号有关的参数,与误差无关;
W(3)(fi,α)=[W(3)(fi,α1),…,W(3)(fi,αk),…,W(3)(fi,αK)]为假设阵列当中只存在阵元位置误差时频点fi上的阵元位置误差扰动矩阵,其中
为只存在阵元位置误差时频点fi上、第k个信号的阵元位置误差扰动矩阵,其中
为第k个信号到达第m个阵元时,由阵元位置误差引入的信源传播时延误差;Δdm为第m个阵元的真实位置与测量位置间的偏差,与信号频率fi无关;
定义阵元位置误差扰动矢量为w(3)(fi)=[Δd2,…,ΔdM]T;Λ(3)(fi)为仅存在阵元位置误差时一个只与原信号有关的参数,与误差无关。
其它步骤和参数与具体实施方式一相同。
具体实施方式三:
本实施方式步骤2中所述的采用期望最大化方法来对w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)进行迭代估计的具体步骤如下:
在期望最大化方法中的E-step步中,首先对 的分布函数进行计算
其中运算符<·>表示求解条件期望;
在期望最大化方法中的M-step步中,分别求取分布函数 对各未知参数的导数,即对 取极值来对各未知参数求解;
分别令以上的导数为0,即可求得第p次迭代时各个未知参数的估计值
其中(p)代表迭代次数;
式(33)中
为矩阵第r1行、r2列的元素,其中tr[·]表示求迹运算;
式(38)中
O(1)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(39)
为中间变量;
Ξ(1)(fi)
(40)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)
为中间变量;
式(33)中
式(41)中,(Ψ(1)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
式(34)中
为矩阵第r1行、r2列的元素;
式(42)中
O(2)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(43)
为中间变量;
Ξ(2)(fi)
(44)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)
为中间变量;
式(34)中
式(45)中,(Ψ(2)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
式(35)中
为矩阵第r1行、r2列的元素;
式(46)中
O(3)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(47)
为中间变量;
Ξ(3)(fi)
(48)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)
为中间变量;
式(35)中
式(49)中,(Ψ(3)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
由于直接利用式(33)~(35)计算w(1)(fi)、w(2)(fi)、w(3)(fi)比较复杂,因此可将式(38)~(49)代入式(33)~(35)中对等式进行化简并对w(1)(fi)、w(2)(fi)、w(3)(fi)求解;
当迭代若干步后,w(1)(fi)、w(2)(fi)、w(3)(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、电源、晶振和复位。
数字信号处理器采用TexasInstruments(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再利用这些结果进行阵列校正。其中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、电源、晶振和复位。
数字信号处理器采用AnalogDeviceInstruments(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再利用这些结果进行阵列校正。其中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、电源、晶振和复位。
数字信号处理器采用AnalogDeviceInstruments(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再利用这些结果进行阵列校正。其中SRAM3-10负责存储数据,JTAG3-11负责对DSP(3-1~3-6)进行调试,电源负责整体供电,晶振负责提供时钟,复位负责提供复位信号。
Claims (4)
1.基于空域稀疏优化的宽带信号超分辨测向误差估计方法,其特征在于包括下述步骤:
步骤1:建立同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号模型:
当阵列中同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时,阵列输出表示为
其中,“·”表示Hadamard积,即两个矩阵对应位置的元素相乘,最后得到与原矩阵相同维数的矩阵,A″′(fi,α)为同时存在以上三种误差时频点fi上的阵列流型矩阵
A”'(fi,α)=[a”'(fi,α1),…,a”'(fi,αk),…,a”'(fi,αK)](13)
A(fi,α)=[a(fi,α1),…,a(fi,αk),…,a(fi,αK)]为理想情况下频点fi上的阵列流型矩阵,a(fi,αk)为理想情况下频点fi上第k个信号的阵列导向矢量;
A”'(fi,α)与理想情况下频点fi上的阵列流型矩阵的关系为
A”'(fi,α)=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)(14)
相应的同时存在以上三种误差时频点fi上第k个信号的阵列导向矢量为
a”'(fi,αk)=W(1)(fi)W(2)(fi)W(3)(fi,αk)·a(fi,αk)(15)
式(12)中,S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵;N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵,均值为0,方差为μ2(fi);则同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时频点fi上的接收信号协方差矩阵
R″′(fi)=E{X″′(fi)(X″′(fi))H},i=1,2,…,J(16)
W(1)(fi)为假设阵列当中只存在阵元间互耦误差时的阵列扰动矩阵,
其中,cq(fi)表示间距为q、信号频率为fi时的阵元间的互耦系数,q=1,2,…,Q;
定义阵元间互耦扰动矢量为w(1)(fi)=[c1(fi),…,cQ(fi)]T;Λ(1)(fi)为仅存在阵元间互耦误差时一个只与原信号有关的参数;
W(2)(fi)为假设阵列当中只存在阵列通道幅相不一致性误差时的幅相阵列扰动矩阵,
W(2)(fi)=diag([W1(fi),…,Wm(fi),…,WM(fi)]T)(18)
其中
为信号频点fi上第m路通道的幅相不一致性误差,ρm(fi)、分别为在信号频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差;
定义频点fi上阵列通道幅相不一致扰动矢量为Λ(2)(fi)为仅存在阵列通道幅相不一致性误差时频点fi上一个只与原信号有关的参数;
W(3)(fi,α)=[W(3)(fi,α1),…,W(3)(fi,αk),…,W(3)(fi,αK)]为假设阵列当中只存在阵元位置误差时频点fi上的阵元位置误差扰动矩阵,其中
为只存在阵元位置误差时频点fi上、第k个信号的阵元位置误差扰动矩阵,其中
为第k个信号到达第m个阵元时,由阵元位置误差引入的信源传播时延误差,△dm为第m个阵元的真实位置与测量位置间的偏差;
定义频点fi上的阵元位置误差扰动矢量为w(3)(fi)=[△d2,…,△dM]T;Λ(3)(fi)为仅存在阵元位置误差时一个只与原信号有关的参数;
步骤2:对同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号参数进行估计:
首先将搜索空间划分为若干离散的角度网格L表示信号可能到达的L个方向,得出频点fi上阵列流型矩阵的稀疏表示
其中, 为频点fi上第l个稀疏信号的阵列导向矢量,同时得出阵元位置误差扰动矩阵的稀疏表示 其中 为频点fi上、第l个稀疏信号的阵元位置误差扰动矩阵,为第l个稀疏信号到达第m个阵元时,由阵元位置误差扰动引入的信源传播时延误差,
相应的获得同时存在以上三种误差时频点fi上阵列流型矩阵的稀疏表示
其中, 为同时存在以上三种误差时频点fi上、第l个稀疏信号对应的阵列导向矢量,则得出同时存在以上三种误差时频点fi上的阵列输出信号的稀疏表示
其协方差矩阵为
式(22)中 为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)的高斯分布;
根据式(22),同时存在以上三种误差时频点fi的阵列输出信号的概率密度为
IM为M×M维的单位阵;结合式(22)、(24)和(25)得
采用期望最大化方法来对w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)进行迭代估计,得出估计值和对应的得到 以及
2.根据权利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个信号的阵列导向矢量,
a(fi,αk)=[1,exp(-jφk),…,exp(-j(M-1)φk)]T,i=1,2,…,J,(6)
其中,φ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)]Ti=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)=[N1(fi,kp),…,Nm(fi,kp),…,NM(fi,kp)]Ti=1,2,…,J(11)
为频点fi上的第kp次噪声采样矩阵,其中Nm(fi,kp)为第m个阵元在频点fi上得到的第kp次噪声采样值;
步骤1.2:在理想阵列信号模型基础上建立同时含有阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差的阵列信号模型:
当阵列中同时存在阵元间互耦误差、阵列通道幅相不一致性误差、阵元位置误差时,频点fi上的阵列输出表示为
其中,“·”表示Hadamard积,即两个矩阵对应位置的元素相乘,最后得到与原矩阵相同维数的矩阵,A″′(fi,α)为同时存在以上三种误差时频点fi上的阵列流型矩阵
A”'(fi,α)=[a”'(fi,α1),…,a”'(fi,αk),…,a”'(fi,αK)](13)
它与理想情况下频点fi上的阵列流型矩阵的关系为
A”'(fi,α)=W(1)(fi)W(2)(fi)W(3)(fi,α)·A(fi,α)(14)
相应的同时存在以上三种误差时频点fi上第k个信号的阵列导向矢量为
a”'(fi,αk)=W(1)(fi)W(2)(fi)W(3)(fi,αk)·a(fi,αk)(15)
式(12)中S(fi)为信号sk(t)经过傅立叶变换后的信号矢量矩阵,N(fi)为噪声nm(t)经过傅立叶变换后的噪声矢量矩阵,均值为0,方差为μ2(fi);相应的协方差矩阵R″′(fi)为
R″′(fi)=E{X″′(fi)(X″′(fi))H},i=1,2,…,J(16)
W(1)(fi)为假设阵列当中只存在阵元间互耦误差时的阵列扰动矩阵,
其中,cq(fi)表示间距为q、信号频率为fi时的阵元间的互耦系数,q=1,2,…,Q;
定义阵元间互耦扰动矢量为w(1)(fi)=[c1(fi),…,cQ(fi)]T;Λ(1)(fi)为仅存在阵元间互耦误差时一个只与原信号有关的参数;
W(2)(fi)为假设阵列当中只存在阵列通道幅相不一致性误差时的幅相阵列扰动矩阵,
W(2)(fi)=diag([W1(fi),…,Wm(fi),…,WM(fi)]T)(18)
其中
为信号频点fi上第m路通道的幅相不一致性误差,ρm(fi)、分别为在信号频点fi上第m路通道相对于第一路通道的幅度增益和相位偏差,;
定义阵列通道幅相不一致扰动矢量为Λ(2)(fi)为仅存在阵列通道幅相不一致性误差时一个只与原信号有关的参数;
W(3)(fi,α)=[W(3)(fi,α1),…,W(3)(fi,αk),…,W(3)(fi,αK)]为假设阵列当中只存在阵元位置误差时频点fi上的阵元位置误差扰动矩阵,其中
为只存在阵元位置误差时频点fi上、第k个信号的阵元位置误差扰动矩阵,其中
为第k个信号到达第m个阵元时,由阵元位置误差引入的信源传播时延误差;△dm为第m个阵元的真实位置与测量位置间的偏差;
定义阵元位置误差扰动矢量为w(3)(fi)=[△d2,…,△dM]T;Λ(3)(fi)为仅存在阵元位置误差时一个只与原信号有关的参数。
3.根据权利2所述的基于空域稀疏优化的宽带信号超分辨测向误差估计方法,其特征在于步骤2中所述的采用期望最大化方法来对w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)进行迭代估计的具体步骤如下:
在期望最大化方法中的E-step步中,首先对 的分布函数进行计算
其中运算符<·>表示求解条件期望;
在期望最大化方法中的M-step步中,分别求取分布函数 对各未知参数的导数,即对 取极值来对各未知参数求解;
分别令以上的导数为0,即求得第p次迭代时各个未知参数的估计值
其中(p)代表迭代次数;
式(33)中
为矩阵第r1行、r2列的元素,其中tr[·]表示求迹运算;
式(38)中
O(1)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(39)
为中间变量;
Ξ(1)(fi)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)(40)
为中间变量;
式(33)中
式(41)中,(Ψ(1)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
式(34)中
为矩阵第r1行、r2列的元素;
式(42)中
O(2)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(43)
为中间变量;
Ξ(2)(fi)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)(44)
为中间变量;
式(34)中
式(45)中,(Ψ(2)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
式(35)中
为矩阵第r1行、r2列的元素;
式(46)中
O(3)(fi)=Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1X”'(fi)(47)
为中间变量;
Ξ(3)(fi)
=Σ(fi)-Σ(fi)(A”'(fi,Ω))H(μ2(fi)IM+A”'(fi,Ω)Σ(fi)(A”'(fi,Ω))H)-1A”'(fi,Ω)Σ(fi)(48)
为中间变量;
式(35)中
式(49)中,(Ψ(3)(fi))r为中间变量,为M×M维的矩阵,只有在第±r对角线上的元素全为1,其余元素全为0;
将式(38)~(49)代入式(33)~(35)中对等式进行化简并对w(1)(fi)、w(2)(fi)、w(3)(fi)求解;
当迭代若干步后,w(1)(fi)、w(2)(fi)、w(3)(fi)、μ2(fi)和δl(fi)的估计值的变化趋于0,则得出最后的估计值和对应得到 以及
4.根据权利要求1、2或3所述的基于空域稀疏优化的宽带信号超分辨测向误差估计方法,其特征在于所述的ρ1(fi)=1,W1(fi)=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510628112.6A CN105182280B (zh) | 2015-09-28 | 2015-09-28 | 基于空域稀疏优化的宽带信号超分辨测向误差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510628112.6A CN105182280B (zh) | 2015-09-28 | 2015-09-28 | 基于空域稀疏优化的宽带信号超分辨测向误差估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105182280A true CN105182280A (zh) | 2015-12-23 |
CN105182280B CN105182280B (zh) | 2017-10-10 |
Family
ID=54904483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510628112.6A Expired - Fee Related CN105182280B (zh) | 2015-09-28 | 2015-09-28 | 基于空域稀疏优化的宽带信号超分辨测向误差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182280B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107037397A (zh) * | 2017-06-21 | 2017-08-11 | 哈尔滨工业大学 | 一种在波达方向估计中校正多种阵列误差的方法 |
CN108880586A (zh) * | 2018-06-28 | 2018-11-23 | 中国人民解放军战略支援部队信息工程大学 | 一种宽带弱信号增强方法与装置 |
CN111123273A (zh) * | 2019-12-24 | 2020-05-08 | 浙江大学 | 基于贝叶斯压缩感知算法的稀疏阵列优化方法 |
RU2814220C1 (ru) * | 2023-07-04 | 2024-02-28 | Общество с ограниченной ответственностью "Специальный Технологический Центр" (ООО "СТЦ") | Способ обнаружения и оценивания характеристик широкополосных сигналов и устройство для его реализации |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011093805A1 (en) * | 2010-01-27 | 2011-08-04 | Tuncer Engin T | A system and a method for simultaneous position, mutual coupling and gain/phase calibration of antenna arrays |
CN103885058A (zh) * | 2014-02-26 | 2014-06-25 | 西安电子科技大学 | 一种利用稀疏统计特性的机载雷达前视超分辨成像方法 |
CN103954961A (zh) * | 2014-05-13 | 2014-07-30 | 西安电子科技大学 | 步进频雷达低空障碍物超分辨一维成像方法 |
-
2015
- 2015-09-28 CN CN201510628112.6A patent/CN105182280B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011093805A1 (en) * | 2010-01-27 | 2011-08-04 | Tuncer Engin T | A system and a method for simultaneous position, mutual coupling and gain/phase calibration of antenna arrays |
CN103885058A (zh) * | 2014-02-26 | 2014-06-25 | 西安电子科技大学 | 一种利用稀疏统计特性的机载雷达前视超分辨成像方法 |
CN103954961A (zh) * | 2014-05-13 | 2014-07-30 | 西安电子科技大学 | 步进频雷达低空障碍物超分辨一维成像方法 |
Non-Patent Citations (4)
Title |
---|
WANG DING ET AL.: "Sensor array calibration for uniform rectangular array in presence of mutual coupling and sensor gain-and-phase errors", 《JOURNAL OF CENTRAL SOUTH UNIV ERSITY》 * |
YONG YANG ET AL.: "A calibration algorithm of millimeter-wave sparse arrays based on simulated annealing", 《MICROWAVE AND MILLIMETER WAVE TECHNOLOGY (ICMMT)》 * |
刘章孟: "基于信号空域稀疏性的阵列处理理论与方法", 《中国博士学位论文全文数据库 信息科技辑》 * |
王鼎等: "一种利用互耦矩阵稀疏性的阵列误差有源校正改进算法", 《信号处理》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107037397A (zh) * | 2017-06-21 | 2017-08-11 | 哈尔滨工业大学 | 一种在波达方向估计中校正多种阵列误差的方法 |
CN108880586A (zh) * | 2018-06-28 | 2018-11-23 | 中国人民解放军战略支援部队信息工程大学 | 一种宽带弱信号增强方法与装置 |
CN108880586B (zh) * | 2018-06-28 | 2019-11-08 | 中国人民解放军战略支援部队信息工程大学 | 一种宽带弱信号增强方法与装置 |
CN111123273A (zh) * | 2019-12-24 | 2020-05-08 | 浙江大学 | 基于贝叶斯压缩感知算法的稀疏阵列优化方法 |
RU2814220C1 (ru) * | 2023-07-04 | 2024-02-28 | Общество с ограниченной ответственностью "Специальный Технологический Центр" (ООО "СТЦ") | Способ обнаружения и оценивания характеристик широкополосных сигналов и устройство для его реализации |
Also Published As
Publication number | Publication date |
---|---|
CN105182280B (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444810B (zh) | 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法 | |
CN105259550B (zh) | 基于压缩感知的多输入多输出雷达二维角度估计方法 | |
CN107290730A (zh) | 互耦条件下双基地mimo雷达角度估算方法 | |
CN104407335B (zh) | 一种3轴交叉阵列的doa估计方法 | |
CN106980104B (zh) | 用于传感器阵列的信号波达方向自校正方法 | |
CN109782238B (zh) | 一种传感器阵列阵元幅相响应和阵元位置的联合校准方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN111046591B (zh) | 传感器幅相误差与目标到达角度的联合估计方法 | |
CN105242236B (zh) | 宽带信号超分辨测向中的阵元位置误差校正方法 | |
CN104793177B (zh) | 基于最小二乘法的麦克风阵列测向方法 | |
CN105353338B (zh) | 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法 | |
CN105182280B (zh) | 基于空域稀疏优化的宽带信号超分辨测向误差估计方法 | |
Semper et al. | Grid-free direction-of-arrival estimation with compressed sensing and arbitrary antenna arrays | |
CN105182279B (zh) | 基于空域稀疏优化的宽带信号超分辨测向误差校正方法 | |
Gu et al. | Fast and efficient DOA estimation method for signals with known waveforms using nonuniform linear arrays | |
Dai et al. | Enhanced second-order off-grid DOA estimation method via sparse reconstruction based on extended coprime array under impulsive noise | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN104950282B (zh) | 连续域内稀疏重构实现的宽带信号超分辨测向方法及装置 | |
CN105223541B (zh) | 宽带信号超分辨测向中的阵元间互耦误差校正方法 | |
CN104635199A (zh) | 块稀疏重构方法实现的宽带信号快速超分辨测向装置 | |
CN112763972A (zh) | 基于稀疏表示的双平行线阵二维doa估计方法及计算设备 | |
CN105572631A (zh) | 基于多波位联合处理的最大似然目标doa估计方法 | |
Shi et al. | Direction of arrival estimation for a non-ideal acoustic vector hydrophone array | |
Park et al. | Alternating projections gridless covariance-based estimation for DOA | |
Li et al. | DOA estimation for sparse nested MIMO radar with velocity receive sensor array |
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 |
Granted publication date: 20171010 Termination date: 20200928 |
|
CF01 | Termination of patent right due to non-payment of annual fee |