CN105093190B - A kind of design method of orthogonal phase coded signal - Google Patents
A kind of design method of orthogonal phase coded signal Download PDFInfo
- Publication number
- CN105093190B CN105093190B CN201510418718.7A CN201510418718A CN105093190B CN 105093190 B CN105093190 B CN 105093190B CN 201510418718 A CN201510418718 A CN 201510418718A CN 105093190 B CN105093190 B CN 105093190B
- Authority
- CN
- China
- Prior art keywords
- phase
- coded signal
- signal
- autocorrelation
- coded
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/003—Bistatic radar systems; Multistatic radar systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/42—Diversity systems specially adapted for radar
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S2013/0236—Special technical features
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides a kind of design method of orthogonal phase coded signal, can submit the design efficiency of orthogonal phase coded signal, obtain the orthogonal phase coded signal with relatively low autocorrelation performance.This method comprises the following steps:M phase-coded signal is produced, M phase-coded signal has identical code length N, and each phase-coded signal has autocorrelation;The normalized autocorrelation side lobe peak of M phase-coded signal is calculated respectively, is ranked up according to its corresponding normalized autocorrelation side lobe peak by order from small to large;P phase-coded signal before being chosen from M phase-coded signal after sequence;According to the normalized crosscorrelation peak value of any two phase-coded signal in preceding P phase-coded signal, the L phase-coded signals with relatively low normalized crosscorrelation peak value are chosen from preceding P phase-coded signal.
Description
Technical field
The present invention relates to technical field of radar communication, more particularly to a kind of design method of orthogonal phase coded signal.
Background technology
Modern radar requirement still has relatively good target detection, identification and ability of tracking in complicated electromagnetic environment.Example
Such as, MIMO radar can according to system requirements flexible arrangement transmitting antenna and reception antenna, can from multiple angle detection targets and
Target echo has independence, can improve radar to the detection of target, identification performance (see document:Fishler E,Haimovich
A,Blum R S,et al.Spatial diversity in radar models and detection performance
[J] .IEEE Transactions on Signal Processing, 2006,54 (3):823-838.).
Distributive array coherent synthesis radar can complete transmitting-receiving coherent synthesis and obtain N3Maximum system gain (N is single
First radar number), it has power big, and detection accuracy is high with disposing the technical advantages such as maneuverability (see document:Gao H W,Cao
Z,Lu Y B.Basic study and principle validate of distributed aperture
coherence-synthetic radar[C]∥Proc.of the 12th Chinese Radar Conference,2012:
129-134.)。
These radars do not interfere with each other to ensure transmission signal and independent target information are obtained from echo, are required for hair
Penetrate orthogonal waveforms.It can be seen that, the design of orthogonal waveforms is one of key technology of modern radar.
Orthogonal waveforms are often realized with the mode of phase code, traditional design mode be mostly based on signal autocorrelation secondary lobe and
The average energy of cross-correlation or the total optimization design philosophy of peak value, and the multi-purpose modern times such as simulated annealing or genetic algorithm
Intelligent algorithm optimizes search, obtains more satisfactory orthogonal signalling (see document:Deng H.Synthesis of Binary
Sequences with Good Autocorrelation and Crosscorrelation Properties by
Simulated Annealing[J].IEEE Transactions on Aerospace and electronic systems,
1996,32(1):98-107.Liu B,He Z,Zeng J,et al.Polyphase orthogonal code design
For MIMO radar systems [C] ∥ IEEE International Conference on radar, CIE'06,
Shanghai:IEEE,2006:1-4.).When code length is longer or signal way is more, the optimized algorithm in above-mentioned document is calculated
Complexity is very big, greatly constrains the design efficiency and performance of radar emission signal.
The content of the invention
For the shortcoming of above-mentioned technology, it is an object of the invention to provide a kind of design side of orthogonal phase coded signal
Method, it is possible to increase the design efficiency of orthogonal phase coded signal, obtains a kind of quadrature phase with relatively low autocorrelation performance and compiles
Code signal.
The technical thought that the present invention is realized is:Code element forecast model is proposed, the model is first passed through and obtains with low auto-correlation
Property sequence sets, then choose the sequence with minimum average cross-correlation peak value from the sequence sets of low autocorrelation, so obtain
Orthogonal phase coded signal there is relatively low autocorrelation performance.
To reach above-mentioned purpose, embodiments of the invention, which are adopted the following technical scheme that, to be achieved.
A kind of design method of orthogonal phase coded signal, methods described comprises the following steps:
Step 1, M phase-coded signal is produced, M phase-coded signal has identical code length N, and each phase is compiled
Code signal has autocorrelation, wherein, M and N are natural number;
Step 2, the normalized autocorrelation side lobe peak of the M phase-coded signal is calculated respectively, by the M phase
Encoded signal is ranked up according to its corresponding normalized autocorrelation side lobe peak by order from small to large, after being sorted
M phase-coded signal;
Step 3, P phase-coded signal before being chosen from M phase-coded signal after the sequence, before calculating is described
The normalized crosscorrelation peak value of any two phase-coded signal in P phase-coded signal, wherein, P is natural number, and P<M;
Step 4, according to the normalized crosscorrelation peak of any two phase-coded signal in the preceding P phase-coded signal
Value, chooses the L phase-coded signals with relatively low normalized crosscorrelation peak value, institute from the preceding P phase-coded signal
It is the orthogonal phase coded signal designed to state the L phase-coded signals with relatively low normalized crosscorrelation peak value, wherein,
L is natural number, and L<P.
The features of the present invention and further it is improved to:
(1) step 1 specifically includes following sub-step:
(1a) given initiation sequence length is k0Phase-coded signal
Wherein, x1, x2...,Wherein, D=2,3,4 ...,It is the phase number of phase-coded signal, initiation sequence length k to constitute set Φ, D0's
Value is met:N>k0>logDM;
(1b) is with initial phase encoded signalFor initiation sequence, N-k is continuously predicted0It is secondary, obtain the phase that code length is N
Encoded signal
(1c) repeats sub-step (1a) to (1b) M times, obtains the M phase-coded signals with autocorrelation, structure
Into signal collection.
Further, in the sub-step (1a),
Initiation sequence length k0=int (min (2logDM, N/2)), int () represents round computing.
The sub-step (1b) specifically includes following sub-step:
OrderK=k0, using mathematical modelingUse current sequence XkPredict xk+1, build
New sequence Xk+1=[Xk,xk+1] so that new sequence Xk+1Normalized autocorrelation peak value relative to current sequence XkIt is on a declining curve,
Wherein, f is cost function.
The mathematical modelingIn cost function f be specially appointing in following four cost function
Anticipate a kind of cost function:
①
②
③
④
Cost function 1., 2., 3., 4. inWithThe calculating formula of i-th of element be respectively:
I.e.It is XkOne-side autocorrelation, length is k-1,It is Xk+1One-side autocorrelation, length is
k。
(2) step 4 specifically includes following sub-step:
(4a) presses the normalized crosscorrelation peak value of any two phase-coded signal in the preceding P phase-coded signal
The dimension that order constitutes cross-correlation symmetrical matrix C, the cross-correlation symmetrical matrix C is P × P;
(4b) sorts every row element in the cross-correlation symmetrical matrix C by order from small to large respectively, is arranged
Matrix C after sequence ', by the Matrix C after sequence ' in the index that corresponds in cross-correlation matrix C of each element recorded, obtain
To index matrix I, the every corresponding L phase-coded signal of L element before row in index matrix I, to the L phase-coded signal
The normalized crosscorrelation peak value summation of middle any two phase-coded signal obtains P and value E, wherein, the index matrix I's
Dimension is P × P;
(4c) chooses minimum in described P and value E and is used as the orthogonal of final design with corresponding L phase-coded signal is worth
Phase-coded signal.
The present invention has advantages below compared with the design method of conventional orthogonal phase-coded signal:
(1) cost function relative to existing algorithm be mostly based on signal autocorrelation secondary lobe and cross-correlation average energy or
The total optimization design philosophy of peak value, the present invention proposes code element forecast model, does not take in, carries out every time from total optimization
Current best selection, is oriented prediction and obtains the signal collection with low NASP, substantially reduce computation complexity, improve
Design efficiency;
(2) prediction of present invention orientation obtains the signal collection with low NASP, and is searched further in the sequence sets with most
Low NASP signal so that the orthogonal phase coded signal of design has lower autocorrelation performance.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the design method schematic flow sheet of orthogonal phase coded signal provided in an embodiment of the present invention;
Fig. 2 is overall implementation process schematic diagram provided in an embodiment of the present invention;
Fig. 3 is utilization cost function f of the present invention1,f2,f3,f4The 1000 low auto-correlation phase-coded signal separately designed
Middle minimum NASP with code length change curve;
Fig. 4 is the normalized autocorrelation curve and cross-correlation for the four phase encoded signals that 3 code lengths that the present invention is designed are 400
Curve;
Fig. 5 be the average NASP of 3 orthogonal four phases encoded signals that is separately designed with the present invention and genetic algorithm and
Be averaged variation relation curves of the NCP with code length.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
For convenience, it is defined as below first:
NASP:Normalized autocorrelation side lobe peak.
NCP:Normalized crosscorrelation peak value.
The embodiment of the present invention provides a kind of orthogonal phase coded signal design method with relatively low autocorrelation performance, is used for
It is L to produce signal number, and code length is the N orthogonal phase coded signal with relatively low autocorrelation performance, as shown in figure 1, specifically
Comprise the following steps:
Step 1, M phase-coded signal is produced, M phase-coded signal has identical code length N, and each phase is compiled
Code signal has autocorrelation.
Wherein, M and N is natural number.
Step 1 specifically includes following sub-step:
(1a) given initiation sequence length is k0Initial phase encoded signal
Wherein, x1,x2,…,Wherein, D=2,3,4 ...,It is the phase number of the initial phase encoded signal, initiation sequence to constitute set Φ, D
Length k0Value meet:N>k0>logDM。
Exemplary, for four phase encoded signals, D=4, Φ={ 0, pi/2, π, 3 pi/2s }.Initiation sequence length k0Value
Meet:N>k0>logDM, can make k0=int (min (2logDM, N/2)), int () represents round computing.Example
Such as, as generation M=10000, during code length N=200 four phase encoded signal, initiation sequence length k0=int (min (13.3,
100))=13.
(1b) is with initial phase encoded signalFor initiation sequence, N-k is continuously predicted0It is secondary, obtain the phase that code length is N
Encoded signal
ByObtain XNN-k can be decomposed into0Secondary code element predicted operation, and operation can be described as using current sequence every time
XkPredict xk+1, build new sequence Xk+1=[Xk,xk+1] so that the NASP of new sequence is on a declining curve relative to former sequence, wherein
k0≤k≤N-1。
Therefore, X is obtainedNConcretely comprise the following steps:Initial value is first assigned, is madeK=k0,Pass through X every timekPredict xk+1Afterwards,
X is made againk=[Xk,xk+1], k=k+1 ("=" represents assignment operation), so continuous prediction N-k0It is secondary with regard to X can be obtainedN, and XNHave
Low NASP.
It can be seen that, use current sequence XkPredict xk+1It is crucial, its mathematical modeling is:
Wherein, f is cost function, to make Xk+1=[Xk,xk+1] NASP is on a declining curve and fall off rate faster, it is necessary to
Select a suitable cost function.The present invention provides four kinds of cost functions, is designated as f1,f2,f3,f4, it is specifically defined respectively formula
1. -4., cost function selection principle referring to emulation experiment 2 detailed description.
①
②
③
④
In the definition of four cost functions of the aboveThe calculating formula of i-th of element be respectively:
I.e.It is XkOne-side autocorrelation, length is k-1,It isOne-side autocorrelation, length is
k.OrderThenWithEqual length, be k, while there is following recurrence Relation:
Therefore,It can regard asNext state, due toEach element can be obtained with recursion, calculate
Amount is greatly reduced so that the present invention can quickly produce phase-coded signal.Formula 1. in mkActuallyMiddle maximum
Index.Formula 4. inRepresentVariance.f1Show that the NASP of current sequence is maximum in NextState
Reduction, referred to as decresting method;f2Show the autocorrelation peak secondary lobe of minimum renewal sequence, be referred to as minimum secondary lobe method;f3Show minimum
Change the autocorrelation sidelobe energy of renewal sequence, referred to as minimizing energy method;f4Show the side of minimum renewal sequence autocorrelation sidelobe
Difference, referred to as least variance method.Four kinds of methods are based on local optimum principle, and globally optimal solution is obtained by continuous recursion.
(1c) repeats sub-step (1a) to (1b) M times, obtains the M phase-coded signals with autocorrelation, structure
Into signal collection.
Specifically, the signal collection constituted in sub-step (1c) is designated as
Step 2, the normalized autocorrelation side lobe peak of the M phase-coded signal is calculated respectively, by the M phase
Encoded signal is ranked up according to its order of corresponding normalized autocorrelation side lobe peak from small to large, the M after being sorted
Individual phase-coded signal.
Step 2 chooses P (P from the M phase-coded signal that step 1 is produced<M, can typically makeIndividual tool
There is minimum NASP phase code signal:
Specifically, the NASP of M phase-coded signal of step 1 generation is calculated, by this M NASP that sort from small to large
And record the corresponding signal index l of each NASP.Solve l (l=1,2 ..., M) individual phase-coded signalNASP meter
Formula is:
Step 3, P phase-coded signal before being chosen from M phase-coded signal after the sequence, before calculating is described
The normalized crosscorrelation peak value of any two phase-coded signal in P phase-coded signal, wherein, P is natural number, and P<M.
Calculate the NCP of any two phase-coded signal in preceding P phase-coded signal, it is desirable to which NCP, which is obtained, first seeks two phases
The cross-correlation function of position encoded signal, specific formula for calculation is:
And p ≠ q
C(sp,sq, it is i) p-th of signal spWith q-th of signal sqIn moment i cross-correlation function, i is time index.Then
The NCP of p-th of phase-coded signal and q-th of phase-coded signal calculating formula is in P phase-coded signal before solving:
And p ≠ q.
Step 4, according to the normalized crosscorrelation peak of any two phase-coded signal in the preceding P phase-coded signal
Value, chooses the L phase-coded signals with relatively low normalized crosscorrelation peak value from the preceding P phase-coded signal.
The L phase-coded signals with relatively low normalized crosscorrelation peak value are the quadrature phase coding designed
Signal, wherein, L is natural number, and L<P.
Step 4 specifically includes following sub-step:
(4a) presses the normalized crosscorrelation peak value of any two phase-coded signal in the preceding P phase-coded signal
The dimension that order constitutes cross-correlation symmetrical matrix C, the cross-correlation symmetrical matrix C is P × P.
Because P phase-coded signal has low autocorrelation performance, so the main diagonal element of cross-correlation matrix is put
For 0, i.e., in this step without considering the autocorrelation of signal.The cross-correlation symmetrical matrix C concrete forms of structure are as follows:
(4b) sorts every row element in the cross-correlation symmetrical matrix C by order from small to large respectively, is arranged
Matrix C after sequence ', by the Matrix C after sequence ' in the index that corresponds in cross-correlation matrix C of each element recorded, obtain
To index matrix I, the every corresponding L phase-coded signal of L element before row in index matrix I, to the L phase-coded signal
The normalized crosscorrelation peak value summation of middle any two phase-coded signal obtains P and value E, wherein, the index matrix I's
Dimension is P × P.
Specifically, cross-correlation symmetrical matrix C every row element respectively by sorting from small to large, the matrix after being sorted
C ', and retain sequence after Matrix C ' in each element correspond to cross-correlation symmetrical matrix C in index, obtain index matrix
L element L phase-coded signal of correspondence before I, index matrix I often row, the index composition of this L phase-coded signal go to
Measure I1×pObtain the NCP's of any two phase-coded signal in this L phase-coded signal and value E, calculation formula is:
P and value are most obtained at last.
(4c) chooses minimum in described P and value E and is used as the orthogonal of final design with corresponding L phase-coded signal is worth
Phase-coded signal
Because cross-correlation symmetrical matrix C has P rows, so can finally obtain P NCP summing value.Minimum is chosen again to ask
Be worth corresponding L phase-coded signal, then this L phase-coded signal has good cross correlation, is used as final design
L × N orthogonal phase coded signal.
It is mentality of designing based on the present invention, what reference picture 2 can be designed that any number and code length has relatively low auto-correlation
The orthogonal phase coded signal of performance.
Emulation experiment 1:
Simulation parameter is set:Initiation sequence length k0=10, the number M=40000 of the phase-coded signal of generation, selection
Phase-coded signal number P=10000, the phase number D=4 of phase-coded signal, the orthogonal phase coded signal of final design
Number L=4, the code length N=40 of phase-coded signal.
4 quadri-phase code signals that the code length that table 1 lists design is 40, four signals are respectively designated as L1, L2, L3, L4,
0,1,2,3 represent phase 0, pi/2, π, 3 pi/2s respectively in table.Table 2 lists its auto-correlation and cross correlation properties, and leading diagonal is
Normalized autocorrelation side lobe peak NASP, remaining is normalized crosscorrelation peak value NCP.Table 3 lists the present invention and background technology
The design result of method (Deng ' s SA, Liu ' s GA) in listed bibliography, wherein the 2nd~5 row represent respectively this 4
Individual signal NASP and NCP maximum and average value.It is lower that contrast finds that the Polyphase Orthogonal Code signal that the present invention is designed has
NASP, and NCP is close to literature value.Illustrate that the present invention can be designed that the more excellent quadrature phase coding letter of autocorrelation performance
Number.
Emulation experiment 2:
Simulation parameter is set:Initiation sequence length k0=10, code length N=20~1600, phase number D=4, signal number
L=1.Fig. 3 is utilization cost function f of the present invention1,f2,f3,f4It is minimum in the 1000 low auto-correlation sample sequence separately designed
NASP is with the change curve of code length, and abscissa represents code length, and ordinate represents normalized autocorrelation side lobe peak (unit is dB),
It is all feasible and effective to understand these four cost functions, and choosing any of which can obtain with low autocorrelation performance
Signal.As code length N≤100, it is not very big to choose the corresponding design effect difference of four kinds of cost functions;As code length N > 100
When, minimizing energy method (f3) and least variance method (f4) corresponding design effect is substantially better than drop Peak Intensity Method (f1) and minimum secondary lobe
Method (f2), NASP maximum differences can reach 2dB or so.So to select suitable generation according to the code length of actually required orthogonal signalling
Valency function.
Emulation experiment 3:
Simulation parameter is set:Initiation sequence length k0=10, the signal number M=40000 that step 1 is produced, what step 2 was selected
Signal number P=4000, phase number D=4, signal number L=3.The average NASP for three sequences that the present invention is designed for-
25.59dB, average NCP is -21.07dB.The normalized autocorrelation curve of these three signals respectively as shown in Fig. 4 (a)-(c), its
Middle Fig. 4 (a)-(c) abscissa represents delay cell, and ordinate represents that auto-correlation normalizes amplitude (unit is dB), normalization
Cross-correlation curve is respectively as shown in Fig. 4 (d)-(f), and wherein Fig. 4 (d)-(f) abscissa represents delay cell, and ordinate is represented
Cross-correlation normalization amplitude (unit is dB), it is seen that the present invention can be designed that more satisfactory orthogonal signalling.
Emulation experiment 4:
The computation complexity expression formula of forecast model of the present invention is:Lost based on mixing
The computation complexity expression formula of propagation algorithm is:Blending heredity is set to calculate
Method parameter is L=3, D=4, G=1000, S=300, F=200, and common parameters are set:Phase number D=4, signal number L=
3.Emulation obtains the orthogonal signalling that code length N is respectively 100,200,400,800,1600.Emulate inventive algorithm, Reasonable adjustment M
With P value, as a result show that the average NASP of the orthogonal signalling of present invention design equal conditions is low relative to genetic algorithm
0.015, and average NCP quite (absolute difference is less than 0.005).Now record M average value is respectively 50,100,200,
400,400, P corresponding average values are respectively 10,20,100,80,100.To avoid contingency, above-mentioned emulation passes through 100
Secondary independent repeated trials.Each simulation parameter is substituted into, computation complexity is as shown in table 4, it is seen that inventive algorithm efficiency is far above mixed
Hybrid genetic algorithm (raising that averagely there are about 3 orders of magnitude), illustrates that the present invention can quickly design orthogonal phase coded signal.
Emulation experiment 5:
Parameter of the present invention is set to:M=4000~20000, P=500~2000;Genetic algorithm (genetic algorithm and
Iterative code selection algorithm) parameter is set to:S=1000~3000, G=300~400, F=200, wherein G are evolutionary generation, S
For Population Size, F is iterations.Common parameters are set:Phase number D=4, signal number L=3, code length N=100~
1000. this parameter settings can ensure that computation complexity of the present invention is less than genetic algorithm.For traditional pulse pressure sequence, its
NASP with sequence length N increase withRate reduction.The present invention depicts curveWith making comparisons, originally
The average NASP and average NCP of 3 orthogonal quadri-phase code sequences of design are invented with the variation relation curve of code length as shown in figure 5,
Abscissa represents code length, and ordinate represents normalization amplitude (unit is dB), works as N>After 300, average NASP and NCP is substantially pressed
According toRate reduction, it is seen that the average NASP of orthogonal signalling that the present invention is designed is averagely lower than genetic algorithm
2.3dB or so, and average NCP is approached.Illustrate that the orthogonal phase coded signal that the present invention is designed has lower autocorrelation performance.
The Polyphase Orthogonal Code (L=4, N=40, D=4) that the present invention of table 1 is designed
The NASP and NCP (L=4, N=40, D=4) of the road signal of table 2 four
L1 | L2 | L3 | L4 | |
L1 | 0.1346 | 0.1953 | 0.1953 | 0.2016 |
L2 | 0.1953 | 0.1250 | 0.2151 | 0.2305 |
L3 | 0.1953 | 0.2151 | 0.1275 | 0.2151 |
L4 | 0.2016 | 0.2305 | 0.2151 | 0.1275 |
The design results of table 3 (L=4, N=40, D=4)
Max(NASP) | Avg(NASP) | Max(NCP) | Avg(NCP) | |
Deng’s(SA) | 0.1820 | 0.1525 | 0.2121 | 0.1988 |
Liu’s(GA) | 0.1581 | 0.1471 | 0.2305 | 0.2078 |
The present invention | 0.1346 | 0.1286 | 0.2305 | 0.2088 |
The computation complexity of table 4 compares
The foregoing is only a specific embodiment of the invention, but protection scope of the present invention is not limited thereto, any
Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all be contained
Cover within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (4)
1. a kind of design method of orthogonal phase coded signal, it is characterised in that methods described comprises the following steps:
Step 1, M phase-coded signal is produced, M phase-coded signal has identical code length N, and each phase code letter
Number there is autocorrelation, wherein, M and N are natural number;
Wherein, step 1 specifically includes following sub-step:
(1a) given initiation sequence length is k0Phase-coded signal
Wherein,Wherein, D=2,3,4 ...,It is the phase number of the phase-coded signal, initiation sequence length to constitute set Φ, D
k0Value meet:N>k0>logDM;
(1b) is with phase-coded signalFor initiation sequence, N-k is continuously predicted0It is secondary, obtain the phase-coded signal that code length is N
Wherein, the sub-step (1b) specifically includes:
OrderUsing mathematical modelingUse current sequence XkPredict xk+1, build new sequence
Arrange Xk+1=[Xk,xk+1] so that new sequence Xk+1Normalized autocorrelation peak value relative to current sequence XkIt is on a declining curve, its
In, f is cost function;
(1c) repeats sub-step (1a) to (1b) M times, obtains the M phase-coded signals with autocorrelation, constitutes letter
Number collection;
Step 2, the normalized autocorrelation side lobe peak of the M phase-coded signal is calculated respectively, by the M phase code
Signal is ranked up according to its corresponding normalized autocorrelation side lobe peak by order from small to large, the M after being sorted
Phase-coded signal;
Step 3, P phase-coded signal before being chosen from M phase-coded signal after the sequence, calculates described preceding P
The normalized crosscorrelation peak value of any two phase-coded signal in phase-coded signal, wherein, P is natural number, and P<M;
Step 4, according to the normalized crosscorrelation peak value of any two phase-coded signal in the preceding P phase-coded signal,
The L phase-coded signals with relatively low normalized crosscorrelation peak value, the L are chosen from the preceding P phase-coded signal
The individual phase-coded signal with relatively low normalized crosscorrelation peak value is the orthogonal phase coded signal designed, wherein, L is
Natural number, and L<P.
2. the design method of orthogonal phase coded signal according to claim 1, it is characterised in that in the sub-step
In (1a),
Initiation sequence length k0=int (min (2logDM, N/2)), int () represents round computing.
3. the design method of orthogonal phase coded signal according to claim 1, it is characterised in that the mathematical modelingIn cost function f be specially any one cost function in following four cost function:
Cost function 1., 2., 3., 4. inWithThe calculating formula of i-th of element be respectively:
I.e.It is XkOne-side autocorrelation, length is k-1,It is Xk+1One-side autocorrelation, length is k.
4. the design method of orthogonal phase coded signal according to claim 1, it is characterised in that step 4 is specifically included
Following sub-step:
(4a) by the normalized crosscorrelation peak value of any two phase-coded signal in the preceding P phase-coded signal in order
The dimension for constituting cross-correlation symmetrical matrix C, the cross-correlation symmetrical matrix C is P × P;
(4b) sorts every row element in the cross-correlation symmetrical matrix C by order from small to large respectively, obtains after sequence
Matrix C ', by the Matrix C after sequence ' in each element correspond to cross-correlation matrix C in index recorded, obtain rope
Draw L element L phase-coded signal of correspondence before often being gone in matrix I, index matrix I, to appointing in the L phase-coded signal
The normalized crosscorrelation peak value summation of two phase-coded signals of meaning obtains P and value E, wherein, the dimension of the index matrix I
For P × P;
(4c) chooses minimum in described P and value E and is worth corresponding L phase-coded signal as the quadrature phase of final design
Encoded signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510418718.7A CN105093190B (en) | 2015-07-16 | 2015-07-16 | A kind of design method of orthogonal phase coded signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510418718.7A CN105093190B (en) | 2015-07-16 | 2015-07-16 | A kind of design method of orthogonal phase coded signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093190A CN105093190A (en) | 2015-11-25 |
CN105093190B true CN105093190B (en) | 2017-10-13 |
Family
ID=54574055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510418718.7A Active CN105093190B (en) | 2015-07-16 | 2015-07-16 | A kind of design method of orthogonal phase coded signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093190B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676219A (en) * | 2016-01-11 | 2016-06-15 | 桂林电子科技大学 | Quadrature-phase modulation-based MIMO radar three-dimensional imaging method |
CN108363047B (en) * | 2018-02-13 | 2022-05-03 | 西安电子科技大学 | Four-phase code Doppler compensation method based on sliding window interpolation method |
CN108562873B (en) * | 2018-03-30 | 2022-06-10 | 西安电子科技大学 | Radar communication integration signal implementation device based on code insertion fusion |
US11175394B2 (en) | 2019-03-27 | 2021-11-16 | Raytheon Company | Binary to quad phase generator |
CN110988816B (en) * | 2019-11-12 | 2023-09-05 | 西安电子科技大学 | Coherent MIMO radar emission waveform design method based on space-time phase coding |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813768A (en) * | 2009-02-19 | 2010-08-25 | 中国科学院电子学研究所 | Method for designing quadrature phase-coded signal based on hyperchaos mapping |
CN104198996A (en) * | 2014-08-30 | 2014-12-10 | 西安电子科技大学 | High bit rate and low range resolution of low sidelobe phase coded signal design method |
-
2015
- 2015-07-16 CN CN201510418718.7A patent/CN105093190B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813768A (en) * | 2009-02-19 | 2010-08-25 | 中国科学院电子学研究所 | Method for designing quadrature phase-coded signal based on hyperchaos mapping |
CN104198996A (en) * | 2014-08-30 | 2014-12-10 | 西安电子科技大学 | High bit rate and low range resolution of low sidelobe phase coded signal design method |
Non-Patent Citations (2)
Title |
---|
基于混合遗传算法的正交多相码波形优化设计;周畅;《中国优秀硕士学位论文全文数据库信息科技辑》;20110715(第07期);1-57 * |
面向任务的MIMO雷达波形设计方法研究;刘韵佛;《中国博士学位论文全文数据库信息科技辑》;20140115(第01期);1-113 * |
Also Published As
Publication number | Publication date |
---|---|
CN105093190A (en) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093190B (en) | A kind of design method of orthogonal phase coded signal | |
CN102323972B (en) | Method for managing phased array radar resource | |
CN102404756B (en) | Antenna parameter optimizing system based on mobile phone measurement report | |
CN101980298B (en) | Multi-agent genetic clustering algorithm-based image segmentation method | |
CN107290709A (en) | The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde | |
CN105426956A (en) | Ultra-short-period photovoltaic prediction method | |
CN104597434B (en) | Improve the multiframe coherent TBD methods of envelope shift compensation and Fourier Transform of Fractional Order | |
CN105870935A (en) | Clustering algorithm-based radial distribution network reactive power optimization method | |
CN109271695A (en) | Multiple target antenna design method neural network based | |
CN105225000A (en) | A kind of wind power probability model nonparametric probability method optimized based on Fuzzy Ordered | |
CN103207390B (en) | The approximate fractal detection method of target in the sea clutter of FRFT territory | |
Wang et al. | Improved prediction method of PV output power based on optimised chaotic phase space reconstruction | |
CN102789542B (en) | Load distribution method for hydrothermal power system based on multi-target distribution estimation | |
CN102426676B (en) | A kind of feature extracting method of intelligent power consumption strategy | |
CN110059376B (en) | Cylindrical surface conformal array based rapid optimization design method and system | |
CN105842666A (en) | Radar sub-array dividing optimization method based on difference algorithm | |
CN106488482B (en) | Wireless sensor network optimizing method based on multi-Agent evolutionary Algorithm | |
CN103267963A (en) | Multi-criterion radar emission signal design method based on phase modulation | |
CN113242066B (en) | Multi-cell large-scale MIMO communication intelligent power distribution method | |
CN102928827B (en) | Rapid dimension-reducing space-time self-adaption processing method based on PAST (Projection Approximation Subspace Tracking) | |
CN105137398B (en) | The pulse pressure filter optimization method of the anti-repeating jamming of radar based on genetic algorithm | |
CN105068052B (en) | The method for designing of one species zero correlation orthogonal phase coded signal | |
CN104008292A (en) | Rapid and accurate prediction method of ultra wide band electromagnetic pulse response of wide-band antenna | |
CN103954960B (en) | Based on the airborne early warn ing radar moving target detection method of clutter ridge steering vector | |
CN114745299A (en) | Non-invasive load monitoring method based on sequence delay reconstruction CSP convolutional neural network |
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 |