本发明属于信号解调领域，尤其涉及一种基于数据库匹配的水声直扩信号m序列盲估计方法。
背景技术：
扩频通信能在负信噪比条件下正常通信,给侦察和解调带来了很大的难度。为了获取ds信号的更多信息,获取伪码序列成为研究的热点。最常用的是相关法,将噪声建模为白噪声,但当存在相关噪声时,相关法的估计性能将严重下降,对伪码序列的获取也无能为力。相关矩阵分解法,能获取基带ds信号的伪码序列,但没有用到伪码序列的相关性；基于神经网络的多主分量分析法,能提高运算速度,也没有用到伪码序列的相关性,性能很难进一步提高；基于massey算法的线性移位反馈虽然用到了伪码的相关性,但巨大的运算量使之仅停留在理论研究阶段。
传统的m序列估计方法是得到截获信号的三阶相关函数及其相关峰坐标后，利用矩阵斜消变换求出各个相关峰坐标之间的最大公因式，从而得到m序列的估计结果。但由于水声信道中严重的噪声干扰和多径效应，三阶相关函数相关峰的坐标会发生偏移，使得估计结果出现较大误差。本发明从三阶相关函数的相关峰坐标和m序列本原多项式的一一对应关系出发，在无噪声的高斯信道条件下求出不同m序列的本原多项式所对应的相关峰坐标，按阶数的不同将m序列与其对应的相关峰坐标建成数据库。这种方法在保证实时性的前提下可以更好的适应水声信道环境。
中国专利cn105680903a公开了一种周期长短码直扩码分多址信号多伪码估计方法，该方法充分利用矩阵填充的数学模型和m序列三阶相关峰特性，依次实现信号各用户复合码序列、长扰码序列和短扩频码序列的盲估计。
中国专利cn106817326a中公开了一种多用户周期长短码直扩信号的伪码盲估计方法，该方法充分利用fast-ica算法以及m序列三阶相关特性，实现对多用户周期长短码直扩信号的两个伪随机码的盲估计。
中国专利cn105680904a中公开了一种非周期长码直扩信号伪码估计方法，该方法通过循环去均值处理和拟合优度检验更为精确地提取三阶相关峰，实现对非周期长码直扩信号的伪码盲估计。
但以上三种方法在利用m序列三阶相关峰特性时均没有建立完备的数据库，并不是通过数据库匹配的方式得到m序列的估计结果。以上方法在噪声严重、多径时延较大的水声信道中性能恶化严重，并不能很好的适应水声环境的实际需求。
技术实现要素：
本发明的目的在于公开可靠性高，实时性强的一种基于数据库匹配的水声直扩信号m序列盲估计方法。
本发明的目的是这样实现的：
一种基于数据库匹配的水声直扩信号m序列盲估计方法，包含如下步骤：
步骤(1)：求出在某一阶数r下所有m序列的本原多项式：
某一阶数r下的m序列的本原多项式f(x)，满足：
(1)f(x)是既约的，即不能再分解因式；
(2)f(x)可整除xm+1，其中m＝2r-1；
(3)f(x)不能整除xq+1，其中q＜r。
对于某一阶数r，通过遍历的方法得到某一阶数r下所有m序列的本原多项式。
步骤(2)：依次求出每一个m序列的三阶相关函数并记录下对应的所有相关峰坐标：
根据三阶相关函数的定义，m序列的三阶相关函数在其值域内可以表示为：
每一行每一列都有且仅有一个相关峰坐标存在，且这些相关峰坐标与m序列一一对应。记录所有相关峰坐标及其对应的m序列。
步骤(3)：对于一个m序列的三阶相关函数，将相关峰坐标及其对应的m序列一起作为样本存放在数据库，以便完成坐标匹配后直接输出m序列的估计结果。
步骤(4)：改变阶数r的值，重复步骤(1)至步骤(3)，完成所需数据库的建立：
根据实际使用需求，确定阶数r的取值范围。分别对每一阶数r建立完备的数据库：
{(p',q')i}r，i＝1,2,…,n，
确保截获的dsss信号所使用的m序列存在于数据库中。
步骤(5)：求出截获dsss信号的三阶相关函数和所有的相关峰坐标{(ps,qs)}r；
步骤(6)：将截获dsss信号的相关峰坐标与数据库中的样本进行匹配；
匹配程度δ：
上式中，ni是{(ps,qs)}r和{(p',q')i}r两个集合中相同元素的个数，nr是{(p',q')i}r中所有元素的个数。n是对于某一阶数r的本原多项式的最大个数。
分别计算{(ps,qs)}r与相同阶数r的数据库中的样本之间的匹配程度δ。
步骤(7)：将匹配程度δ最高的数据库中的样本作为m序列的估计结果：当δ最大时，对应的数据库中的样本所代表的m序列即为截获dsss信号所使用的m序列的估计结果。
本发明的有益效果为：
本发明简单易行，可靠性高，既能保证非合作通信的实时性，又能很好地适应水下声信道噪声严重，多径时延大的特点。
附图说明
图1为数据库匹配三阶相关函数的m序列盲估计原理图；
图2为数据库结构图。
具体实施方式
下面结合附图来进一步描述本发明：
如图1，一种基于数据库匹配的水声直扩信号m序列盲估计方法，包含如下步骤：
步骤(1)：求出在某一阶数r下所有m序列的本原多项式：
某一阶数r下的m序列的本原多项式f(x)，满足：
(1)f(x)是既约的，即不能再分解因式；
(2)f(x)可整除xm+1，其中m＝2r-1；
(3)f(x)不能整除xq+1，其中q＜r。
对于某一阶数r，通过遍历的方法得到某一阶数r下所有m序列的本原多项式。
步骤(2)：依次求出每一个m序列的三阶相关函数并记录下对应的所有相关峰坐标：
根据三阶相关函数的定义，m序列的三阶相关函数在其值域内可以表示为：
每一行每一列都有且仅有一个相关峰坐标存在，且这些相关峰坐标与m序列一一对应。记录所有相关峰坐标及其对应的m序列。
步骤(3)：对于一个m序列的三阶相关函数，将相关峰坐标及其对应的m序列一起作为样本存放在数据库，以便完成坐标匹配后直接输出m序列的估计结果。
步骤(4)：改变阶数r的值，重复步骤(1)至步骤(3)，完成所需数据库的建立：
如图2，根据实际使用需求，确定阶数r的取值范围。分别对每一阶数r建立完备的数据库：
{(p',q')i}r，i＝1,2,…,n，
确保截获的dsss信号所使用的m序列存在于数据库中。
步骤(5)：求出截获dsss信号的三阶相关函数和所有的相关峰坐标{(ps,qs)}r；
步骤(6)：将截获dsss信号的相关峰坐标与数据库中的样本进行匹配；
匹配程度δ：
上式中，ni是{(ps,qs)}r和{(p',q')i}r两个集合中相同元素的个数，nr是{(p',q')i}r中所有元素的个数。n是对于某一阶数r的本原多项式的最大个数。
分别计算{(ps,qs)}r与相同阶数r的数据库中的样本之间的匹配程度δ。
步骤(7)：将匹配程度δ最高的数据库中的样本作为m序列的估计结果：当δ最大时，对应的数据库中的样本所代表的m序列即为截获dsss信号所使用的m序列的估计结果。
以上所述并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。