阅读说明：本技术 一种高灵敏度导航接收机tow解析方法 (TOW analysis method of high-sensitivity navigation receiver ) 是由 杨克元 张蓬 田嘉 王延光 钱小平 于 2019-09-16 设计创作，主要内容包括：本发明提供了一种高灵敏度导航接收机TOW解析方法,该方法利用导航电文TOW的特点,对多个连续子帧TOW相应位置的比特交错累加,提高了每个比特的信噪比,降低了解调误码率,可应用于高灵敏度导航接收机TOW解调。该方法的应用,解决了低信噪比环境下,导航信号解调误码率过高导致TOW无法正常解调的问题。方法简单易实现,具有良好的嵌入型,只需要在传统的导航接收机帧同步后加上该累加模块即可实现低信噪比下TOW的可靠解调。(The invention provides a TOW analysis method of a high-sensitivity navigation receiver, which utilizes the characteristics of a navigation message TOW to accumulate the bits at corresponding positions of a plurality of continuous sub-frames TOW in a staggered way, improves the signal-to-noise ratio of each bit, reduces the demodulation error rate and can be applied to the TOW demodulation of the high-sensitivity navigation receiver. The method solves the problem that the TOW can not be demodulated normally due to overhigh demodulation error rate of the navigation signals in the environment with low signal-to-noise ratio. The method is simple and easy to implement, has good embedded type, and can realize the reliable demodulation of the TOW under low signal-to-noise ratio only by adding the accumulation module after the frame synchronization of the traditional navigation receiver.)

1. A TOW analysis method of a high-sensitivity navigation receiver is characterized by comprising the following steps: the method comprises the following steps:

1) dividing the TOW into a high bit part and a low bit part, and setting the low bit part of the TOW as b_(N-1)To b₀The high bit part of TOW in the week is b_NTo b₁₆Wherein b is₀Is the lowest order, b_NIs the highest order;

2) store K +1 frame data, denoted F_kK is 0,1, 2.. ·, K; from item F₀When a frame starts, the position of all '0' of the low bit part of the TOW in the week is searched, and the frame data frame number corresponding to the position is marked as L;

3) demodulating the high bit part of the TOW in the week;

4) the TOW2 at the current week is calculated from the demodulation results of the low bit part and the high bit part.

2. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 1, wherein: the specific process of the step 2) is as follows:

21) let the nth bit of the kth frame be denoted b_knSearching the position of the 0 th bit with the 1 st value as 0, and performing superposition operation to obtain the result recorded as sum _ b₀If sum _ b₀Less than 0, k is 0 or b₀The first value is the position of 0, otherwise, k is 1 and b₀The first value is the position of 0, the frame number of the position of bit0 with the first value of 0 is recorded as k₀；

22) Search the 1 st bit 1 st 0 position by k₀As a starting point, calculating to obtain the 1 st '0' position of the 1 st bit, and performing superposition operation to obtain a result which is recorded as sum _ b₁: if sum _ b₁Less than 0, then k₀Is b is₁The first value is the position of 0, otherwise k₀+2 is the position where the first value is 0, b₁The position frame number with the first value 0 is marked as k₁；

23) Search bit 2 for the position of value 0 of 1, in k₁As a starting point, the 2bit is at the 1 st '0' position, and the superposition operation is performed, and the result after the superposition is obtained is recorded as sum _ b₂(ii) a If sum _ b₂Less than 0, then k₁Is b is₂The first value is the position of 0, otherwise k₁+4 is the position where the first value is 0, b₂The position frame number with the first value 0 is marked as k₂；

24) According to the method, the position of the 1 st value of the 3 rd, 4 th, … … th bit, N-1 st bit is '0', the position of the 1 st '0' of the nth bit is searched, the superposition operation is carried out, and the result after the superposition is obtained is recorded as sum _ b_n(ii) a Searching out the low bit part b₀To b_(N-1)The frame data frame number corresponding to the position of all 0 values is denoted as L.

3. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 2, wherein: in the step 21), searchingSearching the position of the 0 th bit 1 with the value of 0, and when the superposition operation is carried out, taking positive when k is even number and taking negative when k is odd number, and recording the result as: sum _ b₀＝b₀₀-b₁₀+b₂₀-b₃₀+b₄₀-b₅₀+......。

4. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 2, wherein: the above-mentioned

sum_b₁＝b_(k0)1+b_(k0+1)1-b_(k0+2)1-b_(k0+3)1+b_(k0+4)1+b_(k0+5)1+......。

5. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 2, wherein: the above-mentioned

sum_b₂＝b_(k1)2+b_(k1+1)2+b_(k1+2)2+b_(k1+3)2-b_(k1+4)2-b_(k1+5)2-b_(k1+6)2-b_(k1+7)2+b_(k1+8)1......。

6. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 2, wherein: the above-mentioned

sum_b_n＝b_(k1)n+b_(k1+1)n+……+b_(k1+2^n-1)n-b_(k1+2^n)n-……，n＝3,4，……，N-1。

7. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 2, wherein: the specific process of the step 3) is as follows:

31) from the lowest bit b of the high bit part_NStarting, accumulation starting from position L, 2^ N<K-L<2^ (N +1), the accumulation length is selected to be 2^ N, the accumulation value is recorded as S0, the data accumulation value between the position 0 and the position L-1 is recorded as S1, the accumulation values of 2^ N +1 to K are recorded as S2, and the sum (N) is S0-S1-S2 if sum (N)>0, then sum (n) 1; otherwise, sum (n) ═ 0;

32) after the Nth bit is demodulated, the demodulation is continued from the (N +1) th bit, and the demodulation b_nWhen a bit in (N +1, N + 2...., 16) is to be 'or' operated on the value that the bit has been previously demodulated, i.e., sum _ or (N) or sum (N +1) or..; if sum _ or is 0, when the current bit is demodulated, the accumulated result between the position 0 and the position L-1 needs to be inverted; otherwise, the accumulated result between the start at position 0 to position L-1 need not be inverted.

8. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 7, wherein: the process of the step 4) is as follows:

41) position week TOW1 values for which the low bit fraction is all zero:

42) the current time of week TOW2 ═ TOW1-L was calculated from TOW 1.

9. The TOW parsing method for the time of week of the high sensitivity navigation receiver according to claim 8, wherein: the TOW1 ═ sum (N) ^ 2^ N + sum (N +1) ^ 2^ (N +1) +. + -. sum (16) × 2^ (16).

10. The time-of-week TOW analysis method for the high-sensitivity navigation receiver according to any one of claims 1-9, wherein: the selection range of N in step 31) is an integer from 1 to 16.

Technical Field

The invention relates to the technical field of navigation receivers, in particular to a time-of-week TOW analysis method for a high-sensitivity navigation receiver.

Background

Compared with a low-orbit satellite, a medium-high orbit satellite mainly uses a side lobe signal of a navigation satellite or a leakage signal from the other side of the earth, so that the received signal power is very weak, and the signal-to-noise ratio received by a satellite-borne navigation receiver is as low as 0dB (the corresponding carrier-to-noise ratio of the navigation signal is 17dBHz) or less. In the signal processing process of the navigation receiver, the time synchronization can be carried out only by correctly demodulating the time of week (TOW), the receiving signal-to-noise ratio of the ground or low-orbit satellite navigation receiver is higher, the time of week TOW demodulation can be directly carried out, but under the low signal-to-noise ratio, the error rate of the time of week TOW demodulation is higher, and the time synchronization is difficult to use. Therefore, it is necessary to reduce the demodulation error rate of TOW in cycles by other means.

The prior art is analyzed, relevant systems such as a universal database, a national defense science and technology information resource service system, a national science and technology book literature center and the like are consulted, in the relevant literature data which can be found at present, the research on the performance of capturing and tracking is more important for a high-sensitivity navigation receiver, most of the research is focused on theoretical simulation, and the demodulation method of the TOW in the week with low signal to noise ratio is not introduced in detail, such as: GNSS for LEO, GEO, HEO and beyond. Feasibility study of GNSS as navigation system to reach the Moon (authors: Vincenzo Capuano et al); standalone GPS L1C/A Receiver for Lunar Missions (author: Vincenzo Capuano et al); design of a top high sensitivity acquisition system for a space GNSS receiver (author: L.Musumici et al); kalman filter-based architecture for robust and high-sensitivity tracking in GNSS receivers (authors: Peral-Rosado et al).

Disclosure of Invention

The technical problem solved by the invention is as follows: the method divides the TOW into a high bit part and a low bit part by utilizing the self characteristics of the TOW in the navigation signal week, finds out the position of all zeros of the low bit part as a reference point for calculating the TOW value in the week, and reduces the complexity of realization; the method of bit interleaving accumulation is adopted, so that the demodulation signal-to-noise ratio of the TOW in the week is improved, and the error probability of the TOW in the week is reduced.

The above object of the present invention is achieved by the following scheme: a TOW analysis method of a high-sensitivity navigation receiver comprises the following steps:

1) dividing the TOW into a high bit part and a low bit part, and setting the low bit part of the TOW as b_(N-1)To b₀The high bit part of TOW in the week is b_NTo b₁₆Wherein b is₀Is the lowest order, b_NIs the highest order;

2) store K +1 frame data, denoted F_kK is 0,1, 2.. ·, K; from item F₀When a frame starts, the position of all '0' of the low bit part of the TOW in the week is searched, and the frame data frame number corresponding to the position is marked as L;

3) demodulating the high bit part of the TOW in the week;

4) the TOW2 at the current week is calculated from the demodulation results of the low bit part and the high bit part.

The specific process of the step 2) is as follows:

21) let the nth bit of the kth frame be denoted b_knSearching the position of the 0 th bit with the 1 st value as 0, and performing superposition operation to obtain the result recorded as sum _ b₀If sum _ b₀Less than 0, k is 0 or b₀The first value is the position of 0, otherwise, k is 1 and b₀The first value is the position of 0, the frame number of the position of bit0 with the first value of 0 is recorded as k₀；

22) Search the 1 st bit 1 st 0 position by k₀As a starting point, calculating to obtain the 1 st '0' position of the 1 st bit, and performing superposition operation to obtain a result which is recorded as sum _ b₁: if sum _ b₁Less than 0, then k₀Is b is₁The first value is the position of 0, otherwise k₀+2 is the position where the first value is 0, b₁The position frame number with the first value 0 is marked as k₁；

23) Search bit 2 for the position of value 0 of 1, in k₁As a starting point, the 2bit is at the 1 st '0' position, and the superposition operation is performed, and the result after the superposition is obtained is recorded as sum _ b₂(ii) a If sum _ b₂Less than 0, then k₁Is b is₂The first value is the position of 0, otherwise k₁+4 is the position where the first value is 0, b₂Position frame sequence with first value 0Number k₂；

24) According to the method, the position of the 1 st value of the 3 rd, 4 th, … … th bit, N-1 st bit is '0', the position of the 1 st '0' of the nth bit is searched, the superposition operation is carried out, and the result after the superposition is obtained is recorded as sum _ b_n(ii) a Searching out the low bit part b₀To b_(N-1)The frame data frame number corresponding to the position of all 0 values is denoted as L.

In step 21), the position of the 0 th bit 1 st value is searched for and is 0, when the superposition operation is performed, the positive value is taken when k is an even number, the negative value is taken when k is an odd number, and the result is recorded as:

sum_b₀＝b₀₀-b₁₀+b₂₀-b₃₀+b₄₀-b₅₀+......。

the above-mentioned

sum_b₁＝b_(k0)1+b_(k0+1)1-b_(k0+2)1-b_(k0+3)1+b_(k0+4)1+b_(k0+5)1+......。

The above-mentioned

sum_b₂＝b_(k1)2+b_(k1+1)2+b_(k1+2)2+b_(k1+3)2-b_(k1+4)2-b_(k1+5)2-b_(k1+6)2-b_(k1+7)2+b_(k1+8)1......。

The above-mentioned

sum_b_n＝b_(k1)n+b_(k1+1)n+……+b_(k1+2^n-1)n-b_(k1+2^n)n-……，n＝3,4，……，N-1。

The specific process of the step 3) is as follows:

31) from the lowest bit b of the high bit part_NStarting, accumulation starting from position L, 2^ N<K-L<2^ (N +1), the accumulation length is selected to be 2^ N, the accumulation value is recorded as S0, the data accumulation value between the position 0 and the position L-1 is recorded as S1, the accumulation values of 2^ N +1 to K are recorded as S2, and the sum (N) is S0-S1-S2 if sum (N)>0, then sum (n) 1; otherwise, sum (n) ═ 0;

32) after the Nth bit is demodulated, the demodulation is continued from the (N +1) th bit, and the demodulation b_nWhen a bit is selected from (N +1, N + 2.. times.16),an or operation is performed on the value that the bit has been demodulated before, i.e., sum _ or _ sum (N) or sum (N +1) or..; if sum _ or is 0, when the current bit is demodulated, the accumulated result between the position 0 and the position L-1 needs to be inverted; otherwise, the accumulated result between the start at position 0 to position L-1 need not be inverted.

The process of the step 4) is as follows:

41) position week TOW1 values for which the low bit fraction is all zero:

42) the current time of week TOW2 ═ TOW1-L was calculated from TOW 1.

The TOW1 ═ sum (N) ^ 2^ N + sum (N +1) ^ 2^ (N +1) +. + -. sum (16) × 2^ (16).

The selection range of N in step 31) is an integer from 1 to 16.

Compared with the prior art, the invention has the following beneficial effects:

(1) the performance is excellent: the method has the advantages that the self characteristics of the TOW of the navigation signal in-week time are utilized, the TOW of the in-week time is divided into a high bit part and a low bit part, and the positions of all zeros of the low bit part are found out to be used as reference points for calculating the TOW value of the in-week time, so that the complexity of implementation is reduced; the demodulation signal-to-noise ratio of the TOW in the week is improved and the error probability of the TOW in the week is reduced by adopting a bit interleaving accumulation method;

(2) the adaptability is good: the method is suitable for the time-of-week TOW demodulation of all currently known navigation constellations and a new navigation constellation which has the same characteristics with the time-of-week TOW of the current navigation signal;

(3) the realization is simple: the method provided by the patent is simple and easy to implement, has a good embedded type, and can realize the reliable demodulation of the TOW in the lower period of the ground signal-to-noise ratio only by adding the accumulation module after the frame synchronization of the traditional navigation receiver.

Drawings

FIG. 1a is a schematic diagram of the division of TOW high ratio bits in a week;

FIG. 1b is a schematic diagram of a bit interleaved accumulation process;

fig. 2 is a processing flow chart of the time-of-week TOW analysis method of the high-sensitivity navigation receiver according to the present invention.

Detailed Description

The invention is described in further detail below with reference to the following figures and specific examples:

as shown in the method flowchart of fig. 2, the multiple access interference suppression method for a spread spectrum receiver based on spectrum sensing of the present invention includes the following steps:

1) dividing the time of week TOW into a high bit portion and a low bit portion, the length of each portion being related to the total accumulated length, assuming that the low bit portion is b_(N-1)To b₀The high bit part is b_NTo b₁₆Wherein b is₀Is the lowest order, b_NIs the highest order, as shown in FIG. 1 a;

2) store K +1 frame data, denoted F_kK is 0,1,2, from F₀At the beginning of the frame, the position of all '0' of the low bit part is found, and for the convenience of description, the nth bit of the kth frame can be represented as b_kn：

(21) Searching the position of the 0 th bit with the 1 st value as 0, taking positive when k is even number and taking negative when k is odd number, performing superposition operation, and recording the result as: sum _ b₀＝b₀₀-b₁₀+b₂₀-b₃₀+b₄₀-b₅₀+.., if sum _ b₀Less than 0, k is 0 or b₀The first value is the position of 0, otherwise, k is 1 and b₀The first value is the position of 0, the frame number of the position of bit0 with the first value of 0 is recorded as k₀，k₀Is 0 or 1;

(22) searching the 1 st bit 1 st 0 position by k₀As a starting point, unlike the 0 th bit, the 1 st bit changes its symbol every two consecutive frames, so that during accumulation, the same symbols are directly accumulated, and different symbols are accumulated after inversion, and the result is recorded as:

sum_b₁＝b_(k0)1+b_(k0+1)1-b_(k0+2)1-b_(k0+3)1+b_(k0+4)1+b_(k0+5)1+.., if sum _ b₁Less than 0, then k₀Is b is₁The first value is the position of 0, otherwise k₀+2 is the position where the first value is 0, b₁The position frame number with the first value 0 is marked as k₁，k₁Is k₀Or k₀+2；

(23) Search bit 2 for the position of value 0 of 1, in k₁As a starting point, the 2 nd bit changes sign every 4 consecutive frames, and the accumulation mode is b₁Similarly, the same symbols are accumulated directly, and different symbols are accumulated after being inverted, as shown in fig. 1b, and the result is recorded as:

sum_b₂＝b_(k1)2+b_(k1+1)2+b_(k1+2)2+b_(k1+3)2-b_(k1+4)2-b_(k1+5)2-b_(k1+6)2-b_(k1+7)2+b_(k1+8)1.., if sum _ b₂Less than 0, then k₁Is b is₂The first value is the position of 0, otherwise k₁+4 is the position where the first value is 0, b₂The position frame number with the first value 0 is marked as k₂，k₂Is k₁Or k₁+4；

(24) By analogy, the low bit part b can be searched₀To b_(N-1)Positions with all values of 0 are marked as L;

(25) in the above process, from b₁At the beginning, the accumulated starting position may not be the stored 0 th frame, so that the data from the 0 th frame to the starting position is not utilized, and in order to improve the utilization efficiency of the stored data and reduce the loss of the signal to noise ratio, the partial data can be accumulated, and the partial data is accumulated to the total accumulated sum after being inverted.

3) Demodulation of TOW high bit part in week:

(31) from the lowest bit b of the high bit part_NStarting, accumulation starting from position L, 2^ N<K-L<2^ (N +1), the accumulation length is 2^ N, the accumulation value is S0, the data accumulation value from the position 0 to the position L-1 is S1, the accumulation values from 2^ N +1 to K are S2, then sum (N) is S0-S1-S2 if sum (N)>0, then sum (n) 1; otherwise, sum (n) is 0.

(32) Demodulation b_nWhen a bit in (N +1, N + 2...., 16), the value that the bit has been demodulated before is subjected to an or operation: note sum _ or ═ sum (N) or sum (N +1) or.. If sum _ or is 0, when the current bit is demodulated, position 0 is onThe accumulated result between the start position L-1 needs to be inverted; otherwise, the accumulated result between the start at position 0 to position L-1 need not be inverted.

(33) The demodulation result according to the high bit part of the TOW in the week is H.

4) According to the demodulation results of the low bit part and the high bit part, the TOW2 in the current week is calculated:

the TOW1 values at the location cycle with low bit fractions all zero are:

TOW1＝sum(N)*2^N+sum(N+1)*2^(N+1)+......+sum(16)*2^(16)

the value of TOW2 at the current week is:

TOW2＝TOW1-L。

the above description is only one embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.

Those skilled in the art will appreciate that the invention may be practiced without these specific details.