CN103499825B - Method for rapidly capturing BDS weak signals of high-orbit spacecraft - Google Patents

Method for rapidly capturing BDS weak signals of high-orbit spacecraft Download PDF

Info

Publication number
CN103499825B
CN103499825B CN201310438975.8A CN201310438975A CN103499825B CN 103499825 B CN103499825 B CN 103499825B CN 201310438975 A CN201310438975 A CN 201310438975A CN 103499825 B CN103499825 B CN 103499825B
Authority
CN
China
Prior art keywords
data
bds
code
coherent accumulation
carried out
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310438975.8A
Other languages
Chinese (zh)
Other versions
CN103499825A (en
Inventor
闻长远
柯颋
王猛
李兴国
高阳
王烁
陆华
岳富占
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Space Star Technology Co Ltd
Original Assignee
Space Star Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Space Star Technology Co Ltd filed Critical Space Star Technology Co Ltd
Priority to CN201310438975.8A priority Critical patent/CN103499825B/en
Publication of CN103499825A publication Critical patent/CN103499825A/en
Application granted granted Critical
Publication of CN103499825B publication Critical patent/CN103499825B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a method for rapidly capturing BDS weak signals of a high-orbit spacecraft. The method includes the steps that NH code phases are removed through 2ms data overlapping and accumulation and autocorrelation of NH codes, weak signal capturing is finished by the cooperation of 20ms coherent accumulation with multiple times of incoherent accumulation, a parallel code phase search strategy based on FFT/IFFT is adopted, and therefore capturing speed is improved. By the adoption of the method, the problem that coherent accumulation time can not be prolonged because the NH code phases of the BDS signals are unknown can be solved, capturing sensitivity of the BDS signals is high, the capturing speed of the BDS signals is high, and the method can be widely applied to a BDS navigation receiver which is used for executing a high-orbit task, a moon exploration return task and other spaceflight tasks, and has wide popularization and application prospect.

Description

A kind of high rail spacecraft BDS weak signal quick capturing method
Technical field
The present invention relates to a kind of high rail spacecraft BDS(triones navigation system) weak signal fast Capture Technique, can be used for exchanging the satellite navigation signals being shaped with the secondary codings such as NH code (Neumann-Hoffman code) and carry out weak signal fast Acquisition, belong to satellite application field.
Background technology
High rail spacecraft all plays an important role in recent years in land and overocean communications, aerological sounding, educational applications, live telecast, disaster early warning etc., particularly geostationary orbit spacecraft, with the stationary characteristic over the ground of its uniqueness, become communication, meteorology, scouting, time service and follow the tracks of and the important track resources of the earth satellite such as data relay and scientific research.The quantity of China to high rail spacecraft demand increases year by year, and the independent navigation of high rail spacecraft has become the space flight new technology of China in the urgent need to development.
GPS (Global Position System) relies on its global, round-the-clock, continuous and high-precision feature, plays an important role in national defence, national security, economic security and social life.China just has the BDS (triones navigation system) of Global coverage ability in early implementation construction, expect the year two thousand twenty and possess basic service ability.High rail spacecraft use navigation neceiver to become the focus of all kinds of engineer applied.
Because high rail spacecraft orbit height is higher than navigation constellation, receiver needs to receive the navigation satellite signal from earth opposite.The earth block and the increase of signal propagation distance causes the factors such as received power reduction, will cause using conventional receiver technology to there is the shortcomings such as visible star number reduces, the navigation signal duration is short, received signal power is faint.
Solving and receiving the faint basic ideas of pilot signal power is adopt multicycle signal energy accumulation, by the Cumulative sum (coherent accumulation, non-coherent accumulation) to multiple 1ms correlation, brings signal to noise ratio (S/N ratio) to promote.But there is NH code in BDS, will in acquisition procedure, the problem causing coherent accumulation length limited because NH code phase is unknown; If only adopt non-coherent accumulation process, now need data to be processed by huger than data during coherent accumulation, reduce acquisition speed simultaneously, thus cause high rail to receive the difficulty of BDS signal.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, a kind of high rail spacecraft BDS weak signal quick capturing method is provided, the present invention can overcome because NH code phase is unknown and limit the problem that coherent accumulation length causes efficiently promoting accumulated process gain, can be used for designing the BDS navigation neceiver with higher acquisition sensitivity level, practical.
Technical solution of the present invention is:
A kind of high rail spacecraft BDS weak signal quick capturing method, comprises step as follows:
(1) obtain the 1ms data under local values of Doppler frequency shift, by its zero padding to 2ms, and FFT conversion is carried out to these data, after getting conjugation, obtain data 1;
(2) carry out continuous print 2ms data buffer storage respectively after the intermediate-freuqncy signal of reception being carried out digital sample, namely 1ms and 2ms carries out buffer memory, and 2ms and 3ms carries out buffer memory, and buffer memory 20 groups, obtains data 2 successively;
(3) the 2ms data of often organizing in the data 2 obtained in step (2) be multiplied with a chip of 20 bit NH codes successively respectively, described NH code is 0,0,0,0,0,1,0,0,1,1,0,1,0,1,0,0,1,1,1,0, obtain 20 groups to be multiplied with NH code after data 3;
(4) coherent accumulation is carried out to 20 groups of 2ms data that step (3) obtains in data 3, obtains the coherent accumulation value that length is 2ms, and carry out FFT operation obtain data 4;
(5) respectively data 4 are multiplied with data 1 successively contraposition, obtain data 5, IFFT conversion is carried out to data 5 again and obtains data 6;
(6) repeatedly non-coherent accumulation is carried out to data 6, obtain 1 group of data 7;
(7) Nai Man-Pierre's Si criterion is utilized to carry out detection threshold judgement to data 7, judge whether successfully lock-on signal, if successful lock-on signal enters step (9), otherwise judge whether 20 NH code phases have been searched for, if do not have search to complete, enter step (8), if searched for the local values of Doppler frequency shift of change to enter step (1);
(8) to NH code ring shift chip of 20 bits, enter step (2), repeat step (2) ~ (7);
(9) end is caught.
The present invention's beneficial effect is compared with prior art:
(1) when the present invention is directed to coherent accumulation, the unknown of NH code phase causes carrying out long-time accumulation problem, adopts and accumulates the overlap of 2ms data, eliminate this impact, extends the coherent accumulation time to 20ms, thus improves the acquisition sensitivity of BDS signal.
(2) present invention employs the mode of overlapping accumulated value being carried out to non-coherent accumulation, improve processing gain further.
(3) the present invention adopts and to walk abreast code phase fast search based on FFT/IFFT, greatly in engineering must improve acquisition speed, reduces the resource loss brought 20msNH code auto-correlation computation.
Accompanying drawing explanation
Fig. 1 is the inventive method process flow diagram;
Fig. 2 is the NH code searching principle schematic of 2ms data overlap of the present invention accumulation;
Fig. 3 the present invention is that input signal is multiplied with NH code schematic diagram;
Fig. 4 the present invention is NH code auto-correlation schematic diagram;
Fig. 5 the present invention is 2ms data overlap accumulated process schematic flow sheet.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is further described in detail.
One provided by the invention high rail BDS weak signal quick capturing method, can be used for design and there is weak signal processing power BDS navigation neceiver, the present invention removes NH code phase by the accumulation of 2ms data overlap, complete the weak signal acquisition procedure that 20ms coherent accumulation coordinates several times non-coherent accumulation, and adopt the parallel code phase search strategy based on FFT/IFFT, improve acquisition speed.
As shown in Figure 1, concrete steps of the present invention are as follows:
(1) obtain the 1ms data under local values of Doppler frequency shift, by its zero padding to 2ms, and FFT conversion is carried out to these data, after getting conjugation, obtain data 1;
(2) carry out continuous print 2ms data buffer storage respectively after the intermediate-freuqncy signal of reception being carried out digital sample, namely 1ms and 2ms carries out buffer memory, and 2ms and 3ms carries out buffer memory, and buffer memory 20 groups, obtains data 2 successively;
(3) the 2ms data of often organizing in the data 2 obtained in step (2) be multiplied with a chip of 20 bit NH codes successively respectively, described NH code is 0,0,0,0,0,1,0,0,1,1,0,1,0,1,0,0,1,1,1,0, obtain 20 groups to be multiplied with NH code after data 3;
(4) the 20 groups of 2ms data obtained in data 3 step (3) carry out coherent accumulation (such as: the 1st group of 2ms data and the 2nd group of 2ms accumulate, result after accumulation is accumulated with the data of the 3rd group of 2ms again, until 20 groups of 2ms data coherency accumulations terminate), obtain the coherent accumulation value that length is 2ms, and carry out FFT operation obtain data 4;
(5) respectively data 4 are multiplied with data 1 successively contraposition, obtain data 5, IFFT conversion is carried out to data 5 again and obtains data 6;
(6) repeatedly non-coherent accumulation is carried out to data 6, obtain 1 group of data 7; Number of times data 6 being carried out to non-coherent accumulation is relevant with acquisition sensitivity, and the acquisition sensitivity of requirement is higher, and it is more to carry out incoherent accumulative number of times;
(7) Nai Man-Pierre's Si criterion is utilized to carry out detection threshold judgement to data 7, judge whether successfully lock-on signal, if successful lock-on signal enters step (9), otherwise judge whether 20 NH code phases have been searched for, if do not have search to complete, enter step (8), if searched for the local values of Doppler frequency shift of change to enter step (1);
(8) to NH code ring shift chip of 20 bits, enter step (2), repeat step (2) ~ (7);
(9) end is caught.
In conjunction with specific engineering construction case, provide technology embodiment and simulation results:
(1) 2ms data overlap accumulative removal NH code phase design
As shown in Figure 2,3, the NH code searching principle of 2ms data overlap accumulation and flow process, be shown in Fig. 2 to the treatment scheme of 2ms data.
Suppose that the phase differential of input signal and NH code sequence is ξ, the result that now continuous print 2ms data are multiplied with NH code is successively: D k=[(1-ξ) N k, N k+1, ξ N k+2] × N k+1
By 20 groups of D kdata are added successively:
Σ D k = ( 1 - ξ ) × Σ k = 1 20 N k × N k + 1 + Σ k = 1 20 N k + 1 × N k + 1 + ξ × Σ k = 1 20 N k + 2 × N k + 1 = ( 1 - ξ ) × R N ( 1 ) + 20 × R N ( 0 ) + ξ × R N ( - 1 )
R in formula n() is autocorrelation value corresponding to the different code phase of NH code.As shown in Figure 4, be the autocorrelation function of NH code, known at code phase 1 and-1 place according to Fig. 4, the autocorrelation value of NH code is 0, and the autocorrelation value at code phase 0 place is maximum.Therefore from above formula computing, by the accumulation to 20 groups of 2ms data, middle 1ms can eliminate the impact of NH code phase completely, and Fig. 5 is shown in by concrete schematic diagram.
Adopt which, the middle 1ms data of 2ms accumulation gained are relevant to local 1ms data, find correlation peak, be equivalent to the gain having carried out 20 coherent accumulations.
The content be not described in detail in instructions of the present invention belongs to the known technology of those skilled in the art.

Claims (1)

1. a high rail spacecraft BDS weak signal quick capturing method, is characterized in that step is as follows:
(1) obtain the 1ms data under local values of Doppler frequency shift, by its zero padding to 2ms, and FFT conversion is carried out to these data, after getting conjugation, obtain data 1;
(2) carry out continuous print 2ms data buffer storage respectively after the intermediate-freuqncy signal of reception being carried out digital sample, namely 1ms and 2ms carries out buffer memory, and 2ms and 3ms carries out buffer memory, and buffer memory 20 groups, obtains data 2 successively;
(3) the 2ms data of often organizing in the data 2 obtained in step (2) be multiplied with a chip of 20 bit NH codes successively respectively, described NH code is 0,0,0,0,0,1,0,0,1,1,0,1,0,1,0,0,1,1,1,0, obtain 20 groups to be multiplied with NH code after data 3;
(4) coherent accumulation is carried out to 20 groups of 2ms data that step (3) obtains in data 3, obtains the coherent accumulation value that length is 2ms, and carry out FFT operation obtain data 4;
(5) respectively data 4 are multiplied with data 1 successively contraposition, obtain data 5, IFFT conversion is carried out to data 5 again and obtains data 6;
(6) repeatedly non-coherent accumulation is carried out to data 6, obtain 1 group of data 7;
(7) Nai Man-Pierre's Si criterion is utilized to carry out detection threshold judgement to data 7, judge whether successfully lock-on signal, if successful lock-on signal enters step (9), otherwise judge whether 20 NH code phases have been searched for, if do not have search to complete, enter step (8), if searched for the local values of Doppler frequency shift of change to enter step (1);
(8) to NH code ring shift chip of 20 bits, enter step (2), repeat step (2) ~ (7);
(9) end is caught;
In step (6), the number of times of non-coherent accumulation carries out to data 6 relevant with acquisition sensitivity.
CN201310438975.8A 2013-09-24 2013-09-24 Method for rapidly capturing BDS weak signals of high-orbit spacecraft Active CN103499825B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310438975.8A CN103499825B (en) 2013-09-24 2013-09-24 Method for rapidly capturing BDS weak signals of high-orbit spacecraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310438975.8A CN103499825B (en) 2013-09-24 2013-09-24 Method for rapidly capturing BDS weak signals of high-orbit spacecraft

Publications (2)

Publication Number Publication Date
CN103499825A CN103499825A (en) 2014-01-08
CN103499825B true CN103499825B (en) 2015-05-27

Family

ID=49865048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310438975.8A Active CN103499825B (en) 2013-09-24 2013-09-24 Method for rapidly capturing BDS weak signals of high-orbit spacecraft

Country Status (1)

Country Link
CN (1) CN103499825B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133225B (en) * 2014-07-29 2016-12-07 上海北伽导航科技有限公司 Big Dipper signal processing method and device
CN104360357A (en) * 2014-11-24 2015-02-18 四川九洲电器集团有限责任公司 Quick Beidou satellite signal capturing method and system based on circulation mode
CN104459731B (en) * 2014-11-27 2017-01-04 上海交通大学 A kind of high rail Satellite Orbit Determination method of the weak GNSS signal of receiving pole
CN104483684B (en) * 2015-01-05 2017-06-23 中国科学院重庆绿色智能技术研究院 A kind of method of fast Acquisition Big Dipper D1 satellite navigation system weak signals
CN104765048B (en) * 2015-04-02 2017-04-12 西安电子科技大学 High-sensitivity Beidou satellite B1I signal capturing method
CN107132554B (en) * 2016-02-26 2019-10-25 大唐半导体设计有限公司 A kind of parallel code phase search device and the method for realizing parallel code phase search
CN105911569B (en) * 2016-04-06 2019-06-14 东南大学 The weak D1 signal acquisition methods of Beidou satellite navigation system B1 frequency band based on sequence of symhols estimation
CN107290761B (en) * 2016-04-13 2020-07-28 大唐半导体设计有限公司 Parallel code phase searching device and method for realizing parallel code phase searching
WO2018107441A1 (en) * 2016-12-15 2018-06-21 深圳开阳电子股份有限公司 Signal capturing method and receiver for satellite navigation system
CN109782311B (en) * 2019-02-15 2020-08-21 航天恒星科技有限公司 Method and apparatus for capturing BDS signals

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928853A (en) * 2012-10-25 2013-02-13 西安华迅微电子有限公司 Method for capturing weak signals of Big Dipper D1 satellite navigation system
CN103091686A (en) * 2011-10-31 2013-05-08 北京大学 Bit synchronization method and device of compatible Beidou satellite navigation system and global position system (GPS)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9000980B2 (en) * 2011-07-13 2015-04-07 Texas Instruments Incorporated GNSS receiver correlating doppler derotation samples with code phases

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103091686A (en) * 2011-10-31 2013-05-08 北京大学 Bit synchronization method and device of compatible Beidou satellite navigation system and global position system (GPS)
CN102928853A (en) * 2012-10-25 2013-02-13 西安华迅微电子有限公司 Method for capturing weak signals of Big Dipper D1 satellite navigation system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高轨环境下BDS弱信号跟踪技术研究;闻长远等;《飞行器测控学报》;20130831;第32卷(第4期);第363-370页 *

Also Published As

Publication number Publication date
CN103499825A (en) 2014-01-08

Similar Documents

Publication Publication Date Title
CN103499825B (en) Method for rapidly capturing BDS weak signals of high-orbit spacecraft
CN104253774B (en) Doppler shift estimating system and method under a kind of high dynamic environment
CN101029923B (en) Software receiver for indoor position and navigation technology of global positioning system
CN104280750B (en) Long code capture systems based on partial matched filter fft algorithm and method
CN103033828B (en) High-sensitivity compass-assisted time servicing device, time service receiver and time service method
CN106291619B (en) A kind of wireless long-period spread spectrum code signal high-performance catching method of aeronautical satellite inter-satellite link
CN101246210B (en) GPS locating method
CN103954977B (en) A kind of GNSS cheating interference cognitive method and system
CN103091687B (en) Big dipper receiver multi-satellite signal joint capturing method
CN102608626A (en) High-sensitivity satellite navigation signal capturing method and system
CN204101732U (en) A kind of GNSS New System signal capture device
CN106646541A (en) Beidou weak signal capture method based on difference correlation integral
CN104765050A (en) Novel Beidou signal secondary acquisition algorithm
CN104459731A (en) Quite-weak GNSS signal receiving high-orbit-satellite orbit positioning method
CN104898136A (en) Capturing method of Beidou second-generation B1 signals and system thereof
CN103543456A (en) Large frequency offset GNSS signal capture method based on segmentation relative combination FFT operation
CN103885071A (en) Satellite signal capturing method and equipment
CN104765052B (en) GEO navigation satellite high-sensitivity carrier tracking method
CN105759288A (en) Feature sequence-based Beidou B1I weak signal acquisition method
CN104880698B (en) Based on the space maneuver object detection method converted apart from frequency domain polynomial-phase
CN103809192A (en) Dynamic correction algorithm of GNSS receiver
CN106291610A (en) A kind of parallel correlation module of compression for GNSS signal compression capture processing means and its implementation
CN102243309A (en) Method and apparatus for restraining cross-correlation interference in GNSS
CN104714241A (en) Rapid GPS bit synchronization method
CN2854622Y (en) System for GPS position at small signal environment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant