CN105093191B - Direct wave recovery method for external radiation source radar - Google Patents
Direct wave recovery method for external radiation source radar Download PDFInfo
- Publication number
- CN105093191B CN105093191B CN201510502616.3A CN201510502616A CN105093191B CN 105093191 B CN105093191 B CN 105093191B CN 201510502616 A CN201510502616 A CN 201510502616A CN 105093191 B CN105093191 B CN 105093191B
- Authority
- CN
- China
- Prior art keywords
- diversity
- signal
- mean square
- square error
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/42—Diversity systems specially adapted for radar
Abstract
The invention discloses a direct wave recovery method for external radiation source radar. The method comprises: S101, receiving, by a reference channel of the external radiation source radar, diversity signals of a plurality of diversity channels through a space diversity method; S102, processing the diversity signals of each diversity channel by adopting a one-path equalizer to obtain multi-path equalizer output signals; S103, summating the multi-path equalizer output signals to obtain direct wave signals; S104, calculating an error term according to the direct wave signals; S105, calculating a mean square error according to the error term, and judging whether the mean square error is smaller than a preset threshold; S106, if the mean square error is not smaller than the preset threshold, returning to S104 to calculate next error term; and S107, if the mean square error is smaller than the preset threshold, iterating the vectors of the equalizers every Ng sampling points. The method is small in calculation quantity and easy to implement in engineering.
Description
Technical field
The present invention relates to Radar Technology field, more particularly to a kind of external illuminators-based radar direct wave restoration methods.
Background technology
Direct wave recovers the key technology for being frequency modulation broadcasting external illuminators-based radar.As system is near cities or mountain area ring
The application in border, serious environment clutter causes external illuminators-based radar direct wave to recover to become extremely difficult.Traditional direct wave is extensive
The main lobe alignment direct-path signal direction of second mining Wave beam forming, using direct wave antenna side lobe the suppression to multipath signal is realized
System.In the case of heavy multipath, multipath clutter is entered from main lobe, causes distorted signals.It is same using constant mould Blind Equalization Technique
Can realize that direct wave is purified, including the multipath conditions for entering main lobe simultaneously.In the case of heavy multipath, constant mould method
Constringency performance is poor.And space diversity blind balance method is adapted to the situation of heavy multipath, strengthen the recovery capability of direct wave.
But with the increase of space diversity dimension, the amount of calculation of algorithm is greatly increased, and is unfavorable for the Project Realization of the method.
The content of the invention
The technical problem to be solved is, there is provided a kind of external illuminators-based radar direct wave restoration methods, the party
Method devises the less interval alternative manner of amount of calculation for the blind balance method of space diversity, realizes strong multipath clutter ring
The recovery of direct wave under border.
In order to solve above-mentioned technical problem, the invention provides a kind of external illuminators-based radar direct wave restoration methods include：
S101：The reference channel of external illuminators-based radar receives the diversity letter of multiple diversity passages using space diversity method
Number；
S102：Processed using balanced device all the way for the diversity signal of each diversity passage, it is equal so as to obtain multichannel
Weighing apparatus output signal；Wherein,In formula, piK () represents the i-th road equalizer output signal, xi(k)
Represent the diversity signal of i-th diversity passage of current time, fiK () is Nf× 1 is vectorial, represents the vector of No. i-th balanced device, fi l
K () represents l-th equalizer coefficients of No. i-th balanced device,I=
0,1 ..., N-1, N represent the number of balanced device；
S103：Multichannel equalizer output signal is sued for peace, direct-path signal is obtained；Wherein,Z (k) represents direct-path signal；
S104：Error term is calculated according to the direct-path signal；Wherein, e (k)=z (k) [γ-z2(k)], in formula,
E (k) represents error term, and γ represents the kurtosis of primary signal a (k) that external illuminators-based radar is launched, and
S105：Mean square error is calculated according to the error term, and judges the mean square error whether less than pre- gating
Limit；Wherein, the mean square error is σ2(k+1)=λ σ2(k)+e(k)e*K (), in formula, λ is forgetting factor, initial value σ2(0)=
0；
S106：If the mean square error is not less than the pre-determined threshold, returns S104 and calculate next error term e (k+1)；
S107：If the mean square error is less than the pre-determined threshold, every NgIndividual sampled point, to per balanced device all the way
Vector is iterated；Wherein, vectorial iterative formula is fi(k+1)=fi(k)+μ·Xi *(k) e (k), in formula, Xi(k)=[xi(k)
xi(k-1) … xi(k-l) … xi(k-Nf+1)]T, XiK () represents the input vector of i-th diversity passage, xi(k-l) represent
L-th diversity signal of i-th diversity passage, μ is step factor.
Implement the present invention, have the advantages that：The external illuminators-based radar direct wave restoration methods of the present invention are for sky
Between diversity blind balance method, devise the less interval alternative manner of amount of calculation, realize under strong multipath clutter environment go directly
The recovery of ripple, it is easy to Project Realization.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is 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 schematic flow sheet of the one embodiment for the external illuminators-based radar direct wave restoration methods that the present invention is provided；
Fig. 2 is the theory diagram of Fig. 1；
Fig. 3 is that the flow process of the emulation verification method of the external illuminators-based radar direct wave restoration methods provided the present invention is illustrated
Figure；
Fig. 4 is the theory diagram of Fig. 3；
Fig. 5 is the theory diagram of S201；
Fig. 6 is SD-CMA algorithm output signal modulus value；
Fig. 7 is the schematic diagram of SD-CMA algorithm range Doppler planes；
Fig. 8 is the IISD-CMA algorithm output signal modulus value at intervals of 2；
Fig. 9 is the range Doppler figure of the IISD-CMA algorithms at intervals of 2
Figure 10 is the IISD-CMA algorithm output signal modulus value at intervals of 3；
Figure 11 is the range Doppler figure of the IISD-CMA algorithms at intervals of 3.
Specific 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 the embodiment of whole.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.
Fig. 1 is the schematic flow sheet of the one embodiment for the external illuminators-based radar direct wave restoration methods that the present invention is provided,
Fig. 2 is the theory diagram of Fig. 1, as depicted in figs. 1 and 2, including：
S101：The reference channel of external illuminators-based radar receives the diversity letter of multiple diversity passages using space diversity method
Number；
S102：Processed using balanced device all the way for the diversity signal of each diversity passage, it is equal so as to obtain multichannel
Weighing apparatus output signal；Wherein,In formula, piK () represents the i-th road equalizer output signal, xiK () represents
The diversity signal of i-th diversity passage of current time, fiK () is Nf× 1 is vectorial, represents the vector of No. i-th balanced device, fi l(k)
L-th equalizer coefficients of No. i-th balanced device are represented,I=0,
1 ..., N-1, N represent the number of balanced device；
S103：Multichannel equalizer output signal is sued for peace, direct-path signal is obtained；Wherein,z
K () represents direct-path signal；
S104：Error term is calculated according to the direct-path signal；Wherein, e (k)=z (k) [γ-z2(k)], in formula,
E (k) represents error term, and γ represents the kurtosis of primary signal a (k) that external illuminators-based radar is launched, and
S105：Mean square error is calculated according to the error term, and judges the mean square error whether less than pre- gating
Limit；Wherein, the mean square error is σ2(k+1)=λ σ2(k)+e(k)e*K (), in formula, λ is forgetting factor, initial value σ2(0)=
0；
S106：If the mean square error is not less than the pre-determined threshold, returns S104 and calculate next error term e (k+1)；
S107：If the mean square error is less than the pre-determined threshold, every NgIndividual sampled point, to per balanced device all the way
Vector is iterated；Wherein, vectorial iterative formula is fi(k+1)=fi(k)+μ·Xi *(k) e (k), in formula, Xi(k)=[xi(k)
xi(k-1) … xi(k-l) … xi(k-Nf+1)]T, XiK () represents the input vector of i-th diversity passage, xi(k-l) represent
L-th diversity signal of i-th diversity passage, μ is step factor.
Below by the correctness of the Computer Simulation signal authentication embodiment of the present invention：
As shown in Figure 3 and Figure 4, the step of emulation verification method includes：
S201：External illuminators-based radar is sent primary signal to multiple diversity passages by space diversity method.
Specifically, as shown in figure 5, the S201 is specifically included：
External illuminators-based radar is sent primary signal to multiple diversity passages by space diversity method, wherein, i-th point
Collection passage on diversity signal beIn formula, a (k) is external illuminators-based radar transmitting
Primary signal, niK () is the noise at i-th diversity channel reception end, hiK () is Lr× 1 channel response vector, and Represent that i-th diversity passage is believed for l-th
Road coefficient.
S202：The external illuminators-based radar direct wave restoration methods provided using the present invention recover direct-path signal.
S203：Clutter cancellation is carried out using least mean square algorithm according to monitoring channel signal and the direct-path signal, is obtained
Signal to after offseting.
Specifically, the S203 is specifically included：
Monitoring channel signal y (k) is produced, wherein,In formula, nsK () is to connect
Receipts machine noise, c is multipath clutter, andLsFor multipath number, g (k) is target echo signal, andfdmFor bistatic Doppler, ξ is bistatic time delay, fsFor sample frequency, MTFor mesh
Mark number；
Adaptive-filtering is carried out to the direct-path signal, wherein, wave filter input recursive vector is Z (k)=[z (k) z
(k-1) … z(k-l) … z(k-Nw+1)]T, z (k-l) is the input recurrence direct-path signal of l-th time delay in formula；
Clutter cancellation is carried out using least mean square algorithm, signal u (k) after being offseted, wherein, u (k)=y (k)-ZH
(k) V (k), in formulaFor weight vector, vlK () is the
L weight coefficient, weight vector iterative formula is V (k+1)=V (k)+η Z (k) u*K (), η is step factor.
S204：Range Doppler plane is built according to the signal after the suppression.
Specifically, the S204 is specifically included：
Range Doppler plane is built according to the signal after the suppression, wherein, the range Doppler plane is：NcFor correlation length, ξ represents time delay, and m represents Doppler frequency shift.
S205：By target location at the range Doppler plane monitoring-network.
Wherein, simulation parameter is specially：System bandwidth is 100kHz, and data transfer rate is 200kHz.It is long that the time of integration is set
Spend for 2s.Balanced device weight vector length is set to 41, and iteration step length is 0.001, the bistatic time delay of two targets and Doppler
Respectively (171km, -300Hz), (276km, 60Hz).Clutter noise ratio is -58dB, and signal to noise ratio (SNR) is set to -3dB.With
Least mean square algorithm (LMS algorithm) the weight vector length that work is offseted is set to 20, and iteration step length is set to 0.001.Space diversity
Dimension is set to 3, and reference channel is provided with 3 diversity subchannels, is shown in Table 1, and by the method for the present invention, that is, is spaced iterative space
Between diversity constant mould blind balance method be designated as IISD-CMA, Conventional spatial diversity constant mould blind balance method is designated as into SD-CMA..
The channel tap coefficient being not zero is listed in table 1.The channel multipath coefficient of monitoring passage is shown in Table the channel c in 1.
5 three footpath channels used by the emulation of table 11
Emulation main contrast SD-CMA is more with the convergence rate in direct wave recovery process of IISD-CMA methods, distance
General Le plane monitoring-network performance etc. compares.SD-CMA method constringency performances are shown in Fig. 6.Modulus value tends to constant and shows that algorithm successfully recovers
Modulus value is the frequency modulation broadcasting reference signal of constant.The fluctuating of modulus value indicates from large to small convergence rate.Range Doppler plane
Target detection performance see Fig. 7, two targets are high-visible.
Iteration interval takes the constringency performance of 2 IISD-CMA methods as shown in Figure 8, and amount of calculation declines 2 compared with without interval arithmetic
Times.It can be seen that the convergence rate of the method is not decreased obviously.The target detection performance of range Doppler plane is shown in Fig. 9, it is seen that with
Corresponding two targets of SD-CMA methods of Fig. 5 are high-visible, and signal to noise ratio is close to.
Change the interval of interval alternative manner, take at intervals of 3, amount of calculation declines 3 times compared with without interval arithmetic.Algorithm output letter
Number modulus value change curve is shown in Figure 10.Compared with Fig. 6, convergence rate is declined slightly.Target detection performance is shown in Figure 11, without significant performance
Loss.
Implement the present invention, have the advantages that：The external illuminators-based radar direct wave restoration methods pin that the present invention is provided
Blind balance method to space diversity, devises the less interval alternative manner of amount of calculation, under realizing strong multipath clutter environment
The recovery of direct wave, it is easy to Project Realization.
It should be noted that herein, term " including ", "comprising" or its any other variant are intended to non-row
His property is included, so that a series of process, method, article or device including key elements not only include those key elements, and
And also include other key elements being not expressly set out, or also include for this process, method, article or device institute inherently
Key element.In the absence of more restrictions, the key element for being limited by sentence "including a ...", it is not excluded that including being somebody's turn to do
Also there is other identical element in the process of key element, method, article or device.
In several embodiments provided herein, it should be understood that disclosed system and method can pass through it
Its mode is realized.For example, system embodiment described above is only schematic, for example, the division of the unit, and only
Only a kind of division of logic function, can there is other dividing mode when actually realizing, such as multiple units or component can be tied
Close or be desirably integrated into another system, or some features can be ignored, or do not perform.It is another, it is shown or discussed
Coupling each other or direct-coupling or communication connection can be the INDIRECT COUPLINGs by some interfaces, device or unit or logical
Letter connection, can be electrical, mechanical or other forms.
Professional further appreciates that, with reference to the unit of each example of the embodiments described herein description
And algorithm steps, can with electronic hardware, computer software or the two be implemented in combination in, in order to clearly demonstrate hardware and
The interchangeability of software, according to function has generally described the composition and step of each example in the above description.These
Function is performed with hardware or software mode actually, depending on the application-specific and design constraint of technical scheme.Specialty
Technical staff can use different methods to realize described function to each specific application, but this realization should not
Think beyond the scope of this invention.
The step of method described with reference to the embodiments described herein or algorithm, directly can be held with hardware, processor
Capable software module, or the combination of the two is implementing.Software module can be placed in random access memory (RAM), internal memory, read-only deposit
Reservoir (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or technology
In field in known any other form of storage medium.
The foregoing description of the disclosed embodiments, enables professional and technical personnel in the field to realize or using the present invention.
Various modifications to these embodiments will be apparent for those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, the present invention
The embodiments shown herein is not intended to be limited to, and is to fit to and principles disclosed herein and features of novelty phase one
The most wide scope for causing.
Claims (1)
1. a kind of external illuminators-based radar direct wave restoration methods, it is characterised in that include：
S101：The reference channel of external illuminators-based radar receives the diversity signal of multiple diversity passages using space diversity method；
S102：Processed using balanced device all the way for the diversity signal of each diversity passage, so as to obtain multichannel balanced device
Output signal；Wherein,In formula, piK () represents the i-th road equalizer output signal, xiK () represents
The diversity signal of i-th diversity passage of current time, fiK () is Nf× 1 is vectorial, represents the vector of No. i-th balanced device,
L-th equalizer coefficients of No. i-th balanced device are represented,N
The number of balanced device is represented, k is discrete-time variable, NfFor fiThe length of (k) vector；
S103：Multichannel equalizer output signal to be sued for peace, and it is average carry out again, obtains direct-path signal；Wherein,Z (k) represents direct-path signal；
S104：Error term is calculated according to the direct-path signal；Wherein, e (k)=z (k) [γ-z2(k)], in formula, e (k)
Error term is represented, γ represents the kurtosis of primary signal a (k) that external illuminators-based radar is launched, and
S105：Mean square error is calculated according to the error term, and judges the mean square error whether less than pre-determined threshold；Its
In, the mean square error is σ2(k+1)=λ σ2(k)+e(k)e*K (), in formula, λ is forgetting factor, initial value σ2(0)=0；
S106：If the mean square error is not less than the pre-determined threshold, returns S104 and calculate next error term e (k+1)；
S107：If the mean square error is less than the pre-determined threshold, every NgIndividual sampled point, to the vector per balanced device all the way
It is iterated；Wherein, vectorial iterative formula isIn formula, Xi(k)=[xi(k)
xi(k-1) … xi(k-l) … xi(k-Nf+1)]T, XiK () represents the input vector of i-th diversity passage, xi(k-l) represent
L-th diversity signal of i-th diversity passage, μ is step factor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510502616.3A CN105093191B (en) | 2015-08-14 | 2015-08-14 | Direct wave recovery method for external radiation source radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510502616.3A CN105093191B (en) | 2015-08-14 | 2015-08-14 | Direct wave recovery method for external radiation source radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105093191A CN105093191A (en) | 2015-11-25 |
| CN105093191B true CN105093191B (en) | 2017-05-03 |
Family
ID=54574056
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510502616.3A Active CN105093191B (en) | 2015-08-14 | 2015-08-14 | Direct wave recovery method for external radiation source radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105093191B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106970382B (en) * | 2017-03-22 | 2019-08-02 | 武汉大学 | One kind being based on external illuminators-based radar unmanned plane real-time monitoring system and method |
| CN113395779A (en) | 2017-03-24 | 2021-09-14 | 中兴通讯股份有限公司 | Processing method and device for beam recovery |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7496619B2 (en) * | 2002-06-18 | 2009-02-24 | Vanderbilt University | System and methods of nonuniform data sampling and data reconstruction in shift invariant and wavelet spaces |
| CN101651500B (en) * | 2009-09-18 | 2013-02-27 | 南京信息工程大学 | Improved T/2 fractional spaced coordinate transforming blind equalization method |
| CN101656696B (en) * | 2009-09-18 | 2012-03-07 | 南京信息工程大学 | Frequency-domain small wave blind equalization method based on united combining space-time diversity |
| FR2959895B1 (en) * | 2010-05-07 | 2012-08-03 | Commissariat Energie Atomique | COOPERATIVE NETWORK WITH STRATEGY FOR ADAPTIVE RETRANSMISSION REQUEST |
| CN104777466B (en) * | 2015-03-23 | 2017-05-03 | 河海大学 | Passive bistatic radar signal processing method based on dimension reduction constant modulus blind equalization |
-
2015
- 2015-08-14 CN CN201510502616.3A patent/CN105093191B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105093191A (en) | 2015-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104777466B (en) | Passive bistatic radar signal processing method based on dimension reduction constant modulus blind equalization | |
| CN103760519B (en) | High resolution DOA estimation method and system | |
| CN101388677A (en) | Interference resisting method for communication receiver based on blind signal separation and system thereof | |
| CN105207964B (en) | A kind of underwater sound adaptive decision-feedback equalization method based on single vector sensor | |
| CN108521291B (en) | Radio frequency vortex wavefront optimization device and method | |
| CN102231280A (en) | Frequency-domain blind separation sequencing algorithm of convolutive speech signals | |
| CN101729461A (en) | System and method for eliminating single-frequency interference and multi-frequency interference | |
| CN103346756A (en) | Space-time self-adaptive filtering method and device | |
| CN105093191B (en) | Direct wave recovery method for external radiation source radar | |
| CN107037409B (en) | MIMO radar waveform separation method based on compressed sensing | |
| CN104537257A (en) | Distributed self-adaptation direct positioning method based on time difference | |
| CN103368885B (en) | The fusion method that a kind of frequency domain bidirectional iteration is balanced | |
| CN112910533B (en) | Broadband signal array system with parallel structure | |
| CN104716982A (en) | Robust anti-interference processing method and device of spread frequency system | |
| CN104749562B (en) | Passive bistatic radar signal processing method based on cyclic constant mould blind equalization | |
| Wijenayake et al. | Space-time spectral white spaces in cognitive radio: Theory, algorithms, and circuits | |
| CN105978833A (en) | Improved SAGE channel parameter estimation method | |
| CN103630880B (en) | For the waveform diversity method of array SAR three-dimensional imaging | |
| CN106026974B (en) | A kind of passband global response error constraints airspace matrix filter design method | |
| CN102355435A (en) | Wavelet weighted multi-modulus blind equalization algorithm based on fractional lower order statistics (WT-FLOSWMMA) | |
| Zhao et al. | Receiver-channel based adaptive blind equalization approach for GPS dynamic multipath mitigation | |
| Walnut et al. | Cornerstones of sampling of operator theory | |
| Kim et al. | HFM design for timing synchronization in underwater communications systems | |
| CN105334502A (en) | Radar time space two dimension self-adaption processing method based on FIR filters | |
| CN109085575B (en) | Signal processing method of OFDM radar communication integrated airborne platform system |
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 |