CN106291612B - A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method - Google Patents
A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method Download PDFInfo
- Publication number
- CN106291612B CN106291612B CN201610618004.5A CN201610618004A CN106291612B CN 106291612 B CN106291612 B CN 106291612B CN 201610618004 A CN201610618004 A CN 201610618004A CN 106291612 B CN106291612 B CN 106291612B
- Authority
- CN
- China
- Prior art keywords
- satellite
- wireless signal
- maximum
- inter
- judgment method
- 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
Links
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
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)
Abstract
The invention discloses a kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, its step is：S1：Digital quadrature demodulation is carried out to reception signal, the complex base band signal of I Q two-way composition is obtained after LPF；S2：Resampling is carried out to complex base band signal, obtains detection statistic by coherent integration, detection and non-coherent integration afterwards, the maximum searched in whole detection statistics, and record the lattice site of maximum；After the completion of first run detection, hunting zone is adjusted, performs more wheel detections, and the uniformity of maximum lattice site is tested；S3：Whether succeeded according to assay judgement capture.The present invention have the advantages that adaptation star between time-varying characteristics, improve acquisition sensitivity and realize fast Acquisition.
Description
Technical field
Present invention relates generally to radio communication and field of measuring technique, refers in particular to one kind and is applied to aeronautical satellite inter-satellite link
The high-performance prize judgment method of wireless signal.
Background technology
GPS (Global Navigation Satellite System, GNSS) can be the earth and
The anywhere of terrestrial space provides round-the-clock precision positions and temporal information.In GPS, maintain compared with
High satellite orbit determines that precision and clock correction determine that precision is to ensure that navigation receiving terminal obtains the positioning of big system requirements or awarded
Shi Jingdu key.
Aeronautical satellite precise orbit determination and time synchronized have expedited the emergence of out chain between navigation constellation star to the demand of H_2O maser function
The concept on road.Once establishing contact between aeronautical satellite by inter-satellite link, the space segment of navigation system will no longer be isolated defends
The combination of star, turn into the entirety mutually cooperateed with.By the cooperation of satellite-ground link, the control section of whole navigation system and space
Duan Zhenzheng forms round-the-clock, round-the-clock a Seamless Network, and this provides huge for the service operation management of navigation system
Performance space.Inter-satellite link can realize satellite navigation system in the case where only configuring a small number of monitoring stations by inter-satellite link
Accurate measurement obtains the metrical information of other segmental arcs on track, so as to reach the ability for obtaining Precise Orbit parameter.
Navigation constellation inter-satellite link network is a kind of complicated satellite network, has the acentric spy of typical flattening
Sign, it is a wireless network with certain amount peer node.Navigation constellation inter-satellite link network needs real present time about
Multi-multipoint measurement under the conditions of beam, to complete accurate measurement function, navigational satellite system is broadcast and connect by inter-satellite link
The precise distance measurement between spread spectrum distance measuring signal progress satellite is received, the capture of spread-spectrum signal is a pseudo-code and carrier wave two-dimensional search
Process, the size of hunting zone directly determine the speed degree of signal capture and realize difficulty.When inter-satellite link uses radio
When distance-finding method completes precise distance measurement with time synchronized, the satellite relative motion in different orbital planes is larger, gives measurement signal band
Carry out larger Doppler frequency shift, so as to increase the capture range of signal, certain difficulty is brought to capture, especially for
The limited situation of computing resource on star.When between any two satellites in constellation establish measurement communication link when, interstellar distance and
Doppler's excursion is larger, in addition, wireless signal is often using length between measurement performance and the angle of security, star
Cycle spreading code, if not utilizing prior information, the two-dimensional search lattice point enormous amount of code phase and Doppler are real to capture
Now bring great difficulty.
For navigation constellation, time synchronized and precise orbit determination are the bases of navigation system operation, the satellite in constellation
It is in a high-precision space-time datum.Inter-satellite link system is believed using the ephemeris and clock correction information prediction of link setup satellite
Number arrival time and Doppler, can be and more by delay time search scope control corresponding to code phase within ten musec orders
The hunting zone of general Le is less than 100Hz.The reduction of capture range reduce to a certain degree capture realize difficulty, but still face with
Lower several respects are difficult：
(1) from system application angle, generally require to be switched fast link setup object, cause signal duration between star
It is short, therefore, it is necessary to receiver is rapidly completed capture.
(2) interstellar distance is remote, it is desirable to which acquisition sensitivity is high.
(3) satellite borne equipment process resource anxiety is, it is necessary to fully improve resource utilization.
Therefore, the capture of inter-satellite link wireless signal, substantially in certain hunting zone, spaceborne set using limited
Standby resource realizes weak signal fast Acquisition.
It is straight that classical signal acquisition procedure includes detection statistic generation, peakvalue's checking and threshold judgement, the selection of thresholding
Connect and have influence on acquisition performance, thresholding is too low to increase false-alarm probability, and thresholding selection is too high, can reduce acquisition probability；Therefore,
The determination of thresholding is most important.Under the conditions of CFAR criterion, decision threshold is distributed with noise probability and its statistical value is relevant.Star
Between link for being switched fast the demand of link setup object, cause antenna to point to saltus step, cause the time-varying characteristics of ambient noise, make
Obtain each link setup and be required for statistical noise distribution situation again, however, under the conditions of fast capturing signal, less number can only be relied on
Noise statisticses are carried out according to sample, easily causes the deviation of statistical result, causes thresholding misalignment, so as to have influence on acquisition performance.Separately
Outside, the setting of thresholding improves the energy scale of detection statistic, can reduce acquisition sensitivity, is unfavorable for weak signal capture.Most
The realization of whole catching method also needs to the resource problem in view of satellite borne equipment, optimizes and improves the utilization of resources of catching method
Rate.
The content of the invention
The technical problem to be solved in the present invention is that：For technical problem existing for prior art, the present invention provides one
Kind can improve acquisition sensitivity, the aeronautical satellite inter-satellite link wireless communication for being adapted to weak signal capture, can realizing fast Acquisition
Number high-performance prize judgment method.
In order to solve the above technical problems, the present invention uses following technical scheme：
A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, its step are：
S1：Digital quadrature demodulation is carried out to reception signal, the complex baseband of I Q two-way composition is obtained after LPF
Signal；
S2：Resampling is carried out to complex base band signal, examined afterwards by coherent integration, detection and non-coherent integration
Statistic, the maximum searched in whole detection statistics are surveyed, and records the lattice site of maximum；After the first run terminates, continue
More wheel detections are performed, and the uniformity of maximum lattice site is tested；
S3：Whether succeeded according to assay judgement capture.
As a further improvement on the present invention：In the step S2, more wheel grid searches and peakvalue's checking are performed, i.e.,
Realize that capture result is adjudicated by the consistency check to lattice site corresponding to peak value；After first run search is completed, subsequently search
Rope scope phase place change situation will be set according to caused by Doppler.
As a further improvement on the present invention：In the step S2, the quadrature demodulation refers to drive according to carrier doppler
Dynamic load ripple NCO generates orthogonal local carrier signal；The signal resampling refers to phase-accumulated according to pseudo-code Doppler driving
Device generates signal resampling enable signal.
As a further improvement on the present invention：The detailed process of the step S2 is：
S201：Phase grid search scope is determined according to prior information and its uncertainty；The prior information includes star
Between wireless signal arrival time TOA and arrival rate FOA；
S202：Calculate the detection statistics result of lattice point in hunting zone；
S203：Searching and detecting result maximum, and lattice site corresponding to record；
S204：Searched for if the first run, follow-up " consistency check " stage Search scope is determined according to lattice site, and return
S202, otherwise into S205；
S205：Consistency check is performed to maximum lattice site；The consistency check is：For every wheel maximum pair
The lattice site situation of change answered is counted, and thinks front and rear " consistent " when change is less than 1 chip；
S206：If not completing one-hundred-percent inspection, S202 is returned, is searched for into next round.
As a further improvement on the present invention：In step s3, it is described to confirm that capture result refers to：If continuous N time is examined
As a result " consistent ", then judge acquisition success, otherwise, it is determined that capture failure.
As a further improvement on the present invention：Detection statistics result described in the step S202, using time-domain calculation side
Method or frequency-domain calculations method, result of calculation are receiving sequence and the coherent integration or non-coherent integration results of local pseudo-code sequence.
As a further improvement on the present invention：" consistency check " stage Search scope is described in the step S204
Refer to：Follow-up " consistency check " stage Search scope is set centered on lattice site corresponding to first run peakvalue's checking result.
As a further improvement on the present invention：If one-hundred-percent inspection described in the step S206 is completed to refer to continuous N time inspection
Result " consistent " is tested, then terminates to examine, otherwise, to the end of N takes turns detection, completes one-hundred-percent inspection.
As a further improvement on the present invention：The detailed process of the step S205 is：
S2051：According to the hunting zone of setting, the detection statistic of whole lattice points is calculated, calculating process, which includes being concerned with, to be accumulated
Point, detection and non-coherent integration；
S2052：Maximum is searched in whole detection statistics, and records lattice point phase
S2053：Compared with last round of maximum lattice point phase, defining conformance criteria isIt is if full
The foot criterion adds up to counter m, otherwise resets；
S2054：If m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see whether reach amount of testing, if
N reaches N, terminates acquisition procedure, and because m is not up to M, announces capture failure；
S2055：If n is not up to N, examined into next round.
As a further improvement on the present invention：Using Doppler's prior information to receiving sequence code phases in acquisition procedure
Carry out first-order dynamic compensation；If the precision deficiency of first-order dynamic compensation, is further compensated, the elder generation for compensation using second order dynamic
Testing information includes arriving signal frequency and frequency change rate.
Compared with prior art, the advantage of the invention is that：
1st, aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the invention, avoids decision threshold meter
Calculating misalignment causes acquisition performance to decline, and is advantageous to improve acquisition sensitivity, is adapted to weak signal capture, by adjusting hunting zone,
Resource consumption is not only reduced, and fast Acquisition can be realized.
2nd, aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the invention, it is general in order to reduce false-alarm
Rate, more wheels search of maximum is performed, and compensated by code phase dynamic, the uniformity of code phase true value is kept, in order to enter one
Step improves acquisition search speed, has readjusted the hunting zone in " consistency check " stage.
3rd, aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the invention, avoids traditional thresholding and sentences
The problems such as statistics misalignment, computing resource consumption and the acquisition sensitivity certainly brought reduces, not only fully meets inter-satellite link body
System causes the fast capturing signal demand under noise characteristic time dependant conditions, can also be widely used in all types of spaceborne and ground and expand
Frequency receiver.
Brief description of the drawings
Fig. 1 is the inter-satellite link wireless signal principle schematic diagram of ranging communicating integral.
Fig. 2 is the present invention principle schematic that inter-satellite link wireless signal captures in specific application.
Fig. 3 is the schematic flow sheet that the present invention scans in specific application.
Fig. 4 is the schematic flow sheet that the present invention carries out consistency check in specific application.
Fig. 5 is the principle schematic that the present invention carries out code phase dynamic compensation function in specific application.
Embodiment
The present invention is described in further details below with reference to Figure of description and specific embodiment.
As shown in figure 1, be the inter-satellite link wireless signal structural representation of ranging communicating integral, including measurement branch road with
Communication leg.Wherein, measurement branch road is free of data message, and only including spreading code, communication leg then contains data message, and passes through
Band spectrum modulation, two branch roads are modulated according to UQPSK modes generates wireless signal between star after carrier wave.According to the signal of the structural generation
Model is：
In formula：
j：Represent satellite number；
Ac：Expression is modulated to the ranging spreading code amplitude of each frequency carrier wave I branch roads；
Ap：Expression is modulated to the communications spread code amplitude of each frequency carrier wave Q branch roads；
C：Represent I branch road ranging spreading codes；
P：Represent Q branch road communications spread codes；
Dp：Represent the numeric data code modulated on Q branch road communications spread codes；
f：Represent inter-satellite link carrier frequency；
Represent inter-satellite link measurement channel carrier wave first phase；
Represent inter-satellite link communication channel carrier wave first phase.
Receiver receives wireless signal between star, by tracking measurement branch road pseudo-code or carrier phase, realize H_2O maser and
Time synchronized function, and aid in completing communication leg data-transformation facility, therefore, receiver is captured only for measurement branch road
Processing.
The present invention central principle be：The detection statistics result of whole lattice points is calculated in the region of search, searches maximum
And record lattice site where it；Then testing result is being calculated repeatedly simultaneously in section near lattice point where original maximum
Maximum is searched for, while records lattice site corresponding to often wheel search, eventually through the uniformity of lattice site where maximum
Examine and whether determine acquisition success.
As shown in Figures 2 and 3, aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the invention, its
Concretely comprise the following steps：
S1：Digital quadrature demodulation is carried out to reception signal, the complex baseband of I Q two-way composition is obtained after LPF
Signal；
S2：Carry out resampling to signal, resampling rate is twice of pseudo- bit rate；Afterwards by coherent integration, detection and
Non-coherent integration obtains detection statistic, the maximum searched in whole detection statistics, and records the lattice site of maximum；
After the first run terminates, more wheel detections are continued executing with, and the uniformity of maximum lattice site is tested；
S3：Whether succeeded according to assay judgement capture.
The present invention performs more wheel grid searches and peakvalue's checking in step s 2, by lattice site corresponding to peak value
Consistency check realize capture result judgement.Moreover, in order to improve acquisition speed, after first run search is completed, subsequent searches
Scope phase place change situation will be set according to caused by Doppler, in the case where ensuring acquisition performance, be avoided as far as possible more
Remaining search.
In step s 2, quadrature demodulation refers to that generating orthogonal local carrier according to carrier doppler driving carrier wave NCO believes
Number.Signal resampling refers to drive phase accumulator generation signal resampling enable signal according to pseudo-code Doppler.Coherent integration
Refer to be required according to resource and capture time etc., time domain approach or frequency domain method can be used, serial structure can be used
Parallel organization can be used.Detection refers to that " envelope detection ", " quadratic detection " or " differential detection " can be used.
In the present embodiment, above-mentioned steps S2 detailed process is：
S201：Phase grid search scope is determined according to prior information and its uncertainty；The prior information includes star
Between wireless signal arrival time (TOA) and arrival rate (FOA)；
S202：Calculate the detection statistics result of lattice point in hunting zone；The detection statistics result, time-domain calculation can be used
Method or frequency-domain calculations method, result of calculation are receiving sequence and the coherent integration or non-coherent integration knot of local pseudo-code sequence
Fruit；
S203：Searching and detecting result maximum, and lattice site corresponding to record；
S204：Searched for if the first run, follow-up " consistency check " stage Search scope is determined according to lattice site, and return
S202, otherwise into S205；" consistency check " the stage Search scope refers to：With lattice corresponding to first run peakvalue's checking result
Follow-up " consistency check " stage Search scope is set centered on point position；
S205：Consistency check is performed to maximum lattice site；The consistency check is：For every wheel maximum pair
The lattice site situation of change answered is counted, and thinks front and rear " consistent " when change is less than 1 chip；
S206：If not completing one-hundred-percent inspection, S202 is returned, is searched for into next round.The one-hundred-percent inspection is completed to refer to
If continuous N time assay " consistent ", terminates to examine, otherwise, to the end of N takes turns detection, one-hundred-percent inspection is completed.
In step s3, it is described to confirm that capture result refers to：If continuous N time assay " consistent ", judgement are captured as
Work(, otherwise, it is determined that capture failure.
As shown in figure 4, in the present embodiment, above-mentioned steps S205 detailed process is：
S2051：According to the hunting zone of setting, the detection statistic of whole lattice points is calculated, calculating process, which includes being concerned with, to be accumulated
Point, detection and non-coherent integration；
S2052：Maximum is searched in whole detection statistics, and records lattice point phase
S2053：Compared with last round of maximum lattice point phase, defining conformance criteria isIt is if full
The foot criterion adds up to counter m, otherwise resets；
S2054：If m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see whether reach amount of testing, if
N reaches N, terminates acquisition procedure, and because m is not up to M, announces capture failure；
S2055：If n is not up to N, examined into next round.
In above process, conformance criteria refers to：Cause the factor of lattice point phase place change to include code during two-wheeled search
The change of phase and the fuzziness of systemic resolution, due to using half chip as phase search lattice point, therefore, system resolution
There is the fuzziness of a phase lattice point in rate, i.e., in theoryIn addition, limit code-phase during two-wheeled search
The maximum changing value of position is 1/4 chip, i.e. 1/2 phase lattice point so thatStill set up.Therefore, withAs conformance criteria, on condition that code phase change is less than 1/4 chip during two-wheeled is searched for, and this needs to lead to
Code Doppler effect correction is crossed to be ensured.
Hunting zone is：When the first round searches for, set according to time of arrival (toa) TOA constant interval, if TOA
Excursion is [TOA-△TOA,TOA+△TOA], wherein, TOARepresent TOA, △ T that estimation obtainsOARepresent estimated bias.It will catch
The initial time obtained is arranged to TOA+△TOA, then, time interval corresponding to phase search is [0,2 △ TOA].So set
Advantage is that signal has just reached when starting search, avoids the false-alarm under no signal condition.
After in the first round, search terminates, into " consistency check " stage, subsequent searches section will all be searched for according to the first run
Peak value lattice site set, center isCode phase during the determination of hunting zone needs to consider search changes,
And false-alarm probability can be influenceed, in the case of code Doppler effect correction is carried out, code phase change is only several chips.If setting is searched
Rope scope be less than current catching method search unit, can tune to it is identical with search unit, although there is unnecessary search,
Realize the multiplexing of searching resource, it is possible to increase resource utilization.
Further, to ensure that lattice site true value has definite uniformity, believed in acquisition procedure using Doppler's priori
Breath carries out first-order dynamic compensation to receiving sequence code phases；If the precision deficiency of first-order dynamic compensation, can use second order dynamic to mend
Repay, the prior information for compensation includes arriving signal frequency and frequency change rate.
As shown in figure 5, present invention code phase dynamic compensation function principle schematic in concrete application example, including single order
Code phase compensation control and the compensation control of second order code phase.First compensation phase control refers to：Added up first according to frequency control word
Digital phase value is obtained, in phase value across all moment, generates complete cycle pulse, the equispaced of pulse is corresponding with frequency control word
Cycle is consistent.The control signal that complete cycle pulse will be taken out as data choosing, realizes the resampling to data.In second order compensation control
Frequency modulation control word is added, on the basis of frequency control word, realizes the real-Time Tracking Control to frequency and its change, because
This, it is more accurate to be controlled than first compensation phase.
Phase expression formula during first compensation phase controls is：
In formula, E represents frequency control word,Initial phase value is represented, L represents phase accumulator length, frequency control word
Computational methods be：
In formula, [] represents to ask whole operation, fcRepresent the resampling frequency after compensation, fsRepresent data sampling clock frequency.
fcEqual to twice of reception spread-spectrum code rate, Doppler shift caused by speed of related movement between star is contained, can be according to arrival
Signal frequency FOA is calculated.
When using first compensation phase control, if the f that prior information is broughtcEvaluated error is △ fc, accordingly, FREQUENCY CONTROL
The deviation △ E of word are：
If sampling pulse number corresponding to two-wheeled search time interval is N, then, because of Doppler frequency difference during two-wheeled search
Caused phase place change is △ EN, works as satisfaction：
△EN>2L-1 (5)
That is code phase change may then cause more than 1/4 chipSo as to cause missing inspection.Such case is then
Need to control using second order compensation, although add resource consumption compared to first compensation phase, but the degree of accuracy is higher.
Phase expression formula during second order compensation controls is：
In formula, F represents frequency modulation control word, and M represents frequency accumulator length.The computational methods of frequency modulation control word are：
Wherein, f 'cEqual to receiving, contain between star caused by relative motion acceleration twice of spreading code percentage speed variation
Doppler frequency changes, and can be calculated according to arriving signal frequency change rate MOA.
Based on the above method of the present invention, the performance evaluation after specific implementation is the following aspects：
(1) false-alarm probability；
If total searching times are N, consistency check thresholding is M, that is, needs continuous N time to meet phase uniform condition, recognize
Determine acquisition success, first round search lattice point quantity is L0, the lattice point quantity that search is often taken turns into " consistency check " stage is Lc。
Generation false-alarm can be divided into two kinds of situations, the first be in the range of " consistency check " stage Search not include true code phase and
Its neighboring lattice points, conditional probability in this case are defined as Pf1。
If the probability that the wheel of cut-off n-th passes through " consistency check " is Pf1(n), can be described as (n-M+1, n-M+2, L, n)
Wheel meets conformance criteria, and the n-th-M wheels do not meet conformance criteria, met unanimously in the absence of continuous N wheel in 1 to n-M-1 wheel
Property criterion.Under these conditions, if n-M+1 does not meet conformance criteria, and other conditions are constant, then when n+1 wheels meet one
Cause property criterion, then end n+1 wheels and in addition, meeting conformance criteria when n-M takes turns, but not will be caused by " consistency check "
When cut-off n-M wheels pass through " consistency check ", if n+1 wheels meet conformance criteria, cut-off n+1 wheels will be by the way that " uniformity is examined
Test ".Obtained using above-mentioned relation：
In formula, C represents combination, works as n<During 2M, Pf1(n-M)=0, then, according to (8), obtain
Pf1(n+1)=Pf1(n)n≥M+1 (9)
Comprehensive (8), (9) obtain, Pf1(n+1)≤Pf1(n), then, false-alarm probability meets in the case of the first
Wherein, Pf1(M) M wheels are represented by by the probability of consistency check before representing：
Pf1(M+1) represent that the first round does not meet the subsequent M wheels of conformance criteria by the probability of consistency check, can represent
For：
Therefore,
The first situation can be described as, and the maximum phase of first run search and the difference of true code phase are more thanThat
, when confirming " consistency check " stage Search scope, true code phase and its neighboring lattice points are foreclosed.Occur the
A kind of probability of situation is：
In formula, Pd3Represent that any one in true code phase and its two neighboring lattice point turns into the probability of maximum.
Second of situation is to include true code phase and its neighboring lattice points, reference in the range of " consistency check " stage Search
Computational methods in the case of the first, conditional probability in this case meet：
And the probability that second of situation occurs can be expressed as：
So, total false-alarm probability is expressed as：
(17)
(2) acquisition probability；
Definition is P by the probability that " consistency check " and result are true code phase or neighboring lattice points by the n-th wheeld
(n), then total acquisition probability is：
Wherein, Pd(M) M wheels are represented by by the probability of consistency check before representing：
Pd(M)=(Pd3)M+1 (19)
So, total acquisition probability meets：
Pd>(Pd3)M+1 (20)
(3) capture time；
If first run search time is T0, the search time often taken turns into " consistency check " stage is Tc, then, during capture
Between meet：
T0+MTc≤T0+nTc≤T0+NTc (21)
Parameter calculates：
Under normal circumstances, to improve acquisition probability under the conditions of certain false-alarm probability as criterion.First, determined according to (17) formula
N,M,Lc,Pd3.In the parameter combination for meeting false-alarm probability condition, less M and larger P is chosen as far as possibled3, so as to
Improve acquisition probability.However, Pd3It is relevant with the factor such as the time of integration, input signal carrier-to-noise ratio and noise profile characteristic, it is larger
Pd3Need to increase post detection integration or reduce acquisition sensitivity, if acquisition sensitivity is certain, in minimum carrier-to-noise ratio
Under the conditions of, according to noise actual distribution characteristic to Pd3Monte Carlo simulation is carried out, and finally determines post detection integration.
In a concrete application example, requiring false-alarm probability less than 10-5And acquisition probability is more than 95% condition
Lower progress parameter selection.Take N=4, M=3, Lc=60, Pd3=0.992, then, according to (17), false-alarm probability meets Pf<2×
10-6, according to (20), acquisition probability Pd>96.8%.
The above is only the preferred embodiment of the present invention, protection scope of the present invention is not limited merely to above-described embodiment,
All technical schemes belonged under thinking of the present invention belong to protection scope of the present invention.It should be pointed out that for the art
For those of ordinary skill, some improvements and modifications without departing from the principles of the present invention, the protection of the present invention should be regarded as
Scope.
Claims (10)
- A kind of 1. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, it is characterised in that step is：S1：Digital quadrature demodulation is carried out to reception signal, the complex baseband letter of I Q two-way composition is obtained after LPF Number；S2：Resampling is carried out to complex base band signal, obtains detection system by coherent integration, detection and non-coherent integration afterwards Metering, the maximum searched in whole detection statistics, and record the lattice site of maximum；After the first run terminates, continue executing with More wheel detections, and the uniformity of maximum lattice site is tested；S3：Whether succeeded according to assay judgement capture.
- 2. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 1, its feature exist In in the step S2, execution takes turns grid searches and peakvalue's checking, that is, passes through one to lattice site corresponding to peak value Cause property, which is examined, realizes capture result judgement；After first run search is completed, subsequent searches scope will the phase according to caused by Doppler Situation of change is set.
- 3. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 1, its feature exist In in the step S2, the quadrature demodulation refers to that generating orthogonal local carrier according to carrier doppler driving carrier wave NCO believes Number；The signal resampling refers to drive phase accumulator generation signal resampling enable signal according to pseudo-code Doppler.
- 4. the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 1 or 2 or 3, its It is characterised by, the detailed process of the step S2 is：S201：Phase grid search scope is determined according to prior information and its uncertainty；The prior information includes nothing between star Line time of arrival (toa) TOA and arrival rate FOA；S202：Calculate the detection statistics result of lattice point in hunting zone；S203：Searching and detecting result maximum, and lattice site corresponding to record；S204：Searched for if the first run, follow-up " consistency check " stage Search scope is determined according to lattice site, and return S202, otherwise into S205；S205：Consistency check is performed to maximum lattice site；The consistency check is：For corresponding to every wheel maximum Lattice site situation of change is counted, and thinks front and rear " consistent " when change is less than 1 chip；S206：If not completing one-hundred-percent inspection, S202 is returned, is searched for into next round.
- 5. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 4, its feature exist In in step s3, described to confirm that capture result refers to：If continuous N time assay " consistent ", judges acquisition success, no Then, capture failure is judged.
- 6. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 4, its feature exist In detection statistics result described in the step S202, using time-domain calculation method or frequency-domain calculations method, result of calculation is to connect Receive sequence and the coherent integration or non-coherent integration results of local pseudo-code sequence.
- 7. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 4, its feature exist In " consistency check " stage Search scope refers to described in the step S204：With lattice corresponding to first run peakvalue's checking result Follow-up " consistency check " stage Search scope is set centered on point position.
- 8. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 4, its feature exist In, if one-hundred-percent inspection described in the step S206 is completed to refer to continuous N time assay " consistent ", terminate to examine, otherwise, To the end of N takes turns detection, one-hundred-percent inspection is completed.
- 9. aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 4, its feature exist In the detailed process of the step S205 is：S2051：According to the hunting zone of setting, the detection statistic of whole lattice points is calculated, calculating process includes coherent integration, inspection Ripple and non-coherent integration；S2052：Maximum is searched in whole detection statistics, and records lattice point phaseS2053：Compared with last round of maximum lattice point phase, defining conformance criteria isShould if meeting Criterion adds up to counter m, otherwise resets；S2054：If m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see whether reach amount of testing, if n reaches To N, terminate acquisition procedure, and because m is not up to M, announce capture failure；S2055：If n is not up to N, examined into next round.
- 10. the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method according to claim 1 or 2 or 3, its It is characterised by, first-order dynamic compensation is carried out to receiving sequence code phases using Doppler's prior information in acquisition procedure；If one The precision deficiency of rank dynamic compensation, is further compensated using second order dynamic, and the prior information for compensation includes arriving signal frequency Rate and frequency change rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610618004.5A CN106291612B (en) | 2016-07-29 | 2016-07-29 | A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610618004.5A CN106291612B (en) | 2016-07-29 | 2016-07-29 | A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106291612A CN106291612A (en) | 2017-01-04 |
| CN106291612B true CN106291612B (en) | 2018-04-06 |
Family
ID=57663614
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610618004.5A Active CN106291612B (en) | 2016-07-29 | 2016-07-29 | A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106291612B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106908810B (en) * | 2017-01-13 | 2019-08-09 | 北京空间飞行器总体设计部 | It is a kind of how long code complex navigation signal phase consistency calibration method |
| CN107820311A (en) * | 2017-12-05 | 2018-03-20 | 中国科学院云南天文台 | For the quick early warning system of solar radio burst interference navigational communications event |
| CN109470884B (en) * | 2018-10-08 | 2020-08-07 | 浙江大学 | Inter-satellite relative angle measurement system and method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101718874A (en) * | 2009-12-14 | 2010-06-02 | 浙江大学 | Method and device for fast capturing signal in high-sensitivity navigational satellite receiver |
| CN102262233A (en) * | 2011-07-15 | 2011-11-30 | 王雪 | Carrier tracking method and carrier tracking loop of GNSS (Global Navigation Satellite System) receiver |
| CN102520423A (en) * | 2011-12-06 | 2012-06-27 | 成都金本华科技有限公司 | Circuit for rapidly capturing long period pseudo random spread spectrum code of satellite navigation receiver and capturing method thereof |
-
2016
- 2016-07-29 CN CN201610618004.5A patent/CN106291612B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101718874A (en) * | 2009-12-14 | 2010-06-02 | 浙江大学 | Method and device for fast capturing signal in high-sensitivity navigational satellite receiver |
| CN102262233A (en) * | 2011-07-15 | 2011-11-30 | 王雪 | Carrier tracking method and carrier tracking loop of GNSS (Global Navigation Satellite System) receiver |
| CN102520423A (en) * | 2011-12-06 | 2012-06-27 | 成都金本华科技有限公司 | Circuit for rapidly capturing long period pseudo random spread spectrum code of satellite navigation receiver and capturing method thereof |
Non-Patent Citations (1)
| Title |
|---|
| 导航星座星间链路信号自适应捕获方法;李献斌等;《系统工程与电子技术》;20150131;第37卷(第1期);正文第12-14页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106291612A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| DE10084224B4 (en) | Method for determining position from GPS signals | |
| CN101099089B (en) | Method and apparatus for increasing coherent integration length while receiving positioning signal | |
| CN102033236B (en) | Position and speed combined estimation method for satellite navigation | |
| CN104714244A (en) | Multi-system dynamic PPP resolving method based on robust self-adaption Kalman smoothing | |
| JP2011247758A (en) | Position calculation method and position calculation device | |
| CN1987515B (en) | Method and system for realizing GPS positioning under weak signal environment | |
| CN104459740A (en) | High-precision position differential positioning method of positioning terminal | |
| CN104536016A (en) | GNSS new-system signal capturing device and method | |
| CN106291612B (en) | A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method | |
| CN103592662A (en) | Carrier wave tracking method and loop for GPS signal receiver | |
| CN106199666A (en) | A kind of positioning and tracing method forwarding GNSS signal based on terminal | |
| DE102015103605B4 (en) | Distance measuring method and apparatus | |
| CN102565825B (en) | Received signal strength fiduciary level decision maker, method and code phase error calculation method | |
| CN101435866A (en) | Carrier wave tracking system and method of GPS coarse / capturing code signal | |
| CN113281786A (en) | Full-airspace multistage trigger type GNSS space signal quality monitoring and evaluating method | |
| CN102590834A (en) | Method of determining adequacy and adequacy determining device | |
| DE102005010782A1 (en) | Method and device for processing GPS signals and GPS receivers | |
| O’Driscoll et al. | Impact of extended coherent integration times on weak signal RTK in an ultra-tight receiver | |
| Di Lecce et al. | Neural technologies for increasing the GPS position accuracy | |
| CN103339526A (en) | Apparatus and method | |
| CN102162853B (en) | Reception signal integrating method and receiver | |
| Xue et al. | A new adaptive weighted fusion algorithm for BDS and UWB positioning system | |
| CN111060939B (en) | Method and device for optimizing density of virtual grid and positioning system | |
| CN104199062B (en) | Antenna array carrier wave combined tracking method based on Kinematic Decomposition | |
| CN114900247B (en) | Intelligent testing device for double pseudo code bidirectional comparison |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |