CN107192992A - Method, calibrator, radar installations and the system of radar calibration - Google Patents

Method, calibrator, radar installations and the system of radar calibration Download PDF

Info

Publication number
CN107192992A
CN107192992A CN201710605439.0A CN201710605439A CN107192992A CN 107192992 A CN107192992 A CN 107192992A CN 201710605439 A CN201710605439 A CN 201710605439A CN 107192992 A CN107192992 A CN 107192992A
Authority
CN
China
Prior art keywords
mrow
radar
mfrac
signal
msub
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
Application number
CN201710605439.0A
Other languages
Chinese (zh)
Other versions
CN107192992B (en
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.)
Beijing Ju Net Technology Co., Ltd.
Original Assignee
Beijing Cape Interspace Technology Ltd
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 Beijing Cape Interspace Technology Ltd filed Critical Beijing Cape Interspace Technology Ltd
Publication of CN107192992A publication Critical patent/CN107192992A/en
Application granted granted Critical
Publication of CN107192992B publication Critical patent/CN107192992B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses the method for radar calibration, calibrator, radar installations and system.Wherein, radar calibration device, including:Reception antenna, the pulse signal Srx suitable for receiving radar;Demodulating unit, suitable for extracting intermediate-freuqncy signal Sdc from Srx;Signal processing unit, suitable for generating frequency shift signal and delay time signal based on intermediate-freuqncy signal Sdc;Transmitter unit, suitable for being mixed and amplifieroperation to frequency shift signal and delay time signal, to generate signal to be launched;Transmitting antenna, suitable for launching the time delay and frequency shift signal, so that radar receives the signal.

Description

