CN106301753B - Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system - Google Patents
Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system Download PDFInfo
- Publication number
- CN106301753B CN106301753B CN201610608799.1A CN201610608799A CN106301753B CN 106301753 B CN106301753 B CN 106301753B CN 201610608799 A CN201610608799 A CN 201610608799A CN 106301753 B CN106301753 B CN 106301753B
- Authority
- CN
- China
- Prior art keywords
- variable
- pseudo
- plural
- random number
- function
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000009792 diffusion process Methods 0.000 claims abstract description 37
- 238000000605 extraction Methods 0.000 claims abstract description 8
- 230000002596 correlated effect Effects 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims abstract description 7
- 230000000875 corresponding effect Effects 0.000 claims description 32
- 239000000284 extract Substances 0.000 claims description 6
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 239000013589 supplement Substances 0.000 claims description 4
- 230000000739 chaotic effect Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 2
- 150000001875 compounds Chemical group 0.000 claims 14
- 238000004891 communication Methods 0.000 abstract description 3
- 239000002131 composite material Substances 0.000 abstract 1
- 238000005259 measurement Methods 0.000 abstract 1
- 230000008859 change Effects 0.000 description 10
- 230000009182 swimming Effects 0.000 description 9
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000005311 autocorrelation function Methods 0.000 description 4
- 229910002056 binary alloy Inorganic materials 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 201000003373 familial cold autoinflammatory syndrome 3 Diseases 0.000 description 3
- 238000012886 linear function Methods 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
- H04L9/001—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols using chaotic signals
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
- H04L9/06—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols the encryption apparatus using shift registers or memories for block-wise or stream coding, e.g. DES systems or RC4; Hash functions; Pseudorandom sequence generators
- H04L9/065—Encryption by serially and continuously modifying data stream elements, e.g. stream cipher systems, RC4, SEAL or A5/3
- H04L9/0656—Pseudorandom key sequence combined element-for-element with data sequence, e.g. one-time-pad [OTP] or Vernam's cipher
- H04L9/0662—Pseudorandom key sequence combined element-for-element with data sequence, e.g. one-time-pad [OTP] or Vernam's cipher with particular pseudorandom sequence generator
Landscapes
- Engineering & Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a kind of multi-state variable space-time chaos plural number pseudo-noise code generator realization method and systems, the present invention constructs a series of plural variable for being arranged and being intercoupled by sequences and constitutes plural pseudo-noise code generator, first, different nonlinear function effects are respectively adopted to the current state value of plural variable real and imaginary parts current location and deviation post, pass through adding using diffusion coefficient as weight, subtract, multiply, or the hybrid operation removed, the plural pseudo-random number sequence being distributed at any time is generated by state iteration, then, binaryzation and mould two and required ranging code is obtained respectively from correlated variables real and imaginary parts tap extraction real number pseudo-random number sequence.It the composite can be widely applied to satellite navigation system, it can also be used to various range-measurement systems, communication system, Broadcast and TV system, control system etc..
Description
Technical field
The present invention relates to technical field of satellite navigation, specially a kind of multi-state variable space-time chaos plural number pseudo noise code hair
Raw device realization method and system.
Background technique
Current four large satellite navigation system in the world have GPS (the Global Positioning System) satellite in the U.S.
Navigation system, GLONASS (Global Navigation Satellite System) satellite navigation system, European Union of Russia
Galileo satellite navigation system and China Beidou (BeiDou) satellite navigation system.Except GLONASS uses frequency division multiple access
Outside (FDMA, Frequency Division Multiple Access) communication mode, remaining satellite navigation system is all made of code
Divide multiple access (CDMA, Code Division Multiple Access) communication mode.Ranging code used by them is divided into civilian
Two class of bigness scale code and military accurate code can be carried out target coarse positioning using the satellite navigation system of bigness scale code, and use accurate code
Satellite navigation system can carry out target with high precision positioning.
The pseudo-noise code generator that satellite navigation system uses is divided into two kinds, and one is current satellite navigation systems to make
Binary system pseudo-noise code generator, ranging code use one or several certain series linear feedbacks by linear method
Shift register first passes through one section of short binary sequence to initialization of register, is then generated by register shift.Bigness scale code is realized
C/A (the Coarse Acquisition Code) code of aspect GPS L1 signal is linearly moved totally by two parallel 10 grades (20 grades)
Bit register generates, and code length is 1023；The bigness scale code of Galileo E1 signal is produced by two parallel linear shift registers
Raw to truncate combined M code realization, code length is 4092；The linear displacement that the ranging code of GLONASS is 9 grades by a maximum length
Register generates (M sequence), and code length is 511；Ranging code C at BeiDou B1I, B2I signalB1ICode and CB2ICode is by two
11 grades of parallel (totally 22 grades) linear shift registers generate, and code length is 2046.Accurate code realizes that aspect only has GPS to provide
Implementation method is generated using the linear shift register of two parallel 12 grades (24 grades) totally respectively.Due to certain length
Pseudo noise code needs the shift register of certain series to generate by displacement, leads to the generally existing complexity of ranging code generated
It is low, safety is poor, code length is fixed and it is shorter, number of codes is limited the disadvantages of, and shift register also need to undertake linear feedback and
Satellite Phase shares out the work, and makes that its structure is complicated.
Another kind is the real number pseudo-noise code generator discussed on document, it is to be mixed by nonlinear method using space-time
Ignorant system can provide that large number of, irrelevant, class is random and determining reproducible signal to the sensitive dependence of initial value,
These signals have the characteristics that pseudo-randomness, aperiodicity, long-term unpredictability and ergodic, and space-time chaos are unidirectional
Coupled map lattice model is applied to GPS system, designs 20 grades of space-time chaos real number pseudo noise codes with one group of 20 space lattice
Generator replaces two 10 grades of parallel linear shift registers of GPS with it.It is first initial with same real number plaid matching point variable state
Change, the real number of state distribution at any time is then generated 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 export real number pseudo random number,
Binaryzation and mould two and be ranging code.Nonlinear function is by variable, parameter (power, position number, the work including variable of variable
For the nonlinear strength of variable weight) and integer constant item composition, variable, the parameter of variable and integer constant item are known as nonlinear function
Parameter.Nonlinear function generally uses univariate polynomial form, and (i.e. low order side's variable is multiplied by negative non-linear strong
Degree adds integer constant) f (x)=- δ x2+ 1, δ are the nonlinear strength of variable, and diffusion coefficient is (as nonlinear function threshold
Weight) be real number field [0,1] in real number, act on the diffusion system of the nonlinear function threshold of different lattice point state variables
Several summations is 1, and different nonlinear function thresholds can only carry out add operation between each other.Real number pseudo-noise code generator
The ranging code code length that current binary system pseudo-noise code generator can be overcome to generate is fixed and the shorter and limited disadvantage of number of codes,
But due to using a kind of nonlinear function and generally quadratic power function effect, and data precision is 10-2, make the ranging generated
Code complexity is not high, safety is not strong, and use space lattice point number is excessive, and each lattice point state variable once only generates a kind of reality
Number pseudo random number.Furthermore to real number pseudo random number binaryzation, currently used processing method is mainly Invariant Distribution method, adjacent value ratio
Compared with method and adjusting thresholds and sequence truncation method: Invariant Distribution method needs the real number pseudo random number generated to each lattice point state variable
Probability statistical analysis is carried out, to determine a reference value again to pseudo-random number sequence binaryzation；Neighbour's value comparison method passes through plaid matching dotted state
The real number pseudo random number consecutive value that variable generates, which compares, carries out pseudo-random number sequence binaryzation；Adjusting thresholds and sequence truncation method are first
It is new sequence that the real number pseudo random number for first needing the dotted state variable of plaid matching to generate intercepts by a certain percentage, then flat in new sequence again
Mean value nearby searches for optimal threshold, realizes former pseudo-random number sequence binaryzation according to threshold value.
Summary of the invention
In view of the above technical problems, high performance satellite navigational range is generated by nonlinear method the present invention provides a kind of
The multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system of code.
In order to solve the above technical problems, the present invention adopts the following technical scheme that:
One, a kind of multi-state variable space-time chaos plural number pseudo-noise code generator implementation method, comprising:
S1 building plural number pseudo-noise code generator G, G be made of a series of plural variables { x (i)+y (i) j } (i=1,
2 ... N, N are positive integer), { x (i) } and { y (i) } respectively constitutes pseudo-noise code generator G1、G2, { x (i) } and { y (i) } is one
The state variable that sequential series are arranged and intercoupled respectively；I indicates plural variable position serial number, and include in G is to have
Imitate plural variable；
S2 optimizes plural number pseudo-noise code generator G and obtains minimum series, using minimum series as the series N of G；
S3 constitutes G in the forward or backward supplement extension plural variable of effective plural variable0, extend plural
Variable number is greater than or equal to preset position offset maximum value；
S4 constructs nonlinear function, will act on the state of effective plural variable real part current location and deviation post
The nonlinear function of value is denoted as real part current location nonlinear function and real part deviation post nonlinear function respectively, will act on
The nonlinear function of the state value of effective plural variable imaginary part current location and deviation post is denoted as imaginary part present bit respectively
Set nonlinear function and imaginary part deviation post nonlinear function；Wherein:
The building of real part current location nonlinear function specifically:
Real part current location nonlinear function is by LL12The function and LL of the value containing different parameters of a difference power13A difference
The variable of the value containing different parameters of power is constituted；It is power with the nonlinear strength that each function is born, takes LL11A function weighted sum obtains
Real part first function item, to remaining (LL12-LL11) a function weighted sum obtains real part second function item；It is born with each variable
Nonlinear strength is power, to LL13A variable weighting is summed to obtain the first variable item of real part；The first variable item of real part is divided by real part second
Function item, in addition real part first function item and the first real constant item, gained multinomial, that is, real part current location nonlinear function；
The building of real part deviation post nonlinear function specifically:
Real part deviation post nonlinear function is by LL21The variable of the value containing different parameters of a difference power is constituted；With each change
The negative nonlinear strength of amount is power, to LL21A variable weighting is summed to obtain the second variable item of real part, the second variable item of real part and second
Real constant item is added gained multinomial, that is, real part deviation post nonlinear function；
The building of imaginary part current location nonlinear function specifically:
Imaginary part current location nonlinear function is by LL32The function and LL of the value containing different parameters of a difference power33A difference
The variable of the value containing different parameters of power is constituted；It is power with the nonlinear strength that each function is born, takes LL31A function weighted sum obtains
Imaginary part first function item, to remaining (LL32-LL31) a function weighted sum obtains imaginary part second function item；It is born with each variable
Nonlinear strength is power, to LL33A variable weighting is summed to obtain the first variable item of imaginary part；Imaginary part first function item subtracts imaginary part first
The product of variable item and imaginary part second function item, and third real constant item is added, gained multinomial, that is, imaginary part current location is non-thread
Property function；
The building of imaginary part deviation post nonlinear function specifically:
Imaginary part deviation post nonlinear function is by LL41The variable of the value containing different parameters of a difference power is constituted；With each change
The negative nonlinear strength of amount is power, to LL41A variable weighting is summed to obtain the second variable item of imaginary part, the second variable item of imaginary part and the 4th
Real constant item is added gained multinomial, that is, imaginary part deviation post nonlinear function；
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 amount that position offset, that is, position number increases or decreases, state translational movement, that is, variable
The amount that state value increases or decreases；The parameter of variable includes power, position number, position offset and the state translation of variable
Amount；
LL12、LL13、LL21、LL32、LL33、LL41It is the integer greater than 0, value sets itself as needed；LL11For not
Greater than LL12Positive integer, LL31For no more than LL32Positive integer, LL11And LL31It is worth sets itself as needed；
S5 parameter initialization and use pseudo-random number sequence or the sequence of real numbers being made of different real numbers initialization are effective
The state value of plural variable and extension plural variable, effective plural variable and extension plural variable are referred to as
For plural variable；
S6 is right respectively using multiple groups different real part current location nonlinear function and real part deviation post nonlinear function
The current state value of effective plural variable real part current location and deviation post is acted on respectively, obtains real part threshold；
Using the different imaginary part current location nonlinear function of multiple groups and imaginary part deviation post nonlinear function respectively to effective plural shape
The current state value of state variable imaginary part current location and deviation post is acted on respectively, obtains imaginary part threshold；It is based on diffusion
Number carries out real part threshold and imaginary part threshold comprising add, subtracting, the hybrid operation of at least one of multiplication and division operation respectively, logical
It crosses state iteration and generates the plural pseudo-random number sequence being distributed at any time；
This step further comprises:
Real part state iteration, specifically:
It is acted on respectively using the different real part current location nonlinear function of multiple groups and real part deviation post nonlinear function
Current state value in effective plural variable real part current location and deviation post obtains real part current location nonlinear function
Value and real part deviation post Nonlinear function；
Using diffusion coefficient as weight, respectively to real part current location Nonlinear function and the non-linear letter of real part deviation post
Numerical value is weighted and averaged, and obtains the first real part threshold and the second real part threshold；To the first real part threshold and the second real part
Threshold is added, is subtracted each other, is multiplied or division operation, obtains the first hybrid operation threshold；
Take part real part current location Nonlinear function that arithmetic average is asked to obtain the first average value；Remaining real part current location
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining real part current location that company multiplies value
The quantity of functional value obtains the second average value, and the first average value and the second average value subtract each other
Take part real part deviation post Nonlinear function that arithmetic average is asked to obtain third average value, remaining real part deviation post
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining real part deviation post that company multiplies value
The quantity of functional value obtains the 4th average value, and third average value and the 4th average value subtract each otherFirst hybrid operation threshold
Divided byObtain the current position state value of k+1 moment real part；
Imaginary part state iteration, specifically:
It is acted on respectively using the different imaginary part current location nonlinear function of multiple groups and imaginary part deviation post nonlinear function
Current state value in effective plural variable imaginary part current location and deviation post obtains imaginary part current location nonlinear function
Value and imaginary part deviation post Nonlinear function；
Using diffusion coefficient as weight, respectively to imaginary part current location Nonlinear function and the non-linear letter of imaginary part deviation post
Numerical value is weighted and averaged, and obtains the first imaginary part threshold and the second imaginary part threshold；To the first imaginary part threshold and the second imaginary part
Threshold is added, is subtracted each other, is multiplied or division operation, obtains the second hybrid operation threshold；
Take part imaginary part current location Nonlinear function that arithmetic average is asked to obtain the 5th average value, remaining imaginary part current location
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining imaginary part current location that company multiplies value
The quantity of functional value obtains the 6th average value, and the 5th average value and the 6th average value are added
Take part imaginary part deviation post Nonlinear function that arithmetic average is asked to obtain the 7th average value, all imaginary part deviation posts
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by all imaginary part deviation posts that company multiplies value
The quantity of functional value obtains the 8th average value, and the 7th average value and the 8th average value subtract each otherSecond hybrid operation threshold
Divided byObtain k+1 moment imaginary part current position state value；
The plural pseudo-random number sequence or rearranged plural number puppet that S7 is obtained using the currently active plural variable
Random number sequence modifies G0The state value of middle extension plural variable, or weight is carried out each other using these state values modified
Combination nova arrangement；Then, next effective plural variable in G is read, step S6 is executed；When effective plural shapes all in G
The equal completion status iteration of state variable executes step S8；
S8 is respectively from G1And G2The real part pseudo-random number sequence and imaginary part puppet that middle correlated variables tap extraction is distributed at any time with
Machine Number Sequence is denoted as the first pseudo random number and the second pseudo random number respectively, will be each in the first pseudo random number and the second pseudo random number
Real number random number is in chronological order respectively compared with corresponding a reference value, and if more than a reference value, real number pseudo random number value 1 is no
Then value 0 can be obtained the first pseudo noise code and the second pseudo noise code of binaryzation；A reference value is respectively to describe the first pseudorandom
The statistic of each real number pseudo random number median size in several and the second pseudo random number；
First pseudo noise code and the second pseudo noise code are carried out two He of mould by S9, obtain ranging code.
Optimize plural number pseudo-noise code generator G in above-mentioned steps S2 and obtain minimum series, specifically:
2.1 initialization take N=1；
2.2 from G1N number of state variable in take k1A variable tap, from remainder (N-k1) a state variable takes k2A variable
Tap seeks number of combinations It respectively indicates from G1N, N-k1K is taken respectively in a state variable1、
k2The number of combinations of a variable tap, k1And k2It is [1, N] interior integer；
2.3 from G2N number of state variable in take k3A variable tap, from remainder (N-k3) a state variable takes k4A variable
Tap seeks number of combinations It respectively indicates from G2N, N-k3K is taken respectively in a state variable3、
k4The number of combinations of a variable tap, k3And k4It is [1, N] interior integer；
2.4 from G1N number of state variable in take k5A variable tap, from G2N number of state variable in take k6A variable is taken out
Head seeks number of combinations It respectively indicates from G1And G2N number of state variable in take k respectively5、k6A change
Measure the number of combinations of tap, k5And k6It is [1, N] interior integer；
2.5 seek the sum of 2.2~2.4 gained number of combinations of sub-step, if the sum of number of combinations is greater than satellite navigation system used and defends
Star sum, current N value are minimum series, are terminated；Otherwise, sub-step 2.6 is executed；
N is added 1 by 2.6, executes sub-step 2.2.
In step S5, effective plural variable is initialized using pseudo-random number sequence and extends the shape of plural variable
State value, specifically:
Two linear pseudo-noise code generators are constructed respectively, are denoted as the first linear pseudo-noise code generator and the second linear puppet
Random code generator；
Respectively drive the first linear pseudo-noise code generator and the second linear pseudo-noise code generator and from related register
Tap output third pseudo noise code and the 4th pseudo noise code；
Different real numbers are respectively set to by third pseudo noise code and the 4th pseudo noise code 0,1, are changed into third pseudorandom
Several and the 4th pseudo random number, gained pseudo-random number sequence will guarantee that effective plural variable is in chaos working condition, if institute
Pseudo-random number sequence cannot be guaranteed chaos working condition, the diffusion coefficient and function and/or variable of nonlinear function need to be adjusted
Nonlinear strength；
Third pseudo random number and the 4th pseudo random number are the original state of the respectively real and imaginary parts of plural variable
Value.
In step S5, effective plural variable and extension plural number are initialized using the sequence of real numbers being made of different real numbers
The state value of state variable, the sequence of real numbers will guarantee that effective plural variable is in chaos working condition, if real number
Sequence cannot be guaranteed chaos working condition, need to adjust nonlinear function diffusion coefficient and function and/or variable it is non-linear
Intensity.
A reference value described in step S8 can be used weight method, ranking method or median method and obtain.
In the weight method, the first pseudo random number and the second pseudo random number are carried out respectively:
The sum of all real number pseudo random numbers in pseudo random number are sought, is denoted as and is worth；Each real number pseudo random number square is removed respectively
With with sum again after value i.e. weight and, by weight and as a reference value.
In the ranking method, a reference value of the first pseudo random number and the second pseudo random number be respectively the first pseudo random number and
The median that each real number pseudo random number obtains after size sorts in second pseudo random number, using median as a reference value.
In the median method, the first pseudo random number and the second pseudo random number are carried out respectively: asking real in pseudo random number
The maximum value of number pseudo random numbers and the difference of minimum value, are denoted as that extreme value is poor, 0.5 times of extreme value difference and minimum value are summed, by gained
A reference value is used as with value.
Two, a kind of multi-state variable space-time chaos plural number pseudo-noise code generator realizes system, comprising:
(1) plural pseudo-noise code generator constructs module, and for constructing plural pseudo-noise code generator G, G is by a series of multiple
Number state variable { x (i)+y (i) j } constitute (i=1,2 ... N, N are positive integer), { x (i) } and { y (i) } respectively constitute puppet with
Machine code generator G1、G2, { x (i) } and { y (i) } is a series of state variable that sequences are arranged and intercoupled respectively；I indicates multiple
Number state variable position serial number, include in G is effective plural variable；
(2) optimization module obtains minimum series for optimizing plural pseudo-noise code generator G, using minimum series as G's
Series N；
(3) expansion module, for being constituted in the forward or backward supplement extension plural variable of effective plural variable
G0, extension plural variable number is more than or equal to preset position offset maximum value；
(4) nonlinear function constructs module, for constructing nonlinear function；
The nonlinear function building module further comprises (4-1) real part current location nonlinear function building mould
Block, (4-2) real part deviation post nonlinear function building module, (4-3) imaginary part current location nonlinear function building module and
(4-4) imaginary part deviation post nonlinear function constructs module；Wherein:
(4-1) real part current location nonlinear function constructs module, for the building of real part current location nonlinear function,
Specifically:
Real part current location nonlinear function is by LL12The function and LL of the value containing different parameters of a difference power13A difference
The variable of the value containing different parameters of power is constituted；It is power with the nonlinear strength that each function is born, takes LL11A function weighted sum obtains
Real part first function item, to remaining (LL12-LL11) a function weighted sum obtains real part second function item；It is born with each variable
Nonlinear strength is power, to LL13A variable weighting is summed to obtain the first variable item of real part；The first variable item of real part is divided by real part second
Function item, in addition real part first function item and the first real constant item, gained multinomial, that is, real part current location nonlinear function；
(4-2) real part deviation post nonlinear function constructs module, for the building of real part deviation post nonlinear function,
Specifically:
Real part deviation post nonlinear function is by LL21The variable of the value containing different parameters of a difference power is constituted；With each change
The negative nonlinear strength of amount is power, to LL21A variable weighting is summed to obtain the second variable item of real part, the second variable item of real part and second
Real constant item is added gained multinomial, that is, real part deviation post nonlinear function；
(4-3) imaginary part current location nonlinear function constructs module, for the building of imaginary part current location nonlinear function,
Specifically:
Imaginary part current location nonlinear function is by LL32The function and LL of the value containing different parameters of a difference power33A difference
The variable of the value containing different parameters of power is constituted；It is power with the nonlinear strength that each function is born, takes LL31A function weighted sum obtains
Imaginary part first function item, to remaining (LL32-LL31) a function weighted sum obtains imaginary part second function item；It is born with each variable
Nonlinear strength is power, to LL33A variable weighting is summed to obtain the first variable item of imaginary part；Imaginary part first function item subtracts imaginary part first
The product of variable item and imaginary part second function item, and third real constant item is added, gained multinomial, that is, imaginary part current location is non-thread
Property function；
(4-4) imaginary part deviation post nonlinear function constructs module, for the building of imaginary part deviation post nonlinear function,
Specifically:
Imaginary part deviation post nonlinear function is by LL41The variable of the value containing different parameters of a difference power is constituted；With each change
The negative nonlinear strength of amount is power, to LL41A variable weighting is summed to obtain the second variable item of imaginary part, the second variable item of imaginary part and the 4th
Real constant item is added gained multinomial, that is, imaginary part deviation post nonlinear function；
Wherein, the parameter of function includes for working frequency, the power of function, the range value of function, the phase of function, position
Serial number, position offset and state translational movement；The amount that position offset, that is, position number increases or decreases, state translational movement become
The amount that amount state value increases or decreases；The parameter of variable includes power, position number, position offset and the state translation of variable
Amount；
LL12、LL13、LL21、LL32、LL33、LL41It is the integer greater than 0, value sets itself as needed；LL11For not
Greater than LL12Positive integer, LL31For no more than LL32Positive integer, LL11And LL31It is worth sets itself as needed；
(5) initialization module, for parameter initialization and using pseudo-random number sequence or the reality being made of different real numbers
The state value of the effective plural variable of number sequence row initialization and extension plural variable, effective plural variable and extension
Plural variable is referred to as plural variable；
(6) state iteration module, real part current location nonlinear function and real part bits of offset for using multiple groups different
Nonlinear function is set respectively to carry out the current state value of effective plural variable real part current location and deviation post respectively
Effect, obtains real part threshold；Use the different imaginary part current location nonlinear function of multiple groups and the non-linear letter of imaginary part deviation post
Number respectively acts on the current state value of effective plural variable imaginary part current location and deviation post respectively, obtains imaginary part
Threshold；Based on diffusion coefficient, real part threshold and imaginary part threshold are carried out respectively comprising adding, subtracting, at least one of multiplication and division
The hybrid operation of operation generates the plural pseudo-random number sequence being distributed at any time by state iteration；
The state iteration module further comprises (6-1) real part state iteration module and (6-2) imaginary part state iteration
Module, in which:
(6-1) real part state iteration module, for the real part current location nonlinear function and real part different using multiple groups
Deviation post nonlinear function is respectively acting on the current state of effective plural variable real part current location and deviation post
Value, obtains real part current location Nonlinear function and real part deviation post Nonlinear function；
Using diffusion coefficient as weight, respectively to real part current location Nonlinear function and the non-linear letter of real part deviation post
Numerical value is weighted and averaged, and obtains the first real part threshold and the second real part threshold；To the first real part threshold and the second real part
Threshold is added, is subtracted each other, is multiplied or division operation, obtains the first hybrid operation threshold；
Take part real part current location Nonlinear function that arithmetic average is asked to obtain the first average value；Remaining real part current location
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining real part current location that company multiplies value
The quantity of functional value obtains the second average value, and the first average value and the second average value subtract each other
Take part real part deviation post Nonlinear function that arithmetic average is asked to obtain third average value, remaining real part deviation post
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining real part deviation post that company multiplies value
The quantity of functional value obtains the 4th average value, and third average value and the 4th average value subtract each otherFirst hybrid operation threshold
Divided byObtain the current position state value of k+1 moment real part；
(6-2) imaginary part state iteration module, for the imaginary part current location nonlinear function and imaginary part different using multiple groups
Deviation post nonlinear function is respectively acting on the current state of effective plural variable imaginary part current location and deviation post
Value, obtains imaginary part current location Nonlinear function and imaginary part deviation post Nonlinear function；
Using diffusion coefficient as weight, respectively to imaginary part current location Nonlinear function and the non-linear letter of imaginary part deviation post
Numerical value is weighted and averaged, and obtains the first imaginary part threshold and the second imaginary part threshold；To the first imaginary part threshold and the second imaginary part
Threshold is added, is subtracted each other, is multiplied or division operation, obtains the second hybrid operation threshold；
Take part imaginary part current location Nonlinear function that arithmetic average is asked to obtain the 5th average value, remaining imaginary part current location
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by remaining imaginary part current location that company multiplies value
The quantity of functional value obtains the 6th average value, and the 5th average value and the 6th average value are added
Take part imaginary part deviation post Nonlinear function that arithmetic average is asked to obtain the 7th average value, all imaginary part deviation posts
Nonlinear function connects again after being multiplied respectively with corresponding diffusion coefficient and multiplies, and it is non-linear divided by all imaginary part deviation posts that company multiplies value
The quantity of functional value obtains the 8th average value, and the 7th average value and the 8th average value subtract each otherSecond hybrid operation threshold
Divided byObtain k+1 moment imaginary part current position state value；
(7) judgment module, for the plural pseudo-random number sequence using the acquisition of the currently active plural variable or through weight
The plural pseudo-random number sequence modification G newly arranged0The state value of middle extension plural variable, or the state using these modifications
Value carries out reconfiguring arrangement each other；Then, next effective plural variable in G is read, state iteration module is gone to；
When effective equal completion status iteration of plural variable all in G, binarization block is gone to；
(8) binarization block, for respectively from G1And G2The real part pseudorandom that middle correlated variables tap extraction is distributed at any time
Number Sequence and imaginary part pseudo-random number sequence, are denoted as the first pseudo random number and the second pseudo random number respectively, by the first pseudo random number and
Each real number random number is in chronological order respectively compared with corresponding a reference value in second pseudo random number, if more than a reference value, real number
Pseudo random number value 1, otherwise value 0, can be obtained the first pseudo noise code and the second pseudo noise code of binaryzation；A reference value difference
For the statistic of each real number pseudo random number median size in the first pseudo random number of description and the second pseudo random number；
(9) ranging code obtains module, for the first pseudo noise code and the second pseudo noise code to be carried out two He of mould, obtains ranging
Code.
The ranging code that the present invention can overcome current binary system pseudo-noise code generator and real number pseudo-noise code generator to generate
Existing all technological deficiencies, are greatly improved the performance of satellite navigation system ranging code, the present invention provides satellite navigation systems
High-performance ranging code of new generation realizes technology, can substitute present satellites navigation system ranging code realization technology comprehensively.
Compared to the prior art, the invention has the advantages that and the utility model has the advantages that
Compared with existing binary system pseudo-noise code generator and real number pseudo-noise code generator realize technology, the present invention is excellent
Point are as follows:
(1) it can get the high pseudo noise code of complexity
The present invention is acted on iteratively using the nonlinear function collective effect of a variety of different powers by hybrid operation
Pseudo noise code is generated, pseudo noise code complexity obtained is high.
(2) it can get the pseudo noise code of high safety
The present invention generates the high pseudo noise code of complexity, can fully ensure that the high security of pseudo noise code obtained.
(3) pseudo noise code code length is not limited by series
What the present invention generated is the random code being distributed at any time, maximum code length and pseudo-noise code generator using series without
It closes, and up to endless.
(4) series spy is few
The present invention passes through plural implementation, can reduce pseudo-noise code generator series to greatest extent.
(5) it is more to generate pseudo noise code pattern
The parameter that pseudorandom pattern of the present invention includes by real number, the nonlinear function of initialization plural variable states value
The parameters precision of such as diffusion coefficient, nonlinear strength determines that these parameters precisions of receiving end are at least 10-5, the present invention can produce
Raw pseudorandom pattern is at least 105×LA, L is code length.
(6) pseudo noise code characteristic can be adjusted at any time
It, can be by adjusting the diffusion coefficient, function and/or change of nonlinear function when pseudo noise code pseudo-randomness is destroyed
The nonlinear strength of amount restores.
Detailed description of the invention
Fig. 1 is the idiographic flow schematic diagram of the method for the present invention；
Fig. 2 is the flow diagram of binaryzation in embodiment；
Fig. 3 is in embodiment according to G1The real part pseudo-random number sequence that 2nd state variable obtains；
Fig. 4 is the pseudo noise code for using weight method binaryzation to obtain real part pseudo-random number sequence shown in Fig. 3；
Fig. 5 is the pseudo noise code for using ranking method binaryzation to obtain real part pseudo-random number sequence shown in Fig. 3；
Fig. 6 is the pseudo noise code for using median method binaryzation to obtain real part pseudo-random number sequence shown in Fig. 3；
Fig. 7 is in embodiment using the ranging code obtained after weight method binary conversion treatment；
Fig. 8 is in embodiment using the ranging code obtained after ranking method binary conversion treatment；
Fig. 9 is in embodiment using the ranging code obtained after median method binary conversion treatment；
Figure 10 is the auto-correlation function of ranging code shown in Fig. 7；
Figure 11 is the auto-correlation function of ranging code shown in Fig. 8；
Figure 12 is the auto-correlation function of ranging code shown in Fig. 9；
Figure 13 is present system structural schematic diagram.
Specific embodiment
Below with reference to embodiment, with forward direction extension plural variable, present invention is further described in detail, according to Fig. 1
Realize that steps are as follows:
S1: building plural number pseudo-noise code generator G, G N number of plural variable { x by sequentially arranging and intercoupling
(i)+y (i) j } it constitutes, wherein { x (i) } constitutes pseudo-noise code generator G1, { y (i) } constitutes pseudo-noise code generator G2, i table
It giving instructions in reply number state variable position serial number, N indicates series, i=1,2 ... N.
S2: optimizing the series of plural pseudo-noise code generator G, obtains minimum series Noptimal, after minimum series optimizes
Plural pseudo-noise code generator G series.
This step further comprises sub-step:
2.1 initialization take N=1；
2.2 from G1N number of state variable in take k1A variable tap, from remainder (N-k1) state variable takes k2A variable is taken out
Head seeks number of combinations It respectively indicates from G1N, N-k1K is taken respectively in a state variable1、k2
The number of combinations of a variable tap；
2.3 from G2N number of state variable in take k3A variable tap, from remainder (N-k3) variable takes k4A variable tap, is asked
Number of combinations It respectively indicates from G2N, N-k3K is taken respectively in a state variable3、k4A variable
The number of combinations of tap；
2.4 from G1N number of state variable in take k5A variable tap, from G2N number of state variable in take k6A variable is taken out
Head seeks number of combinations It respectively indicates from G1And G2N number of state variable in take k respectively5、k6A change
Measure the number of combinations of tap；
2.5 seek the sum of 2.2~2.4 gained number of combinations of sub-step, if the sum of number of combinations is greater than Chinese Beidou satellite navigation system
System satellite number 35, current N value are minimum series Noptimal, terminate；Otherwise, step 2.6 is executed；
N is added 1 by 2.6, executes step 2.2.
Minimum series is listed in table 1, i.e., plural variable number in table.
Minimum series obtained by 1 embodiment of table
| Plural variable number (unit: a) | 1 | 2 | 3 | 4 |
| Variable tap number of combinations (unit: kind) | 1 | 11 | 51 | >100 |
S3: by before effective plural variable to extension IMA plural variable simultaneously constitutes G0。
The specific embodiment of this step:
It takes than max { Il2,1,Il4,1Big numerical value IM, that is, extend plural variable number；Forward direction extends IMA plural
Variable, the I extendedMA plural variable is denoted as G0。G0Middle extension plural variable position serial number i is [1, IM] in it is whole
It counts, effectively plural variable position serial number i is [1+I in GM,N+IM] in integer.
To guarantee subsequent computational efficiency, IMValue should not be too large, and generally take the positive integer no more than 10.
S4: the nonlinear function of effective plural variable real and imaginary parts and the building of state iterative formula are acted on.
The building of (4-1) nonlinear function:
In formula (1)~(2):
K indicates discrete time coordinate；
I indicates the position number of effective plural variable, i=1+IM,...,N+IM；
xk(i)、yk(i) the current position state value of i-th of plural variable real and imaginary parts of k moment is respectively indicated；
Respectively indicate act on i-th of effective plural variable real part current location,
Real part deviation post, imaginary part current location, imaginary part deviation post nonlinear function,Separately including frequency is f0's
Sin function and cos function, phase are respectively ΦkWith Φ 'k；
l1、l2、l3、l4Respectively indicate nonlinear functionSerial number, that is,Expression acts on
The l of i effective plural variable real part current locations1A nonlinear function,Expression acts on i-th of effectively plural shape
The l of state variable real part deviation post2A nonlinear function,Expression acts on i-th of effective plural variable imaginary part and works as
The l of front position3A nonlinear function,Expression acts on the l of i-th of effective plural variable imaginary part deviation post4
A nonlinear function；
ll11、ll12Respectively indicate nonlinear functionIn include function, variable serial number；ll2Indicate nonlinear functionIn include variable serial number；ll31、ll32Respectively indicate nonlinear functionIn include function, variable serial number；ll4
Indicate nonlinear functionIn include variable serial number；
For nonlinear functionMiddle power isLl11A function,For nonlinear functionMiddle power
ForLl31A function；
LL12、LL13Respectively indicate nonlinear functionMiddle functionWith variable xk(i) quantity；LL21Indicate non-linear
FunctionMiddle variable xk(i) quantity；LL32、LL33Respectively indicate nonlinear functionMiddle functionWith variable yk(i) number
Amount；LL41Indicate nonlinear functionMiddle variable yk(i) quantity；
LL12、LL13、LL21、LL32、LL33、LL41It is the integer greater than 0, value sets itself as needed；LL11For not
Greater than LL12Positive integer, LL31For no more than LL32Positive integer, value sets itself as needed；
Respectively indicate k moment nonlinear functionIn ll11A functionWith ll12A variableNonlinear strength；
Indicate k moment nonlinear functionIn ll2A variableNonlinear strength；
Respectively indicate k moment nonlinear functionIn ll31A functionWith ll32A change
AmountNonlinear strength；
For k moment nonlinear functionIn ll4A variableNonlinear strength；
For nonlinear functionIn ll12A variable xk(i) power；For nonlinear functionIn ll2A change
Measure xk(i) power；For nonlinear functionIn ll32A variable yk(i) power；For nonlinear functionIn
ll4A variable yk(i) power；
Respectively nonlinear functionThe real constant for including, value
Sets itself as needed.
Above-mentioned nonlinear strength is real number, and above-mentioned secondary number formulary is the positive integer greater than 1, and value is voluntarily set as needed
It is fixed；
(4-2) state iterative formula
If the code length that plural number pseudo-noise code generator G is generated is L, the nonlinear function of different parameters value is contained using multiple groups
The current state value of effective plural variable real and imaginary parts current location and deviation post is acted on, is power with diffusion coefficient
Value by the inclusion of add, subtract, the hybrid operation mode of at least one of multiplication and division operation, divide the generation of plural variable at any time
The plural pseudo random number of cloth.
Pass through the state for i-th of effective plural variable (k+1) moment real and imaginary parts that hybrid operation mode obtains
Value xk+1(i)、yk+1(i) state iterative formula is as follows:
In formula (3):
Wherein:
Indicate the k moment theThe state value of the real part of a effective plural variable of deviation post；
Indicate the k moment theThe state value of the imaginary part of a effective plural variable of deviation post；
L12、L22、L32、L42Respectively indicate nonlinear functionQuantity, L12、L22、L32、L42?
For the integer greater than 0, value sets itself as needed；
L11、L21、L31、L41Respectively it is not more than L12、L22、L32、L42Positive integer, value sets itself as needed；
WithRespectively indicate k moment nonlinear functionWithThe position offset for including；
Respectively indicate nonlinear functionWithThe state translational movement for including；
Respectively k moment nonlinear functionWithThe state translational movement for including；
Respectively indicate the k moment
Nonlinear functionThreshold；
Respectively nonlinear functionIn the diffusion system at k moment
Number is real number；
S5: the initialization of parameter and plural variable states value.
(5-1) parameter initialization
Parameter includes plural variable number N (i.e. series), code length L, extension plural variable number IM, nonlinear function
Quantity, nonlinear function parameter, diffusion coefficient.Nonlinear function is the multinomial for including different power functions and/or variable
Form, parameter further comprise the parameter and real constant item of function and/or variable, function and/or variable, wherein the parameter of function
For working frequency, the power of function, the range value of function, the phase of function, position number, position offset, state translational movement,
The parameter of variable is power, position number, position offset, the state translational movement of variable.
In the present invention, function power is the integer greater than 1, no maximum requirement；The diffusion coefficient of each nonlinear function is equal
For real number.Nonlinear functionSeparately including frequency is f0Sin function and cos function, phase is respectively ΦkWith Φ 'k,
The range value of all functions is numerical value 1, f0For groundwork frequency；Nonlinear functionIt respectively acts on effectively multiple
The nonlinear function of number state variable real and imaginary parts current location.
The present embodiment takes plural variable number N=3, i.e., minimum series N according to table 1optimal, code length L=511, expand
Open up plural variable quantity IM=5.
Act on the quantity L of the nonlinear function of effective plural variable real and imaginary parts current location11=1, L12=
3、L31=1, L32=3；It acts on before effective plural variable real and imaginary parts to the quantity of the nonlinear function of deviation post
L21=L41=1, L22=L42=2.Act on each nonlinear function of real part current locationThe function highest power difference for including
For 2,4,6 powers, each nonlinear function of imaginary part current location is acted onThe function highest power for including is respectively 3,5,7
Power；Act on each nonlinear function before real part to deviation postThe function highest power for including is respectively 2,3 powers, is made
For before imaginary part to each nonlinear function of deviation postThe function highest power for including is respectively 2,3 powers.
The diffusion coefficient for acting on effective plural variable real and imaginary parts current location is respectively Wherein, 1st, 2,3 nonlinear function when respectively indicating moment 1Diffusion coefficient,
Respectively moment 1 when the 1st, 2,3 nonlinear functionDiffusion coefficient.
It acts on before effective plural variable real and imaginary parts and is respectively to the diffusion coefficient of deviation postWherein,Point
Not Wei moment 1 when the 1st, 2 nonlinear functionDiffusion coefficient,Respectively moment 1 when the 1st, 2 it is non-
Linear functionDiffusion coefficient.
Act on the nonlinear function of effective plural variable real part current locationThe highest power function for including exists
Preceding LL11In a function, corresponding nonlinear strength is respectively Remaining LL11The corresponding nonlinear strength of -1 function is 0；LL12-LL11A function is corresponding non-linear
Intensity sets itself；LL13The corresponding nonlinear strength of a variable is 0；Its forward direction deviation post highest power variable is corresponding non-
Linear intensity is respectivelyRemaining LL21- 1 variable is corresponding non-linear strong
Degree is 0；Act on the nonlinear function of current location imaginary partThe highest power function for including is in preceding LL31It is corresponding in a function
Nonlinear strength be respectively Remaining LL31- 1 letter
The corresponding nonlinear strength of number is 0；LL32-LL31The corresponding nonlinear strength sets itself of a function；LL33A variable is corresponding
Nonlinear strength is 0；Its corresponding nonlinear strength of forward direction deviation post highest power variable is respectivelyRemaining LL41The corresponding nonlinear strength of -1 variable is 0；Wherein,Respectively moment 1 when the 1st, 2,3 nonlinear functionHighest power function is corresponding non-
Linear intensity,Respectively moment 1 when the 1st, 2 nonlinear functionHighest power variable is corresponding non-thread
Property intensity,Respectively moment 1 when the 1st, 2,3 nonlinear functionHighest power function pair
The nonlinear strength answered,Respectively moment 1 when the 1st, 2 nonlinear functionHighest power variable pair
The nonlinear strength answered.
Effective plural variable real part current location the 1st, 2,3 nonlinear functionThe real constant for including is respectivelyWithThe state translational movement for including is respectively 1.50005；
Effective plural variable imaginary part current location the 1st, 2,3 nonlinear functionThe real constant for including is respectivelyWithThe state translational movement for including is 2.50005；Have
To the 1st, 2 nonlinear function of deviation post before effect plural variable real partThe real constant for including is respectivelyWithThe state translational movement for including is respectively 1.00001；Effective plural becomes
To the 1st, 2 nonlinear function of deviation post before amount imaginary partThe real constant for including is respectivelyWithThe state translational movement for including is 2.00002.
Forward location offset is respectively Respectively indicate the moment
1st, 2 nonlinear function when 1Corresponding forward location offset,When respectively indicating moment 1 the 1st, 2 it is non-
Linear functionCorresponding forward location offset；Moment 1 is initial time.
f0=1.023Hz, nΦ=16, monodrome section is [0, pi/2], f0Groundwork frequency is provided for the ranging code of generation,
nΦFor to nonlinear functionPhase of the sin function and cos function separately included in monodrome section carries out initial
Phase angle is divided equally.
In the present invention, the setting of the nonlinear strength of the diffusion coefficient and function and/or variable of each nonlinear function will be protected
Card system is in chaos working condition, and series is set as minimum series, the nonlinear function of corresponding current locationPhase
The product that place value is divided for the position number and monodrome section angle of effective plural variable, monodrome section, that is, variable and its function
It is worth one-to-one section, it is n that phase, which always divides number,Φ。
The initialization of (5-2) plural variable states value
Pseudo-random number sequence can be used or sequence of real numbers is initialized.Plural shape is initialized according to pseudo-random number sequence
The state value of state variable can be used linear shift register by related register tap and obtain pseudo noise code；According to real number
Sequence initialization plural variable needs to guarantee that effective plural variable works in chaos state.
The specific implementation process using pseudo-random number sequence initialization plural variable is provided below:
5.1: two 11 grades of linear displacement pseudo-noise code generators are constructed according to multinomial (4) and (5):
Z1(X)=1+X+X7+X8+X9+X10+X11 (4)
Z2(X)=1+X+X2+X3+X4+X5+X8+X9+X11 (5)
5.2: being respectively driven with initial value [1 110111011 1] and [1 010101010 1]
Two 11 grades of linear displacement pseudo-noise code generators generate pseudo noise code, and feed back pseudo noise code from 11 grades of linear displacements respectively
Last shift register of generator extracts 530 pseudo noise codes, is denoted as third pseudo noise code and the 4th pseudo noise code respectively；
5.3: numerical value 1 and 0 in third pseudo noise code is replaced respectively using 0.00085 and 0.00035, using 0.00065 He
0.00015 replaces numerical value 1 and 0 in the 4th pseudo noise code respectively, obtains third pseudo random number and the 4th pseudo random number；With third puppet with
Machine number is real part, the 4th pseudo random number is that imaginary part constitutes plural pseudo-random number sequence；
5.4: the real and imaginary parts of plural variable obtain respectively with position number multiplied by 0.1 and 0.2 with step 5.3
The real and imaginary parts values of 530 plural pseudo-random number sequences be separately summed, i.e. real and imaginary parts of plural variable
Initial value.
S6: using formula (3) state is carried out to effective plural variable each in plural pseudo-noise code generator G respectively and changed
In generation, generates 530 seconds plural pseudo-random number sequence being distributed at any time.
S7: to G0In respectively extend plural variable, with its position number multiplied by 530 obtained with step S6 after 0.001
The value of real part of the plural pseudo-random number sequence in position is added, and obtains the real part state value of the extension plural variable；Again with its position-order
It number is added multiplied by after 0.003 with the imaginary values for 530 plural pseudo-random number sequences that step S6 is obtained, obtains the extension plural
The imaginary part state value of variable realizes G0The modification of middle extension plural variable states value；Then it reads next effectively multiple in G
Number state variable executes step S6.All effective equal completion status iteration of plural variable in G execute step S8.
S8: delay 5 seconds, to avoid initial non-chaotic oscillator signal, first respectively from G12nd, the 3rd variable extraction length
The real part pseudo-random number sequence of one group of time state distribution in 511 seconds, is shown in Fig. 3, and binaryzation, sees Fig. 4~6, and mould two and obtains
First pseudo noise code；Again from G21st variable and the 3rd variable extract the imaginary part pseudorandom of 511 seconds time states of length distribution
Number Sequence, binaryzation and mould two and be the second pseudo noise code, binaryzation process is shown in Fig. 2；
It, will be from G in the present invention1The real part pseudo-random number sequence of middle correlated variables tap extraction being distributed at any time is denoted as
One pseudo random number, will be from G2The imaginary part pseudo-random number sequence of middle correlated variables tap extraction being distributed at any time be denoted as second it is pseudo- with
Machine number.The binaryzation of pseudo random number can be used weight method, ranking method or median method and obtain:
(1) weight method
First pseudo random number and the second pseudo random number are carried out respectively:
The sum of all real number pseudo random numbers in pseudo random number are sought, is denoted as and is worth；By square difference of each real number pseudo random number
Sum again divided by with after value, be denoted as weight and, by weight and as a reference value.By each real number pseudo random number respectively with a reference value ratio
Compared with being greater than a reference value, real number pseudo random number value 1, otherwise value 0, to obtain the first pseudo random number and the second pseudo random number
Corresponding binaryzation pseudo noise code, is denoted as the first pseudo noise code and the second pseudo noise code.
(2) ranking method
First pseudo random number and the second pseudo random number are carried out respectively:
Real number pseudo random number in pseudo random number in magnitude order, is taken into median, using median as a reference value.It will be each
Real number pseudo random number is respectively compared with a reference value, if more than a reference value, real number pseudo random number value 1, and otherwise value 0, thus
To the first pseudo random number and the corresponding binaryzation pseudo noise code of the second pseudo random number, it is denoted as the first pseudo noise code and the second pseudorandom
Code.
(3) median method
First pseudo random number and the second pseudo random number are carried out respectively:
The maximum value and minimum value of real number pseudo random number in pseudo random number, then the difference of maximizing and minimum value are asked, is denoted as
Extreme value is poor, and 0.5 times of extreme value difference and minimum value are summed, extreme value and value are denoted as, and regard extreme value and value as a reference value.By each real number
Pseudo random number is greater than a reference value, real number pseudo random number value 1, otherwise value 0 compared with a reference value, thus obtain first it is pseudo- with
The binaryzation pseudo noise code that machine number and the second pseudo random number are answered, is denoted as the first pseudo noise code and the second pseudo noise code.
S9: carrying out two He of mould for the first pseudo noise code and the second pseudo noise code, forms the pseudo noise code that code length is 511,
That is navigation satellite ranging code, is shown in Fig. 7~9.
The pseudo-randomness of ranging code is evaluated, evaluation result is shown in Table 2~7 and Figure 10~12.Pseudo-randomness can be used
Balance, the distance of swimming and autocorrelation evaluation: numerical value 0 and 1 accounts for the percentage of sum in balance, that is, ranging code, ideally
0 and 1 accounts for the 50% of sum respectively；The percentage of the total number of runs of distance of swimming Zhan of different length, ideal feelings in the distance of swimming, that is, ranging code
The percentage of each total number of runs of length runs Zhan is under conditionWherein, a indicates run length；Autocorrelation, that is, ranging
The δ function characteristic of code auto-correlation function.
The balance of ranging code shown in 2 Fig. 7 of table
| Numerical value | Percent of total (%) |
| 0 | 50.68 |
| 1 | 49.32 |
The balance of ranging code shown in 3 Fig. 8 of table
| Numerical value | Percent of total (%) |
| 0 | 46.18 |
| 1 | 53.82 |
The balance of ranging code shown in 4 Fig. 9 of table
| Numerical value | Percent of total (%) |
| 0 | 48.53 |
| 1 | 51.47 |
The distance of swimming of ranging code shown in 5 Fig. 7 of table
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 45.65 |
| 2 | 24.35 |
| 3 | 14.78 |
| 4 | 3.91 |
| 5 | 2.61 |
| 6 | 1.30 |
The distance of swimming of ranging code shown in 6 Fig. 8 of table
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 44.28 |
| 2 | 34.69 |
| 3 | 13.65 |
| 4 | 4.06 |
| 5 | 2.21 |
| 6 | 1.11 |
The distance of swimming of ranging code shown in 7 Fig. 9 of table
| Run length | Account for distance of swimming sum percentage (%) |
| 1 | 55.08 |
| 2 | 26.56 |
| 3 | 15.74 |
| 4 | 1.64 |
| 5 | 0.66 |
| 6 | 0.33 |
Claims (8)
1. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method, characterized in that include:
S1 building plural number pseudo-noise code generator G, G be made of a series of plural variables { x (i)+y (i) j }, { x (i) } and
{ y (i) } respectively constitutes pseudo-noise code generator G1、G2, { x (i) } and { y (i) } is that a series of sequences are arranged and intercoupled respectively
State variable；I expression plural variable position serial number, i=1,2 ... N, N are positive integer；Include in G is effectively multiple
Number state variable；
S2 optimizes plural number pseudo-noise code generator G and obtains minimum series, using minimum series as the series N of G；
S3 constitutes G in the forward or backward supplement extension plural variable of effective plural variable0, extend plural variable
Number is greater than or equal to preset position offset maximum value；
S4 building is respectively acting on effective plural variable real and imaginary parts current location and deviation post current state value
Nonlinear function, specifically:
Act on real part or imaginary part current location and deviation post current state value nonlinear function be one group of identical power or
The function and/or compound function and/or variable of the value containing different parameters of different powers, with the non-thread of function and/or compound function
Property intensity be power, a part of function and/or compound function weighted sum are obtained to the first function item of real part or imaginary part, by remaining letter
Several and/or compound function weighted sum obtains the second function item of real part or imaginary part；It is power with the nonlinear strength of variable, by one
Variation per minute weighted sum obtains the first variable item of real part or imaginary part, and second that surplus variable weighted sum is obtained real part or imaginary part becomes
Quantifier；By first function item, second function item, the first variable item and the second variable item by predetermined manner carry out include add, subtract,
The hybrid operation of at least two operations in multiplication and division, along with corresponding real constant item, gained multinomial acts on real part or imaginary part
The nonlinear function of current location and deviation post current state value；
S5 parameter initialization and use pseudo-random number sequence or the sequence of real numbers being made of different real numbers initialization effectively plural number
The state value of state variable and extension plural variable, effective plural variable and extension plural variable are referred to as multiple
Number state variable；
S6 is using the multiple groups nonlinear function of building respectively to effective plural variable real and imaginary parts current location and offset
The current state value of position is acted on, and real part threshold and imaginary part threshold are obtained；It is respectively real part with different diffusion coefficients
The weighting coefficient of threshold and imaginary part threshold includes by predetermined manner respectively to real part threshold and imaginary part threshold
Add, subtract, in multiplication and division at least two operations hybrid operation, the plural pseudo random number that is distributed at any time is generated by state iteration
Sequence；
The plural pseudo-random number sequence or rearranged plural pseudorandom that S7 uses the currently active plural variable to obtain
Number Sequence modifies G0The state value of middle extension plural variable, or group again is carried out each other using these state values modified
Close arrangement；Then, next effective plural variable in G is read, step S6 is executed；When effective plurals all in G become
Equal completion status iteration is measured, step S8 is executed；
S8 is respectively from G1And G2The real part pseudo-random number sequence and imaginary part pseudo random number that middle correlated variables tap extraction is distributed at any time
Sequence is denoted as the first pseudo random number and the second pseudo random number respectively, by each real number in the first pseudo random number and the second pseudo random number
Random number is in chronological order respectively compared with corresponding a reference value, and if more than a reference value, otherwise real number pseudo random number value 1 takes
Value 0 can be obtained the first pseudo noise code and the second pseudo noise code of binaryzation；A reference value be respectively describe the first pseudo random number and
The statistic of each real number pseudo random number median size in second pseudo random number；
Choosing the right moment for first pseudo noise code and the second pseudo noise code is as follows with length:
Delay 5 seconds, to avoid initial non-chaotic oscillator signal, first respectively from G12nd, the 3rd variable extract length 511 seconds one
The real part pseudo-random number sequence of group time state distribution, and binaryzation, and mould two and obtain the first pseudo noise code；Again from G21st
A variable and the 3rd variable extract the imaginary part pseudo-random number sequence of 511 seconds time states of length distribution, binaryzation and mould two and are
Second pseudo noise code；
First pseudo noise code and the second pseudo noise code are carried out two He of mould by S9, obtain ranging code.
2. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
Step S2 specifically:
From G1N number of state variable in take k1A variable tap, from G1Remaining state variable take k2A variable tap, by default side
Formula is combined to obtain the first number of combinations；From G2N number of state variable in take k3A variable tap, from G2Remaining state variable take
k4A variable tap is combined to obtain the second number of combinations by predetermined manner；From G1N number of state variable in take k5A variable tap,
From G2N number of state variable in take k6A variable tap is combined to obtain third number of combinations by predetermined manner；If the first combination
The sum of number, the second number of combinations and third number of combinations are greater than the satellite navigation system population of satellite used, and the number of combinations is corresponding effectively
Plural variable number is the plural pseudo-noise code generator minimum series obtained.
3. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
In step S4, nonlinear function is by the function containing different parameters value of one group of identical power or different power and/or multiple
Close function and/or variable, as the nonlinear strength and real constant structure of each function and/or compound function and/or variable weight
At multinomial, the parameter that wherein function and/or compound function include be working frequency, the power of function, function range value,
The phase of function, position number, position offset, state translational movement, the parameter that variable includes are power, the position-order of variable
Number, position offset, state translational movement, real number precision is 10 in parameter-5, and the real and imaginary parts of plural variable are non-linear
The phase difference for the function and/or compound function that function includes is orthogonal, using nonlinear strength as the letter of weight in nonlinear function
The operation for 1. adding, subtract, multiplying or removing is carried out between several and/or compound function and/or variable, 2. presets hybrid operation.
4. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
In step S5, effective plural variable is initialized using pseudo-random number sequence and extends the state of plural variable
Value, specifically:
Two linear pseudo-noise code generators are constructed respectively, are denoted as the first linear pseudo-noise code generator and the second linear pseudorandom
Code generator；
Respectively drive the first linear pseudo-noise code generator and the second linear pseudo-noise code generator and from related register tap
Export third pseudo noise code and the 4th pseudo noise code；
Be respectively set to different real numbers in third pseudo noise code and the 4th pseudo noise code 0,1, be changed into third pseudo random number and
4th pseudo random number, gained pseudo-random number sequence will guarantee that effective plural variable is in chaos working condition, if gained is pseudo-
Random number sequence cannot be guaranteed chaos working condition, need to adjust the diffusion coefficient and function and/or compound function of nonlinear function
And/or the nonlinear strength of variable；
Third pseudo random number and the 4th pseudo random number are the initial state value of the respectively real and imaginary parts of plural variable.
5. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
It is effectively multiple using the sequence of real numbers being made of different real numbers or the initialization of the sequence of real numbers by rearranging in step S5
The state value of number state variable and extension plural variable, the sequence of real numbers will guarantee that effective plural variable is in
Chaos working condition needs to adjust the diffusion coefficient of nonlinear function if gained sequence of real numbers cannot be guaranteed chaos working condition
With the nonlinear strength of function and/or compound function and/or variable.
6. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
In step S6,1. the nonlinear function threshold progress for acting on effective plural variable real and imaginary parts is added, subtracted,
Multiply or remove operation, 3. using diffusion coefficient as the operation for add, subtract, multiplying or removing of weight and 4. 2. preset hybrid operation, with expansion
Dissipate the default hybrid operation that coefficient is weight.
7. multi-state variable space-time chaos plural number pseudo-noise code generator implementation method as described in claim 1, it is characterized in that:
A reference value described in step S8 is obtained using weight method, ranking method or median method；
In weight method, the first pseudo random number and the second pseudo random number are carried out respectively:
The sum of all real number pseudo random numbers in pseudo random number are sought, is denoted as and is worth；Each real number pseudo random number square respectively divided by with
Sum again after value i.e. weight and, by weight and as a reference value；
In ranking method, a reference value of the first pseudo random number and the second pseudo random number is respectively each real number pseudo random number by big float
The median obtained after sequence, using median as a reference value；
In median method, the first pseudo random number and the second pseudo random number are carried out respectively:
The difference for asking the maximum value of real number pseudo random number and minimum value in pseudo-random number sequence, is denoted as that extreme value is poor, by 0.5 times of extreme value
Difference is summed with minimum value, regard gained and value as a reference value.
8. multi-state variable space-time chaos plural number pseudo-noise code generator realizes system, characterized in that include:
(1) plural pseudo-noise code generator constructs module, and for constructing plural pseudo-noise code generator G, G is by a series of plural shapes
State variable { x (i)+y (i) j } is constituted, and { x (i) } and { y (i) } respectively constitutes pseudo-noise code generator G1、G2, { x (i) } and { y
(i) } a series of state variable for arranging for sequences and intercoupling respectively；I expression plural variable position serial number, i=1,
2 ... N, N are positive integer；Include in G is effective plural variable；
(2) optimization module obtains minimum series for optimizing plural pseudo-noise code generator G, using minimum series as the series of G
N；
(3) expansion module, for constituting G in the forward or backward supplement extension plural variable of effective plural variable0, expand
It opens up plural variable number and is greater than or equal to preset position offset maximum value；
(4) nonlinear function constructs module, is respectively acting on effective plural variable real and imaginary parts present bit for constructing
The nonlinear function with deviation post current state value is set, specifically:
Act on real part or imaginary part current location and deviation post current state value nonlinear function be one group of identical power or
The function and/or compound function and/or variable of the value containing different parameters of different powers, with the non-thread of function and/or compound function
Property intensity be power, a part of function and/or compound function weighted sum are obtained to the first function item of real part or imaginary part, by remaining letter
Several and/or compound function weighted sum obtains the second function item of real part or imaginary part；It is power with the nonlinear strength of variable, by one
Variation per minute weighted sum obtains the first variable item of real part or imaginary part, and second that surplus variable weighted sum is obtained real part or imaginary part becomes
Quantifier；By first function item, second function item, the first variable item and the second variable item by predetermined manner carry out include add, subtract,
The hybrid operation of at least two operations in multiplication and division, along with corresponding real constant item, gained multinomial acts on real part or imaginary part
The nonlinear function of current location and deviation post current state value；
(5) initialization module, for parameter initialization and using pseudo-random number sequence or the real number sequence being made of different real numbers
The state value of the effective plural variable of row initialization and extension plural variable, effective plural variable and extension plural number
State variable is referred to as plural variable；
(6) state iteration module, for use the multiple groups nonlinear function of building respectively to effective plural variable real part and
The effect of the current state value of imaginary part current location and deviation post, obtains real part threshold and imaginary part threshold；With different diffusions
Coefficient is respectively the weighting coefficient of real part threshold and imaginary part threshold, is pressed respectively to real part threshold and imaginary part threshold default
Mode carry out include add, subtract, in multiplication and division at least two operations hybrid operation, generated by state iteration and be distributed at any time
Plural pseudo-random number sequence；
(7) judgment module, plural pseudo-random number sequence for being obtained using the currently active plural variable or is again arranged
The plural pseudo-random number sequence of column modifies G0The state value of middle extension plural variable, or the state value phase using these modifications
It carries out reconfiguring arrangement between mutually；Then, next effective plural variable in G is read, state iteration module is executed；Work as G
In all effective equal completion status iteration of plural variable, execute binarization block；
(8) binarization block, for respectively from G1And G2The real part pseudorandom number sequence that middle correlated variables tap extraction is distributed at any time
Column and imaginary part pseudo-random number sequence, are denoted as the first pseudo random number and the second pseudo random number, by the first pseudo random number and second respectively
Each real number random number is in chronological order respectively compared with corresponding a reference value in pseudo random number, if more than a reference value, real number puppet with
Machine number value 1, otherwise value 0, can be obtained the first pseudo noise code and the second pseudo noise code of binaryzation；A reference value is respectively to retouch
State the statistic of each real number pseudo random number median size in the first pseudo random number and the second pseudo random number；
Choosing the right moment for first pseudo noise code and the second pseudo noise code is as follows with length:
Delay 5 seconds, to avoid initial non-chaotic oscillator signal, first respectively from G12nd, the 3rd variable extract length 511 seconds one
The real part pseudo-random number sequence of group time state distribution, and binaryzation, and mould two and obtain the first pseudo noise code；Again from G21st
A variable and the 3rd variable extract the imaginary part pseudo-random number sequence of 511 seconds time states of length distribution, binaryzation and mould two and are
Second pseudo noise code；(9) ranging code obtains module, for the first pseudo noise code and the second pseudo noise code to be carried out two He of mould, obtains
Ranging code.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610608799.1A CN106301753B (en) | 2016-07-29 | 2016-07-29 | Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610608799.1A CN106301753B (en) | 2016-07-29 | 2016-07-29 | Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106301753A CN106301753A (en) | 2017-01-04 |
| CN106301753B true CN106301753B (en) | 2019-03-19 |
Family
ID=57663091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610608799.1A Expired - Fee Related CN106301753B (en) | 2016-07-29 | 2016-07-29 | Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106301753B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109861743B (en) * | 2019-03-10 | 2021-03-23 | 西安电子科技大学 | Data safety transmission method based on Beidou short message positioning information transmission protocol |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101621287A (en) * | 2009-08-10 | 2010-01-06 | 太原理工大学 | Real random code generating device based on chaos laser and generating method thereof |
| CN104614732A (en) * | 2014-12-19 | 2015-05-13 | 上海电机学院 | Chaotic light ranging technology based anti-collision signal processing system and method |
| US9465582B1 (en) * | 2010-11-18 | 2016-10-11 | The Boeing Company | Significant random number generator |
-
2016
- 2016-07-29 CN CN201610608799.1A patent/CN106301753B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101621287A (en) * | 2009-08-10 | 2010-01-06 | 太原理工大学 | Real random code generating device based on chaos laser and generating method thereof |
| US9465582B1 (en) * | 2010-11-18 | 2016-10-11 | The Boeing Company | Significant random number generator |
| CN104614732A (en) * | 2014-12-19 | 2015-05-13 | 上海电机学院 | Chaotic light ranging technology based anti-collision signal processing system and method |
Non-Patent Citations (2)
| Title |
|---|
| 一种利用混沌同步的测控系统测距方法;余金峰等;《电讯技术》;20150730;第55卷(第7期);全文 |
| 基于重叠的P码频域直接捕获法研究;肖模艳，何志彪;《湖北第二师范学院学报》;20130228;第30卷(第2期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106301753A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101394266B (en) | Method for generating variable parameter chaos signal and chaos secret communication system | |
| CN103698785B (en) | Based on the satellite-signal cycle-slip detection and repair method of passing through order limit learning machine | |
| CN100562762C (en) | The method and system that is used for positioning global positioning system correlated peak signal | |
| CN102081359A (en) | DSP Builder-based time-varying delay hyperchaos digital circuit design method and circuit | |
| CN103197912B (en) | A kind of uniform random number generation method for deep space communication protocol code | |
| CN109725294A (en) | Radar array sparse optimization method based on improved adaptive GA-IAGA | |
| CN106291616B (en) | Space-time chaos vector pseudo-noise code generator offset carrier modulator approach and system | |
| CN108886380A (en) | Safe lane measurement | |
| CN106301753B (en) | Multi-state variable space-time chaos plural number pseudo-noise code generator realization method and system | |
| CN107632323A (en) | Radar echo signal processing method for Ionospheric F detection | |
| CN109039522B (en) | Chaos sequence-based method for optimizing balance of spreading codes | |
| CN106911384A (en) | Suitable for the building method of the two-dimentional light address code of visible ray OCDMA communications | |
| Barratt | A free group functor for stable homotopy | |
| Sonnenschein | Algebraic approach to solve tt dilepton equations | |
| CN106817326B (en) | Pseudo code blind estimation method of multi-user long-period short code direct sequence spread spectrum signal | |
| Li et al. | A chaos-based pseudo random number generator using timing-based reseeding method | |
| CN106772476B (en) | Multivariable space-time chaos satellite navigation high-performance ranging code method and system | |
| CN1697366B (en) | Controllable method for generating chaos cipher sequence of natural number, and core of digit chip | |
| CN101807215A (en) | Method for designing chip for real-time decomposition of mixed pixel of hyper-spectral image | |
| CN108833156A (en) | A kind of appraisal procedure and system of the simulation performance index for power telecom network | |
| CN114489567A (en) | Parallel pseudo-random code generator and method | |
| CN107547106A (en) | The long short code method of estimation of the long short code direct sequency-code division multiple access signal of multi tate | |
| Rui et al. | A novel ranging code based on improved logistic map chaotic sequences | |
| CN102799413B (en) | A kind of random digit generation method and device | |
| Jingye et al. | GPS C/A code signal simulation based on MATLAB |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190319 |