基于混合聚类算法的直扩信号伪码周期及序列估计方法与流程

1.本发明涉及通信技术领域，具体涉及一种基于混合聚类算法的直扩信号伪码周期及序列估计方法。


背景技术：

2.针对直扩信号伪码周期及序列估计问题，现有的矩阵分解算法，及一些搜索迭代算法，算法运算量都比较大，估计速度慢，且对噪声敏感，对低信噪比下的直扩信号伪码周期及序列无法准确、快速地获取，难以满足应用需求。


技术实现要素：

3.本发明所要解决的技术问题在于避免上述背景技术中的不足之处而提供一种新的易于实现的直扩信号伪码周期及序列估计方法。
4.本发明所采取的技术方案是：
5.一种基于混合聚类算法的直扩信号伪码周期及序列估计方法，包括以下步骤：
6.(1)计算信号二次功率谱，搜索谱峰的位置，用匹配法统计出谱峰的间距，即得到伪码周期；
7.(2)将信号按照伪码周期进行数组排列，通过混合聚类算法对数组进行聚类处理，估计出信息序列；
8.(3)将估计得到的信息序列进行积累，并去均值，然后判决得到伪码序列。
9.进一步的，步骤(1)具体包括以下步骤：
10.(101)计算信号二次功率谱，即先计算信号的功率谱，然后将信号的功率谱再次计算一次功率谱，即得到信号的二次功率谱，并搜索谱峰位置，将相邻位置两两相减，计算出谱峰的间距值dis＝{dis1,dis2,
…
,disn}，n为谱峰间距值的个数；
11.(102)从谱峰的间距值中，采用匹配法统计出正确的间距；具体为：根据码片速率rs计算过采样率pps，pps＝fs/rs，fs为信号采样率，并根据已知的k种常见伪码周期p＝{p1,p2,
…
,pi,
…
pk}，计算pps*p，得到伪码采样周期可能存在的初始值d＝{d1,d2,
…
di,
…
dk}，di＝pps*pi；
12.(103)计算谱峰间距dis的所有元素与各初始值d的相似度，相似度为计算每个dis值与d的余数，即mod(disj,di),j＝1,2,
…
,n,i＝1,2,
…
,k，将disj归于余数最小值对应的di的集合，完成计算后，统计每个集合所包含的元素个数，元素最多且大于一定数量的集合所对应的初始值di，即为伪码采样周期，pi则为对应伪码周期估计值，用pn表示，pn＝pi。
13.进一步的，步骤(2)具体包括以下步骤：
14.(201)找到信号最佳起始采样点，按过采样率pps进行抽取，得到同步后的信号，并按照伪码周期进行数组排列，得到数组x＝(x1,x2,
…
xi,
…
xn)
t
，其中xi＝(x
i1
,x
i2
,
…
,x
ipn
)，n为信息码片数量，pn为伪码周期，每一行为一个信息码片，包含一个完整周期的伪码序列；
15.(202)对层次聚类法及kmeans算法进行改进，结合信息序列的特点，形成新的混合
聚类算法，对数组进行聚类处理，具体方法为：
16.首先伪码序列为 1和-1，因此数组x中的对象通过聚类分为两类，定义为c1和c2；计算数组x中每两个对象的欧氏距离，即m(i,j)＝dist(xi,xj)，i,j＝1,2,
…
,n；选取距离值最大的两个对象，分别作为c1的聚类中心u1和c2的聚类中心u2；
17.计算x中所有对象与u1和u2的欧氏距离dist(xi，uk)，i＝1，2，
…
,n,k＝1,2；若dist(xi,u1)小于dist(xi,u2)，则xi属于c1，否则，xi属于c2；
18.完成聚类后，将c1和c2对应的样本分别赋值为 1和-1，即对数组x＝(x1,x2,
…
xi,
…
xn)
t
进行处理，xi若属于c1，则赋值为 1，若属于c2，则赋值为-1，由此得到一个由 1和-1组成的长度为n的序列，即为信息序列。
19.进一步的，步骤(3)具体包括以下步骤：
20.(301)对已完成聚类处理的信息码片向量作同相相加，反相相减的累加处理；具体为：将所有属于c1的xi相加，得到s1，将所有属于c2的xi相加，得到s2，再用s1-s2，得到向量m，m＝(m1,m2,
…mpn
)，即完成累加处理；
21.(302)对序列m去均值处理，再进行判决，大于0的为 1，小于0的为-1，即得到伪码序列。
22.本发明与现有技术相比具有如下优点：
23.1、本发明的伪码周期估计算法，对伪码周期进行匹配并统计，相比于其它搜索算法，算法运算量小，能在低信噪比下准确、快速地完成伪码周期估计。
24.2、本发明采用的混合聚类算法进行码元信息序列估计，避免了层次聚类法合并或分裂点的选择及反复的迭代过程，也避免了kmeans算法由于初始中心的选择错误，导致会形成局部最优解或迭代运算时间过长。该算法对象的分组更加合理，能得到更加稳定的聚类结果，且降低了计算复杂度，节约了聚类运算的时间，提高了算法运算速度，更利于工程实时检测估计。
附图说明
25.图1是本发明的工作过程示意图。
具体实施方式
26.参照图1，本发明基于混合聚类算法的直扩信号伪码周期及序列估计方法，工作过程包括：匹配法估计伪码周期1、混合聚类算法估计信息序列2、估计伪码序列3。具体包括以下步骤：
27.(1)计算信号二次功率谱，搜索谱峰的位置，用匹配法统计出谱峰的间距，即得到伪码周期；
28.具体包括以下步骤：
29.(101)计算信号二次功率谱，即先计算信号的功率谱，然后将信号的功率谱再次计算一次功率谱，即得到信号的二次功率谱，并搜索谱峰位置，将相邻位置两两相减，计算出谱峰的间距值dis＝{dis1,dis2,
…
,disn}，n为谱峰间距值的个数；
30.(102)谱峰的间距，即为伪码周期。在信噪比低时，部分谱峰已经被噪声淹没，导致检测失败，也可能检测到虚假的谱峰，因此估算的谱峰之间的距离也存在虚警，对应估计出
的伪码周期的准确率也随之下降。本发明采用匹配法从谱峰的间距值中，统计出正确的间距；具体为：根据码片速率rs计算过采样率pps，pps＝fs/rs，fs为信号采样率，并根据已知的k种常见伪码周期p＝{p1,p2,
…
,pi,
…
pk}，计算pps*p，得到伪码采样周期可能存在的初始值d＝{d1,d2,
…
di,
…
dk}，di＝pps*pi；
31.(103)计算谱峰间距dis的所有元素与各初始值d的相似度，相似度为计算每个dis值与d的余数，即mod(disj,di),j＝1,2,
…
,n,i＝1,2,
…
,k，将disj归于余数最小值对应的di的集合，完成计算后，统计每个集合所包含的元素个数，元素最多且大于一定数量的集合所对应的初始值di，即为伪码采样周期，pi则为对应伪码周期估计值，用pn表示，pn＝pi。
32.(2)得到伪码周期后，将信号按照伪码周期进行数组排列，通过混合聚类算法对数组进行聚类处理，估计出信息序列；
33.步骤(2)具体包括以下步骤：
34.(201)将同步后的信号，即找到信号最佳起始采样点，按过采样率pps进行抽取，得到同步后的信号，并按照伪码周期进行数组排列，得到数组x＝(x1,x2,
…
xi,
…
xn)
t
，其中xi＝(x
i1
,x
i2
,
…
,x
ipn
)，n为信息码片数量，pn为伪码周期，每一行为一个信息码片，包含一个完整周期的伪码序列；
35.(202)对层次聚类法及kmeans算法进行改进，结合信息序列的特点，形成新的混合聚类算法，对数组进行聚类处理，具体方法为：
36.首先伪码序列为 1和-1，因此数组x中的对象通过聚类分为两类，定义为c1和c2；计算数组x中每两个对象的欧氏距离，即m(i,j)＝dist(xi,xj)，i,j＝1,2,
…
,n；选取距离值最大的两个对象，分别作为c1的聚类中心u1和c2的聚类中心u2；
37.计算x中所有对象与u1和u2的欧氏距离dist(xi，uk)，i＝1，2，
…
,n,k＝1,2；若dist(xi,u1)小于dist(xi,u2)，则xi属于c1，否则，xi属于c2；
38.完成聚类后，将c1和c2对应的样本分别赋值为 1和-1，即对数组x＝(x1,x2,
…
xi,
…
xn)
t
进行处理，xi若属于c1，则赋值为 1，若属于c2，则赋值为-1，由此得到一个由 1和-1组成的长度为n的序列，即为信息序列。
39.(3)将估计得到的信息序列进行积累，并去均值，然后判决得到伪码序列。
40.具体包括以下步骤：
41.(301)对已完成聚类处理的信息码片向量作同相相加，反相相减的累加处理；具体为：将所有属于c1的xi相加，得到s1，将所有属于c2的xi相加，得到s2，再用s1-s2，得到向量m，m＝(m1,m2,
…mpn
)，即完成累加处理；
42.(302)对序列m去均值处理，再进行判决，大于0的为 1，小于0的为-1，即得到伪码序列。