CN109901166A - A method of it quickly generates and simulation ISAR echo data - Google Patents

A method of it quickly generates and simulation ISAR echo data Download PDF

Info

Publication number
CN109901166A
CN109901166A CN201910245260.8A CN201910245260A CN109901166A CN 109901166 A CN109901166 A CN 109901166A CN 201910245260 A CN201910245260 A CN 201910245260A CN 109901166 A CN109901166 A CN 109901166A
Authority
CN
China
Prior art keywords
echo
isar
scattering center
frequency
target
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.)
Pending
Application number
CN201910245260.8A
Other languages
Chinese (zh)
Inventor
尹灿斌
李智
徐灿
方宇强
林财永
许洁平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Original Assignee
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peoples Liberation Army Strategic Support Force Aerospace Engineering University filed Critical Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority to CN201910245260.8A priority Critical patent/CN109901166A/en
Publication of CN109901166A publication Critical patent/CN109901166A/en
Pending legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention provides a kind of methods quickly generated with simulation ISAR echo data, comprising: reads in the scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter and deposits into register;Register parameters are read parallel by parallel computing module, different scattering center echo datas are completed to calculate, it is quickly generated from frequency domain and the ISAR solution line tune echo data that envelope goes tilting processing with residual video phase correction process is completed, and different scattering center echo datas are summed and are exported to RAM memory port.The present invention is disposably unified in frequency domain and realizes that the operation of solution line tune, envelope go tilting processing, residual video phase correction etc. are crucial to calculate and simulation process, pass through the thousands of target scattering center echo datas of the frequency information parallel computation that can uniquely divide, substantially increase the speed that ISAR echo data is generated with simulated, meet the higher analogue simulation demand of real-time, while enormously simplifying the hardware system structure of analogue echoes process and simulator.

Description

A method of it quickly generates and simulation ISAR echo data
Technical field
The present invention relates to radar system modelings and (half is in kind/in kind) imitation technology field, while still belonging to signal Process field more particularly to a kind of method quickly generated with simulation ISAR echo data.
Background technique
Inverse Synthetic Aperture Radar (Inverse Synthetic Aperture Radar, ISAR) space situation awareness, Play an important role in the military and civilians such as air-defense anti-missile and civil aviation control, foreign countries had developed ground, bank base, sea base and The multi-platform ISAR system such as airborne.Wherein, the ISAR system of the U.S. and Germany is most representative.
1970, under DARPA approval, america's MIT Lincoln laboratory is built on Kwajalein to be completed and puts into make With first remote high performance wideband imaging radar ALCOR, radar works in C-band, bandwidth 512MHz, Range resolution Rate is 50 centimetres, and when the same year China emits Dongfanghong satellite, the radar is once by booster rocket one-dimensional range profile sequence Its dimension information is extrapolated in analysis.1973, ALCOR was once observed imaging to Skylab satellite, and helped NASA analysis should The damaged condition of satellite, convenient for the subsequent maintenance of the satellite.
Since ALCOR antenna size is relatively small, it is only capable of that low earth orbital target is observed and is imaged.1978, Lincoln laboratory upgrades the Haystack radar that the sixties in last century builds, and working frequency range rises to X-band, bandwidth Reach 1GHz, distance resolution is 25 centimetres, can be observed to the geo-synchronous orbit satellite target outside about 36000 kilometers And imaging.Later in order to adapt to the high resolution observations to microsatellite target, Lincoln laboratory again to Haystack radar into It has gone upgrading, has been allowed to the ability to work for having W-waveband (92GHZ-100GHZ), bandwidth reaches 8GHz, which is thus referred to as Ultra wide band satellite imagery radar (Haystack Ultra wideband Satellite Imaging Radar, HUSIR).By The ISAR radar system development history of above-mentioned U.S.'s Lincoln laboratory can be seen that, ALCOR, Haystack backup radar, MMW radar and The bandwidth of HUSIR system incrementally increases, and resolution ratio gradually increases, and radar carrier frequency is also in rising trend.Can more intuitively it show Big bandwidth and high-resolution can get stronger target identification and descriptive power out.Therefore, high frequency, big bandwidth and high-resolution are The trend of ISAR system and function development.
The tracking of German high-frequency physical and Radar Technology research institute and imaging radar (Tracking And Imaging Radar, TIRA) it is another more famous ground ISAR system, when which works in L-band, tracked for narrow-band impulse Radar;Work in Ku wave band be used to be imaged, at present bandwidth oneself reach 2.1GHz.By ISAR imaging results, can be obtained The features such as kinematic parameter, shape, size, quality and the material property of extraterrestrial target, to improve the space situation awareness of radar Ability.
Russia, Canada and Japan and other countries are also owned by the ISAR system of oneself.The boat embarked such as Russia Its monitoring system (Space Surveillance System, SSS) and space junk radar detection net, Kyoto Univ Japan MUF (Middle and Upper Atmosphere) radar, aerospace research Uchinoura radar and Usuda radar etc. It can be used to be observed extraterrestrial target.
Bank base Cobra Dane and carrier-borne Cobra in addition to ground ISAR system, in U.S.'s ballistic missile defense system Judy radar system, two radars are built by Raytheon Co., and Lincoln laboratory is responsible for technical support and the data processing in later period, Cobra Judy is then made of a phased array S-band radar and parabola antenna X-band radar.Lincoln laboratory was in 1999 Development completes the dual-purpose Cobra Gemini of land-sea and moves radar, is mainly used to track tactical ballistic missile and seen It surveys.
In addition to the ISAR system of above-mentioned ground and sea base, the airborne Ocean Master400 of France and German joint study Radar, Raytheon Co., the U.S. are 5 system of AN/APS-137B (V) that USN develops, the AN/ that Faxforward Inc., the U.S. develops APS-147 surveillance radar and airborne multimode surveillance radar device APS-143B (V) 3 " ocean eye " and american lockheed ALQ-147 multifunction radar of Martin Corporation development etc. all has ISAR imaging function.There are also the U.S. F-35 and F-22 etc. the 4th For AN/APG-77, AN/APG-79 and AN/APG-81FF equipped in fighter plane etc. can also to aerial target carry out ISAR at Picture, and have certain noncooperative target recognition capability.
ISAR can not only obtain the one-dimensional of target, two and three dimensions high-resolution imaging as a result, also fixed using image The processing such as mark, interpretation and identification obtain the characteristic informations abundant such as size and the structure of target, anti-to Situation Awareness, air defense in sky It leads and plays an important role in the military and civilians fields such as civil aviation control.Current ISAR imaging h ardware equipment and basic theory The research of aspect has become mature, but measured data resource is valuable, and the background of Military Application is strong, greatly limits its measured data and exists Application range in target identification work.
For this problem, the ISAR echo generation based on radar echo simulator is come into being with analogue technique, and is taken Obtained very big progress.Currently, various to be emerged one after another based on software, material object or half simulator in kind, ISAR is alleviated significantly Measured data resource is valuable, and military background is strong, is not easy to the problem promoted and applied.But as target becomes increasingly complex, ISAR The resolution ratio of imaging is higher and higher, and signal bandwidth is increasing, and traditional broadband ISAR analogue echoes are also faced with very big choose War.
The distance resolution of ISAR image is related with the bandwidth of radar emission signal, and transmitted signal bandwidth is bigger, ISAR figure The distance resolution of picture is higher.The working frequency range of ISAR is higher, and wavelength is shorter, under the conditions of same target observation corner, The azimuth resolution of ISAR imaging is higher.Meanwhile the working frequency range of ISAR is higher, to the texture information descriptive power of target Also stronger, more conducively identify observed object.Therefore, developing higher working frequency range and the ISAR system of more large signal bandwidth is skill The main trend of art development.It is external at present that oneself has W frequency range the ISAR even relevant report of THz frequency range ISAR.Such as Lincoln laboratory The HUSIR system of construction works in W-waveband, and bandwidth reaches 8GHz, can get the Precise imaging of observed object as a result, having very strong Target decomposition and recognition capability.THz frequency range is between microwave and far infrared frequency range, closer to light in electromagnetic spectrum It learns, thus, for ISAR imaging results closer to optical imagery effect, the textural characteristics such as details of target are more abundant.Meanwhile There is microwave frequency band round-the-clock, all weather operations can be again by THz frequency range ISAR, in terms of anti-stealthy and fine motion Also there is huge application potential.But relative to traditional low band ISAR, the development of high band ISAR is also faced with Problem.Firstly, higher working frequency range means the movement and phase that radar is higher to the required precision of motion compensation, to target Error is more sensitive, and traditional imaging algorithm based on low-frequency range will be not directly applicable among high band ISAR imaging.Its Secondary, due to high band ISAR, especially constraint of the THz radar by current device level, transmitting signal is deposited with transceiver channel In biggish nonlinear phase error, the coherence of the Range compress processing and azimuth dimension of target echo will affect, and then reduce The quality of ISAR imaging.Therefore, for high band ISAR imaging contexts, it is necessary to which the accurate effective phase error compensation of research is calculated Method promotes the high-resolution imaging ability of high band ISAR.With the upsurge and satellite of the space development of our times various countries, critical The rapid development of the space platforms manufacturing technology such as spacecraft, the value volume and range of product of extraterrestrial target increases year by year, and its size It is then smaller and smaller.The space junk target for caused by the reasons such as the disintegration of part satellite, failure and collision needs to carry out them Cataloguing and effective management, threaten to avoid to in-orbit extraterrestrial target.Meanwhile with the extensive application of miniature self-service device, Serious influence is caused on the flight safety in low latitude, needs that such miniature self-service device is detected and managed.Due to above-mentioned Small size weak signal target extraterrestrial target backscattering coefficient Radar Cross Section it is minimum, moreover, the wave band of ISAR it is wider, Operating distance farther out, and needs the modulation in view of electromagnetic propagation environment to influence, should in this way under the conditions of certain transmission power The noise of class target echo is relatively low, will challenge to existing ISAR translational compensation with imaging algorithm.Therefore, for low letter The research than lower ISAR imaging algorithm of making an uproar is a research hotspot and difficulties for current ISAR imaging field.Also, this is asked The research of topic will promote detection and imaging capability of the radar to Weak target, and then improve the situational awareness of ISAR.
More stringent requirements are proposed to ISAR echo simulator for above-mentioned technical need and latest development.
In recent years, point target simulation and its simulator technology gradually tend to be mature.To inspection ISAR system performance and respectively The superiority and inferiority of kind imaging algorithm is of great significance to.On the basis of simulating to ideal point target, people are more wished to Imaging of the ISAR to complicated maneuvering target is simulated, to measure some system indexs.During complex target analogue echoes, need Establish complicated Electromagnetic scattering model.To improve accuracy, also commonly used in simulated target backscattering coefficient Kirchhoff model can calculate more accurate facet according to the geometrical relationship of small plane unit and incident electromagnetic wave The backscattering coefficient of unit.
ISAR echo semi-physical object simulating platform includes two hardware, software parts, by ISAR analogue echoes data producer, D/A transformation and orthogonal modulation composition.Wherein ISAR analogue echoes data producer is the key that generate analogue echoes data, is needed It is special to develop, require it with rapid computations, real-time parallel processing, the ability for meeting exchange data bandwidth on hardware.In and Frequently, it is then to have borrowed the material object for emitting data feedback channel in ISAR receiver that up-conversion and pumping signal, which generate part,.
On software, ISAR Echo searching generator will be according to aerial or the space parameters of target motion, rail being modeled Mark, spatial observation are several how to carry out necessary analog parameter calculating, and follows strictly the basic principle and echo of ISAR detection Signal receive process, which calculate, to be realized.ISAR analogue echoes parameter calculating mainly have running track, satellite platform, the figure of the earth, The factors such as earth rotation, target scattering also need to calculate the phase compensated, complete convolution algorithm, time delay, system thermal noise and estimate The functions such as meter, to obtain close to being really coherent video echo-signal of the target in actual moving process.
The classical signals system that current ISAR is used is linear FM signal, and the major way of echo reception is that solution line is mediated Reason.Solution line tune is the special pulse compress mode of one kind of linear FM signal.Not only operation is simple for it, but also can letter Change equipment.Currently, this method is widely used to ISAR.
It should be pointed out that although solution line mediates reason, the two or some difference identical with matched filtering basic principle 's.In order to correct understanding solution line tune method and its complex target broadband ISAR analogue echoes realization in bring problem, below Reason is mediated to ISAR reception solution line and makees some detailed description.
Assuming that transmitting signal is
WhereinfcFor center frequency, TpFor pulsewidth, γ is frequency modulation rate,
For the fast time, m is integer, T pulse repetition period, tm=mT is the slow time.Solution line tune is with delay Fixed, centre frequency and the identical LFM signal of frequency modulation rate are measured as reference signal and echo and makees difference frequency processing.
If reference distance is Rref, then reference signal are as follows:
T in formularefFor the pulsewidth of reference signal, it is usually larger than Echo width.Assuming that certain point target arrives radar Distance is Rt, the echo signal that radar receivesFor
The signal of line tune is solved as shown in Fig. 1 a to Fig. 1 c, if RΔ=Rt-Rref, then difference frequency exports are as follows:
Discussion is limited in (R in a pulse periodΔFor constant), then above formula is frequency and RΔDirectly proportional single-frequency arteries and veins Punching.Assume that observed range is [Rref-Δr/2,Rref+ Δ r/2], it depicts at range both sides of the edge in Fig. 1 a to Fig. 1 c Echo.Wherein ordinate is frequency, in Fig. 1 (a) in addition to reference signal, there is two remote, close echoes.Reference signal with return Wave conjugate multiplication, echo becomes simple signal, and its frequency is directly proportional to the range difference of echo and reference signal, and by Fig. 1 (b) Know fi=-2 γ RΔ/c.Fourier transform is made to the signal after solution line tune, can obtain corresponding to the sinc of each target in frequency domain The burst pulse of shape, pulse width 1/Tp, and pulse position and RΔIt is directly proportional.Obviously, the narrow pulse signal pair of frequency domain is transformed to The frequency domain resolution answered is 1/Tp, utilize the corresponding relationship f of frequency and distancei=-2 γ RΔ/ c, can obtain distance resolution is ρr= c/(2γTp)=c/ (2B), corresponding temporal resolution are 1/B, this is completely the same with the result of matched filtering.
Since the result of solution line tune pulse compression shows frequency domain, this method is otherwise known as " time-frequency conversion pulse pressure Contracting ".It, need to be multiplied by coefficient-c/ (2 γ) from frequency-domain transform to distance (relative to reference point).If Δ r is certain, after solving line tune Frequency range be [- γ Δ r/c, γ Δ r/c], signal maximum bandwidth be 2 γ Δ r/c=2B Δ r/ (cTp)=B Δ r/Rp, Middle RpFor TpCorresponding distance.As it can be seen that ratios delta r/RpSmaller, signal bandwidth is also smaller relative to original modulating bandwidth, This ratio is sometimes small to 1/tens or even more than one percent in spotlight mode SAR and ISAR, by taking ISAR as an example, flies The length of machine one kind target is generally less than 100 meters, and corresponding time width is several microseconds of zero point, and the broadband signal of big time width is general More than tens microseconds, so as to which signal band is reduced to only several megahertzs, to follow-up equipment (in especially from several hundred megahertzs Put and converted with A/D) it can simplify very much.It should be noted that the reduction of this frequency band is to be brought using the lengthening of time as cost, The signal in the short time is handled with the long time.
As it can be seen that the echo expression formula after solution line tune is more complicated, phase term there are three it.RΔWith slow time tmBecome Change, makes corresponding frequency fiChange with two phases other in echo expression formula, and it is this transformation and different scatterings Center is relevant, it is meant that the corresponding frequency f of different scattering centersiAnd other two phases are all different.Obviously , the phase change of Section 2 makes echo generate Doppler, this is the required normal information of imaging, and Section 3 is solution line tune institute Exclusive, referred to as video residual phase (RVP), it will lead to doppler information and unwanted change occur.
Independent by latter two phase term for solving line tune echo-signal writes out:
In a short time, it is assumed that RΔVariation be approximately linear (high-order term can be ignored), i.e. RΔ=RΔ0+ Vrtm, andBy RΔWithIt substitutes into, obtains
It can thus be concluded that Doppler
B in formulaΔ0=γ TΔ0(and TΔ0=2RΔ0/ c), i.e., target is R relative to the distance of reference pointΔ0When, after solving line tune The frequency of signal.In fact, the above results can be by that will solve the time-domain signal after line tune to the fast time (using the time of reference point as base It is quasi-) it is obtained as Fourier transform:
Above formula shows after solving line tune pulse compression, is 1/T in the narrow pulse width of frequency domainp, frequency displacement is -2 γ RΔ/ c, separately It is outer there are two and RΔRelated phase term (Doppler and RVP).
By Fig. 1 (b) as it can be seen that by solution line tune, chirp becomes single-frequency, and frequency and distance negative (when For reference point) it is directly proportional, this is people's needs.But it can also be seen that each pure-tone pulse misalignment in time, between them There is certain time shift (- RΔ/ c=fi/ γ), and time shift is directly proportional to the solution frequency of line tune.According to the property of Fourier transform, when The time shift in domain is equivalent to frequency-domain linear phase and moves, and which results in the generations of RVP, can be made by the waveform to Fig. 1 (b) Dispersion delay process enables delay and fiDirectly proportional (=fi/ γ), the echo of the different distance in Fig. 1 (b) can be corrected, and make It is perfectly aligned in time, to eliminate RVP.
In practical applications, line tune afterpulse misalignment in time is solved, in addition to generating RVP, after also resulting in pulse pressure Secondary lobe is higher.The frequency waveform that the time domain impulse of rectangle passes through Fourier transform is sinc function, and the secondary lobe near main lobe is phase When high, it is necessary to be weighted processing and be inhibited.Solution line mediates reason can only be in time domain weighting, when all pulses are equal in time When alignment, each pulse can unify to weight, and pulse pressure result secondary lobe is lower;And the pulse being staggered on the time shown in Fig. 1 (b), it can only To Tpr=Tp+2RΔr/ c makees unified time weight, suitable to intermediate signal weighting, and the signal at both ends is improper.Unless 2 Δ r/c < < Tp, otherwise the peak sidelobe ratio of pulse pressure output is extremely difficult to actual requirement.For suppressed sidelobes, it is necessary to try by Each signal plays point alignment, as shown in Fig. 1 (c), then makees unified weighting processing.Because the time delay of time-domain signal is equivalent in frequency Domain changes to Fig. 1 (c) multiplied by linear phase factor, from Fig. 1 (b), and the time delay of signal is directly proportional to difference frequency, i.e.,
F in formulaΔFor difference frequency, γ is the frequency modulation rate of LFM signal.Treatment process by Fig. 1 (b) to Fig. 1 (c) is as shown in Figure 2. The dotted line left side is that time delay adjusts process in figure, is weighting pulse pressure process on the right of dotted line.
After solving line tune pulse compression, distance samples rate are as follows:
According to multiple sampling thheorem, the maximum band width of signal is equal to sample frequency, then has
fs=γ Δ T=2 γ Δ r/c
Due to number of samples M and sampling time TprAnd sample rate fsRelationship be
M=Tprfs
Then maximum observed range is
Δ r=M δr
Work as Tpr=TpWhen, solution line tune Range compress post-sampling rate is equal to resolution ratio:
In conclusion if the echo for generating compound movement target is simulated according to classical ISAR solution line tune mode, in addition to must The target scattering characteristics of palpus calculate and simulation is outer, it is necessary to the treatment process of the complexity such as oblique, meaning are gone by echo solution line tune, envelope Taste reference signal generator, frequency mixer, RVP compensation phase generator and a series of matched filtering must be provided on hardware Device etc. must complete the signal processing of the complexity such as Fourier transformation, convolution algorithm on software.On the one hand, it considerably increases The hardware device and complexity of echo simulator;On the other hand the error of calculation processing is also increased.Seriously affect broadband The quality and benefits of ISAR target echo simulation.
Classical ISAR echo simulator hardware system composition is as shown in Figure 3.
Specifically from the point of view of the step of analogue echoes decomposes, there is following ask in the realization of classical simulation device or analogue technique Topic:
(1) complex target scattering unit is numerous, and the corresponding delay parameter of different scattering units is different, and scattering coefficient is not Together, corresponding solution line tune output frequency parameter is different, needs to calculate separately;
(2) traditional analog technology needs to generate wide-band LFM echo-signal and mating generation wide-band linearity tune first Frequency reference signal, when signal bandwidth is larger, operand is very big;
(3) traditional analog technology needs to complete wide-band LFM echo-signal and the conjugation of linear frequency modulation reference signal multiplies Method operation and be filtered with RVP compensation operation, further increase the complexity of simulation;
(4) since complex target scattering unit is numerous, and the analog parameter that different scattering units are related to is entirely different, therefore (1) all calculating in-(3) are required to implement to all scattering units, considerably increase the operand of simulation, reduce generation The speed of echo data, it is difficult to meet the needs of simulating in real time;
(5) as the demand simulated to high-resolution imaging increases severely, broadband, the simulation demand of ultra-broadband signal are inevitable increasingly Greatly, the calculation amount and simulated pressure that conventional echo simulation faces will be increasing.
Summary of the invention
Recall to that wave simulation is computationally intensive, and the analogue simulation time is long for classical ISAR solution line, simulator is complicated, simulation essence It is poor to spend, it is difficult to the technical issues of meeting requirement of real-time higher imaging simulation application demand, especially, classical ISAR echo Signal processing that signal time domain simulation method needs to be directed to by more complicated signal processing (echo solution line tune, Envelope goes tilting etc.) operand is bigger, and especially in the higher imaging emulation of requirement of real-time, traditional time domain is returned Wave generating algorithm is very time-consuming, and there are hardware systems complicated, simulation error technical problem rambunctious, the invention proposes A method of it quickly generates and is disposably unified in frequency domain using the technical solution with simulation ISAR echo data and realizes solution line Crucial calculating and the simulation processes such as tilting processing, residual video phase correction are gone in tune operation, envelope, pass through the frequency that can uniquely divide The thousands of target scattering center echo datas of information parallel computation substantially increase the speed that ISAR echo data is generated with simulated Degree, meets the higher analogue simulation demand of real-time, while enormously simplifying the hardware system of analogue echoes process and simulator Structure.
The present invention provides a kind of methods quickly generated with simulation ISAR echo data, comprising:
The scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter are read in and deposit to In register;
Register parameters are read parallel by parallel computing module, are completed different scattering center echo datas and are calculated, frequency Domain, which quickly generates, is completed the ISAR solution line tune echo data that envelope goes tilting processing and residual video phase correction process, and will Different scattering center echo datas are summed and are exported to RAM memory port.
Further, described by the scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter It reads in and deposits the step into register and include:
Step 1: the scattering center parameter of needed simulated target is read in and is deposited into register:
According to the scattering point range information R of complex target to be simulatedtAnd solution line tune reference signal reference distance information Rref, determine corresponding frequency parameter after scattering point solution line tune:
And corresponding echo delay:
Wherein, γ is the chirp rate of the linear FM signal of radar transmitting, and c is free space Electromagnetic Wave Propagation speed Degree, the i.e. light velocity.
Step 2: the time series parameters of needed simulated target are read in and are deposited into register:
According to sample frequency fsAnd the corresponding echo of target enrolls window initial time TsAnd end time Tf, determine Sampling time sequence
Step 3: the spectrum sequence parameter of needed simulated target is read in and is deposited into register:
Corresponding spectral sample sequence:
Wherein, total sampling number is Ns, value is identical as echo samples length of time series.
Further, register parameters are read by computing module parallel, different scattering center echo datas is completed and calculates Generate be completed envelope go it is tilting processing and residual video phase correction process ISAR solution line tune echo data the step of include:
According to scattering center echo scattering coefficient σt, radar pulse width Tp, during spectral sample sequence f, each component respectively scatter Corresponding frequency f after heart solution line tunet, radar carrier frequency fcAnd echo delay tdt, generate the corresponding frequency domain solved after line tune Echo-signal sdc, mathematic(al) representation are as follows:
sdctTpsinc[Tp(f-ft)]exp(-j2πfctdt)
Wherein, j indicates imaginary unit,The expression-form of sinc function sinc are as follows:
Further, different scattering center echo datas are summed and is exported to the step of RAM memory port and include:
It can be quickly completed in such a way that parallel computation is disposably accumulated in distributed thousands of scatterings of collection general objective Heart solution line tune echo data quickly generates:
It is expressed in the form of matrix operation are as follows:
By above-mentioned matrix sdcmIt sums by row, obtains the solution line tune echo data s that each component frequency domain generatesdcm, described will return Wave number is stored according to component separation is pressed;Wherein, σtmnFor the corresponding scattering coefficient of n-th scattering center of m-th of component, ftmnFor The corresponding solution line tune signal frequency of n-th of scattering center of m-th of component, tdtmnFor n-th of scattering center pair of m-th of component The echo delay answered.
The present invention also provides a kind of echo simulators, including quickly generating and simulating ISAR number of echoes described one kind According to method, be solidificated in programmable chip with software mode, realize that ISAR echo data quickly generating and simulates.
Further, programmable chip includes fpga chip.
The beneficial effects of the present invention are:
Technical solution provided by the invention is related to a kind of for broadband Inverse Synthetic Aperture Radar echo simulator Fast simulation With the new technique solution for generating complex target echo.Power-driven plane, surface vessel, space suitable for model-driven are defended The high-resolution Inverse Synthetic Aperture Radar analogue echoes and imaging simulation of the non-cooperation compound movement target of magnitude can effectively solve inverse conjunction The problems such as complicated, computationally intensive, computational efficiency is low at simulation system when the radar simulation complex target echo of aperture.The present invention proposes A kind of method quickly generated with simulation ISAR echo data, driven completely by mathematical model, can guarantee analog-quality In the case of, high efficiency, the original echo simulation for fast implementing complex target broadband Inverse Synthetic Aperture Radar, and the number of echoes generated Understand the envelope that line is recalled to wave and must be completed according to synchronously completing and go tilting processing, residual video phase correction process etc., it can be ensured that The high quality of compound movement target, the simulation of high efficiency inverse synthetic aperture radar imaging modeling and simulation.
Detailed description of the invention
Fig. 1 a is target low coverage, long distance echo and the relationship between frequency and time schematic diagram for solving line tune reference signal;Fig. 1 b is directly to solve Gained echo-signal relationship between frequency and time schematic diagram after line tune;Fig. 1 c is after solution line tune and to implement envelope and remove tilting processing and residual video Echo relationship between frequency and time schematic diagram after phasing;
Fig. 2, which is shown, recalls to the flow diagram that wave implementation envelope goes tilting processing to solution line;
Fig. 3 show classical ISAR echo simulator hardware system composition schematic diagram;
Fig. 4 show ISAR echo simulator hardware system composition schematic diagram;
Fig. 5 show simulation objectives model schematic;
Frequency domain echo data (vertical slice is one-dimensional range profile) schematic diagram that Fig. 6 is quickly generated;
Fig. 7 show corresponding time domain echo data schematic diagram;
Fig. 8 show the target ISAR image schematic diagram by imaging.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work It encloses.
Embodiment one
Classical ISAR solves line tune simulated radar echo, is substantially the method for a kind of time domain generation and analogue echo.Because of mesh Mark scattering center enormous amount, each scattering center time domain parameter is more, difference is big, and classical time-domain-simulation simulation needs node-by-node algorithm, Point-by-point superposition, operand are big;And analogue simulation has to pass through reference signal and generates, the multiplying of echo solution line tune, goes envelope oblique The signal processing of the complexity such as frequency-domain transform, the multiplication of frequency domain compensation function, frequency domain to time domain inverse transformation is set, operand is further Increase, the calculating link of simulation is more, and control errors are difficult, and the hsrdware requirements and complexity of simulator are high, it is difficult to meet the following high score Distinguish the demand of radar echo simulation.
Recall to that wave simulation is computationally intensive, and the analogue simulation time is long for classical ISAR solution line, simulator is complicated, simulation essence It is poor to spend, it is difficult to the problem of meeting requirement of real-time higher imaging simulation application demand, the embodiment of the present invention one propose it is a kind of fastly Fast-growing is at the method with simulation ISAR echo data, using new technical solution, is disposably unified in frequency domain and realizes solution line tune behaviour Work, envelope go crucial calculating and the simulation processes such as tilting processing, residual video phase correction, pass through the frequency information that can uniquely divide The echo data of the thousands of target scattering centers of parallel computation substantially increases the speed that ISAR echo data is generated with simulated Degree, meets the higher analogue simulation demand of real-time, while enormously simplifying the hardware system of analogue echoes process and simulator Structure.
The embodiment of the present invention one provides a kind of method quickly generated with simulation ISAR echo data, comprising:
The scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter are read in and deposit to In register;
Register parameters are read parallel by parallel computing module, are completed different scattering center echo datas and are calculated, from Frequency domain, which quickly generates, is completed the ISAR solution line tune echo data that envelope goes tilting processing and residual video phase correction process, and Different scattering center echo datas are summed and are exported to RAM memory port.
What the embodiment of the present invention one reached has the beneficial effect that:
(1) frequency domain quickly generates echo data, greatly improves operation efficiency;Overcome time domain simulation algorithm target scattering center The problems such as parameter is more, difference is big realize that thousands of target scattering center echo datas are parallel with the frequency information that can uniquely divide Simulation;
(2) envelope is completed while generating echo data and goes tilting processing, and when avoiding, frequency domain is mutually converted and classical calculation Method frequency domain residual video phase (RVP) compensation, reduces operand;
(3) the frequency domain echo data generated are exactly one-dimensional range profile;
(4) time solution line tune echo data is the inverse Fourier transform for generating frequency domain echo data;
(5) analogue echoes process is simplified, is referred back to without first generating wide-band LFM echo and wide-band LFM Wave, without the conjugation multiplying carried out between the two;Simulator hardware system structure is simplified accordingly, is not necessarily to broadband line Property frequency modulation echo-signal generator, reference signal generator, frequency mixer, RVP compensation phase generator and matched filter.
It is described that the scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter are read in and posted The step deposited into register includes:
Step 1: the scattering center parameter of needed simulated target is read in and is deposited into register:
According to the scattering point range information R of complex target to be simulatedtAnd solution line tune reference signal reference distance information Rref, determine corresponding frequency parameter after scattering point solution line tune:
And corresponding echo delay:
Wherein, γ is the chirp rate of the linear FM signal of radar transmitting, and c is free space Electromagnetic Wave Propagation speed Degree, the i.e. light velocity.
Step 2: the time series parameters of needed simulated target are read in and are deposited into register:
According to sample frequency fsAnd the corresponding echo of target enrolls window initial time TsAnd end time Tf, determine Sampling time sequence
Step 3: the spectrum sequence parameter of needed simulated target is read in and is deposited into register:
Corresponding spectral sample sequence:
Wherein, total sampling number is Ns, value is identical as echo samples length of time series.
Register parameters are read parallel by computing module, are completed different scattering center echo datas calculating generations and are completed Envelope goes the step of ISAR solution line tune echo data (one-dimensional range profile) of tilting processing and residual video phase correction process to wrap It includes:
According to scattering center echo scattering coefficient σt, pulse width Tp, spectral sample sequence f, each scattering center solution of each component Corresponding frequency f after line tunet, radar carrier frequency fcAnd echo delay tdt, generate the corresponding frequency domain echo solved after line tune Signal sdc, mathematic(al) representation are as follows:
sdctTpsinc[Tp(f-ft)]exp(-j2πfctdt)
Wherein, j indicates imaginary unit,The expression-form of sinc function sinc are as follows:
Solidification computing module is programmed with fpga chip, which completes echo data meter by reading register parameters It calculates and exports data to RAM memory port.
Different scattering center echo datas are summed and are exported to the step of RAM memory port and include:
In view of target scattering center is more, node-by-node algorithm is than relatively time-consuming, in such a way that parallel computation is disposably accumulated Quickly generating for the distributed collection thousands of a scattering center solution line tune echo datas of general objective can be quickly completed:
It is expressed in the form of matrix operation are as follows:
By above-mentioned matrix sdcmIt sums by row, obtains the solution line tune echo data s that each component frequency domain generatesdcm, described will return Wave number is stored according to component separation is pressed;Wherein, σtmnFor the corresponding scattering coefficient of n-th scattering center of m-th of component, ftmnFor The corresponding solution line tune signal frequency of n-th of scattering center of m-th of component, tdtmnFor n-th of scattering center pair of m-th of component The echo delay answered.
As shown in figure 4, the present invention also provides a kind of echo simulator, including quickly generating and simulating described one kind The method of ISAR echo data, is solidificated in programmable chip with software mode, realize ISAR echo data quickly generate with Simulation.Programmable chip includes FPGA chip.
Singualtion may be implemented in entire echo simulator, in a manner of software radio (such as using single programmable chip FPGA it) realizes.
Using the above method, solution line tune ISAR echo data can be quickly generated, and the echo data generated has been completed Envelope goes tilting processing, may be directly applied to imaging.The test result carried out on common PC machine shows based on this hair Bright ISAR solves line tune echo data frequency domain Fast Generation, and echo data generation process is quick, substantially increases emulation meter Calculate efficiency.Compared to conventional Time-domain simulation algorithm, calculating speed improves 30 to 50 times.
The frequency domain echo number of simulation objectives model, 1024 pulses quickly generated is set forth in Fig. 5, Fig. 6, Fig. 7, Fig. 8 According to (one-dimensional range profile), corresponding time domain echo data and its imaging results after imaging.
The beneficial effects of the present invention are:
Technical solution provided by the invention is related to a kind of for broadband Inverse Synthetic Aperture Radar echo simulator Fast simulation With the new technique solution for generating complex target echo.Power-driven plane, surface vessel, space suitable for model-driven are defended The high-resolution Inverse Synthetic Aperture Radar analogue echoes and imaging simulation of the non-cooperation compound movement target of magnitude can effectively solve inverse conjunction The problems such as complicated, computationally intensive, computational efficiency is low at simulation system when the radar simulation complex target echo of aperture.The present invention proposes A kind of method quickly generated with simulation ISAR echo data, driven completely by mathematical model, can guarantee analog-quality In the case of, high efficiency, the original echo simulation for fast implementing complex target broadband Inverse Synthetic Aperture Radar, and the number of echoes generated Understand the envelope that line is recalled to wave and must be completed according to synchronously completing and go tilting processing, residual video phase correction process etc., it can be ensured that The high quality of compound movement target, the simulation of high efficiency inverse synthetic aperture radar imaging modeling and simulation.
The serial number of the above embodiments of the invention is only for description, does not represent the advantages or disadvantages of the embodiments.
It should be noted that for the various method embodiments described above, for simple description, therefore, it is stated as a series of Combination of actions, but those skilled in the art should understand that, the present invention is not limited by the sequence of acts described because According to the present invention, some steps may be performed in other sequences or simultaneously.Secondly, those skilled in the art should also know It knows, the embodiments described in the specification are all preferred embodiments, and related actions and modules is not necessarily of the invention It is necessary.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain Lid is within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.

