CN106772476A - Multivariable space-time chaos satellite navigation high-performance ranging code method and system - Google Patents
Multivariable space-time chaos satellite navigation high-performance ranging code method and system Download PDFInfo
- Publication number
- CN106772476A CN106772476A CN201611170030.2A CN201611170030A CN106772476A CN 106772476 A CN106772476 A CN 106772476A CN 201611170030 A CN201611170030 A CN 201611170030A CN 106772476 A CN106772476 A CN 106772476A
- Authority
- CN
- China
- Prior art keywords
- function
- real
- imaginary part
- pseudo
- random number
- 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.)
- Granted
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
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
Abstract
The invention discloses a kind of multivariable space-time chaos satellite navigation high-performance ranging code method and system, the present invention constructs multivariable vector pseudo-noise code generator, first, the multigroup nonlinear function containing different parameters value is respectively adopted to the current state value of the current location and deviation post of multivariable plural vector real part and each component of imaginary part to act on, added by weights of diffusion coefficient respectively, subtract, multiply, or the hybrid operation for removing, produced with the plural pseudo-random number sequence of Annual distribution by state iteration, secondly from real part and imaginary part correlated components tap extraction real number pseudo-random number sequence, binaryzation and mould two and obtain required high-performance ranging code.The composite can be widely applied to satellite navigation system, it can also be used to the aspect such as various range-measurement systems, communication system, Broadcast and TV system, control system.
Description
Technical field
The present invention relates to technical field of satellite navigation, specially multivariable space-time chaos satellite navigation high-performance ranging code reality
Existing method and system.
Background technology
Four large satellite navigation system have GPS (the Global Positioning System) satellite in the U.S. in the world at present
Navigation system, GLONASS (Global Navigation Satellite System) satellite navigation system, the European Union of Russia
Galileo satellite navigation system and China the Big Dipper (BeiDou) satellite navigation system.Except GLONASS uses frequency division multiple access
Outside (FDMA, Frequency Division Multiple Access) communication mode, remaining satellite navigation system uses code
Divide multiple access (CDMA, Code Division Multiple Access) communication mode.The ranging code that they are used is divided into civilian
Bigness scale code and the military class of accurate code two, target coarse positioning can be carried out using the satellite navigation system of bigness scale code, and uses accurate code
Satellite navigation system can carry out target with high precision positioning.
The ranging code that present satellites navigation system is used is the binary system pseudo noise code that simple primary is used using linear approach
It is initial to register that one or several certain series linear feedback shift registers of generator first pass through one section of short binary sequence
Change, then produced by register shift, the ranging code generally existing complexity that this method is produced is low, security is poor, code length is solid
It is fixed and shorter, the shortcomings of number of codes is limited, and shift register also needs to undertake linear feedback and Satellite Phase shares out the work, and
Make its complex structure.
Another kind is the real number pseudo-noise code generator discussed in existing document, and it is by nonlinear method, by space-time
Chaos unidirectional couplings Map Lattices model construction real number pseudo-noise code generator is first initial with same real number plaid matching point variable state
Change, the real number being distributed with time state is then produced under Nonlinear Mapping kinetic function (referred to as " nonlinear function ") effect
Pseudo random number, then the lattice point variant time state Distribution Value is obtained through related lattice point variable tap, and real number pseudo random number is exported,
Binaryzation and mould two and be ranging code.Nonlinear function typically takes f (x)=- δ x2+ 1, δ are nonlinear strength, and it is by non-linear
Intensity, variable, the power of variable and integer constant a composition, nonlinear strength, variable, the power of variable and constant term are referred to as non-thread
The parameter of property function.Real number pseudo-noise code generator can overcome the ranging code code that current binary system pseudo-noise code generator is produced
The fixed and shorter and limited shortcoming of number of codes long, but due to using a kind of nonlinear function and generally quadratic power function work
With wherein nonlinear function is made up of nonlinear strength, variable, integer constant, and real number precision is 10-2, answer the ranging code of generation
The not high, security of miscellaneous degree is not strong, and excessive using space lattice number, and each lattice point state variable once only produces a kind of real number pseudo-
Random number.
Additionally, existing patent discusses plural pseudo-noise code generator and vector pseudo-noise code generator produces satellite navigation
System high-performance ranging code realizes technology, and they are all to use nonlinear method to contain difference by being built based on single argument
The function of parameter, compound function or variable, using multigroup nonlinear function act on state variable or component current location and
The current state value of deviation post, makes it produce the plural pseudo random number being distributed with time state, then related to imaginary part by real part
Variable or component tap output real number pseudo random number and mould two change to obtain aeronautical satellite ranging code, and the technology can be entered with comprehensively solve two
All technological deficiencies that the aeronautical satellite ranging code that pseudo-noise code generator processed and real number pseudo-noise code generator are produced is present, entirely
Face lifting satellite navigation system range finding code performance.But this technology is realized using based on single argument, it is not met by producing
The need for raw complexity satellite navigation system ranging code technology higher.
The content of the invention
In order to solve the above-mentioned technical problem, the invention provides it is a kind of can General Promotion aeronautical satellite ranging code complexity
Multivariable space-time chaos satellite navigation high-performance ranging code method and system.
The technical scheme that the method for the present invention is used is：A kind of multivariable space-time chaos satellite navigation high-performance ranging code
Implementation method, comprises the following steps：
Step 1：Build multivariable vector pseudo-noise code generator G；
Step 2：Judge whether to produce sin type ranging codes；
If so, then performing following step 5；
If it is not, then performing following step 3；
Step 3：Judge whether to produce cos type ranging codes；
If so, then performing following step 5；
If it is not, then performing following step 4；
Step 4：Produce sin types and cos type ranging codes；
Step 5：Extension plural component G0；
Step 6：Build multivariable nonlinearity function；
Step 7：In parameter initialization and the sequence of real numbers initialization G constituted using pseudo-random number sequence or by different real numbers
Effective plural component and G0The state value of middle extension plural component；
Step 8：Effective plural component in G is taken, is produced with the plural pseudo random number of Annual distribution by state iteration
Sequence；
Step 9：Whether effectively plural component takes in judging G；
If so, then performing following step 10；
If it is not, then changing extension component state value, and turn round execution step 8；
Step 10：Real part pseudo random number and imaginary part pseudo random number are taken respectively；
Step 11：By real part pseudo random number and imaginary part pseudo random number binaryzation and the He of mould two is carried out, obtain high-performance range finding
Code.
The technical scheme that device of the invention is used is：A kind of multivariable space-time chaos satellite navigation high-performance ranging code
System is realized, including multivariable vector pseudo-noise code generator builds module, the first judge module, expansion module, non-linear letter
Number builds module, initialization module, state iteration module, the second judge module, real part pseudo random number and imaginary part pseudo random number and carries
Modulus block, binarization block, the computing module of mould two；
The multivariable vector pseudo-noise code generator builds module, for building multivariable vector pseudo-noise code generator
G；
First judge module, for performing following judgements and operation；
Judge whether to produce sin type ranging codes, if then performing method and step 5, if otherwise performing method and step 3；Judge
Whether cos type ranging code is produced, if then performing method and step 5；If otherwise performing method and step 4；Produce sin types and cos types
Ranging code；
The expansion module, for extending plural component G0；
The nonlinear function builds module, for building multivariable nonlinearity function；
The initialization module, the real number constituted for parameter initialization and using pseudo-random number sequence or by different real numbers
Effective plural component and G in sequence initialization G0The state value of middle extension plural component；
The state iteration module, for taking in G effectively plural component, is produced with Annual distribution by state iteration
Plural pseudo-random number sequence；
Second judge module, for judging whether effectively plural component takes in G；If so, then performing method
Step 10；If it is not, then changing extension component state value, and turn round execution method and step 8；
The real part pseudo random number and imaginary part pseudo random number extraction module, it is pseudo- for taking real part pseudo random number and imaginary part respectively
Random number；
The binarization block, for by real part pseudo random number and imaginary part pseudo random number binaryzation；
The computing module of the mould two, for the real part pseudo random number and imaginary part pseudo random number after binaryzation to be carried out into mould two
With acquisition high-performance ranging code.
The present invention can overcome the ranging code that current binary system pseudo-noise code generator and real number pseudo-noise code generator are produced
Exist all technological deficiencies, while can General Promotion plural number pseudo-noise code generator and vector pseudo-noise code generator produce
The complexity of high-performance ranging code.
Compared to the prior art, the invention has the advantages that and beneficial effect：
(1) complexity pseudo noise code high can be obtained；
Using the different nonlinear function collective effects of multivariable, acted on by hybrid operation iteratively produce it is pseudo- with
Machine code, the pseudo noise code complexity for being obtained is high.
(2) pseudo noise code of high safety can be obtained；
The pseudo noise code that complexity can be produced high, so as to fully ensure that the high security of pseudo noise code.
(3) pseudo noise code code length is not limited by series；
The random code with Annual distribution for producing, its maximum code length is unrelated using series with pseudo-noise code generator, and can
Up to endless.
(4) series is special few；
By plural implementation, pseudo-noise code generator series can be to greatest extent reduced.
(5) produce pseudo noise code pattern many；
The parameter that pseudorandom pattern of the present invention is included by the real number of initialization plural component state value, nonlinear function
Such as diffusion coefficient, the parameters precision of nonlinear strength determine that receiving terminal these parameters precisions are at least 10-5, the present invention can produce
Raw pseudorandom pattern is at least 105×L×(nx+ny) individual, L is code length.
(6) pseudo noise code characteristic can at any time be adjusted；
When pseudo noise code pseudo-randomness is destroyed, can by adjust nonlinear strength, diffusion coefficient, fine setting function work
Working frequency, phase etc. are recovered.
Brief description of the drawings
Fig. 1 is the method idiographic flow schematic diagram that the present invention is implemented；
Fig. 2 is G in the embodiment of the present invention1The pseudo-random number sequence that 3rd intermediateness component is obtained using sin functions；
Fig. 3 is G in the embodiment of the present invention1The pseudo-random number sequence that 3rd intermediateness component is obtained using cos functions；
Fig. 4 is G in the embodiment of the present invention1The pseudo random number that 3rd intermediateness component is obtained using sin and cos functions
Sequence；
Fig. 5 is using the pseudorandom obtained after ranking method treatment in the embodiment of the present invention to pseudo-random number sequence shown in Fig. 2
Code；
Fig. 6 is using the pseudorandom obtained after ranking method treatment in the embodiment of the present invention to pseudo-random number sequence shown in Fig. 3
Code；
Fig. 7 is using the pseudorandom obtained after ranking method treatment in the embodiment of the present invention to pseudo-random number sequence shown in Fig. 4
Code；
Fig. 8 is the sin type high-performance ranging codes of acquisition in the embodiment of the present invention；
Fig. 9 is the cos type high-performance ranging codes of acquisition in the embodiment of the present invention；
Figure 10 is sin the and cos type high-performance ranging codes of acquisition in the embodiment of the present invention；
Figure 11 is the auto-correlation function of ranging code shown in Fig. 8 in the embodiment of the present invention；
Figure 12 is the auto-correlation function of ranging code shown in Fig. 9 in the embodiment of the present invention；
Figure 13 is the auto-correlation function of ranging code shown in Figure 10 in the embodiment of the present invention；
Figure 14 is the system structure diagram of the embodiment of the present invention.
Specific embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, below in conjunction with the accompanying drawings and embodiment is to this hair
It is bright to be described in further detail, it will be appreciated that implementation example described herein is merely to illustrate and explain the present invention, not
For limiting the present invention.
See a kind of multivariable space-time chaos satellite navigation high-performance ranging code implementation method that Fig. 1, the present invention are provided, use
In sin type high-performance ranging codes, including：
S1 builds multivariable vector pseudo-noise code generator G, G and is made up of single plural vector X+Yj, in If+N+|Ib
| its component of dimensional linear space isReferred to as multivariable plural component,WithRespectively a series of order arrangements and the state component for intercoupling respectively；I represents multivariable plural
The position number of component, i=1,2 ..., If+N+|Ib|, If, N be positive integer, IbIt is negative integer, IfTo spreading number, I before representingb
Represent backward spreading number；ixRepresent the sequence number of multivariable plural component real component x, iyRepresent multivariable plural point
The sequence number of amount imaginary y, ix=1,2 ..., nx,iy=1,2 ..., ny,nxAnd nyThe respectively quantity of state component x and y,
Preceding If×(nx+ny) and rear | Ib|×(nx+ny) individual multivariable plural component be designated as extend plural component, its position-order
Number difference i ∈ [1, If] and i ∈ [If+N+1,If+N+|Ib|]；Middle N × (nx+ny) individual multivariable plural component has been designated as
Effect plural component, its position number i ∈ [If+1,If+ N], effective plural componentWithRespectively constitute
Pseudo-noise code generator G1、G2；
S2 extension plural components constitute G0, extension plural number of components is more than or equal to default position offset
The maximum of absolute value；
S3 builds and is respectively acting on effective plural component real part and imaginary part current location and deviation post current state
The nonlinear function of value, specially：
Will act on the non-linear letter of the current state value of effective plural component real part current location and deviation post
Number is designated as real part nonlinear function respectively, will act on the current of effective plural component imaginary part current location and deviation post
The nonlinear function of state value is designated as imaginary part nonlinear function respectively；Wherein：
The structure of real part nonlinear function is specially：
Real part nonlinear functionBy k moment real part vectors ΧkNxIndividual componentN is built respectivelyx
The sin functions and n of the current location containing different parameters value of individual different powersxThe forward direction containing different parameters value of individual different powers
With the variable of reserve migration position；It is power with the nonlinear strength that each sin functions are born, takes IxIndividual sin functions weighted sum obtains reality
Portion's first function, to remaining nx-IxIndividual sin functions weighted sum obtains real part second function；With each variable bear it is non-linear
Intensity is power, to nxThe variable weighting of individual deviation post is sued for peace to obtain the variable item of real part first；Real part first function is divided by real part
Plus real part second function and the first real constant, gained multinomial is real part nonlinear function to first variable item；
Constructed real part nonlinear function is used for acting on effective plural component real part current location and bits of offset
The current state value put；
The structure of imaginary part nonlinear function is specially：
Imaginary part nonlinear functionBy k moment imaginary part vector YskNyIndividual componentN is built respectivelyyIt is individual
The sin functions and n of the current location containing different parameters value of different powersyThe forward direction containing different parameters value of individual different powers and
The variable of reserve migration position；It is power with the nonlinear strength that each sin functions are born, takes nyIndividual sin functions weighted sum obtains imaginary part
First function；It is power with the nonlinear strength that each variable is born, to IyThe mutually multiplied imaginary part first of variable weighting of individual deviation post
Variable item, to remaining ny-IyIndividual variable weighting is sued for peace to obtain the variable item of imaginary part second；The variable item of imaginary part first subtracts imaginary part second
Variable item is divided by, plus the second real constant, gained multinomial is imaginary part nonlinear function after imaginary part first function；
Constructed imaginary part nonlinear function is used for acting on effective plural component imaginary part current location and bits of offset
The current state value put；
Represent by k moment real part vectors ΧkThe serial number l of structure1Real part nonlinear function,Represent by k
Moment imaginary part vector YkThe serial number l of structure2Imaginary part nonlinear function,
nxRepresent the quantity of effective plural component x, nyRepresent the quantity of effective plural component y；IxFor non-linear
FunctionMiddle part number of components, it is no more than nxPositive integer, IyIt is nonlinear functionMiddle part number of components, it is little
In nyPositive integer, its value sets itself as needed；
Wherein, the parameter of function includes working frequency, the power of function, the range value of function, the phase of function, position-order
Number, position offset and state translational movement；The parameter of variable includes power, position number, position offset and the state of variable
Translational movement；Position offset is the amount that position number is increased or decreased, and state translational movement is that variable states value is increased or decreased
Amount；
In S4 parameter initializations and the sequence of real numbers initialization G constituted using pseudo-random number sequence or by different real numbers effectively
Plural component and G0The state value of middle extension plural component；
S5 is worked as to effective plural component real part current location and deviation post respectively using real part nonlinear function
Preceding state value is acted on, and obtains real part threshold；Using imaginary part nonlinear function to effective plural component imaginary part present bit
Put and acted on the current state value of deviation post, obtain imaginary part threshold；Based on diffusion coefficient, to real part threshold and imaginary part
Threshold is carried out adding respectively, subtracted, multiplication and division or comprising adding, subtract, in multiplication and division at least two computings hybrid operation, changed by state
In generation, is produced with the plural pseudo-random number sequence of Annual distribution；
This step is further included：
Real part state iteration, specially：
The current of effective plural component real part current location and deviation post is acted on using real part nonlinear function
State value, obtains real part Nonlinear function；
With diffusion coefficient as weights, to the k momentIndividual real part Nonlinear functionSummation is weighted, first is obtained
Real part threshold；RemainingIndividual real part Nonlinear functionMultiplication is weighted, the second real part threshold is obtained；Institute
There is real part Nonlinear functionArithmetic is sued for peace to obtain the 3rd real part threshold；Second real part threshold is made divided by the 3rd real part
With value subsequent time real part current position state value is obtained plus the first real part threshold；
Imaginary part state iteration, specially：
The current of effective plural component imaginary part current location and deviation post is acted on using imaginary part nonlinear function
State value, obtains imaginary part Nonlinear function；
With diffusion coefficient as weights, to the k momentIndividual imaginary part Nonlinear functionSummation is weighted, first is obtained
Imaginary part threshold；RemainingIndividual imaginary part Nonlinear functionMultiplication is weighted, the second imaginary part threshold is obtained；Institute
There is imaginary part Nonlinear functionArithmetic is sued for peace to obtain the 3rd imaginary part threshold；First imaginary part threshold subtracts the 3rd imaginary part work
With value subsequent time imaginary part current position state value is obtained divided by the second imaginary part threshold；
Plural pseudo-random number sequence or sequence of real numbers that S6 is obtained using current each effective plural component, change G0In
The state value of each extension plural component, or carry out reconfiguring arrangement each other using these state values changed；Then,
Next group of effective plural component in G is read, step S5 is performed to next group of effective plural component；Have when all in G
The effect equal completion status iteration of plural component, performs step S7；
S7 is respectively from G1And G2Middle correlated components tap extraction is designated as respectively with the real number pseudo-random number sequence of Annual distribution
The pseudo random number of real part first and the pseudo random number of imaginary part first；
S8 by each real number random number in the pseudo random number of real part first and the pseudo random number of imaginary part first in chronological order respectively with
A reference value compares, if being more than a reference value, the real number random number value 1, otherwise value 0, you can obtain binaryzation first is pseudo-
Random code and the second pseudo noise code；A reference value is respectively real number in the description pseudo random number of real part first and the pseudo random number of imaginary part first
The statistic of pseudo random number median size；
First pseudo noise code and the second pseudo noise code are carried out the He of mould two by S9, obtain the 3rd good pseudorandom of pseudo-randomness
Code, i.e. sin types high-performance ranging code；If the pseudo-randomness of sin type high-performance ranging codes is destroyed, then need to adjust non-thread
Property function diffusion coefficient, function and/or variable nonlinear strength, finely tune function working frequency, phase, finely tune a reference value
Again the pseudo random number of binaryzation real part first and the pseudo random number of imaginary part first.
See a kind of multivariable space-time chaos satellite navigation high-performance ranging code implementation method that Fig. 1, the present invention are provided, use
In cos type high-performance ranging codes, including：
S1~S2, with step S1~S2 in sin type high-performance ranging code methods；
S3 builds and is respectively acting on effective plural component real part and imaginary part current location and deviation post current state
The nonlinear function of value, specially：
Will act on the non-linear letter of the current state value of effective plural component real part current location and deviation post
Number is designated as real part nonlinear function respectively, will act on the current of effective plural component imaginary part current location and deviation post
The nonlinear function of state value is designated as imaginary part nonlinear function respectively；Wherein：
The structure of real part nonlinear function is specially：
Real part nonlinear functionBy k moment real part vectors ΧkNxIndividual componentN is built respectivelyx
The cos functions and n of the current location containing different parameters value of individual different powersxThe forward direction containing different parameters value of individual different powers
With the variable of reserve migration position；It is power with the nonlinear strength that each cos functions are born, takes IxIndividual cos functions weighted sum obtains reality
The function item of portion the 3rd, to remaining nx-IxIndividual cos functions weighted sum obtains the function item of real part the 4th；With each variable bear it is non-linear
Intensity is power, to nxThe variable weighting of individual deviation post is sued for peace to obtain the variable item of real part second；The function item of real part the 3rd is divided by real part
Plus the function item of real part the 4th and the 3rd real constant, gained multinomial is real part nonlinear function to second variable item；
Constructed real part nonlinear function is used for acting on effective plural component real part current location and bits of offset
The current state value put；
The structure of imaginary part nonlinear function is specially：
Imaginary part nonlinear functionBy k moment imaginary part vector YskNyIndividual componentN is built respectivelyyIt is individual
The cos functions and n of the current location containing different parameters value of different powersyThe forward direction containing different parameters value of individual different powers and
The variable of reserve migration position；It is power with the nonlinear strength that each cos functions are born, takes nyIndividual cos functions weighted sum obtains imaginary part
Second function；It is power with the nonlinear strength that each variable is born, to IyThe mutually multiplied imaginary part the 3rd of variable weighting of individual deviation post
Variable item, to remaining ny-IyIndividual variable weighting is sued for peace to obtain the variable item of imaginary part the 4th；Imaginary part ternary subtracts imaginary part the 4th
Variable item is divided by, plus the 4th real constant, gained multinomial is imaginary part nonlinear function after imaginary part second function；
Constructed imaginary part nonlinear function is used for acting on effective plural component imaginary part current location and bits of offset
The current state value put；
S4~S8, with step S4~S8 in sin type high-performance ranging code methods；
First pseudo noise code and the second pseudo noise code are carried out the He of mould two by S9, obtain the 4th good pseudorandom of pseudo-randomness
Code, i.e. cos types high-performance ranging code；If the pseudo-randomness of cos type high-performance ranging codes is destroyed, then need to adjust non-thread
Property function diffusion coefficient, function and/or variable nonlinear strength, finely tune function working frequency, phase, finely tune a reference value
Again the pseudo random number of binaryzation real part first and the pseudo random number of imaginary part first.
See a kind of multivariable space-time chaos satellite navigation high-performance ranging code implementation method that Fig. 1, the present invention are provided, use
In sin and cos type high-performance ranging codes, including：
S1~S2, with step S1~S2 in sin type high-performance ranging code methods；
S3 builds and is respectively acting on effective plural component real part and imaginary part current location and deviation post current state
The nonlinear function of value, specially：
Will act on the non-linear letter of the current state value of effective plural component real part current location and deviation post
Number is designated as real part nonlinear function respectively, will act on the current of effective plural component imaginary part current location and deviation post
The nonlinear function of state value is designated as imaginary part nonlinear function respectively；Wherein：
The structure of real part nonlinear function is specially：
Real part nonlinear functionBy k moment real part vectors ΧkNxIndividual componentN is built respectivelyx
The sin functions and n of the current location containing different parameters value of individual different powersxThe forward direction containing different parameters value of individual different powers
With the variable of reserve migration position；It is power with the nonlinear strength that each sin functions are born, takes IxIndividual sin functions weighted sum obtains reality
The function item of portion the 5th, to remaining nx-IxIndividual sin functions weighted sum obtains the function item of real part the 6th；With each variable bear it is non-linear
Intensity is power, to nxThe variable weighting of individual deviation post sues for peace to obtain real part ternary item；The function item of real part the 5th is divided by real part
Plus the function item of real part the 6th and the 5th real constant, gained multinomial is real part nonlinear function to ternary；
Constructed real part nonlinear function is used for acting on effective plural component real part current location and bits of offset
The current state value put；
The structure of imaginary part nonlinear function is specially：
Imaginary part nonlinear functionBy k moment imaginary part vector YskNyIndividual componentN is built respectivelyyIt is individual
The cos functions and n of the current location containing different parameters value of different powersyThe forward direction containing different parameters value of individual different powers and
The variable of reserve migration position；It is power with the nonlinear strength that each cos functions are born, takes nyIndividual cos functions weighted sum obtains imaginary part
3rd function item；It is power with the nonlinear strength that each variable is born, to IyThe mutually multiplied imaginary part the 5th of variable weighting of individual deviation post
Variable item, to remaining ny-IyIndividual variable weighting is sued for peace to obtain the variable item of imaginary part the 6th；The variable item of imaginary part the 5th subtracts imaginary part the 6th
Variable item is divided by, plus the 6th real constant, gained multinomial is imaginary part nonlinear function after the function item of imaginary part the 3rd；
Constructed imaginary part nonlinear function is used for acting on effective plural component imaginary part current location and bits of offset
The current state value put；
S4~S8, with step S4~S8 in sin type high-performance ranging code methods；
First pseudo noise code and the second pseudo noise code are carried out the He of mould two by S9, obtain the 5th good pseudorandom of pseudo-randomness
Code, i.e. sin and cos types high-performance ranging code；If the pseudo-randomness of sin and cos type high-performance ranging codes is destroyed, then need
The nonlinear strength of diffusion coefficient, function and/or the variable of nonlinear function is adjusted, working frequency, the phase of function is finely tuned,
Finely tune a reference value pseudo random number of binaryzation real part first and the pseudo random number of imaginary part first again.
In the step of above-mentioned three kinds of methods S2, extension plural component is used to aid in multivariable vector pseudo noise code to occur
The state value iteration of each plural component of position skew in device, its numerical value should be greater than or equal to default plural component
The maximum of position offset absolute value.
In the step of above-mentioned three kinds of methods S4, effective plural in the sequence of real numbers initialization G be made up of different real numbers
Component and G0The state value of middle extension plural component, described sequence of real numbers will ensure in G at effectively plural component
In chaos working condition, if sequence of real numbers is it cannot be guaranteed that chaos working condition is, it is necessary to adjust the diffusion coefficient of nonlinear function, letter
The nonlinear strength of number and/or variable.
In the step of above-mentioned three kinds of methods S4, using effectively plural component and G in pseudo-random number sequence initialization G0
The state value of middle extension plural component, specially：
Two linear pseudo-noise code generators are built respectively, are designated as the first linear pseudo-noise code generator and the second linear puppet
Random code generator；
The first linear pseudo-noise code generator and the second linear pseudo-noise code generator are driven respectively and from related register
The pseudo noise code of tap output the 6th and the 7th pseudo noise code；
Different decimals are respectively set to by the 6th pseudo noise code and the 7th pseudo noise code 0,1, are changed into the second pseudorandom
Number and the 3rd pseudo random number；If gained pseudo random number is it cannot be guaranteed that effectively plural component needs tune in chaos working condition
The nonlinear strength of the diffusion coefficient, function and/or variable of whole nonlinear function；
Second pseudo random number and the 3rd pseudo random number are the real part of respectively each plural component and the initial value of imaginary part.
In the step of above-mentioned three kinds of methods S8, a reference value is obtained using ranking method, i.e.,：The pseudo random number of real part first and imaginary part
The a reference value of the first pseudo random number is respectively wherein real number pseudo random number through size sequence gained median.
In the present invention, if gained sin types, cos types, sin and cos type high-performance ranging code pseudo-randomness be not good, it is necessary to
The nonlinear strength of Tuning function and/or variable, the diffusion coefficient of nonlinear function, or finely tune a reference value binaryzation reality again
The pseudo random number of portion first and the pseudo random number of imaginary part first, or finely tune working frequency, the phase of function.
See Figure 14, the present invention provides a kind of multivariable space-time chaos satellite navigation high-performance ranging code and realizes system, uses
In sin type high-performance ranging codes, including：
(1) multivariable vector pseudo-noise code generator builds module, for building multivariable vector pseudo-noise code generator
G, G are made up of single plural vector X+Yj, in If+N+|Ib| its component of dimensional linear space is Referred to as multivariable plural component,WithRespectively a series of order arrangements and respectively phase
The state component of mutual coupling；I represents the position number of multivariable plural component, i=1,2 ..., If+N+|Ib|, If, N be
Positive integer, IbIt is negative integer, IfTo spreading number, I before representingbRepresent backward spreading number；ixRepresent multivariable plural component reality
The sequence number of partial amount x, iyRepresent the sequence number of multivariable plural component imaginary y, ix=1,2 ..., nx,iy=1,
2,...,ny,nxAnd nyThe respectively quantity of state component x and y, preceding If×(nx+ny) and rear | Ib|×(nx+ny) individual multivariable answer
Number state component is designated as extending plural component, and its position number is respectively i ∈ [1, If] and i ∈ [If+N+1,If+N+|Ib
|]；Middle N × (nx+ny) individual multivariable plural component is designated as effective plural component, its position number i ∈ [If+1,
If+ N], effective plural componentWithRespectively constitute pseudo-noise code generator G1、G2；
(2) expansion module, G is constituted for extending plural component0, extension plural number of components is more than or equal to pre-
If position offset absolute value maximum；
(3) nonlinear function builds module, is respectively acting on effective plural component real part for structure and imaginary part is worked as
Front position and the nonlinear function of deviation post current state value, specially：Act on real part or imaginary part current location and skew
The nonlinear function of position current state value builds one group of not homogeneous respectively based on a group component of real part or imaginary part vector
The sin functions and/or variable containing different parameters value of side, are power with the negative nonlinear strength of function, by a part of letter of real part
Number or imaginary part function weighted sum obtain the first function of real part or imaginary part, and real part survival function weighted sum is obtained into the of real part
Two function items；It is power with the negative nonlinear strength of variable, by real part variable or a part of variable weighting summation of imaginary part or multiplication
The first variable item of real part or imaginary part is obtained, imaginary part surplus variable weighted sum is obtained the second variable item of imaginary part；By first function
, second function, the first variable item and the second variable item include adding by predetermined manner, subtract, at least two fortune in multiplication and division
The hybrid operation of calculation, along with correspondence real constant, gained multinomial is to act on real part or imaginary part current location and bits of offset
Put the nonlinear function of current state value；
(4) initialization module, the real number constituted for parameter initialization and using pseudo-random number sequence or by different real numbers
Effective plural component and G in sequence initialization G0The state value of middle extension plural component；
(5) state iteration module, for using multigroup different real part nonlinear function to effective plural component reality
The current state value of portion current location and deviation post is acted on, and obtains real part threshold；Use multigroup different imaginary part non-thread
Property the effective plural component imaginary part current location of function pair and the current state value of deviation post acted on, obtain imaginary part effect
Value；Real part threshold and imaginary part threshold are carried out adding respectively based on diffusion coefficient, are subtracted, multiplication and division or comprising adding, subtract, in multiplication and division
The hybrid operation of at least two computings, is produced with the plural pseudo-random number sequence of Annual distribution by state iteration；
(6) judge module, for the plural pseudo-random number sequence that is obtained using each current each effective plural component or
Sequence of real numbers, changes G0In each extension plural component state value, or carried out each other using these state values changed
Reconfigure arrangement；Then, next group of effective plural component in G is read, next group of effective plural component is gone to
State iteration module；When all effective equal completion status iteration of plural component in G, the extraction of real number pseudo-random number sequence is gone to
Module；
(7) real number pseudo-random number sequence extraction module, for respectively from G1And G2Middle correlated components tap extraction is with the time point
The real number pseudo-random number sequence of cloth, is designated as the pseudo random number of real part first and the pseudo random number of imaginary part first respectively；
(8) binarization block, for by each real number random number in the pseudo random number of real part first and the pseudo random number of imaginary part first
Compare with a reference value respectively in chronological order, if being more than a reference value, the real number random number value 1, otherwise value 0, you can obtain
Obtain first pseudo noise code and the second pseudo noise code of binaryzation；A reference value is respectively the description pseudo random number of real part first and imaginary part the
The statistic of real number pseudo random number median size in one pseudo random number；
(9) computing module of mould two, for the first pseudo noise code and the second pseudo noise code to be carried out into the He of mould two, obtains pseudo-randomness
The 3rd good pseudo noise code, i.e. sin types high-performance ranging code.
See Figure 14, the present invention provides a kind of multivariable space-time chaos satellite navigation high-performance ranging code and realizes system, uses
In cos type high-performance ranging codes, including：
(1) multivariable vector pseudo-noise code generator builds module, with multivariable vector pseudo noise code in said system four
Generator builds module；
(2) expansion module, with expansion module in said system four；
(3) nonlinear function builds module, is respectively acting on effective plural component real part for structure and imaginary part is worked as
Front position and the nonlinear function of deviation post current state value, specially：Act on real part or imaginary part current location and skew
The nonlinear function of position current state value builds one group of not homogeneous respectively based on a group component of real part or imaginary part vector
The cos functions and/or variable containing different parameters value of side, are power with the negative nonlinear strength of function, by a part of letter of real part
Number or imaginary part function weighted sum obtain the function item of real part the 3rd or imaginary part second function, and real part survival function weighted sum is obtained
The function item of real part the 4th；It is power with the negative nonlinear strength of variable, by real part variable or a part of variable weighting summation of imaginary part
Or the mutually multiplied variable item of real part second or imaginary part ternary, imaginary part surplus variable weighted sum is obtained into the variable of imaginary part the 4th
；By second function item, the 3rd function item, the 4th function item, the second variable item, ternary item and the 4th variable item by default
Mode include adding, subtracted, in multiplication and division at least two computings hybrid operation, along with correspondence real constant, gained multinomial
Act on the nonlinear function of real part or imaginary part current location and deviation post current state value；
(4) initialization module, with initialization module in said system four；
(5) state iteration module, with state iteration module in said system four；
(6) judge module, with judge module in said system four；
(7) real number pseudo-random number sequence extraction module, with real number pseudo-random number sequence extraction module in said system four；
(8) binarization block, with binarization block in said system four；
(9) computing module of mould two, for the first pseudo noise code and the second pseudo noise code to be carried out into the He of mould two, obtains pseudo-randomness
The 4th good pseudo noise code, i.e. cos types high-performance ranging code.
See Figure 14, the present invention provides a kind of multivariable space-time chaos satellite navigation high-performance ranging code and realizes system, uses
In sin and cos type high-performance ranging codes, including：
(1) multivariable vector pseudo-noise code generator builds module, with multivariable vector pseudo noise code in said system four
Generator builds module；
(2) expansion module, with expansion module in said system four；
(3) nonlinear function builds module, is respectively acting on effective plural component real part for structure and imaginary part is worked as
Front position and the nonlinear function of deviation post current state value, specially：Act on real part or imaginary part current location and skew
The nonlinear function of position current state value builds one group of not homogeneous respectively based on a group component of real part or imaginary part vector
Sin the or cos functions and/or variable containing different parameters value of side, with the negative nonlinear strength of function to weigh, by real part one
Function or imaginary part function weighted sum is divided to obtain the 3rd function item of the function item of real part the 5th or imaginary part, by the weighting of real part survival function
Sue for peace to obtain the function item of real part the 6th；It is power with the negative nonlinear strength of variable, real part variable or a part of variable of imaginary part is added
Power summation or mutually multiplied real part ternary or the variable item of imaginary part the 5th, imaginary part the 6th is obtained by imaginary part surplus variable weighted sum
Variable item；3rd function item, the 5th function item, the 6th function item, ternary, the 5th variable item and the 6th variable item are pressed
Predetermined manner include adding, subtracted, in multiplication and division at least two computings hybrid operation, along with correspondence real constant, institute is much
Item formula is the nonlinear function for acting on real part or imaginary part current location and deviation post current state value；
(4) initialization module, with initialization module in said system four；
(5) state iteration module, with state iteration module in said system four；
(6) judge module, with judge module in said system four；
(7) real number pseudo-random number sequence extraction module, with real number pseudo-random number sequence extraction module in said system four；
(8) binarization block, with binarization block in said system four；
(9) computing module of mould two, for the first pseudo noise code and the second pseudo noise code to be carried out into the He of mould two, obtains pseudo-randomness
The 5th good pseudo noise code, i.e. sin and cos types high-performance ranging code.
The present invention is described in further detail with reference to sin and cos type high-performance ranging code embodiments, according to Fig. 1
Realize that step is as follows：
S1 builds multivariable vector pseudo-noise code generator G, G and is made up of single plural vector X+Yj, in If+N+|Ib
| its component of dimensional linear space isReferred to as multivariable plural component,WithRespectively a series of order arrangements and the state component for intercoupling respectively；I represents multivariable plural
The position number of component, i=1,2 ..., If+N+|Ib|, If, N be positive integer, IbIt is negative integer, IfTo spreading number, I before representingb
Represent backward spreading number；ixRepresent the sequence number of multivariable plural component real component x, iyRepresent multivariable plural point
The sequence number of amount imaginary y, ix=1,2 ..., nx,iy=1,2 ..., ny,nxAnd nyThe respectively quantity of state component x and y,
Preceding If × (nx+ny) and rear | Ib|×(nx+ny) individual multivariable plural component be designated as extend plural component, its position-order
Number be respectively i ∈ [1, If] and i ∈ [If+N+1,If+N+|Ib|]；Middle N × (nx+ny) individual multivariable plural component is designated as
Effective plural component, its position number is i ∈ [If+1,If+ N], effective plural componentWithDifference structure
Into pseudo-noise code generator G1、G2；
S2 extension plural components constitute G0, extension plural number of components is more than or equal to default position offset
The maximum of absolute value；
To ensure follow-up computational efficiency, IfAnd IbValue is unsuitable excessive, typically takes the positive integer that its absolute value is not more than 10.
S3：Act on the nonlinear function and the structure of state iterative formula of effective plural component real part and imaginary part；
The structure of (3-1) nonlinear function：
In formula (1)：
K represents discrete time coordinate；
I represents the position number of effective plural component, i=If+1,If+2,...,If+N,IfPositive integer is with N；
ixRepresent the sequence number of effective plural component x, ix=1,2 ..., nx, nxRepresent effective plural component x's
Quantity；
iyRepresent the sequence number of effective plural component y, iy=1,2 ..., ny, nyRepresent effective plural component y's
Quantity；
Represent that k moment position number is effective plural component real part and imaginary part sequence number of i respectively
Respectively ixAnd iyComponent x and y state value；
The real part that acts on effective plural component and imaginary part current location and deviation post are represented respectively
Nonlinear function,Sin functions and cos functions are included respectively, and phase is respectivelyWith Point
The working frequency of the function not included is f0；
l1、l2Nonlinear function is represented respectivelySequence number, i.e.Expression acts on effective plural component reality
The l in portion1Individual nonlinear function,Expression acts on the l of effective plural component imaginary part2Individual nonlinear function；
IxIt is nonlinear functionMiddle part number of components, it is no more than nxPositive integer, IyIt is nonlinear functionMiddle part
Divide number of components, it is no more than nyPositive integer, its value sets itself as needed；
K moment nonlinear functions are represented respectivelyIn the serial number i that includesxSin functions and component
Nonlinear strength；K moment nonlinear functions are represented respectivelyIn the serial number i that includesyComponentWith
The nonlinear strength of cos functions；
Represent nonlinear functionIn include serial number ixSin functions range value；Represent non-linear letter
In numberComprising serial number iyCos functions range value；
Represent nonlinear functionMiddle serial number ixSin functions power；Represent nonlinear functionIn
Serial number ixComponentPower；Represent nonlinear functionMiddle serial number iyComponentPower；Represent
Nonlinear functionMiddle serial number iyCos functions power；
Nonlinear function is represented respectivelyWithIn include component forward and reverse position skew
Amount；
Represent k moment nonlinear functionsMiddle serial number ixSin functions state translational movement；Represent k
Moment nonlinear functionMiddle serial number ixComponentState translational movement；Represent k moment nonlinear functions
Middle serial number iyComponentState translational movement；Represent k moment nonlinear functionsMiddle serial number iyCos letters
Number state translational movement；
Nonlinear function is represented respectivelyWithComprising real constant, its value sets itself as needed.
(3-2) state iterative formula
If the code length that multivariable vector plural number pseudo-noise code generator G is produced is L, the non-of different parameters is contained using multigroup
Linear function acts on the current state value of effective plural component real part and imaginary part current location and deviation post, to spread
Coefficient be weights by comprising add, subtract, the hybrid operation mode of at least one of multiplication and division computing, make the generation of plural component
With the plural pseudo random number of Annual distribution.
In the real part and imaginary part of the effective plural component in (k+1) moment current location obtained by hybrid operation mode
Between state xk+1(i)、yk+1(i) and each state componentIterative formula it is as follows：
In formula (2)：
Real part and imaginary part intermediateness xk+1(i)、yk+1I () is respectively intended to preserve the non-linear letter in (k+1) moment current location
Number acts on the threshold of each state component；
It is the serial number i of (k+1) moment current location serial number ixReal part state component；
It is the serial number i of (k+1) moment current location serial number iyImaginary part state component；
Represent with componentThe nonlinear function of structureQuantity；Represent with componentThe non-linear letter for building
NumberQuantity；
The positive integer more than 0 is, its value sets itself as needed；To be not more thanPositive integer,To be not more thanPositive integer, its value sets itself as needed；
Represent the k moment with component respectivelyThe serial number l of structure1Nonlinear functionWith with componentBuild
Serial number l2Nonlinear functionDiffusion coefficient, be real number.
S4：The initialization of parameter and plural variable states value；
(4-1) parameter initialization
Parameter includes effective plural number of components N × (nx+ny) (N is series), code length L, nonlinear function quantity, non-
Linear function parameter, diffusion coefficient.Nonlinear function is the polynomial form for including different power functions and/or variable, ginseng
Number further includes the parameter and real constant of function and/or variable, function and/or variable, and the wherein parameter of function is work
Frequency, the power of function, the range value of function, the phase of function, position number, position offset, state translational movement, variable
Parameter is power, position number, position offset, the state translational movement of variable.
In the present invention, function power is the positive integer more than 1, no maximum requirement；The diffusion coefficient of each nonlinear function
It is real number.Nonlinear functionSin functions and cos functions are included respectively, and the working frequency of function is f0, phase difference
ForWithNonlinear functionRespectively act on plural component real part and imaginary part current location and
The nonlinear function of deviation post.
The present embodiment, takes series N=3, code length L=511, acts on the non-of effective plural component real part and imaginary part
The quantity of linear functionBuild the state point of the nonlinear function of real part and imaginary part
Amount x and y quantity is respectively nx=3, ny=3, the range value of all sin functions and cos functionsWorking frequency f0=1.023Hz, real part and each component of imaginary part
Forward direction extension plural number of components is 5, and backward extension plural number of components is 4.
The nonlinear function of effective plural component real part current location and deviation post is acted on respectively by containing difference
The sin functions and component of parameterComposition, the quantity of sin functions is Ix=nx=3, componentQuantity be nx=3, wherein the
One nonlinear functionComprising sin function parameters be：Secondary number formulary is respectivelyWithNonlinear strength is respectively State is put down
Shifting amount is respectively Comprising componentParameter
For：Secondary number formulary is respectivelyWithNonlinear strength is respectively WithPosition offset is respectivelyWithShape
State translational movement is respectivelyWith Comprising it is real often
NumberSecond nonlinear functionComprising sin function parameters be：Secondary number formulary is respectivelyWithNonlinear strength is respectively State translational movement is respectively Comprising componentParameter be：Secondary number formulary is respectivelyWithIt is non-linear
Intensity is respectivelyWithPosition offset is respectivelyWithState translational movement is respectively
WithComprising real constant3rd nonlinear functionComprising sin letters
Counting parameter is：Secondary number formulary is respectivelyWithNonlinear strength is respectively State translational movement is respectively Comprising componentParameter be：Secondary number formulary is respectivelyWithNonlinear strength is respectively
WithPosition offset is respectively WithState is put down
Shifting amount is respectively With Comprising real constant
The nonlinear function of effective plural component imaginary part current location and deviation post is acted on respectively by containing difference
The cos functions and component of parameterComposition, the quantity of cos functions is ny=3, componentQuantity be Iy=2 and ny=3, wherein
First nonlinear functionComprising cos function parameters be：Secondary number formulary is respectivelyWithNonlinear strength is respectively State
Translational movement is respectively Comprising component's
Parameter is：Secondary number formulary is respectivelyWithNonlinear strength is respectively WithPosition offset is respectively WithState translational movement is respectively With Comprising real constantSecond nonlinear functionComprising cos functions
Parameter is：Secondary number formulary is respectivelyWithNonlinear strength is respectively State translational movement is respectively Comprising componentParameter be：Secondary number formulary is respectively WithNonlinear strength is respectively WithPosition offset is respectively WithState is translated
Amount is respectively With Comprising real constant3rd nonlinear functionComprising cos function parameters be：Secondary number formulary is respectively WithNonlinear strength is respectively State translational movement is respectively Comprising componentParameter be：Secondary number formulary is respectively WithIt is non-linear
Intensity is respectively WithPosition offset is respectively WithState translational movement is respectively With Comprising real constant
Act on the phase of the sin functions that the nonlinear function of effective plural component real part current location is includedFor position number is multiplied byAct on what the nonlinear function of effective plural component imaginary part current location was included
The phase of cos functionsFor position number is multiplied by
Act on the nonlinear function of effective plural component real part current location and deviation postCorresponding diffusion system
Number is respectivelyWith
Act on the nonlinear function of effective plural component imaginary part current location and deviation postCorresponding diffusion
Coefficient is respectivelyWith
The initialization of (4-2) plural component state value
Can be initialized using pseudo-random number sequence or by the sequence of real numbers that different real numbers are constituted, and to be ensured effectively multiple
Number state component serves are in chaos state.The state value of plural component, pseudorandom are initialized according to pseudo-random number sequence
Number Sequence can be obtained using linear shift register by related register tap.
The specific implementation that plural component is initialized using the sequence of real numbers being made up of different real numbers is provided below
Journey.
Each plural component real part position number is multiplied by 0.1, then with moment and 0.00001 product addition, structure
Into value of real part；Each plural component imaginary part position number is multiplied by 0.2 again, then with moment and 0.00001 product addition,
Constitute imaginary values.Value of real part and imaginary values constitute each plural component initialization value.
S5：State is carried out respectively to each effective plural component in plural pseudo-noise code generator G according to formula (2) to change
In generation, produce the plural pseudo-random number sequence with 530 of Annual distribution；
State iteration of the present invention is based on nonlinear function weights and acts on realization, and formula (2) is shown in the effect of nonlinear function weights
Shown hybrid operation mode.
S6：To G0In each extension plural component, 530 obtained with S5 are multiplied by after 0.001 with its position number multiple
The value of real part of number pseudo-random number sequence is added, and obtains the real part state value of the extension plural component；Multiplied with its position number again
It is added with the imaginary values after 0.003 with 530 plural pseudo-random number sequences that S5 is obtained, obtains the void of the extension plural component
Portion's state value, realizes G0In each extension plural component state value modification；Then next group of effective plural in G is read
Component, performs step S5.All effective equal completion status iteration of plural component in G, perform step S7.
During effective plural component carries out state iterative process in G, the iterative value of each plural component is distinguished
It is stored in the intermediateness component of real part and imaginary part, state value tap extraction is carried out in order to subsequent step.
S7：Postpone 5 seconds, to avoid initial non-chaos oscillator signal, from G12nd, the 3rd intermediateness component tap carry
One group of pseudo random number of real part first of time state distribution of length 511 second is taken, Fig. 2~4 is seen, from G21st, the 3rd middle shape
511 seconds pseudo random numbers of imaginary part first of time state distribution of state component tap extraction length.
S8 by each real number random number in the pseudo random number of real part first and the pseudo random number of imaginary part first in chronological order respectively with
A reference value compares, if being more than a reference value, the real number random number value 1, otherwise value 0, you can obtain binaryzation first is pseudo-
Random code and the second pseudo noise code, are shown in Fig. 5~7；A reference value be respectively description the pseudo random number of real part first and the puppet of imaginary part first with
The statistic of real number pseudo random number median size in machine number；
First pseudo noise code and the second pseudo noise code are carried out the He of mould two by S9, obtain the 5th good pseudorandom of pseudo-randomness
Code, i.e. sin and cos types high-performance ranging code, are shown in Fig. 8~10.
The pseudo-randomness evaluation of ranging code is shown in Table 1~6 and Figure 11~13.Pseudo-randomness can using balance, the distance of swimming and
Autocorrelation is evaluated.Balance is the percentage that numerical value 0 and 1 accounts for sum in ranging code, ideally 0 and 1 accounts for sum respectively
50%.The distance of swimming is that the distance of swimming of different length in ranging code accounts for the percentage of total number of runs, ideally each length runs
The percentage for accounting for total number of runs isWherein, a represents run length.Autocorrelation is ranging code auto-correlation function
δ function characteristics.
The balance of ranging code shown in the Fig. 8 of table 1
| Numerical value | Percent of total (%) |
| 0 | 49.90 |
| 1 | 50.10 |
The distance of swimming of ranging code shown in the Fig. 8 of table 2
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 47.54 |
| 2 | 29.51 |
| 3 | 13.11 |
| 4 | 5.33 |
| 5 | 2.05 |
| 6 | 1.23 |
The balance of ranging code shown in the Fig. 9 of table 3
| Numerical value | Percent of total (%) |
| 0 | 50.68 |
| 1 | 49.32 |
The distance of swimming of ranging code shown in the Fig. 9 of table 4
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 49.17 |
| 2 | 20.25 |
| 3 | 16.94 |
| 4 | 7.02 |
| 5 | 3.31 |
| 6 | 1.65 |
The balance of ranging code shown in the Figure 10 of table 5
| Numerical value | Percent of total (%) |
| 0 | 55.19 |
| 1 | 44.81 |
The Figure 10 of table 6 shows the distance of swimming of ranging code
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 46.25 |
| 2 | 25.42 |
| 3 | 13.75 |
| 4 | 7.92 |
| 5 | 2.92 |
| 6 | 0.83 |
It should be appreciated that the part that this specification is not elaborated belongs to prior art.
It should be appreciated that the above-mentioned description for preferred embodiment is more detailed, therefore can not be considered to this
The limitation of invention patent protection scope, one of ordinary skill in the art is not departing from power of the present invention under enlightenment of the invention
Profit requires under protected ambit, can also make replacement or deform, each falls within protection scope of the present invention, this hair
It is bright scope is claimed to be determined by the appended claims.
Claims (10)
1. a kind of multivariable space-time chaos satellite navigation high-performance ranging code implementation method, it is characterised in that comprise the following steps：
Step 1：Build multivariable vector pseudo-noise code generator G；
Step 2：Judge whether to produce sin type ranging codes；
If so, then performing following step 5；
If it is not, then performing following step 3；
Step 3：Judge whether to produce cos type ranging codes；
If so, then performing following step 5；
If it is not, then performing following step 4；
Step 4：Produce sin types and cos type ranging codes；
Step 5：Extension plural component G0；
Step 6：Build multivariable nonlinearity function；
Step 7：In parameter initialization and the sequence of real numbers initialization G constituted using pseudo-random number sequence or by different real numbers effectively
Plural component and G0The state value of middle extension plural component；
Step 8：Effective plural component in G is taken, is produced with the plural pseudo-random number sequence of Annual distribution by state iteration；
Step 9：Whether effectively plural component takes in judging G；
If so, then performing following step 10；
If it is not, then changing extension component state value, and turn round execution step 8；
Step 10：Real part pseudo random number and imaginary part pseudo random number in G are taken respectively；
Step 11：By real part pseudo random number and imaginary part pseudo random number binaryzation and the He of mould two is carried out, obtain high-performance ranging code.
2. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In：Build multivariable vector pseudo-noise code generator G, G described in step 1 to be made up of single plural vector X+Yj, in If+
N+|Ib| its component of dimensional linear space isReferred to as multivariable plural divides
Amount,WithRespectively a series of order arrangements and the state component for intercoupling respectively；I represents multivariable plural
The position number of component, i=1,2 ..., If+N+|Ib|, If, N be positive integer, IbIt is negative integer, IfTo spreading number, I before representingb
Represent backward spreading number；ixRepresent the sequence number of multivariable plural component real component x, iyRepresent multivariable plural point
The sequence number of amount imaginary y, ix=1,2 ..., nx, iy=1,2 ..., ny, nxAnd nyRespectively the quantity of state component x and y, preceding
If×(nx+ny) and rear | Ib|×(nx+ny) individual multivariable plural component be designated as extend plural component, its position number
Difference i ∈ [1, If] and i ∈ [If+ N+1, If+N+|Ib|]；Middle N × (nx+ny) individual multivariable plural component is designated as effectively
Plural component, its position number i ∈ [If+ 1, If+ N], effective plural componentWithRespectively constitute it is pseudo- with
Machine code generator G1、G2。
3. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In：Plural component is extended described in step 5 and constitutes G0, it is inclined that extension plural number of components is more than or equal to default position
The maximum of shifting amount absolute value.
4. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In：
Multivariable nonlinearity function is built described in step 6, for for sin type ranging codes, it implements process is：Build
The nonlinear function of effective plural component real part or imaginary part current location and deviation post current state value is respectively acting on,
The nonlinear function of real part or imaginary part current location and deviation post current state value is specially acted on to be sweared with real part or imaginary part
The current location containing different parameters value and partially of one group of identical power or different power is built based on one group component of amount respectively
The sin functions and/or sin compound functions and/or variable of pan position, with each sin functions and/or the negative non-thread of sin compound functions
Property intensity for power, take a part of sin functions of real part and/or sin compound functions or imaginary part sin functions and/or sin compound functions add
Power computing obtains real part or imaginary part first function, sin functions and/or sin compound function ranking operations to real part remainder
Obtain real part second function；It is power with the nonlinear strength that each variable is born, takes real part variable or a part of variable weighting fortune of imaginary part
Real part or the variable item of imaginary part first are calculated to obtain, the variable weighting computing to imaginary part remainder obtains the variable item of imaginary part second；By first
Function item, second function, the first variable item and the second variable item include adding, subtracted by predetermined manner, at least two in multiplication and division
The hybrid operation of kind of computing, along with correspondence real constant, gained multinomial is to act on real part or imaginary part current location and partially
The nonlinear function of pan position current state value；
Multivariable nonlinearity function is built described in step 6, for for cos type ranging codes, it implements process is：Build
The nonlinear function of effective plural component real part or imaginary part current location and deviation post current state value is respectively acting on,
The nonlinear function of real part or imaginary part current location and deviation post current state value is specially acted on to be sweared with real part or imaginary part
The current location containing different parameters value and partially of one group of identical power or different power is built based on one group component of amount respectively
The cos functions and/or cos compound functions and/or variable of pan position, with each cos functions and/or the negative non-thread of cos compound functions
Property intensity for power, take a part of cos functions of real part and/or cos compound functions or imaginary part cos functions and/or cos compound functions add
Power computing obtains the function item of real part the 3rd or imaginary part second function, and the cos functions and/or cos to real part remainder are combined letter
Number ranking operation obtains the function item of real part the 4th；It is power with the nonlinear strength that each variable is born, takes real part variable or an imaginary part part
Variable weighting computing obtains the variable item of real part second or imaginary part ternary, and imaginary part remainder variable weighting computing is obtained into imaginary part
4th variable item；By second function item, the 3rd function item, the 4th function item, the second variable item, ternary item and the 4th variable
Include add by predetermined manner, subtract, in multiplication and division at least two computings hybrid operation, along with correspondingly real constant item, institute
It is the nonlinear function for acting on real part or imaginary part current location and deviation post current state value to obtain multinomial；
Multivariable nonlinearity function is built described in step 6, for for sin types ranging code and cos type ranging codes, its specific reality
Now process is：Structure is respectively acting on effective plural component real part or imaginary part current location and deviation post current state value
Nonlinear function, specially act on the nonlinear function of real part or imaginary part current location and deviation post current state value with
Built respectively based on one group component of real part or imaginary part vector one group of identical power or different power containing different parameters value
Current location and sin the or cos functions and/or sin or cos compound functions and/or variable of deviation post, with each sin or cos letters
Number and/or the negative nonlinear strength of sin or cos compound functions are power, take a part of sin or cos functions of real part and/or sin or
Cos compound functions or imaginary part sin or cos function and/or sin or cos compound function ranking operations obtain the function item of real part the 5th or
The function item of imaginary part the 3rd, sin or cos functions and/or sin or cos compound function ranking operations to real part remainder obtain reality
The function item of portion the 6th；It is power with the nonlinear strength that each variable is born, takes real part variable or a part of variable weighting computing of imaginary part is obtained
Real part ternary or the variable item of imaginary part the 5th, the variable weighting computing to imaginary part remainder obtain the variable item of imaginary part the 6th；
By the 3rd function item, the 5th function item, the 6th function item, ternary, the 5th variable item and the 6th variable item by default side
Formula include adding, subtracted, in multiplication and division at least two computings hybrid operation, along with correspondence real constant, gained multinomial is
Act on the nonlinear function of real part or imaginary part current location and deviation post current state value；
In step 6, multivariable nonlinearity function is divided based on a group component of multivariable plural component real part or imaginary part
One group of identical power or the different current locations containing different parameters value of power and the function of deviation post are not built and/or are answered
Close function and/or variable, as each function and/or the nonlinear strength and real constant structure of compound function and/or variable weights
Into multinomial, the parameter that wherein function and/or compound function are included be working frequency, the power of function, the range value of function,
The phase value of function, position number, position offset, state translational movement, the parameter that variable is included are power, the position-order of variable
Number, position offset, state translational movement, in parameter real number precision be 10-5, it is with nonlinear strength in multivariable nonlinearity function
1. and 2. computing carried out between the function and/or compound function and/or variable of weights, wherein being 1. the fortune for adding, subtract, multiplying or removing
Calculate, be 2. default hybrid operation.
5. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In：
It is each effectively multiple in the sequence of real numbers initialization G for constituted or rearranged by pseudo-random number sequence or different real numbers in step 7
Number state component and G0In each extension plural component state value, described pseudo-random number sequence or sequence of real numbers will ensure
Each effective plural component is in chaos working condition in G, if pseudo-random number sequence or sequence of real numbers are it cannot be guaranteed that chaos work
Make state, it is necessary to nonlinear strength, the diffusion coefficient of nonlinear function of Tuning function and/or compound function and/or variable.
6. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In the process that implements of step 8 is：Effective plural component real part is worked as respectively using real part multivariable nonlinearity function
The current state value of front position and deviation post is acted on, and obtains real part threshold；Use imaginary part multivariable nonlinearity function pair
The current state value of effective plural component imaginary part current location and deviation post is acted on, and obtains imaginary part threshold；To reality
Portion's threshold and imaginary part threshold carry out 1., 2., 3. or 4. computing, by state iteration produce the plural number with Annual distribution pseudo- with
Machine Number Sequence；1. it is wherein the computing for adding, subtract, multiplying or removing, is 3. with diffusion coefficient as weights 2. to preset hybrid operation
Plus, the computing that subtracts, multiply or remove, be 4. the default hybrid operation with diffusion coefficient as weights.
7. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 1, its feature exists
In：Change extension component state value described in step 9, be the plural pseudorandom obtained using current each effective plural component
Number Sequence or sequence of real numbers, are respectively modified G0In each extension plural component state value, or using the state of these modifications
Value carries out reconfiguring arrangement each other.
8. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 2, its feature exists
In the process that implements of step 10 is：Respectively from G1And G2Middle correlated components tap extraction with Annual distribution real number pseudorandom
Number Sequence, is designated as the pseudo random number of real part first and the pseudo random number of imaginary part first respectively.
9. multivariable space-time chaos satellite navigation high-performance ranging code implementation method according to claim 2, its feature exists
In by real part pseudo random number and imaginary part pseudo random number binaryzation described in step 11, it implements process is：By real part first
Each real number random number compares with a reference value respectively in chronological order in pseudo random number and the pseudo random number of imaginary part first, if being more than benchmark
It is worth, then the real number random number value 1, otherwise value 0, you can obtain first pseudo noise code and the second pseudo noise code of binaryzation；
A reference value is respectively real number pseudo random number median size in the description pseudo random number of real part first and the pseudo random number of imaginary part first
Statistic.
10. a kind of multivariable space-time chaos satellite navigation high-performance ranging code realizes system, it is characterised in that：Including multivariable arrow
Amount pseudo-noise code generator build module, the first judge module, expansion module, nonlinear function build module, initialization module,
State iteration module, the second judge module, real part pseudo random number and imaginary part pseudo random number extraction module, binarization block, mould two
Computing module；
The multivariable vector pseudo-noise code generator builds module, for building multivariable vector pseudo-noise code generator G；
First judge module, for performing following judgements and operation；
Judge whether to produce sin type ranging codes, if then performing method and step 5, if otherwise performing method and step 3；Judge whether
Cos type ranging codes are produced, if then performing method and step 5；If otherwise performing method and step 4；Produce sin types and the range finding of cos types
Code；
The expansion module, for extending plural component G0；
The nonlinear function builds module, for building multivariable nonlinearity function；
The initialization module, the sequence of real numbers constituted for parameter initialization and using pseudo-random number sequence or by different real numbers
Effective plural component and G in initialization G0The state value of middle extension plural component；
The state iteration module, for taking in G effectively plural component, answering with Annual distribution is produced by state iteration
Number pseudo-random number sequence；
Second judge module, for judging whether effectively plural component takes in G；If so, then performing method and step
10；If it is not, then changing extension component state value, and turn round execution method and step 8；
The real part pseudo random number and imaginary part pseudo random number extraction module, for taking real part pseudo random number and imaginary part pseudorandom respectively
Number；
The binarization block, for by real part pseudo random number and imaginary part pseudo random number binaryzation；
The computing module of the mould two, for the real part pseudo random number and imaginary part pseudo random number after binaryzation to be carried out into the He of mould two, obtains
Obtain high-performance ranging code.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611170030.2A CN106772476B (en) | 2016-12-16 | 2016-12-16 | Multivariable space-time chaos satellite navigation high-performance ranging code method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611170030.2A CN106772476B (en) | 2016-12-16 | 2016-12-16 | Multivariable space-time chaos satellite navigation high-performance ranging code method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106772476A true CN106772476A (en) | 2017-05-31 |
| CN106772476B CN106772476B (en) | 2019-01-04 |
Family
ID=58893331
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611170030.2A Active CN106772476B (en) | 2016-12-16 | 2016-12-16 | Multivariable space-time chaos satellite navigation high-performance ranging code method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772476B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117406249A (en) * | 2023-12-14 | 2024-01-16 | 成都安则优科技有限公司 | Anti-unmanned aerial vehicle satellite navigation equipment and data playback method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101441265A (en) * | 2008-12-17 | 2009-05-27 | 北京航空航天大学 | Method for capturing navigation satellite signal by using chaos system |
| CN102156288A (en) * | 2008-12-17 | 2011-08-17 | 北京航空航天大学 | Method for acquiring navigation satellite signal by using chaotic system |
| CN103988094A (en) * | 2011-08-16 | 2014-08-13 | 欧洲太空署 | A navigation system using spreading codes based on pseudo-random noise sequences |
| WO2016027727A1 (en) * | 2014-08-20 | 2016-02-25 | ソニー株式会社 | Reception device and reception method for global navigation satellite system, and program |
-
2016
- 2016-12-16 CN CN201611170030.2A patent/CN106772476B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101441265A (en) * | 2008-12-17 | 2009-05-27 | 北京航空航天大学 | Method for capturing navigation satellite signal by using chaos system |
| CN102156288A (en) * | 2008-12-17 | 2011-08-17 | 北京航空航天大学 | Method for acquiring navigation satellite signal by using chaotic system |
| CN103988094A (en) * | 2011-08-16 | 2014-08-13 | 欧洲太空署 | A navigation system using spreading codes based on pseudo-random noise sequences |
| WO2016027727A1 (en) * | 2014-08-20 | 2016-02-25 | ソニー株式会社 | Reception device and reception method for global navigation satellite system, and program |
Non-Patent Citations (1)
| Title |
|---|
| 洪婷 等: "《时空混沌伪随机序列及其在C/A码中的应用》", 《计算机仿真》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117406249A (en) * | 2023-12-14 | 2024-01-16 | 成都安则优科技有限公司 | Anti-unmanned aerial vehicle satellite navigation equipment and data playback method |
| CN117406249B (en) * | 2023-12-14 | 2024-03-08 | 成都安则优科技有限公司 | Anti-unmanned aerial vehicle satellite navigation equipment and data playback method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106772476B (en) | 2019-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Adams | On the nonexistence of elements of Hopf invariant one | |
| CN106291616A (en) | Space-time chaos vector pseudo-noise code generator offset carrier modulator approach and system | |
| CN100562762C (en) | The method and system that is used for positioning global positioning system correlated peak signal | |
| Krips et al. | Hybrid model of fixed and floating point numbers in secure multiparty computations | |
| CN101813768A (en) | Method for designing quadrature phase-coded signal based on hyperchaos mapping | |
| CN106772476B (en) | Multivariable space-time chaos satellite navigation high-performance ranging code method and system | |
| US4706211A (en) | Digital multiplying circuit | |
| Chervyakov et al. | The effective neural network implementation of the secret sharing scheme with the use of matrix projections on FPGA | |
| CN106301753B (en) | Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system | |
| Chu et al. | Design of FSM-based function with reduced number of states in integral stochastic computing | |
| CN106597493B (en) | Nonlinear function drives space-time chaos satellite navigation ranging code realization method and system | |
| Holman et al. | Fermion masses in E8× E′ 8 superstring theories | |
| Lapkova | Class number one problem for real quadratic fields of a certain type | |
| CN101127575A (en) | An equably distributed random number generator and its generation method | |
| CN111190571B (en) | Modular multiplication circuit based on binary domain expansion and control method thereof | |
| Isbell | On the multiplicative semigroup of commutative ring | |
| Giuliani et al. | Critical “dimension” in shell model turbulence | |
| Jagadeeshkumar et al. | A novel design of low power and high speed hybrid multiplier | |
| Cook Jr et al. | Evaluation of the H-function integral and of probability distributions of functions of independent random variables | |
| CN115543911B (en) | Method for calculating computing capacity of heterogeneous computing equipment | |
| Barkalov et al. | Design of fpga-based mealy fsms with counters | |
| Choi et al. | Exact formulas for traces of singular moduli of higher level modular functions | |
| CN102710355B (en) | Method for developing and designing robust self-synchronizing dissipative coupling complex network | |
| Liu | A high quality PN sequence generator based on chaotic maps | |
| Ismail et al. | Dynamics of fractional and double-humped logistic maps versus the conventional one |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |