一种同步直扩信号的源数估计方法技术

本发明专利技术提供了一种同步直扩信号的源数估计方法，通过将每个阵元的接收信号分段重构，使得重构后的观测矩阵可以重新结合经典的信源数目估计方法，避免了在波达方向重合时经典方法失效的问题，并提出了一种新的判决准则，能够改善比值法的性能，进一步提高了信源估计的准确率，通过重构观测矩阵，避免了导向矩阵秩亏；最后提出了用特征值比值的导数作为判决函数，改善了比值法的效果，提高了信源估计的准确率。本发明专利技术实现了在信源波达方向重合时的信源数估计，该方法同时适用于白噪声和有色噪声背景，在较低的信噪比下也具有良好的估计性能，同时该方法在欠定的情况下也能够稳定地实现信源数估计。现信源数估计。现信源数估计。

【技术实现步骤摘要】
一种同步直扩信号的源数估计方法


[0001]本专利技术涉及信号处理领域，具体为一种判断源信号数目的参数估计方法。

技术介绍

[0002]阵列信号处理作为信号处理领域的一个重要分支，其在雷达、语音处理、无线通信、射电天文和医学诊断等领域有着十分重要且广泛的应用。相比使用单个传感器，利用阵列具有获得更多的空间信息和更高的信号增益等优势，因此阵列信号处理技术近几十年来一直在蓬勃发展。
[0003]在阵列信号处理中信源数目往往假设已知，许多经典的DOA方法如MUSIC和ESPRIT等需要依据信源数目才能够正确地划分信号子空间和噪声子空间，从而估计出目标的方向；在语音盲源分离领域，经典的盲源分离方法如Fast
‑
ICA、JADE等也需要根据信源数目才能够成功地分离出各个源信号；在直扩信号系统参数估计中，伪码周期和伪码序列等估计方法往往也需要信源数目这一先验信息；信源数目的不准估计通常会直接影响以上方法的性能结果。
[0004]在众多信源数估计方法中，基于信息论的Akaike信息论准则(AIC)和最小描述长度准则(MDL)是最为常见的方法，但是该方法不适用于色噪声背景，Hsien
‑
Tsai Wu等提出的盖世圆盘法(GDE)根据盖世圆盘的半径划分信号和噪声子空间，能够实现色噪声下的信源数估计，但低信噪比下表现不佳。A.P.Liavas等基于对数函数重新拟合特征值减少了噪声的离散程度，能够在色噪声小快拍下估计信源数。D.T.Hoang利用对数特征值构建LogECNet分类网络进一步提高了估计精度，但是在非合作情况下获取数据困难使得其应用受到限制。李鑫凯等结合平滑秩方法和盖氏圆方法研究了多用户周期和非周期长码的信源数估计问题。以上对信源数估计的研究大多都没有考虑信源之间波达方向很近时的情况，此时导向矩阵不满足列满秩的条件，会严重影响源数估计方法的性能，当信源之间的波达方向相同时，经典的信源估计方法会失效。

技术实现思路