Method, calibrator, radar installations and the system of radar calibration
Technical field
The present invention relates to method, calibrator, radar installations and the system of remote sensing fields, more particularly to radar calibration.
Background technology
In remote sensing using microwave it is most basic the reason for be that they can provide different information.Composed using electromagnetic wave (EM) Microwave region, we can obtain the information complementary with optical remote sensing --- generally speaking, influence target and microwave interactive As a result physical parameter and its characteristic is different from other electromagnetic waves.In addition, microwave also has many good qualities:Some wavelength The top layer of cloud layer even dry ground, sand can be penetrated (penetration depth is up to several meters sometimes).And except can be distant from passive microwave Thermal radiation information is obtained in sense image, meanwhile, user can also obtain aggressive mode image by providing signal source.It is this to survey Amount can be obtained at any time, it is not necessary to by context sources (such as sun).In atmospheric remote sensing, microwave is relative to its all band (such as infrared) also has other advantages, you can pass through the selection of microwave wavelength so that other particles, the effect of particulate around observation medium It should can be ignored.
Radar to be such as imaged, navigate, remote sensing and the miscellaneous application of global location have long been used for military affairs With non-military purpose.Synthetic aperture radar (" SAR ") is that mainly radar system (" SLAR ") is checked in relevant aviation or satellite side, should Radar system carrys out electronics using the flight path of mobile platform (that is, delivery vehicle, such as aircraft or satellite) disposed thereon SAR The very big antenna of simulation or aperture, and generation high-resolution remote sensing image.At present, existing radar system is square in image calibration Face need to be improved.
Therefore, the present invention proposes a kind of new radar calibration scheme.
The content of the invention
Therefore, the present invention provides a kind of radar calibration scheme, effectively solve it is above-mentioned at least one problem.
According to an aspect of the present invention there is provided a kind of radar calibration device, including:Reception antenna, suitable for receiving radar Pulse signal Srx;Demodulating unit, suitable for extracting intermediate-freuqncy signal Sdc from Srx;Signal processing unit, suitable for based on intermediate-freuqncy signal Sdc generates frequency shift signal and delay time signal;Transmitter unit, suitable for behaviour is mixed and amplified to frequency shift signal and delay time signal Make, to generate transmitting time delay to be launched and frequency shift signal;Transmitting antenna, suitable for transmitting time delay and frequency shift signal, so as to radar Receive the signal.
Alternatively, in the calibrator according to the present invention, demodulating unit is suitable to filter carrier wave from Srx according to following manner Signal, radar and the intermediate-freuqncy signal Sdc of calibrator distance are included to extract:
Alternatively, in the calibrator according to the present invention, the signal processing unit includes:Frequency shift unit, suitable for generation The frequency shift signal of the bright spot of orientation in the radar image of correspondence radar generation;Time delay elements, suitable for generation correspondence radar The delay time signal of the bright spot of the radar image middle-range descriscent of generation.
According to a further aspect of the invention there is provided a kind of radar installations, including:Reception antenna, thunder is come from suitable for receiving Up to the time delay and frequency shift signal of calibrator;Wave filter, suitable for eliminating time delay and frequency based on Wigner-Ville distribution and Hankel Time frequency signal interferes with each other in shifting signal;Geometry correction unit, suitable for generating radar image according to time delay and frequency shift signal, and According to the coordinate pair of radar image lieutenant colonel on schedule, the radar image carries out geometric correction;Radiant correction unit, suitable for according to time delay With the energy balane radiant correction coefficient of frequency shift signal, and in radar image pixel carry out radiant correction.
Alternatively, in the radar installations according to the present invention, the time delay and frequency shift signal received is suitable to pass through following sides Formula is represented:
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one when being moved on radar installations under speed u The slow time function.
Alternatively, in the radar installations according to the present invention, geometry correction unit is suitable to according to following manner to radar map As carrying out geometric correction:The default spacing conduct in radar image bright spot is produced according to the frequency shift unit and the time delay elements Precision reference control point carries out geometric correction, wherein default spacing is by setting the frequency shift amount and delay volume of radar calibration device come really It is fixed.
Alternatively, according to the present invention radar installations in, radiant correction unit be suitable to according to following manner come according to when Prolong the energy balane radiant correction COEFFICIENT K with frequency shift signal:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazRespectively the oblique distance and orientation of image are divided Resolution, σTFor the radar cross section of object, α is radar incidence angle.
Alternatively, in the radar installations according to the present invention, wave filter is suitable to according to following manner come based on Wigner- Ville is distributed and Hankel eliminates interfering for time frequency signal:
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported, Produce the high-resolution simultaneously in time delay and frequency displacement;
Multiple signals are disassembled using Hankel conversions, each single signal is other separated, then input to WVD In, to eliminate interfering for time frequency signal, and the linear FM characteristic in stick signal.
According to a further aspect of the invention there is provided a kind of radar system, including according to the present invention radar calibration device and Radar installations.
According to a further aspect of the invention there is provided a kind of radar calibration method, suitable for being performed in radar calibration device.Should Method includes:Receive the pulse signal Srx of radar;Intermediate-freuqncy signal Sdc is extracted from Srx;Based on intermediate-freuqncy signal Sdc, generation frequency Shifting signal and delay time signal;Frequency shift signal and delay time signal are mixed and amplifieroperation, with generate time delay to be launched and Frequency shift signal;Launch time delay and frequency shift signal, so that radar receives the time delay and frequency shift signal.
According to a further aspect of the invention there is provided a kind of method that radar installations is performed, including:Receive and come from radar school The time delay and frequency shift signal of quasi- device;Time frequency signal in time delay and frequency shift signal is eliminated based on Wigner-Ville distribution and Hankel Interfere with each other;According to time delay and frequency shift signal generation radar image, and according to the coordinate pair of the radar image lieutenant colonel on schedule thunder Geometric correction is carried out up to image;Radiant correction unit, suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal, And radiant correction is carried out to pixel in radar image.
, can be along range direction to received radar by radar calibration device according to the radar calibration technology of the present invention Pulse signal carries out time delay operations, and carries out frequency displacement operation to radar pulse signal along azimuth direction.So, radar installations root The original image generated according to time delay and frequency shift signal can include two-dimension netted calibration point.Each calibration point is to be based on frequency displacement Operation and time delay operations and the result that occurs.On this basis, radar installations can improve radiation according to netted calibration point Correction and the precision of geometric correction.
Brief description of the drawings
In order to realize above-mentioned and related purpose, some illustrative sides are described herein in conjunction with following description and accompanying drawing Face, these aspects indicate the various modes of principles disclosed herein that can put into practice, and all aspects and its equivalent aspect It is intended to fall under in the range of theme claimed.The following detailed description by being read in conjunction with the figure, the disclosure it is above-mentioned And other purposes, feature and advantage will be apparent.Throughout the disclosure, identical reference generally refers to identical Part or element.
Fig. 1 shows the application scenarios schematic diagram of radar system 100 according to some embodiments of the invention;
Fig. 2 shows the schematic diagram of radar calibration device 200 according to some embodiments of the invention;
Fig. 3 shows the schematic diagram of demodulating unit 220 according to an embodiment of the invention;
Fig. 4 shows the schematic diagram of radar installations 400 according to some embodiments of the invention;And
Fig. 5 shows the schematic diagram of radar image according to an embodiment of the invention.
Embodiment
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in accompanying drawing Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here Limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure Complete conveys to those skilled in the art.
Fig. 1 shows the application scenarios schematic diagram of radar system 100 according to some embodiments of the invention.
As shown in figure 1, radar installations 110 (rader) is deployed in such as satellite, aircraft on aircraft (not shown).Allusion quotation Type, radar installations 110 is, for example, synthetic aperture radar (SAR).Based on the three-dimensional system of coordinate (x, y, z) shown in Fig. 1, radar Device is h relative to the flying height on ground.Radar installations is moved with speed u along x-axis forward direction, and with incident angle α and along distance Launch pulse signal in direction.Here pulse signal is, for example, chirp pulse signal.Chirp (Linear Frequency Modulation, LFM) it is the signal that a frequency changes over time, it is referred to as chirping of birds communication signal (chirp signal).The frequency of signal can rise (up chirp) over time during modulation or decline (down Chirp change) and linearly.In one embodiment, the mode of signal frequency change is as follows:
ft=f0+aTt
f0Represent initial frequency, aTRepresent the rate of change (chirp rate) of frequency.
According to one embodiment of the invention, the mathematic(al) representation of LFM pulse signals is:
S (t)=exp (j π art2)
The scanning area (swath) of radar installations 110 includes location point P.P position point is suitable to thunder of the arrangement according to the present invention Up to calibrator 120 (TF-IRC).Here, the position coordinates (longitude and latitude) of radar calibration device 120 is determined.Radar calibration device 120 can To carry out time delay operations to received radar pulse signal along range direction, and along azimuth direction to radar pulse signal Carry out frequency displacement operation.So, the original image that radar installations 110 is generated according to time delay and frequency shift signal can include two-dimensional mesh The calibration point of shape.Each calibration point is the result occurred based on frequency displacement operation and time delay operations.On this basis, radar installations 110 can improve the precision of radiant correction and geometric correction according to netted calibration point.With reference to Fig. 2 to according to the present invention Radar calibration device it is illustrative.
Fig. 2 shows the schematic diagram of radar calibration device 200 according to some embodiments of the invention.
As shown in Fig. 2 radar calibration device 200 includes reception antenna 210, demodulating unit 220, signal processing unit 230, hair Penetrate unit 240 and transmitting antenna 250.Here, reception antenna 210 both can be discrete or multiplexing, and the present invention is right This does not do excessive limitation.After reception antenna 210 receives radar pulse signal Srx, demodulating unit 220 is suitable to carry from Srx Take intermediate-freuqncy signal Sdc.Signal processing unit 230 can generate two paths of signals based on Sdc.Wherein, all the way to pass through time delay operations Signal St, another road be by frequency displacement operate signal Sp.According to one embodiment of the invention, the quilt of signal processing unit 230 It is configured to include frequency shift unit and time delay elements.Frequency shift unit is suitable to orientation in the radar image of generation correspondence radar generation Bright spot frequency shift signal.Time delay elements are suitable to the time delay of the bright spot of the radar image middle-range descriscent of generation correspondence radar generation Signal.Here, signal processing unit 230 can be application specific integrated circuit (ASIC), can also be realized based on fpga chip, this In repeat no more.For example, by the FPGA of programmable, radar system 100 is able to produce the school of different frequency displacement and time shift distances On schedule.Frequency displacement and time shift distance due to generation can be adjusted voluntarily, are very helpful for the geometric correction on image.Simultaneously Radar calibration device 200 is difficult the characteristic of decay based on digital data transmission so that the power output of frequency displacement and time shift point is more steady It is fixed.
Transmitter unit 240 is suitable to St and Sp being mixed and amplifieroperation, is then launched by transmitting antenna 250 and returned To the signal of radar installations.
Fig. 3 shows the schematic diagram of demodulating unit 220 according to an embodiment of the invention.Demodulating unit 220 is following Formula carries out signal transacting.
Fig. 4 shows the schematic diagram of radar installations 400 according to some embodiments of the invention.As shown in figure 4, radar installations 400 include reception antenna 410, wave filter 420, geometry correction unit 430 and radiant correction unit 440.Reception antenna 410, is fitted In receiving time delay from radar calibration device and frequency shift signal.Wave filter 420 be suitable to based on Wigner-Ville distribution and Time frequency signal interferes with each other in Hankel elimination time delays and frequency shift signal.Geometry correction unit 430 is suitable to according to time delay and frequency Shifting signal generates radar image, and the radar image carries out geometric correction according to the coordinate pair of radar image lieutenant colonel on schedule.Radiation Unit 440 is corrected suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal, and to pixel in radar image Carry out radiant correction.
Fig. 5 shows radar image schematic diagram according to an embodiment of the invention.As shown in figure 5, figure include it is multiple Bright spot.Each bright spot represents a calibration point (calibrator for example as shown in Figure 2 carries out frequency displacement and time delay operations and produced), i.e., Determine the radar cross section (Radar cross section, RCS) in place.So, although Fig. 5 radar image has not been Whole earth's surface image (in other words, radar can not quick obtaining to whole radar image coordinate), radar installations can pass through school (its gray value and coordinate are determined) to carry out radiant correction and geometric correction to whole radar image on schedule.Wherein, geometric correction Operation be mainly based upon the coordinate value of calibration point, according to the signal delay time of calibration point, longitude and latitude is done to whole radar image Registration process is spent to realize.
It is illustrative to geometry and radiant correction operation with reference to Fig. 1,4 and 5.
Firstly the need of explanation, the time delay received and frequency shift signal of radar installations are as follows:
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one and moved down on radar installations 100 in speed u The function of slow time when dynamic.
The phase of wherein orientation is exactlyDistance to phase be The side of being Position chirp rate.
In an embodiment in accordance with the invention, geometry correction unit 430 is entered with specific reference to following manner to radar image Row geometric correction:The default spacing in radar image (Fig. 5) bright spot is produced as precision ginseng according to frequency shift unit and time delay elements Examine control point and carry out geometric correction.Wherein, spacing is preset to determine by setting the frequency shift amount and delay volume of radar calibration device.
In addition, can be derived from the backward energy mathematical formulae that radar installations 400 is received by above-mentioned radar equation:
Wherein,
Pr:The echo power that SAR system is received
Pt:The power of SAR system transmitting
G:The gain of antenna
λ:The wavelength of radar wave
G:The distance of antenna and object
Ls:The loss of air in transmitting procedure
σ:The radar cross section (RCS) of object
For the region (Homogeneous region) of a piece of homogeneous atural object, RCS can be write as:
α:Incidence angle of the incidence wave to object
ρrg:Oblique distance directional resolution
ρaz:Azimuth direction resolution ratio
σ0:Backscattering coefficient
Thus, it is possible to which backward energy equation is as follows:
The method for calculating the backscattering coefficient of homogeneous area is as follows.Known radar section is put in the region of a homogeneous Product is σTObject, and the image energy of object be PrT, and homogeneous area is AHHomogeneous region, received energy For PrH.Then the backscattering coefficient of this homogeneous area can be obtained so
For homogeneous area AH
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images.Homogeneous region and known radar Sectional area σTThe distance of object to radar antenna have to close or equal, can just obtain close, correct correction coefficient K, if oblique distance gap is excessive on image, may set up multiple known RCS values along oblique distance direction to calculate object, finally Obtain the correction coefficient of multiple different oblique distances.
By above-mentioned formula, we can obtain radar image backscattering conversion coefficient:
PrT:The reflection power of object
ρrg:Oblique distance directional resolution
ρaz:Azimuth direction resolution ratio
σT:The radar cross section of object
α:Radar incidence angle of the incidence wave to object
Calculate K value modes as follows:The image of Digital Synthetic Aperture Radar includes the imaging contexts of homogeneous area and object. One piece of area is first cut out from image for AHHomogenous area.N is contained in the homogenous areaHIndividual pixel, and ensure this homogeneous Without other objects in region.Now, the reflection general power expression of the homogenous area is:
αij:Quantized value (digital number) on image
Power density expression is:
Then cut out from image and contain N in another piece of region containing object and background, this regionpIndividual pixel, And the reflection general power expression for containing object and background is:
It can finally calculate and comprise only the reflection power of object and be:
The general power that can calculate object by above formula is PrT, add ρrg、ρazIt can be learnt in image, and object Radar cross section σTAnd incident angle α is, it is known that therefore correction coefficient K can be calculated.
In an embodiment in accordance with the invention, radiant correction unit 440 be suitable to according to following manner come according to time delay and The energy balane radiant correction COEFFICIENT K of frequency shift signal:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazRespectively the oblique distance and orientation of image are divided Resolution, σTFor the radar cross section of object, α is radar incidence angle.
Wigner-Ville (WVD) distributions can all reach high-resolution over time and frequency.But for synthesis hole For the radar of footpath, time delay and frequency shift signal be in the range of radar beam from return relevant of each individual goal and.This will Interference is formed between multiple signals.Fortunately, Fourier's Bezier expansion can be decomposed preferably compared to WVD analyses These multiple signals.
Radar installations 400 also includes being used for eliminating based on the Novel Filter that Wigner-Ville distribution and Hankel are converted The interference of time frequency signal.
The operation principle of wave filter is described as follows:
The power spectrum P (ω) for describing a time domain communication signal s (t) is the auto-correlation function of this communication signal in itself Fourier transform, its expression is as follows:
Wherein R (τ) is time domain communication signal s (t) and the auto-correlation function of time independence, is expressed as follows:
However, this power spectrum only shows the composition of communication signal frequency in this time range, and without Faxian The composition for showing frequency how is distributed in this time range, that is, frequency is how to change over time.Therefore Redefine its expression with auto-correlation function R (t, the τ) of time correlation as follows:
The R (τ) of (1) formula is replaced with R (t, τ) formula of (3) formula, following formula is can obtain
After the replacement of auto-correlation function, power spectrum has become the function with time correlation, then by arranging, can Define WVD as follows:
Then the mathematic(al) representation of LFM communication signals is:
S (t)=exp (j π art2) (6)
Then WVD is:
And the time delay and frequency shift signal for adding Doppler frequency can be write as following formula:
Wherein fDopplerFor Doppler frequency
(9) formula of utilization makees correlation output with (10) formula and derives such as following formula:
Wherein TaFor generated time
By (11), formula is known that Scorr(toffset) maximum must be in toffsetShift time (offset time) isWhen have maximum, therefore it is recognised that Doppler frequency is relative to the displacement of transfer point Formula is shown below:
In an embodiment in accordance with the invention, when wave filter 420 is based on Wigner-Ville distribution and Hankel eliminations The course of work interfered of frequency signal is as follows.Utilize the corresponding Doppler frequency of maximum of the WVD correlation functions exported Relative to the displacement of transfer point, the high-resolution simultaneously in time delay and frequency displacement is produced.Disassembled using Hankel conversions Multiple signals, each single signal is other separated, then input into WVD, to eliminate interfering for time frequency signal, And the linear FM characteristic in stick signal.
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported, Very high resolution ratio can be produced in time delay and frequency displacement simultaneously.Because calibrator signal is multiple signal, i.e., multiple time delays Point and many frequency displacement points, can be interfered in WVD times and frequency.The interference of time frequency signal can be disassembled using Hankel conversions Multiple signals, it is each single signal is other separated, then input into WVD, you can the problem of to solve distracter, and Linear FM characteristic in signal can also completely retain.
In addition, the radar cross section (RCS) for providing maximum according to the amplifying power of the radar calibration device of the present invention can be full Sufficient radar system can receive the RCS of maximum.It is therefore assumed that the dynamic range that radar system can detect maximum is backscattering system Number ρ0 max, radar system resolution ratio is A, then the maximum radar cross section that radar system can be detected is:
σDR_max0 max×A (13)
The radar cross section of active radar calibration device is represented by:
When the sectional area of active radar calibration device is equal to the maximum radar cross section that may detect that, we can assess The gain amplifier G of active radar calibration deviceampFor:
Attenuation range:
The attenuation range of the program-controlled attenuator of design system, is according to the dynamic range that can be detected in radar system.
Peak power output:
When designing active radar calibration device, it is necessary to its peak power output is considered, to avoid intraware Beyond working range, so first understanding the communication signal intensity that active radar calibration device receives maximum first, relation is as follows:
P in above formulatIt is the transmission power of SAR system, G is the gain of SAR system transmitting antenna, RnearBe it is airborne to ground The distance of table low coverage point.Have receive communication signal size, we combine 1.25 and 1.26 formulas can learn we maximum Power output is:
Therefore, the power that the component of radar calibration device of the invention can bear, will be more than peak power output, this Sample just can ensure that component operational is normal.
Time delay
In SAR width Wide, temporal delay also has limitation.If delay exceedes greatly very much covering scope, that is in shadow Do not show not come out as upper, so time difference τ of the scope of our time delays no more than low coverage point to remote point yet1, such as following formula It is shown:
Radar system echo communication signal emulator:
Radar system echo communication signal emulator is with the maximum difference of active radar calibration device before, emulator When doing echo simulation to radar system indoors, not via antenna but be directly connected to inside loop because radar system Communication signal out is very high for rate, so communication signal must first pass through high power attenuator, decaying to can be with emulator work The scope of work.Pad value is as follows:
Attenuation range be also the dynamic range according to radar system as foundation, the part of time delay because being echo Emulator, so being the echo communication signal between emulating communication signal from low coverage point to remote point, so time delay τ2Such as Under:
In the specification that this place is provided, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention Example can be put into practice in the case of these no details.In some instances, known method, knot is not been shown in detail Structure and technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the disclosure and help to understand one or more of each inventive aspect, exist Above in the description of the exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention:It is i.e. required to protect The application claims of shield are than the feature more features that is expressly recited in each claim.More precisely, as following As claims reflect, inventive aspect is all features less than single embodiment disclosed above.Therefore, abide by Thus the claims for following embodiment are expressly incorporated in the embodiment, wherein each claim is in itself It is used as the separate embodiments of the present invention.
Those skilled in the art should be understood the module or unit or group of the equipment in example disclosed herein Part can be arranged in equipment as depicted in this embodiment, or alternatively can be positioned at and the equipment in the example In different one or more equipment.Module in aforementioned exemplary can be combined as a module or be segmented into addition multiple Submodule.
Those skilled in the art, which are appreciated that, to be carried out adaptively to the module in the equipment in embodiment Change and they are arranged in one or more equipment different from the embodiment.Can be the module or list in embodiment Member or component be combined into a module or unit or component, and can be divided into addition multiple submodule or subelement or Sub-component.In addition at least some in such feature and/or process or unit exclude each other, it can use any Combination is disclosed to all features disclosed in this specification (including adjoint claim, summary and accompanying drawing) and so to appoint Where all processes or unit of method or equipment are combined.Unless expressly stated otherwise, this specification (including adjoint power Profit is required, summary and accompanying drawing) disclosed in each feature can or similar purpose identical, equivalent by offer alternative features come generation Replace.
Although in addition, it will be appreciated by those of skill in the art that some embodiments described herein include other embodiments In included some features rather than further feature, but the combination of the feature of be the same as Example does not mean in of the invention Within the scope of and form different embodiments.For example, in the following claims, times of embodiment claimed One of meaning mode can be used in any combination.
In addition, be described as herein can be by the processor of computer system or by performing for some in the embodiment Method or the combination of method element that other devices of the function are implemented.Therefore, with for implementing methods described or method The processor of the necessary instruction of element forms the device for implementing this method or method element.In addition, device embodiment Element described in this is the example of following device:The device is used to implement as in order to performed by implementing the element of the purpose of the invention Function.
As used in this, unless specifically stated so, come using ordinal number " first ", " second ", " the 3rd " etc. Description plain objects are merely representative of the different instances for being related to similar object, and are not intended to imply that the object being so described must Must have the time it is upper, spatially, in terms of sequence or given order in any other manner.
Although describing the present invention according to the embodiment of limited quantity, above description, the art are benefited from It is interior it is clear for the skilled person that in the scope of the present invention thus described, it can be envisaged that other embodiments.Additionally, it should be noted that The language that is used in this specification primarily to readable and teaching purpose and select, rather than in order to explain or limit Determine subject of the present invention and select.Therefore, in the case of without departing from the scope and spirit of the appended claims, for this Many modifications and changes will be apparent from for the those of ordinary skill of technical field.For the scope of the present invention, to this The done disclosure of invention is illustrative and not restrictive, and it is intended that the scope of the present invention be defined by the claims appended hereto.

