[发明专利]基于DFT的GNSS导航电文数据解调方法、装置及GNSS接收终端

权利要求书

1.一种基于DFT的GNSS导航电文数据解调方法，其特征在于，包括：

在进行导航电文数据解调之前，对预检测积分结果进行各搜索频点的DFT运算，得到各搜索频点上的DFT运算结果；

判断是否为第一次进行导航电文数据解调；

若是第一次进行导航电文数据解调，则对本次DFT运算结果进行筛选、旋转相位和存储；

若不是第一次进行导航电文数据解调，则读取上一次进行导航电文数据解调后存储的DFT运算结果，并选择上次和本次的频率搜索范围的组合方式，采用猜比特方式计算各组合的能量，选择能量最大的组合对应的电文结果作为本次导航电文数据解调结果；

所述对预检测积分结果进行各搜索频点的DFT运算，得到各搜索频点上的DFT运算结果包括：按照如下公式得到各搜索频点上的DFT运算结果：

DFT_out,i＝Σ_{n＝0,1,…,N-1}(I_n+k+j*Q_n+k)*exp(-j*2*π*Δf_i*n*T₀),i＝1,2,…,M

其中，预检测积分结果为：PDI_out＝DFT_in＝I_k+j*Q_k；I_k和Q_k为正交两路的预检测积分结果；I_k为第k个积分结果，单个积分时间为T₀秒；j满足j*j＝-1；Δf_i为DFT搜索的频率，i＝1,2,…,M，M为总共需要搜索的频点数；其中，用N个预检测积分结果计算一组DFT，总共的预检测积分时间T＝N*T₀小于等于单比特电文数据的宽度，且积分过程未跨越两比特电文数据的边界；

所述对本次DFT运算结果进行筛选、旋转相位和存储包括：遍历b，选择使得(I_b²+Q_b²)最大的几个频点的DFT运算结果，b＝1,2,…,M；对本次DFT运算结果按照如下公式进行旋转相位：I_b⁺+j*Q_b⁺＝(I_b+j*Q_b)*exp(j*2*π*Δf_i*T)；将旋转后的本次DFT运算结果进行存储；

所述选择上次和本次的频率搜索范围的组合方式包括：在上次的DFT运算结果中选择最值得信任的M^-个频点的结果，记为I_a^-+j*Q_a^-，其中，a取1,2,…,M中的M^-个结果，M^-小于等于M；所有组合的数量为M^-*M，其中，本次有M个频点的DFT运算结果，记为I_b+j*Q_b，b＝1,2,…,M；

所述采用猜比特方式计算各组合的能量，选择能量最大的组合对应的电文结果作为本次导航电文数据解调结果包括：遍历a、b和D的所有值，使得[(I_a^-+D*I_b)²+(Q_a^-+D*Q_b)²]有最大值，其中，D为本次解调的电文比特和上次解调的电文比特相比是否翻转，D取1或-1，D为1的时候代表没有翻转，则本次解调的电文结果沿用上次的结果；D为-1的时候代表有翻转，则本次解调的电文结果为上次解调电文结果取反的值；将本次解调的电文结果作为本次导航电文数据解调结果。

2.根据权利要求1所述的方法，其特征在于，Δf_i的取值范围为：[-1/T₀/2Hz,1/T₀/2Hz]。