CN111289962A - Multi-factor ship target radar echo simulation method - Google Patents

Multi-factor ship target radar echo simulation method Download PDF

Info

Publication number
CN111289962A
CN111289962A CN202010127454.0A CN202010127454A CN111289962A CN 111289962 A CN111289962 A CN 111289962A CN 202010127454 A CN202010127454 A CN 202010127454A CN 111289962 A CN111289962 A CN 111289962A
Authority
CN
China
Prior art keywords
echo
ship
radar
target
intensity
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
CN202010127454.0A
Other languages
Chinese (zh)
Other versions
CN111289962B (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.)
China Shipbuilding Pengli Nanjing Atmospheric And Ocean Information System Co ltd
Cssc Pengli Nanjing Technology Group Co ltd
Original Assignee
Csic Nanjing Pride Technology Group Co ltd
CSIC Pride Nanjing Atmospheric and Oceanic Information System Co 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 Csic Nanjing Pride Technology Group Co ltd, CSIC Pride Nanjing Atmospheric and Oceanic Information System Co Ltd filed Critical Csic Nanjing Pride Technology Group Co ltd
Priority to CN202010127454.0A priority Critical patent/CN111289962B/en
Publication of CN111289962A publication Critical patent/CN111289962A/en
Application granted granted Critical
Publication of CN111289962B publication Critical patent/CN111289962B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/66Radar-tracking systems; Analogous systems
    • 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/4052Means for monitoring or calibrating by simulation of echoes
    • 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/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418Theoretical aspects

Landscapes

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

Abstract

The invention discloses a multi-factor ship target radar echo simulation method, which comprises the following steps of 1, setting target quantity and 6 factor simulation of each target: center coordinates (xn, yn), ship length Ln, width Wn, echo intensity An, motion direction d (t), and motion velocity v (t); step 2, adding a central coordinate random error Sn; step 3, adding a random error Wn of the motion speed; step 4, establishing a correction motion curve; step 5, coordinate conversion; step 6, adding an intensity random error Qn; and 7, solving the real echo intensity Tn. The invention introduces 3 random errors of elements, so that the simulation method can simulate the radar echo of a single ship or a plurality of ships in a motion state as truly as possible, can be used for simulating various complex ship motion states during the research and test of a ship radar target tracking algorithm, greatly reduces the intensity of radar outfield test and ship running test, and greatly reduces the development cost.

Description

Multi-factor ship target radar echo simulation method
Technical Field
The invention relates to the field of radar simulation methods, in particular to a multi-factor ship target radar echo simulation method.
Background
The research and test of ship radar target tracking algorithm usually need to organize one or more ships, complete complex actions such as crossing, turning, pursuing, meeting and the like under various water area environments according to requirements, erect radars on the shore, scan and record radar echoes of ship targets through long-time radars, in order to ensure the richness of obtained data and the reliability of data samples, the actual measurement often needs a large amount of manpower and material resources for several months or several years. At present, no effective radar echo simulator can effectively and more truly simulate ship radar echoes.
The existing simulation methods are all simulated by a rectangular or polygonal color block with uniform color distribution. The method actually ignores the most important factors (environmental factors) in radar echo, and when the echo is used for testing a tracking algorithm, the information such as the coordinates, the speed and the like of the target is accurate and error-free. In a real situation, the echo shape and intensity distribution of a ship are distributed with strong center and weak periphery, and are influenced by the environment, and the distribution cannot be accurately described by a numerical analysis method, so that in fact, information such as coordinates and speed of the ship is continuously changed, and the change has certain randomness. In more complex conditions, when multiple ships intersect, turn, overtake, meet and the like, and even collide, the shapes of radar echoes of the ships are changed more complexly, and the conditions are usually difficult to test through real ships, so that the algorithm research of ship target tracking is greatly limited by a measuring method.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a multi-element ship target radar echo simulation method aiming at the defects of the prior art, wherein the multi-element ship target radar echo simulation method simulates the complex ship motion state through software, combines the radar principle and some characteristics of radar echo existing in actual detection, and provides real approaching test data for the research of a radar tracking algorithm.
In order to solve the technical problems, the invention adopts the technical scheme that:
a multi-element ship target radar echo simulation method comprises the following steps.
Step 1, setting target number and simulation of 6 elements of each target: the number n of target ships to be tracked is set, and 6 key elements are set for radar echoes of each target ship. Wherein, 6 key elements are respectively: center coordinates (xn, yn), ship length Ln, width Wn, echo intensity An, motion direction d (t), and motion velocity v (t). Then, the set 6 key elements are simulated in a rectangular coordinate system. Where n =1,2, … …, t denotes the movement time.
Step 2, adding a center coordinate random error Sn: and (3) adding the center coordinate random error Sn to the center coordinates (Xn, Yn) simulated in the step 1 to form corrected center coordinates (Xn, Yn), wherein Xn = Xn Sn and Yn = Yn Sn.
Step 3, adding a random error Wn of the motion speed: and (4) adding the random error Wn of the motion speed to the motion speed V (t) simulated in the step 1 to form a corrected motion speed Wn V (t).
Step 4, establishing a correction motion curve: and establishing a corrected motion curve of the reaction environment change factor according to the corrected central coordinates (Xn, Yn) in the step 2 and the corrected motion speed Wn V (t) in the step 3.
And 5, coordinate conversion: and (3) converting the rectangular coordinate system where the echo intensity An in the step (1) is located into a polar coordinate system.
Step 6, adding intensity random error Qn: for echo intensity An in a polar coordinate system, the echo intensity An has radial echo intensity An1 and tangential echo intensity An 2; adding a radial echo intensity error Qn1 to the radial echo intensity An1 to form a corrected radial echo intensity Qn1 An 1; adding a tangential echo intensity error Qn2 to the tangential echo intensity An2 to form a corrected radial echo intensity Qn2 An 2; the radial echo intensity error is used for simulating the influence of a range side lobe on the radar echo shape when two coordinate radars are detected; the tangential echo intensity error is used for simulating the influence of the azimuth side lobe on the radar echo shape when the two-coordinate radar is detected.
And 7, solving the real echo intensity Tn: the vector sum of the modified radial echo intensities Qn1 An1 and the modified radial echo intensities Qn2 An2 constitutes the true echo intensity Tn.
In step 2, the center coordinate random error Sn = Rand (min 1, max 1) x K1, wherein K1 represents a random influence reduction coefficient; min1 represents the minimum random error of the center coordinate, and max1 represents the maximum random error of the center coordinate.
K1 is related to the distance between the target ship and the radar station and the shape of the ship, and the closer the target ship is to the radar station, the smaller K1 is; the simpler the ship shape, the smaller K1; min1= -Ln/2, max1= Ln/2.
In step 3, the random error Wn =1+ Rand (min 2, max 2) x K2 of the movement speed, where K2 represents the random influence reduction coefficient, min2 represents the minimum random error of the movement speed, and max2 represents the maximum random error of the movement speed.
K2 is related to the degree to which the speed of the target vessel is influenced by the environment, the greater the influence degree, the greater K2; min2= -1, max2= 1.
And selecting a corresponding distribution function type according to the ship length Ln, the width Wn and the ship type simulated in the step 1 by the radial echo intensity error Qn1 and the tangential echo intensity error Qn 2. Wherein the distribution function types include normal distribution, positive bias distribution, negative bias distribution, multimodal distribution, and chi-squared distribution.
The invention has the following beneficial effects:
1. the complex ship motion state is simulated by software, and the radar principle and some characteristics of radar echo in actual detection are combined to provide real-approaching test data for the research of a radar tracking algorithm.
2. The random error Sn of the central coordinate, the random error Wn of the movement speed and the random error Qn of the echo intensity are adopted to simulate the detection errors of a real radar when the radar is influenced by factors such as the detection electromagnetic environment of the radar, the storm of the water surface, the weather condition and the like. Meanwhile, the three random factors are represented by three random numbers with different value ranges. The selection of the value range indicates the severity of the environmental factors under the simulation condition. The noise environment is indicated to be good or bad compared to the radar.
Drawings
Fig. 1 shows a flow chart of a multi-element ship target radar echo simulation method of the invention.
Fig. 2 shows the effect comparison of the influence of different radial echo intensity distribution functions and tangential echo intensity distribution functions on the echo intensity.
Figure 3 shows the simulated echo effect in the case of a radial, cut-away intensity distribution.
Figure 4 shows the effect of radar echo simulation on a moving target.
Fig. 5 shows a simulation diagram of 6 key elements in the overtaking state.
Fig. 6 shows a schematic diagram when the distribution function type is normal distribution.
Fig. 7 shows a schematic diagram of the distribution function type with a negative bias distribution.
Fig. 8 shows a schematic diagram of the distribution function type as a forward biased distribution.
Fig. 9 shows a simulation diagram of 6 key elements in the cross state.
Fig. 10 shows a simulation diagram of 6 key elements in the circling state.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific preferred embodiments.
In the description of the present invention, it is to be understood that the terms "left side", "right side", "upper part", "lower part", etc., indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and that "first", "second", etc., do not represent an important degree of the component parts, and thus are not to be construed as limiting the present invention. The specific dimensions used in the present example are only for illustrating the technical solution and do not limit the scope of protection of the present invention.
As shown in fig. 1, a multi-element ship target radar echo simulation method includes the following steps.
Step 1, setting target number and simulation of 6 elements of each target: the number n of target ships to be tracked is set, and 6 key elements are set for radar echoes of each target ship. Wherein, 6 key elements are respectively: center coordinates (xn, yn), ship length Ln, width Wn, echo intensity An, motion direction d (t), and motion velocity v (t).
The invention uses rectangle to simulate the ship echo, so the length Ln and the width Wn of the target ship become two important factors. Then, the set 6 key elements are simulated in a rectangular coordinate system. Where n =1,2, … … denotes a ship number, and t denotes a movement time.
The motion direction and the motion speed are all functions changing along with time, and can be a regular curve or any motion track customized by a user, but it is emphasized that the motion direction and the motion speed are smooth curves in consideration of the reality of the ship motion; v (t) and D (t) together form the motion curve of the ship target.
And 2, adding a central coordinate random error Sn.
Any ship, even if anchored and parked, continuously changes in the scanning echo of each circle of the radar, because the shape of the echo obtained by each circle of the radar is different, the coordinates of a target are calculated when the radar echo is processed, each circle slightly changes, and in order to describe the change, a central coordinate random error is introduced.
After the central coordinate is influenced by random factors, the motion direction after calculation naturally changes randomly, so that the random factors of the motion direction are not considered. The length and width of the target are closely related to the intensity distribution function in the next calculation and are not considered here.
And (3) adding the center coordinate random error Sn to the center coordinates (Xn, Yn) simulated in the step 1 to form corrected center coordinates (Xn, Yn), wherein Xn = Xn Sn and Yn = Yn Sn.
The random error of the center coordinates Sn = Rand (min 1, max 1) K1, wherein K1 represents a random influence reduction coefficient, K1 is related to the distance of a target ship from a radar station and the shape of the ship, and the closer the target ship is to the radar station, the smaller K1 is; the simpler the ship shape, the smaller K1; here, it is preferably 0.1; min1 represents the minimum random error of the center coordinates, max1 represents the maximum random error of the center coordinates, min1= -Ln/2, max1= Ln/2.
Since the change in direction is calculated from the change in coordinates, the direction of motion d (t) is expressed by the following equation:
D(t)=arctan[(Yn(t+2)-Yn(t+1))/(Xn(t+2)-Xn(t+1))]-arctan[(Yn(t+1)-Yn)/(Xn(t+1)-Xn)])
in the above formula, if Xn and Yn are the corrected central coordinates of the nth target ship at time t, Xn (t +1) and Yn (t +1) are the corrected central coordinates of the nth target ship at time t + 1; xn (t +2) and Yn (t +2) are the corrected center coordinates of the nth target vessel at time t + 2.
Step 3, adding a random error Wn of the motion speed: the moving speed of the ship is easily influenced by the environment such as water flow, wind speed and the like, therefore, the real speed of the ship with fixed speed is always changed along with the conditions such as water speed, wind speed and the like, and in order to describe the change, a random speed error Wn is introduced.
And (4) adding the random error Wn of the motion speed to the motion speed V (t) simulated in the step 1 to form a corrected motion speed Wn V (t).
The above-mentioned random error Wn =1+ Rand (min 2, max 2) × K2 of the movement speed, where min2 represents the minimum random error of the movement speed and max2 represents the maximum random error of the movement speed. K2 represents a random influence reduction factor, K2 is related to the degree to which the speed of the target vessel is influenced by the environment, the greater the influence degree, the greater K2. In general, the speed influence of the environment on ships varies from ship to ship, and for a yacht, K2 is smaller and may be set to 0.1 or less, and for a kayak, it may be set to 0.5 or more.
The random motion speed error Wn is usually between 0 and 2, so min2= -1, and max2= 1.
Step 4, establishing a correction motion curve: and establishing a corrected motion curve of the reaction environment change factor according to the corrected central coordinates (Xn, Yn) in the step 2 and the corrected motion speed Wn V (t) in the step 3.
And 5, coordinate conversion: and (3) converting the rectangular coordinate system where the echo intensity An in the step (1) is located into a polar coordinate system.
Step 6, adding intensity random error Qn: for echo intensities An in a polar coordinate system, there are radial echo intensities An1 and tangential echo intensities An 2.
The intensity of each real radar echo is not invariable all the time, and random errors Qn of the target intensity are introduced along with the influences of the change of the equivalent reflection area of a moving target, the environmental change and the like.
The radial echo intensity error Qn1 is added to the radial echo intensity An1 to form a modified radial echo intensity Qn1 An 1. The tangential echo intensity error Qn2 is added to the tangential echo intensity An2 to form a modified radial echo intensity Qn2 An 2. The radial echo intensity error is used for simulating the influence of a range side lobe on the radar echo shape during the detection of the two-coordinate radar. The tangential echo intensity error is used for simulating the influence of the azimuth side lobe on the radar echo shape when the two-coordinate radar is detected.
And 7, solving the real echo intensity Tn: the vector sum of the modified radial echo intensities Qn1 An1 and the modified radial echo intensities Qn2 An2 constitutes the true echo intensity Tn.
And selecting a corresponding distribution function type according to the ship length Ln, the width Wn and the ship type simulated in the step 1 by the radial echo intensity error Qn1 and the tangential echo intensity error Qn 2. Wherein the distribution function types include normal distribution, positive bias distribution, negative bias distribution, multimodal distribution, and chi-squared distribution. In fig. 2, the left side shows four intensity profiles, and the right side shows the effect of the echo intensity profile corresponding to the four intensity profiles. The intensity distribution function can simulate a real ship as much as possible, for example, the 4 th component distribution function in fig. 2 can simulate the situation that the echo of the bow and the stern is weak, and two piles of goods with strong reflection are loaded on the ship (such as an LNG ship, a container ship and a special operation ship provided with a crane).
FIG. 3 shows a ship rectangular radar echo in which the radial echo intensity distribution function, the tangential echo intensity distribution function and the random error of the echo target intensity are considered simultaneously, and a rectangular radar echo in which any random factor and intensity distribution are not considered completely. The comparison of the two figures shows that the echo simulation effect of the device is closer to reality.
As shown in fig. 4, after the motion function of the target is combined, the multi-element ship target echo simulation method provided by the invention provides a more real radar video echo for the research of a ship target tracking algorithm.
The three random factors Sn, Wn and Qn are represented by three random numbers with different value ranges. The selection of the value range indicates the severity of the environmental factors under the simulation condition. The noise environment is indicated to be good or bad compared to the radar.
The present invention will be described in detail with reference to the following, meeting, crossing and circling states as specific embodiments.
1. Overtaking meeting state
The two states are similar, for example, the tracking is taken as an example, the state simulates that the A ship is in front of the B ship and is 200 meters away from the front and the back, the A ship moves in the same direction at a speed higher than that of the B ship, the A ship is required to be within 15 meters away from the B ship when the A ship passes through the B ship, and the simulation is finished after the A ship passes through the B ship by 200 meters, as shown in FIG. 5.
Let A, B be yachts 5m long and 2 m wide, with the other parameters shown in FIG. 5. Initially, the motion curves of the two vessels can be expressed as:
YA=15;XA=100+10t;
YB=30;XB=300+5t;
considering the random error of the motion direction detection, namely the random error Sn of the central coordinate, taking Rand (-2.5, 2.5) as random data in the range of (-2.5, 2.5), considering that the speed boat is small, the influence of the random error factor in the direction can be reduced, and the random influence is reduced to 0.1, so that the motion curves of the two boats can be expressed as follows:
YA=15+ Rand(-2.5,2.5)*0.1;XA=100+ Rand(-2.5,2.5)*0.1+10t;
YB=30+ Rand(-2.5,2.5)*0.1;XB=300+ Rand(-2.5,2.5)*0.1+5t。
considering the random error of the motion speed detection and then introducing the random error Wn of the motion speed, the motion curves of the two ships can be expressed as follows:
YA=15+ Rand(-2.5,2.5)*0.1;XA=100+ Rand(-2.5,2.5)*0.1+10*(1+ Rand(-1,1)*0.1)*t;
YB=30+ Rand(-2.5,2.5)*0.1;XB=300+ Rand(-2.5,2.5)*0.1+5*(1+ Rand(-1,1)*0.1)*t。
after the coordinate conversion is completed, the influence of the echo intensity distribution on the radial direction and the tangential direction is considered, and since the influence of the change of the echo in the Y direction of the radar on the tracking when two ships approach about 15 meters is concerned under the test condition, the radial intensity distribution function Qn1 can be ignored, and only the tangential intensity distribution function Qn2 is considered. Due to bilateral symmetry of the yacht, a standard normal distribution function can be selected as a tangential intensity distribution function. (there are many alternative distribution functions including, but not limited to, normal distribution, positive bias distribution, negative bias distribution, multi-peak distribution, chi-square distribution, etc., where the distribution function morphology can be changed by adjusting parameters, for example, a normal distribution, which can be selected to be "flat" for sand carriers, and "steep" for LNG). The preferred distribution is as shown in fig. 6-8.
After a rectangular ship is processed by an intensity distribution function in the length direction and the width direction or the radial direction and the tangential direction after coordinate conversion, the echo shape of the rectangular ship is not rectangular any more after processing but is closer to the echo of a real ship because the intensity distribution function is selected according to different ship shapes.
2. Cross state
The state A ship and the state B ship run on a 90-degree vertical route. Sailing in the south and east directions, respectively, at the same speed. Requiring a distance of less than 15 meters at the crossover.
Suppose that ship a is a river crossing yacht with the length of 5 meters and the width of 2 meters, ship B is a cargo ship in a river channel (the middle section of the front section of the ship is lower, and the cockpit is higher at the back of the ship), and other parameters are as shown in fig. 9. Initially, the motion curves of the two vessels can be expressed as:
XA=190;YA=200-10t;
YB=100;XB=100+10t;
considering the random error of the motion direction detection, the direction stability of the yacht is good, the influence of the random error factor (namely, the random error Sn of the central coordinate) in the direction can be reduced, the random influence is reduced to 0.1, and then the motion curves of the two ships can be expressed as follows:
XA=190+ Rand(-2.5,2.5)*0.1;YA=200+ Rand(-2.5,2.5)*0.1-10t;
YB=100+ Rand(-2.5,2.5)*0.1;XB=100+ Rand(-2.5,2.5)*0.1+10t;
considering the random error of the motion speed detection and then introducing the random error Wn of the motion speed, the motion curves of the two ships can be expressed as follows:
XA=190+ Rand(-2.5,2.5)*0.1;YA=200+ Rand(-2.5,2.5)*0.1-10*(1+ Rand(-1,1)*0.1)t;
YB=100+ Rand(-2.5,2.5)*0.1;XB=100+ Rand(-2.5,2.5)*0.1+10*(1+ Rand(-1,1)*0.1)t;
after the coordinate conversion is completed, the influence of the echo intensity distribution on the radial direction and the tangential direction is considered, and under the test condition of intersection, the influence of the change of the radar echo at the intersection point on tracking is concerned when two ships are close to about 15 meters, so that the ship A can only consider the tangential echo intensity distribution Qn2, and the ship B can only consider the radial echo intensity distribution function Qn 1. At the moment, because the ship A is a yacht, the standard normal distribution can be selected as a tangential intensity distribution function; the B-ship is a cargo ship and a "flat" negative bias profile is selected as a function of the radial intensity profile.
3. State of revolution
A ship, a speed boat 5 meters long and 2 meters wide, takes a circle with a diameter of about 200 meters as a path and has a speed of 5 m/s.
Initially, the motion curves of the two vessels can be expressed as:
X=100*cos(5t/200/pi*360);
Y=100*sin(5t/200/pi*360);
considering the random error of direction and speed, the motion curves of the two ships can be expressed as:
X=100*cos(5t/200/pi*360)*(1+ Rand(-1,1)*0.1)+ Rand(-2.5,2.5)*0.1;
Y=100*sin(5t/200/pi*360)*(1+ Rand(-1,1)*0.1)+ Rand(-2.5,2.5)*0.1。
after the coordinate transformation is completed, the influence of the echo intensity distribution on the radial direction and the tangential direction is considered, at the moment, the distribution on the radial direction and the tangential direction is considered, and the standard normal distribution can be selected as the intensity distribution function in both directions.
Although the preferred embodiments of the present invention have been described in detail, the present invention is not limited to the details of the embodiments, and various equivalent modifications can be made within the technical spirit of the present invention, and the scope of the present invention is also within the scope of the present invention.

Claims (6)

1. A multi-element ship target radar echo simulation method is characterized by comprising the following steps: the method comprises the following steps:
step 1, setting target number and simulation of 6 elements of each target: setting the number n of target ships to be tracked, and setting 6 key elements for the radar echo of each target ship; wherein, 6 key elements are respectively: center coordinates (xn, yn), ship length Ln, width Wn, echo intensity An, motion direction d (t), and motion velocity v (t); then, simulating the set 6 key elements under a rectangular coordinate system; wherein n =1,2, … …, t denotes movement time;
step 2, adding a center coordinate random error Sn: adding a center coordinate random error Sn to the center coordinates (Xn, Yn) simulated in the step 1 to form corrected center coordinates (Xn, Yn), wherein Xn = Xn Sn and Yn = Yn Sn;
step 3, adding a random error Wn of the motion speed: adding a random error Wn of the motion speed to the motion speed V (t) simulated in the step 1 to form a corrected motion speed Wn V (t);
step 4, establishing a correction motion curve: establishing a corrected motion curve for obtaining the reaction environment change factor according to the corrected central coordinates (Xn, Yn) in the step 2 and the corrected motion speed Wn V (t) in the step 3;
and 5, coordinate conversion: converting the rectangular coordinate system where the echo intensity An in the step 1 is into a polar coordinate system;
step 6, adding intensity random error Qn: for echo intensity An in a polar coordinate system, the echo intensity An has radial echo intensity An1 and tangential echo intensity An 2; adding a radial echo intensity error Qn1 to the radial echo intensity An1 to form a corrected radial echo intensity Qn1 An 1; adding a tangential echo intensity error Qn2 to the tangential echo intensity An2 to form a corrected radial echo intensity Qn2 An 2; the radial echo intensity error is used for simulating the influence of a range side lobe on the radar echo shape when two coordinate radars are detected; the tangential echo intensity error is used for simulating the influence of a direction side lobe on the radar echo shape when the two-coordinate radar is detected;
and 7, solving the real echo intensity Tn: the vector sum of the modified radial echo intensities Qn1 An1 and the modified radial echo intensities Qn2 An2 constitutes the true echo intensity Tn.
2. The multi-element marine target radar echo simulation method of claim 1, wherein: in step 2, the center coordinate random error Sn = Rand (min 1, max 1) x K1, wherein K1 represents a random influence reduction coefficient; min1 represents the minimum random error of the center coordinate, and max1 represents the maximum random error of the center coordinate.
3. The multi-element marine target radar echo simulation method of claim 2, wherein: k1 is related to the distance between the target ship and the radar station and the shape of the ship, and the closer the target ship is to the radar station, the smaller K1 is; the simpler the ship shape, the smaller K1; min1= -Ln/2, max1= Ln/2.
4. The multi-element marine target radar echo simulation method of claim 1, wherein: in step 3, the random error Wn =1+ Rand (min 2, max 2) x K2 of the movement speed, where K2 represents the random influence reduction coefficient, min2 represents the minimum random error of the movement speed, and max2 represents the maximum random error of the movement speed.
5. The multi-element marine target radar echo simulation method of claim 4, wherein: k2 is related to the degree to which the speed of the target vessel is influenced by the environment, the greater the influence degree, the greater K2; min2= -1, max2= 1.
6. The multi-element marine target radar echo simulation method of claim 1, wherein: selecting a corresponding distribution function type according to the ship length Ln, the width Wn and the ship type simulated in the step 1 by the radial echo intensity error Qn1 and the tangential echo intensity error Qn 2; wherein the distribution function types include normal distribution, positive bias distribution, negative bias distribution, multimodal distribution, and chi-squared distribution.
CN202010127454.0A 2020-02-28 2020-02-28 Multi-factor ship target radar echo simulation method Active CN111289962B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010127454.0A CN111289962B (en) 2020-02-28 2020-02-28 Multi-factor ship target radar echo simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010127454.0A CN111289962B (en) 2020-02-28 2020-02-28 Multi-factor ship target radar echo simulation method

Publications (2)

Publication Number Publication Date
CN111289962A true CN111289962A (en) 2020-06-16
CN111289962B CN111289962B (en) 2022-02-11

Family

ID=71028350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010127454.0A Active CN111289962B (en) 2020-02-28 2020-02-28 Multi-factor ship target radar echo simulation method

Country Status (1)

Country Link
CN (1) CN111289962B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2503352A2 (en) * 2011-03-23 2012-09-26 Fujitsu Ten Limited Calculation device for radar apparatus, radar apparatus and calculation method
CN103197299A (en) * 2013-03-25 2013-07-10 南京信息工程大学 Extraction and quantitative analysis system of weather radar radial wind information
CN104237882A (en) * 2014-02-20 2014-12-24 中船重工鹏力(南京)大气海洋信息系统有限公司 Scanning system achieving multiple ocean radar scanning modes and control method
CN108051789A (en) * 2017-12-06 2018-05-18 上海无线电设备研究所 Moving-target SAR Imaging Simulations test System and method under a kind of sea background

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2503352A2 (en) * 2011-03-23 2012-09-26 Fujitsu Ten Limited Calculation device for radar apparatus, radar apparatus and calculation method
CN103197299A (en) * 2013-03-25 2013-07-10 南京信息工程大学 Extraction and quantitative analysis system of weather radar radial wind information
CN104237882A (en) * 2014-02-20 2014-12-24 中船重工鹏力(南京)大气海洋信息系统有限公司 Scanning system achieving multiple ocean radar scanning modes and control method
CN108051789A (en) * 2017-12-06 2018-05-18 上海无线电设备研究所 Moving-target SAR Imaging Simulations test System and method under a kind of sea background

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YING ZHANG等: ""Performance analysis of moving ship target signal integration and detection in geosynchronous synthetic aperture radar"", 《THE JOURNAL OF ENGINEERING》 *
朱雪慧等: ""基于动态RCS的舰船雷达回波仿真与分析"", 《雷达科学与技术》 *

Also Published As

Publication number Publication date
CN111289962B (en) 2022-02-11

Similar Documents

Publication Publication Date Title
CN111708038B (en) Unmanned ship laser radar point cloud data correction method based on attitude sensor and GNSS
CN105487061B (en) The multicharacteristic information fusion method of target data association
CN108957437B (en) Radar target tracking method based on environmental situation analysis
CN105445718B (en) A kind of DOA estimation method of the distributed overloading warship over-the-horizon radar based on array reconfiguration
CN114910915A (en) Multi-mode imaging method for underwater target of side scan sonar
CN114942428B (en) Multi-target tracking hybrid radar data association algorithm
CN106872970A (en) A kind of multiple target based on differential evolution becomes data transfer rate tracks of device and its method
CN113093159B (en) Multi-beam sounding error improved model design method
CN111289962B (en) Multi-factor ship target radar echo simulation method
CN113156418A (en) Monte Carlo simulation-based radar target tracking precision prediction method
CN109490868B (en) Offshore target motion analysis method based on distributed vertical line array
CN111090094B (en) Method and system for measuring dual-beam angle of pulse Doppler radar and storage medium
CN112731564A (en) Intelligent thunder forecasting method based on Doppler weather radar data
Brekke et al. A long simulation scenario for evaluation of multi-target tracking methods
CN112731400B (en) Method and system for estimating target vector speed of marine vessel
Chen et al. Research on AIS and radar information fusion method based on distributed Kalman
US5805101A (en) Process for determining the position of a radar target
CN113534130B (en) Multi-station radar multi-target data association method based on sight angle
CN111846144B (en) Virtual-real fusion unmanned surface vessel collision avoidance test method
Davide et al. Identification of vessel class with LSTM using kinematic features in maritime traffic control
CN115372956B (en) Hybrid system radar track starting method based on forward and reverse logic Hough transformation
CN115508837B (en) Mobile multi-base sonar target positioning and error analysis method based on UUV cluster
CN115392117B (en) High-frame-rate fuzzy-free acoustic navigation method for underwater high-speed maneuvering platform
Beschnidt et al. “Virtual Waterway”–a traffic simulation environment for inland and coastal waterways
CN113608432B (en) Ship-borne two-degree-of-freedom servo system high dynamic loop parameter adjusting method and system

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 211153 No.32 Changqing street, Jiangning Development Zone, Nanjing City, Jiangsu Province

Patentee after: China Shipbuilding Pengli (Nanjing) Atmospheric and Ocean Information System Co.,Ltd.

Patentee after: CSSC Pengli (Nanjing) Technology Group Co.,Ltd.

Address before: 211153 No.32 Changqing street, Jiangning Development Zone, Nanjing City, Jiangsu Province

Patentee before: CSIC PRIDE (NANJING) ATMOSPHERE MARINE INFORMATION SYSTEM Co.,Ltd.

Patentee before: CSIC NANJING PRIDE TECHNOLOGY GROUP CO.,LTD.