CN105044681B - The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high - Google Patents
The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high Download PDFInfo
- Publication number
- CN105044681B CN105044681B CN201510305913.9A CN201510305913A CN105044681B CN 105044681 B CN105044681 B CN 105044681B CN 201510305913 A CN201510305913 A CN 201510305913A CN 105044681 B CN105044681 B CN 105044681B
- Authority
- CN
- China
- Prior art keywords
- phase
- coded signal
- main lobe
- mismatched filter
- signal
- 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/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
Abstract
The invention belongs to Radar Technology field, a kind of mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high is disclosed.The optimization method comprises the following steps：Determine phase-coded signal s；According to the length N of phase-coded signal s and mismatched filter hh, build the object function on mismatched filter h；Build the bound for objective function on mismatched filter h；According to bound for objective function, the object function of mismatched filter h is solved；Determine mismatched filter.The present invention can reduce the distance side lobe level after phase-coded signal pulse compression.
Description
Technical field
The invention belongs to Radar Technology field, specifically a kind of low range resolution ratio phase-coded signal of code check high
Mismatched filter optimization method, for reducing the distance side lobe level after phase-coded signal mismatch filter.
Background technology
Pulse compression technique solves the contradiction between the operating distance of radar and range resolution ratio, is obtained in field of radar
It is widely applied.But pulse compression signal generally has distance side lobe higher, and distance side lobe higher is unfavorable for thunder
Effectively detected up to target, particularly the Faint target detection under multiple target or strong clutter background, the main lobe pole of weak signal target
Easily flooded by the distance side lobe of strong target echo and cause false dismissal, influence the detection performance of radar.For example, the distance side lobe of high-amplitude
False-alarm may be triggered, the small-power target that the distance side lobe of high-power target echo can be flooded at range cell nearby.
In order to reduce the distance side lobe of phase-coded signal, generally using the pulse compression signal of design low sidelobe come real
It is existing.In existing various phase-coded signals, Barker code signal is optimal biphase coded signal, and it has minimum peak side-lobe
Level, but, the Baud Length up to 13 of Barker code signal is unsatisfactory for actual application requirement.Compared with biphase coded signal,
Multiphase Barker code signal has lower peak sidelobe.Frank code, P1 yards, P2 yards, P3 yards and P4 yards be several typical cases
Polyphase code signal, they are all the polyphase code signals as derived from linear FM signal.P (n, k) code signal is by non-linear tune
Polyphase code signal derived from frequency signal, compared with the polyphase code signal as derived from linear FM signal, P (n, k) code signal has
Lower peak sidelobe.
Existing phase-coded signal, under many circumstances, its peak sidelobe is still unsatisfactory for actual application will
Ask, and existing phase-coded signal Timed automata be equal to phase-coded signal Baud Length, phase-coded signal
Bandwidth be equal to the inverse of code element time width.The phase-coded signal that time width is fixed can be entered by increasing its Baud Length
One step reduction distance side lobe.But for a wideband radar system, its signal transacting bandwidth is fixed, cannot now be passed through
It is simple to increase Baud Length to reduce distance side lobe, there is presently no what is proposed for this problem based on phase-coded signal
The method for designing of mismatched filter.
The content of the invention
For above-mentioned existing methods shortcoming, it is an object of the invention to propose a kind of mismatch filter of phase-coded signal
The optimization method of device, can reduce the distance side lobe level after phase-coded signal pulse compression.
Realizing the technical thought of the object of the invention is：The phase-coded signal of the low range resolution ratio of code check high is given, in phase
In the case that the bandwidth of position encoded signal is constant, with snr loss as constraints, led to minimizing the phase-coded signal
Cross the peak sidelobe of mismatched filter output and approach desired main lobe and be shaped as object function, optimization is actually needed
Mismatched filter.
To reach above-mentioned purpose, the embodiment of the present invention is adopted the following technical scheme that and is achieved.
A kind of mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high, it is characterised in that bag
Include following steps：
(1) phase-coded signal s is determined；
(2) the length N of setting mismatched filter hh, and according to the length of the phase-coded signal s and mismatched filter h
Nh, build the object function on the mismatched filter h；
(3) bound for objective function on the mismatched filter h is built；
(4) according to the bound for objective function, the object function of the mismatched filter h is solved；
(5) the mismatched filter h is determined.
The characteristics of technical solution of the present invention and it is further modified to：
Step (1) specifically includes following sub-step：
(1a) setting first phase encoded signal s1Baud Length N1；
(1b) encodes s according to the first phase1Baud Length N1, determine the Baud Length of the phase-coded signal s
Ns, wherein, Ns=b × N1, b is the increase multiple of the Baud Length of phase-coded signal s, and b values are integer；
(1c) sets main lobe width controlled quentity controlled variable M=fix (δ × b) of the phase-coded signal s, and δ is empirical coefficient, is taken
In the range of (0.5~1.0), b is the increase multiple of the Baud Length of phase-coded signal s to value, and b values are integer；
The Baud Length N of (1d) according to the phase-coded signal ss, the main lobe width control of the phase-coded signal s
Amount M, builds object function and constraints on the phase-coded signal s；
(1e) solves the phase-coded signal s according to the object function and the constraints.
Step (2) specifically includes following sub-step：
(2a) phase-coded signal s is ρ by the pulse compression vector after mismatched filter hc, wherein,
(2b) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, to each distance side lobe
Modulus of access, obtains distance side lobe modulus value vector ρ respectively：
Wherein, k ∈ [1, (Ns+Nh)/2-M-1], l ∈ [(Ns+Nh)/2+M+1, Ns+Nh-1]；
(2c) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, the actual main lobe of output
BmFor
Wherein, t ∈ [(Ns+Nh)/2-M, (Ns+Nh)/2+M]；
(2d) determines expectation main lobe bs of the phase-coded signal s by output after mismatched filter hm, wherein, desired master
Valve can be rectangle, the main lobe of sinc functions, or do not increase the main lobe of the phase-coded signal before he number.This example
Described in expect main lobe bmThe main lobe for being shaped as the first phase encoded signal, by the first phase encoded signal s1's
The interval equably discrete column vector for turning to 2 × M+1 dimensions of the corresponding functional value in main lobe part, and using the column vector as expectation
Main lobe bmValue；
(2e) is according to the distance side lobe modulus value vector ρ, the actual main lobe BmWith the expectation main lobe bm, build target
Function
Step (3) specifically includes following sub-step：
The symmetrical zero padding of phase-coded signal s head and the tail is extended to length for N by (3a)hPhase-coded signal
(3b) is N by lengthhPhase-coded signalSpin upside down and take conjugation and obtain matched filter
Wherein, ()*Represent the conjugation of ().
(3c) is according to the mismatched filter h, the matched filterObtain the mesh on the mismatched filter h
The constraints of scalar functions is
Further, step (1d) specifically includes sub-step：
A () calculates the kth of phase-coded signal s1Individual distance side lobe ρk1：
Wherein, k1=M+1, M+2 ..., (N are shifteds- 1), M is the main lobe width controlled quentity controlled variable of phase-coded signal s, Jk1For
Matrix is slided, matrix J is slidedk1Form be：
In formula, ()TTransposition is represented, 0 represents full null matrix, and I represents unit matrix, the dimension of the subscript representing matrix of 0 and I
Number；
The distance side lobe ρ of (b) according to phase-coded signal sk1, obtain the peak sidelobe PSL of distance side lobes, PSLs
=max | ρk1|；Peak sidelobe PSL according to distance side lobes, the object function for obtaining phase-coded signal s is：
PSLs, wherein, p is the phase vectors of phase-coded signal s；
C () determines the expectation main lobe b of phase-coded signal sn, wherein, desired main lobe can be rectangle, sinc functions
Main lobe, or do not increase the main lobe of the phase-coded signal before he number.Main lobe b is expected described in this examplenShape
It is the main lobe of the first phase encoded signal, by the first phase encoded signal s1The corresponding functional value area in main lobe part
Between the equably discrete column vector for turning to 2 × M+1 dimensions, take m-th element therein to the 1st element as expecting main lobe bn's
Value；
D () is according to the expectation main lobe bn, the constraints for obtaining phase-coded signal s is：
0≤p (k)≤2 π, k=1,2 ..., Ns
Wherein, bn(k2) represent and expect main lobe bn2 elements of kth, Jk2To slide matrix, k2=1,2 ..., M, M are shifted
It is the main lobe width controlled quentity controlled variable of phase-coded signal s, p (k) represents k-th element of phase vectors p, k=1,2 ..., Ns, γ
It is constant, γ represents that the main lobe of phase-coded signal s approaches expectation main lobe bnDegree, value be 0.01~0.2.
The present invention has the following advantages that compared with prior art.The mismatched filter of present invention design is ensureing phase code
While signal bandwidth is constant, phase-coded signal Baud Length is increased, therefore had during optimization design bigger
The free degree, so as to reduce further the distance side lobe level of phase-coded signal.
Brief description of the drawings
Fig. 1 is a kind of mismatched filter of low range resolution ratio phase-coded signal of code check high provided in an embodiment of the present invention
The schematic flow sheet of optimization method；
Fig. 2 be mismatched filter provided in an embodiment of the present invention phase-coded signal is carried out the main lobe after mismatch filter with
Phase-coded signal carries out the main lobe contrast schematic diagram after matched filtering, and abscissa represents relative time delay, and unit is us, ordinate
Expression amplitude, unit is dB；
Fig. 3 be mismatched filter provided in an embodiment of the present invention phase-coded signal is carried out the result after mismatch filter with
Phase-coded signal carries out the Comparative result schematic diagram after matched filtering, and abscissa represents relative time delay, and unit is us, ordinate
Expression amplitude, unit is dB；
Fig. 4 is that mismatched filter provided in an embodiment of the present invention carries out the main lobe phase after mismatch filter to phase-coded signal
Position carries out the contrast schematic diagram of the main lobe phase after matched filtering with phase-coded signal, and abscissa represents relative time delay, unit
It is us, ordinate represents phase, and unit is degree.
Specific embodiment
Reference picture 1, it is of the invention to realize that step is as follows.
Step 1, determines phase-coded signal s.
Determine that phase-coded signal s specifically includes following sub-step：
(1a) setting first phase encoded signal s1Baud Length N1。
(1b) encodes s according to first phase1Baud Length N1, determine the Baud Length N of phase-coded signal ss。
Wherein, Ns=b × N1, b is the increase multiple of the Baud Length of phase-coded signal s, and b values are integer, and Nh-
NsIt is even number, Nh+NsAlso it is even number.
Main lobe width controlled quentity controlled variable M=fix (δ × b) of (1c) setting phase-coded signal s.
Wherein, δ is empirical coefficient, and in the range of (0.5-1.0), usual value is 0.8 to general value；B represents phase
The increased multiple of transmission of symbols speed of encoded signal s, namely phase-coded signal s Baud Lengths increase multiple；fix(·)
Represent downward round numbers.
Here value formula M=fix (δ × b) for determining the main lobe width controlled quentity controlled variable M of phase-coded signal s is a warp
Test formula.In practice, the M values for being obtained according to the empirical equation there may be the bandwidth that significantly changes phase-coded signal s
Situation, now needs the value artificially to phase-coded signal s main lobe width controlled quentity controlled variables M to be adjusted, and the process of M values adjustment is abided by
Follow following rule：Increase the value of phase-coded signal s main lobe width controlled quentity controlled variables M, the bandwidth of phase-coded signal s reduces；Subtract
The value of small phase-coded signal s main lobe width controlled quentity controlled variables M, the bandwidth of phase-coded signal s increases.Therefore, it is determined that main lobe
Compromise is needed to consider that main lobe approaches the bandwidth of the degree and phase-coded signal s for expecting main lobe during the value of width control system amount M
Size.
The Baud Length N of (1d) according to phase-coded signal ss, the main lobe width controlled quentity controlled variable M of phase-coded signal s, structure
Object function and constraints on phase-coded signal s.
Object function and constraints on phase-coded signal s are built, following sub-step is specifically included：
A () calculates 1 distance side lobe ρ of kth of phase-coded signal sk1,
Wherein, k1=M+1, M+2 ..., (N are shifteds- 1), M is the main lobe width controlled quentity controlled variable of phase-coded signal s, NsRepresent
The Baud Length of phase-coded signal s, Jk1To slide matrix, matrix J is slidedk1Concrete form be：
In formula, ()TTransposition is represented, 0 represents full null matrix, and I represents unit matrix, the dimension of the subscript representing matrix of 0 and I
Number.
The distance side lobe ρ of (b) according to phase-coded signal sk1, obtain the peak sidelobe PSL of distance side lobes, PSLs
=max | ρk1|；Peak sidelobe PSL according to distance side lobes, the object function for obtaining phase-coded signal s is：
PSLs, wherein, p is the phase vectors of phase-coded signal s；Wherein, max is represented and is taken maximum, and min is represented and taken minimum value, |
| represent Modulus of accesss
C () determines the expectation main lobe b of phase-coded signal sn。
Wherein, desired main lobe can be rectangle, the main lobe of sinc functions, or do not increase the phase before he number
The main lobe of encoded signal.Main lobe b is expected described in this examplenThe main lobe for being shaped as the first phase encoded signal, by institute
State first phase encoded signal s1The interval equably discrete column vector for turning to 2 × M+1 dimensions of the corresponding functional value in main lobe part,
M-th element therein is taken to the 1st element as expectation main lobe bnValue.
D () is according to expectation main lobe bn, the constraints for obtaining phase-coded signal s is：
0≤p (k)≤2 π, k=1,2 ..., Ns,
Wherein, s.t. represents constraints, | | represent Modulus of access, bn(k2) represent and expect main lobe bn2 elements of kth,
Jk2To slide matrix, k2=1 is shifted, 2 ..., M, M are the main lobe width controlled quentity controlled variable of phase-coded signal s, NsRepresent phase code
The Baud Length of signal s, p (k) represents k-th element of the phase vectors p of phase-coded signal s, k=1,2 ..., Ns, γ is
Constant, γ represents that the main lobe of phase-coded signal s approaches expectation main lobe bnDegree, value be 0.01~0.2.
In theory, 0≤p of phase restriction (k)≤2 π, k=1,2 ..., NsCan ignore, because exp function cycles are 2 π,
But in actual solution procedure, phase restriction ensure that optimized algorithm is preferably restrained, therefore, be given in constraints
Phase restriction.
Constant γ is an empirical value, is found in emulation, when the shape for being shaped as sinc function main lobes of main lobe is expected,
γ values 0.1 just can make the main lobe shape after phase-coded signal s pulse compressions approach the shape for expecting main lobe well.Such as
Fruit reduces the value of constant γ, and the main lobe after phase-coded signal s pulse compressions can preferably approach the shape for expecting main lobe,
But the peak sidelobe of distance side lobe can be raised.If the value of increase constant γ, the peak sidelobe of distance side lobe
Can reduce, but the main lobe after phase-coded signal s pulse compressions is deteriorated to expecting the Approximation effect of main lobe.Therefore, it is determined that
Compromise is needed to consider during the value of constant γ.
(1e) solves the phase-coded signal s according to the object function and the constraints.
Under constraints, p norm optimization Algorithm for Solving object functions, the phase-coded signal s after being optimized are used.
Step 2, the length N of setting mismatched filter hh, and according to phase-coded signal s and the length of mismatched filter h
Nh, build the object function on mismatched filter h.
Mismatched filter h is the mismatched filter of phase-coded signal s.
The length N of setting mismatched filter hh, and according to the Baud Length N of phase-coded signal ssWith main lobe width control
Amount M, and mismatched filter h length Nh, the object function for solving mismatched filter h is built, specifically include following sub-step：
(2a) phase-coded signal s is ρ by the pulse compression vector after mismatched filter hc,
Wherein,Represent convolution, []TRepresent transposition, i ∈ [1, Ns+Nh- 1], NsRepresent that the code element of phase-coded signal s is long
Degree, NhRepresent the length of mismatched filter h.
(2b) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, to each distance side lobe
Modulus of access, obtains distance side lobe modulus value vector ρ respectively：
Wherein, k ∈ [1, (Ns+Nh)/2-M-1], l ∈ [(Ns+Nh)/2+M+1, Ns+Nh-1]；
In phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, extract the distance on the main lobe left side
The left comer code vector L of secondary lobe is L=[1 ..., m ..., (Ns+Nh)/2-M-1] and main lobe on the right of distance side lobe right corner code
Vectorial R is R=[(Ns+Nh)/2+M+1 ..., n ..., Ns+Nh- 1], wherein, m ∈ [1, (Ns+Nh)/2-M-1], n ∈ [(Ns+
Nh)/2+M+1, Ns+Nh, and m and n are positive integer -1]；After remembering the distance side lobe Modulus of access of output, secondary lobe modulus value vector is obtained：
Wherein, | | represent Modulus of access, k ∈ [1, (Ns+Nh)/2-M-1], l ∈ [(Ns+Nh)/2+M+1, Ns+Nh-1]。
(2c) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, the actual main lobe of output
BmFor
Wherein, t ∈ [(Ns+Nh)/2-M, (Ns+Nh)/2+M]。
(2d) determines expectation main lobe bs of the phase-coded signal s by output after mismatched filter hm, wherein, desired master
Valve can be rectangle, the main lobe of sinc functions, or do not increase the main lobe of the phase-coded signal before he number.This example
Middle expectation main lobe bmBe shaped as first phase encoded signal s1Main lobe, by first phase encoded signal s1Main lobe part pair
The interval equably discrete column vector for turning to 2 × M+1 dimensions of the functional value answered, and using column vector as expecting main lobe bmValue.
(2e) is according to distance side lobe modulus value vector ρ, actual main lobe BmWith expectation main lobe bm, build object function
Wherein, h is that length is NhMismatched filter, | | | |∞Infinite Norm is represented, α is the arithmetic number being manually set,
For Sidelobe Suppression effect and the main lobe approximation ratio of compromising.If the value of α is larger, phase-coded signal s is by mismatch filter
Main lobe after device h can preferably approach the shape of expectation main lobe, but can raise the peak sidelobe of distance side lobe, if
The value of α is smaller, and the peak sidelobe of distance side lobe can be reduced, but phase-coded signal s is by after mismatched filter h
Main lobe to expect main lobe Approximation effect be deteriorated, therefore, it is determined that arithmetic number α value when need compromise consider.
Step 3, builds the bound for objective function on mismatched filter h.
The symmetrical zero padding of phase-coded signal s head and the tail is extended to length for N by (3a)hPhase-coded signal
(3b) is N by lengthhPhase-coded signalSpin upside down and take conjugation and obtain matched filter
Wherein, ()*Represent the conjugation of ().
(3c) is according to mismatched filter h, matched filterObtain the constraint of the object function on mismatched filter h
Condition is：
Wherein, s.t. represents constraints, | | | |22 norms are represented, ε is that length is NhPhase-coded signalMistake
Distribution coefficient, if the value of mismatching ε is larger, phase-coded signal s is by the peak sidelobe meeting after mismatched filter h
Reduce, but snr loss can raise, and if the value of ε is smaller, snr loss can diminish, but phase-coded signal s is passed through
The peak sidelobe crossed after mismatched filter h can be raised；Snr loss that can be according to actual needs is adjusted, generally
Value is 0 ＜ ε ＜ 10.
Step 4, according to bound for objective function, solves the object function of mismatched filter h：
The value of energy weight coefficient α belongs in the range of [0,1].
Step 5, determines mismatched filter h.
In constraintsUnder, solve object function using the convex optimization tool bag cvx of MatlabObtain mismatched filter h.
Effect of the invention can be further illustrated with reference to emulation experiment.
(1) simulation parameter
First phase encoded signal s1Baud Length N1=64, transmission of symbols speed increases multiple b=4, then phase code
The Baud Length N of signal ss=256；The length N of mismatched filterhThe main lobe width controlled quentity controlled variable M of=512, phase-coded signal s
=4, that expects main lobe is shaped as Baud Length NsMain lobe shape after=256 phase-coded signal s pulse compressions, weight coefficient
α=0.3, mismatching ε=2.25.
(2) emulation content
Emulation 1, designed phase encoded signal s.
Baud Length N according to phase-coded signal ssWith main lobe width controlled quentity controlled variable M, by minimizing the phase code
The peak sidelobe of signal s and approach and expect that main lobe is shaped as criterion, using p norm optimization algorithms obtain low distance side lobe
Phase-coded signal s.
Emulation 2, matched filtering and mismatch filter are carried out to phase-coded signal s respectively.
In simulations, setting Baud Length N1=64 phase-coded signal s1With Baud Length Ns=256 phase code
Signal s has identical time width 15.36us, and sample frequency is 100MHz.
A () phase-coded signal s carries out the emulation of the main lobe characteristic of output signal after matched filtering and mismatch filter respectively
With analysis.
Matched filtering and mismatch filter are carried out respectively to phase-coded signal s, by the main lobe result Modulus of access of output signal
After be drawn as X-Y scheme, as shown in Fig. 2 its center line-o- represents the main lobe characteristic of output signal after phase-coded signal s matched filterings
Schematic diagram, main lobe characteristic schematic diagrames of the line-*-expression phase-coded signal s by output signal after mismatched filter h.
As shown in Figure 2, after phase-coded signal s matched filterings output signal main lobe, i.e. ,-o- the lines in figure, with position compile
The main lobe of output signal, i.e.-* in figure-line, essentially coincide after code signal s mismatch filters；Range value is at -3dB in the figure,
The width of the main lobe of output signal is equal after phase-coded signal s matched filterings and mismatch filter.Therefore, the inventive method design
The mismatched filter for obtaining is constant by the bandwidth of output signal after mismatch filter to phase-coded signal s.
B () phase-coded signal s carries out the emulation of the amplitude characteristic of output signal after matched filtering and mismatch filter respectively
With analysis.
Matched filtering and mismatch filter are carried out respectively to phase-coded signal s, will be drawn after the result Modulus of access of output signal
Into X-Y scheme, as shown in figure 3, after where the dotted line signifies that phase-coded signal s matched filterings output signal amplitude characteristic figure, it is real
Line represents amplitude characteristic figures of the phase-coded signal s by output signal after mismatched filter h.
From the figure 3, it may be seen that the maximum amplitude value of output signal is -36.98dB, phase after phase-coded signal s matched filterings
Encoded signal s by output signal after mismatched filter h maximum amplitude value for the distance after -42.92dB, with matched filtering by
Valve is compared, reduce 5.94dB, snr loss -0.085dB；Therefore, the mismatched filter h that present invention design is obtained is being protected
While card phase-coded signal s bandwidth is constant, by increasing the Baud Length of phase-coded signal s, phase is reduce further
The distance side lobe level of position encoded signal s.
C () phase-coded signal s carries out the main lobe phase characteristic of output signal after matched filtering and mismatch filter respectively
Simulation and analysis.
Take out phase-coded signal s carries out the phase of the main lobe part of output signal after matched filtering and mismatch filter respectively
Position, changes into angle value, is drawn as X-Y scheme, as shown in figure 4, being exported after where the dotted line signifies that phase-coded signal s matched filterings
The phase characteristic figure of the main lobe of signal, solid line represents phase-coded signal s by the main lobe of output signal after mismatched filter h
Phase characteristic figure.
As shown in Figure 4, after phase-coded signal s matched filterings output signal main lobe phase excursion for [-
0.164 °, 0.164 °], solid line represents changes of the phase-coded signal s by the main lobe phase of output signal after mismatched filter h
Scope is [- 7.20710-9°, -7.20710-9°], and phase-coded signal s is by output signal after mismatched filter h
Main lobe phase compared with the main lobe phase of output signal after phase-coded signal s matched filterings, maximum difference can ignore not
Meter.Therefore, phase-coded signal s is defeated after the mismatched filter h that the inventive method design is obtained ensure that pulse compression
The main lobe phase for going out signal is basically unchanged.
Claims (4)
1. the mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high, it is characterised in that including following
Step：
Step 1, determines phase-coded signal s；
Step 2, the length N of setting mismatched filter hh, and according to the phase-coded signal s and the length of the mismatched filter h
Degree Nh, build the object function on the mismatched filter h；
Wherein, step 2 specifically includes following sub-step：
(2a) phase-coded signal s is ρ by the pulse compression vector after mismatched filter hc,i∈[1,Ns+Nh- 1], NsIt is the Baud Length of phase-coded signal s；
(2b) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, each distance side lobe is distinguished
Modulus of access, obtains distance side lobe modulus value vector ρ：
Wherein, k ∈ [1, (Ns+Nh)/2-M-1], l ∈ [(Ns+Nh)/2+M+1,Ns+Nh-1]；
(2c) is in phase-coded signal s by the pulse compression vector ρ after mismatched filter hcIn, the actual main lobe B of outputmFor
Wherein, t ∈ [(Ns+Nh)/2-M,(Ns+Nh)/2+M]；
(2d) determines expectation main lobe bs of the phase-coded signal s by output after mismatched filter hm, wherein, the expectation main lobe bm
The main lobe for being shaped as the phase-coded signal s, the corresponding functional value in the main lobe part of the phase-coded signal s is interval
The equably discrete column vector for turning to 2 × M+1 dimensions, and using the column vector as expecting main lobe bmValue；
(2e) is according to the distance side lobe modulus value vector ρ, the actual main lobe BmWith the expectation main lobe bm, build on described
The object function of mismatched filter hThe value of energy weight coefficient α belongs to the scope of [0,1]
It is interior；
Step 3, builds the bound for objective function on the mismatched filter h；
Step 4, according to the bound for objective function, solves the object function of the mismatched filter h；
Step 5, determines the mismatched filter h.
2. the mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high according to claim 1,
Characterized in that, step 1 specifically includes following sub-step：
(1a) setting first phase encoded signal s1Baud Length N1；
(1b) encodes s according to the first phase1Baud Length N1, determine the Baud Length N of the phase-coded signal ss,
Wherein, Ns=b × N1, b is the increase multiple of the Baud Length of phase-coded signal s, and b values are integer；
(1c) sets main lobe width controlled quentity controlled variable M=fix (δ × b) of the phase-coded signal s, and δ values are in (0.5~1.0)
In the range of, b is the increase multiple of the Baud Length of phase-coded signal s, and b values are integer；
The Baud Length N of (1d) according to the phase-coded signal ss, the main lobe width controlled quentity controlled variable M of the phase-coded signal s,
Build object function and constraints on the phase-coded signal s；
(1e) solves the phase-coded signal s according to the object function and constraints of the phase-coded signal s.
3. the mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high according to claim 1,
Characterized in that, step 3 specifically includes following sub-step：
The symmetrical zero padding of phase-coded signal s head and the tail is extended to length for N by (3a)hPhase-coded signal
(3b) is N by lengthhPhase-coded signalSpin upside down and take conjugation and obtain matched filter
Wherein, ()*Represent the conjugation of ()；
(3c) is according to the mismatched filter h, the matched filterObtain the object function on the mismatched filter h
Constraints beWherein, ε is that length is NhPhase-coded signalMismatching.
4. the mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high according to claim 2,
Characterized in that, sub-step (1d) also specifically includes following sub-step：
A () calculates 1 distance side lobe ρ of kth of phase-coded signal sk1：
Wherein, k1=M+1, M+2 ..., (N are shifteds- 1), M is the main lobe width controlled quentity controlled variable of phase-coded signal s, Jk1To slide
Matrix, slides matrix Jk1Form be：
In formula, ()TTransposition is represented, 0 represents full null matrix, and I represents unit matrix, the dimension of the subscript representing matrix of 0 and I；
The distance side lobe ρ of (b) according to phase-coded signal sk1, obtain the peak sidelobe PSL of distance side lobes, PSLs=max |
ρk1|；Peak sidelobe PSL according to distance side lobes, the object function for obtaining phase-coded signal s is：Its
In, p is the phase vectors of phase-coded signal s；
C () determines the expectation main lobe b of phase-coded signal sn, wherein, the expectation main lobe bnBe shaped as sinc function main lobes
Shape, the interval equably discrete column vector for turning to 2 × M+1 dimensions of the corresponding functional value in sinc function main lobes part takes wherein
M-th element to the 1st element as expecting main lobe bnValue；
D () is according to the expectation main lobe bn, the constraints for obtaining phase-coded signal s is：
0≤p (k)≤2 π, k=1,2 ..., Ns,
Wherein, bn(k2) represent and expect main lobe bn2 elements of kth, Jk2To slide matrix, k2=1 is shifted, 2 ..., M, M are phase
The main lobe width controlled quentity controlled variable of position encoded signal s, p (k) represents k-th element of phase vectors p, k=1,2 ..., Ns, γ is normal
Number, γ represents that the main lobe of phase-coded signal s approaches expectation main lobe bnDegree, value be 0.01~0.2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510305913.9A CN105044681B (en) | 2015-06-05 | 2015-06-05 | The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510305913.9A CN105044681B (en) | 2015-06-05 | 2015-06-05 | The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105044681A CN105044681A (en) | 2015-11-11 |
| CN105044681B true CN105044681B (en) | 2017-06-20 |
Family
ID=54451358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510305913.9A Active CN105044681B (en) | 2015-06-05 | 2015-06-05 | The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105044681B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106093876B (en) * | 2016-07-19 | 2018-11-09 | 西安电子科技大学 | The orthogonal wide main lobe phase encoded signal design method of distributed MIMO radar |
| CN107942327A (en) * | 2017-10-20 | 2018-04-20 | 西安电子科技大学 | Single channel HRWS SAR imaging methods based on impulse phase coding |
| CN110471042B (en) * | 2019-08-05 | 2021-05-18 | 西安电子科技大学 | Design method of wide-main-lobe multi-pulse phase coding signal |
| CN113489573B (en) * | 2021-09-07 | 2021-11-12 | 中国人民解放军国防科技大学 | Bistatic radar communication integrated waveform design method, computer device and medium |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5786788A (en) * | 1996-10-08 | 1998-07-28 | Raytheon Company | Radar system and method for reducing range sidelobes |
| CN103344944B (en) * | 2013-07-02 | 2015-05-13 | 西安电子科技大学 | Radar pulse compression filter optimization design method applied to random signal waveforms |
| CN103969633B (en) * | 2014-03-28 | 2016-05-18 | 西安电子科技大学 | In clutter, detect the grading design method of target MIMO radar emission waveform |
| CN104198996B (en) * | 2014-08-30 | 2016-08-24 | 西安电子科技大学 | The low sidelobe phase-coded signal method for designing of the low range resolution ratio of high code check |
-
2015
- 2015-06-05 CN CN201510305913.9A patent/CN105044681B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105044681A (en) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105044681B (en) | The mismatched filter optimization method of the low range resolution ratio phase-coded signal of code check high | |
| CN104198996B (en) | The low sidelobe phase-coded signal method for designing of the low range resolution ratio of high code check | |
| CN105044682B (en) | Wide-main-lobe phase coding signal and mismatch filter combined optimization method | |
| CN103344944B (en) | Radar pulse compression filter optimization design method applied to random signal waveforms | |
| CN106772457B (en) | A kind of satellite navigation anti-interference method based on empty time-frequency framework | |
| CN104865558B (en) | The combined optimization method of phase-coded signal and mismatched filter based on p norms | |
| CN105974376B (en) | A kind of SAR radio frequency interferences suppressing method | |
| CN103197300B (en) | Real-time processing method for cancellation of direct wave and clutter of external radiation source radar based on graphic processing unit (GPU) | |
| CN102608588A (en) | Broadband sub-matrix adaptive beamforming method based on sub-band decomposition | |
| CN103885045B (en) | Based on the circulation associating Adaptive beamformer method of Subarray partition | |
| CN102355273B (en) | Digital channelizing method and circuit | |
| CN103439691B (en) | Method for broadband networking radar to restrict narrow-band interference | |
| CN106093877B (en) | Orthogonal width main lobe phase encoded signal and mismatched filter combined optimization method | |
| CN106526569A (en) | Broadband MIMO radar sparse spectrum waveform design method based on alternative iteration | |
| CN108874743A (en) | A kind of adding window method controlling frequency spectrum principal subsidiary lobe characteristic | |
| CN105337587A (en) | Non maximum extraction system integrated filter bank construction method based on DFT | |
| CN106093876A (en) | The orthogonal wide main lobe phase coding Design of Signal method of distributed MIMO radar | |
| CN104242980B (en) | The design of a kind of Sub-1G RF front-end circuit based on RF energy measuring and parameter regulation means | |
| CN105301572A (en) | Multi-peak low sidelobe phase encode signal and mismatched filter joint design method | |
| CN104536013A (en) | Weight calculation method for nulling antenna of satellite navigation receiver | |
| CN112034434A (en) | Radar radiation source identification method based on sparse time-frequency detection convolutional neural network | |
| CN105044685A (en) | Ultra low sidelobe pulse transformer design method of linear frequency modulation pulse signal | |
| CN102710284A (en) | Two-dimensional self-adaptive filtering narrow-band interference suppression device | |
| Zhao et al. | Block NLMS cancellation algorithm and its real-time implementation for passive radar | |
| CN103630880B (en) | For the waveform diversity method of array SAR three-dimensional imaging |
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 |