CN105093190B - A kind of design method of orthogonal phase coded signal - Google Patents

A kind of design method of orthogonal phase coded signal Download PDF

Info

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
Application number
CN201510418718.7A
Other languages
Chinese (zh)
Other versions
CN105093190A (en
Inventor
杨明磊
陈伯孝
张威
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201510418718.7A priority Critical patent/CN105093190B/en
Publication of CN105093190A publication Critical patent/CN105093190A/en
Application granted granted Critical
Publication of CN105093190B publication Critical patent/CN105093190B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/35Details of non-pulse systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/003Bistatic radar systems; Multistatic radar systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/42Diversity systems specially adapted for radar
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S2013/0236Special 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

A kind of design method of orthogonal phase coded signal
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.
CN201510418718.7A 2015-07-16 2015-07-16 A kind of design method of orthogonal phase coded signal Active CN105093190B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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