阅读说明：本技术 一种高灵敏卫星导航信号捕获方法 (High-sensitivity satellite navigation signal capturing method ) 是由 汪平平 党群 郭华昌 田雪涛 于 2019-11-07 设计创作，主要内容包括：本发明涉及一种高灵敏卫星导航信号捕获方法,对接收信号进行滤波、下变频、数字采样后,再与中频载波混频,然后抽取滤波处理、混频、分组后构成二叉树,对二叉树进行寻优,对最优累加结果进行FFT后和各卫星PRN码的FFT值的共轭值进行相乘再进行IFFT处理,选出最大值即为多普勒、码相位。本发明通过对二叉树进行寻优,利用比特信息,可以实现不同场景下微弱卫星信号的快速捕获。(The invention relates to a high-sensitivity satellite navigation signal capturing method, which comprises the steps of filtering, down-converting and digitally sampling a received signal, mixing the signal with an intermediate frequency carrier, extracting, filtering, mixing and grouping to form a binary tree, optimizing the binary tree, multiplying an optimal accumulation result by a conjugate value of an FFT value of each satellite PRN code after FFT, and then performing IFFT processing to select a maximum value, namely Doppler and a code phase. The invention can realize the rapid capture of weak satellite signals in different scenes by optimizing the binary tree and utilizing bit information.)

1. A high-sensitivity satellite navigation signal acquisition method is characterized by comprising the following steps:

step 1: filtering, down-converting and digitally sampling the received radio frequency signal to obtain a baseband signal _ rec, and mixing the signal _ rec and a local intermediate frequency carrier local _ if to obtain a complex signal _ if;

step 2: performing Q-time CIC decimation filtering processing on the signal _ if, and enabling the sampling rate of the signal to be equal to the intermediate frequency sampling rate f_sIs reduced to f_s', get signal _ dif; q times is more than 2 times of bandwidth;

and step 3: mixing signal _ dif and local _ freq _ bin to obtain data, wherein the local _ freq _ bin is Doppler compensation;

and 4, step 4: grouping data, every f_s'/1000 points are a group, N groups are formed, a binary tree with the depth of floor (N/20) +1 is formed, the degree of each internal node is 2, two different state paths are represented, and the accumulated metric value is acc (k);

and 5: optimizing the binary tree, initializing first, the priority queue containing the start node n₀Initializing the cumulative metric value acc (n)₀)＝0；

Step 6: continuous optimizerThe node n with the maximum accumulated metric value in the priority queue_kOutputting, wherein the node has the best path in the traditional Viterbi gridding graph;

and 7: output node n of the optimizer_kInserting a priority queue, the contents of the insertion including n_kThe time, state and forward pointer of the device; if the node is not a terminal node, n is also added_kIs connected to the successor node n_k+1Inserting into a priority queue structure, inserting node n_k+1The cumulative metric of (c) is acc (n)_k+1)＝acc(n_k)+d(n_k,k+1)a_kAfter inserting the node into the priority queue, returning to the step 6; stopping the optimizing if the trellis has been expanded to the terminal node;

and 8: accumulating the data by using the priority queue to obtain data b_1,k,b_2,k,…b_n,k(ii) a The accumulation mode is as follows:

set up b_1,0,b_2,0,…b_n,kThe initial value is 0, and the initial value is 0,

m, d, wherein k is 1,2_k∈{-1,1}；

And step 9: performing FFT on the optimal accumulation result to obtain rec _ FFT;

step 10: generating FFT values of PRN codes of each satellite in advance by using a local code table, and performing conjugation processing on the FFT values to obtain a signal local _ FFT;

step 11: multiplying the values obtained in the steps 9 and 10 to obtain rst _ fft;

step 12: performing IFFT processing on the rst _ fft to obtain a result, storing the result into the acqu _ mat, and returning to the step 3 until all local _ freq _ bins are traversed to obtain the acqu _ mat;

step 13: and comparing acqu _ mat to obtain the maximum value, wherein the corresponding value is Doppler and code phase.

Technical Field

The invention relates to the technical field of satellite navigation, in particular to a weak signal processing method based on a receiver.

Background

In complex environments such as urban canyons, forests, indoor environments and the like, GNSS signals are shielded, serious attenuation can be generated, multipath and interference exist, the signal power is lower than the working range of a normal GPS receiver, and the usability and the positioning accuracy of the GNSS receiver are greatly reduced. High-sensitivity GNSS reception technology has become a hotspot in navigation technology research.

The high-sensitivity technology mainly improves the processing gain of signals through coherent integration and non-coherent integration, and extracts GNSS signals from background noise. The selection of the coherent integration time is limited by data bit inversion and Doppler frequency difference, if the data bit value, the edge position and the real-time Doppler frequency difference can be obtained for code Doppler compensation, long-time coherent integration can be carried out, and the signal gain is improved.

At present, weak signal acquisition mainly comprises coherent detection, incoherent detection and differential detection. The coherent detection capture algorithm can provide the maximum gain under the condition of white gaussian noise. However, the maximum coherent integration time is affected by crystal oscillator stability, bit flipping, doppler, and other factors; the non-coherent accumulation is insensitive to phase changes caused by navigation data and doppler frequency, but there is also a "squaring loss", which is more severe as C/N0 decreases.

Psiaki proposes a half-bit alternating capture technique. Dividing the received 20ms data block into two sections, ensuring that the data of one section is completely in the same bit, respectively performing coherent integration on the two sections of data, and then further performing non-coherent accumulation on the result of the coherent integration. Ziedan presents a weak signal acquisition algorithm based on navigation data prediction, but the accuracy of data prediction and the complexity of the algorithm are still not effectively solved.

In addition to navigation data flipping, long coherent integration time leads to increased frequency bin search, resulting in significantly increased frequency domain search times, limiting coherent integration time.

Disclosure of Invention

Technical problem to be solved

In order to overcome the defects and shortcomings of the existing satellite weak signal capturing technology, the invention provides a high-sensitivity satellite navigation signal capturing method which is suitable for the American GPS and the 'Beidou' system in China and can realize the rapid capturing of weak satellite signals in different scenes.

Technical scheme

A high-sensitivity satellite navigation signal acquisition method is characterized by comprising the following steps:

step 1: filtering, down-converting and digitally sampling the received radio frequency signal to obtain a baseband signal _ rec, and mixing the signal _ rec and a local intermediate frequency carrier local _ if to obtain a complex signal _ if;

step 2: performing Q-time CIC decimation filtering processing on the signal _ if, and enabling the sampling rate of the signal to be equal to the intermediate frequency sampling rate f_sIs reduced to f_s', get signal _ dif; q times is more than 2 times of bandwidth;

and step 3: mixing signal _ dif and local _ freq _ bin to obtain data, wherein the local _ freq _ bin is Doppler compensation;

and 4, step 4: grouping data, every f_s'/1000 points are a group, N groups are formed, a binary tree with the depth of floor (N/20) +1 is formed, the degree of each internal node is 2, two different state paths are represented, and the accumulated metric value is acc (k);

and 5: optimizing the binary tree, initializing first, the priority queue containing the start node n₀Initializing the cumulative metric value acc (n)₀)＝0；

Step 6: the optimizer continuously converts the node n with the maximum accumulated metric value in the priority queue_kOutputting, wherein the node has the best path in the traditional Viterbi gridding graph;

and 7: output node n of the optimizer_kInserting a priority queue, the contents of the insertion including n_kThe time, state and forward pointer of the device; if the node is not a terminal node, n is also added_kIs connected to the successor node n_k+1Inserting into a priority queue structure, inserting node n_k+1The cumulative metric of (c) is acc (n)_k+1)＝acc(n_k)+d(n_k,k+1)a_kAfter inserting the node into the priority queue, returning to the step 6; stopping the optimizing if the trellis has been expanded to the terminal node;

and 8: accumulating the data by using the priority queue to obtain data b_1,k,b_2,k,…b_n,k(ii) a The accumulation mode is as follows:

set up b_1,0,b_2,0,…b_n,kThe initial value is 0, and the initial value is 0,