Claims (6)

1. a kind of method quickly generated with simulation ISAR echo data characterized by comprising
The scattering center parameter of needed simulated target, time series parameters, spectrum sequence parameter are read in and deposited to deposit In device;
Register parameters are read parallel by parallel computing module, are completed different scattering center echo datas and are calculated, from frequency domain It quickly generates and the ISAR solution line tune echo data that envelope goes tilting processing and residual video phase correction process is completed, and will not It sums and is exported to RAM memory port with scattering center echo data.
2. the method as described in claim 1, which is characterized in that it is described by the scattering center parameter of needed simulated target, when Between sequential parameter, spectrum sequence parameter is read in and the step deposited into register includes:
Step 1: the scattering center parameter of needed simulated target is read in and is deposited into register:
According to the scattering point range information R of complex target to be simulatedtAnd solution line tune reference signal reference distance information Rref, really Determine corresponding frequency parameter after scattering point solution line tune:
And corresponding echo delay:
Wherein, γ is the chirp rate of the linear FM signal of radar transmitting, and c is free space propagation velocity of electromagnetic wave, i.e., The light velocity.
Step 2: the echo samples time series parameters of needed simulated target are read in and are deposited into register:
According to sample frequency fsAnd the corresponding echo of target enrolls window initial time TsAnd end time Tf, determine sampling Time series
Step 3: the spectrum sequence parameter of needed simulated target is read in and is deposited into register:
Corresponding spectral sample sequence:
Wherein, total sampling number is Ns, value is identical as echo samples length of time series.
3. the method as described in claim 1, which is characterized in that read register parameters parallel by computing module, complete not With the calculating of scattering center echo data, generates and the ISAR that envelope goes tilting processing and residual video phase correction process is completed Solve line tune echo data the step of include:
According to the scattering coefficient σ of scattering centert, radar pulse width Tp, spectral sample sequence f, each scattering center solution line of each component Corresponding frequency f after tunet, radar carrier frequency fcAnd echo delay tdt, generate corresponding solution line frequency modulation domain echo-signal sdc, Its mathematic(al) representation are as follows:
sdctTpsinc[Tp(f-ft)]exp(-j2πfctdt)
Wherein, j indicates imaginary unit,The expression-form of sinc function sinc are as follows:
4. the method as described in claim 1, which is characterized in that different scattering center echo datas are summed and are exported to RAM The step of port memory includes:
The distributed collection thousands of a scattering center solutions of general objective can be quickly completed in such a way that parallel computation is disposably accumulated Line tune echo data quickly generates:
It is expressed in the form of matrix operation are as follows:
By above-mentioned matrix sdcmIt sums by row, obtains the solution line tune echo data s that each component frequency domain generatesdcm, by the number of echoes It is stored according to component separation is pressed;Wherein, σtmnFor the corresponding scattering coefficient of n-th of scattering center of m-th of component, ftmnIt is m-th The corresponding solution line tune signal frequency of n-th of scattering center of component, tdtmnN-th of scattering center for m-th of component is corresponding Echo delay.
5. a kind of echo simulator, which is characterized in that including quickly generating and simulating one kind described in one of claim 1-4 The method of ISAR echo data, is solidificated in programmable logic chip with software mode, realizes the fast fast-growing of ISAR echo data At with simulation.
6. method as claimed in claim 5, which is characterized in that programmable chip includes fpga chip.
CN201910245260.8A 2019-03-28 2019-03-28 A method of it quickly generates and simulation ISAR echo data Pending CN109901166A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910245260.8A CN109901166A (en) 2019-03-28 2019-03-28 A method of it quickly generates and simulation ISAR echo data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910245260.8A CN109901166A (en) 2019-03-28 2019-03-28 A method of it quickly generates and simulation ISAR echo data

Publications (1)

Publication Number Publication Date
CN109901166A true CN109901166A (en) 2019-06-18

Family

ID=66953138

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910245260.8A Pending CN109901166A (en) 2019-03-28 2019-03-28 A method of it quickly generates and simulation ISAR echo data

Country Status (1)

Country Link
CN (1) CN109901166A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112182062A (en) * 2020-09-15 2021-01-05 中国人民解放军63660部队 Multi-target radar networking measurement data matching and cataloguing method
CN112882031A (en) * 2021-01-14 2021-06-01 中山大学 ISAR imaging method, device and storage medium of discrete frequency coding waveform
CN113030963A (en) * 2021-02-04 2021-06-25 中国人民解放军陆军工程大学 Bistatic ISAR sparse high-resolution imaging method combining residual phase elimination
CN114545959A (en) * 2022-02-24 2022-05-27 中国人民解放军战略支援部队航天工程大学 Remote sensing satellite platform control based on flutter information and image correction method thereof

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4450447A (en) * 1981-07-17 1984-05-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Synthetic aperture radar target simulator
CN101295019A (en) * 2008-05-23 2008-10-29 北京航空航天大学 Simulation method of linear frequency modulation continuous wave synthetic aperture radar video signal
CN101592733A (en) * 2009-07-01 2009-12-02 电子科技大学 A kind of parallel real-time imaging processing method for inverse synthetic aperture radar
CN101710177A (en) * 2009-12-22 2010-05-19 电子科技大学 Multi-target imaging method for inverse synthetic aperture radar
CN102401898A (en) * 2011-08-25 2012-04-04 北京理工大学 Quantified simulation method for forest remote sensing data of synthetic aperture radar
CN103176170A (en) * 2013-02-06 2013-06-26 中国科学院电子学研究所 SAR (synthetic aperture radar) echo simulating method based on GPU (graphics processing unit) parallel computing
CN104122551A (en) * 2014-07-12 2014-10-29 西安电子科技大学 Two-dimensional unitary ESPRIT based ISAR horizontal calibration method
CN104965197B (en) * 2015-06-30 2017-12-12 南京理工大学 Radar echo signal simulator based on FPGA
CN107765226A (en) * 2017-09-18 2018-03-06 北京空间飞行器总体设计部 A kind of SAR satellite radars analogue echoes method, system and medium
CN108459308A (en) * 2018-03-29 2018-08-28 北京环境特性研究所 A kind of analogue echoes method and device based on time-varying RCS data
CN108828541A (en) * 2018-08-01 2018-11-16 北京润科通用技术有限公司 A kind of radar echo simulation system
CN109343060A (en) * 2018-12-07 2019-02-15 电子科技大学 ISAR imaging method and system based on deep learning time frequency analysis

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4450447A (en) * 1981-07-17 1984-05-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Synthetic aperture radar target simulator
CN101295019A (en) * 2008-05-23 2008-10-29 北京航空航天大学 Simulation method of linear frequency modulation continuous wave synthetic aperture radar video signal
CN101592733A (en) * 2009-07-01 2009-12-02 电子科技大学 A kind of parallel real-time imaging processing method for inverse synthetic aperture radar
CN101710177A (en) * 2009-12-22 2010-05-19 电子科技大学 Multi-target imaging method for inverse synthetic aperture radar
CN102401898A (en) * 2011-08-25 2012-04-04 北京理工大学 Quantified simulation method for forest remote sensing data of synthetic aperture radar
CN103176170A (en) * 2013-02-06 2013-06-26 中国科学院电子学研究所 SAR (synthetic aperture radar) echo simulating method based on GPU (graphics processing unit) parallel computing
CN104122551A (en) * 2014-07-12 2014-10-29 西安电子科技大学 Two-dimensional unitary ESPRIT based ISAR horizontal calibration method
CN104965197B (en) * 2015-06-30 2017-12-12 南京理工大学 Radar echo signal simulator based on FPGA
CN107765226A (en) * 2017-09-18 2018-03-06 北京空间飞行器总体设计部 A kind of SAR satellite radars analogue echoes method, system and medium
CN108459308A (en) * 2018-03-29 2018-08-28 北京环境特性研究所 A kind of analogue echoes method and device based on time-varying RCS data
CN108828541A (en) * 2018-08-01 2018-11-16 北京润科通用技术有限公司 A kind of radar echo simulation system
CN109343060A (en) * 2018-12-07 2019-02-15 电子科技大学 ISAR imaging method and system based on deep learning time frequency analysis

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
XU Y , ET AL: ". Broadband electromagnetic scattering echo generation and its ISAR imaging simulation", 《 INTERNATIONAL CONFERENCE ON MICROWAVE & MILLIMETER WAVE TECHNOLOGY》 *
李卫星: "宽带雷达数据采集与实时成像处理技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
杨义先等编著: "《2006北京地区高校研究生学术交流会 通信与信息技术会议论文集 上》", 31 December 2006, 北京邮电大学出版社 *
雷杰等: "一种基于钟摆模型的舰船目标成像方法", 《电子与信息学报》 *
黄琳琳等: "动态目标宽带雷达回波频域模拟及成像仿真", 《航空兵器》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112182062A (en) * 2020-09-15 2021-01-05 中国人民解放军63660部队 Multi-target radar networking measurement data matching and cataloguing method
CN112182062B (en) * 2020-09-15 2022-09-13 中国人民解放军63660部队 Multi-target radar networking measurement data matching and cataloguing method
CN112882031A (en) * 2021-01-14 2021-06-01 中山大学 ISAR imaging method, device and storage medium of discrete frequency coding waveform
CN112882031B (en) * 2021-01-14 2023-09-05 中山大学 ISAR imaging method, device and storage medium of discrete frequency coding waveform
CN113030963A (en) * 2021-02-04 2021-06-25 中国人民解放军陆军工程大学 Bistatic ISAR sparse high-resolution imaging method combining residual phase elimination
CN114545959A (en) * 2022-02-24 2022-05-27 中国人民解放军战略支援部队航天工程大学 Remote sensing satellite platform control based on flutter information and image correction method thereof
CN114545959B (en) * 2022-02-24 2024-01-12 中国人民解放军战略支援部队航天工程大学 Remote sensing satellite platform control and image correction method based on flutter information

Similar Documents

Publication Publication Date Title
CN109901166A (en) A method of it quickly generates and simulation ISAR echo data
CN105259552B (en) A kind of synthetic aperture radar image-forming method and apparatus based on NLFM signal
Xin et al. Signal processing for digital beamforming FMCW SAR
EP2650695A1 (en) Imaging method for synthetic aperture radar in high squint mode
Saeedi et al. Synthetic aperture radar imaging using nonlinear frequency modulation signal
Cimmino et al. Efficient spotlight SAR raw signal simulation of extended scenes
Belcher et al. Theory and simulation of ionospheric effects on synthetic aperture radar
Vu et al. Fast time-domain algorithms for UWB bistatic SAR processing
CN109270528B (en) One-station fixed type double-station SAR imaging method based on full-analytic distance model
CN111505639A (en) Synthetic aperture radar wide-width sparse imaging method based on variable repetition frequency sampling mode
CN105445704A (en) Radar moving object inhibition method in SAR image
Spreng et al. UWB near-field MIMO radar: Calibration, measurements and image reconstruction
CN111190151A (en) System parameter design and emission power optimization method of multi-mode small satellite SAR in scanning mode
Mohr et al. Design and generation of stochastically defined, pulsed FM noise waveforms
US9134415B2 (en) Wideband waveform synthesis using frequency jump burst-type waveforms
CN109143236B (en) Bistatic bunching SAR large-scene imaging method suitable for complex flight trajectory
Wang et al. Squint spotlight SAR raw signal simulation in the frequency domain using optical principles
CN113484859A (en) Two-dimensional super-resolution radar imaging method based on fusion technology
Chen et al. Forward looking imaging of airborne multichannel radar based on modified iaa
Liu et al. Efficient and precise frequency-modulated continuous wave synthetic aperture radar raw signal simulation approach for extended scenes
Kim et al. SAR raw data simulation for multiple‐input multiple‐output video synthetic aperture radar using beat frequency division frequency modulated continuous wave
Ahmed et al. Enhanced azimuth resolution for spaceborne interrupted FMCW sar through spectral analysis
Grove MIMO radar Systems and Algorithms-Imperfections and Calibration
Rogers et al. The synthetic aperture radar trans-ionospheric radio propagation simulator (SAR-TIRPS)
Bączyk et al. Moving target imaging in multistatic passive radar

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190618

RJ01 Rejection of invention patent application after publication