CN109471080B - High-speed platform radar echo signal simulation system based on simulink - Google Patents
High-speed platform radar echo signal simulation system based on simulink Download PDFInfo
- Publication number
- CN109471080B CN109471080B CN201811328332.7A CN201811328332A CN109471080B CN 109471080 B CN109471080 B CN 109471080B CN 201811328332 A CN201811328332 A CN 201811328332A CN 109471080 B CN109471080 B CN 109471080B
- Authority
- CN
- China
- Prior art keywords
- module
- radar
- echo
- scene
- sub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
-
- 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
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 a simulink-based high-speed platform radar echo signal simulation system, which mainly solves the problems of low simulation precision and echo efficiency in the prior art. The scheme is as follows: different functional modules are designed in simulink, wherein a PRT synchronization module generates a pulse trigger signal to control time sequence synchronization; the radar pulse transmitting module simulates and transmits radar pulse baseband signals; the radar track import module outputs radar real-time position and motion information; the beam center calculation module calculates and outputs real-time beam center coordinates; the sub-scene intercepting module intercepts the imaging sub-scene according to the beam center; the antenna directional diagram module calculates and updates the antenna gain of each point in the sub-scene; the system function calculation module calculates and updates the system function of the echo in real time; the echo generating module convolutes to generate an original echo signal. The invention improves the efficiency and the precision of echo simulation and is used for simulating radar echo signals under different radar working modes, environments and target scenes.
Description
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a radar echo signal simulation system which can be used for designing a radar system and verifying and evaluating the performance of an SAR system.
Background
The synthetic aperture radar SAR is a modern high-resolution microwave imaging radar and is widely applied to various fields, and in the research and development process of the SAR system, radar echo data plays an extremely important role in the research and development of an imaging algorithm, the design of radar system parameters and the verification and evaluation of the performance of the SAR system. However, the SAR system is very complex, and in practical application, various known, complex, even unknown, extreme situations and external conditions will be encountered, so that in the development and development stages of the whole system, a large amount of echo data needs to be combined to design and correct parameters of each system, and an algorithm is optimized and improved. If the huge echo data volume is obtained only by mounting the SAR on an actual flight device, such as an airplane and a satellite, actual measurement is carried out, the cost and the safety are both very big problems, and errors are easy to occur. The appearance and the development of the SAR echo simulation technology bring great convenience to the research and development of the SAR system.
The application of the echo simulation technology greatly reduces the research and development cost of the SAR system, and in the research and development process of the SAR, each module can be tested and debugged at any time without waiting for the mounting test after the complete machine is manufactured.
At present, SAR echo simulation at home and abroad has achieved a lot of achievements, and the SAR echo simulation system comprises a pure theoretical model, a semi-physical simulation system, computer software simulation and a pure hardware platform for realizing FPGA and DPS. However, different models and systems have disadvantages, or the calculated amount is too large, especially for a large scene, or the hardware system is too large, and the echo simulation support for a natural scene is poor.
In 1978, v.h.kaupp and j.c.holtaman et al, the university of kansas, developed a Ku radar simulator named RIS, which is based on a point scattering model and can simulate a variety of different types of scenes, but is not widely used due to limited backscattering coefficient data.
In 2004, mori et al proposed a multi-working-mode SAR echo simulator based on a time domain algorithm, which can simulate original echo signals under various non-ideal conditions, but has a large calculation amount under the condition of an excessively large scene.
In 2006, the university of qinghua provided an SAR original echo signal simulation method based on an inverse wave number domain algorithm, the method obtains a scene backscattering coefficient by processing an optical image, and then obtains an original echo signal by inversion of the wave number domain algorithm.
In 2010, a new SAR echo simulator RTS-RF is developed by Mistral of American company, the system has various functions, the man-machine interaction is convenient, but the hardware system is too large, and the echo simulation support for a natural scene is poor.
Disclosure of Invention
The invention aims to provide a high-speed platform radar echo signal simulation system based on simulink aiming at the defects of the radar echo simulation technology, so as to simulate radar echo signals under different radar working modes, different external environments and different target scenes, get rid of a huge hardware system on the premise of ensuring the precision, reduce the calculated amount and improve the echo simulation efficiency.
In order to achieve the above object, the technical solution of the present invention is to generate different modules by simulink, wherein the generated modules include:
the PRT synchronization module is used for finishing the generation of synchronization pulses, outputting the synchronization pulses to each module and synchronizing the clock of the whole system;
the radar pulse signal transmitting module is used for receiving the PRT synchronous pulse input, generating a radar pulse transmitting signal and outputting the radar pulse transmitting signal to the echo generating module;
the radar track importing module is used for reading radar track file information under the control of PRT synchronous pulse and respectively outputting the read radar track information to the beam center computing module, the sub-scene intercepting module, the antenna directional diagram computing module and the system function computing module;
the beam center calculation module is used for updating pixel point coordinate information of a beam center in a scene graph according to the input radar and the current mode marking bit under the control of the PRT synchronous pulse, and respectively outputting the beam center coordinate information to the sub-scene interception module, the antenna directional diagram calculation module and the system function calculation module;
the sub-scene intercepting module is used for intercepting an imaged sub-scene from a guided large scene according to different beam irradiation modes of different modes of the radar, input beam center coordinates and radar track information under the control of the PRT synchronous pulse, and respectively outputting data of the sub-scene to the antenna directional diagram calculating module and the system function calculating module;
the antenna directional pattern calculation module is used for calculating the antenna gain of each point in the sub-scene according to the input radar track information, the beam center coordinate information, the sub-scene size, the receiving beam offset angle, the distance direction and azimuth direction resolution ratio and the antenna beam width under the control of the PRT synchronous pulse, and outputting the result to the system function calculation module;
the system function calculation module is used for calculating a system function of an echo by utilizing a concentric circle algorithm according to an input radar position, a beam center coordinate, sub-scene data, an antenna directional diagram and initial range gate information of the echo under the control of a PRT synchronous pulse, and outputting the system function to the echo generation module;
the echo generating module is used for convolving the radar pulse transmitting signal with a system function, generating an original echo signal under the control of a PRT synchronous pulse, and outputting a result to the delay and range gate gating module;
and the time delay and range gate gating module is used for outputting echo data in a 'streaming' mode under the control of the PRT synchronous pulse, so that the radar echo signal stream is better simulated.
The invention has the following advantages:
1. because the invention is based on simulink to carry out modularization and hierarchical design, the system is relatively simple in construction and huge in function, namely radar echo signals under different radar working modes, different external environments and different target scenes can be conveniently and flexibly simulated according to the input SAR mode identifier, the imported radar track information file and the target scene;
2. the invention adopts the concentric circle algorithm to calculate the system function of the echo, so that a huge hardware system can be got rid of on the premise of ensuring the precision, the calculation amount is reduced, and the echo simulation efficiency is higher.
Drawings
FIG. 1 is a block diagram of a radar echo signal simulation system;
FIG. 2 is a flow diagram of a simulation system for generating radar return signals in simulink;
FIG. 3 is a flow chart for simulating a radar echo signal using the system of the present invention;
FIG. 4 is a schematic illustration of point target imaging using an imaging algorithm from a front side view;
FIG. 5 is a geometric model of the echo system function generated by the system function calculating module according to the present invention by using a concentric circle algorithm.
Detailed Description
The invention is described in further detail below with reference to the attached drawing figures:
at present, SAR echo simulation at home and abroad has achieved a lot of achievements, and the SAR echo simulation system comprises a pure theoretical model, a semi-physical simulation system, computer software simulation and a pure hardware platform for realizing FPGA and DPS. The invention utilizes simulink to carry out modular design on a radar echo signal simulation system.
The simulink is used as a visual simulation tool in MATLAB, provides an integrated environment of dynamic system modeling, simulation and comprehensive analysis, does not need to write a lengthy program, can construct a complex system through simple and intuitive interface operation, and is widely used for complex simulation and design of control theory and digital signal processing.
The process of the invention for system construction using simulink is shown in fig. 2. Firstly, defining basic parameters required by a seeker system through an M script, then constructing a model of each module in a simulink, setting input and output ports to realize functions of each module, calling an S function to package each module, generating a parameter setting interface, and finally connecting corresponding input and output ports in each module to generate a radar echo signal simulation system, wherein the system is shown in figure 1.
Referring to fig. 1, the system for simulating a radar echo signal generated by simulink of the present invention includes a PRT synchronization module 1, a radar pulse signal emission module 2, a radar track introduction module 3, a beam center calculation module 4, a sub-scene interception module 5, an antenna directional pattern calculation module 6, a system function calculation module 7, an echo generation module 8, and a delay and range gate gating module 9, wherein:
the PRT synchronization module 1 is used for finishing the generation of synchronization pulse, outputting the synchronization pulse to each module and synchronizing the clock of the whole system;
the radar pulse signal transmitting module 2 is used for receiving PRT synchronous pulse input, generating a radar pulse transmitting signal and outputting the radar pulse transmitting signal to the echo generating module 8;
the radar track importing module 3 is used for reading radar track file information under the control of PRT synchronous pulses and respectively outputting the read radar track information to the beam center calculating module 4, the sub-scene intercepting module 5, the antenna directional diagram calculating module 6 and the system function calculating module 7;
the beam center calculation module 4 is used for updating the pixel point coordinate information of the beam center in the scene graph according to the input radar and the current mode marking bit under the control of the PRT synchronous pulse, and outputting the beam center coordinate information to the sub-scene interception module 5, the antenna directional pattern calculation module 6 and the system function calculation module 7;
the sub-scene intercepting module 5 is used for intercepting an imaged sub-scene from a guided large scene according to different beam irradiation modes of different modes of the radar, input beam center coordinates and radar track information under the control of a PRT synchronous pulse, and outputting data of the sub-scene to the antenna directional diagram calculating module 6 and the system function calculating module 7;
an antenna directional pattern calculation module 6, configured to calculate antenna gains of each point in a sub-scene according to input radar track information, beam center coordinate information, sub-scene size, received beam offset angle, distance direction and azimuth direction resolution, and antenna beam width under control of the PRT synchronization pulse, and output a result to a system function calculation module 7;
the system function calculation module 7 is configured to calculate a system function of an echo by using a concentric circle algorithm according to an input radar position, a beam center coordinate, sub-scene data, an antenna directional diagram, and initial range gate information of the echo under the control of a PRT synchronization pulse, and output the system function to the echo generation module 8;
the echo generating module 8 is used for convolving the radar pulse transmitting signal with a system function, generating an original echo signal under the control of a PRT synchronous pulse, and outputting a result to the delay and range gate gating module 9;
and the time delay and range gate gating module 9 is used for outputting echo data in a 'streaming' mode under the control of the PRT synchronous pulse, so that the radar echo signal stream is better simulated.
Referring to fig. 3, the process of simulating a radar echo signal using the system of the present invention is as follows:
in the process 1, the PRT synchronization module 1 reads radar working frequency, sampling frequency, emission signal pulse width, emission signal bandwidth, pulse repetition time and emission pulse number in various modes from the simulation interactive interface, outputs PRT synchronization signals PRTTrigger to each module, and outputs a marker FrameTrigger of the radar working mode to the beam center calculation module.
And 2, receiving PRT synchronous pulse input by the radar pulse signal transmitting module 2, detecting a rising edge, generating a baseband linear frequency modulation signal LFM when the rising edge comes, and outputting N according to a pulse repetition period PRT a Sending the LFM pulse signals to an echo generating module 8;
the expression for generating the baseband chirp signal LFM is as follows:
wherein the content of the first and second substances,
t is the fast time, t m Is a slow time, T p For the pulse width of the LFM signal, f c Is the carrier frequency, k r Frequency is adjusted linearly;
and 3, receiving PRT synchronous pulse input trigger by the radar track leading-in module 3, detecting a rising edge, reading track file information when the rising edge comes, updating a radar position, a speed, an acceleration, an incidence angle, a speed vector, a beam projection included angle on the ground and radar position information with errors, and outputting a result to the beam center calculating module 4, the sub-scene intercepting module 5, the antenna directional diagram calculating module 6 and the system function calculating module 7.
And 4, receiving PRT synchronous pulse input trigger by the beam center calculation module 4, detecting a rising edge, inputting the position, the speed, the target position coordinate and the current mode marking bit of the radar when the rising edge comes, calculating the pixel point coordinate information of the beam center in the scene graph, and outputting the result to the sub-scene interception module 5, the antenna directional diagram calculation module 6 and the system function calculation module 7.
And 5, inputting a beam center coordinate and a radar position coordinate by the sub-scene intercepting module 5 according to different beam irradiation modes of different radar modes, synchronously triggering by the PRT, intercepting an imaged sub-scene from a guided large scene, and outputting data to the antenna directional diagram calculating module 6 and the system function calculating module 7.
And 6, inputting radar position information, pixel point coordinate information of a beam center in a large scene, sub-scene size, receiving beam bias angle, distance direction and azimuth direction resolution and antenna beam width by the antenna directional pattern calculation module 6, receiving PRT synchronous pulse input trigger, detecting a rising edge, calculating antenna gain of each point in the sub-scene when the rising edge comes, and outputting a result to the system function calculation module 7.
The calculation of the antenna gain of each point in the sub-scene has two modes:
the first method is to use the directional diagram of a single antenna and the bias angle of the receiving antenna to be 0 for calculation in an SAR imaging mode, namely, firstly, the azimuth angle alpha of the beam center is obtained according to the input radar position information and the coordinates of the beam center on the central pixel point of the scene c And a pitch angle beta c Azimuth angle alpha of target relative to radar RT And a pitch angle beta RT (ii) a Then calculating the antenna gain rcs modulated by the directional diagram of the transmitting antenna 1 ：
rcs 1 ＝abs((sinc(α RT -α c ))*(sinc(β RT -β c ))) <2>
Wherein abs is a function of absolute value, which represents multiplication,
the second is that under the single-pulse mode, four single antennas are used, and corresponding receiving offset angles are set to form a four-antenna directional diagram for calculation,firstly, obtaining an azimuth angle alpha of a beam center according to input radar position information, coordinates of the beam center at a pixel point at the scene center, an azimuth deviation angle +/-Delta alpha of a sub-beam relative sum beam and a pitching deviation angle +/-Delta beta c And a pitch angle beta c Azimuth angle alpha of target relative to radar RT And a pitch angle β RT (ii) a Respectively calculating the antenna gain rcs modulated by the directional diagrams of the transmitting antennas with different offsets 21 、rcs 22 、rcs 23 、rcs 24 ：
rcs 21 ＝abs((sinc(α RT +Δα-α c ))*(sinc(β RT +Δβ-β c ))) <3>
rcs 22 ＝abs((sinc(α RT +Δα-α c ))*(sinc(β RT -Δβ-β c ))) <4>
rcs 23 ＝abs((sinc(α RT -Δα-α c ))*(sinc(β RT +Δβ-β c ))) <5>
rcs 24 ＝abs((sinc(α RT -Δα-α c ))*(sinc(β RT -Δβ-β c ))) <6>
Total antenna gain rcs 2 Comprises the following steps:
rcs 2 ＝rcs 21 *rcs 22 *rcs 23 *rcs 24 <7>
wherein, rcs 21 The antenna gain, rcs, is the azimuth deviation angle of the sub-beam relative to the sum beam, Δ α, and the elevation deviation angle, Δ β 22 The antenna gain with the sub-beams having an azimuth deviation angle of delta alpha and a pitch deviation angle of-delta beta with respect to the sum beam, rcs 23 The antenna gain with the sub-beams having an azimuth deviation angle of- Δ α and a pitch deviation angle of Δ β with respect to the sum beam, rcs 24 The antenna gain is obtained by the antenna with the sub-beam relative to the beam and the azimuth deviation angle of-delta alpha and the pitching deviation angle of-delta beta.
A process 7, inputting the position of the radar, the beam center coordinate, the sub-scene data, the antenna directional diagram and the initial distance gate information of the echo by the system function calculation module 7, receiving the input trigger of the PRT synchronous pulse, detecting the rising edge, calculating the system function of the echo by using a concentric circle algorithm when the rising edge comes, and outputting the result to the echo generation module 8;
the principle of the concentric circle algorithm is as follows:
without considering the wavefront curvature, taking the positive side as an example, the distribution of echoes of point targets on a two-dimensional plane like a rectangular array shown in the first graph of fig. 4 will become the result of the second graph after compressing the distance to the pulse, which is caused by the curvature generated by the radar motion. After correcting the overbending, the energy distribution of the point target will be in the same range unit, as shown in the third graph in fig. 4, and then the target can be subjected to azimuth imaging along the azimuth direction, and the imaging result of the point target is obtained, as shown in the last graph in fig. 4.
It can be seen that different point targets, which are located at different distances from the radar, are distributed in different range units, because the distances from different target points to the radar platform are different, which results in different delay times, and the distance sampling frequency f is used s Sampling the distance delay of the target point, wherein the interval of the sampling units is c/2f s C is the light speed, and the echo distribution of each target point is distributed according to the integral multiple relation of the interval of the sampling units; for reaching target points with the same radar range in a scene at different azimuth moments, the integral multiple relation of the sampling units is the same, so that the complex echoes of the target points are accumulated in the same range unit.
Therefore, in order to quickly obtain an echo system function so as to meet the requirement of real-time echo signal generation and maintain the calculation accuracy of the echo signal, points with the same length from a radar are considered to be located on the same distance unit, the points in a scene are accumulated along a concentric circle with the radar as the center of a circle to obtain a one-dimensional range profile of the radar, and then the FFT is used for quickly realizing the generation of the echo system function in a frequency domain, so that multiple points can be simultaneously processed, and the calculation amount is reduced.
According to the above thought, the specific method for calculating the echo system function is as follows:
at each azimuth moment, firstly, the distances R (k) from all point targets in the scene to the radar are calculated, and the distances are compared with the distance sampling unit to obtain the distribution situation of the points on all concentric circles, namely the distribution situation of the points on all concentric circles
Formula (II)<8>In, n k Indicating the position of the distance elements, i.e. the points are distributed on several concentric circles, delta r Is the size of a distance unit, and
as shown in fig. 5, after the concentric circles are distributed, P point targets are shared on a certain concentric circle within the beam irradiation range, according to the formula<9>It can be known that the P scattering points should be distributed in the same range unit, and they can uniformly generate echo signals, and because the azimuth phase information of the echo signals is more sensitive than the range envelope information, the integrity of the azimuth phase information is ensured, i.e. the P scattering points cannot be used as the range envelope<9>The distance approximation calculation is performed so that the azimuth phase signal s (mT; R) of each point B ) Independent calculation is required:
where σ is the gray value of the point target, mT is t m In discrete form of (A), R B Is the closest distance of the radar to the target, λ is the radar operating wavelength, R (mT; R) is the distance of the radar to the scattering point at mT, exp represents an exponential function.
Then, the pair formula<10>Summing to obtain the same concentric circleUpper point target azimuth phase signal s 2 ：
Wherein σ i The gray value of the ith point target on the same concentric circle.
In pair type<11>Summing to obtain all the distance unit data s at the whole moment 3 ：
Formula (II)<12>Where δ is the impulse response function, k denotes the distance falling within the several distance units, P n Indicating the number of point targets in the nth range bin. The number of ground scattering points in different distance units is different, and the difficulty of the accumulation process is different. In pair type<12>Fourier transform FFT is carried out to transform the frequency domain to a frequency modulation item of multiplying the frequency domain by the distance direction, and then inverse Fourier transform IFFT is utilized to transform the frequency domain back to the time domain, so that a system function s of an echo can be obtained 4 (k,mT；R B )：
Wherein f is r Indicating the range-wise frequency.
And 8, inputting a system function and a baseband linear frequency modulation signal LFM by the echo generating module 8, performing convolution operation to generate an original echo signal, receiving PRT synchronous pulse input trigger, detecting a rising edge, and outputting a result to the delay and range gate gating module 9 when the rising edge comes.
In the process 9, the delay and range gate gating module 9 outputs the input original echo signal in a 'stream' manner under the control of the PRT synchronization pulse.
In summary, the simulink-based high-speed platform radar echo signal simulation system has the advantages of being modular and hierarchical in design, being capable of conveniently and flexibly simulating radar echo signals under different radar working modes, different external environments and different target scenes, helping scientific researchers get rid of the limitation of radar equipment conditions, not needing to rely on expensive radar equipment to obtain relevant radar echo data, and being more efficient and convenient than the traditional actual measurement data mode.
While the invention has been described in connection with what is presently considered to be the most practical and preferred embodiment, it is to be understood that the invention is not to be limited to the disclosed embodiment, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.
Claims (4)
1. High-speed platform radar return signal analog system based on simulink generates different modules through simulink, its characterized in that, the module that generates includes:
the PRT synchronization module (1) is used for finishing the generation of synchronization pulses, outputting the synchronization pulses to each module and synchronizing the clock of the whole system;
the radar pulse signal transmitting module (2) is used for receiving PRT synchronous pulse input, generating a radar pulse transmitting signal and outputting the radar pulse transmitting signal to the echo generating module (8);
the radar track importing module (3) is used for reading radar track file information under the control of PRT synchronous pulse, and respectively outputting the read radar track information to the beam center computing module (4), the sub-scene intercepting module (5), the antenna directional diagram computing module (6) and the system function computing module (7);
the beam center calculating module (4) is used for updating pixel point coordinate information of a beam center in a scene graph according to the input radar and the current mode marking bit under the control of a PRT synchronous pulse, and outputting the beam center coordinate information to the sub-scene intercepting module (5), the antenna directional diagram calculating module (6) and the system function calculating module (7);
the sub-scene intercepting module (5) is used for intercepting an imaged sub-scene from a guided large scene according to different beam irradiation modes of different modes of the radar, input beam center coordinates and radar track information under the control of the PRT synchronous pulse, and outputting data of the sub-scene to the antenna directional diagram calculating module (6) and the system function calculating module (7);
an antenna directional pattern calculation module (6) for calculating the antenna gain of each point in the sub-scene according to the input radar track information, the beam center coordinate information, the sub-scene size, the receiving beam bias angle, the distance direction and azimuth direction resolution and the antenna beam width under the control of the PRT synchronous pulse, and outputting the result to a system function calculation module (7);
the system function calculation module (7) is used for calculating a system function of an echo by utilizing a concentric circle algorithm according to an input radar position, a beam center coordinate, sub-scene data, an antenna directional diagram and initial range gate information of the echo under the control of a PRT synchronous pulse, and outputting the system function to the echo generation module (8);
the echo generating module (8) is used for convolving the radar pulse transmitting signal with a system function, generating an original echo signal under the control of a PRT synchronous pulse, and outputting a result to the delay and range gate gating module (9);
and the time delay and range gate gating module (9) is used for outputting echo data in a 'streaming' mode under the control of the PRT synchronous pulse, so that the radar echo signal stream is better simulated.
2. The system of claim 1, wherein different modules are generated by simulink, which is implemented as follows:
defining basic parameters required by a seeker system through an M script, constructing a model of each module in a simulink, and setting input and output ports to realize the functions of each module;
calling an S function to package each module and generating a parameter setting interface;
and connecting corresponding input and output ports in each module.
3. The system according to claim 1, characterized in that the antenna pattern calculation module (6) calculates the antenna gain for each point in the sub-scene, which is implemented as follows:
firstly, according to the input radar position information, the beam center coordinate information, the receiving beam offset angle and the antenna beam width, the azimuth angle alpha of the beam center is obtained c And a pitch angle β c Azimuth angle alpha of target relative to radar RT And a pitch angle β RT ；
And then calculating the antenna gain rcs after the directional diagram modulation of the transmitting antenna:
rcs＝abs((sinc(α RT -α c ))*(sinc(β RT -β c )))， <1>
wherein abs is a function of absolute value, which means multiplication,
4. the system according to claim 1, characterized in that the system function calculation module (7) calculates the system function of the echo using a concentric circle algorithm, which is implemented as follows:
calculating the distances R (k) from all point targets to the radar in the scene, and comparing the distances with a distance sampling unit to obtain the distribution condition of points of all concentric circles, namely obtaining P point targets on a certain concentric circle in total in the beam irradiation range;
distributing the P scattering points in the same range unit, and calculating azimuth phase signal s of each point 1 (mT；R B ):
Wherein σ is the gray value of the point target; mT is the slow time t m Of discrete form of (A), R B Is the closest distance of the radar to the scattering point, λ is the radar operating wavelength, R (mT; R) B ) Is the distance of the radar to the scattering point at time mT;
by adding point targets on the same concentric circle, i.e. pair<2>Summing to obtain a distance sheetData s of elements containing azimuth phase 2 ：
In pair type<3>Summing to obtain all the distance unit data s at the whole moment 3 ：
Wherein, δ is an impact response function, and k is within a few distance units;
in pair type<4>Fourier transform FFT is carried out, the FFT is converted into a frequency domain multiplied by a frequency modulation item in the distance direction, and then inverse Fourier transform IFFT is utilized to convert the frequency domain back into a time domain, so that a system function s of an echo can be obtained 4 (k,mT；R B )：
Wherein k is r Indicating the range chirp, f r Indicating the range-wise 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 CN109471080A (en) | 2019-03-15 |
| CN109471080B true 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) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901165B (en) * | 2019-03-28 | 2023-03-14 | 河南九乾电子科技有限公司 | Satellite-borne SAR echo simulation device and method |
| CN109782241B (en) * | 2019-03-29 | 2021-01-29 | 北京润科通用技术有限公司 | Meteorological radar echo simulation method and system |
| CN110988858B (en) * | 2019-11-11 | 2021-12-07 | 西安空间无线电技术研究所 | 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 |
| CN111624564B (en) * | 2020-05-27 | 2023-06-16 | 北京润科通用技术有限公司 | Radar pitching angle target simulation system and method |
| CN112213721B (en) * | 2020-09-16 | 2023-10-03 | 西安科技大学 | Millimeter wave three-dimensional imaging method for scanning external or internal cylindrical scene facing security inspection |
| CN114442051B (en) * | 2020-11-05 | 2024-05-24 | 北京华航无线电测量研究所 | 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 |
| CN112698280B (en) * | 2020-12-09 | 2024-05-31 | 南京长峰航天电子科技有限公司 | Double-base SAR real-time echo simulation method based on DSP and FPGA architecture |
| CN114325606B (en) * | 2021-11-17 | 2024-06-25 | 西安电子科技大学 | Multi-system agile radar radio frequency echo signal simulation method |
| CN114636983B (en) * | 2022-03-21 | 2024-06-11 | 中国电子科技集团公司第三十八研究所 | SAR target scene generation and simulation device |
| CN114966585B (en) * | 2022-04-12 | 2024-06-25 | 西安电子科技大学 | Radar echo generation method for simulating high-speed motion by adopting near-field low-speed motion |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006022814A1 (en) * | 2006-05-13 | 2007-11-15 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | High-resolution Synthetic Aperture Side View Radar System using Digital Beamforming |
| CN104090277B (en) * | 2014-07-21 | 2016-09-28 | 中国科学院电子学研究所 | A kind of slidingtype circumferential synthetic aperture radar formation method |
| CN107765226B (en) * | 2017-09-18 | 2020-06-09 | 北京空间飞行器总体设计部 | SAR satellite radar echo simulation method, system and medium |
-
2018
- 2018-11-09 CN CN201811328332.7A patent/CN109471080B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109471080A (en) | 2019-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109471080B (en) | High-speed platform radar echo signal simulation system based on simulink | |
| CN109239684B (en) | Radar target echo simulation system based on electromagnetic scattering model and simulation method thereof | |
| CN104007440B (en) | One accelerated decomposition rear orientation projection spot beam SAR formation method | |
| CN101526614B (en) | SAR echo rapid simulation method based on sub-aperture and equivalent scatterer | |
| CN101907704B (en) | Method for evaluating simulation imaging of multi-mode synthetic aperture radar | |
| CN103487802B (en) | Scanning radar angle super-resolution imaging method | |
| CN107765226B (en) | SAR satellite radar echo simulation method, system and medium | |
| CN103941243B (en) | A kind of spinning aircraft based on SAR three-dimensional imaging surveys high method | |
| CN107479038B (en) | High-precision radar target echo real-time simulation method | |
| CN109164428A (en) | Digital radar analogue system and method | |
| CN104614713A (en) | Radar echo signal simulator suitable for onboard radar system | |
| CN102798861B (en) | Bistatic synthetic aperture radar (SAR) imaging method based on optimal image space | |
| CN110208760A (en) | A kind of radar return emulation mode based on time domain up-sampling | |
| Ren et al. | 3D Imaging Algorithm for Down‐Looking MIMO Array SAR Based on Bayesian Compressive Sensing | |
| CN107907870B (en) | Signal generation method for verifying intersection butt joint microwave radar angle measurement function | |
| CN112649806B (en) | MIMO radar near-field three-dimensional imaging method | |
| CN116381629A (en) | Radar large-scale target simulation system and method based on real-time dynamic convolution | |
| Shu et al. | Development of multichannel real-time Hardware-in-the-Loop radar environment simulator for missile-borne Synthetic Aperture Radar | |
| CN111175734A (en) | Radiation source positioning and speed measuring method and system based on virtual equivalent synthetic aperture | |
| CN113447896B (en) | Undulating terrain SAR echo simulation method based on dynamic occlusion judgment | |
| Wang et al. | Study on motion compensation method for W-band UAV MISAR real-time imaging | |
| 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 | |
| CN117148351B (en) | Missile-borne SAR image imaging method and device based on satellite SAR image | |
| CN114636983B (en) | SAR target scene generation and simulation device |
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 |