KR20110048933A - Method and apparatus for acquiring weak gps l1 c/a code modulated signal - Google Patents
Method and apparatus for acquiring weak gps l1 c/a code modulated signal Download PDFInfo
- Publication number
- KR20110048933A KR20110048933A KR1020090105696A KR20090105696A KR20110048933A KR 20110048933 A KR20110048933 A KR 20110048933A KR 1020090105696 A KR1020090105696 A KR 1020090105696A KR 20090105696 A KR20090105696 A KR 20090105696A KR 20110048933 A KR20110048933 A KR 20110048933A
- Authority
- KR
- South Korea
- Prior art keywords
- signal
- digital signal
- code
- integration period
- gps
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
Abstract
In the GPS signal acquisition method, a first digital signal is generated by sampling an input signal corresponding to a GPS signal during a first integration period corresponding to a multiple of a repetition period of a C / A code. The second digital signal is generated by compensating the Doppler effect on the first digital signal during the first integration period. A duplicate carrier corresponding to the second digital signal is generated, and a first complex signal is generated by mixing the second digital signal and the duplicate carrier. The second complex signal is generated by overlapping the first complex signal in units of repetition periods during the first integration period, and a C / A code is obtained through the second complex signal.
Description
The present invention relates to a GPS signal acquisition method and apparatus, and more particularly, to a weak GPS C / A code modulated signal acquisition method and apparatus.
The present invention is derived from a study conducted as part of the IT source technology development project of the Ministry of Knowledge Economy [Task Management Number: 2007-F-040-03, Task name: Indoor and outdoor continuous positioning technology development].
In general, global positioning system (GPS) satellites use two carriers and three binary messages. Two carriers are the L1 carrier and the L2 carrier, and three binary messages are the navigation data message at 50 Hz, the private coarse / acquisition (C / A) code at 1.023 MHz, and the P code at 10.23 MHz. Here the navigation data message contains the system parameters needed to calculate the location. The L1 carrier carries a navigation message modulated with a C / A code and a P code, and the L2 carrier carries a navigation message modulated with a P code. The C / A code and the P code are both pseudo-random codes, but have a repeated property after a certain period of time. The P code repeats after 7 days, but the C / A code repeats after 1 [ms].
Each satellite is assigned a unique C / A code, and the C / A code assigned to each satellite is known to all receivers. The GPS receiver correlates with the C / A code that received its local C / A code to identify satellites in the visible range and to measure the propagation delay time. Since the C / A code is repeated every 1 [ms], the receiver needs an input signal of adequate power for at least 1 [ms] for successful acquisition. This period is called the coherent integration period.
On the other hand, when the GPS signal is weak (for example, a signal having a power of 33 dB-Hz or less) such as a room, a GPS signal, that is, a C / A code cannot be obtained only by the coherent integration period. Acquire C / A codles during the integral period that is a multiple of the runt integral period. However, the magnitude of the Doppler effect can be modeled as a function of the length of the integration period and the Doppler value. The Doppler effect can be ignored for relatively short integration periods such as 1-10 [ms], but in the case of long integration periods, the Doppler effect causes variations in the code phase, making it difficult to obtain C / A codes.
In addition, when the sampling rate used in the analog / digital converter of the GPS receiver is different from the actual sampling rate, a problem of sample fraction occurs. These sample fragments can also be ignored for short integration periods, but for longer integration periods, it is difficult to obtain C / A codes by increasing the deviation of the code phase.
An object of the present invention is to provide a GPS signal acquisition method and apparatus capable of acquiring a weak GPS signal in an indoor environment.
According to an embodiment of the present invention, a method for obtaining a GPS signal of a GPS receiver is provided. The method includes sampling an input signal corresponding to a GPS signal during a first integration period corresponding to a multiple of a repetition period of a C / A code to generate a first digital signal, wherein the first digital signal during the first integration period. Compensating a Doppler effect on the signal to generate a second digital signal, generating a duplicate carrier corresponding to the second digital signal, and mixing the second digital signal and the duplicate carrier to generate a first complex signal. And generating a second complex signal by overlapping the first complex signal in units of the repetition period during the first integration period, and obtaining a C / A code through the second complex signal.
According to another embodiment of the present invention, there is provided an apparatus for obtaining GPS signals including an analog / digital converter, an oscillator, a mixer, a compensator, and a determiner. The analog-to-digital converter samples the input signal corresponding to the GPS signal, and the oscillator generates a duplicate carrier. The mixer generates a complex signal by mixing the digital signal and the replica carrier. The compensator generates the digital signal by compensating the Doppler effect for the output of the analog-to-digital converter during the first integration period corresponding to a multiple of the repetition period of the C / A code, and generates the complex signal during the first integration period. Overlapping in the repeating cycle unit. The determination unit obtains the C / A code based on the output of the compensation unit.
According to another aspect of the present invention, there is provided a method for acquiring a GPS signal of a GPS receiver by sampling an input signal corresponding to a GPS signal during a first integration period corresponding to a multiple of a repetition period of a C / A code. Generating a second digital signal in consideration of a position compensation value according to a sample fragment, an actual number of samples during the first integration period, and a position according to a Doppler effect, from the first digital signal during the first integration period Generating a duplicate carrier corresponding to the second digital signal, generating a first complex signal by mixing the second digital signal and the duplicate carrier, for the first complex signal during the first integration period Generating a second complex signal by overlapping the first complex signal in units of the repetition period while compensating for the difference in points caused by the sample fragments; Multiplying a second complex signal with a local C / A code, accumulating at least a portion of a multiplication result for each first integration period during a second integration period corresponding to a multiple of the first integration period, and peaking from the accumulation result Detecting a point to obtain a C / A code.
DETAILED DESCRIPTION Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those skilled in the art may easily implement the present invention. The present invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. In the drawings, parts irrelevant to the description are omitted in order to clearly describe the present invention, and like reference numerals designate like parts throughout the specification.
Throughout the specification, when a part is said to "include" a certain component, it means that it can further include other components, without excluding other components unless specifically stated otherwise.
1 is a schematic block diagram of a GPS signal acquisition apparatus according to an embodiment of the present invention.
Referring to FIG. 1, the GPS signal acquisition apparatus includes an analog / digital converter (hereinafter referred to as “A / D”)
The A /
The
The
The C /
The
The IFFT
Here, R odd (i) is the multiplication result of the odd-numbered [(2i-1) th] segment among all segments, R even (i) is the multiplication result of the even- numbered [(2i) th] segment, and N is all The number of segments. For example, if the total integration period is 1 second and one segment consists of an integration period of 10 [ms], N is 100, and it is determined that 50 odd-numbered segments and even-numbered segments are accumulated by 50, respectively.
On the other hand, the
As described above, according to an embodiment of the present invention, since a baseband signal of a plurality of coherent integration periods is superposed in one segment, and a result value for the plurality of segments is accumulated, a C / A code is obtained. Weak GPS signals in the environment and the like can be easily obtained. In addition, since the Doppler effect and the sample fragments caused by the long integration period can be compensated, the code phase deviation due to the Doppler effect and the sample fragments can be reduced.
Next, a GPS signal acquisition method according to an embodiment of the present invention will be described in detail. In the following, the coherent integration period corresponding to the repetition period of the C / A code is 1 [ms], and the input signal is processed for an integration period T i of 1 second (= 1000 [ms]). It is assumed that the period of 1 second is divided into 100 segments, and each segment has an integration period T c of 10 [ms].
First, the sample fragment problem and its solution will be described in detail with reference to FIGS. 2 to 6B.
2 is a diagram illustrating a sample fragment problem of the present invention, FIG. 3 is a diagram illustrating a code phase deviation caused by sample fragments for 1 second, and FIG. 4 is a sample fragment compensation according to an embodiment of the present invention. It is a figure which shows. 5A and 5B are diagrams showing the cumulative result of the even and odd segments before compensating the sample fragments, respectively. FIGS. 6A and 6B are the even numbers after compensating the sample fragments according to one embodiment of the present invention. It is a figure which shows the accumulation result of a 1st segment and an odd numbered segment.
Assuming a sampling rate (F s ) of the A /
As shown in FIG. 2, the actual number of
In order to solve such a sample fragment problem, first, the
Here, s has a value from 1 to 100 as an index of a segment, and c 1 is a sample fragment compensation value as given by
In this way, the total number of points read at the current position from the digital signal sampled and output by the A /
Next, the compensator 140 rounds or adds T c × n points to form an n r length vector including a baseband value of 1 [ms]. To this end, the
Here, c 2 is defined as in
Here, ceil () is a function that returns an integer greater than or equal to the argument of ceil ().
Referring to FIG. 4,
Therefore, the
5A, 5B, 6A, and 6B, according to an embodiment of the present invention, when the
Next, the Doppler effect and its solution will be described in detail with reference to FIGS. 7 to 9B.
FIG. 7 is a diagram for explaining the Doppler effect, and FIGS. 8A and 8B are diagrams showing code phases in a case where a positive Doppler effect and a negative Doppler effect occur, and FIGS. 9A and 9B are each an embodiment of the present invention. According to an example, a diagram showing a phase when a positive Doppler effect and a negative Doppler effect are compensated for.
The Doppler offset according to the Doppler effect is modeled as a function of the length of the integration period and the Doppler value. The Doppler effect can be neglected for relatively short integration periods such as 1 to 10 ms. However, for weak signals, such as less than 33 dB-Hz, a long integration period of 1000 times is required depending on the power. This integration, however, requires as many Doppler shifted replica carriers as possible.
According to one embodiment of the present invention, in order to use a non-remarkable Doppler shifted duplicated carrier, it is possible to coherently integrate each segment after dividing the signal into 100 segments as previously assumed. Then all the even-numbered segments may be incoherent integrated and all the odd-numbered segments may be incoherent integrated. In this way, the frequency search spatial resolution can be maintained at 100 Hz (fd) for each segment that is easy to calculate, and sign switching of navigation bits can occur in odd or even segments. That is, odd and even segments may not be involved in sign conversion.
The signal can be extended or reduced by the positive / negative Doppler effect. Therefore, a 1 [ms] signal contains more or fewer samples as Δ samples than n samples. If this Doppler effect is not properly compensated for, an offset of Δ size occurs every millisecond, as shown in FIG. Accumulating such offsets over long integration periods is a large error, which can cause deviations in the code phase estimate in the sign direction of the Doppler value. That is, as shown in FIG. 8A, the positive Doppler effect causes the code phase estimate to decrease with a negative slope as the integral period elapses, and as shown in FIG. As a result, the positive slope increases. This code phase deviation makes peak detection difficult. The following describes in detail how to compensate for this Doppler effect.
First, the
Here, F r is a Doppler shifted copy carrier frequency and is given by Equation (8).
In this case, IF is an intermediate frequency, and f di is a Doppler shift.
Here, f d is 1 / T c = 100, and i is [-d r / f d ... -5 -4 -3 -2 -1 +1 +2 +3 +4 +5 ... + d r / f d ] Either, d r is the Doppler range and has a value of, for example, ± 5000.
First, the
Here, s is an index of the segment and has a value from 1 to 100.
As shown in
Meanwhile, the Doppler compensation value d c is determined by the Doppler compensation point dc p and the index j of the current millisecond, as shown in Equation (11).
Here, the current millisecond index j and the Doppler compensation point dc p are given by Equations 12 and 13, respectively.
Where s is the index of the segment and fo j is the index of the fold in each segment and has a value from 1 to 10.
As shown in Equation 9, the
9A and 9B, it can be seen that the code phase can be kept constant when compensating the Doppler effect as in an embodiment of the present invention.
Next, a GPS signal acquisition method in a GPS signal acquisition apparatus according to an embodiment of the present invention will be described with reference to FIG. 10.
10 is a schematic flowchart of a GPS signal acquisition method according to an embodiment of the present invention.
Referring to FIG. 10, first, the GPS signal obtaining apparatus checks whether an unidentified Doppler shift remains among a plurality of Doppler shifts f di , and if one remains, selects one of the unidentified Doppler shifts and selects a variable. Initialize (S914). That is, the GPS signal obtaining apparatus calculates the Doppler compensation value d c as 0, the pointer p as 1, the index s of the segment as 1, and the Doppler shifted duplicate carrier frequency and the Doppler compensation point, respectively. Initialize with Next, GPS signal acquisition device is to advance the entire integration period (T i) by increasing the index (s) of the segment by 1 (S984) a peak detection process for.
Specifically, the GPS signal obtaining apparatus calculates the sample fragment and the Doppler effect for the corresponding segment as shown in Equation 10 (S920), and considers the sample fragment and the Doppler effect as shown in
Next, the GPS signal acquisition apparatus overlaps the baseband signal in units of a length corresponding to the repetition period of the C / A code (1 [ms] in the previous example) during the T c period for the corresponding segment, and the baseband of 1 [ms]. Generate a signal (S960). Specifically, while increasing the index (fo j ) of the fold in the segment (S962), and compensates for the sample fragment by the sample fragment compensation value (c 2 ) as shown in equation (5) (S964), the Doppler compensation value (d) as shown in equation (11) After calculating c ) (S966), the fold is overlapped with the previous folds to generate a baseband signal of 1 [ms] corresponding to n r points (S968).
The GPS signal obtaining apparatus performs n r point FFT operations on the overlapped baseband signals of 1 [ms] to convert them to the frequency domain (S972), and multiplies the signals converted to the frequency domain by the C / A code (S974).
Subsequently, the GPS signal obtaining apparatus accumulates the multiplication result for the current segment and the multiplication result for the previous segment. In this case, the odd-numbered segment and the even-numbered segment are separated and accumulated (S982). The GPS signal obtaining apparatus determines whether the index s of the current segment is the last index T i / T c (S984), and if not, increases the index s by 1 (S986) and repeats from step S920. If the index s of the current segment is the last index T i / T c , the GPS signal obtaining apparatus outputs a cumulative value of the last odd-numbered segment and an even-numbered segment, and repeats from step S912.
If all the Doppler shifts are confirmed by repeating the above process (S912), the GPS signal acquisition apparatus obtains a C / A code by detecting peak points from the accumulated values (S990).
11A and 11B are diagrams showing cumulative results of even-numbered and odd-numbered segments after compensating sample fragments and Doppler effects, respectively, according to an embodiment of the present invention, and measuring the weak GPS signal of 17 dB-Hz. One result.
As shown in Figs. 11A and 11B, according to one embodiment of the present invention, even a weak GPS signal can easily detect peak points exceeding a threshold.
The embodiments of the present invention described above are not only implemented through the apparatus and the method, but may be implemented through a program for realizing a function corresponding to the configuration of the embodiments of the present invention or a recording medium on which the program is recorded.
Although the embodiments of the present invention have been described in detail above, the scope of the present invention is not limited thereto, and various modifications and improvements of those skilled in the art using the basic concepts of the present invention defined in the following claims are also provided. It belongs to the scope of rights.
1 is a schematic block diagram of a GPS signal acquisition apparatus according to an embodiment of the present invention.
2 is a view for explaining the problem of sample fragments of the present invention.
3 is a diagram showing code phase deviation caused by sample fragments for 1 second.
4 illustrates sample fragment compensation according to an embodiment of the present invention.
5A and 5B are graphs showing the cumulative result of the even and odd segments before compensating the sample fragments, respectively.
6A and 6B are diagrams illustrating cumulative results of even-numbered and odd-numbered segments after compensating sample fragments, respectively, according to one embodiment of the present invention.
7 is a diagram for explaining the Doppler effect.
8A and 8B are diagrams showing code phases when a positive Doppler effect and a negative Doppler effect occur.
9A and 9B are diagrams illustrating phases when the positive Doppler effect and the negative Doppler effect are compensated according to an embodiment of the present invention, respectively.
10 is a schematic flowchart of a GPS signal acquisition method according to an embodiment of the present invention.
11A and 11B are diagrams illustrating cumulative results of even-numbered and odd-numbered segments after compensating sample fragments and Doppler effects, respectively, according to one embodiment of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090105696A KR20110048933A (en) | 2009-11-03 | 2009-11-03 | Method and apparatus for acquiring weak gps l1 c/a code modulated signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090105696A KR20110048933A (en) | 2009-11-03 | 2009-11-03 | Method and apparatus for acquiring weak gps l1 c/a code modulated signal |
Publications (1)
Publication Number | Publication Date |
---|---|
KR20110048933A true KR20110048933A (en) | 2011-05-12 |
Family
ID=44360150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020090105696A KR20110048933A (en) | 2009-11-03 | 2009-11-03 | Method and apparatus for acquiring weak gps l1 c/a code modulated signal |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR20110048933A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160046815A (en) * | 2013-08-27 | 2016-04-29 | 마이크로소프트 테크놀로지 라이센싱, 엘엘씨 | Cloud-offloaded global satellite positioning |
-
2009
- 2009-11-03 KR KR1020090105696A patent/KR20110048933A/en not_active Application Discontinuation
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160046815A (en) * | 2013-08-27 | 2016-04-29 | 마이크로소프트 테크놀로지 라이센싱, 엘엘씨 | Cloud-offloaded global satellite positioning |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2280491B1 (en) | Apparatus and method for estimating a frequency shift and a time shift | |
US8457177B2 (en) | System and method for fast code phase and carrier frequency acquisition in GPS receiver | |
CN109791208B (en) | Method and apparatus for reception of DS/FH signals | |
EP3264622B1 (en) | Low complexity, low power and long range radio receiver | |
US20020025011A1 (en) | System and method for fast code phase and carrier frequency acquisition in GPS receiver | |
JP2016524710A (en) | Receiver and method for direct sequence spread spectrum signals | |
JP2008521012A (en) | Method and apparatus for increasing coherent integration length while receiving positioning signal | |
US7852264B2 (en) | Systems and methods for fast GNSS signals acquisition | |
JP2007322233A (en) | Device for reproducing phase modulation series | |
JP5202540B2 (en) | Method and apparatus for receiving a BOC modulated radio navigation signal | |
JP2006261985A (en) | Receiver for spread spectrum communications | |
CN108562918B (en) | BOC (n, n) ambiguity-free capturing method and device based on correlation shift | |
KR20110048933A (en) | Method and apparatus for acquiring weak gps l1 c/a code modulated signal | |
JP5607606B2 (en) | Multipath detection apparatus and GNSS reception apparatus | |
CN108226968B (en) | Navigation signal rapid capturing method | |
US20060262835A1 (en) | Method and device for estimating a frequency in the acquisition of a spread-spectrum code signal | |
CN107132552B (en) | Parallel code phase searching device and method for realizing parallel code phase searching | |
CN117590446B (en) | Broadband composite navigation signal tracking method and device | |
KR101043689B1 (en) | Method and apparatus for acquiring initial synchronous signal for gps | |
JP4842172B2 (en) | Spread spectrum receiver | |
Tome | Post-Processed Acquisition & Tracking of GPS C/A L1 Signals | |
KR20200017804A (en) | A method for Generating an Unambiguous Correlation Function | |
Ju et al. | An improved CL code acquisition method of L2C signal | |
JP2000324188A (en) | Synchronous generation system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WITN | Withdrawal due to no request for examination |