一种基于小波域滤波码捕获方法
【技术领域】
[0001] 本发明属于导航信号检测技术领域，涉及一种基于小波域滤波码捕获方法。
【背景技术】
[0002] 全球导航卫星系统（GlobalNavigationSatelliteSystem，简称GNSS)是全球性 的位置与时间测定系统，它包括卫星星座、地面监护系统及用户终端设备，可为地球表面、 近地面和地球外空任意地点的用户提供全天候、实时、高精度的三维位置、速度以及精密的 时间信息。GPS系统、Galilea系统和GLONASS系统都属于全球导航卫星系统的范畴。
[0003]GNSS信号是经过直接序列扩频调制的扩频信号，对经过直接扩频调制的卫星导航 信号的捕获是GNSS系统需要解决的首要问题。串行捕获方法是一种信号捕获技术，但此 方法采用时域计算接受码和本地码相关，一次计算只能得到一个码相位的相关值，且计算 量大，捕获耗时较长。为弥补串行捕获方法的自身缺陷，出现了并行捕获方法。并行捕获 方法是将时域计算接收码和本地码相关转化为频域方法计算，且频域有相应的快速傅里叶 变换（FastFourierTransformation，FFT)。并行捕获方法虽节省了大量的捕获时间，但 由于需要捕获接收GNSS弱信号，需要大数据量的相关积累。为了更加快速、准确的计算接 收码和本地码相关，部分匹配滤波（partialmatchedfilter，PMF)方法应运而生，其中具 有代表性的方法是基于非相干积累的部分匹配滤波（partialmatchedfilterbasedon non-coherentaccumulation，PMF_NCA)方法。
[0004] 基于非相干积累的部分匹配滤波方法的基本过程可以描述为：将进行下变频后的 基带信号进行部分匹配滤波，然后进行非相干积累，最后将积累的判决量进行阈值判定。采 用非相干积累有两点好处：一是提高了信噪比进而提高了检测概率；二是克服了导航数据 调制所带来的比特翻转问题。但是非相干积累采用将信号平方进行积累的方式，这种是积 累方式也将噪声幅度也给平方了，所以非相干积累提高检测概率的能力是有限的。
[0005] 由于非相干积累增大了噪声幅度，限制了提高检测概率的能力，因此，需要探寻一 种克服现有基于非相干积累的部分匹配滤波方法的途径，从而提高信号的信噪比，将会更 加有利于对GNSS弱信号的捕获。
【发明内容】
[0006] 本发明所要解决的技术问题是针对上述现有技术提供一种基于小波域滤波码捕 获方法，该方法能有效地滤除噪声，能更大程度地提高信号的信噪比，进而提高检测概率的 能力。
[0007] 本发明上述技术问题所采用的技术方案为：一种基于小波域滤波码捕获方法，其 特征在于包括如下步骤：
[0008] 步骤1、对接收的中频信号进行下变频处理，得到基带信号s>);
[0009] 步骤2、将基带信号s>)分为Q个数据段，各数据段表示为(η)，…，& (η)，…， (η)，分别将各个数据段进行快速傅里叶变换，然后再与本地码c(n+m)的快速傅里叶变 换的复共辄进行相乘，得到的乘积再分别取快速傅里叶逆变换，获得与各个数据段相对应 的相关函数，用札(η)、!?>)、& (η)分别表示基带信号s>)每一数据段对应的相关函数：
[0011] 其中，FFT0表示取信号快速傅里叶变换；IFFT0表示求信号的快速傅里叶逆变 换；conj()表示求信号的复共辄；
[0012] 步骤3、对步骤2获得的相关函数进行分段处理获取同步相关值向量；
[0013] 步骤4、对步骤3获得的同步相关值向量进行小波变换，获得小波系数；
[0014] 步骤5、对步骤4获得的小波系数进行低通滤波；
[0015] 步骤6、对低通滤波后的小波系数进行小波反变换得到重构信号；
[0016] 步骤7、对重构信号进行非相干积累，并构造判决统计量；
[0017] 步骤8、门限判决。
[0018] 所述基带信号的获取方法为：
[0019] 对接收的中频信号r(η)进行下变频处理，得到基带信号s(n):
[0021] 其中，sjn)为基带信号；P(n)为中频信号r(n)的噪声，且为复高斯噪声，实部、虚 部均服从均值为〇,方差为σ2的高斯分布；为信号频率，Af为搜索步长，i为搜索次数； r(n)为接收的中频信号；h(l)为信道的衰落因子；c(n+p)为接收的C/A码，C/A码为全球导 航卫星系统发出的一种伪随机码；P为初始相位；η= 0,. . .，N-1 ;N为采样点数；
[0022] 所述步骤3中对步骤2获得的相关函数进行分段处理获取同步相关值向量的方法 为：
[0023] 在Rq(n)中，设0<η<Ρ，其中PXQ=N，N为采样点数，当本地码c(n+m)与接收 的C/A码c(n+p)同步时，即m=p，则相关函数的同步相关值向量s为：
[0024] s=[Ri(p) ...Rq (p) ...Rq (p)]
[0025] 其中，Rq(p)为第qms同步时的相关值，即
[0027] 其中，q= 1，. . .，Q，
：Co为lms的点数；h⑴为信道的衰落因子，f。为基 带信号的频率;W(q)为基带信号噪声且为复高斯噪声，实部、虚部都服从均值为0,方差为Co*σ2的高斯分布。
[0028] 所述步骤4中对步骤3获得的同步相关值向量进行小波变换获得小波系数的方法 为：
[0029] 对同步相关值向量s进行小波变换：
[0031] 其中，A;(q)djq)表示小波基函数；J为小波分解的级数；sr(l)，sr(2)，· · ·，sr(k) ，...，S(J+1)为与小波基函数对应的小波系数。
[0032] 所述步骤6获得重构信号的方法为：
[0033] 设低频段的小波系数\(1)进行低通滤波后的值为\(1)';
[0034] 对低频段的小波系数\(1)进行低通滤波后的值\(1)'进行小波反变换得到的 重构信号为：
[0036] 其中AJ(q)为低频小波基函数；q= 1，. . .，Q。
[0037] 所述步骤7的方法为：
[0038] 对重构信号sf(q)进行非相干积累，并构造判决统计量Z:
[0040] 其中Co为1ms的点数；〇 2为噪声方差；Q为数据段数。
[0041] 所述步骤8中的门限判决的方法为：设γ为判决门限，判决门限γ为预先设定的 经验值，或者通过设先设定的虚警概率求取获得，将判决门限γ与决统计量Ζ进行比较，若 决统计量Ζ大于判决门限γ，表示接收端检测到信号，若决统计量Ζ小于判决门限γ则表 示接收端没有检测到信号。
[0042] 本发明还可以包括如下步骤：步骤9、根据构造的判决统计量Ζ得出总的虚警概率 Pfa 为：
[0044] 其中P表示带搜索C/A码相位数；γ为判决门限，判决门限γ为预先设定的经验 值，或者通过设先设定的虚警概率求取获得，Ρ{Ζ〈γ}表示当Ζ小于判决门限γ的概率。
[0045] 与现有技术相比，本发明的优点在于：
[0046] 1、本发明对获取的信号进行频域相关处理，使得本方法能一次相关计算所有码相 位对应的相关值。
[0047] 2、本方法引入了小波域滤波处理信号的方式，保证了滤波后的有用信号能量没有 减少且使滤波后噪声信号的能量和方差减小，提高了检测概率。
【附图说明】
[0048] 图1为基于小波域滤波码捕获方法的处理流程示意图。
【具体实施方式】
[0049] 以下结合附图实施例对本发明作进一步详细描述。
[0050] 本实施例提供了一种基于小波域滤波码捕获方法，从图1可以看出，接收的中频 信号依次经过下变频单元、频域相关处理、分段处理、小波变换、低通滤波、小波反变换、非 相干积累单元和门限判决等处理步骤，详细为：
[0051] 步骤1、对接收的中频信号进行下变频处理，得到基带信号s>);
[0052] 对接收的中频信号r(η)进行下变频处理，得到基带信号sjn):
[0054] 其中，sjn)为基带信号；P(n)为中频信号r(n)的噪声，且为复高斯噪声，实部、虚 部均服从均值为〇,方差为σ2的高斯分布；fi为信号频率，△f为搜索步长，i为搜索次数； r(n)为接收的中频信号；h(l)为信道的衰落因子；c(n+p)为接收的C/A码，C/A码为全球导 航卫星系统发出的一种伪随机码；P为初始相位；η= 0,. . .，N-1 ;N为采样点数
[0055] 步骤2、频域相关处理：将基带信号\(η)分为Q个数据段，各数据段表示为 srl(n)，…，srq(n)，…，Sr(j(n)，分别将各个数据段进行快速傅里叶变换，然后再与本地码 c(n+m)的快速傅里叶变换的复共辄进行相乘，得到的乘积