m, d, wherein k is 1,2_k∈{-1,1}；

And step 9: performing FFT on the optimal accumulation result to obtain rec _ FFT;

step 10: generating FFT values of PRN codes of each satellite in advance by using a local code table, and performing conjugation processing on the FFT values to obtain a signal local _ FFT;

step 11: multiplying the values obtained in the steps 9 and 10 to obtain rst _ fft;

step 12: performing IFFT processing on the rst _ fft to obtain a result, storing the result into the acqu _ mat, and returning to the step 3 until all local _ freq _ bins are traversed to obtain the acqu _ mat;

step 13: and comparing acqu _ mat to obtain the maximum value, wherein the corresponding value is Doppler and code phase.

Advantageous effects

According to the high-sensitivity satellite navigation signal capturing method, the binary tree is optimized, and bit information is utilized, so that weak satellite signals in different scenes can be rapidly captured.

Drawings

FIG. 1 is a flow chart of signal acquisition according to the present invention;

FIG. 2 is a data packet diagram of the present invention;

FIG. 3 is a diagram of a binary tree structure according to the present invention.

Detailed Description

The invention will now be further described with reference to the following examples and drawings:

the technical solution of the present invention is further explained by taking GPS signals as an example through a specific implementation manner with reference to the accompanying drawings.

(1) Filtering, down-converting and digitally sampling a received radio frequency signal of 99ms to obtain a baseband signal _ rec, and mixing the signal _ rec and a local intermediate frequency carrier local _ if to obtain a complex signal _ if;

(2) CIC (common information center) extraction filtering processing is carried out on the signal _ if, and the sampling rate of the signal is reduced to 2.046M from the intermediate frequency sampling rate of 16.384M to obtain signal _ dif;

(3) mixing signal _ dif with local _ freq _ bin, wherein the local _ freq _ bin is used for fine Doppler compensation in a 1KHz range to obtain data;

(4) grouping data, wherein each 2046 points form one group, the total number of the groups is 99, the starting point is k, k is 1,2.. 19, each 20 groups form one block, a binary tree with the depth of 5 is formed, the degree of each node in the binary tree is 2, two different state paths [1, -1] are represented, and the path accumulation metric value is acc (k);

(5) optimizing the binary tree, initializing first, the priority queue containing the start node n₀Initializing the cumulative metric value acc (n)₀)＝0；

(6) The optimizer continuously converts the node n with the maximum accumulated metric value in the priority queue_kAt the output, this node has the best path in the conventional Viterbi trellis diagram.

(7) Output node n of the optimizer_kInserting a priority queue, the contents of the insertion including n_kTime of day, state, and forward pointer. If the node is not a terminal node, n is also added_kIs connected to the successor node n_k+1Inserting into a priority queue structure, inserting node n_k+1The cumulative metric of (c) is acc (n)_k+1)＝acc(n_k)+d(n_k,k+1)a_kAfter inserting the node into the priority queue, returning to the step (6); if the trellis has expanded to the last underlying node, then the optimization is stopped.

(8) Accumulating the data by using the priority queue to obtain data b_1,k,b_2,k,…b_n,k. The accumulation mode is as follows:

set up b_1,0,b_2,0,…b_n,kThe initial value is 0, and the initial value is 0,

m, d, wherein k is 1,2_k∈{-1,1}；

Obtaining an accumulated value rec _ data;

(9) and performing FFT on the optimal accumulation result rec _ data to obtain rec _ FFT.

(10) FFT values of PRN codes of each satellite are generated in advance by using a local code table, and conjugation processing is carried out on the FFT values to obtain a signal local _ FFT.

(11) Multiplying the values obtained in the steps (9) and (10) to obtain the rst _ fft.

(12) And (3) performing IFFT processing on the rst _ fft to obtain a result, storing the result into acqu _ mat (: k), returning to the step (3) until all local _ freq _ bin is traversed to obtain acqu _ mat, returning to the step (4), and returning to the step (13) when k is 19.

(13) And comparing acqu _ mat to obtain the maximum value, wherein the corresponding value is the capture Doppler and the code phase.