CN113721218A - Heterogeneous radar multiband bandwidth synthesis method and system - Google Patents
Heterogeneous radar multiband bandwidth synthesis method and system Download PDFInfo
- Publication number
- CN113721218A CN113721218A CN202111098818.8A CN202111098818A CN113721218A CN 113721218 A CN113721218 A CN 113721218A CN 202111098818 A CN202111098818 A CN 202111098818A CN 113721218 A CN113721218 A CN 113721218A
- Authority
- CN
- China
- Prior art keywords
- echo signal
- pole
- coherent
- band
- echo
- 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.)
- Granted
Links
- 238000001308 synthesis method Methods 0.000 title description 6
- 230000001427 coherent effect Effects 0.000 claims abstract description 68
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 6
- 238000000265 homogenisation Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000003786 synthesis reaction Methods 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
Images
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
-
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a method and a system for synthesizing a multi-band bandwidth of a heterologous radar. The method comprises the steps of measuring echo signals of different frequency bands of the same target by using different radars; determining a DE model of echo signals of corresponding frequency bands according to the echo signals of different frequency bands; taking any echo signal as a reference echo signal, solving the pole and the amplitude of an all-pole DE model of the echo signal of the residual frequency band by adopting an improved TLS-ESPRIT algorithm, estimating and compensating non-phase parameters contained in the pole and the amplitude, and performing coherent registration on non-coherent multi-segment echo signals; and reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration, and reconstructing a full-frequency-band echo signal by combining an improved TLS-ESPRIT algorithm. The method and the device can improve the accuracy of the data reconstruction of the heterogeneous radar multiband.
Description
Technical Field
The invention relates to the field of radar data reconstruction, in particular to a heterogeneous radar multiband bandwidth synthesis method and system.
Background
The existing heterogeneous radar multiband data reconstruction, namely multiband fusion broadband coherent processing, mainly comprises two main methods, namely a spectrum estimation method and a sparse reconstruction method. Most of spectrum estimation methods are difficult to estimate model orders and sensitive to noise, and a sparse reconstruction method based on compressed sensing has the possibility of failure as a method for reconstructing signals with high probability due to more iteration times and larger calculated amount.
Therefore, a method or system is needed to solve the above problems.
Disclosure of Invention
The invention aims to provide a heterogeneous radar multiband bandwidth synthesis method and system, which can improve the accuracy of data reconstruction of the heterogeneous radar multiband.
In order to achieve the purpose, the invention provides the following scheme:
a heterogeneous radar multiband bandwidth synthesis method comprises the following steps:
measuring echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
determining attenuation index (DE) models of echo signals of corresponding frequency bands according to the echo signals of different frequency bands; the DE model is an all-pole DE model;
taking any echo Signal as a reference echo Signal, solving the pole and the amplitude of an all-pole DE model of the echo signals of the residual frequency band by adopting an improved Total Least square-rotation invariant technical Signal parameter estimation (TLS-ESPRIT) algorithm, Estimating and compensating non-phase parameters contained in the pole and the amplitude, and enabling the non-coherent multi-section echo signals to be coherent and registered;
and reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration, and reconstructing a full-frequency-band echo signal by combining an improved TLS-ESPRIT algorithm.
Optionally, determining a DE model of the echo signal of the corresponding frequency band according to the echo signals of different frequency bands specifically includes:
carrying out spline interpolation homogenization treatment on echo signals of different frequency bands;
and determining a DE model of the echo signals of the corresponding frequency bands according to the processed echo signals of different frequency bands and a geometric diffraction theory.
Optionally, the method for performing coherent registration on multiple echo signals with non-coherent segments includes:
using formulas
Carrying out all-pole DE model coherent registration;
wherein the content of the first and second substances,
for the coherent registered all-pole DE model, s2For the unregistered echo signal, j is the imaginary unit j ═ sqrt (-1), N2Number of frequency samples, n, corresponding to the coherent registered echo2＝0,1,…,N2-1，
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
Optionally, the reconstructing an echo signal by using a Hankel matrix and the echo signal after the coherent registration and reconstructing a full-band echo signal by combining an improved TLS-ESPRIT algorithm specifically includes:
determining parameters of an all-band all-pole DE model by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
according to the parameters of the full-band all-pole DE model, using a formula
Reconstructing a full-band echo signal;
wherein Q is the number of target scattering centers, N is the length of a full-band signal, and diIs the amplitude coefficient, P, of the ith scattering centeriN-1, which is the pole of the ith scattering center, 0, 1.
A heteroradar multiband bandwidth synthesis system, comprising:
the echo signal acquisition module is used for measuring echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
the model determining module is used for determining a DE model of the echo signal of the corresponding frequency band according to the echo signals of different frequency bands; the DE model is an all-pole DE model;
the coherent registration module is used for solving the pole and the amplitude of an all-pole DE model of the echo signals of the residual frequency band by adopting an improved TLS-ESPRIT algorithm by taking any echo signal as a reference echo signal, estimating and compensating non-coherent parameters contained in the pole and the amplitude, and performing coherent registration on non-coherent multi-section echo signals;
and the full-band echo signal reconstruction module is used for reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration and reconstructing the full-band echo signal by combining an improved TLS-ESPRIT algorithm.
Optionally, the echo signal attenuation index and model determining module specifically includes:
the spline interpolation homogenization processing unit is used for carrying out spline interpolation homogenization processing on the echo signals of different frequency bands;
and the attenuation index and model determining unit of the echo signals is used for determining a DE model of the echo signals of the corresponding frequency bands according to the processed echo signals of different frequency bands and a geometric diffraction theory.
Optionally, the coherent registration module specifically includes:
a coherent registration unit for utilizing the formula
Carrying out all-pole DE model coherent registration;
wherein the content of the first and second substances,
for the coherent registered all-pole DE model, s2For the unregistered echo signal, j is the imaginary unit j ═ sqrt (-1), N2Number of frequency samples, n, corresponding to the coherent registered echo2＝0,1,…,N2-1，
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
Optionally, the full-band echo signal reconstructing module specifically includes:
the echo processing data determining unit of the all-band all-pole is used for determining all-band all-pole DE model parameters by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
full-band echo signal reconstruction unit for formula
Reconstructing a full-band echo signal;
wherein Q is the number of target scattering centers, N is the length of a full-band signal, and diIs the amplitude coefficient, P, of the ith scattering centeriN-1, which is the pole of the ith scattering center, 0, 1.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
according to the method and the system for synthesizing the multi-band bandwidth of the heterologous radar, an improved TLS-ESPRIT algorithm is adopted to obtain the attenuation index of an echo signal of a residual frequency band and the pole and the amplitude of a model; the full-band echo signal is reconstructed based on the Hankel matrix and the improved TLS-ESPRIT algorithm, so that the problems of distance image blurring caused by insufficient bandwidth in space target monitoring and identification under the condition of low signal-to-noise ratio, large number of targets in the missile target range or the missile early warning system, different sizes and the like and high resolution difficulty are effectively solved. The improved TLS-ESPRIT-based heterogeneous radar multiband fusion coherent processing method has the characteristics of high precision, high efficiency and good noise resistance. Furthermore, the accuracy of the data reconstruction of the heterogeneous radar multiband is effectively improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without inventive labor.
Fig. 1 is a schematic flow chart of a method for synthesizing a multiband bandwidth of a radar system according to the present invention;
fig. 2 is a schematic structural diagram of a hetero-radar multi-band bandwidth synthesis system provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a heterogeneous radar multiband bandwidth synthesis method and system, which can improve the accuracy of data reconstruction of the heterogeneous radar multiband.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a schematic flow chart of a method for synthesizing a multi-band bandwidth of a radar diversity system according to the present invention, as shown in fig. 1, the method for synthesizing a multi-band bandwidth of a radar diversity system according to the present invention includes:
s101, measuring echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
namely, two radars satisfy the principle of adjacent registration; the proximity registration principle is that the spacing between radars is smaller than the radar-to-target radial distance, and the spacing between radars is much smaller than the radar-to-target radial distance.
S102, determining attenuation indexes and DE models of echo signals of corresponding frequency bands according to the echo signals of different frequency bands; the DE model is an all-pole DE model;
s102 specifically comprises the following steps:
carrying out spline interpolation homogenization treatment on echo signals of different frequency bands;
and determining attenuation indexes and models of the echo signals of the corresponding frequency bands according to the processed echo signals of different frequency bands and a geometric diffraction theory.
The specific establishing process of the attenuation index and the model of the echo signal is as follows:
two radars working independently are set, and the initial carrier frequencies under adjacent registration are respectively f1,f2The number of frequency samples is N1，N2The frequency sampling interval of both is delta f, the blank band width between the initial carrier frequencies is delta B, and f2＝f1+ Δ B. The baseband echoes of two radars of a stationary target consisting of Q scattering points can be represented as:
the two above formulas are respectively s1And s2The Geometric Theory of Diffraction (GTD) model, where j is sqrt (-1), model parameters
Collectively, the characteristics of the Q scattering centers described above are shown. A. theiIs the amplitude of the scattering center, riIs the distance of the scattering center from a reference point, αiThe diffraction term coefficient (shown in Table 1) is an integral multiple of 0.5. f. ofn＝f1+n·Δf,fn'＝f2+ n'. DELTA.f, c is the propagation velocity of the electromagnetic wave in air, here approximately the speed of light. Gamma ray1And gamma2Are respectively s1And s2The non-coherent quantities of the fixed phase term and the linear phase term existing in between.
In a relative bandwidth of not large
When is driven by s1For example, the GTD model may be approximated as follows:
will f isn、fn'And (3) substituting the expressions into the expressions (1) and (2) to obtain attenuation indexes of corresponding echoes and the expressions (4) and (5) of the models, and equivalent to an all-pole model.
Wherein P is1i、P2iAre respectively s1,s2The pole of (a); d1i、d2iAre respectively s1,s2The pole of (a) corresponds to the amplitude coefficient. When two radars are arranged adjacently, the target scattering center intensity, namely the amplitude, can be normalized on a frequency domain. And it is the phase relationship of the signals that mainly affects the imaging quality, so it can be assumed that a is belowi＝1。
TABLE 1
S103, with any echo signal as a reference echo signal, solving the pole and the amplitude of an all-pole DE model of the echo signal of the residual frequency band by adopting an improved total least square-rotation invariant technology signal parameter estimation TLS-ESPRIT algorithm, estimating and compensating non-phase parameters contained in the pole and the amplitude, and enabling the non-phase-coherent multi-section echo signals to be in coherent registration;
s103 specifically comprises the following steps:
using formulas
Carrying out all-pole DE model coherent registration;
wherein the content of the first and second substances,
for the coherent registered all-pole DE model, s2For the unregistered echo signal, j is the imaginary unit j ═ sqrt (-1), N2Number of frequency samples, n, corresponding to the coherent registered echo2＝0,1,…,N2-1，
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
The process of solving the pole and the amplitude of the all-pole DE model of the echo signals of the residual frequency bands by adopting the improved TLS-ESPRIT algorithm comprises the following steps:
definitions x (N) ═ s (N), N ═ 1,2, … N, y (N) ═ s (N + 1). Is provided with
Wherein Q is the number of scattering centers and N is the signal length; introduction of M0X 1-dimensional vector:
X(n)＝[s(n),s(n+1),…,s(n+M0-1)]T
n＝1,2,…N；
Y(n)＝[s(n+1),s(n+2),…,s(n+M0)]T
n＝1,2,…N；
then, the estimated values of the autocorrelation matrix and the cross-correlation matrix thereof can be respectively expressed as:
by using
Direct pair autocorrelation function RXXCarrying out singular value decomposition;
determining the effective order of the model by utilizing the regular correlation technique, and taking the diagonal matrix sigma corresponding to Q main left singular valuesZ1，UZ1And VZ1The matrix is formed by Q main left and right singular vectors, and the product of the vectors is used as a signal subspace. SigmaZ2Diagonal matrix, U, composed of non-dominant singular valuesZ2And VZ2The matrices formed by the non-dominant left and right singular vectors are multiplied to form a noise subspace.
According to the obtained sigma, the minimum singular value is taken and recorded as sigma2Constructing a covariance matrix Cxy＝Rxy-σ2I'；
Using obtained sigmaZ1、
Cxy、VZ1Building a matrix bundle
Solving the generalized eigenvalue of the matrix beam is the pole PiAnalogically determining the pole in each echo segment, and combining the equations
That is, the estimated value of the linear phase term can be obtained
And (3) amplitude estimation process:
according to the obtained model pole
By passing
Find out the
Individual pole and its corresponding amplitude coefficient
In the formula
Is made by
The Vandermonde matrix generated by each pole.
When two radars collect echo signals of the same target under coherent registration, non-coherent quantities between the two radars can be summarized into a linear phase term contained in a pole and a fixed phase term contained in amplitude.
Wherein, γ1' is a fixed phase term, γ2' is a linear phase term. p is a radical of1i、p2iAre respectively s1,s2The pole of (a); d1i、d2iAre respectively s1,s2The pole of (a) corresponds to the amplitude coefficient.
And S104, reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration, and reconstructing a full-band echo signal by combining an improved TLS-ESPRIT algorithm.
S104 specifically comprises the following steps:
determining a full-band all-pole model as follows:
wherein j ═ sqrt (-1), AiFor echo signal amplitude, Q represents the number of scattering centers of the target, αiFor the diffraction term coefficients (as shown in table 1),
is the frequency sampling interval, B is the full-band signal bandwidth, N is the full-band signal length, f0A carrier frequency is initiated for the signal. r isi＝Ri-Rref，RiRadial distance, R, from the i-th scattering center to the radarrefIs a reference distance;
determining parameters of an all-band all-pole DE model by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
according to the parameters of the full-band all-pole DE model, using a formula
Reconstructing a full-band echo signal;
and a vacant frequency band exists between the high frequency band and the low frequency band after the coherent registration, and a new matrix needs to be constructed in order to continue to use the improved TLS-ESPRIT algorithm to complete parameter estimation of the full-frequency band all-pole model. Therefore, an improved TLS-ESPRIT all-band all-pole model solving method based on the Hankel matrix is provided. According to the method, a Hankel matrix is introduced during multi-band echo signal processing, and the influence of noise on signals is greatly reduced by using the special properties of a special structure of the Hankel matrix.
Rearrangement into a form with a Hankel matrix, namely:
where L is often 1/3 the signal length. Composed new matrix
Echo processing data as full band all poles;
reconstructing a full-band echo signal by adopting an improved TLS-ESPRIT algorithm according to the echo processing data of the full-band all-pole and the full-band all-pole model;
and obtaining estimation data of each sub-band and the vacant frequency band:
in order to reduce estimation errors, raw measurement data may be used as much as possible. Thus, the full-band bandwidth extrapolated fusion signal can be expressed as:
the Hankel matrix greatly reduces the influence of noise on signals by utilizing the special property of a special structure of the Hankel matrix; second, the improved TLS-ESPRIT corrects for noise present in both signal subspaces by perturbing the signal subspaces of both data matrices simultaneously and keeping the square of the perturbation norm to a minimum.
Therefore, the full-band echo signal estimated by the algorithm is closer to the original target echo signal acquired by an ideal full-band radar (the initial carrier frequency of a low-band radar is taken as the initial carrier frequency of the full band, and the total synthesized bandwidth is taken as the actual bandwidth of the full-band radar). The method is embodied in that the estimation and compensation precision of the multiband phase to be registered (a linear phase term and a fixed phase term) and the extrapolation data is high, and the calculation efficiency is high.
Fig. 2 is a schematic structural diagram of a hetero-radar multi-band bandwidth synthesis system provided by the present invention, and as shown in fig. 2, the hetero-radar multi-band bandwidth synthesis system provided by the present invention includes:
an echo signal acquiring module 201, configured to measure echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
an attenuation index and model determining module 202 for echo signals, configured to determine a DE model of echo signals of corresponding frequency bands according to echo signals of different frequency bands; the DE model is an all-pole DE model;
the coherent registration module 203 is configured to use any echo signal as a reference echo signal, solve the pole and amplitude of the all-pole DE model of the echo signal of the remaining frequency band by using an improved TLS-ESPRIT algorithm, estimate and compensate non-coherent parameters included in the pole and amplitude, and perform coherent registration on non-coherent multiple echo signals;
and the full-band echo signal reconstruction module 204 is configured to reconstruct an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration, and reconstruct a full-band echo signal by combining with an improved TLS-ESPRIT algorithm.
The echo signal attenuation index and model determining module 202 specifically includes:
the spline interpolation homogenization processing unit is used for carrying out spline interpolation homogenization processing on the echo signals of different frequency bands;
and the attenuation index and model determining unit of the echo signals is used for determining an ED model of the echo signals of the corresponding frequency band according to the processed echo signals of different frequency bands and a geometric diffraction theory.
The coherent registration module 203 specifically includes:
a coherent registration unit for utilizing the formula
Carrying out all-pole DE model coherent registration;
wherein
For the coherent registered all-pole DE model, s1For the unregistered echo signal, j is the imaginary unit j ═ sqrt (-1), N2Number of frequency samples, n, corresponding to the coherent registered echo2＝0,1,…,N2-1， 
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
The full-band echo signal reconstruction module 204 specifically includes:
the echo processing data determining unit of the all-band all-pole is used for determining all-band all-pole DE model parameters by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
full-band echo signal reconstruction unit for formula
Reconstructing a full-band echo signal;
wherein Q is the number of target scattering centers, N is the length of a full-band signal, and diIs the amplitude coefficient of the ith scattering center, PiN-1, which is the pole of the ith scattering center, 0, 1.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (8)
1. A method for synthesizing a multiband bandwidth of a heterologous radar, comprising:
measuring echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
determining attenuation indexes and DE models of echo signals of corresponding frequency bands according to the echo signals of different frequency bands; the DE model is an all-pole DE model;
taking any echo signal as a reference echo signal, solving the pole and the amplitude of an all-pole DE model of the echo signal of the residual frequency band by adopting an improved total least square-rotation invariant technology signal parameter estimation TLS-ESPRIT algorithm, estimating and compensating non-phase parameters contained in the pole and the amplitude, and performing coherent registration on non-coherent multi-segment echo signals;
and reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration, and reconstructing a full-frequency-band echo signal by combining an improved TLS-ESPRIT algorithm.
2. The method according to claim 1, wherein the determining the DE model of the echo signals of the corresponding frequency band according to the echo signals of the different frequency bands specifically comprises:
carrying out spline interpolation homogenization treatment on echo signals of different frequency bands;
and determining a DE model of the echo signals of the corresponding frequency bands according to the processed echo signals of different frequency bands and a geometric diffraction theory.
3. The method according to claim 1, wherein any echo signal is used as a reference echo signal, an improved overall least squares-rotation invariant technology signal parameter estimation TLS-ESPRIT algorithm is used to solve the pole and amplitude of an all-pole DE model of the echo signals of the remaining frequency bands, and non-coherent parameters included in the pole and amplitude are estimated and compensated, so that non-coherent multi-segment echo signals are subjected to coherent registration, specifically comprising:
using formulas
Carrying out all-pole DE model coherent registration;
wherein the content of the first and second substances,
for the coherent registered all-pole DE model, s2For the unregistered echo signal, j is the imaginary unit j ═ sqrt(-1)，N2Number of frequency samples, n, corresponding to the coherent registered echo2＝0,1,…,N2-1，
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
4. The method as claimed in claim 1, wherein the reconstructing echo signals of the reference echo signal and the echo signals after coherent registration by using a Hankel matrix and reconstructing full-band echo signals by combining with an improved TLS-ESPRIT algorithm specifically comprises:
determining parameters of an all-band all-pole DE model by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
according to the parameters of the full-band all-pole DE model, using a formula
Reconstructing a full-band echo signal;
wherein Q is the number of target scattering centers, N is the length of a full-band signal, and diIs the amplitude coefficient, P, of the ith scattering centeriN-1, which is the pole of the ith scattering center, 0, 1.
5. A heteroradar multiband bandwidth synthesis system, comprising:
the echo signal acquisition module is used for measuring echo signals of different frequency bands of the same target by using different radars; different radar initial carrier frequencies are different; the distance between the radars is smaller than the radial distance from the radars to the target; the difference of the visual angles between the radars is between 5 and 10 degrees;
the model determining module is used for determining a DE model of the echo signal of the corresponding frequency band according to the echo signals of different frequency bands; the DE model is an all-pole DE model;
the coherent registration module is used for solving the pole and the amplitude of an all-pole DE model of the echo signals of the residual frequency band by adopting an improved TLS-ESPRIT algorithm by taking any echo signal as a reference echo signal, estimating and compensating non-coherent parameters contained in the pole and the amplitude, and performing coherent registration on non-coherent multi-section echo signals;
and the full-band echo signal reconstruction module is used for reconstructing an echo signal by using a Hankel matrix according to the reference echo signal and the echo signal after the coherent registration and reconstructing the full-band echo signal by combining an improved TLS-ESPRIT algorithm.
6. The hetero-radar multi-band bandwidth synthesis system of claim 5, wherein the echo signal attenuation index and model determination module specifically comprises:
the spline interpolation homogenization processing unit is used for carrying out spline interpolation homogenization processing on the echo signals of different frequency bands;
and the attenuation index and model determining unit of the echo signals is used for determining a DE model of the echo signals of the corresponding frequency bands according to the processed echo signals of different frequency bands and a geometric diffraction theory.
7. The hetero-radar multi-band bandwidth synthesis system of claim 5, wherein the coherent registration module specifically comprises:
a coherent registration unit for utilizing the formula
Carrying out all-pole DE model coherent registration;
wherein the content of the first and second substances,
for the coherent registered all-pole DE model, s2For the unregistered echo signal, j is the imaginary unit j ═ sqrt (-1), N2Are coherentFrequency sampling number n corresponding to the registered echo2＝0,1,…,N2-1，
Is an estimate of the fixed phase term,
is an estimate of the linear phase term.
8. The hetero-radar multi-band bandwidth synthesis system of claim 5, wherein the full-band echo signal reconstruction module specifically comprises:
the echo processing data determining unit of the all-band all-pole is used for determining all-band all-pole DE model parameters by combining a Hankel matrix and a TLS-ESPRIT algorithm according to the reference echo signal and the echo signal after the coherent registration;
full-band echo signal reconstruction unit for formula
Reconstructing a full-band echo signal;
wherein Q is the number of target scattering centers, N is the length of a full-band signal, and diIs the amplitude coefficient, P, of the ith scattering centeriN-1, which is the pole of the ith scattering center, 0, 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111098818.8A CN113721218B (en) | 2021-09-18 | 2021-09-18 | Heterogeneous radar multiband bandwidth synthesis method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111098818.8A CN113721218B (en) | 2021-09-18 | 2021-09-18 | Heterogeneous radar multiband bandwidth synthesis method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113721218A true CN113721218A (en) | 2021-11-30 |
| CN113721218B CN113721218B (en) | 2023-12-08 |
Family
ID=78684367
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111098818.8A Active CN113721218B (en) | 2021-09-18 | 2021-09-18 | Heterogeneous radar multiband bandwidth synthesis method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113721218B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060077091A1 (en) * | 2004-08-06 | 2006-04-13 | Sparta, Inc. | Super-resolution based on frequency domain interferometric processing of sparse multi-sensor measurements |
| US20060238411A1 (en) * | 2004-05-28 | 2006-10-26 | Time Domain Corporation | System and method for spatially diverse radar signal processing |
| EP2533069A1 (en) * | 2011-06-10 | 2012-12-12 | Sony Corporation | Signal processing unit and method |
| CN103412306A (en) * | 2013-07-10 | 2013-11-27 | 北京航空航天大学 | Ultra wideband coherent processing method based on multi-station and multi-band bandwidth fusion |
| CN103675783A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A broadband multiband imaging coherent processing method |
| JP2018132523A (en) * | 2017-02-13 | 2018-08-23 | 国立大学法人電気通信大学 | Discrete frequency synthesis rader device, and distance estimation method and program |
| CN109143151A (en) * | 2018-08-16 | 2019-01-04 | 南京航空航天大学 | The uniform surface battle array tensor reconstructing method and signal source locating method of part array element damage |
| CN111308438A (en) * | 2020-03-06 | 2020-06-19 | 中国人民解放军海军航空大学 | Echo data scattering feature extraction method and system |
| CN111610505A (en) * | 2020-06-19 | 2020-09-01 | 中国人民解放军空军工程大学航空机务士官学校 | Airborne bistatic radar STAP algorithm based on exponential form time-varying weighting |
| US10866304B1 (en) * | 2018-09-12 | 2020-12-15 | Neural Propulsion Systems, Inc. | Signal detection and denoising systems |
| CN112596051A (en) * | 2020-12-11 | 2021-04-02 | 中国科学院空天信息创新研究院 | Multi-band fusion algorithm based on improved relaxation algorithm |
-
2021
- 2021-09-18 CN CN202111098818.8A patent/CN113721218B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060238411A1 (en) * | 2004-05-28 | 2006-10-26 | Time Domain Corporation | System and method for spatially diverse radar signal processing |
| US20060077091A1 (en) * | 2004-08-06 | 2006-04-13 | Sparta, Inc. | Super-resolution based on frequency domain interferometric processing of sparse multi-sensor measurements |
| EP2533069A1 (en) * | 2011-06-10 | 2012-12-12 | Sony Corporation | Signal processing unit and method |
| CN103675783A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A broadband multiband imaging coherent processing method |
| CN103412306A (en) * | 2013-07-10 | 2013-11-27 | 北京航空航天大学 | Ultra wideband coherent processing method based on multi-station and multi-band bandwidth fusion |
| JP2018132523A (en) * | 2017-02-13 | 2018-08-23 | 国立大学法人電気通信大学 | Discrete frequency synthesis rader device, and distance estimation method and program |
| CN109143151A (en) * | 2018-08-16 | 2019-01-04 | 南京航空航天大学 | The uniform surface battle array tensor reconstructing method and signal source locating method of part array element damage |
| US10866304B1 (en) * | 2018-09-12 | 2020-12-15 | Neural Propulsion Systems, Inc. | Signal detection and denoising systems |
| CN111308438A (en) * | 2020-03-06 | 2020-06-19 | 中国人民解放军海军航空大学 | Echo data scattering feature extraction method and system |
| CN111610505A (en) * | 2020-06-19 | 2020-09-01 | 中国人民解放军空军工程大学航空机务士官学校 | Airborne bistatic radar STAP algorithm based on exponential form time-varying weighting |
| CN112596051A (en) * | 2020-12-11 | 2021-04-02 | 中国科学院空天信息创新研究院 | Multi-band fusion algorithm based on improved relaxation algorithm |
Non-Patent Citations (6)
| Title |
|---|
| SHUO ZHANG 等: "Multiband Radar Signal Phase Reference Synthesis Based on Matrix Pencil Method", 2021 IEEE 4TH ADVANCED INFORMATION MANAGEMENT, COMMUNICATES, ELECTRONIC AND AUTOMATION CONTROL CONFERENCE, pages 1271 - 1277 * |
| SHU-YU ZHENG 等: "GTD Model Parameters Estimation Based on Improved LS-ESPRIT Algorithm", 2019 PHOTONICS & ELECTROMAGNETICS RESEARCH SYMPOSIUM - FALL, pages 2282 - 2289 * |
| 张小宽;郑舒予;奚之飞;葛启超;宗彬锋;: "基于改进LS-ESPRIT算法的GTD模型参数估计与RCS重构", 电子与信息学报, no. 10, pages 2493 - 2499 * |
| 张硕 等: "多频带带宽相参合成技术研究现状", 《兵器装备工程学报》, vol. 42, no. 8, pages 147 - 152 * |
| 张静;徐政;牛卢璐;: "TLS-ESPRIT在电力系统信号高精度频谱估计中应用", 电力自动化设备, no. 05, pages 48 - 51 * |
| 李尚生 等: "基于Hankel矩阵改进TLS-ESPRIT算法的散射中心参数提取及RCS重构", 系统工程与电子技术, vol. 43, no. 1, pages 62 - 71 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113721218B (en) | 2023-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sun et al. | Advanced signal processing methods for ground-penetrating radar: Applications to civil engineering | |
| CN111142105A (en) | ISAR imaging method for complex moving target | |
| Li et al. | Scaled Radon-Wigner transform imaging and scaling of maneuvering target | |
| CN108037494B (en) | Radar target parameter estimation method under impulse noise environment | |
| CN105807241B (en) | A kind of exponential signal denoising method using prior information | |
| CN105137373A (en) | Index signal de-noising method | |
| CN102901966B (en) | Moving-target instantaneous distance Doppler imaging method based on detuning and keystone processing | |
| CN112255595A (en) | Sea clutter data preprocessing method based on simulated airborne measurement | |
| Wang et al. | Prior-guided deep interference mitigation for FMCW radars | |
| CN112596051A (en) | Multi-band fusion algorithm based on improved relaxation algorithm | |
| Lagunas et al. | Compressive sensing for through wall radar imaging of stationary scenes using arbitrary data measurements | |
| CN113721218A (en) | Heterogeneous radar multiband bandwidth synthesis method and system | |
| CN109521411A (en) | A kind of detection method of range extension target | |
| Chen et al. | Examination of Antenna Calibration Methodologies in an Extrapolation Range | |
| CN115453530B (en) | Double-base SAR filtering back projection two-dimensional self-focusing method based on parameterized model | |
| Sandström | Implementation of FMCW radar at low frequencies | |
| CN113189570B (en) | Array signal processing method and system based on complex domain compressed sensing | |
| CN115469286A (en) | Super-resolution angle measurement method based on millimeter wave automobile radar minimum redundancy MIMO array | |
| CN114325700A (en) | Satellite-borne multi-channel SAR moving target imaging method | |
| CN113655443A (en) | Low-frequency band SAR radio frequency interference suppression method for broadband digital television signal | |
| Tran et al. | An experimental study of radar tomographic imaging in a multi-bistatic scenario | |
| CN114609606B (en) | Parameter estimation method and device of target object and electronic equipment | |
| CN111024305B (en) | Method for detecting vacuum degree by using THz signal | |
| Pardini et al. | Tomosar vertical profiles at multiple sar frequencies: An experimental analysis | |
| JP2690606B2 (en) | Sound source number determination method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |