CN102208930A - 快速信号子空间估计方法 - Google Patents
快速信号子空间估计方法 Download PDFInfo
- Publication number
- CN102208930A CN102208930A CN2010101352052A CN201010135205A CN102208930A CN 102208930 A CN102208930 A CN 102208930A CN 2010101352052 A CN2010101352052 A CN 2010101352052A CN 201010135205 A CN201010135205 A CN 201010135205A CN 102208930 A CN102208930 A CN 102208930A
- Authority
- CN
- China
- Prior art keywords
- matrix
- signal
- observation data
- vector
- iteration
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及数字信号处理技术领域。本发明提供了一种快速信号子空间估计方法,所述方法包括:根据观测数据模型,获得接收天线阵元的初始观测信号;构造新的观测数据,对所述新的观测数据进行迭代;根据所述迭代结果构造转换矩阵;以及,获得所述转换矩阵的列向量。本发明的方法在数据层次上进行多级维纳滤波前向迭代,得到前向迭代过程的转换矩阵,进而利用Householder矩阵的酉性质快速地计算该转换矩阵,从而获得信号子空间的一组标准基向量;所述方法借助了多级维纳滤波前向迭代的Householder变换,不仅进一步降低了计算复杂度,还保证了算法在有限精度下的数值稳健性。
Description
技术领域
本发明涉及数字信号处理技术领域,特别涉及一种降低信号子空间估计计算复杂度并保证数值稳健性的方法。
背景技术
随着大规模集成电路和数字信号处理技术的发展,信号子空间估计问题已经成为许多领域的关键问题之一,如相控阵雷达的波达方向(DOA)估计和降秩空时自适应处理(STAP)、移动通信的多用户检测(Multi-user Detection)等。传统的信号子空间估计方法通常是对观测数据协方差矩阵进行特征值分解(EVD),对应大特征值的特征向量张成信号子空间,其余特征向量张成噪声子空间。然而由于特征值分解的计算复杂度非常高,数量级为O(M3),其中M为观测数据的维数,并且需要计算观测数据协方差矩阵,计算复杂度为O(M2N),其中N为形成协方差矩阵的采样支持长度。因此利用特征值分解估计信号子空间的计算复杂度为O(M3+M2N)。然而在很多实际应用中,特别是相控阵雷达和第三代移动通信系统,观测数据维数M通常很大,因此巨大的计算量往往限制了实时应用。
学术界和产业界相继提出了一系列的信号子空间快速估计方法来避免计算特征值分解,如G.Xu提出了基于Lanczos的快速子空间估计方法,其计算复杂度为O(M2P+M2N),其中P为信号子空间的维数。当P远小于M时,该方法的计算量得到一定程度减轻,然而其仍然需要计算协方差矩阵。另外一种方法是基于多级维纳滤波(MSWF)前向迭代来估计信号子空间。多级维纳滤波前向迭代可以在数据层次(data-level)或者协方差矩阵层次(power-level)实现。协方差矩阵层次的前向迭代与G.Xu的方法的计算量一致,而数据层次的前向迭代,特别是基于相关相减结构的阻塞矩阵选择,能够大大降低其计算复杂度,然而其数值稳健性无法得到保证,尤其在有限精度运算中。
发明内容
本发明的目的是克服上述缺陷,提出一种数值稳健的快速信号子空间估计方法。
为实现上述目的,本发明的快速子空间估计方法包括:
获得天线阵元的观测数据模型;
构造新的观测数据,对所述新的观测数据进行迭代;
根据所述迭代结果构造转换矩阵;以及
获得所述转换矩阵的列向量。
进一步地,所述方法在所述构造新的观测数据,对所述新的观测数据进行迭代的步骤之前还包括:
根据所述观测数据模型获得天线阵元接收信号数据;
根据所述天线阵元接收信号数据构造新的期望信号和新的观测数据;
其中,所述观测数据模型为:
x(k)=[x1(k),x2(k),…,xM(k)]T,k=1,2,…,N
=A(θ)s(k)+n(k)
=[a(θ1),a(θ2),…,a(θP)][s1(k),s2(k),…,sP(k)]T
+[n1(k),n2(k),…,nM(k)]T
其中,xM(k)代表第M个天线阵元的信号数据;s(k),A(θ)和n(k)分别表示信号矢量、导向矩阵以及噪声矢量;N为采样支持长度,a(θi)代表波达方向为θi的信号导向矢量,并且其具有如下形式:
进一步地,所述对所述新的观测数据进行迭代的步骤通过多级维纳滤波前向迭代计算实现;
进一步地,所述迭代结果为:归一化互相关矢量hi和阻塞矩阵QB,i;
进一步地,所述构造后的转换矩阵为:
其中,ki为具有单位范数,任意相位的旋转因子;hi为归一化互相关矢量;QBi为阻塞矩阵;qi为Householder矩阵Qi的第一列;
进一步地,所述获得所述转换矩阵的列向量的步骤通过按照多级维纳滤波前向迭代的Householder矩阵后向连乘来实现;
进一步地,所述新的期望信号为天线参考阵元的观测数据,即d0(k)=x1(k);所述新的观测数据为:x0(k)=[x2(k),x3(k),…,xM(k)]T,其中,xM(k)代表第M个天线阵元的信号数据;
进一步地,所述多级维纳滤波前向迭代计算在所述新的观测数据的数据层次上进行。
本发明的方法不需要计算观测数据协方差矩阵或者特征值分解,而是在数据层次上进行多级维纳滤波前向迭代,得到前向迭代过程的转换矩阵,进而利用Householder矩阵的酉性质快速地计算该转换矩阵,从而获得信号子空间的一组标准基向量;所述方法借助了多级维纳滤波前向迭代的Householder变换,不仅进一步降低了计算复杂度,还保证了算法在有限精度下的数值稳健性。
附图说明
图1是本发明具体实施例的ULA天线阵观测数据模型示意图;
图2是本发明具体实施例的基于Householder变换的多级维纳滤波前向迭代的示意图;
图3是M=20,K=19,N=4K时本发明具体实施例与传统算法的乘法计算量比较示意图;
图4是M=40,K=39,N=4K时本发明具体实施例与传统算法的乘法计算量比较示意图;
具体实施方式
本发明提出的快速子空间估计方法,结合附图和实施例说明如下。
S101:获得天线阵元的观测数据模型;
M个单元的均匀天线阵(ULA)接收到P个窄带远场信号,在第k个快拍时刻,M×1维的复观测数据可以表示为:
x(k)=[x1(k),x2(k),…,xM(k)]T,k=1,2,…,N
=A(θ)s(k)+n(k) (1)
=[a(θ1),a(θ2),…,a(θP)][s1(k),s2(k),…,sP(k)]T
+[n1(k),n2(k),…,nM(k)]T
其中,xM(k)代表第M个天线阵元的信号数据;s(k),A(θ)和n(k)分别表示信号矢量、导向矩阵以及噪声矢量;N为采样支持长度,a(θi)代表波达方向为θi的信号导向矢量,并且其具有如下形式:
首先,对于接收到的信号,假设其为窄带信号,具有不同的波达方向θi,并且所有信号为零均值的联合平稳高斯过程,即E{s(k)}=0;此外,对于噪声,假设其为时域白、零均值的平稳高斯过程,即E{n(k)}=0;并且有:信号之间不相关,信号和噪声之间不相关,不同阵元的噪声之间亦不相关,即有E{s(k)sH(k)}=Rs,E{s(k)nH(k)}=0,E{n(k)nH(k)}=σ2I。由此,式(1)构成了天线阵元接收信号的基本观测模型。
如图1所示为本发明实施例的天线阵接收数据的观测数据模型;
S102:构造新的观测数据,对所述新的观测数据进行迭代;
更具体地:
S1021:根据所述观测数据模型获得天线阵元接收信号数据;
S1022:根据所述天线阵元接收信号数据构造新的期望信号和新的观测数据;
S1023:对所述新的观测数据进行迭代。
本实施例采用基于Householder变换的多级维纳滤波前向迭代,所述迭代过程基于数据层次实现,避免了观测数据协方差矩阵的计算,从而降低了计算复杂度。
在具体实施过程中,为了能够执行多级维纳滤波前向迭代,在进行迭代计算之前,必须通过各个阵元接收的观测数据得到先验的本地期望信号(也称为训练信号)d0(k)以及输入数据向量x0(k);在本实施例中,将d0(k)设置为天线参考阵元的观测数据x1(k)(第一个阵元数据),即d0(k)=x1(k);而多级维纳滤波前向迭代的输入数据向量x0(k)则表示为x0(k)=[x2(k),x3(k),…,xM(k)]T,迭代计算在此新的观测数据基础上进行。
经过多次的前向迭代后,可以得到每级互相关矢量ri、范数归一化互相关矢量hi、Householder矩阵Qi、Householder向量qi和阻塞矩阵QBi,它们之间存在如下的数学关系,其中ki是具有单位范数的复相位旋转因子,即,即||ki||2=1:
hi=ri/||ri||2 (3)
QB,ihi=QB,iqi=0
在基于Householder变换的多级维纳滤波前向迭代中,利用hi、||ki||2=1以及的关系,可计算出各级的阻塞矩阵QB,i。
如图2所示,在具体实施过程中,可以通过以下方法完成基于Householder变换的多级维纳滤波前向迭代:
for i=1:P
hi=ri/δi
ki=-sign(real(hi,1))
vi=hi-kiδiui (4)
end
从i=1到P执行上述过程,共进行P次迭代。
S103:根据所述迭代的结果构造转换矩阵;
假设多级维纳滤波前向迭代进行P次后得到的转换矩阵由TP表示,它直接作用于输入的期望信号d0(k)和输入数据向量x0(k),从而可得到一组期望信号向量d(k)=[d0(k),d1(k),…,dP-1(k)]T,并且期望信号向量d(k)的协方差矩阵为三对角矩阵,即其中Rx和Rd分别表示输入数据向量x0(k)和输出期望信号向量d(k)的协方差矩阵,且Rd为三对角矩阵。基于步骤S102计算得到的归一化互相关矢量hi和阻塞矩阵QB,i、图2所示的多级维纳滤波前向迭代过程以及前述的||ki||2=1,转换矩阵TP可以表示如下:
其中,其中,ki为具有单位范数,任意相位的旋转因子;hi为归一化互相关矢量;QB,i为阻塞矩阵;qi为Householder矩阵Qi的第一列。
TP的列向量张成了信号子空间,并且这组列向量是信号子空间的一组标准正交基,所以,通过计算所述转换矩阵并获得其列向量即可实现信号子空间估计。
在基于Householder变换的多级维纳滤波前向迭代中,Householder矩阵Qi具有严格的酉性质,即因此酉变换减小了有限精度运算下的舍入误差在前向迭代过程中的传播,TP列向量(即信号子空间的标准正交基)的正交性得到了非常好的保证,理论上可以证明其正交性误差非常接近机器的运算精度,即有限精度运算下的正交性误差与2-ε在同一个数量级,其中ε为定点数中表示小数部分的比特长度,因此子空间估计的数值稳健性得到显著增强。
S104:获得所述转换矩阵的列向量;
本实施例采用转换矩阵TP的快速计算方法获得其列向量,从而大大减小了信号子空间估计的计算量。利用Householder矩阵Qi的酉性质从数学上可以证明:转换矩阵TP是矩阵Q的前P列,并且Q可表示为如下一系列矩阵的连乘:
其中矩阵的(k-1)×(k-1)主子式为单位矩阵Ik-1。计算矩阵有两种方法:前向连乘和后向连乘。在矩阵的后向连乘中(即第一步先计算,再计算,依此类推),中间矩阵的(k-1)×(k-1)主子式逐渐消失。利用该性质,本实施例的转换矩阵TP的快速计算方法为:不按照多级维纳滤波前向迭代的顺序前向连乘来计算矩阵(即第一步先计算,再计算,依此类推),而是通过其后向连乘来进一步降低计算量。
具体实施过程中,可以通过以下方法完成转换矩阵TP的快速计算:
初始化:Q=I
连乘:for i=(P-1):1
end
截取:TP=Q(:,1:P)
其中,K=M-1,I为单位矩阵;Qi为MSWF前向迭代过程中的Householder矩阵;vi为Householder向量;βi为系数因子。
通过上述方法得到转换矩阵TP的计算结果,从而获得其列向量,所述列向量即为估计信号子空间的一组标准正交基。
传统的基于相关相减结构的多级维纳滤波亦可得到信号子空间的一组标准正交基,并且这组基向量即为范数归一化互相关矢量hi。然而在基于相关相减结构的多级维纳滤波前向迭代中,阻塞矩阵Bi不具有酉性质,即将导致有限精度运算下的舍入误差通过多级维纳滤波前向迭代而传播,舍入误差产生积累,最终使得信号子空间基向量的正交性弱化,算法产生数值不稳健。本发明给出的方法中,数值稳健性通过Householder变换得到非常好的保证,因为Householder变换具有酉变换性质,多级维纳滤波的前向迭代不传播舍入误差,并且信号子空间基向量的正交性得到保证,非常接近于机器的运算精度。此外,通过步骤S104的转换矩阵TP的快速计算方法,其计算量亦可以得到进一步降低。
表1所示为本发明的计算方法、传统相关相减结构的多级维纳滤波方法以及基于公式(5)直接计算转换矩阵TP这三种方法的乘法计算量比较(实际应用中通常不考虑加法计算量),其中K=M-1。表1的乘法计算量公式难以明显比较三种方法的计算量。
表1
如图3、4所示为具体实施过程中,当M=20或40(K=19或39)以及N=4K时本发明的方法与其余两种传统方法的乘法计算量比较。由图3和4可以看出,本发明提出的信号子空间估计方法,可以大大减小计算量,尤其对于天线阵接收信号较多、多级维纳滤波前向迭代次数大的环境。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (10)
1.一种快速信号子空间估计方法,其特征在于,所述方法包括如下步骤:
(1)获得天线阵元的观测数据模型;
(2)构造新的观测数据,对所述新的观测数据进行迭代;
(3)根据所述迭代结果构造转换矩阵;以及
(4)获得所述转换矩阵的列向量。
2.如权利要求1所述的方法,其特征在于,在步骤(1)和(2)之间还包括如下步骤:
根据所述观测数据模型获得天线阵元的初始接收信号数据;
根据所述天线阵元的初始接收信号数据构造新的期望信号和新的观测数据。
4.如权利要求1所述的方法,其特征在于,所述步骤(2)通过多级维纳滤波前向迭代计算实现。
5.如权利要求1所述的方法,其特征在于,所述迭代结果为:归一化互相关矢量hi和阻塞矩阵QB,i。
6.如权利要求1所述的方法,其特征在于,所述构造后的转换矩阵为:
其中,κi为具有单位范数,任意相位的旋转因子;hi为归一化互相关矢量;QB,i为阻塞矩阵;qi为Householder矩阵Qi的第一列。
7.如权利要求1所述的方法,其特征在于,所述获得所述转换矩阵的列向量的步骤通过按照多级维纳滤波前向迭代的Householder矩阵后向连乘来实现。
8.如权利要求2所述的方法,其特征在于,所述新的期望信号为天线参考阵元的观测数据,即d0(k)=x1(k)。
9.如权利要求2所述的方法,其特征在于,所述新的观测数据为:x0(k)=[x2(k),x3(k),…,xM(k)]T,其中,xM(k)代表第M个天线阵元的信号数据。
10.如权利要求4所述的方法,其特征在于,所述多级维纳滤波前向迭代计算在所述新的观测数据的数据层次上进行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010135205.2A CN102208930B (zh) | 2010-03-30 | 2010-03-30 | 快速信号子空间估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010135205.2A CN102208930B (zh) | 2010-03-30 | 2010-03-30 | 快速信号子空间估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102208930A true CN102208930A (zh) | 2011-10-05 |
CN102208930B CN102208930B (zh) | 2014-01-22 |
Family
ID=44697611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010135205.2A Active CN102208930B (zh) | 2010-03-30 | 2010-03-30 | 快速信号子空间估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102208930B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005038A (zh) * | 2015-06-18 | 2015-10-28 | 许昌学院 | 一种改进的声矢量阵相干源doa估计算法 |
CN105158751A (zh) * | 2015-08-29 | 2015-12-16 | 许昌学院 | 一种声矢量阵快速方位估计方法 |
CN107179521A (zh) * | 2017-06-01 | 2017-09-19 | 电子科技大学 | 一种低复杂度的鲁棒子空间估计方法 |
CN108226849A (zh) * | 2018-01-30 | 2018-06-29 | 电子科技大学 | 一种基于gpu的子空间快速求解方法 |
CN110719123A (zh) * | 2019-09-21 | 2020-01-21 | 广东石油化工学院 | 一种利用子空间优化理论的plc信号重构方法和系统 |
CN114721015A (zh) * | 2022-04-08 | 2022-07-08 | 中山大学 | 一种gnss接收机盲稳健stap波束形成方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271021A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种快速得到信号子空间的信号处理方法 |
WO2008139865A1 (ja) * | 2007-05-09 | 2008-11-20 | Nec Corporation | Mmse等化回路、受信装置、及び通信システム、並びにその方法及びプログラム |
CN101592721A (zh) * | 2009-06-19 | 2009-12-02 | 中国人民解放军理工大学 | 基于特征值重构的相干信号到达角估计方法 |
CN101626354A (zh) * | 2009-05-20 | 2010-01-13 | 高拓讯达(北京)科技有限公司 | 多径信道的噪声估计方法和装置 |
CN101662321A (zh) * | 2008-08-25 | 2010-03-03 | 上海无线通信研究中心 | 认知无线mimo系统中基于零子空间技术的次级预编码发送方法 |
-
2010
- 2010-03-30 CN CN201010135205.2A patent/CN102208930B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271021A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种快速得到信号子空间的信号处理方法 |
WO2008139865A1 (ja) * | 2007-05-09 | 2008-11-20 | Nec Corporation | Mmse等化回路、受信装置、及び通信システム、並びにその方法及びプログラム |
CN101662321A (zh) * | 2008-08-25 | 2010-03-03 | 上海无线通信研究中心 | 认知无线mimo系统中基于零子空间技术的次级预编码发送方法 |
CN101626354A (zh) * | 2009-05-20 | 2010-01-13 | 高拓讯达(北京)科技有限公司 | 多径信道的噪声估计方法和装置 |
CN101592721A (zh) * | 2009-06-19 | 2009-12-02 | 中国人民解放军理工大学 | 基于特征值重构的相干信号到达角估计方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005038A (zh) * | 2015-06-18 | 2015-10-28 | 许昌学院 | 一种改进的声矢量阵相干源doa估计算法 |
CN105158751A (zh) * | 2015-08-29 | 2015-12-16 | 许昌学院 | 一种声矢量阵快速方位估计方法 |
CN107179521A (zh) * | 2017-06-01 | 2017-09-19 | 电子科技大学 | 一种低复杂度的鲁棒子空间估计方法 |
CN107179521B (zh) * | 2017-06-01 | 2019-06-18 | 电子科技大学 | 一种低复杂度的鲁棒子空间估计方法 |
CN108226849A (zh) * | 2018-01-30 | 2018-06-29 | 电子科技大学 | 一种基于gpu的子空间快速求解方法 |
CN108226849B (zh) * | 2018-01-30 | 2019-10-18 | 电子科技大学 | 一种基于gpu的子空间快速求解方法 |
CN110719123A (zh) * | 2019-09-21 | 2020-01-21 | 广东石油化工学院 | 一种利用子空间优化理论的plc信号重构方法和系统 |
CN114721015A (zh) * | 2022-04-08 | 2022-07-08 | 中山大学 | 一种gnss接收机盲稳健stap波束形成方法及装置 |
CN114721015B (zh) * | 2022-04-08 | 2022-12-09 | 中山大学 | 一种gnss接收机盲稳健stap波束形成方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102208930B (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102208930B (zh) | 快速信号子空间估计方法 | |
Xu et al. | Distributed estimation over sensor networks based on distributed conjugate gradient strategies | |
CN105335615B (zh) | 一种低复杂度的二维角度和极化参数联合估计方法 | |
Mao et al. | An Improved DOA Estimation Algorithm Based on Wavelet Operator. | |
Zheng et al. | Two-dimensional DOA estimation for coprime planar array: A coarray tensor-based solution | |
Li et al. | Distributed MIMO radar based on sparse sensing: Analysis and efficient implementation | |
CN101483280B (zh) | 一种稳健波束合成器权值求解方法 | |
Zheng et al. | Robust and efficient adaptive beamforming using nested subarray principles | |
CN114884841A (zh) | 基于高阶统计和非均匀阵列的欠定参数联合估计方法 | |
CN115421119A (zh) | 一种基于收发翻转互质mimo雷达结构的doa估计方法 | |
Pan et al. | Direction-of-arrival estimation with ULA: A spatial annihilating filter reconstruction perspective | |
CN109283486A (zh) | 基于广义复数相关熵的相干分布式非圆信号角度参数联合估计方法 | |
CN107645460A (zh) | 实值平行因子分解的多径参数估算方法 | |
CN109507634A (zh) | 一种任意传感器阵列下的基于传感算子的盲远场信号波达方向估计方法 | |
CN115913408A (zh) | 天线相位误差校正方法、装置、基站和存储介质 | |
Chen et al. | DOA and noncircular phase estimation of noncircular signal via an improved noncircular rotational invariance propagator method | |
CN106534028B (zh) | 基于压缩感知的ofdm信道跟踪方法 | |
CN104777492A (zh) | 卫星导航系统抗干扰处理方法及装置 | |
Jeon et al. | Superresolution TOA estimation with computational load reduction | |
Li et al. | DOA Estimation Based on Sparse Reconstruction via Acoustic Vector Sensor Array under Non-uniform Noise | |
CN109283485A (zh) | 一种基于共变的相干分布式非圆信号角度参数联合估计方法 | |
Shi et al. | Direct Position Determination of Non-circular Signals For Distributed Antenna Arrays: Optimal Weight And Polynomial Rooting Approach | |
Leite et al. | Partially-Calibrated Sparse Subarrays for Direction-of-Arrival Estimation with Increased Degrees of Freedom | |
Pinto et al. | Multi-step knowledge-aided iterative MUSIC for direction finding using nested arrays | |
CN114785383B (zh) | 一种基于智能超表面辅助1比特adc通信系统的导频图案设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |