一种基于Chirp相关峰值位置偏差趋势估计采样偏差的方法与流程

一种基于chirp相关峰值位置偏差趋势估计采样偏差的方法
技术领域
1.本发明涉及无线扩频通信系统的同步技术领域，尤其设计一种基于chirp相关峰值位置偏差趋势估计采样偏差的方法。


背景技术：

2.扩频码的同步是指到达接收机的扩频编码信号与本地参考扩频信号在码的图案位置和码时钟速率在时间上都是准确一致的，如果不一致就有了码元同步的偏移。扩频码未完全同步时，扩频信号无法解扩，将导致信息传输的失败。
3.在目前的无线通信接收机中，中频数字化的软件无线电系统结构如图1所示，结构中的本地振荡器提供的频率源，使ad采样器对接收到的模拟中频信号进行固定频率的采样将其转换成数字信号。但由于发射机和接收机的采样时钟是不同晶振产生的，因此，接收机对接收到的信号进行采样时，会产生采样偏差(也称定时偏差)，使得接收机采样得到的信号不是采样在最佳采样点上，如图2所示，从而产生扩频码的解扩失败。
4.以上背景技术内容的公开仅用于辅助理解本发明的构思及技术方案，其并不必然属于本专利申请的现有技术，在没有明确的证据表明上述内容在本专利申请的申请日已经公开的情况下，上述背景技术不应当用于评价本技术的新颖性和创造性。


技术实现要素：

5.为解决上述技术问题，本发明提出一种基于chirp相关峰值位置偏差趋势估计采样偏差的方法，充分利用chirp信号解扩信息，利用较少的资源开销，联合实现了ppm的估计。
6.为了达到上述目的，本发明采用以下技术方案：
7.本发明公开了一种基于chirp相关峰值位置偏差趋势估计采样偏差的方法，包括：采用接收数据与本地基础图样的时域圆周相关来对接收的chirp信号进行解扩，以查找解扩的相关峰值位置；利用解扩的相关峰值位置，通过迭代方式对峰值位置的移动趋势进行分析以实现对采样偏差的估计。
8.优选地，其中的chirp信号是通过对基础图样的循环移位来表示不同的扩频符号。
9.优选地，其中解扩的相关峰值位置与扩频符号存在映射关系。
10.优选地，在通过迭代方式对峰值位置的移动趋势进行分析以实现对采样偏差的估计之前还包括：对连续接收的多个相关峰值位置求平均，以估计得到有效参考峰值位置，并根据有效参考峰值位置，来分析峰值位置的移动趋势。
11.优选地，通过迭代方式对相关峰值位置的移动趋势进行分析以实现对采样偏差的估计具体包括：设定偏差阈值，对相关峰值位置的移动趋势的合理性进行筛选，通过迭代方式进行趋势累加得到累计偏移取值，再将累计偏移取值转化为采样偏差。
12.优选地，将累计偏移取值转化为采样偏差具体包括：将累计偏移取值转换为ppm，并结合chirp信号的扩频因子sf和采样率，转换成采样点补偿值。
13.优选地，采用接收数据与本地基础图样的时域圆周相关来对接收的chirp信号进行解扩，以查找解扩的相关峰值位置具体包括：将接收的chirp信号rx与本地的基础图样local利用fft和ifft实现时域的圆周相关值pv，查找主峰值pvmax及对应的相关峰值位置pdmax，具体实现公式为：pv＝ifft(fft(rx).*conj(local))，pvmax＝max(abs(pv).^2)＝abs(pv[pdmax].^2)。
[0014]
优选地，其中相关峰值位置的计算包含整数部分和小数部分，整数部分为直接的峰值索引pdint，小数部分pdfrac为基于主峰值和左右次峰值的二分插值所得：
[0015][0016]
其中，pvmaxedge0、pvmaxedge1分别为左右次峰值；
[0017]
并更新相关峰值位置：pdmax＝pdint+pdfrac。
[0018]
优选地，利用解扩的相关峰值位置，通过迭代方式对峰值位置的移动趋势进行分析以实现对采样偏差的估计具体包括：
[0019]
先对连续接收的m个峰值满足
‑
pdwin<pdmax[k]<pdwin，k＝1,
…
m的相关峰值位置进行存储，对存储的m个峰值求平均得到初始的有效峰值位置real_calc，首次成功符号数值下标：real_start＝m/2+1；
[0020]
将当前峰值设置为当前的有效峰值calc＝pdmax，符号索引计数idx＝cntnum，计算位置累计偏移sum_trendvalue＝calc
‑
real_calc，符号索引sum_trendidx＝idx
‑
real_start；
[0021]
在确定有效起始位置后，通过迭代方式对后续的连续峰值位置信息pdmax，进行偏差合理性判定：abs(pdmax
‑
calc)<(cntnum
‑
idx)*deltasthr，其中deltasthr为连续符号的偏移判决门限，对于满足偏差策略的峰值，累计入位置偏移sum_trendvalue及符号索引sum_trendidx，并更新当前峰值为calc，其对应的符号索引为idx，作为下次迭代偏差合理性判决的参考信息；
[0022]
将累计位置偏移sum_trendvalue转换为ppm，并结合chirp的扩频因子sf和采样率，转换为采样点补偿值。
[0023]
与现有技术相比，本发明的有益效果在于：本发明提出的基于chirp相关峰值位置偏差趋势估计采样偏差的方法，对所有接收到的chirp信号的解扩，都与本地基础图样的时域圆周相关实现，从而简化了本地的序列存储，有效的降低了解扩的运算量和时间；进一步充分利用了chirp信号的解扩峰值位置信息，使系统将解扩和ppm的估计有效结合，从而减少了计算资源。
附图说明
[0024]
图1是软件无线电的接收机硬件框图；
[0025]
图2是采样时钟偏差示意图；
[0026]
图3是本发明优选实施例的基于特殊chirp相关峰值位置偏差趋势估计采样偏差的方法流程图；
[0027]
图4是4进制各符号时频图；
[0028]
图5是本发明具体实施例的采样偏差工作流程图。
具体实施方式
[0029]
下面对照附图并结合优选的实施方式对本发明作进一步说明。
[0030]
本发明优选实施例公开了一种基于特殊chirp相关峰值位置偏差趋势估计采样偏差的方法，该方法包括以下步骤：(1)该特殊chirp信号是通过对基础图样的循环移位来表示不同的扩频符号；(2)对chirp信号的解扩采用接收数据与本地基础图样的时域圆周相关实现，其解扩的相关峰值位置与扩频符号存在映射关系；(3)该采样偏差的估计方案充分利用解扩信息
‑
相关峰值位置，通过迭代的对峰值位置的移动趋势进行分析实现采样偏差的估计；(4)该采样偏差估计对峰值位置的趋势判定分三步走：有效参考峰值位置的估计
‑
>迭代的进行峰值位置的合理性判定并进行趋势的累加存储
‑
>将累加趋势转换为ppm，并结合chirp的扩频因子sf和采样率，转换为采样点补偿值。本发明充分利用chirp信号的解扩峰值位置信息，使系统将解扩和ppm的估计有效结合，从而减少了计算资源。
[0031]
在本发明优选实施例的无线扩频通信系统中，构造了一种特殊的chirp信号，以某一个符号的chirp信号作为基础图样，通过对基础图样的循环移位实现对不同符号的扩频映射。在rx端，对所有接收到的chirp信号的解扩，都与本地基础图样的时域圆周相关实现，因此简化了本地的序列存储，有效的降低了解扩的运算量和时间，且其解扩的相关峰值位置与扩频符号存在映射关系。具体地，本发明优选实施例公开的基于特殊chirp相关峰值位置偏差趋势估计采样偏差的方法具有以下特征：
[0032]
1)构造了一种特殊的chirp信号，其是对基础图样的循环移位产生，其解扩简化了收端的本地序列的存储；
[0033]
2)采样偏差估计是基于对特殊chirp信号的解扩峰值位置进行趋势分析实现采样偏差的估计；
[0034]
3)该估计方法需要预设偏差阈值，进行峰值趋势合理性筛选，迭代的进行趋势累加，最终将偏移取值转化为ppm偏差。
[0035]
其中，对于不同符号的chirp信号，都是通过对基础图样的循环移位后的chirp信号表示；在rx端，对chirp信号的解扩采用接收数据与本地基础图样的时域圆周相关实现，其解扩的相关峰值位置与扩频符号存在映射关系。
[0036]
将接收chirp信号rx与本地的基础图样local利用fft和ifft实现时域的圆周相关值pv，查找主峰值pvmax及对应的峰值位置pdmax，具体实现见下式(1)、(2)：
[0037]
pv＝ifft(fft(rx).*conj(local))
ꢀꢀ
(1)
[0038]
pvmax＝max(abs(pv).^2)＝abs(pv[pdmax].^2)
ꢀꢀ
(2)
[0039]
其中，峰值位置的计算包含整数部分和小数部分，整数部分为直接的峰值索引pdint，小数部分pdfrac为基于主峰和左右次峰(pvmaxedge0、pvmaxedge1)的二分插值所得，具体见公式(3)。此处理实现减少chirp符号的个数的同时，提升估计精度；
[0040]
pvmaxedge0＝pv[pdint
‑
1]，pvmaxedge1＝pv[pdint+1]
[0041][0042]
更新峰值位置：pdmax＝pdint+pdfrac。
[0043]
通过求平均确定参考有效初始峰值位置，即在启动ppm的估计算法后，先对连续接收的m个峰值满足
‑
pdwin<pdmax[k]<pdwin(k＝1,
…
m)的相关峰值进行存储。对存储的m个峰值求平均得到初始的有效峰值位置real_calc，首次成功符号数值下标：real_start＝m/2+1。
[0044]
将当前峰值设置为当前的有效峰值calc＝pdmax，符号索引计数idx＝cntnum，计算位置累计偏移sum_trendvalue＝calc
‑
real_calc，符号索引sum_trendidx＝idx
‑
real_start。
[0045]
在确定有效起始位置后，迭代的对后续的连续峰值位置信息pdmax，进行偏差合理性判定：abs(pdmax
‑
calc)<(cntnum
‑
idx)*deltasthr。其中deltasthr为连续符号的偏移判决门限，对于满足偏差策略的峰值，累计入位置偏移sum_trendvalue及符号索引sum_trendidx，并更新当前峰值为calc，其对应的符号索引为idx，作为下次迭代偏差合理性判决的参考信息。
[0046]
将累计位置偏移sum_trendvalue转换为ppm，并结合chirp的扩频因子sf和采样率，转换为采样点补偿值。
[0047]
下述以具体实施例对本发明优选实施例公开的基于特殊chirp相关峰值位置偏差趋势估计采样偏差的方法作进一步的说明。
[0048]
在本具体实施例中，先设定调频信号的基础图样由up
‑
chirp和down
‑
chirp组合而成，如公式(4)所示。这里设定不同符号的扩频图样为基础图样(符号0)的向右移位循环。当然本发明的基础图样不限于此图样。
[0049][0050]
其中，a为信号的幅度，μ为线性扩频斜率，分别为第一段和第二段的起始相位，t
c
为一个符号周期；以m＝4进制为例，各个符号的时频图见图4：图4中，将符号长度t
c
等分为4段，每段长t
step
，即t
step
为移动步长。符号1相对于符号0移动t
step
，符号2相对于符号1移动t
step
，符号3相对于符号2移动t
step
，其他进制以此类推。
[0051]
本实施例中，设定系统的基本参数：带宽bw＝125khz，扩频因子sf＝7，符号周期tc＝1.004ms，频率步进f
step
＝bw/2^sf＝976.56hz；当sym＝28时，该线性扩频信号为将分为三段，其中第一段的起始频率f
01
＝sym*fstep＝54.684khz，第二段的起始频率f
02
＝0hz，第三段的起始频率为bw，即sym28的线性扩频信号为：
[0052][0053]
下面详细说明采样偏差估计模块的工作步骤，其实现流程见图5。
[0054]
本实施例中，设置系统配置：射频rf＝470mhz，采样倍率sampleratio＝4，符号采样点samplenum＝2^sf*sampleratio*1.5＝768，采样点的偏移门限deltasstart＝64(放大32倍)，deltasarea＝32，本实施例以添加了50ppm的jitter为例，假定本实施例中因此进入的ptdmax＝
[0055]
{
‑
4,3,10,13,19,17,22,
‑
23,21,24,11,20,20,11,43,23,9,43,30,26,41,35,19,52,21,37,41,36,44,54}
[0056]
s1：初始化采样偏差估计的符号计数cntnum＝0，判定当前峰值ptdmax是否满足峰值位置的合理性。利用以下条件：
[0057]
ptdmax∈[
‑
deltasstart,deltasstart]
[0058]
如果满足条件则进入步骤s2，同时复位寄存器sum＝0、sum_index＝0、calc＝0、index＝0、real_calc＝0、real_start＝0、buf数组初始化；否则，继续筛选，如果cntnum计数为cntnum＝＝10？仍未找到符合条件的buf[cntnum]，则认为估计失败。
[0059]
这里假定依次进入ptdmax＝{40,45,
‑
4,3,10,13,19,17}，则ptdmax＝{45,46}时不满足条件，则舍弃，当ptdmax＝
‑
4时，满足条件，此时cntnum＝3；并初始化所有的参数，并进入步骤s2。
[0060]
s2：确定峰值位置的有效性，从步骤s1中确定合理范围的cntnum作为起始值，依次将峰值位置ptdmax更新buf数组数据，按如下方式更新：
[0061]
buf[k]＝buf[k
‑
1]；k＝1,2,3,4；其中，buf[0]＝ptdmax。
[0062]
计算delatasstate9buf数组中各个数据之间的差值，这里筛选出首个有效值，筛选条件如下：
[0063]
tmp[k]＝abs(buf[k]
‑
buf[k
‑
1])，且满足tmp[k]∈(0，deltasarea)
[0064]
当首次满足上述条件时，计算buf[0]～buf[4]的均值记为初始有效值real_calc，首次成功数值下标real_start＝cntnum
‑
2；
[0065]
记录当前的有效数值：calc＝ptdmax；下标计数index＝cntnum。
[0066]
计算sum数值：sum+＝calc
‑
real_calc；
[0067]
计算下标sum_index数值：sum_index+＝index
‑
real_start；
[0068]
然后进入步骤s3，否则继续筛选；如果一直不满足此条件，当计数器cntnum＝＝32时，认为估计失败，复位计数器cntnum；
[0069]
本实施例中，buff首次存满后的数据为buf[5]＝{19,13,10,3,
‑
4}，则tmp[4]＝{6,3,7,7}，都满足系统要求；因此real_calc＝8，real_start＝3，并记录当前的有效数值calc＝19，下标计数index＝5；sum＝11，sum_index＝2。
[0070]
s3：采用迭代处理方式，确定之后的峰值有效性，并进行偏差累计。其有效性的判决机制如下：
[0071]
计算当前p路峰值位置ptdmax和上一个记录的有效峰值位置calc的差值：
[0072]
tmp1＝abs(ptdmax
‑
calc)
[0073]
计算当前计数器cntnum与上一个记录的有效峰值位置下标计数器index的差值：
[0074]
tmp2＝cntnum
‑
index
[0075]
判断是否满足以下条件：
[0076]
tmp1≤tmp2
×
deltasarea
[0077]
如果满足条件，更新最新峰值有效值calc和峰值下标：
[0078]
calc＝ptdmax
[0079]
index＝cntnum
[0080]
并同时更新sum和sum_index：
[0081]
sum+＝ptdmax
‑
real_calc
[0082]
sum_index+＝cntnum
‑
real_start
[0083]
重复以上步骤直至cntnum＝＝cntnumthr，这里设定cntnumthr＝30；
[0084]
这里设定本次迭代完成后sum＝440；sum_index＝330。
[0085]
s4：计算采样偏差估计值：
[0086]
ppm_est_calc＝sum/sum_index/samplenum/放大倍率
[0087]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
＝440/336/768/32*1e6
[0088]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
＝53.28ppm
[0089]
对采样偏差的补偿处理，可以采用插值滤波器实现，或调整buffer的取值索引，进行修正，假定本方案采用buffer的调整方案，则每个符号的调整采样点个数为：
[0090]
ppm_comp＝sum/sum_index放大倍率＝440/336/32＝0.0409
[0091]
设定symidx＝0，symidx++；remain_ppm＝0；
[0092]
当floor(symidx*ppm_comp+remain_ppm)＝＝1时，对缓存的接收数据，向前取一个采样点，进行采样偏差的补偿，并更新残留偏差：
[0093]
remain_ppm＝symidx*ppm_comp+remain_ppm
‑
1，symidx＝0。
[0094]
从本实施例中可以看出，本方案能有效地估计采样偏差，估计误差为3.28ppm，且可通过调整cntnumthr以满足不同系统的估计误差要求。
[0095]
本发明具体实施例中充分利用了chirp信号的解扩峰值位置信息，使系统将解扩和ppm的估计有效结合，从而减少了计算资源。
[0096]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明的保护范围。