CN103197284B - Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance - Google Patents
Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance Download PDFInfo
- Publication number
- CN103197284B CN103197284B CN201310092762.4A CN201310092762A CN103197284B CN 103197284 B CN103197284 B CN 103197284B CN 201310092762 A CN201310092762 A CN 201310092762A CN 103197284 B CN103197284 B CN 103197284B
- Authority
- CN
- China
- Prior art keywords
- phi
- burst
- correlation
- doppler tolerance
- cost function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a radar waveform design method restrained by autocorrelation, orthogonality and Doppler tolerance. Firstly, a cost function with stable properties and with overall consideration of the autocorrelation, the cross correlation and the Doppler tolerance of wave forms is designed. Secondly, by using a Greedy random search optimization algorithm, repeated iterations are carried out until a system no longer receives any phase change, and the wave form design meeting the requirements of the full polarization measurement of radar is obtained.
Description
Technical field
The present invention relates to radar waveform design field, particularly relate to the radar waveform method for designing under a kind of auto-correlation, the constraint of orthogonal and doppler tolerance.
Background technology
In Electromagnetic Wave Propagation process, on a certain point of fixity in space, the time dependent mode of spatial orientation of electric field intensity is called polarization.Target is equivalent to a Polarization changer, and its polarization information can characterize the full detail of its scattering properties, and therefore the utilization of polarization information can improve the target detection of radar, tracking and recognition performance.Polarization technique is adopted to also help ground/ocean clutter cancellation.The prerequisite that polarization information utilizes is that radar possesses polarization measurement ability.Instantaneous Polarimetry requires that horizontal polarization and vertical polarization transmitted waveform have good auto-correlation and cross-correlation (orthogonal) performance.Autocorrelation performance is undesirable will cause pulse compression after distance side lobe high, be unfavorable for detect Small object; The undesirable receiving cable isolation that will cause of orthogonal performance reduces, and produces cross jamming.
The form of orthogonal waveforms has a lot, comprises frequency orthogonal, phase encoding orthogonal etc.The PARSAX radar of Delft Polytechnics of Holland development adopts the Continuous Wave with frequency modulation system of bistatic, has instantaneous Polarimetry ability.Horizontal polarization (H) transmission channel adopts the LFM waveform of positive frequency modulation slope, and vertical polarization (V) transmission channel adopts the LFM waveform of negative frequency modulation slope.Because positive frequency modulation slope LFM signal and negative frequency modulation slope LFM signal are just accurate orthogonal, go tiltedly process to create cross aisle interference, be presented as linear FM signal item.Timed automata is larger, and cross aisle interference is more serious.The Hai Deng of University of New Orleans in 2004, using signal autocorrelation side-lobe energy and cross-correlation energy as cost function, adopts simulated annealing to be optimized combination to frequency coding sequence.The people such as the Wang Dunyong of Airforce Radar institute improve the cost function of waveform on the basis that Hai Deng works, and add autocorrelation sidelobe peak value and cross-correlation peak value, and genetic algorithm is carried out wave sequence optimizing as optimizing algorithm.
But adopt during phase modulation waveform and face Doppler sensitivity problem, for high-speed target, the signal to noise ratio (S/N ratio) that the Doppler shift that echo exists makes pulse pressure export sharply reduces.Therefore there is a need to when design object function increase doppler tolerance constraint condition.
Desirable transmitted waveform is that respective matched filtering output will have enough low secondary lobe (autocorrelation performance), pairwise orthogonal is also wanted to require (their cross correlation) to be applicable to instantaneous Polarimetry between waveform, simultaneously each waveform should have enough large doppler tolerance (matched filter as pulse compression can not be too responsive to frequency displacement) as far as possible, but does not also have a kind of waveform design method to consider the three aspect constraint of autocorrelation, cross correlation and doppler tolerance before simultaneously.
Summary of the invention
The invention provides a kind of waveform design method of the autocorrelation, cross correlation and the doppler tolerance performance that consider waveform, the method is by the weight of 3 elements above-mentioned in Readjusting cost function and use Greedy Stochastic search optimization algorithm, designs the various phase modulation waveforms satisfied the demands.
The waveform design method taking into account waveform autocorrelation, cross correlation and doppler tolerance performance of the present invention comprises: 1) first design a kind of cost function, if total L burst in orthogonal signal code character, each burst code length N, encoding phase value number is M, each chip in each burst is subcode, then l burst shows by formula (1);
{s
l(n)=exp[jφ
l(n)],n=1,2,---,N},l=1,2,---,L (1)
Encoding phase span is
The phasing matrix of whole orthogonal signal code character can be expressed as
Cost function is:
E=λ
ac× auto-correlation cost+λ
cc× cross-correlation cost+λ
dt× doppler tolerance cost (10)
Wherein λ
acrepresent the weighting coefficient of auto-correlation constraint, λ
ccrepresent the weighting coefficient of interrelational constraint, λ
dtbe the weighting coefficient of doppler tolerance constraint, doppler tolerance cost is expressed as E
dt:
Wherein B is signal bandwidth, f
dtfor the doppler tolerance of signal, what target setting speed was more than or equal to 200m/s is high-speed target, and what target velocity was less than 200m/s is slower-velocity target;
2) according to the requirement of radar Polarimetry, in cost function, auto-correlation, the weighting coefficient that orthogonal and doppler tolerance three performances are corresponding is set;
3) with the formula (10) are cost function, optimized algorithm is utilized to be optimized phasing matrix in encoding phase span, finally obtain the phasing matrix value that least cost function value is corresponding, and then be met the optimum waveform of radar Polarimetry requirement.
Beneficial effect of the present invention:
1. have devised and consider waveform autocorrelation, cross correlation and doppler tolerance three cost function of stable performance of requiring of aspect
In Waveform Design, this 3 individual character retrains by self-corresponding weighting coefficient each in cost function.The relative value changing each weighting coefficient in cost function can adjust the performance of waveform, and weighting coefficient is larger, shows more to stress corresponding performance.
2. utilize Greedy Stochastic search optimization algorithm, the orthogonal phase modulation wave sequence meeting radar Polarimetry demand of phase and random length arbitrarily can be designed expeditiously: after initialization is carried out to phasing matrix, exchange phase place disturbance in effective valued space of each subcode of facies-suite, the knots modification Δ E of cost function before and after calculation perturbation.Only accept disturbance as Δ E < 0.Through successive ignition, until obtain the Waveform Design meeting radar Polarimetry demand when system no longer accepts any phase change.
Accompanying drawing explanation
Fig. 1 is the Waveform Design process flow diagram with Greedy algorithm;
Embodiment
In order to understand technical scheme of the present invention better, below in conjunction with drawings and the specific embodiments, the present invention is described in further detail.
The invention provides a kind of radar waveform method for designing, flow process as shown in Figure 1, basic ideas design to consider waveform autocorrelation, cross correlation and doppler tolerance three cost function of stable performance of requiring of aspect, this 3 individual character retrains by self-corresponding weighting coefficient each in cost function, and the performance of waveform controls by the relative value changing each weighting coefficient in cost function; Then Greedy Stochastic search optimization algorithm design is utilized to go out the orthogonal phase modulation wave sequence meeting radar Polarimetry demand of phase and random length arbitrarily.Specific as follows:
First, cost function is determined.
If total L burst in orthogonal signal code character, each sequence code length N, encoding phase value number is M, and each chip in each burst is subcode.Then l burst shows by formula (1);
{s
l(n)=exp[jφ
l(n)],n=1,2,---,N},=1,2,---,L (1)
Encoding phase span is
The phasing matrix of whole orthogonal signal code character can be expressed as
The autocorrelation function of orthogonal signal and cross correlation function should meet following two conditions:
With
Wherein A (s
l, k) represent in l burst be spaced apart k subcode autocorrelation function, C (s
p, s
q, k) represent p burst and q burst be spaced apart k subcode cross correlation function;
Associating code character signal matrix can obtain
Utilize conventional method, auto-correlation cost and cross-correlation cost can be expressed as
with
Doppler tolerance is defined as frequency separation during ambiguity function peak value decline 3dB or 6dB, is obtained by the mode of Greedy random search.The doppler tolerance of note signal is f
dt, the doppler tolerance key element in cost function is expressed with following formula:
Wherein λ
dtbe the weighting coefficient of doppler tolerance constraint, B is signal bandwidth, and what target setting speed was more than or equal to 200m/s is high-speed target, and what target velocity was less than 200m/s is slower-velocity target.For high-speed target detection, when doppler tolerance is larger, cost is less.For slower-velocity target, owing to not needing to pursue large doppler tolerance, and doppler tolerance is less means that velocity resolution is higher, and be conducive to distinguishing low-speed motion target and clutter, more hour cost is less therefore to make doppler tolerance.
Consider the requirement of autocorrelation, cross correlation, doppler tolerance three aspect, cost function is designed to
Wherein λ
acrepresent the weighting coefficient of auto-correlation constraint, λ
ccrepresent the weighting coefficient of interrelational constraint, λ
dtrepresent the weighting coefficient of doppler tolerance.Change the performance that each weighting coefficient relative value can adjust waveform.Weighting coefficient is larger, shows more to stress corresponding performance.
Secondly, adopt Greedy algorithm to be optimized phase sequence, algorithm flow as shown in Figure 1.
As previously mentioned, if total L burst in orthogonal signal code character, each sequence code length N, encoding phase number is M, and whole block signal matrix can be expressed as:
1) according to the requirement of radar Polarimetry, in cost function, set auto-correlation, the weighting coefficient that orthogonal and doppler tolerance three performances are corresponding, and in encoding phase span, random initializtion is carried out to block signal matrix;
2) for the phase of the 1st subcode of the 1st burst
1(1), replace by the residue encoding phase value except intialization phase value in encoding phase span successively, and calculate the cost function value after each replacement, by φ corresponding for least cost function value
1(1) phase value is fixed up as a result;
3) all the other N-1 subcodes of the 1st burst are repeated to the operation of step 2;
4) all the other L-1 bursts repeat and the 1st operation that sequence is identical, obtain the phasing matrix value that least cost function value in whole burst is corresponding, namely meet the optimum waveform of radar Polarimetry requirement.
This optimization method is specially: after carrying out initialization to phasing matrix, to phase place disturbance in effective valued space of each subcode of each burst, the knots modification Δ E of cost function before and after calculation perturbation, only accept disturbance as Δ E < 0, until the phase place of each subcode of each burst no longer disturbance be namely met the optimum waveform that radar Polarimetry requires.
Greedy algorithm operation efficiency is very high.Owing to only accepting the phase perturbation that cost function reduces, be therefore likely absorbed in local minimum, but the waveform of requirement can be met by Multi simulation running.
The improvement optimizing rear waveform performance is illustrated below in conjunction with simulation example:
Greedy algorithm is adopted to be optimized Frank code phase sequences.The Frank code of 44 phase 16 subcodes is used, bandwidth B=1MHz, sample frequency f in emulation
s=20MHz, doppler tolerance is defined as the frequency separation of ambiguity function peak value decline 3dB.Before optimizing, the performance of each sequence of Frank code is in table 1.For the weighting coefficient λ of the doppler tolerance factor in cost function
acbe respectively the situation of 1 and 10, the waveform performance after optimization is in table 2,3;
Original auto-correlation PSL(dB) | -17.0927 | -17.0927 | -21.0721 | -17.0927 | ---- | ---- |
Original cross-correlation peak value (dB) | -23.5218 | 23.5218 | -18.0618 | -22.9226 | -19.0849 | -19.0849 |
Original doppler tolerance (kHz) | 268.38 | 267.71 | 267.71 | 268.38 | ---- | ---- |
The waveform performance of the original Frank code of table 1
Auto-correlation PSL(dB after optimizing) | -17.0927 | -17.0927 | -21.0721 | -17.0927 | ---- | ---- |
Cross-correlation peak value (dB) after optimizing | -23.5218 | -23.5218 | -18.0618 | -22.9226 | -19.0849 | -19.0849 |
Doppler tolerance (kHz) after optimizing | 283.69 | 267.71 | 267.71 | 268.38 | ---- | ---- |
Waveform performance (the λ of the Frank code after table 2 optimization
ac=1 λ
cc=1 λ
dt=1)
Auto-correlation PSL(dB after optimizing) | -12.9430 | -12.9430 | -10.1030 | -14.0824 | ---- | ---- |
Cross-correlation peak value (dB) after optimizing | -20.8279 | -21.8469 | -18.0618 | -21.1394 | -18.1291 | -19.0849 |
Doppler tolerance (kHz) after optimizing | 297.79 | 297.79 | 290.34 | 278.73 | ---- | ---- |
Waveform performance (the λ of the Frank code after table 3 optimization
ac=1 λ
cc=1 λ
dt=10)
Waveform performance before and after being optimized by contrast, can find that the weight increasing the doppler tolerance factor can improve the doppler tolerance optimizing waveform.
To sum up, waveform design method provided by the invention adopts Greedy random search algorithm to carry out waveform optimization design, by the weighting of element in cost function, the autocorrelation of waveform, cross correlation and doppler tolerance performance can be taken into account, and the weight of three can be regulated arbitrarily, thus the demand (orthogonal performance) of channel separation can be met, the demand that weak target measures (the low sidelobe performance of autocorrelation function) and high-speed target measurement (doppler tolerance performance) can be taken into account again.
In sum, these are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (2)
1. the radar waveform method for designing under auto-correlation, the constraint of orthogonal and doppler tolerance, is characterized in that,
1) a kind of cost function is first designed, if total L burst in orthogonal signal code character, each burst code length N, encoding phase value number is M, each chip in each burst is subcode, then l burst shows by formula (1);
{s
i(n)=exp[jφ
i(n)],n=1,2,---,N),l=1,2,---,L (1)
Encoding phase span is
The phasing matrix of whole orthogonal signal code character can be expressed as
The autocorrelation function of orthogonal signal and cross correlation function should meet following two conditions:
With
Wherein A (s
l, k) represent in l burst be spaced apart k subcode autocorrelation function, C (s
p, s
q, k) represent p burst and q burst be spaced apart k subcode cross correlation function;
Cost function is:
E=λ
ac× auto-correlation cost+λ
cc× cross-correlation cost+λ
dt× doppler tolerance cost (6)
Wherein λ
acrepresent the weighting coefficient of auto-correlation constraint, λ
ccrepresent the weighting coefficient of interrelational constraint, λ
dtbe the weighting coefficient of doppler tolerance constraint, doppler tolerance cost is expressed as E
dt:
(7)
Wherein B is signal bandwidth, f
dtfor the doppler tolerance of signal, what target setting speed was more than or equal to 200m/s is high-speed target, and what target velocity was less than 200m/s is slower-velocity target;
2) according to the requirement of radar Polarimetry, in cost function, auto-correlation, the weighting coefficient that orthogonal and doppler tolerance three performances are corresponding is set;
3) with the formula (6) are cost function, optimized algorithm is utilized to be optimized phasing matrix in encoding phase span, finally obtain the phasing matrix value that least cost function value is corresponding, and then be met the optimum waveform of radar Polarimetry requirement.
2. the radar waveform method for designing that a kind of auto-correlation as claimed in claim 1, the constraint of orthogonal and doppler tolerance are lower, it is characterized in that, employing Greedy optimized algorithm is optimized, and is specially:
One, in encoding phase span, initialization is carried out to each matrix element in phasing matrix;
Two, for the phase of the 1st subcode of the 1st burst
1(1), replace by the residue encoding phase value except intialization phase value in encoding phase span successively, and calculate the cost function value after each replacement, by φ corresponding for least cost function value
1(1) phase value is fixed up as a result;
Three, all the other N-1 subcodes of the 1st burst are repeated to the operation of step 2;
Four, all the other L-1 bursts repeat and the 1st operation that sequence is identical, obtain the phasing matrix value that whole burst least cost function value is corresponding, namely meet the optimum waveform of radar Polarimetry requirement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310092762.4A CN103197284B (en) | 2013-03-21 | 2013-03-21 | Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310092762.4A CN103197284B (en) | 2013-03-21 | 2013-03-21 | Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103197284A CN103197284A (en) | 2013-07-10 |
CN103197284B true CN103197284B (en) | 2015-02-25 |
Family
ID=48719961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310092762.4A Expired - Fee Related CN103197284B (en) | 2013-03-21 | 2013-03-21 | Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103197284B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698751B (en) * | 2014-01-11 | 2016-02-10 | 西安电子科技大学 | The multistation radar waveform method for designing of orthogonal low Doppler sidelobne |
CN106019237B (en) * | 2016-06-23 | 2021-07-30 | 哈尔滨工业大学(威海) | Radar LFM composite waveform design method |
CN106291475A (en) * | 2016-07-20 | 2017-01-04 | 南京御达电信息技术有限公司 | A kind of for radar waveform design optimization method |
CN106443595B (en) * | 2016-09-05 | 2019-03-05 | 电子科技大学 | A kind of cognition radar waveform design method of anti-instantaneous forwarding slice reconstruct interference |
CN109204390B (en) * | 2018-09-29 | 2021-03-12 | 交控科技股份有限公司 | Train control method based on deep learning |
CN112162243B (en) * | 2020-08-28 | 2024-04-30 | 西安电子科技大学 | Method for generating MIMO radar orthogonal waveform based on DDPG model |
CN112068122B (en) * | 2020-09-09 | 2023-02-17 | 中国航空工业集团公司雷华电子技术研究所 | Doppler compensation method of phase coding signal based on tracking information |
CN112965034B (en) * | 2021-04-13 | 2023-11-17 | 哈尔滨工业大学 | Method for improving Doppler tolerance of slow time phase coded signal of sky-wave radar |
CN113567935B (en) * | 2021-07-20 | 2023-09-19 | 电子科技大学 | High-speed target detection method based on multi-sub-pulse processing |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540187A (en) * | 2010-12-13 | 2012-07-04 | 电子科技大学 | Orthogonal waveform designing method for formation flying satellites SAR (synthetic aperture radar) |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012111141A1 (en) * | 2011-02-18 | 2012-08-23 | 三菱電機株式会社 | Passive radar device |
-
2013
- 2013-03-21 CN CN201310092762.4A patent/CN103197284B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540187A (en) * | 2010-12-13 | 2012-07-04 | 电子科技大学 | Orthogonal waveform designing method for formation flying satellites SAR (synthetic aperture radar) |
Non-Patent Citations (3)
Title |
---|
Optimization of Orthogonal Discrete Frequency-Coding Waveform Based on Modified Genetic Algorithm for MIMO Radar;Bo Liu等;《Communications, Circuits and Systems, 2007. ICCCAS 2007. International Conference on》;20071231;966-970 * |
Xiongjun Fu等.Channel Isolation and Range Sidelobe Suppression for Instantaneous Full-polarization Radar.《International Conference on Information Science and Technology,March 26-28, 2011 Nanjing, Jiangsu, China》.2011, * |
多输入多输出合成孔径雷达关键技术研究;王力宝;《万方学位论文》;20110920;1-158 * |
Also Published As
Publication number | Publication date |
---|---|
CN103197284A (en) | 2013-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103197284B (en) | Radar wave form design method restrained by autocorrelation, orthogonality and Doppler tolerance | |
CN102540187B (en) | Orthogonal waveform designing method for formation flying satellites SAR (synthetic aperture radar) | |
CN106772224A (en) | A kind of L-type array estimating two-dimensional direction-of-arrival algorithm of use time frequency analysis | |
CN109782243A (en) | Incomplete data MIMO radar angle estimating method based on block Hankel matrix filling | |
CN101241180A (en) | Orthonormal discrete frequency coding design method possessing relative low self correlation performance | |
CN110535537B (en) | Underwater communication and detection integrated method | |
CN109814073A (en) | A kind of MTD radar ambiguity solution speed-measuring method | |
CN106342236B (en) | Based on the airborne radar clutter suppression method of difference beam | |
CN105116384A (en) | Radar orthogonal waveform design method based on frequency modulation and phase modulation of chaotic sequence | |
CN104076351A (en) | Phase-coherent accumulation detection method for high-speed high maneuvering target | |
CN101975939A (en) | Compressive sensing theory-based Doppler ambiguity-resolution processing method | |
CN102928828B (en) | Phase difference estimation method of distributed radar based on orthogonal waveforms | |
CN101702017A (en) | Multi-input multi-output radar waveform design and processing method | |
CN106443595A (en) | Cognition radar waveform design method for resisting instantaneous transmitting slice reconstruction interference | |
CN105675986A (en) | Estimation of arrival angle of narrowband frequency modulation signal based on time frequency analysis during data loss | |
CN102879783A (en) | Sparse detection frequency signal-based inverse synthetic aperture radar (ISAR) imaging method | |
CN104181509A (en) | Incoherent scattering radar signal processing method based on frequency hopping and polyphase alternating codes | |
CN105510895B (en) | A kind of signal processing method of random pulses intervening sequence | |
CN101334469A (en) | Wind profile radar clutter suppression method based on fraction order Fourier transform | |
CN110412514A (en) | A kind of linear frequency modulation continuous wave waveform optimization method under MIMO system | |
CN104991239B (en) | A kind of distance side lobe suppressing method based on pulse compression radar | |
CN109239675A (en) | A kind of distributed spaceborne radar LFM waveform determining method based on same bandwidth difference chirp rate | |
CN104808188A (en) | High-speed stealth target detection method of polynomial Hough Fourier transform | |
CN106597381A (en) | Full coherent full polarization MIMO radar four-channel integrated target detecting method | |
CN102385054A (en) | TPRF method for avoiding velocity measurement obscureness of Doppler weather radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150225 Termination date: 20170321 |
|
CF01 | Termination of patent right due to non-payment of annual fee |