[0005]为了克服现有技术的不足，本专利技术提供一种同步直扩信号的源数估计方法。针对信源波达方向相同时同步直扩信号的信源数估计问题，本专利技术提出了一种分段重构的方法，通过将每个阵元的接收信号分段重构，使得重构后的观测矩阵可以重新结合经典的信源数目估计方法，避免了在波达方向重合时经典方法失效的问题，并提出了一种新的判决准则，能够改善比值法的性能，进一步提高了信源估计的准确率。
[0006]针对同步直扩信号信源数估计方法在信源之间的波达方向很近时性能下降甚至失效的问题，本专利技术提出了一种针对同步直扩信号的信源数估计方法，通过重构观测矩阵，避免了导向矩阵秩亏；最后提出了用特征值比值的导数作为判决函数，改善了比值法的效果，提高了信源估计的准确率。本专利技术利用直扩信号扩频码之间的不相关性，实现了在信源波达方向重合时的信源数估计，该方法同时适用于白噪声和有色噪声背景，在较低的信噪
比下也具有良好的估计性能，同时该方法在欠定(阵元数小于信源数)的情况下也能够稳定地实现信源数估计。
[0007]针对在同步直扩系统中，当信源之间的波达方向很近时造成信源数估计方法性能下降甚至失效的问题，本专利技术提出了针对同步直扩信号的信源数估计方法。
[0008]本专利技术解决其技术问题所采用的技术方案包括如下步骤：
[0009]步骤一：利用阵列得到接收信号X(t)，并将第r个阵元的接收信号单独提取出来，记为x
r
；
[0010]步骤二：将x
r
按照码元长度分段得到
[0011][0012]其中L是x
r
的长度，L0是码元长度，表示x
r
的第i段数据；
[0013]步骤三：将组成重构矩阵X
r
：
[0014][0015]其中p是分段数目，满足p＝[L/L0]；
[0016]步骤四：计算重构矩阵X
r
的协方差矩阵，得到重构协方差矩阵R
r
：
[0017]R
r
＝E(X
r
X
rH
)
ꢀꢀꢀ
(3)
[0018]其中H代表共轭转置；
[0019]步骤五：对重构协方差矩阵R
r
取平均值，得到平均协方差矩阵R
x
：
[0020][0021]其中M为阵元个数；
[0022]步骤六：对平均协方差矩阵R
x
进行特征分解：
[0023]R
x
＝U∑U
H
ꢀꢀꢀ
(5)
[0024]其中U是特征值对应的特征向量矩阵，∑是由特征值组成的对角阵；
[0025]步骤七：将特征值矩阵∑中的特征值λ按照从大到小排列，从大到小排列后计算特征值序列的比值导数，比值导数γ
i
按照式(7)计算：
[0026][0027]其中，i＝1,2,
…
,p
‑
2；
[0028]步骤八：计算比值导数γ
i
取最大时i的取值，该值等于信源数目的估计值k；
[0029][0030]最终实现信源波达方向相同时的同步直扩信号的信源数目估计。
[0031]所述步骤四中，重构协方差矩阵R
r
为：
[0032][0033]式中，R
s
为源信号的协方差矩阵，为重构噪声协方差矩阵，H代表共轭转置，I代
表单位矩阵，A表示分段信号的导向矩阵；其中：
[0034][0035][0036]式(4)S
p
表示经过截取后第p段源信号矩阵，fs代表采样率，式(11)中a
K
为第k个信源的信号幅度，θ
i
是第i个源信号的方向角，λ是源信号的波长，重构后的协方差矩阵中，第一项A是由子矩阵[a
1 a2ꢀ…ꢀ
a
K
]组成的块对角阵，利用初等列变换容易证明A的秩为K，第二项是S
p
的协方差矩阵R
s
，写成分块矩阵的形式：
[0037][0038]由于不同信号源的码元之间是不相关的，因此每一个子矩阵都是对角阵，而每一行的非零元素表示同一个源信号的不同码元之间的相关，非零元素的模相同，因此经过初等列变换只有第一列的子矩阵不等于单位阵，此时矩阵R
s
的秩等于子矩阵的秩等于K，得到rank(R
r
)≤K，实际上矩阵R
r
在无噪声的条件下有rank(R
r
)＝K，且R
r
为p阶Hermitian矩阵，对其进行特征分解能够得到p
‑
k重为0的特征值和k个较大的特征值。
[0039]本专利技术的有益效果在于本专利技术提出了一种同步直扩信号的信源数估计方法，通过重构观测矩阵，避免了导向矩阵秩亏，能够在所有信源本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种同步直扩信号的源数估计方法，其特征在于包括下述步骤：步骤一：利用阵列得到接收信号X(t)，并将第r个阵元的接收信号单独提取出来，记为x
r
；步骤二：将x
r
按照码元长度分段得到得到其中L是x
r
的长度，L0是码元长度，表示x
r
的第i段数据；步骤三：将组成重构矩阵X
r
：其中p是分段数目，满足p＝[L/L0]；步骤四：计算重构矩阵X
r
的协方差矩阵，得到重构协方差矩阵R
r
：R
r
＝E(X
r
X
rH
)
ꢀꢀꢀꢀ
(3)其中H代表共轭转置；步骤五：对重构协方差矩阵R
r
取平均值，得到平均协方差矩阵R
x
：其中M为阵元个数；步骤六：对平均协方差矩阵R
x
进行特征分解：R
x
＝U∑U
H
ꢀꢀꢀꢀ
(5)其中U是特征值对应的特征向量矩阵，∑是由特征值组成的对角阵；步骤七：将特征值矩阵∑中的特征值λ按照从大到小排列，从大到小排列后计算特征值序列的比值导数，比值导数γ
i
按照式(7)计算：其中，i＝1,2,
…
,p
‑
2；步骤八：计算比值导数γ
i
取最大时i的取值，该值等于信源数目的估计值k；最...

【专利技术属性】
技术研发人员：申晓红，王奕成，王海燕，姜喆，马高悦，赵重达，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：