Claims (17)

1. a kind of radar calibration device, including:
Reception antenna, the pulse signal Srx suitable for receiving radar;
Demodulating unit, suitable for extracting intermediate-freuqncy signal Sdc from Srx;
Signal processing unit, suitable for generating frequency shift signal and delay time signal based on intermediate-freuqncy signal Sdc;
Transmitter unit, suitable for being mixed and amplifieroperation to frequency shift signal and delay time signal, with generate time delay to be launched and Frequency shift signal;
Transmitting antenna, suitable for launching the time delay and frequency shift signal, so that radar receives the signal.
2. calibrator as claimed in claim 1, wherein, the demodulating unit is suitable to filter load from Srx according to following manner Ripple signal, radar and the intermediate-freuqncy signal Sdc of calibrator distance are included to extract:
<mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>pf</mi> <mi>c</mi> </msub> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>pa</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>+</mo> <msub> <mi>j&amp;pi;a</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3. calibrator as claimed in claim 1, wherein, the signal processing unit includes:
Frequency shift unit, the frequency shift signal suitable for the bright spot of orientation in the radar image of generation correspondence radar generation;
Time delay elements, the delay time signal suitable for the bright spot of the radar image middle-range descriscent of generation correspondence radar generation.
4. a kind of radar installations, including:
Reception antenna, suitable for receiving time delay and frequency shift signal from radar calibration device;
Wave filter, suitable for eliminating the mutual of time frequency signal in time delay and frequency shift signal based on Wigner-Ville distribution and Hankel Interference;
Geometry correction unit, suitable for generating radar image according to time delay and frequency shift signal, and according to radar image lieutenant colonel on schedule The coordinate pair radar image carries out geometric correction;And
Radiant correction unit, suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal, and in radar image Pixel carries out radiant correction.
5. radar installations as claimed in claim 4, wherein, the time delay and frequency shift signal received is suitable to pass through following manner table Show:
<mrow> <msub> <mi>s</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;sigma;</mi> <mi>exp</mi> <mo>&amp;lsqb;</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> <mi>&amp;lambda;</mi> </mfrac> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>&amp;times;</mo> <mi>exp</mi> <mo>{</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mi>a</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>}</mo> </mrow>
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one slow when being moved on radar installations under speed u The function of slow time.
6. radar installations as claimed in claim 5, wherein, the geometry correction unit is suitable to according to following manner to radar map As carrying out geometric correction:
Precision reference control is used as in the default spacing of radar image bright spot according to the frequency shift unit and time delay elements generation System point carries out geometric correction, wherein default spacing is determined by setting the frequency shift amount and delay volume of radar calibration device.
7. radar installations as claimed in claim 5, wherein, the radiant correction unit be suitable to according to following manner come according to when Prolong the energy balane radiant correction COEFFICIENT K with frequency shift signal:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
<mrow> <mi>K</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mi>T</mi> </mrow> </msub> <msub> <mi>&amp;rho;</mi> <mrow> <mi>r</mi> <mi>g</mi> </mrow> </msub> <msub> <mi>&amp;rho;</mi> <mrow> <mi>a</mi> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;sigma;</mi> <mi>T</mi> </msub> <mi>sin</mi> <mi>&amp;alpha;</mi> </mrow> </mfrac> </mrow>
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazThe respectively oblique distance and azimuth resolution of image, σTFor the radar cross section of object, α is radar incidence angle.
8. radar installations as claimed in claim 5, wherein, the wave filter is suitable to according to following manner come based on Wigner- Ville is distributed and Hankel eliminates interfering for time frequency signal:
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported, produce While the high-resolution in time delay and frequency displacement;
Multiple signals are disassembled using Hankel conversions, each single signal is other separated, then input into WVD, with Just interfering for time frequency signal, and the linear FM characteristic in stick signal are eliminated.
9. a kind of radar system, including:
At least one radar calibration device as described in claim 1-4;
Such as any one of claim 6-8 radar installations.
10. a kind of radar calibration method, suitable for being performed in radar calibration device, this method includes:
Receive the pulse signal Srx of radar;
Intermediate-freuqncy signal Sdc is extracted from Srx;
Based on intermediate-freuqncy signal Sdc, generation frequency shift signal and delay time signal;
Frequency shift signal and delay time signal are mixed and amplifieroperation, to generate time delay to be launched and frequency shift signal;And
Launch time delay and frequency shift signal, so that radar receives the time delay and frequency shift signal.
11. method as claimed in claim 10, wherein, it is described that carrier signal is filtered from Srx, with extract include radar with The step of intermediate-freuqncy signal Sdc of calibrator distance, includes:
Operation is filtered according to following formula,
<mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>pf</mi> <mi>c</mi> </msub> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>pa</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> 2
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mi>&amp;phi;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>4</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>+</mo> <msub> <mi>j&amp;pi;a</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
12. method as claimed in claim 10, wherein, it is described to be based on intermediate-freuqncy signal Sdc, generation frequency shift signal and delay time signal The step of include:
The frequency shift signal of the bright spot of orientation in the radar image of generation correspondence radar generation;
The delay time signal of the bright spot of the radar image middle-range descriscent of generation correspondence radar generation.
13. a kind of method that radar installations is performed, including:
Receive time delay and frequency shift signal from radar calibration device;
Interfering with each other for time frequency signal in time delay and frequency shift signal is eliminated based on Wigner-Ville distribution and Hankel;
According to time delay and frequency shift signal generation radar image, and the radar image enters according to the coordinate pair of radar image lieutenant colonel on schedule Row geometric correction;And
Radiant correction unit, suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal, and in radar image Pixel carries out radiant correction.
14. method as claimed in claim 13, wherein, the time delay and frequency shift signal received is suitable to pass through following manner table Show:
<mrow> <msub> <mi>s</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;sigma;</mi> <mi>exp</mi> <mo>&amp;lsqb;</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> <mi>&amp;lambda;</mi> </mfrac> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>&amp;times;</mo> <mi>exp</mi> <mo>{</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mi>a</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>t</mi> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>}</mo> </mrow>
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one slow when being moved on radar installations under speed u The function of slow time.
15. method as claimed in claim 13, wherein, it is described that the operation that radar image carries out geometric correction is included:
Carried out according to the default spacing of bright spot in time delay and the corresponding radar image of frequency shift signal as Precision reference control point several What is corrected, wherein default spacing is determined by setting the frequency shift amount and delay volume of radar calibration device.
16. method as claimed in claim 13, wherein, it is described according to time delay and the energy balane radiant correction system of frequency shift signal Number K:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
<mrow> <mi>K</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mi>T</mi> </mrow> </msub> <msub> <mi>&amp;rho;</mi> <mrow> <mi>r</mi> <mi>g</mi> </mrow> </msub> <msub> <mi>&amp;rho;</mi> <mrow> <mi>a</mi> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;sigma;</mi> <mi>T</mi> </msub> <mi>sin</mi> <mi>&amp;alpha;</mi> </mrow> </mfrac> </mrow>
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazThe respectively oblique distance and azimuth resolution of image, σTFor the radar cross section of object, α is radar incidence angle.
17. method as claimed in claim 13, wherein, it is described that time-frequency letter is eliminated based on Wigner-Ville distribution and Hankel Number include the step of interfering:
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported, produce While the high-resolution in time delay and frequency displacement;
Multiple signals are disassembled using Hankel conversions, each single signal is other separated, then input into WVD, with Just interfering for time frequency signal, and the linear FM characteristic in stick signal are eliminated.
CN201710605439.0A 2016-07-24 2017-07-24 Radar calibration method, calibrator, radar device and system Active CN107192992B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610589957 2016-07-24
CN2016105899573 2016-07-24

Publications (2)

Publication Number Publication Date
CN107192992A true CN107192992A (en) 2017-09-22
CN107192992B CN107192992B (en) 2020-06-16

Family

ID=59883923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710605439.0A Active CN107192992B (en) 2016-07-24 2017-07-24 Radar calibration method, calibrator, radar device and system

Country Status (1)

Country Link
CN (1) CN107192992B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108196253A (en) * 2017-12-27 2018-06-22 北京航空航天大学 A kind of azimuth ambiguity suppression method for synthetic aperture radar
CN109444836A (en) * 2018-12-28 2019-03-08 中国人民解放军63891部队 Radar simulation model bearing calibration based on measured data
CN111459086A (en) * 2020-03-30 2020-07-28 深圳市科楠科技开发有限公司 System and method for realizing scaler control and data processing
CN114114178A (en) * 2021-12-10 2022-03-01 南京邮电大学 Calibration device for radar and video image
CN116295158A (en) * 2023-05-17 2023-06-23 四川华恒升科技发展有限公司 Instrument and method for measuring sludge depth based on airborne double-frequency depth finder

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3935573A (en) * 1974-04-08 1976-01-27 Rca Corporation Doppler radar sensor calibrator
CN101069103A (en) * 2004-02-27 2007-11-07 霍尼韦尔国际公司 RF channel calibration for non-linear FM waveforms
CN102121984A (en) * 2010-12-20 2011-07-13 南京鹏力科技有限公司 Calibration source for ground wave radar signal
CN102608582A (en) * 2012-02-02 2012-07-25 北京航空航天大学 Carrier-borne full-coherent phased-array radar calibrator
CN102798844A (en) * 2011-05-23 2012-11-28 中国科学院电子学研究所 Method for measuring input output delay signal isolation of same port

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3935573A (en) * 1974-04-08 1976-01-27 Rca Corporation Doppler radar sensor calibrator
CN101069103A (en) * 2004-02-27 2007-11-07 霍尼韦尔国际公司 RF channel calibration for non-linear FM waveforms
CN102121984A (en) * 2010-12-20 2011-07-13 南京鹏力科技有限公司 Calibration source for ground wave radar signal
CN102798844A (en) * 2011-05-23 2012-11-28 中国科学院电子学研究所 Method for measuring input output delay signal isolation of same port
CN102608582A (en) * 2012-02-02 2012-07-25 北京航空航天大学 Carrier-borne full-coherent phased-array radar calibrator

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108196253A (en) * 2017-12-27 2018-06-22 北京航空航天大学 A kind of azimuth ambiguity suppression method for synthetic aperture radar
CN108196253B (en) * 2017-12-27 2019-09-20 北京航空航天大学 A kind of azimuth ambiguity suppression method for synthetic aperture radar
CN109444836A (en) * 2018-12-28 2019-03-08 中国人民解放军63891部队 Radar simulation model bearing calibration based on measured data
CN109444836B (en) * 2018-12-28 2023-02-28 中国人民解放军63891部队 Radar simulation model correction method based on measured data
CN111459086A (en) * 2020-03-30 2020-07-28 深圳市科楠科技开发有限公司 System and method for realizing scaler control and data processing
CN111459086B (en) * 2020-03-30 2023-08-29 深圳市科楠科技开发有限公司 System and method for realizing scaler control and data processing
CN114114178A (en) * 2021-12-10 2022-03-01 南京邮电大学 Calibration device for radar and video image
CN116295158A (en) * 2023-05-17 2023-06-23 四川华恒升科技发展有限公司 Instrument and method for measuring sludge depth based on airborne double-frequency depth finder
CN116295158B (en) * 2023-05-17 2023-08-08 四川华恒升科技发展有限公司 Instrument and method for measuring sludge depth based on airborne double-frequency depth finder

Also Published As

Publication number Publication date
CN107192992B (en) 2020-06-16

Similar Documents

Publication Publication Date Title
CN107192992A (en) Method, calibrator, radar installations and the system of radar calibration
KR102196734B1 (en) Apparatus and Method for Jamming in Synthetic Aperture Radar
Cheney et al. Fundamentals of radar imaging
CN109471080B (en) High-speed platform radar echo signal simulation system based on simulink
CN108226927B (en) SAR imaging method based on weighted iteration minimum sparse Bayesian reconstruction algorithm
CN104166126B (en) A kind of simulated radar echo method for continuous wave radar
CN107656253B (en) Electromagnetic vortex synthetic aperture radar echo signal simulation method and device
CN104635230B (en) Method for MIMO (multi-input multi-output)-SAR (synthetic aperture radar) near field measurement imaging azimuth side lobe suppression
CN111025246B (en) Simulation system and method for composite scene imaging of sea surface and ship by using stationary orbit SAR
CN107765226A (en) A kind of SAR satellite radars analogue echoes method, system and medium
CN109507666A (en) The sparse frequency band imaging method of ISAR based on off-network variation bayesian algorithm
CN112558497B (en) Anti-interference digital simulation method and system for radar altimeter
CN102854505A (en) Weighting sparse-driven self-focusing SAR (Synthetic Aperture Radar) imaging method
CN102721965B (en) Omni-directional side lobe suppression method for synthetic aperture radar imaging processing
CN104215937A (en) Airborne radar space time adaptation processing method based on environment dynamic perception
Wang et al. Multichannel wideband synthetic aperture radar for ice sheet remote sensing: Development and the first deployment in Antarctica
Li et al. Azimuth super-resolution of forward-looking imaging based on Bayesian learning in complex scene
CN104020465B (en) External illuminators-based radar angle-measuring method based on eight unit small-bore circle array antennas
RU2562614C1 (en) Method of simulating radar targets
CN117491978A (en) Satellite-borne laser radar echo simulation method and device based on point cloud data
CN106772368B (en) The super-resolution three-D imaging method of more random frequency radar arrays
KR102570575B1 (en) Test system of seeker for seeking multiple target and method of thereof
CN110376586B (en) Distributed MIMO radar moving target detection method based on chromatography principle
RU186130U1 (en) MULTIFUNCTIONAL RADAR TARGET SIMULATOR
CN116609742B (en) SAR echo simulation and interference assessment method and system based on real-time map inversion

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhu Wenju

Inventor before: Zhu Wenju

Inventor before: Liu Ning

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20200519

Address after: 100176 313A, floor 3, building 4, yard 4, Yongchang Middle Road, economic and Technological Development Zone, Daxing District, Beijing

Applicant after: Beijing Ju Net Technology Co., Ltd.

Address before: 100000 Beijing Haidian District Meiyuan malianwa No. 6 Building 1 storey commercial section 113

Applicant before: BEIJING TUOPU SPACE TECHNOLOGY Co.,Ltd.

GR01 Patent grant
GR01 Patent grant