一种直扩信号码片宽度估计方法技术

本发明专利技术公开了一种直扩信号码片宽度估计方法，利用直扩信号伪码游程分布特性，即长度为1的游程出现概率最高，而长度为1的伪码对应直扩信号的一个码片，通过对出现概率最大的基带信号过零点间隔搜索获得码片宽度的估计。本发明专利技术提出的方法提高了码片宽度估计精度，适应低信噪比条件且运算量较小，具有较高的实际应用价值。

【技术实现步骤摘要】

本专利技术属于通信和信号处理
，涉及一种直扩信号码片宽度估计方法，用于非合作条件下直扩信号码片宽度的估计。
技术介绍
直接序列扩频(以下简称直扩)技术在水声通信与探测领域已得到广泛应用，其特点是通过伪随机序列将信号频带展宽进行传输，接收端通过解扩处理获得扩频增益。而对于非合作方，直扩信号单位频点信噪比低，其截获检测与参数估计技术是通信与水声对抗中的研究热点与难点。直扩信号码片宽度估计是非合作截获处理的重要部分，并且是直扩信号参数估计的难点。现有的直扩信号码片宽度估计方法主要特征有：1基于功率谱零点的方法利用直扩信号功率谱左右零点间隔与码片宽度倒数的对应关系进行估计，存在功率谱零点位置估计精度较低、低信噪比下性能严重下降的问题；2基于周期谱估计的方法需要较大的运算量，并且通过多次累积来提高估计精度和适应低信噪比。
技术实现思路
技术问题：本专利技术提供一种估计结果精度高、适应低信噪比条件且运算量较小的直扩信号码片宽度估计方法。技术方案：本专利技术的直扩信号码片宽度估计方法，包括以下步骤：步骤1：获取离散数据序列x(n),n＝0,1,2,…,N-1：从传感器接收N个采样点的实时采集数据或提取N个采样点存储的现成数据作为待处理的数据序列x(n),n＝0,1,2,…,N-1，其中N为采样点数；步骤2：估计所述离散数据序列x(n),n＝0,1,2,…,N-1的载波频率fc；步骤3：提取基带信号y(n)：利用所述步骤2估计出的载波频率fc，对数据序列x(n)进行下移频处理，即将数据序列x(n)搬移到基带，得到基带信号y(n)；步骤4：统计所述基带信号y(n)过零点之间的间隔J(k),k＝1,2,…,K-1，K为基带信号y(n)的过零点数；步骤5：对所述过零点间隔J(k)从小到大进行排序，得到Js(k),k＝1,2,…,K；步骤6：设定过零点间隔允许偏差ε，将Js(k)分为L段：划分依据为每段中最后一个过零间隔与该段开始过零间隔之间的相对偏差小于ε，其中ε＜0.5；步骤7：计算Js(k)分段后的每段中包含的过零点间隔的个数N(l),l＝1,2,…,L，以及均值m(l),l＝1,2,…,L；步骤8：找出N(l),l＝1,2,…,L中的最大值Nmax(lmax)，则第lmax段中过零点间隔出现概率最大；步骤9：根据下式计算码片宽度其中m(lmax)为Js(k)第lmax段的均值，fs为采样频率。进一步的，本专利技术方法中，步骤1中，采样点数N覆盖多个码元，即包括多个码片宽度。本专利技术方法利用伪码的游程分布特性，通过基带信号最大概率跳变间隔的搜索，获得码片宽度的估计结果，并且具有估计精度高、适应低信噪比条件且运算量较小的特性。有益效果：本专利技术与现有技术相比，具有以下优点：1.利用直扩信号伪码游程分布特性，即长度为1的游程出现概率最高，长度为1的伪码对应直扩信号的一个码片，通过对出现概率最大的过零点间隔搜索即可得到码片宽度的估计。本专利技术方法原理简单，并且通过多个对应码片宽度的过零点间隔进行平均得到码片宽度估计结果，提高了估计精度；2.相较于功率谱零点估计方法，本专利技术方法具有更高的估计精度和更低的信噪比下限；3.相较于基于周期谱估计的方法，本专利技术方法计算量小，能够实时处理大量数据，实用性强。现有伪码宽度估计技术常采用功率谱零点法和周期谱估计方法。功率谱零点法通过估计直扩信号功率谱零点位置获得信号带宽的估计值，根据带宽与码片宽度的关系估计出码片宽度。功率谱零点法算法简单，但估计精度较低。本专利技术方法通过多个码片宽度估计的平均，提高了估计精度。周期谱估计方法需要估计信号的周期谱，通过对幅值谱的谱峰搜索得到码片宽度的估计。周期谱方法估计精度高，但原理复杂，数字实现时计算复杂度高、计算速度慢，实时系统应用难度大。本专利技术方法保证高估计精度的同时，算法简单，计算量小，实时处理能力强。附图说明图1所示为本专利技术的流程图。图2所示为实施例1中无噪声时过零点间隔统计结果。图3所示为实施例1中-5dB信噪比下过零点间隔统计结果。图4所示为实施例1中不同信噪比条件下本专利技术方法与接收信号零点法码片宽度估计均方根误差统计结果。具体实施方式下面结合实施例和说明书附图对本专利技术作进一步的说明。本说明书及其附图中给出了本专利技术的较佳的实施例，但本专利技术可以以许多不同的形式来实现，并不限于本说明书所描述的实施例。相反地，提供这些实施例的目的是使对本专利技术的公开内容的理解更为透彻全面。步骤1：获取离散数据序列x(n),n＝0,1,2,…,N-1：从传感器接收N个采样点的实时采集数据或提取存储的现成数据作为待处理的数据序列x(n),n＝0,1,2,…,N-1，其中N为采样点数。在步骤1中，采样数据可以通过数据采集器采集。步骤2：估计所述数据序列x(n),n＝0,1,2,…,N-1的载波频率fc。在步骤2中，载波频率的估计属于现有技术，例如《一种水声二相调制直扩信号载频估计方法》专利技术专利。步骤3：提取基带信号y(n)：利用所述步骤2估计出的载波频率fc，对数据序列x(n)进行下移频处理，即将数据序列x(n)搬移到基带，得到基带信号y(n)。步骤4：统计所述基带信号y(n)过零点之间的间隔J(k),k＝1,2,…,K-1，K为基带信号y(n)的过零点数；。步骤5：对所述过零点间隔J(k)进行排序，得到Js(k),k＝1,2,…,K。在步骤5中排序可以按照从小到大排列，也可以按照从大到小排列。步骤6：设定过零点间隔允许偏差ε，将Js(k)分为L段：划分依据为每段中最后一个过零间隔与该段开始过零间隔之间的相对偏差小于ε，其中ε＜0.5。在步骤6中，ε的设定需要小于0.5，否则会将两个游程对应的过零点间隔划分为一组，也不能太小，以避免将相同游程对应的过零点间隔分为两组，一般可设定为0.2。步骤7：计算Js(k)分段后的每段中包含的过零点间隔的个数N(l),l＝1,2,…,L，以及均值m(l),l＝1,2,…,L；步骤8：找出N(l),l＝1,2,…,L中的最大值Nmax(lmax)，则第lmax段中过零点间隔出现概率最大；步骤9：根据下式计算码片宽度其中m(lmax)为Js(k)第lmax段的均值，fs为采样频率。本专利技术的原理是利用直扩信号伪码的游程分布特性，即长度为1的游程出现概率最高，而长度为1的伪码对应直扩信号的一个码片宽度，反应在基带信号上则为出现概率最大的过零点间隔，因此通过对出现概率最大的过零点间隔搜索和过零点间隔平均则可得到码片宽度的估计。本专利技术原理简单，并且充分利用了多个过零点间隔的平均，可以有效提高估计精度。实施例1仿真信号参数为：采样频率fs＝50KHz，载波频率f0＝5.0kHz，伪码阶数为5，采用二相调制方式，一个码片内的载波周期数为8，则码片宽度为1.6ms，离散采样点数为80点，信号幅度A0＝1，叠加零均值高斯白噪声，方差σ2的大小由信噪比SNR决定：图2和图3分别为无噪声和信噪比为-5dB两种情况下基带信号过零点间隔从小到大排序后的结果。从图2可以看出过零点间隔从小到大排序后，过零点间隔出现概率服从伪码游程分布特性，采样长度为80的间隔出现概率最大，根据最大概率获得离散采样长度为80，除以采样频率后可得到1.6ms的码片宽度估计结果。存在噪声后，过零点间本文档来自技高网...

【技术保护点】
一种直扩信号码片宽度估计方法，其特征在于，该方法包括以下步骤：步骤1：获取离散数据序列x(n),n＝0,1,2,…,N‑1：从传感器接收N个采样点的实时采集数据或提取N个采样点存储的现成数据作为待处理的数据序列x(n),n＝0,1,2,…,N‑1，其中N为采样点数；步骤2：估计所述离散数据序列x(n),n＝0,1,2,…,N‑1的载波频率fc；步骤3：提取基带信号y(n)：利用所述步骤2估计出的载波频率fc，对数据序列x(n)进行下移频处理，即将数据序列x(n)搬移到基带，得到基带信号y(n)；步骤4：统计所述基带信号y(n)过零点之间的间隔J(k),k＝1,2,…,K‑1，K为基带信号y(n)的过零点数；步骤5：对所述过零点间隔J(k)从小到大进行排序，得到Js(k),k＝1,2,…,K；步骤6：设定过零点间隔允许偏差ε，将Js(k)分为L段：划分依据为每段中最后一个过零间隔与该段开始过零间隔之间的相对偏差小于ε，其中ε＜0.5；步骤7：计算Js(k)分段后的每段中包含的过零点间隔的个数N(l),l＝1,2,…,L，以及均值m(l),l＝1,2,…,L；步骤8：找出N(l),l＝1,2,…,L中的最大值Nmax(lmax)，则第lmax段中过零点间隔出现概率最大；步骤9：根据下式计算码片宽度其中m(lmax)为Js(k)第lmax段的均值，fs为采样频率。...

【技术特征摘要】
1.一种直扩信号码片宽度估计方法，其特征在于，该方法包括以下步骤：步骤1：获取离散数据序列x(n),n＝0,1,2,…,N-1：从传感器接收N个采样点的实时采集数据或提取N个采样点存储的现成数据作为待处理的数据序列x(n),n＝0,1,2,…,N-1，其中N为采样点数；步骤2：估计所述离散数据序列x(n),n＝0,1,2,…,N-1的载波频率fc；步骤3：提取基带信号y(n)：利用所述步骤2估计出的载波频率fc，对数据序列x(n)进行下移频处理，即将数据序列x(n)搬移到基带，得到基带信号y(n)；步骤4：统计所述基带信号y(n)过零点之间的间隔J(k),k＝1,2,…,K-1，K为基带信号y(n)的过零点数；步骤5：对所述过零点间隔J(k)...

【专利技术属性】
技术研发人员：王晓燕，方世良，李帝水，姚帅，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32