CN106291612A - A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method - Google Patents

Abstract

The invention discloses a kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, the steps include: that S1: docking is collected mail and number carry out digital quadrature demodulation, obtain after low-pass filtering I the complex base band signal that constitutes of Q two-way；S2: complex base band signal carries out resampling, obtains detection statistic through coherent integration, detection and non-coherent integration afterwards, searches for the maximum in whole detection statistic, and records the lattice site of maximum；After the first run has detected, adjust hunting zone, perform to take turns detection more, and the concordance of maximum lattice site is tested；S3: the most successful according to assay judgement capture.The present invention has time-varying characteristics between adaptation star, improves acquisition sensitivity and realize the advantages such as fast Acquisition.

Description

A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method

Technical field

Present invention relates generally to radio communication and field of measuring technique, refer in particular to one and be applicable 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 relatively High satellite orbit determines that precision and clock correction determine that precision is to ensure that navigation receives terminal and obtains the location of big system requirements or award The key of Shi Jingdu.

The demand of H_2O maser function has been expedited the emergence of out chain between navigation constellation star by aeronautical satellite precise orbit determination and time synchronized The concept on road.Setting up contact between aeronautical satellite upon inter-satellite link, the space segment of navigation system will be no longer isolated defending The combination of star, becomes a most collaborative entirety.By the cooperation of satellite-ground link, the control section of whole navigation system and space Duan Zhenzheng defines round-the-clock a, Seamless Network for round-the-clock, this be navigation system service operation management provide huge Performance space.Inter-satellite link can realize satellite navigation system in the case of only configuring minority monitoring station by inter-satellite link Accurate measurement obtains the metrical information of other segmental arc on track, thus reaches to obtain the ability of Precise Orbit parameter.

Navigation constellation inter-satellite link network is the satellite network that a class is complicated, has the acentric spy of typical flattening Levy, be a wireless network with some peer node.Navigation constellation inter-satellite link network needs real present time about Multi-multipoint measurement under the conditions of bundle, for completing accurate measurement function, navigational satellite system is broadcast by inter-satellite link and connects Receiving spread spectrum distance measuring signal and carry out the precise distance measurement between satellite, the capture of spread-spectrum signal is a pseudo-code and carrier wave two-dimensional search Process, the size of hunting zone directly determines the speed degree of signal capture and realizes 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 relatively big, gives and measures signal band Carry out bigger Doppler frequency shift, thus increase the capture range of signal, bring certain difficulty to capture, especially for Resource-constrained situation is calculated on star.Between two satellites any in constellation set up measure communication link time, interstellar distance and Doppler's excursion is relatively big, it addition, from measurement performance and the angle of safety, between star, wireless signal often uses length Cycle spread spectrum code, if not utilizing the two-dimensional search lattice point enormous amount of prior information, code phase and Doppler, gives capture real Now bring great difficulty.

For navigation constellation, time synchronized and precise orbit determination are the bases that navigation system is run, the satellite in constellation It is in a high-precision space-time datum.Inter-satellite link system utilizes ephemeris and the clock correction information prediction letter of link setup satellite Number time of advent and Doppler, it is possible to by delay time search scope control corresponding for code phase within ten musec order, and many The hunting zone of general Le is less than 100Hz.The reduction of capture range to a certain degree reduce capture realize difficulty, but still face with Lower several respects difficulty:

(1) from system application angle, generally require and be switched fast link setup object, cause signal duration between star It is short, accordingly, it would be desirable to receiver is rapidly completed capture.

(2) interstellar distance is remote, it is desirable to acquisition sensitivity is high.

(3) satellite borne equipment processes resource anxiety, needs fully to improve resource utilization.

Therefore, the capture of inter-satellite link wireless signal, substantially in certain hunting zone, utilize limited spaceborne set Standby resource realizes weak signal fast Acquisition.

Classical signal acquisition procedure includes detection statistic generation, peakvalue's checking and threshold judgement, and choosing of thresholding is straight Connecing and have influence on acquisition performance, the too low false-alarm probability that can increase of thresholding, thresholding selects too high, can reduce acquisition probability；Therefore, The determination of thresholding is most important.Under the conditions of CFAR criterion, decision threshold is relevant with noise probability distribution and statistical value thereof.Star Between link for being switched fast the demand of link setup object, cause antenna point to saltus step, cause the time-varying characteristics of background noise, make Obtain each link setup and be required for statistical noise distribution situation again, but, under the conditions of fast capturing signal, less number can only be relied on Carry out noise statistics according to sample, easily cause the deviation of statistical result, cause thresholding misalignment, thus have influence on acquisition performance.Separately Outward, the setting of thresholding improves the energy scale of detection statistic, can reduce acquisition sensitivity, is unfavorable for that weak signal captures.? The realization of whole catching method also needs to the resource problem in view of satellite borne equipment, optimizes and improve the utilization of resources of catching method Rate.

Summary of the invention

The technical problem to be solved in the present invention is that the technical problem existed for prior art, and the present invention provides one Kind can improve acquisition sensitivity, be suitable for weak signal capture, be capable of the aeronautical satellite inter-satellite link wireless communication of fast Acquisition Number high-performance prize judgment method.

For solve above-mentioned technical problem, the present invention by the following technical solutions:

A kind of aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, the steps include:

S1: docking is collected mail and number is carried out digital quadrature demodulation, obtain after low-pass filtering I the complex baseband that constitutes of Q two-way Signal；

S2: complex base band signal is carried out resampling, is examined through coherent integration, detection and non-coherent integration afterwards Survey statistic, search for the maximum in whole detection statistic, and record the lattice site of maximum；After the first run terminates, continue Perform to take turns detection more, and the concordance of maximum lattice site is tested；

S3: the most successful according to assay judgement capture.

As a further improvement on the present invention: in described step S2, many wheel grid searches and peakvalue's checking are performed, i.e. Adjudicate by the consistency check of lattice site corresponding to peak value is realized capture result；After completing first run search, follow-up search The phase place situation of change caused according to Doppler is set by rope scope.

As a further improvement on the present invention: in described step S2, described quadrature demodulation refers to drive according to carrier doppler Dynamic load ripple NCO generates orthogonal local carrier signal；Described signal resampling refers to drive phase-accumulated according to pseudo-code Doppler Device generates signal resampling and enables signal.

As a further improvement on the present invention: the detailed process of described step S2 is:

S201: determine phase place grid search scope according to prior information and uncertainty thereof；Described prior information includes star Between the wireless signal TOA time of advent and arrival rate FOA；

S202: the detection statistics result of lattice point in calculating hunting zone；

S203: searching and detecting result maximum, and record the lattice site of correspondence；

S204: if first run search, determine follow-up " consistency check " stage Search scope according to lattice site, and return S202, otherwise enters S205；

S205: maximum lattice site is performed consistency check；Described consistency check is: for often taking turns maximum pair The lattice site situation of change answered is added up, before and after change is i.e. thought less than 1 chip " unanimously "；

S206: if not completing one-hundred-percent inspection, returns S202, enters next round search.

As a further improvement on the present invention: in step s3, described confirmation capture result refers to: if continuous N time inspection Result " unanimously ", then judge acquisition success, otherwise, it is determined that capture unsuccessfully.

As a further improvement on the present invention: detection statistics result described in described step S202, time-domain calculation side is used Method or frequency-domain calculations method, result of calculation is receiving sequence and the coherent integration of local pseudo-code sequence or non-coherent integration results.

As a further improvement on the present invention: " consistency check " stage Search scope described in described step S204 is Refer to: centered by the lattice site that first run peakvalue's checking result is corresponding, set follow-up " consistency check " stage Search scope.

As a further improvement on the present invention: if one-hundred-percent inspection described in described step S206 completes to refer to continuous N time inspection Test result " unanimously ", then terminate inspection, otherwise, treat that N wheel detection terminates, complete one-hundred-percent inspection.

As a further improvement on the present invention: the detailed process of described step S205 is:

S2051: according to the hunting zone set, calculate the detection statistic of whole lattice point, calculating process includes relevant long-pending Point, detection and non-coherent integration；

S2052: search for maximum in whole detection statistic, and record lattice point phase place

S2053: compare with last round of maximum lattice point phase place, definition conformance criteria isIf it is full Enumerator m is added up by this criterion of foot, otherwise resets；

S2054: if m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see and whether reach amount of testing, if N reaches N, terminates acquisition procedure, and owing to m is not up to M, announces to capture unsuccessfully；

S2055: if n is not up to N, enters next round inspection.

As a further improvement on the present invention: utilize Doppler's prior information to receiving sequence code phases in acquisition procedure Carry out first-order dynamic compensation；If the precision that first-order dynamic compensates is not enough, use second order dynamic compensation further, for the elder generation compensated The information of testing includes arriving signal frequency and frequency change rate.

Compared with prior art, it is an advantage of the current invention that:

1, the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the present invention, it is to avoid decision threshold meter Calculating misalignment causes acquisition performance to decline, and is conducive to improving acquisition sensitivity, is suitable for weak signal capture, by adjusting hunting zone, Not only reduce resource consumption, and be capable of fast Acquisition.

2, the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the present invention, general in order to reduce false-alarm Rate, performs the search of taking turns of maximum more, and by code phase dynamic compensation, keeps the concordance of code phase true value, in order to enter one Step improves acquisition search speed, has readjusted the hunting zone in " consistency check " stage.

3, the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the present invention, it is to avoid traditional thresholding is sentenced The statistics misalignment certainly brought, calculate the problem such as resource consumption and acquisition sensitivity reduction, the most fully meet inter-satellite link body Making plays the fast capturing signal demand under noise characteristic time dependant conditions, it is also possible to be widely used in all types of spaceborne and ground expansion Frequency receiver.

Accompanying drawing explanation

Fig. 1 is the inter-satellite link wireless signal structural principle schematic diagram of range finding communicating integral.

Fig. 2 is present invention principle schematic of link radio signals capture between concrete application clock star.

Fig. 3 is the schematic flow sheet that the present invention scans for when concrete application.

Fig. 4 is the schematic flow sheet that the present invention carries out consistency check when concrete application.

Fig. 5 is the principle schematic that the present invention carries out code phase dynamic compensation function when concrete application.

Detailed description of the invention

Below with reference to Figure of description and specific embodiment, the present invention is described in further details.

As it is shown in figure 1, for range finding communicating integral inter-satellite link wireless signal structural representation, including measure branch road with Communication leg.Wherein, measuring branch road and do not contain data message, only include spreading code, communication leg then contains data message, and passes through Band spectrum modulation, two branch roads generate wireless signal between star after modulating carrier wave according to UQPSK mode.Signal according to this structural generation Model is:

In formula:

J: represent satellite number；

A_c: represent the range finding spread spectrum code amplitude being modulated to each frequency carrier wave I branch road；

A_p: represent the communications spread code amplitude being modulated to each frequency carrier wave Q branch road；

C: represent I branch road range finding spread spectrum code；

P: represent Q branch road communications spread code；

D_p: represent the numeric data code of modulation on Q branch road communications spread code；

F: represent inter-satellite link carrier frequency；

Represent that inter-satellite link measures 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, it is achieved H_2O maser and Time synchronized function, and assisted communication leg data-transformation facility, therefore, receiver captures only for measuring branch road Process.

The central principle of the present invention is: calculate the detection statistics result of whole lattice point in the region of search, searches maximum And record its place lattice site；Then near the lattice point of original maximum place, testing result is repeatedly calculated also in interval Search maximum, records the lattice site that often wheel search is corresponding, eventually through the concordance of maximum place lattice site simultaneously Whether inspection determines acquisition success.

As shown in Figures 2 and 3, the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method of the present invention, its Concretely comprise the following steps:

S1: docking is collected mail and number is carried out digital quadrature demodulation, obtain after low-pass filtering I the complex baseband that constitutes of Q two-way Signal；

S2: signal carries out resampling, resampling rate is pseudo-bit rate twice；Afterwards through coherent integration, detection and Non-coherent integration obtains detection statistic, searches for the maximum in whole detection statistic, and records the lattice site of maximum； After the first run terminates, continue executing with and take turns detection more, and the concordance of maximum lattice site is tested；

S3: the most successful according to assay judgement capture.

The present invention performs to take turns grid search and peakvalue's checking in step s 2, by lattice site corresponding to peak value Consistency check realize capture result judgement.And, in order to improve acquisition speed, after completing first run search, subsequent searches The phase place situation of change caused according to Doppler is set by scope, in the case of guaranteeing acquisition performance, avoids many as far as possible Remaining search.

In step s 2, quadrature demodulation refers to drive carrier wave NCO to generate orthogonal local carrier letter according to carrier doppler Number.Signal resampling refers to that driving phase accumulator to generate signal resampling according to pseudo-code Doppler enables signal.Coherent integration Refer to, according to the requirement such as resource and capture time, time domain approach or frequency domain method to be used, serial structure can be used also Parallel organization can be used.Detection refers to use " envelope detection ", " quadratic detection " or " differential detection ".

In the present embodiment, the detailed process of above-mentioned steps S2 is:

S201: determine phase place grid search scope according to prior information and uncertainty thereof；Described prior information includes star Between the wireless signal time of advent (TOA) and arrival rate (FOA)；

S202: the detection statistics result of lattice point in calculating hunting zone；Described detection statistics result, can use time-domain calculation Method or frequency-domain calculations method, result of calculation is receiving sequence and the coherent integration of local pseudo-code sequence or non-coherent integration knot Really；

S203: searching and detecting result maximum, and record the lattice site of correspondence；

S204: if first run search, determine follow-up " consistency check " stage Search scope according to lattice site, and return S202, otherwise enters S205；Described " consistency check " stage Search scope refers to: with the lattice that first run peakvalue's checking result is corresponding Point sets follow-up " consistency check " stage Search scope centered by position；

S205: maximum lattice site is performed consistency check；Described consistency check is: for often taking turns maximum pair The lattice site situation of change answered is added up, before and after change is i.e. thought less than 1 chip " unanimously "；

S206: if not completing one-hundred-percent inspection, returns S202, enters next round search.Described one-hundred-percent inspection completes to refer to If continuous N time assay " unanimously ", then terminate inspection, otherwise, treat that N wheel detection terminates, complete one-hundred-percent inspection.

In step s3, described confirmation capture result refers to: if continuous N time assay " unanimously ", then judge to be captured as Merit, otherwise, it is determined that capture unsuccessfully.

As shown in Figure 4, in the present embodiment, the detailed process of above-mentioned steps S205 is:

S2051: according to the hunting zone set, calculate the detection statistic of whole lattice point, calculating process includes relevant long-pending Point, detection and non-coherent integration；

S2052: search for maximum in whole detection statistic, and record lattice point phase place

S2053: compare with last round of maximum lattice point phase place, definition conformance criteria isIf it is full Enumerator m is added up by this criterion of foot, otherwise resets；

S2054: if m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see and whether reach amount of testing, if N reaches N, terminates acquisition procedure, and owing to m is not up to M, announces to capture unsuccessfully；

S2055: if n is not up to N, enters next round inspection.

In above process, conformance criteria refers to: the factor causing lattice point phase place to change includes code during two-wheeled search The change of phase place and the fuzziness of systemic resolution, owing to using half chip as phase search lattice point, therefore, system resolution There is the fuzziness of a phase place lattice point in rate, the most in theoryIt addition, limit code-phase during two-wheeled search The maximum changing value of position is 1/4 chip, i.e. 1/2 phase place lattice point so thatStill set up.Therefore, withAs conformance criteria, premise is that during two-wheeled is searched for, code phase change is less than 1/4 chip, and this needs logical Cross code Doppler effect correction to be ensured.

Hunting zone is: during first round search, be set according to the constant interval of time of arrival (toa) TOA, if TOA Excursion is [T_OA-△T_OA,T_OA+△T_OA], wherein, T_OARepresent TOA, the △ T estimating to obtain_OARepresent estimated bias.To catch The initial time obtained is set to T_OA+△T_OA, then, the time interval that phase search is corresponding is [0,2 △ T_OA].So arrange Advantage is, when starting search, signal arrives the most, it is to avoid the false-alarm under no signal condition.

After search terminates in the first round, entering " consistency check " stage, subsequent searches interval all will be searched for according to the first run Peak value lattice site be set, center isThe determination of hunting zone needs the code phase change considered during search, And can affect false-alarm probability, in the case of carrying out yard Doppler effect correction, code phase change is only several chips.If set searches Rope scope, less than the search unit of current catching method, can tune to identical with search unit, although occur in that unnecessary search, but Achieve the multiplexing of searching resource, it is possible to increase resource utilization.

Further, for guaranteeing that lattice site true value has definite concordance, acquisition procedure utilize Doppler's priori believe Cease and reception sequence code phases is carried out first-order dynamic compensation；If the precision that first-order dynamic compensates is not enough, second order can be used dynamically to mend Repay, include arriving signal frequency and frequency change rate for the prior information compensated.

As it is shown in figure 5, present invention code phase dynamic compensation function principle schematic in concrete application example, including single order Code phase compensates and controls and the compensation control of second order code phase.First compensation phase controls to refer to: first add up according to frequency control word Obtaining digital phase value, at phase value across the moment in week, generate complete cycle pulse, the equispaced of pulse is corresponding with frequency control word Cycle is consistent.Complete cycle pulse will be selected, as data, the control signal taken out, it is achieved the resampling to data.During second order compensation controls Add frequency modulation control word, on the basis of frequency control word, it is achieved that to frequency and the real-Time Tracking Control of change thereof, because of This, more accurate than first compensation phase controls.

Phase place expression formula during first compensation phase controls is:

In formula, E represents frequency control word,Representing initial phase value, L represents phase accumulator length, frequency control word Computational methods be:

$E=\[\frac{f_c\·2^L}{f_s}\]---\left(3\right)$

In formula, whole operation, f are asked in [] expression_cRepresent the resampling frequency after compensating, f_sRepresent data sampling clock frequency. f_cEqual to the twice of reception spread-spectrum code rate, contain the Doppler shift that speed of related movement between star causes, can be according to arrival Signal frequency FOA is calculated.

When using first compensation phase to control, if the f that prior information is brought_cEstimation difference is △ f_c, accordingly, FREQUENCY CONTROL The deviation △ E of word is:

$\mathrm{\Δ}E=\[\frac{{\mathrm{\Δf}}_c\·2^L}{f_s}\]---\left(4\right)$

If the sampling pulse number of two-wheeled interval search time correspondence is N, then, because of Doppler frequency difference during two-wheeled search The phase place caused is changed to △ EN, when meeting:

△EN>2^L-1 (5)

I.e. code phase varies more than 1/4 chip, then may causeThus cause missing inspection.This situation is then Needing to use second order compensation to control, although comparing first compensation phase to add resource consumption, but accuracy is higher.

Phase place 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:

$F=\[\frac{f_c^{\′}\·2^M}{f_s}\]---\left(7\right)$

Wherein, f '_cEqual to receiving the twice spreading code percentage speed variation, contain what relative motion acceleration between star caused Doppler frequency changes, and can be calculated according to arriving signal frequency change rate MOA.

Said method based on the present invention, the performance evaluation after being embodied as is the following aspects:

(1) false-alarm probability；

If total searching times is N, consistency check thresholding is M, i.e. needs continuous N time to meet phase place uniform condition, recognizes Determining acquisition success, first round search lattice point quantity is L₀, it is L that " consistency check " stage of entrance often takes turns the lattice point quantity of search_c。 Occur false-alarm can be divided into two kinds of situations, the first be do not include in the range of " consistency check " stage Search true code phase and Its neighboring lattice points, conditional probability in this case is defined as P_f1。

If it is P that the probability by " consistency check " is taken turns in cut-off n-th_f1(n), can be described as (n-M+1, n-M+2, L, n) Wheel all meets conformance criteria, and the n-th-M wheel does not meets conformance criteria, there is not continuous N wheel and meet consistent in 1 to n-M-1 wheel Property criterion.Under these conditions, if n-M+1 does not meets conformance criteria, and other conditions are constant, then when n+1 wheel meets one Cause property criterion, then cut-off n+1 wheel, it addition, meet conformance criteria when n-M takes turns but will not resulted in by " consistency check " When cut-off n-M wheel is by " consistency check ", if n+1 wheel meets conformance criteria, cut-off n+1 wheel will be by " concordance is examined Test ".Above-mentioned relation is utilized to obtain:

$\begin{array}{c}{P}_{f1}(n+1)=\left(\frac{C_{L_c-3}^1}{C_3^1}{P}_{f1}\right(n\left)\frac{C_3^1}{C_{L_c}^1}\right)+\left(\frac{C_3^1}{C_{L_c-3}^1}\frac{C_{L_c-3}^1}{C_3^1}{P}_{f1}\right(n\left)\frac{C_3^1}{C_{L_c}^1}\right)(1-P_{f1}(n-M\left)\right)\\ =\left(\frac{C_{L_c-3}^1}{C_{L_c}^1}{P}_{f1}\right(n\left)\right)+\left(\frac{C_3^1}{C_{L_c}^1}{P}_{f1}\right(n\left)\right)(1-P_{f1}(n-m\left)\right)\end{array}---\left(8\right)$

In formula, C represents combination, when n < during 2M, P_f1(n-M)=0, then, according to (8), obtain

P_f1(n+1)=P_f1(n)n≥M+1 (9)

Comprehensively (8), (9) obtain, P_f1(n+1)≤P_f1(n), then, in the case of the first, false-alarm probability meets

$P_{f1}=\underset{n=M}{\overset{n=N}{\&Sigma;}}{P}_{f1}\left(n\right)\&le;P_{f1}\left(M\right)+(N-M)P_{f1}(M+1)---\left(10\right)$

Wherein, P_f1(M) the M wheel probability by consistency check before representing, is represented by:

$P_{f1}\left(M\right)={\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M---\left(11\right)$

P_f1(M+1) represent that the first round does not meets the conformance criteria M subsequently wheel probability by consistency check, can represent For:

$P_{f1}\left(M\right)=\frac{C_{L_c-3}^1}{C_{L_c}^1}{\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M---\left(12\right)$

Therefore,

$\begin{array}{c}{P}_{f1}\&le;{\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M+(N-M)\frac{C_{L_c-3}^1}{C_{L_c}^1}{\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M\\ <(N-M+1){\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M\end{array}---\left(13\right)$

The first situation can be described as, and the maximum phase place of first run search is more than with the difference of true code phaseThat , when confirming " consistency check " stage Search scope, true code phase and neighboring lattice points thereof are foreclosed.Occur the The probability of a kind of situation is:

$(1-P_{d3})\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}---\left(14\right)$

In formula, P_d3Represent in true code phase and adjacent two lattice points thereof that any one becomes the probability of maximum.

The second situation is to include true code phase and neighboring lattice points thereof, reference in the range of " consistency check " stage Search Computational methods in the case of the first, conditional probability in this case meets:

$P_{f2}<(N-M+1){(1-P_{d3})}^M{\left(\frac{C_3^1}{C_{L_c-3}^1}\right)}^M---\left(15\right)$

And occur the probability of the second situation to be expressed as:

$1-(1-P_{d3})\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}---\left(16\right)$

So, total false-alarm probability is expressed as:

$\begin{array}{c}{P}_f=(1-P_{d3})\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}{P}_{f1}+(1-(1-P_{d3}\left)\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}\right)P_{f2}\\ <(N-M+1)\left(\right(1-P_{d3})\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}{\left(\frac{C_3^1}{C_{L_c}^1}\right)}^M+(1-\left(1-P_{d3}\right)\frac{C_{L_0-3-L_c}^1}{C_{L_0-3}^1}\left){\left(1-P_{d3}\right)}^M{\left(\frac{C_3^1}{C_{L_c-3}^1}\right)}^M\right)\\ <(N-M+1)(1-P_{d3}){\left(\frac{C_3^1}{C_{L_c-3}^1}\right)}^M\end{array}$

(2) acquisition probability；

Define that to take turns by " consistency check " and probability that result is true code phase or neighboring lattice points by n-th be P_d (n), then total acquisition probability is:

$P_d=\underset{n=M}{\overset{n=N}{\&Sigma;}}{P}_d\left(n\right)>P_d\left(M\right)---\left(18\right)$

Wherein, P_d(M) the M wheel probability by consistency check before representing, is represented by:

P_d(M)=(P_d3)^M+1 (19)

So, total acquisition probability meets:

P_d>(P_d3)^M+1 (20)

(3) capture time；

If first run search time being T₀, the search time that " consistency check " stage of entrance often takes turns is T_c, then, during capture Between meet:

T₀+MT_c≤T₀+nT_c≤T₀+NT_c (21)

Parameter calculates:

Under normal circumstances, under the conditions of certain false-alarm probability, acquisition probability is improved as criterion.First, determine according to (17) formula N,M,L_c,P_d3.In the parameter meeting false-alarm probability condition combines, choose less M and bigger P as far as possible_d3Such that it is able to Improve acquisition probability.But, P_d3Relevant, more greatly with factors such as the time of integration, input signal carrier-to-noise ratio and noise profile characteristics P_d3Need increase post detection integration or reduce acquisition sensitivity, if acquisition sensitivity is certain, then in minimum carrier-to-noise ratio Under the conditions of, according to noise actual distribution characteristic to P_d3Carry out Monte Carlo simulation, and finally determine post detection integration.

In a concrete application example, requiring that false-alarm probability is less than 10^-5And the condition that acquisition probability is more than 95% Under carry out parameter selection.Take N=4, M=3, L_c=60, P_d3=0.992, then, according to (17), false-alarm probability meets P_f<2× 10^-6, according to (20), acquisition probability P_d> 96.8%.

Below being 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, should be regarded as the protection of the present invention Scope.

Claims (10)

1. an aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method, it is characterised in that step is:

S1: docking is collected mail and number is carried out digital quadrature demodulation, obtain after low-pass filtering I the complex baseband letter that constitutes of Q two-way Number；

S2: complex base band signal is carried out resampling, obtains detection system through coherent integration, detection and non-coherent integration afterwards Metering, searches for the maximum in whole detection statistic, and records the lattice site of maximum；After the first run terminates, continue executing with Take turns detection more, and the concordance of maximum lattice site is tested；

S3: the most successful according to assay judgement capture.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 1, its feature exists In, in described step S2, perform many wheel grid searches and peakvalue's checking, i.e. by lattice site corresponding to peak value The inspection of cause property realizes capture result judgement；After completing first run search, the phase place that subsequent searches scope will cause according to Doppler Situation of change is set.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 1, its feature exists In, in described step S2, described quadrature demodulation refers to drive carrier wave NCO to generate orthogonal local carrier letter according to carrier doppler Number；Described signal resampling refers to that driving phase accumulator to generate signal resampling according to pseudo-code Doppler enables signal.

4. according to the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method described in claim 1 or 2 or 3, its Being characterised by, the detailed process of described step S2 is:

S201: determine phase place grid search scope according to prior information and uncertainty thereof；Described prior information includes nothing between star Line time of arrival (toa) TOA and arrival rate FOA；

S202: the detection statistics result of lattice point in calculating hunting zone；

S203: searching and detecting result maximum, and record the lattice site of correspondence；

S204: if first run search, determine follow-up " consistency check " stage Search scope according to lattice site, and return S202, otherwise enters S205；

S205: maximum lattice site is performed consistency check；Described consistency check is: corresponding for often taking turns maximum Lattice site situation of change is added up, before and after change is i.e. thought less than 1 chip " unanimously "；

S206: if not completing one-hundred-percent inspection, returns S202, enters next round search.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 4, its feature exists In, in step s3, described confirmation capture result refers to: if continuous N time assay " unanimously ", then judge acquisition success, no Then, it is determined that capture unsuccessfully.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 4, its feature exists In, detection statistics result described in described step S202, use time-domain calculation method or frequency-domain calculations method, result of calculation is for connecing Receive sequence and the coherent integration of local pseudo-code sequence or non-coherent integration results.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 4, its feature exists In, " consistency check " stage Search scope described in described step S204 refers to: with the lattice that first run peakvalue's checking result is corresponding Point sets follow-up " consistency check " stage Search scope centered by position.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 4, its feature exists In, if one-hundred-percent inspection described in described step S206 completes to refer to continuous N time assay " unanimously ", then terminate to check, otherwise, Treat that N wheel detection terminates, complete one-hundred-percent inspection.

Aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method the most according to claim 4, its feature exists In, the detailed process of described step S205 is:

S2051: according to the hunting zone set, calculating the detection statistic of whole lattice point, calculating process includes coherent integration, inspection Ripple and non-coherent integration；

S2052: search for maximum in whole detection statistic, and record lattice point phase place

S2053: compare with last round of maximum lattice point phase place, definition conformance criteria isShould if meeting Enumerator m is added up by criterion, otherwise resets；

S2054: if m reaches M, terminate acquisition procedure, and confirm acquisition success, otherwise, see and whether reach amount of testing, if n reaches To N, terminate acquisition procedure, and owing to m is not up to M, announce to capture unsuccessfully；

S2055: if n is not up to N, enters next round inspection.

10. according to the aeronautical satellite inter-satellite link wireless signal high-performance prize judgment method described in claim 1 or 2 or 3, its It is characterised by, utilizes Doppler's prior information to carry out first-order dynamic compensation to receiving sequence code phases in acquisition procedure；If one The precision of rank dynamic compensation is not enough, uses second order dynamic compensation, the prior information for compensating to include arriving signal frequency further Rate and frequency change rate.