A missile seeker-oriented multi-resolution hybrid simulation method and device

Abstract

The application discloses a kind of multi-resolution hybrid simulation method and device for seeker, the method comprises: obtaining the initial pulse description word data of radar;Obtain the relative distance of radar and seeker, and obtain pulse arrival time;According to relative distance and pulse arrival time, it is judged whether radar satisfies radar interception condition, if it is satisfied, it is obtained that transmitting antenna gain, if it is not satisfied, it is returned to the step of obtaining the initial pulse description word data of radar;Obtain spatial link loss and receiving antenna gain, generate seeker receiving pulse description word data;According to the signal pulse fusion processing of seeker receiving pulse description word data, obtain mixed intermediate frequency sampling signal stream;According to mixed intermediate frequency sampling signal stream, obtain measurement pulse description word data stream;According to signal sorting of measurement pulse description word data stream, obtain sorting result.The application can improve simulation speed and simulation fidelity, and can be widely applied in seeker simulation technical field.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seeker simulation, and in particular to a multi-resolution hybrid simulation method and device for a seeker. BACKGROUND

[0002] At present, research and development of a seeker simulation system are mainly divided into signal level and function level. The signal level simulation is a high-resolution simulation method, which is committed to reflecting the signal flow process and processing process of a real system internal zero intermediate frequency and below, and can reflect the real system situation realistically. However, the signal level simulation has problems of simulation speed and memory space occupation, and is difficult to adapt to the needs of large-scale system simulation or digital twinning, especially when simulating a wideband complex electromagnetic signal environment, the effect is very poor. The function level simulation is a low-resolution simulation method, which aims to realize system functions, does not require to reflect the signal change realistically, and usually only involves some spatial energy calculation and threshold detection calculation based on probability, so the simulation speed is high and the memory occupation is small. However, the function level simulation needs to rely on subjective judgment of people to make decisions in many aspects, and needs to set some parameters artificially, so the subjectivity is strong, the simulation result is easy to be affected, and it does not have strong persuasiveness, so it is not often used when verifying and evaluating system performance. SUMMARY

[0003] Therefore, the embodiments of the present application mainly aim to provide a multi-resolution hybrid simulation method and device for a seeker, so as to solve at least one of the problems in the prior art, and the present application can improve the simulation speed and simulation fidelity.

[0004] To achieve the above-mentioned purpose, one aspect of the embodiments of the present application provides a multi-resolution hybrid simulation method for a seeker, comprising the following steps:

[0005] Obtaining initial pulse description word data of a radar;

[0006] According to the initial pulse description word data, obtaining a relative distance between the radar and a seeker, and obtaining a pulse arrival time;

[0007] According to the relative distance and the pulse arrival time, judging whether the radar meets a radar interception condition, if yes, obtaining a transmitting antenna gain, and if no, returning to the step of obtaining the initial pulse description word data of the radar;

[0008] Obtaining a spatial link loss and a receiving antenna gain, and generating seeker receiving pulse description word data according to the transmitting antenna gain, the spatial link loss and the receiving antenna gain;

[0009] According to the seeker receiving pulse description word data, a plurality of signal pulses are fused to obtain a mixed intermediate frequency sampling signal stream;

[0010] According to the mixed intermediate frequency sampling signal stream, a measurement pulse description word data stream is obtained;

[0011] According to the measurement pulse description word data stream, signal sorting is performed to obtain a sorting result.

[0012] In some embodiments, according to the initial pulse description word data, the relative distance between the radar and the seeker is obtained, and the pulse arrival time is obtained, including the following steps:

[0013] According to the latitude and longitude of the radar of the initial pulse description word data, the relative distance between the radar and the seeker is obtained;

[0014] According to the relative distance and the speed of light, the pulse arrival time is obtained.

[0015] In some embodiments, according to the relative distance and the pulse arrival time, it is judged whether the radar satisfies the radar interception condition, if it satisfies, the transmitting antenna gain is obtained, if it does not satisfy, the step of obtaining the initial pulse description word data of the radar is returned, including the following steps:

[0016] According to the relative distance, a space domain interception condition is constructed;

[0017] According to the pulse arrival time, a time domain interception condition is constructed;

[0018] The frequency scanning range, intermediate frequency bandwidth and switching frequency band speed of the seeker are obtained;

[0019] According to the frequency scanning range, the intermediate frequency bandwidth, the switching frequency band speed and the pulse arrival time, a frequency domain interception condition is constructed;

[0020] It is judged whether the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, if the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, the transmitting antenna gain is obtained, if the radar does not satisfy the space domain interception condition, the time domain interception condition or the frequency domain interception condition, the step of obtaining the initial pulse description word data of the radar is returned;

[0021] The radar interception condition includes the space domain interception condition, the time domain interception condition and the frequency domain interception condition.

[0022] In some embodiments, the spatial link loss and the receiving antenna gain are obtained, including the following steps:

[0023] acquire a wavelength of the electromagnetic wave signal, and acquire signal power attenuation caused by meteorological conditions;

[0024] obtain a path loss according to the relative distance and the wavelength of the electromagnetic wave signal;

[0025] obtain the spatial link loss according to the signal power attenuation and the path loss;

[0026] acquire a relative position of the thunder to the seeker;

[0027] convert the relative position to obtain a radar signal incidence direction;

[0028] preset an antenna pointing direction;

[0029] acquire an included angle between the antenna pointing direction and the radar signal incidence direction;

[0030] obtain the receiving antenna gain through an antenna pattern function according to the included angle.

[0031] In some embodiments, the step of obtaining the mixed intermediate frequency sampling signal stream by fusing a plurality of signal pulses according to the seeker received pulse description word data comprises the following steps:

[0032] signal sampling is performed on adjacent times of a plurality of the signal pulses according to the seeker received pulse description word data to obtain a plurality of pulse sampling signal streams and corresponding pulse sampling time streams;

[0033] the current pulse sampling time stream of the current signal pulse is taken as a total time stream, and the current pulse sampling signal stream of the current signal pulse is taken as a total signal stream;

[0034] a time relationship between the current signal pulse and a next signal pulse is acquired;

[0035] the next signal pulse and the current signal pulse are fused according to the total time stream, the total signal stream and the time relationship to obtain a fused pulse, the fused pulse is taken as the current signal pulse, and the step of taking the current pulse sampling time stream of the current signal pulse as the total time stream and taking the current pulse sampling signal stream of the current signal pulse as the total signal stream is returned until all the signal pulses are fused to obtain the mixed intermediate frequency sampling signal stream.

[0036] In some embodiments, the step of fusing the next signal pulse and the current signal pulse according to the total time stream, the total signal stream and the time relationship comprises the following steps:

[0037] when the time relationship is that adjacent pulses are not aliasing, adding a next said pulse sampling time stream to the total time stream;

[0038] when the time relationship is that adjacent pulses are not aliasing, adding a next said pulse sampling signal stream to the total signal stream;

[0039] when the time relationship is that adjacent pulses are partially aliasing, adding a next said pulse sampling time stream to the total time stream and removing repeated time values;

[0040] when the time relationship is that adjacent pulses are partially aliasing, zero padding the end of the total signal stream to the length of the total time stream, performing a backward traversal operation on the total signal stream, and adding a next said pulse sampling signal stream to the total signal stream until the length of the next said pulse sampling signal stream is exceeded;

[0041] when the time relationship is that adjacent pulses are completely aliasing, performing a backward traversal operation on a next said pulse sampling signal stream, and adding the next said pulse sampling signal stream to the total signal stream at the corresponding time.

[0042] In some embodiments, the method further comprises the following steps:

[0043] superimposing the mixed intermediate frequency sampling signal stream and a thermal noise signal stream to generate a composite pulse signal;

[0044] performing parameter measurement on the composite pulse signal to generate the measurement pulse descriptor word data stream.

[0045] In some embodiments, the method further comprises the following steps:

[0046] performing threshold detection on the composite pulse signal by envelope detection to obtain a detection result and a detection threshold;

[0047] comparing the detection result and the detection threshold, and if the detection result is greater than the detection threshold, taking the time corresponding to the rising edge of the composite pulse signal as a pulse arrival time;

[0048] obtaining a falling edge time of the composite pulse signal;

[0049] obtaining a pulse width according to the falling edge time and the pulse arrival time;

[0050] obtaining a center time of the composite pulse signal according to the pulse width and the pulse arrival time;

[0051] the signal amplitude value at the center time as a pulse amplitude;

[0052] a signal carrier frequency is measured by a frequency domain threshold;

[0053] an azimuth angle of arrival is measured by a four-antenna amplitude comparison direction finding method;

[0054] The measurement pulse description word data stream includes the pulse time of arrival, the pulse width, the pulse amplitude, the signal carrier frequency, and the azimuth angle of arrival.

[0055] To achieve the above object, another aspect of the embodiment of the present application provides a kind of mixed simulation device for seeker, the device includes:

[0056] The first module is used to obtain the initial pulse description word data of radar;

[0057] The second module is used to obtain the relative distance between the radar and seeker according to the initial pulse description word data, and obtain pulse time of arrival;

[0058] The third module is used to judge whether the radar satisfies radar interception condition according to the relative distance and the pulse time of arrival, if satisfies, obtain transmitting antenna gain, if not satisfies, return the step of the initial pulse description word data of radar is obtained;

[0059] The fourth module is used to obtain spatial link loss and receiving antenna gain, and generate seeker receiving pulse description word data according to the transmitting antenna gain, the spatial link loss and the receiving antenna gain;

[0060] The fifth module is used to fuse processing according to the seeker receiving pulse description word data for several signal pulses, and obtain mixed intermediate frequency sampling signal stream;

[0061] The sixth module is used to obtain measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream;

[0062] The seventh module is used to carry out signal sorting according to the measurement pulse description word data stream, and obtain sorting result.

[0063] To achieve the above object, another aspect of the embodiment of the present application provides an electronic device, the electronic device includes memory and processor, the memory stores computer program, the processor executes the computer program and realizes the above-mentioned one kind of mixed simulation method for seeker of multi-resolution.

[0064] To achieve the above object, another aspect of the embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the above-mentioned missile seeker-oriented multi-resolution hybrid simulation method.

[0065] To achieve the above object, another aspect of the embodiment of the present application provides a computer program product or computer program, which comprises computer instructions stored in a computer readable storage medium. A processor of a computer device can read the computer instructions from the computer readable storage medium, and the processor executes the computer instructions to make the computer device execute the above-mentioned missile seeker-oriented multi-resolution hybrid simulation method.

[0066] The embodiment of the present application has at least the following beneficial effects: the present application provides a missile seeker-oriented multi-resolution hybrid simulation method and device, the scheme acquires initial pulse descriptor word data of a radar; according to the initial pulse descriptor word data, the relative distance between the radar and a seeker is acquired, and pulse arrival time is acquired; according to the relative distance and the pulse arrival time, it is judged whether the radar satisfies radar interception conditions, if yes, the transmitting antenna gain is acquired, if not, the step of acquiring the initial pulse descriptor word data of the radar is returned; the spatial link loss and the receiving antenna gain are acquired, according to the transmitting antenna gain, the spatial link loss and the receiving antenna gain, the seeker receiving pulse descriptor word data is generated; according to the seeker receiving pulse descriptor word data, a plurality of signal pulses are fused to obtain a mixed intermediate frequency sampling signal stream; according to the mixed intermediate frequency sampling signal stream, measurement pulse descriptor word data stream is acquired; according to the measurement pulse descriptor word data stream, signal sorting is performed to obtain a sorting result, which can improve the simulation speed and simulation fidelity. BRIEF DESCRIPTION OF DRAWINGS

[0067] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative effort.

[0068] Figure 1 is a flow chart of a missile seeker-oriented multi-resolution hybrid simulation method provided by the embodiment of the present application;

[0069] Figure 2 is a schematic diagram of adjacent pulse conditions provided by the embodiment of the present application;

[0070] Figure 3is a signal arrival time measurement schematic diagram provided by an embodiment of the present application;

[0071] Figure 4 is a composite simulation system simulation scheme framework schematic diagram provided by an embodiment of the present application;

[0072] Figure 5 is a hardware structure schematic diagram of an electronic device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0073] In order to make the purpose, technical scheme and advantages of the present application clearer, the present application is further described in detail below in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not intended to limit the present application. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of the present application, but are only examples of devices and methods consistent with some aspects of the embodiments of the present application as described in the appended claims.

[0074] It should be noted that although the functional modules are divided in the system schematic diagram, and the logical sequence is shown in the flowchart, in some cases, the steps shown or described can be performed in a manner different from the module division in the system or the sequence in the flowchart. The terms "first / S100", "second / S200" in the specification and claims and the above drawings can be used in this paper to describe various concepts, but unless specifically stated, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, the first information can also be referred to as the second information, and similarly, the second information can also be referred to as the first information, without departing from the scope of the embodiments of the present application. Depending on the context, the word "if" as used herein can be interpreted as "when" or "when" or "in response to determining".

[0075] The terms "at least one", "multiple", "each", "any" and the like used in the present application include one, two or more than two, multiple includes two or more than two, each refers to each of the corresponding multiple, and any refers to any one of the multiple.

[0076] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as understood by those skilled in the art to which the present application belongs. The terms used herein are only for the purpose of describing the embodiments of the present application and are not intended to limit the present application.

[0077] Radar seeker, also known as radio seeker, is a kind of detection device. There are three ways to verify and evaluate the working efficiency of the seeker system: formula derivation, field experiment and computer simulation. Pure mathematical formula is used to analyze the system statically and qualitatively, which has low cost but it is difficult to consider all factors and the result is not persuasive enough; field experiment is the most real way to reflect the performance of the seeker system, but it has high cost, long preparation time and is limited by the real environment, so it is difficult to simulate the influence of different battlefield environments on the performance of the system, and the flexibility and reusability are poor, so it is generally used only in the last stage of system development test; in recent years, the performance of computers has improved rapidly, and the computer simulation capability has been greatly enhanced, so computer simulation has gradually become the main means of system performance verification and evaluation. Computer simulation consumes less resources, is flexible and convenient to use, and has strong reusability, and thanks to the continuous development of modern electronic information equipment towards digitization, the fidelity of computer simulation can also be guaranteed. Therefore, the research on the simulation modeling of the seeker system can effectively promote the development of the related field of the seeker system.

[0078] At present, the research and development of the seeker simulation system are mainly divided into signal level and function level. The signal level simulation can truly reflect the flow process and processing steps of the internal signal of the system, and the simulation result is persuasive, but it has the problems of slow simulation speed and large memory occupation, which is difficult to adapt to the needs of large-scale system simulation or digital twinning, especially when simulating the wideband complex electromagnetic signal environment. The function level simulation is a low resolution simulation method, which aims to realize the function of the system and does not require to truly reflect the change of the signal, usually only involves the calculation of some spatial energy and the threshold detection calculation based on probability, so the simulation speed is high and the memory occupation is small. However, the function level simulation needs to rely on subjective judgment of people to make decisions in many aspects, and some parameters need to be set artificially, so the setting of some parameters is disturbed by subjective factors, and the simulation result is easy to be affected, so the persuasion is not enough, and the real processing and transformation of the signal cannot be reflected, so the function level simulation is not often used in system performance verification and evaluation. As can be seen, the signal level simulation and the function level simulation have their own advantages and disadvantages, but neither of them can effectively simulate the seeker system in the wideband electromagnetic environment.

[0079] Therefore, as shown in the Figure 1 embodiments of the present application provide a multi-resolution hybrid simulation method for a seeker, which can include but is not limited to steps S100 to S700:

[0080] In step S100, initial pulse description word data of a radar is obtained;

[0081] In step S200, the relative distance between the radar and the seeker is obtained according to the initial pulse description word data, and the pulse arrival time is obtained.

[0082] Step S300: Based on the relative distance and the pulse arrival time, determine whether the radar meets the radar interception conditions. If it does, obtain the transmitting antenna gain; otherwise, return to the step of obtaining the initial pulse descriptor data of the radar.

[0083] Step S400: Obtain the spatial link loss and the receiving antenna gain; generate seeker receive pulse descriptor data based on the transmitting antenna gain, the spatial link loss, and the receiving antenna gain.

[0084] Step S500: Based on the pulse description word data received by the seeker head, several signal pulses are fused to obtain a mixed intermediate frequency sampling signal stream;

[0085] Step S600: Obtain the measurement pulse description word data stream based on the hybrid intermediate frequency sampling signal stream;

[0086] Step S700: Perform signal sorting based on the measurement pulse description word data stream to obtain the sorting result.

[0087] In step S100 of some embodiments, initial pulse descriptor data of the radar is generated based on radar parameter information. This initial pulse descriptor data may include, but is not limited to, antenna scan type, scan parameters, signal intra-pulse modulation type, modulation parameters, radar signal transmit power, radar pulse width, repetition rate type, repetition rate parameters, radar latitude and longitude, altitude, radar transmit antenna gain, and radar half-power beamwidth. Optionally, the radar parameter information may be user-preset, loaded from a local file, or obtained from a server; no limitation is imposed here.

[0088] In some embodiments, step S200 may include, but is not limited to, steps S210 to S220:

[0089] Step S210: Based on the latitude and longitude of the radar in the initial pulse description word data, obtain the relative distance between the radar and the seeker.

[0090] Step S220: Based on the relative distance and the speed of light, the arrival time of the pulse is obtained.

[0091] In step S210 of some embodiments, the latitude and longitude of the radar that obtains the initial pulse description word data are acquired. The relative distance R between the radar and the seeker can be calculated from the latitude and longitude. Next, the relative azimuth angle between the radar and the seeker can be calculated. Before calculating the relative azimuth angle, a coordinate system transformation is required. Optionally, the coordinate system transformation process is as follows: Geodetic coordinate system (L,B,H) -> Geocentric rectangular coordinate system (X,Y,Z) -> North-East coordinate system (X... N ,YN ,Z N )。 Then the relative bearing angle of the radar to the seeker can be calculated in the North-East coordinate system, and the expression of the relative bearing angle is:

[0092] direction = arctan (X N / Z N )

[0093] wherein, direction represents the relative bearing angle of the radar to the seeker.

[0094] In step S220 of some embodiments, according to the relative distance of the radar to the seeker and the speed of light, the time delay of the radar to the seeker can be calculated. By adding the time delay to the pulse emission time, the pulse arrival time can be obtained. Exemplarily, the calculation formula of the time delay of the radar to the seeker is:

[0095] TimeDelay = R / c

[0096] The calculation formula of the pulse arrival time is:

[0097] TOA = TOE + TimeDelay

[0098] wherein, TimeDelay represents the time delay of the radar to the seeker; R represents the relative distance of the radar to the seeker; c represents the speed of light; TOA represents the pulse arrival time; and TOE represents the pulse emission time.

[0099] In some embodiments, step S300 can include but is not limited to steps S310 to S350:

[0100] Step S310, constructing a space domain interception condition according to the relative distance;

[0101] Step S320, constructing a time domain interception condition according to the pulse arrival time;

[0102] Step S330, acquiring the frequency scanning range of the seeker, the intermediate frequency bandwidth, and the switching frequency range speed;

[0103] Step S340, constructing a frequency domain interception condition according to the frequency scanning range, the intermediate frequency bandwidth, the switching frequency range speed, and the pulse arrival time;

[0104] Step S350, judging whether the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, if the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, acquiring the transmitting antenna gain, if the radar does not satisfy the space domain interception condition, the time domain interception condition or the frequency domain interception condition, returning to the step of acquiring the initial pulse description word data of the radar.

[0105] The radar interception condition includes the space domain interception condition, the time domain interception condition and the frequency domain interception condition.

[0106] In step S310 of some embodiments, the height of the seeker and the height of the radar antenna of the radiation source are acquired, and the maximum visual distance of the seeker to the radiation source radar signal detection is calculated by the height of the seeker and the height of the radar antenna of the radiation source. By comparing the relative distance and the maximum visual distance, it can be judged whether the seeker can receive the pulse signal transmitted by the radar. The maximum visual distance refers to the maximum visual distance of the radar to the target due to the influence of the curvature of the earth's surface. When the relative distance of the radar to the seeker exceeds the maximum visual distance, the pulse signal transmitted by the radar will be blocked by the ground. Exemplarily, the calculation formula of the maximum visual distance is as follows:

[0107]

[0108] In the formula, R smax (km) is the maximum visual distance; h a (m) is the height of the seeker; h t (m) is the height of the radar antenna of the radiation source.

[0109] According to the relative distance and the maximum visual distance, the space domain interception condition can be constructed, which is as follows:

[0110]

[0111] When the relative distance is less than or equal to the maximum visual distance, it means that the space domain screening is passed, that is, the seeker can receive the pulse signal transmitted by the radar; when the relative distance is greater than the maximum visual distance, it means that the space domain screening is not passed, that is, the seeker cannot receive the pulse signal transmitted by the radar.

[0112] In step S320 of some embodiments, the working start time and the working range of the seeker can be pre-set, and the time domain interception condition can be constructed after the pulse arrival time is calculated. According to the time domain interception condition, the time domain screening of the pulse can be performed to judge whether the pulse arrival time is within the working time range of the seeker. Exemplarily, the time domain interception condition for time domain screening is as follows:

[0113] t beginTOA < t begin + duration

[0114] Wherein, duration represents the working time length of the seeker; t begin represents the working start time of the seeker.

[0115] In steps S330 to S340 of some embodiments, the frequency scanning range of the seeker is [fre min ,fre max ], the intermediate frequency bandwidth is bw, and the switching frequency range speed is vec seconds / time. For a pulse with a pulse carrier frequency of rf arriving at time t, the frequency domain interception condition of the seeker receiving the pulse is:

[0116]

[0117] In the formula,

[0118] Wherein, is the minimum receiving frequency of the seeker.

[0119] In step S350 of some embodiments, the radar is subjected to space domain interception, time domain interception, and frequency domain interception, i.e., it is judged whether the radar meets the space domain interception condition, the time domain interception condition, and the frequency domain interception condition. If the radar does not meet one of the space domain interception condition, the time domain interception condition, or the frequency domain interception condition, the step of acquiring the initial pulse descriptor word data of the radar is returned, i.e., the step of generating the initial pulse descriptor word data of the next radar is entered. If the radar meets the space domain interception condition, the time domain interception condition, and the frequency domain interception condition, the next step of acquiring the transmitting antenna gain is entered. On the basis of the initial pulse descriptor word data storing the signal transmitting power, the transmitting antenna gain can be determined according to the antenna scanning type and specific parameters of the radar. The radar antenna rotates in a certain rule during the working process, and the angle between the radar antenna pointing direction and the radar receiver pointing direction determines the transmitting antenna gain of the radar signal received by the seeker, which can be calculated by the antenna pattern function.

[0120] For each type of antenna, the antenna pattern function is different. Exemplarily, taking a Gaussian antenna as an example, the mathematical expression of the main lobe pattern is:

[0121]

[0122] In the formula, θ b is the one-way half-power beamwidth (rad); θ1 is the offset angle corresponding to the pattern.

[0123] The mathematical expression of the pattern other than the main lobe is:

[0124]

[0125] wherein k2 = 1.3916 / sin(0.5θ b ).

[0126] After obtaining the above directional diagram function, the antenna gain can be calculated by the included angle between the radar antenna pointing direction and the radar receiver pointing direction. The antenna gain is the maximum antenna gain multiplied by the antenna directional diagram function, i.e.

[0127] G(θ1) = G max f(θ1)

[0128] wherein G(θ1) is the antenna gain; G max is the maximum antenna gain.

[0129] In some embodiments, the steps of obtaining the spatial link loss and the receiving antenna gain can include, but are not limited to, steps S410 to S480:

[0130] Step S410, obtaining the wavelength of the electromagnetic wave signal and the signal power attenuation caused by the weather condition;

[0131] Step S420, obtaining the path loss according to the relative distance and the wavelength of the electromagnetic wave signal;

[0132] Step S430, obtaining the spatial link loss according to the signal power attenuation and the path loss;

[0133] Step S440, obtaining the relative position of the radar to the seeker head;

[0134] Step S450, transforming the relative position to obtain the radar signal incident direction;

[0135] Step S460, pre-setting the antenna pointing direction;

[0136] Step S470, obtaining the included angle between the antenna pointing direction and the radar signal incident direction;

[0137] Step S480, obtaining the receiving antenna gain through the antenna directional diagram function according to the included angle.

[0138] In steps S410 to S430 of some embodiments, the spatial link loss can be calculated by the pulse carrier frequency, signal propagation distance, weather condition and specific parameters. The spatial link loss is the sum of the path loss caused by signal propagation and the signal power attenuation caused by weather condition. During the propagation of electromagnetic signal in the atmosphere, the signal will be affected by weather and will be attenuated. The loss caused by various weather conditions can be obtained by the corresponding empirical formula. Exemplarily, the calculation formula of path loss is as follows:

[0139]

[0140] In the formula, λ is the wavelength of the electromagnetic signal; and R is the relative distance between the radar and the seeker.

[0141] In steps S440 to S480 of some embodiments, the receiving antenna gain of the corresponding antenna channel is calculated. The radar signal incidence direction can be obtained by converting the relative azimuth of the radar to the seeker, and the conversion formula is as follows:

[0142] θ=(θ0-180°)%360°

[0143] In the formula, θ is the radar signal incidence direction; θ0 is the relative azimuth of the radar to the seeker; and % is the modulo operation.

[0144] In the simulation system, it can be assumed that the directions of the four antennas are 0°, 90°, 180° and 270°, and all are Gaussian antennas. Then the receiving antenna gain can be determined according to the angle between the radar signal incidence direction and the antenna. Assuming that the antenna direction is Then the angle angle between the antenna direction and the radar signal incidence direction is as follows:

[0145]

[0146] After the angle is determined, the receiving antenna gain can be determined by the same method as that for obtaining the transmitting antenna gain.

[0147] In some embodiments, the transmitting antenna gain, the spatial link loss and the receiving antenna gain jointly determine the characteristics of the signal received by the seeker, thereby affecting the generation of the pulse description word data. Optionally, in order to more comprehensively describe this process, other factors such as the modulation mode of the signal, the pulse width, the reflection characteristics of the target, etc. can also be considered. Then the pulse description word data received by the seeker is generated, and the structure of the signal pulse description word group arriving at a certain reconnaissance antenna channel can be as follows:

[0148] PDW=[TOA,PW,Pr,RF,mod type,mod para]

[0149] In the formula, PDW represents a seeker received pulse description word data; TOA represents a pulse arrival time; PW represents a pulse width; Pr represents a pulse received power; RF represents a pulse corresponding carrier frequency; mod type represents a pulse intra-pulse modulation type; and mod para represents a pulse intra-pulse modulation parameter.

[0150] In step S500 of some embodiments, a function level simulation-signal level simulation conversion is performed, a simulation acceleration method for a wideband electromagnetic environment is proposed, and in view of the problem that a computer memory is difficult to bear a large bandwidth and full-time intermediate frequency sampling signal stream, it is proposed to divide the full-time intermediate frequency sampling signal stream based on a time window, and only simulate the intermediate frequency sampling signal stream near the pulse signal.

[0151] In some embodiments, step S500 can include but is not limited to steps S510 to S540:

[0152] Step S510, according to the seeker received pulse description word data, signal sampling is performed on the adjacent time of a plurality of signal pulses, and a plurality of pulse sampling signal streams and corresponding pulse sampling time streams are obtained;

[0153] Step S520, taking the current pulse sampling time stream of the current signal pulse as a total time stream, and taking the current pulse sampling signal stream of the current signal pulse as a total signal stream;

[0154] Step S530, obtaining the time relationship between the current signal pulse and the next signal pulse;

[0155] Step S540, according to the total time stream, the total signal stream and the time relationship, performing fusion processing on the next signal pulse and the current signal pulse to obtain a fusion pulse, taking the fusion pulse as the current signal pulse, returning to the step of taking the current pulse sampling time stream of the current signal pulse as a total time stream and taking the current pulse sampling signal stream of the current signal pulse as a total signal stream, until all the signal pulses are fused to obtain the mixed intermediate frequency sampling signal stream.

[0156] In step S510 of some embodiments, according to each seeker received pulse description word data, an initial intermediate frequency sampling signal stream near the corresponding signal pulse is generated. In order to ensure the authenticity of subsequent detection and measurement, the time of 1us left and right of the signal pulse is optionally sampled. Then the initial intermediate frequency sampling signal stream can be obtained as follows:

[0157] TimeStream = [[SampleTime-1], [SampleTime-2], [SampleTime-3]...]

[0158] SignalStream = [[SampleSignal-1], [SampleSignal-2], [SampleSignal-3]…]

[0159] Wherein, TimeStream is a collection of pulse sampling time streams; SignalStream is a collection of pulse sampling signal streams; [SampleTime-1], [SampleTime-2], [SampleTime-3]… are each pulse sampling time stream; [SampleSignal-1], [SampleSignal-2], [SampleSignal-3]… are each pulse sampling signal stream.

[0160] In step S520 of some embodiments, the basic idea of the intermediate frequency sampling signal stream fusion is to take the sampling time stream and the sampling signal stream of the first signal pulse as the total time stream and the total signal stream, and then fuse the subsequent signal pulses one by one. Optionally, the current pulse sampling time stream of the current signal pulse is taken as the total time stream, and the current pulse sampling signal stream of the current signal pulse is taken as the total signal stream, which lays a foundation for the subsequent fusion of signal pulses one by one.

[0161] In step S530 of some embodiments, since the pulse description word information has been sorted according to the pulse arrival time in the generation of the seeker receiving pulse description word data by default, it can be ensured that the starting time of the current signal pulse is before the starting time of the subsequent signal pulse, and then the time relationship between the current signal pulse and the next signal pulse has the following three kinds: the current signal pulse and the next signal pulse are not aliasing, the current signal pulse and the next signal pulse are partially aliasing, and the next signal pulse is completely aliasing in the current signal pulse. Then, as shown in the adjacent pulse case, there are three kinds of adjacent pulse cases: adjacent pulse aliasing, adjacent pulse partial aliasing, and adjacent pulse complete aliasing. Figure 2

[0162] In some embodiments, the fusion processing of the next signal pulse and the current signal pulse according to the total time stream, the total signal stream, and the time relationship can include but is not limited to steps S541 to S545:

[0163] Step S541, when the time relationship is adjacent pulse aliasing, adding the next pulse sampling time stream to the total time stream;

[0164] Step S542, when the time relationship is adjacent pulse aliasing, adding the next pulse sampling signal stream to the total signal stream;

[0165] ​Step S543, when the time relationship is adjacent pulse partial aliasing, add the next pulse sampling time stream to the total time stream, and remove the repeated time values;

[0166] Step S544, when the time relationship is adjacent pulse partial aliasing, pad 0 to the end of the total signal stream to make the total signal stream equal to the total time stream, perform a backward traversal operation on the total signal stream, and add the next pulse sampling signal stream to the total signal stream until the length of the next pulse sampling signal stream is exceeded.

[0167] Step S545, when the time relationship is adjacent pulse complete aliasing, perform a backward traversal operation on the next pulse sampling signal stream, and add the next pulse sampling signal stream to the total signal stream at the corresponding time.

[0168] In steps S541 to S542 of some embodiments, if the start time of the next signal pulse is greater than the end time of the current signal pulse, the adjacent pulses are not aliased, and the sampling time stream of the next signal pulse is directly added to the total time stream, and the sampling signal stream of the next signal pulse is added to the total signal stream.

[0169] In steps S543 to S544 of some embodiments, if the start time of the next signal pulse is less than the end time of the current signal pulse, but the end time exceeds the end time of the current signal pulse, the adjacent pulses are partially aliased, the total time stream is added to the next pulse sampling time stream, the repeated time values are removed, 0 is added to the end of the total signal stream to make the total signal stream equal to the total time stream, and finally, the next signal pulse sampling signal stream is added to the total signal stream from the end of the total signal stream to the beginning in a backward traversal operation until the length of the next pulse sampling signal stream is exceeded.

[0170] In step S545 of some embodiments, if the end time of the next signal pulse is less than the end time of the current signal pulse, the adjacent pulses are completely aliased, the total time stream remains unchanged, the next signal pulse sampling signal stream is traversed backward, and the signal values of the next signal pulse sampling signal stream are added to the signal values of the total signal stream at the corresponding time. The corresponding time here refers to the corresponding part of the next signal pulse sampling time stream and the total time stream, so that the signal values of the next signal pulse can be superimposed on the signal values of the total signal stream corresponding to the total time stream.

[0171] By continuously fusing the sampling time stream and the sampling signal stream of adjacent pulses, a total time stream and signal stream can be finally formed as a mixed intermediate frequency sampling signal stream of the radiation source, which is used in the subsequent link of analog seeker detection and measurement of signal parameters. It can be seen that this simulation acceleration method fuses the isolated pulse sampling signal stream according to the sampling time stream, which can effectively reproduce the pulse aliasing situation. In the simulation system design, this method is used for simulation, which improves the simulation speed under the premise of ensuring the fidelity of the simulation system.

[0172] In some embodiments, step S600 can include but is not limited to steps S610 to S620:

[0173] Step S610, superimposing the mixed intermediate frequency sampling signal stream and the thermal noise signal stream to generate a composite pulse signal;

[0174] Step S620, performing parameter measurement on the composite pulse signal to generate the measurement pulse descriptor word data stream.

[0175] In step S610 of some embodiments, the reconnaissance receiver module superimposes the received mixed intermediate frequency sampling signal stream of the radiation source and the thermal noise signal stream, and generates a signal stream after superposition as a composite pulse signal for subsequent parameter measurement.

[0176] In step S620 of some embodiments, the signal detection and measurement model is used to measure the time domain, frequency domain and spatial domain parameters of the mixed intermediate frequency sampling signal stream of the radar radiation source in the composite pulse signal, and through digital signal processing of the mixed intermediate frequency sampling signal of the radar radiation source, the signal time of arrival (TOA), pulse width (PW), pulse amplitude (PA), signal carrier frequency (RF) and direction of arrival (DOA) are measured. Angle) and other parameters form a measurement pulse descriptor word data stream, which lays a foundation for subsequent sorting and identification processing.

[0177] In some embodiments, step S620 can include but is not limited to steps S621 to S628:

[0178] Step S621, threshold detection is performed on the composite pulse signal by envelope detection method to obtain detection results and detection threshold;

[0179] Step S622, comparing the detection results with the detection threshold, if the detection results are greater than the detection threshold, the rising edge corresponding time of the composite pulse signal is taken as the pulse arrival time;

[0180] Step S623, obtaining the falling edge time of the composite pulse signal;

[0181] Step S624, obtaining pulse width according to the falling edge time and the pulse arrival time;

[0182] Step S625, obtaining the center time of the composite pulse signal according to the pulse width and the pulse arrival time;

[0183] Step S626, taking the signal amplitude value at the center time as pulse amplitude;

[0184] Step S627, measuring signal carrier frequency through frequency domain threshold;

[0185] Step S628, measuring arrival azimuth through four-antenna amplitude comparison direction finding method;

[0186] The measurement pulse description word data stream includes the pulse arrival time, the pulse width, the pulse amplitude, the signal carrier frequency and the arrival azimuth.

[0187] In steps S621 to S622 of some embodiments, the signal arrival time is measured, referring to Figure 3 , Figure 3 The figure shows the measurement of signal arrival time for a seeker. By the method of envelope detection, the input composite pulse signal is threshold detected to obtain the detection result and the detection threshold. The detection result and the detection threshold are compared. If the detection result is greater than the detection threshold, it means that the composite pulse signal is detected, and the rising edge time of the composite pulse signal is taken as the pulse arrival time TOA. Each time the rising edge of the composite pulse signal arrives, the counting result of the time counter is accumulated with the last pulse arrival time to form a new pulse arrival time and output as a parameter. Then the time counter is emptied and the counting is restarted until the arrival of the next pulse rising edge.

[0188] In steps S623 to S624 of some embodiments, the pulse width is measured. The estimation of the pulse width is similar to the estimation of the pulse arrival time, both of which use the counter to time, but the falling edge PFE is added to the timing, and the calculation formula is:

[0189]

[0190] In the formula, PFE is the falling edge time; Δt is the minimum value of quantization; INT is the rounding operation. Then the pulse width PW is the quantization value of the time between the falling edge time PFE and the pulse arrival time TOA, with Δt as the unit.

[0191] In steps S625 to S626 of some embodiments, the pulse amplitude is measured, and the measurement of the pulse amplitude PA utilizes the measured pulse arrival time and pulse width. For example, the center time of the composite pulse signal is calculated according to the measured pulse arrival time and pulse width, and the signal amplitude value at the center time is taken as the measurement value of the pulse amplitude. Optionally, the calculation formula of the center time is:

[0192]

[0193] where T center is the center time. The center time refers to the middle time value of a pulse signal.

[0194] In steps S627 to S628 of some embodiments, the frequency domain parameter signal carrier frequency RF is measured using the form of a frequency domain threshold, and the radar signal arrival azimuth DOA is measured using the four-antenna amplitude comparison direction finding method.

[0195] In step S700 of some embodiments, according to the generated measurement pulse description word data stream, the radar signal can be sorted to obtain a sorting result. For example, in the radar signal sorting process, the OPTICS clustering algorithm is used for pre-sorting and the PRI transformation method is used for main sorting. In the pre-sorting stage, the radar signal arrival azimuth DOA and the frequency domain parameter signal carrier frequency RF can be used as the basis for clustering to preliminarily sort the radar signal. The purpose of pre-sorting is to dilute the pulse stream so that the subsequent processing (including but not limited to matching deduction of known radar radiation sources and preliminary signal classification) is simpler. In the main sorting stage, the PRI transformation method further sorts the pulses by processing the pulse arrival time TOA. By setting a detection threshold, the PRI transformation value exceeding the threshold can be identified to estimate the PRI value. This method can effectively suppress the influence of sub-harmonics and improve the accuracy of PRI estimation. Through the application of the OPTICS clustering algorithm and the PRI transformation method in radar signal sorting, the radar signal can be effectively classified and identified.

[0196] Among them, OPTICS (Ordering Points To Identify the Clustering Structure) is a density-based clustering algorithm, which does not need to specify the number of clusters in advance, can identify clusters of any shape, and can identify noise points. The algorithm processes sample points in the data set by constructing two queues, an ordered queue O and a result queue R. The ordered queue O stores the core objects and their density direct objects, and is arranged in ascending order of reachable distance; the result queue R stores the output order of the sample points. The algorithm flow includes steps such as selecting core objects, calculating reachable distance, updating queues, etc., until all points are processed or there is no core point.

[0197] PRI transform method (Pulse Repetition Interval Transform) is a secondary processing method based on pulse arrival time, which is used to estimate the pulse repetition period (PRI) value. By weighting the autocorrelation function of the signal with a phase factor, a PRI transform spectrum can be obtained, and a peak value will appear at the PRI value of the signal, thereby achieving the estimation of the PRI. Discrete PRI transform estimates the PRI value by dividing the PRI value range into multiple parts (called PRI bins) and calculating the PRI transform value in each bin.

[0198] As shown in Figure 4 , taking the composite simulation system loaded with the seeker-oriented multi-resolution hybrid simulation method of the embodiment of the application as an example, hybrid simulation is performed, and the following steps are included:

[0199] Step 1: Generate the pulse description word data of the radar end according to the radar parameter information set by the user or loaded from a file, including antenna scanning type, scanning parameter, signal intra-pulse modulation type, modulation parameter, radar signal transmission power, radar pulse width, frequency type, frequency parameter, latitude and longitude of the radar, height, radar transmission antenna gain, radar half-power beam width, etc.

[0200] Step 2: Calculate the relative distance between the radar and the seeker from the latitude and longitude, and then calculate the relative azimuth angle between the radar and the seeker. Before calculating the relative azimuth angle, the coordinate system needs to be converted to the North-East coordinate system. Then the relative azimuth angle of the radar to the seeker can be calculated in the North-East coordinate system.

[0201] Step 3: First, calculate the time delay of the radar to the seeker, and add the time delay to the time of pulse emission (TOE) to obtain the time of arrival (TOA).

[0202] Step 4: Perform spatial interception, time domain interception, and frequency domain interception on the radar, and the specific steps are as follows:

[0203] Step 1: Spatial interception is performed: after the relative distance of the radiation source radar to the seeker is calculated, spatial screening can be performed by comparing the relative distance with the maximum range of the seeker for the radar signal detection to determine whether the seeker can receive the pulse signal transmitted by the radar. If not, the subsequent process of the radar pulse description word generation is directly skipped, and the next radar pulse description word generation step is entered.

[0204] Step 2: Time domain interception is performed: in the simulation system, the user sets the working start time and working range of the seeker, and can perform time domain screening on the pulse after calculating the pulse arrival time to determine whether the pulse arrival time is within the working time range of the seeker. If the pulse arrival time is not within the working time range of the seeker, the subsequent process of the radar pulse description word generation is directly skipped, and the next radar pulse description word generation step is entered.

[0205] Step 3: Frequency domain interception is performed: the frequency scanning range, intermediate frequency bandwidth, and frequency band switching speed of the seeker receiver are obtained, and the frequency domain interception condition is constructed according to the frequency scanning range, intermediate frequency bandwidth, and frequency band switching speed and the pulse arrival time. If the frequency domain interception condition is not met, the subsequent process of the radar pulse description word generation is directly skipped, and the next radar pulse description word generation step is entered.

[0206] Step 5: Based on the signal transmission power stored in the pulse description word, the antenna transmission gain is determined according to the antenna scanning type and specific parameters of the radar, and the antenna transmission gain can be calculated by the antenna pattern function.

[0207] Step 6: The path loss caused by signal propagation and the signal power attenuation caused by weather conditions are added to obtain the spatial link loss.

[0208] Step 7: The calculation of the corresponding antenna channel receiving antenna gain is performed. In the simulation system, it is assumed that the directions of the four antennas are 0°, 90°, 180°, and 270°, and they are all Gaussian type antennas. The receiving antenna gain is determined according to the angle between the radar signal incident direction and the antenna. After the angle is determined, the same method as the determination of the transmitting antenna gain is used to determine the receiving antenna gain.

[0209] Step 8: The seeker receives the pulse description word data to obtain the signal pulse description word composition structure arriving at a certain reconnaissance antenna channel.

[0210] Step nine: functional level simulation-signal level simulation conversion. According to the pulse description word information, the corresponding pulse vicinity intermediate frequency sampling signal stream is generated. After generating a plurality of signal pulse vicinity sampling signal streams according to the pulse description word information, the sampling signal streams are fused, and through continuously fusing the sampling time streams and the sampling signal streams of adjacent pulses, a total time stream and a signal stream can be finally formed as a mixed intermediate frequency sampling signal stream of the radiation source.

[0211] Step ten: the reconnaissance receiver module superimposes the received mixed intermediate frequency sampling signal stream of the radiation source and the thermal noise signal stream.

[0212] Step eleven: time domain, frequency domain and space domain parameter measurement. Through the signal detection and measurement model, the time domain, frequency domain and space domain parameters of the radar mixed intermediate frequency sampling signal stream are measured, and through the digital signal processing of the radar mixed intermediate frequency signal, the signal arrival time, pulse width, pulse amplitude, signal carrier frequency and arrival direction angle and other parameters are measured, and the measurement pulse description word data stream is formed for subsequent sorting and identification processing. The specific steps are as follows:

[0213] Step 1: signal arrival time measurement. The input signal is envelope detected, and the detection result and the detection threshold are compared. If it is greater than the threshold, it is considered that the pulse is detected, and the pulse rising edge corresponding time is taken as the pulse arrival time. When each pulse rising edge comes, the counting result of the time counter is added to the previous pulse arrival time to form a new pulse arrival time and is taken as a parameter output. Then the time counter is emptied and the counting is restarted until the next pulse rising edge comes.

[0214] Step 2: pulse width measurement. The estimation of pulse width is similar to the estimation of pulse arrival time, both of which use the counter to time, but the falling edge is added to the timing, and the pulse width is in the unit of the quantized minimum value, and the quantized value of the time interval between the falling edge time and the rising edge time of the pulse.

[0215] Step 3: pulse amplitude measurement. According to the measured pulse arrival time and pulse width, the pulse amplitude is measured, the center time of the pulse is calculated, and the signal amplitude at the center time is taken as the measurement value of the pulse amplitude.

[0216] Step 4: frequency domain parameter signal carrier frequency measurement using frequency domain threshold.

[0217] Step 5: four-antenna amplitude comparison direction finding method is used to measure the radar signal arrival azimuth.

[0218] Step twelve: performing sorting identification. The sorting identification is divided into two steps, namely pre-sorting and main sorting. The formed measurement pulse description word data stream is pre-sorted by an OPTICS clustering algorithm, the pre-sorted measurement pulse description word data stream is main-sorted by using a PRI transformation method, and a sorting result is output.

[0219] The embodiment of the present application also provides a seeker-oriented multi-resolution hybrid simulation device, which can realize the above-mentioned seeker-oriented multi-resolution hybrid simulation method.

[0220] The first module is used for acquiring initial pulse description word data of a radar.

[0221] The second module is used for acquiring a relative distance between the radar and a seeker according to the initial pulse description word data, and acquiring a pulse arrival time.

[0222] The third module is used for judging whether the radar satisfies a radar interception condition according to the relative distance and the pulse arrival time, acquiring a transmitting antenna gain if the radar satisfies the radar interception condition, and returning to the step of acquiring the initial pulse description word data of the radar if the radar does not satisfy the radar interception condition.

[0223] The fourth module is used for acquiring a spatial link loss and a receiving antenna gain, and generating seeker reception pulse description word data according to the transmitting antenna gain, the spatial link loss and the receiving antenna gain.

[0224] The fifth module is used for performing fusion processing on a plurality of signal pulses according to the seeker reception pulse description word data, and obtaining a hybrid intermediate frequency sampling signal stream.

[0225] The sixth module is used for acquiring a measurement pulse description word data stream according to the hybrid intermediate frequency sampling signal stream.

[0226] The seventh module is used for performing signal sorting according to the measurement pulse description word data stream, and obtaining a sorting result.

[0227] It can be understood that the contents in the above method embodiments are all applicable to the present device embodiment, the present device embodiment specifically realizes the same functions as the above method embodiments, and achieves the same beneficial effects as the above method embodiments.

[0228] The embodiment of the present application also provides an electronic device, which comprises a processor and a memory, the memory stores a computer program, and the processor realizes the above-mentioned seeker-oriented multi-resolution hybrid simulation method when executing the computer program. The electronic device can be any intelligent terminal, such as a tablet computer or a vehicle-mounted computer.

[0229] It can be understood that the contents in the above method embodiments are all applicable to the present device embodiments, the present device embodiments specifically implement the functions same as those of the above method embodiments, and achieve the same beneficial effects as those of the above method embodiments.

[0230] Reference Figure 5 , Figure 5 The hardware structure of the electronic device of another embodiment is illustrated, and the electronic device comprises:

[0231] The processor 801 can be implemented in the form of a general-purpose CPU (Central Processing Unit), a microprocessor, an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits, and is used to execute related programs to implement the technical solutions provided by the present embodiment.

[0232] The memory 802 can be implemented in the form of a ROM (Read Only Memory), a static storage device, a dynamic storage device, or a RAM (Random Access Memory). The memory 802 can store an operating system and other application programs. When the technical solutions provided by the present embodiment are implemented by software or firmware, the related program codes are stored in the memory 802 and are called and executed by the processor 801 to implement a kind of multi-resolution hybrid simulation method for guide head.

[0233] The input / output interface 803 is used to realize information input and output.

[0234] The communication interface 804 is used to realize the communication interaction between the present device and other devices. The communication can be realized by a wired manner (for example, USB, network cable, etc.), or can be realized by a wireless manner (for example, mobile network, WIFI, Bluetooth, etc.).

[0235] The bus 805 is used to transmit information between various components (for example, the processor 801, the memory 802, the input / output interface 803, and the communication interface 804) of the device.

[0236] The processor 801, the memory 802, the input / output interface 803, and the communication interface 804 are connected to each other through the bus 805 to realize the communication connection between them in the device.

[0237] The present embodiment further provides a computer readable storage medium, which stores a computer program. The computer program is executed by a processor to implement the above-mentioned multi-resolution hybrid simulation method for guide head.

[0238] It can be understood that the contents in the above method embodiments are applicable to the present storage medium embodiments, the present storage medium embodiments specifically implement the functions same as the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0239] The embodiment of the present application also provides a computer program product or a computer program, which comprises computer instructions stored in a computer readable storage medium. A processor of a computer device can read the computer instructions from the computer readable storage medium, and the processor executes the computer instructions, so that the computer device executes the aforementioned multi-resolution hybrid simulation method for a seeker.

[0240] The multi-resolution hybrid simulation method for a seeker and the device of the embodiment of the present application have the following advantages:

[0241] 1. The embodiment of the present application adopts different simulation granularities according to the characteristics and importance of different steps in simulation, proposes a hybrid seeker simulation method combining function level and signal level, builds a hybrid seeker simulation system combining function level and signal level, and realizes the balance between simulation speed and simulation fidelity.

[0242] 2. The function level modeling of a traditional complex electromagnetic signal environment directly generates pulse description word data at a receiving end of a reconnaissance receiver, and cannot reflect the influence of signal propagation in space. The embodiment of the present application improves the generation method of pulse description word data of the function level modeling of a complex electromagnetic signal environment. The construction of pulse description word data starts from a radar transmitting end, and considers the influences of signal types, intra-pulse modulation modes, carrier frequency information, transmitting antenna gains, spatial link losses and other influences in the signal transmitting and propagation process, so that the pulse description word data reaching the seeker receiver is more realistic.

[0243] 3. The embodiment of the present application divides the full-time intermediate frequency sampling signal stream based on a time window, only simulates the intermediate frequency sampling signal stream near the pulse signal, then fuses the mixed pulse, effectively reduces the simulation memory occupation without affecting the simulation result, and also speeds up the speed of the subsequent signal processing link.

[0244] To sum up, in the complex signal environment modeling process, the embodiment of the application improves the generation method of the pulse description word data by using the function level simulation, constructs the pulse description word from the radar transmitting end, and reflects the changes of the signal in the transmitting, propagating and receiving processes in the pulse description word through calculation, so that the pulse description word data reaching the receiving end is more realistic; in the seeker modeling process, the signal level simulation is used, the core modules of the seeker, i.e., the signal detection and parameter measurement module and the signal sorting module, are researched, and the algorithm suitable for the broadband electromagnetic signal environment is selected for simulation modeling; for the interface problem of different granularity simulation models, the signal flow generation method based on the pulse description word-sampling signal mapping relationship is used, and the effective information transmission problem between different granularity modules is solved.

[0245] In some alternative embodiments, the functions / operations mentioned in the block diagrams can not occur in the order mentioned in the operation diagram. For example, depending on the functions / operations involved, two blocks shown in succession can actually be executed substantially concurrently or the blocks can sometimes be executed in reverse order. Furthermore, the embodiments presented and described in the flow diagrams of the application are only provided by way of example, with the purpose of providing a more comprehensive understanding of the technology. The disclosed methods are not limited to the operations and logical flows presented herein. Alternative embodiments are contemplated in which the order of various operations is changed and in which sub-operations described as part of a larger operation are independently executed.

[0246] Furthermore, although the present application is described in the context of functional modules, it is to be understood that one or more of the functions and / or features described can be integrated in a single physical device and / or software module, or one or more functions and / or features can be implemented in separate physical devices or software modules. It is also to be understood that detailed discussion of the actual implementation of each module is unnecessary to an understanding of the present application. Rather, the actual implementation is within the routine of an engineer's knowledge given the property, functionality and internal relationships of the various functional modules disclosed herein. Accordingly, the present application is not limited to the specific details of the functional modules described herein. In addition, it is to be understood that the specific concepts disclosed are merely illustrative and are not intended to limit the scope of the present application, which is defined by the appended claims and their equivalents.

[0247] If the functions are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application or the parts of the technical solutions that essentially contribute to the prior art or the parts of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0248] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered a list of executable instructions for implementing logic functions, and can be specifically embodied in any computer-readable medium for use by an instruction execution system, apparatus, or device, such as a computer-based system, a system including a processor, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions, or in conjunction with these instructions execution systems, apparatuses, or devices. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport programs for use by an instruction execution system, apparatus, or device, or in conjunction with these instruction execution systems, apparatuses, or devices.

[0249] More specific examples (a non-exhaustive list) of the computer-readable medium include the following: an electrical connection having one or more wires (electrical devices), a portable computer diskette (magnetic devices), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium can even be paper or other suitable medium on which the program can be printed, as the program can be electronically obtained, for example, by optical scanning of the paper or other medium, followed by editing, interpreting, or otherwise processing, if necessary, in other suitable ways to be electronically obtained, and then stored in the computer memory.

[0250] It should be understood that aspects of the application can be implemented in hardware, software, firmware or a combination thereof. In the above embodiments, various steps or methods can be implemented in software or firmware which is stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, and in another embodiment, any of the following technologies, known in the art, can be used: a combination of discrete logic circuits having logic gates for implementing logic functions upon an application of data signals, application specific integrated circuits having logic gates, field programmable gate arrays (FPGA), or other components, in combination or as the case can be.

[0251] In the description of the present specification, the description of the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" and the like means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Also, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0252] Although embodiments of the present application have been shown and described, it would be appreciated by those skilled in the art that changes, modifications, alternatives and variations to these embodiments could be made without departing from the principles and spirit of the application, the scope of which is defined by the claims and their equivalents.

[0253] The above is a specific description of the preferred embodiments of the present application, but the present application is not limited to the described embodiments, and those skilled in the art can make various equivalent modifications or replacements without departing from the spirit of the present application, and these equivalent modifications or replacements are included in the scope defined by the claims of the present application.

Claims

1. A missile seeker oriented multi-resolution hybrid simulation method, characterized in that, The method comprises the following steps: acquiring initial pulse descriptor word data of a radar; acquiring relative distance between the radar and a seeker according to the initial pulse descriptor word data, and acquiring pulse arrival time; judging whether the radar satisfies radar interception conditions according to the relative distance and the pulse arrival time, acquiring transmitting antenna gain if the radar satisfies the radar interception conditions, and returning to the step of acquiring the initial pulse descriptor word data of the radar if the radar does not satisfy the radar interception conditions; acquiring spatial link loss and receiving antenna gain, and generating seeker receiving pulse descriptor word data according to the transmitting antenna gain, the spatial link loss and the receiving antenna gain; performing fusion processing on a plurality of signal pulses according to the seeker receiving pulse descriptor word data to obtain a mixed intermediate frequency sampling signal stream; acquiring a measured pulse descriptor word data stream according to the mixed intermediate frequency sampling signal stream; performing signal sorting according to the measured pulse descriptor word data stream to obtain a sorting result.

2. The method according to claim 1, wherein, The step of acquiring the relative distance between the radar and the seeker according to the initial pulse descriptor word data comprises the following steps: obtaining the relative distance between the radar and the seeker according to the latitude and longitude of the radar in the initial pulse descriptor word data; obtaining the pulse arrival time according to the relative distance and the speed of light.

3. The method of claim 1, wherein, The step of judging whether the radar satisfies radar interception conditions according to the relative distance and the pulse arrival time, acquiring transmitting antenna gain if the radar satisfies the radar interception conditions, and returning to the step of acquiring the initial pulse descriptor word data of the radar if the radar does not satisfy the radar interception conditions comprises the following steps: constructing a space domain interception condition according to the relative distance; constructing a time domain interception condition according to the pulse arrival time; acquiring a frequency scanning range of the seeker, an intermediate frequency bandwidth and a switching frequency range speed of the seeker; constructing a frequency domain interception condition according to the frequency scanning range, the intermediate frequency bandwidth, the switching frequency range speed and the pulse arrival time; judging whether the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, acquiring transmitting antenna gain if the radar satisfies the space domain interception condition, the time domain interception condition and the frequency domain interception condition, and returning to the step of acquiring the initial pulse descriptor word data of the radar if the radar does not satisfy the space domain interception condition, the time domain interception condition or the frequency domain interception condition; wherein the radar interception conditions comprise the space domain interception condition, the time domain interception condition and the frequency domain interception condition.

4. The method of claim 1, wherein, The step of acquiring the spatial link loss and the receiving antenna gain comprises the following steps: acquiring an electromagnetic wave signal wavelength and signal power attenuation caused by meteorological conditions; obtaining a path loss according to the relative distance and the electromagnetic wave signal wavelength; obtaining the spatial link loss according to the signal power attenuation and the path loss; acquiring a relative direction between the radar and the seeker; obtaining a radar signal incident direction by converting the relative direction; pre-setting an antenna pointing direction; acquiring an angle between the antenna pointing direction and the radar signal incident direction; obtaining the receiving antenna gain by an antenna directional diagram function according to the angle.

5. The missile seeker oriented multi-resolution hybrid simulation method according to claim 1, wherein, The method comprises the following steps of: According to the pulse description word data received by the seeker head, a plurality of signal pulses are fused to obtain a mixed intermediate frequency sampling signal stream. According to the pulse description word data received by the seeker head, a plurality of signal pulses are fused to obtain a mixed intermediate frequency sampling signal stream. The current pulse sampling time stream of the current signal pulse is taken as a total time stream, and the current pulse sampling signal stream of the current signal pulse is taken as a total signal stream. The time relationship between the current signal pulse and the next signal pulse is obtained.

6. The missile seeker-oriented multi-resolution hybrid simulation method according to claim 5, wherein, According to the total time stream, the total signal stream and the time relationship, the next signal pulse and the current signal pulse are fused to obtain a fused pulse, and the fused pulse is taken as the current signal pulse. The method comprises the following steps of: When the time relationship is adjacent pulse non-aliasing, the next pulse sampling time stream is added to the total time stream. When the time relationship is adjacent pulse non-aliasing, the next pulse sampling signal stream is added to the total signal stream. When the time relationship is adjacent pulse partial aliasing, the next pulse sampling time stream is added to the total time stream, and the repeated time values are removed. When the time relationship is adjacent pulse partial aliasing, 0 is added to the end of the total signal stream to make the total signal stream equal to the total time stream in length, the total signal stream is traversed from back to front, and the next pulse sampling signal stream is added to the total signal stream until the length of the next pulse sampling signal stream is exceeded.

7. The method of claim 1, wherein, When the time relationship is adjacent pulse complete aliasing, the next pulse sampling signal stream is traversed from back to front, and the next pulse sampling signal stream is added to the total signal stream at the corresponding time. The method comprises the following steps of: The mixed intermediate frequency sampling signal stream is superimposed with a thermal noise signal stream to generate a composite pulse signal.

8. The method according to claim 7, wherein, The composite pulse signal is measured to generate the measurement pulse description word data stream. The method comprises the following steps of: The composite pulse signal is threshold detected by envelope detection to obtain a detection result and a detection threshold. The detection result is compared with the detection threshold, and if the detection result is greater than the detection threshold, the rising edge time of the composite pulse signal is taken as a pulse arrival time. The falling edge time of the composite pulse signal is obtained. According to the falling edge time and the pulse arrival time, a pulse width is obtained; According to the pulse width and the pulse arrival time, a center time of the composite pulse signal is acquired; A signal amplitude value at the center time is taken as a pulse amplitude; A signal carrier frequency is measured through a frequency domain threshold; An arrival azimuth is measured through a four-antenna amplitude comparison direction finding method. The measurement pulse description word data stream includes the pulse arrival time, the pulse width, the pulse amplitude, the signal carrier frequency, and the arrival azimuth.

9. A missile seeker oriented multi-resolution hybrid simulation apparatus, characterized by, The method comprises: A first module is configured to acquire initial pulse description word data of a radar; A second module is configured to acquire a relative distance between the radar and a seeker according to the initial pulse description word data, and acquire a pulse arrival time; A third module is configured to judge whether the radar satisfies a radar intercepting condition according to the relative distance and the pulse arrival time, acquire a transmitting antenna gain if the radar satisfies the radar intercepting condition, and return to the step of acquiring the initial pulse description word data of the radar if the radar does not satisfy the radar intercepting condition; A fourth module is configured to acquire a spatial link loss and a receiving antenna gain, and generate seeker receiving pulse description word data according to the transmitting antenna gain, the spatial link loss, and the receiving antenna gain; A fifth module is configured to fuse a plurality of signal pulses according to the seeker receiving pulse description word data, and obtain a mixed intermediate frequency sampling signal stream; A sixth module is configured to acquire a measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream; A seventh module is configured to perform signal sorting according to the measurement pulse description word data stream, and obtain a sorting result.

10. An electronic device, comprising: The method comprises a processor and a memory; The memory is configured to store a program; The processor is configured to execute the program to implement the method in any one of claims 1 to 8.

Images