CN109471080A - High speed platform radar simulated radar echo system based on simulink - Google Patents
High speed platform radar simulated radar echo system based on simulink Download PDFInfo
- Publication number
- CN109471080A CN109471080A CN201811328332.7A CN201811328332A CN109471080A CN 109471080 A CN109471080 A CN 109471080A CN 201811328332 A CN201811328332 A CN 201811328332A CN 109471080 A CN109471080 A CN 109471080A
- Authority
- CN
- China
- Prior art keywords
- radar
- module
- scene
- echo
- computing module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
Abstract
The high speed platform radar simulated radar echo system based on simulink that the invention discloses a kind of, mainly solves the problems, such as prior art simulation precision and echo low efficiency.Its scheme is: different function module is designed in simulink, wherein PRT synchronization module generates pulse triggering signal, and control sequential is synchronous；Radar pulse transmitting module analog transmissions radar pulse baseband signal；Radar track import modul exports radar real time position and motion information；Beam center computing module calculates and exports real time beam centre coordinate；Sub-scene interception module intercepts imaging sub-scene according to beam center；Antenna direction module calculates and updates the antenna gain of each point in sub-scene；System function computing module calculates in real time and updates the system function of echo；Echo generation module convolution generates original echoed signals.The present invention improves analogue echoes efficiency and precision, for simulating the radar echo signal under different radar operation modes, environment, target scene.
Description
Technical field
The invention belongs to Radar Technology field, in particular to a kind of radar echo signal simulation system can be used for radar system
System design, the verifying and assessment of SAR system performance.
Background technique
Synthetic aperture radar SAR is a kind of modern times high-resolution microwave imaging radar, has been widely used in every field,
And in the R&D process of SAR system, radar return data are for the research of imaging algorithm, the design of radar system parameters, SAR
The verifying and assessment of system performance have extremely important effect.But SAR system is sufficiently complex, will meet in practical application
Face various known, complicated or even unknown, extreme situations and external condition, thus, research and development and development stage in whole system
It needs that a large amount of echo data is combined to be designed each system parameter, correct, algorithm is optimized, is improved.It is so huge
If big echo data amount, which is only relied on, is mounted to actual flight instruments for SAR, obtained as aircraft, satellite are surveyed,
Cost and safety are all very big problems, and are easy to appear error.And the appearance and development of SAR echo simulation analogue technique,
Very big convenience is brought for the research and development and development of SAR system.
The application of analogue echoes technology greatly reduces the research and development cost of SAR system, in the R&D process of SAR, with
When modules can be tested and be debugged, and do not have to until complete machine be made after carry out carry test, moreover, passing through back
Wave simulation technology can obtain various scene radar returns, can largely improve the reliability of SAR.
Currently, SAR analogue echoes have obtained many achievements both at home and abroad, there is pure theory model, semi-physical object simulating system, calculates
The simulation of machine software, pure hardware platform realize FPGA and DPS.But different models, system all Shortcomings place or calculation amount mistake
It is excessively huge in particular for large scene or hardware system greatly, and the analogue echoes of natural scene are supported poor.
1978, V.H.Kaupp and J.C.Holtaman of University of Kansas et al. developed a kind of Ku thunder for being named as RIS
Up to simulator, which is based on target scattering Model, can simulate a variety of different types of scenes, but because of backscattering coefficient number
According to being limited, it is not widely used.
2004, Mori et al. proposed a kind of multi-operation mode SAR echo simulator based on Time-Domain algorithm, can simulate
Original echoed signals under a variety of non-idealities, but in the case where scene is excessive, calculation amount is larger.
2006, Yu Mingcheng of Tsinghua University et al. proposed a kind of SAR original echoed signals based on inverse Wavenumber Domain Algorithms
Then analogy method, this method pass through Wavenumber Domain Algorithms inverting by obtaining scene backscattering coefficient to optical image security
It obtains original echoed signals, though the simulation efficiency of this method is higher, deviates these unreasonablys in platform shake, motion profile
In the case of thinking, simulation precision is lower.
2010, u s company Mistral developed a new SAR echo simulator RTS-RF, and system function is various,
Human-computer interaction is convenient, but hardware system is excessively huge, and supports the analogue echoes of natural scene poor.
Summary of the invention
It is an object of the invention to propose that one kind is based on for deficiency existing for above-mentioned radar echo simulation technology
The high speed platform radar simulated radar echo system of simulink, with simulate different radar operation modes, different external environment,
Radar echo signal under different target scene, and under the premise of guaranteeing precision, huge hardware system is got rid of, reduces and calculates
Amount, improves the efficiency of analogue echoes.
To achieve the above object, the technical scheme is that generating disparate modules by simulink, which is characterized in that
Module generated includes:
PRT synchronization module generates for completing lock-out pulse, and is output to modules, to the clock of whole system into
Row synchronizes；
Radar pulse signal transmitting module, for receiving the input of PRT lock-out pulse, generation radar pulse emits signal, and
Radar pulse transmitting signal is exported to echo generation module；
Radar track import modul, under the control of PRT lock-out pulse, reading radar trail file information, and will
The radar trace information of reading exports respectively to beam center computing module, sub-scene interception module, antenna radiation pattern and calculates mould
Block and system function computing module；
Beam center computing module, under the control of PRT lock-out pulse, according to the radar of input and present mode mark
Show position, updates pixel coordinate information of the beam center in scene figure, and beam center coordinate information is exported respectively to son
Scene interception module, antenna radiation pattern computing module and system function computing module；
Sub-scene interception module, for being shone according to the different beams of radar different mode under the control of PRT lock-out pulse
Beam center coordinate, the radar trace information for penetrating mode and input intercept the sub-scene of imaging from the large scene of importing, and will
The data of the sub-scene are exported respectively to antenna radiation pattern computing module and system function computing module；
Antenna radiation pattern computing module, under the control of PRT lock-out pulse, according to the radar trace information of input,
Beam center coordinate information, sub-scene size, receive wave beam angle of eccentricity, distance to and azimuth resolution, antenna beam width
Degree, calculates the antenna gain of each point in sub-scene, and result is exported to system function computing module；
System function computing module, under the control of PRT lock-out pulse, according in the radar site of input, wave beam
Heart coordinate, sub-scene data, the starting distance door information of antenna radiation pattern, echo, using concentric algorithm, calculate echo is
System function, and the system function is exported to echo generation module；
Echo generation module, for radar pulse transmitting signal and system function to be carried out convolution, and in PRT lock-out pulse
Control under, generate original echoed signals, and by result export to delay and range-gate selection module；
Delay and range-gate selection module, for making echo data in a manner of " stream " under the control of PRT lock-out pulse
Output, thus preferably guinea pig echo-signal stream.
The present invention has the advantage that
1. the present invention carries out modularization, Hierarchical Design by being then based on simulink, system building is relatively easy, and function
Can be huge, it can be convenient, flexible according to the SAR mode identifier of input, the radar trace information file of importing and target scene
Simulation goes out the radar echo signal under different radar operation modes, different external environments, different target scene；
2. system function of the present invention due to being calculated echo using concentric algorithm can be put under the premise of guaranteeing precision
Huge hardware system is taken off, calculation amount is reduced, there is higher analogue echoes efficiency.
Detailed description of the invention
Fig. 1 is radar echo signal simulation system block diagram；
Fig. 2 is the flow chart that radar echo signal simulation system is generated in simulink；
Fig. 3 is the flow chart using present system guinea pig echo-signal；
Fig. 4 is the schematic diagram for carrying out point target imaging under positive side view using imaging algorithm；
Fig. 5 is that system function computing module utilizes concentric algorithm generation Echo System letter geometry of numbers mould in the present invention
Type.
Specific embodiment
Present invention is further described in detail referring to the drawings:
Currently, SAR analogue echoes have obtained many achievements both at home and abroad, there is pure theory model, semi-physical object simulating system, calculates
The simulation of machine software, pure hardware platform realize FPGA and DPS.The present invention is using simulink to radar echo signal simulation system
Carry out modularized design.
Simulink is as a Visual Simulation Tools in MATLAB, Modelling of Dynamic System that it is provided, emulation and comprehensive
The integration environment analyzed is closed, without writing interminable program, is by what the interface operation of simple, intuitive can construct complexity
System, is widely used in the complex simulation and design of control theory and Digital Signal Processing.
The present invention is as shown in Figure 2 using the process that simulink carries out system building.Firstly, passing through the definition guiding of M script
Then basic parameter required for head system constructs the model of modules in simulink, input, output port is arranged,
The function of realizing modules, recalls S function and is packaged to modules, generates parameter setting interface, finally will be each
Corresponding input, output port are connected in module, generate radar echo signal simulation system, as shown in Figure 1.
Referring to Fig.1, the radar echo signal simulation system that the present invention utilizes simulink to generate, including PRT synchronization module
1, radar pulse signal transmitting module 2, radar track import modul 3, beam center computing module 4, sub-scene interception module 5,
Antenna radiation pattern computing module 6, system function computing module 7, echo generation module 8 and delay and range-gate selection module 9,
Wherein:
PRT synchronization module 1 generates for completing lock-out pulse, and is output to modules, to the clock of whole system into
Row synchronizes；
Radar pulse signal transmitting module 2, for receiving the input of PRT lock-out pulse, generation radar pulse emits signal, and
Radar pulse transmitting signal is exported to echo generation module 8；
Radar track import modul 3, under the control of PRT lock-out pulse, reading radar trail file information, and will
The radar trace information of reading is exported respectively to beam center computing module 4, sub-scene interception module 5, antenna radiation pattern and is calculated
Module 6 and system function computing module 7；
Beam center computing module 4, under the control of PRT lock-out pulse, according to the radar and present mode of input
Position is indicated, updates pixel coordinate information of the beam center in scene figure, and beam center coordinate information is exported to subfield
Scape interception module 5, antenna radiation pattern computing module 6 and system function computing module 7；
Sub-scene interception module 5, under the control of PRT lock-out pulse, according to the different beams of radar different mode
Radiation modality and the beam center coordinate of input, radar trace information, intercept the sub-scene of imaging from the large scene of importing, and
The data of the sub-scene are exported to antenna radiation pattern computing module 6 and system function computing module 7；
Antenna radiation pattern computing module 6, under the control of PRT lock-out pulse, according to the radar trace information of input,
Beam center coordinate information, sub-scene size, receive wave beam angle of eccentricity, distance to and azimuth resolution, antenna beam width
Degree, calculates the antenna gain of each point in sub-scene, and result is exported to system function computing module 7；
System function computing module 7, under the control of PRT lock-out pulse, according in the radar site of input, wave beam
Heart coordinate, sub-scene data, the starting distance door information of antenna radiation pattern, echo, using concentric algorithm, calculate echo is
System function, and the system function is exported to echo generation module 8；
Echo generation module 8, for radar pulse transmitting signal and system function to be carried out convolution, and in the synchronous arteries and veins of PRT
Under the control of punching, original echoed signals are generated, and result is exported to delay and range-gate selection module 9；
Delay and range-gate selection module 9, for making echo data with the side of " stream " under the control of PRT lock-out pulse
Formula output, thus preferably guinea pig echo-signal stream.
Referring to Fig. 3, the process using present system guinea pig echo-signal is as follows:
Process 1, PRT synchronization module 1 read radar operating frequency, sampling frequency under various modes from simulation interactive interface
Rate, transmitting signal pulsewidth, transmitted signal bandwidth, pulse-recurrence time and transmitting pulse number, export PRT synchronization signal
The identifier FrameTrigger of radar operation mode is exported to beam center and is calculated mould to modules by PRTTrigger
Block.
Process 2, radar pulse signal transmitting module 2 receive the input of PRT lock-out pulse, detect rising edge, when rising edge comes
Temporarily, base band linear FM signal LFM is generated, and exports N according to pulse repetition period PRTaA LFM pulse signal is raw to echo
At module 8；
The generation base band linear FM signal LFM expression formula is as follows:
Wherein,T is fast time, tmFor slow time, TpFor LFM signal pulsewidth, fcFor carrier wave
Frequency, krFor linear frequency modulation rate；
Process 3, radar track import modul 3 receive the input triggering of PRT lock-out pulse, detect rising edge, when rising edge comes
Temporarily, trail file information is read, radar site, speed, acceleration, incidence angle, velocity vector and wave beam is updated and is thrown on ground
Shadow angle and radar site information with error, and result is exported to beam center computing module 4, sub-scene interception module
5, antenna radiation pattern computing module 6 and system function computing module 7.
Process 4, beam center computing module 4 receive the input triggering of PRT lock-out pulse, detect rising edge, when rising edge comes
Temporarily, position, speed, target location coordinate, the present mode mark position for inputting radar, calculate beam center in scene figure
Pixel coordinate information, and result is exported to sub-scene interception module 5, antenna radiation pattern computing module 6 and system function meter
Calculate module 7.
Process 5, sub-scene interception module 5 input beam center according to the different beams radiation modality of radar different mode
Coordinate, radar site coordinate, and triggered by PRT is synchronous, intercept out the sub-scene being imaged from the large scene of importing, then by data
It exports to antenna radiation pattern computing module 6 and system function computing module 7.
Process 6, antenna radiation pattern computing module 6 input the pixel of radar site information, beam center in large scene
Coordinate information, sub-scene size, receive wave beam angle of eccentricity, distance to and azimuth resolution, antenna beamwidth, and receive
The input triggering of PRT lock-out pulse, detects rising edge, when rising edge temporarily, to calculate the antenna gain of each point in sub-scene, and will
As a result it exports to system function computing module 7.
The antenna gain of each point in the calculating sub-scene, there are two types of modes:
The first is to be counted under SAR imaging pattern using the directional diagram and receiving antenna offset angle of single antenna for 0
It calculates, i.e., the side of beam center is first obtained according to the coordinate of the radar site information of input, beam center imago vegetarian refreshments in the scene
Parallactic angle αcWith pitch angle βc, azimuth angle alpha of the target with respect to radarRTWith pitch angle βRT；It is calculated again by transmitting antenna directional diagram tune
Antenna gain rcs after system1:
rcs1=abs ((sinc (αRT-αc))*(sinc(βRT-βc))) <2>
Wherein, abs is ABS function, and * indicates to be multiplied,
Second is under single pulse mode using four single antennas, and the corresponding offset angle that receives to be arranged and forms four
Antenna radiation pattern is calculated, i.e., first according to the coordinate of the radar site information of input, beam center imago vegetarian refreshments in the scene
And beamlet is opposite and azimuth deviation angle ± Δ α, the pitch deviation angle ± Δ β of wave beam, obtains the orientation of beam center
Angle αcWith pitch angle βc, azimuth angle alpha of the target with respect to radarRTWith pitch angle βRT；The transmitting by different biasings is calculated separately again
The modulated antenna gain rcs of antenna radiation pattern21、rcs22、rcs23、rcs24:
rcs21=abs ((sinc (αRT+Δα-αc))*(sinc(βRT+Δβ-βc))) <3>
rcs22=abs ((sinc (αRT+Δα-αc))*(sinc(βRT-Δβ-βc))) <4>
rcs23=abs ((sinc (αRT-Δα-αc))*(sinc(βRT+Δβ-βc))) <5>
rcs24=abs ((sinc (αRT-Δα-αc))*(sinc(βRT-Δβ-βc))) <6>
Total antenna gain rcs2Are as follows:
rcs2=rcs21*rcs22*rcs23*rcs24 <7>
Wherein, rcs21It is that beamlet is opposite and the azimuth deviation angle of wave beam is Δ α, the day that pitch deviation angle is Δ β
Line gain, rcs22It is that beamlet is opposite and the azimuth deviation angle of wave beam is Δ α, pitch deviation angle is the-antenna increasing of Δ β
Benefit, rcs23It is that beamlet is opposite and the azimuth deviation angle of wave beam is-Δ α, the antenna gain that pitch deviation angle is Δ β,
rcs24Be that beamlet is opposite and the azimuth deviation angle of wave beam be-Δ α, pitch deviation angle for-Δ β antenna gain.
Process 7, system function computing module 7 input the position of radar, beam center coordinate, sub-scene data, antenna side
The starting distance door information of Xiang Tu, echo, and the input triggering of PRT lock-out pulse is received, rising edge is detected, when rising edge arrives
When, using concentric algorithm, the system function of echo is calculated, and result is exported to echo generation module 8；
The concentric algorithm principle is as follows:
In the case where not considering inclination of wave front, example, distribution of the echo of point target on two-dimensional surface are considered as with positive side
The rectangle battle array as shown in the first width of Fig. 4, by distance to pulse compression after will become the second width figure as a result, this is because radar
Movement produces bending.After correcting bending, the Energy distribution of the point target can in identical distance unit, such as in 4 the
Shown in three width figures, at this moment orientation imaging can be carried out to target along orientation, obtain the imaging results of the point target, such as 4 most
Shown in the latter figure.
As it can be seen that different point targets, since the distance to radar is different, they be will be distributed in different distance unit,
This is because the distance of different target points to radar platform is different, the difference of delay time is caused, with distance samples frequency
fsThe range delay of target point is sampled, c/2f is divided between sampling units, c is the light velocity, the echo point of each target point
Cloth is by being divided into what integral multiple relation was distributed between sampling unit；Make at the different orientation moment for reaching radar in scene
With for identical target point, their sampling unit integral multiple relation is identical, and therefore, they will be accumulated in plural echo
In identical distance unit, it is conceivable that, at some orientation moment, it is distributed in using radar platform as the identical concentric of origin
Since the operating distance to radar is identical, what their energy should add up is distributed in identical distance unit for point target on circle
It is interior.
As a result, in order to quickly obtain Echo System function, to meet the requirement of real-time echo signal generation, while to keep
The computational accuracy of echo-signal, it is contemplated that the point apart from radar equal length is located in the same distance unit, first by scene
In point add up along by the concentric circles in the center of circle of radar, obtain the one-dimensional range profile of radar, then using FFT in frequency domain
The generation of Echo System function is fast implemented, multiple spot can be handled simultaneously in this way, reduce operand.
According to the above thought, calculating Echo System function, the specific method is as follows:
At each orientation moment, first have to calculate all point target in scene to radar distance R (k), and this
It is compared apart from same distance sampling unit, obtains the distribution situation put on all concentric circles, i.e.,
In formula<8>, nkIndicate the position of distance unit, i.e. which concentric circles the point is distributed on, δrIt is big for distance unit
It is small, and
As shown in figure 5, within the scope of beam, P is shared on some concentric circles after obtaining the distribution situation of concentric circles
A point target, according to formula<9>it is recognised that this P scattering point should be distributed in identical distance unit, they can unify
Echo-signal is generated, and since the orientation phase information ratio of echo-signal is more sensitive apart from envelope information, so, it is ensured that side
The integrality of position phase information, that is, it is the same from envelope to be unable to image distance, is carried out with formula<9>apart from approximate calculation, so, the side of each point
Position phase signal s (mT；RB) need to independently calculate:
Wherein, σ is the gray value of point target, and mT is tmDiscrete form, RBIt is minimum distance of the radar to target, λ is thunder
Up to operation wavelength, R (mT；It R) is distance of the mT moment radar to scattering point, exp indicates exponential function.
Later, summed to obtain the point target orientation phase signal s on identical concentric circles to formula<10>2:
Wherein, σiFor the gray value of i-th of point target on identical concentric circles.
It sums to formula<11>, obtains all distance unit data s of entire moment3:
In formula<12>, δ is impulse Response Function, which distance unit k expression falls in, PnIndicate that n-th of distance is single
Point target number in member.The number of its ground scatter point is different in different distance unit, the difficulty of cumulative process also not phase
Together.To formula<12>carry out Fourier transformation FFT it change to frequency domain multiplied by distance to frequency modulation item, recycle inverse Fourier transform
IFFT, which changes back to time domain, can be obtained by the system function s of echo4(k,mT；RB):
Wherein, frIndicate distance to frequency.
Process 8,8 input system function of echo generation module and base band linear FM signal LFM carry out convolution algorithm,
Generate original echoed signals, and receive the input triggering of PRT lock-out pulse, detect rising edge, when rising edge come it is interim, result is defeated
Out to delay and range-gate selection module 9.
Process 9, delay and range-gate selection module 9 are under the control of PRT lock-out pulse, by the original echoed signals of input
It is exported in a manner of " stream ".
To sum up, the high speed platform radar simulated radar echo system proposed in this paper based on simulink, have modularization,
The design of stratification can be simulated easily and flexibly under different radar operation modes, different external environments, different target scene
Radar echo signal helps scientific research personnel to get rid of the limitation of radar equipment condition, needs not rely on expensive radar equipment and obtains
Correlation radar echo data, it is more efficient, more convenient than traditional measured data mode.
Above disclosed is only a preferred embodiment of the present invention, it is clear that the power of the present invention cannot be limited with this
Sharp range, therefore equivalent changes made in accordance with the claims of the present invention still fall within the range that the present invention is covered.
Claims (4)
1. the high speed platform radar simulated radar echo system based on simulink is to generate disparate modules by simulink,
It is characterized in that, module generated includes:
PRT synchronization module (1) generates for completing lock-out pulse, and is output to modules, carries out to the clock of whole system
It is synchronous；
Radar pulse signal transmitting module (2) generates radar pulse and emits signal, and will for receiving the input of PRT lock-out pulse
Radar pulse transmitting signal is exported to echo generation module (8)；
Radar track import modul (3) under the control of PRT lock-out pulse, reading radar trail file information, and will be read
The radar trace information taken is exported respectively to beam center computing module (4), sub-scene interception module (5), antenna radiation pattern meter
Calculate module (6) and system function computing module (7)；
Beam center computing module (4), under the control of PRT lock-out pulse, according to the radar of input and present mode mark
Show position, updates pixel coordinate information of the beam center in scene figure, and beam center coordinate information is exported to sub-scene
Interception module (5), antenna radiation pattern computing module (6) and system function computing module (7)；
Sub-scene interception module (5), for being shone according to the different beams of radar different mode under the control of PRT lock-out pulse
Beam center coordinate, the radar trace information for penetrating mode and input intercept the sub-scene of imaging from the large scene of importing, and will
The data of the sub-scene are exported to antenna radiation pattern computing module (6) and system function computing module (7)；
Antenna radiation pattern computing module (6), under the control of PRT lock-out pulse, according to radar trace information, the wave of input
Beam center coordinate information, sub-scene size, receive wave beam angle of eccentricity, distance to and azimuth resolution, antenna beamwidth,
The antenna gain of each point in sub-scene is calculated, and result is exported to system function computing module (7)；
System function computing module (7), under the control of PRT lock-out pulse, according to the radar site of input, beam center
Coordinate, sub-scene data, the starting distance door information of antenna radiation pattern, echo, using concentric algorithm, the system for calculating echo
Function, and the system function is exported to echo generation module (8)；
Echo generation module (8), for radar pulse transmitting signal and system function to be carried out convolution, and in PRT lock-out pulse
Control under, generate original echoed signals, and by result export to delay and range-gate selection module (9)；
Delay and range-gate selection module (9), for making echo data in a manner of " stream " under the control of PRT lock-out pulse
Output, thus preferably guinea pig echo-signal stream.
2. system according to claim 1, which is characterized in that generate disparate modules by simulink, be accomplished by
Basic parameter required for guiding head system is defined by M script, then constructs the model of modules in simulink,
Setting input, output port, realize the function of modules；
It calls S function to be packaged modules, generates parameter setting interface；
Input corresponding in modules, output port are connected.
3. system according to claim 1, which is characterized in that antenna radiation pattern computing module (6) calculates each in sub-scene
The antenna gain of point, is accomplished by
First according to the radar site information of input, beam center coordinate information, reception wave beam angle of eccentricity and antenna beam width
Degree, obtains the azimuth angle alpha of beam centercWith pitch angle βc, azimuth angle alpha of the target with respect to radarRTWith pitch angle βRT；
It is calculated again by the modulated antenna gain rcs of transmitting antenna directional diagram:
Rcs=abs ((sinc (αRT-αc))*(sinc(βRT-βc))),<1>
Wherein, abs is ABS function, and * indicates to be multiplied,
4. system according to claim 1, which is characterized in that system function computing module (7) utilizes concentric algorithm, meter
The system function for calculating echo, is accomplished by
All point target in scene is calculated to the distance R (k) of radar, and this is compared apart from same distance sampling unit,
The distribution situation of the point of all concentric circles is obtained, i.e., within the scope of beam, obtains sharing P mesh on some concentric circles
Mark；
This P scattering point is distributed in identical distance unit, the orientation phase signal s of each point is calculated1(mT；RB):
Wherein, σ is the gray value of point target；MT is slow time tmDiscrete form, RBIt is minimum distance of the radar to scattering point, λ
For radar operation wavelength, R (mT；RB) it is distance of the mT moment radar to scattering point；
Point target on identical concentric circles is added up, i.e., is summed to formula<2>, obtain some distance unit includes side
The data s of position phase2:
It sums to formula<3>, obtains all distance unit data s of entire moment3:
Wherein, δ is impulse Response Function, which distance unit k expression falls in；
Fourier transformation FFT is carried out to formula<4>, it change to frequency domain multiplied by distance to frequency modulation item, recycle inverse Fourier to become
It changes IFFT and changes back to time domain, the system function s of echo can be obtained4(k,mT；RB):
Wherein, krIndicate distance to frequency modulation rate, frIndicate distance to frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811328332.7A CN109471080B (en) | 2018-11-09 | 2018-11-09 | High-speed platform radar echo signal simulation system based on simulink |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811328332.7A CN109471080B (en) | 2018-11-09 | 2018-11-09 | High-speed platform radar echo signal simulation system based on simulink |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109471080A true CN109471080A (en) | 2019-03-15 |
| CN109471080B CN109471080B (en) | 2022-12-02 |
Family
ID=65672311
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811328332.7A Active CN109471080B (en) | 2018-11-09 | 2018-11-09 | High-speed platform radar echo signal simulation system based on simulink |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109471080B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782241A (en) * | 2019-03-29 | 2019-05-21 | 北京润科通用技术有限公司 | A kind of meteorological radar echo analogy method and system |
| CN109901165A (en) * | 2019-03-28 | 2019-06-18 | 河南九乾电子科技有限公司 | The simulator and analogy method of satellite-borne SAR echo |
| CN110726977A (en) * | 2019-11-29 | 2020-01-24 | 中国舰船研究设计中心 | Ship radar performance evaluation method in interference environment |
| CN110988858A (en) * | 2019-11-11 | 2020-04-10 | 西安空间无线电技术研究所 | High-precision distance measurement method and system for dual-beam microwave landing radar |
| CN111323761A (en) * | 2020-03-20 | 2020-06-23 | 北京华力创通科技股份有限公司 | Echo system function construction method and device and echo simulator |
| CN111624564A (en) * | 2020-05-27 | 2020-09-04 | 北京润科通用技术有限公司 | Radar pitch angle target simulation system and method |
| CN112213721A (en) * | 2020-09-16 | 2021-01-12 | 西安科技大学 | Millimeter wave three-dimensional imaging method for scanning outer or inner scenes of cylinder for security inspection |
| CN112578350A (en) * | 2020-12-02 | 2021-03-30 | 西安电子科技大学 | Airborne SAR interference effect simulation method under high-energy microwave interference |
| CN112698280A (en) * | 2020-12-09 | 2021-04-23 | 南京长峰航天电子科技有限公司 | Bistatic SAR real-time echo simulation method based on DSP and FPGA architecture |
| CN114442051A (en) * | 2020-11-05 | 2022-05-06 | 北京华航无线电测量研究所 | High-fidelity missile-borne radar echo simulation method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090109086A1 (en) * | 2006-05-13 | 2009-04-30 | Gerhard Krieger | High-Resolution Synthetic Aperture Side View Radar System Used By Means of Digital Beamforming |
| CN104090277A (en) * | 2014-07-21 | 2014-10-08 | 中国科学院电子学研究所 | Method for imaging sliding circumferential synthetic aperture radar |
| CN107765226A (en) * | 2017-09-18 | 2018-03-06 | 北京空间飞行器总体设计部 | A kind of SAR satellite radars analogue echoes method, system and medium |
-
2018
- 2018-11-09 CN CN201811328332.7A patent/CN109471080B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090109086A1 (en) * | 2006-05-13 | 2009-04-30 | Gerhard Krieger | High-Resolution Synthetic Aperture Side View Radar System Used By Means of Digital Beamforming |
| CN104090277A (en) * | 2014-07-21 | 2014-10-08 | 中国科学院电子学研究所 | Method for imaging sliding circumferential synthetic aperture radar |
| CN107765226A (en) * | 2017-09-18 | 2018-03-06 | 北京空间飞行器总体设计部 | A kind of SAR satellite radars analogue echoes method, system and medium |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901165A (en) * | 2019-03-28 | 2019-06-18 | 河南九乾电子科技有限公司 | The simulator and analogy method of satellite-borne SAR echo |
| CN109901165B (en) * | 2019-03-28 | 2023-03-14 | 河南九乾电子科技有限公司 | Satellite-borne SAR echo simulation device and method |
| CN109782241A (en) * | 2019-03-29 | 2019-05-21 | 北京润科通用技术有限公司 | A kind of meteorological radar echo analogy method and system |
| CN109782241B (en) * | 2019-03-29 | 2021-01-29 | 北京润科通用技术有限公司 | Meteorological radar echo simulation method and system |
| CN110988858A (en) * | 2019-11-11 | 2020-04-10 | 西安空间无线电技术研究所 | High-precision distance measurement method and system for dual-beam microwave landing radar |
| CN110726977A (en) * | 2019-11-29 | 2020-01-24 | 中国舰船研究设计中心 | Ship radar performance evaluation method in interference environment |
| CN111323761B (en) * | 2020-03-20 | 2022-04-08 | 北京华力创通科技股份有限公司 | Echo system function construction method and device and echo simulator |
| CN111323761A (en) * | 2020-03-20 | 2020-06-23 | 北京华力创通科技股份有限公司 | Echo system function construction method and device and echo simulator |
| CN111624564A (en) * | 2020-05-27 | 2020-09-04 | 北京润科通用技术有限公司 | Radar pitch angle target simulation system and method |
| CN112213721A (en) * | 2020-09-16 | 2021-01-12 | 西安科技大学 | Millimeter wave three-dimensional imaging method for scanning outer or inner scenes of cylinder for security inspection |
| CN112213721B (en) * | 2020-09-16 | 2023-10-03 | 西安科技大学 | Millimeter wave three-dimensional imaging method for scanning external or internal cylindrical scene facing security inspection |
| CN114442051A (en) * | 2020-11-05 | 2022-05-06 | 北京华航无线电测量研究所 | High-fidelity missile-borne radar echo simulation method |
| CN112578350B (en) * | 2020-12-02 | 2022-04-19 | 西安电子科技大学 | Airborne SAR interference effect simulation method under high-energy microwave interference |
| CN112578350A (en) * | 2020-12-02 | 2021-03-30 | 西安电子科技大学 | Airborne SAR interference effect simulation method under high-energy microwave interference |
| CN112698280A (en) * | 2020-12-09 | 2021-04-23 | 南京长峰航天电子科技有限公司 | Bistatic SAR real-time echo simulation method based on DSP and FPGA architecture |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109471080B (en) | 2022-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109471080A (en) | High speed platform radar simulated radar echo system based on simulink | |
| CN104614713B (en) | A kind of ship that is suitable for carries the radar echo signal simulator of radar system | |
| CN109164428A (en) | Digital radar analogue system and method | |
| CN103207387B (en) | Method for quickly simulating airborne phased array pulse Doppler (PD) radar clutter | |
| CN101581779A (en) | Method for generating three-dimensional imaging original echoed signals of chromatography synthetic aperture radars | |
| CN107765226A (en) | A kind of SAR satellite radars analogue echoes method, system and medium | |
| CN105467369B (en) | A kind of target echo simulation method and apparatus | |
| CN107783092A (en) | To rcs measurement system and method behind near field based on Chain relation | |
| CN104597444A (en) | Microwave gaze high resolution imaging method based on intensity association | |
| Singh et al. | Simulation of the radar cross-section of dynamic human motions using virtual reality data and ray tracing | |
| US6069582A (en) | Method and apparatus for synthesizing multi-channel radar or sonar data | |
| CN102798861B (en) | Bistatic synthetic aperture radar (SAR) imaging method based on optimal image space | |
| CN106707253A (en) | Networking radar and networking jammer countermeasure test device and method in test room | |
| Ehrman et al. | Automated target recognition using passive radar and coordinated flight models | |
| CN116381629A (en) | Radar large-scale target simulation system and method based on real-time dynamic convolution | |
| CN109738890A (en) | A method of distance figure is generated based on missile-borne Bistatic SAR range Doppler image | |
| Bączyk et al. | Moving target imaging in multistatic passive radar | |
| Rouffet et al. | Digital twin: A full virtual radar system with the operational processing | |
| Abadpour | Modeling Backscattering Behavior of Vulnerable Road Users Based on High-Resolution Radar Measurements | |
| CN106646409A (en) | SAR echo signal simulation method based on quasi-double-station model | |
| Jun et al. | Target Echo Simulation and Implementation Technology Based on Missile-Borne SAR Scene Matching Guidance | |
| Pan et al. | Impacts of space-time-frequency synchronization errors onwideband target echo characteristics of bistatic/multistatic radar | |
| Lu et al. | Model-based radar system simulation and verification | |
| Liu et al. | Bistatic FMCW SAR raw signal simulator for extended scenes | |
| Bokov et al. | Reduction of the discrete error of a digital delay simulation method for reflected radar signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |