Multi-resolution hybrid simulation method and device for seeker

Abstract

The invention discloses a seeker-oriented multi-resolution hybrid simulation method and device. The method comprises the following steps: acquiring initial pulse description word data of a radar; the relative distance between the radar and the seeker is obtained, and the pulse arrival time is obtained; according to the relative distance and the pulse arrival time, judging whether the radar meets a radar interception condition, if so, obtaining the gain of the transmitting antenna, and if not, returning to the step of obtaining the initial pulse description word data of the radar; acquiring space link loss and receiving antenna gain, and generating seeker receiving pulse description word data; performing fusion processing on the signal pulses according to pulse description word data received by the seeker to obtain a mixed intermediate frequency sampling signal flow; acquiring a measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream; and carrying out signal sorting according to the measurement pulse description word data stream to obtain a sorting result. The method can improve the simulation speed and simulation verisimilitude, and can be widely applied to the technical field of seeker simulation.

Description

Technical Field

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

[0002] At present, the research and development of the seeker simulation system is mainly divided into signal level and function level. Signal level simulation is a high-resolution simulation method, which is committed to reflecting the signal flow process and processing process of zero intermediate frequency and below in the real system, and can realistically reflect the real system situation. However, signal level simulation has problems with simulation speed and memory space occupation, and it is difficult to adapt to the needs of large-scale system simulation or digital twins, especially when simulating broadband complex electromagnetic signal environments. The effect is very poor; functional level simulation is a low-resolution simulation method, which aims to realize system functions and does not require the true reflection of signal changes. It usually only involves some spatial energy calculations and probability-based threshold detection calculations, so the simulation speed is high and the memory occupancy is small. However, in many links, functional level simulation requires human subjective judgment to make decisions, and some parameters need to be set manually. Therefore, it is highly subjective, the simulation results are easily affected, and it is not very convincing. Therefore, it is often not used in system performance verification and evaluation. Summary of the invention

[0003] In view of this, the main purpose of the embodiments of the present invention is to provide a multi-resolution hybrid simulation method and device for a guide head, in order to solve at least one of the problems in the prior art. The present invention can improve the simulation speed and simulation realism.

[0004] To achieve the above object, an embodiment of the present invention provides a multi-resolution hybrid simulation method for a seeker, comprising the following steps:

[0005] Get the initial pulse description word data of the radar;

[0006] 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;

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

[0008] Acquire a space link loss and a receiving antenna gain, and generate a seeker receiving pulse description word data according to the transmitting antenna gain, the space link loss and the receiving antenna gain;

[0009] According to the pulse description word data received by the seeker, 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, obtaining a measurement pulse description word data stream;

[0011] Signal sorting is performed according to the measurement pulse description word data stream to obtain a sorting result.

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

[0013] Obtaining the relative distance between the radar and the seeker according to the latitude and longitude of the radar in the initial pulse description word data;

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

[0015] In some embodiments, judging whether the radar satisfies the radar interception condition according to the relative distance and the pulse arrival time, and obtaining the transmitting antenna gain if the condition is satisfied, and returning to the step of obtaining the initial pulse description word data of the radar if the condition is not satisfied, comprises the following steps:

[0016] According to the relative distance, constructing airspace interception conditions;

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

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

[0019] Constructing frequency domain interception conditions according to the frequency scanning range, the intermediate frequency bandwidth, the speed of switching frequency bands, and the pulse arrival time;

[0020] Determine whether the radar satisfies the air domain interception condition, the time domain interception condition and the frequency domain interception condition; if the radar satisfies the air domain interception condition, the time domain interception condition and the frequency domain interception condition, obtain the transmit antenna gain; if the radar does not satisfy the air domain interception condition, the time domain interception condition or the frequency domain interception condition, return to the step of obtaining the initial pulse description word data of the radar;

[0021] Among them, the radar interception condition includes the air domain interception condition, the time domain interception condition and the frequency domain interception condition.

[0022] In some embodiments, the obtaining of the spatial link loss and the receiving antenna gain comprises the following steps:

[0023] Obtain the wavelength of electromagnetic wave signals and the signal power attenuation caused by meteorological conditions;

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

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

[0026] Acquiring the relative position of the radar to the seeker;

[0027] Converting the relative orientation to obtain the incident direction of the radar signal;

[0028] Preset antenna pointing direction;

[0029] Obtaining the angle between the antenna pointing direction and the incident direction of the radar signal;

[0030] According to the angle, the receiving antenna gain is obtained through an antenna pattern function.

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

[0032] According to the pulse description word data received by the seeker, the adjacent time of a plurality of the signal pulses is sampled to obtain a plurality of pulse sampling signal streams and corresponding pulse sampling time streams;

[0033] Taking the current pulse sampling time flow of the current signal pulse as the total time flow, and taking the current pulse sampling signal flow of the current signal pulse as the total signal flow;

[0034] Obtaining the time relationship between the current signal pulse and the next signal pulse;

[0035] According to the total time flow, the total signal flow and the time relationship, the next signal pulse is fused with the current signal pulse to obtain a fused pulse, and the fused pulse is used as the current signal pulse. The step of returning to the step of using the current pulse sampling time flow of the current signal pulse as the total time flow and the current pulse sampling signal flow of the current signal pulse as the total signal flow is performed until all the signal pulses are fused to obtain the mixed intermediate frequency sampling signal flow.

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

[0037] When the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling time stream to the total time stream;

[0038] When the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling signal stream to the total signal stream;

[0039] When the time relationship is a partial aliasing of adjacent pulses, adding the next pulse sampling time stream to the total time stream, and removing repeated time values;

[0040] When the time relationship is partial aliasing of adjacent pulses, the end of the total signal stream is padded with 0 until the total signal stream is 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 exceeds;

[0041] When the time relationship is that adjacent pulses are completely overlapped, 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 moment.

[0042] In some embodiments, obtaining the measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream comprises the following steps:

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

[0044] Parameter measurement is performed on the composite pulse signal to generate the measured pulse description word data stream.

[0045] In some embodiments, the performing parameter measurement on the composite pulse signal to generate the measured pulse description word data stream comprises the following steps:

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

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

[0048] Obtaining the 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] According to the pulse width and the pulse arrival time, obtaining the center time of the composite pulse signal;

[0051] Taking the signal amplitude value at the central moment as the pulse amplitude;

[0052] Measure the signal carrier frequency through frequency domain threshold;

[0053] The arrival azimuth is measured by the four-antenna amplitude comparison direction finding method;

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

[0055] To achieve the above object, another aspect of an embodiment of the present invention provides a multi-resolution hybrid simulation device for a seeker, the device comprising:

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

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

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

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

[0060] The fifth module is used to perform fusion processing on a plurality of signal pulses according to the pulse description word data received by the seeker to obtain a mixed intermediate frequency sampling signal stream;

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

[0062] The seventh module is used to perform signal sorting according to the measurement pulse description word data stream to obtain a sorting result.

[0063] To achieve the above-mentioned purpose, another aspect of an embodiment of the present invention provides an electronic device, which includes a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, it implements the multi-resolution hybrid simulation method of a directional guide head described above.

[0064] To achieve the above objectives, another aspect of an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the multi-resolution hybrid simulation method for a directional guide head described above is implemented.

[0065] To achieve the above object, another aspect of the embodiment of the present invention provides a computer program product or a computer program, the computer program product or the computer program includes computer instructions, and the computer instructions are 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 directional guide head.

[0066] The embodiments of the present invention include at least the following beneficial effects: The present invention provides a multi-resolution hybrid simulation method and device for a seeker, which obtains the initial pulse description word data of the radar; obtains the relative distance between the radar and the seeker and the pulse arrival time according to the initial pulse description word data; determines whether the radar meets the radar interception condition according to the relative distance and the pulse arrival time, and obtains the transmitting antenna gain if it meets the condition, and returns to the step of obtaining the initial pulse description word data of the radar if it does not meet the condition; obtains the space link loss and the receiving antenna gain, and generates the seeker receiving pulse description word data according to the transmitting antenna gain, the space link loss and the receiving antenna gain; performs fusion processing on several signal pulses according to the seeker receiving pulse description word data to obtain a mixed intermediate frequency sampling signal stream; obtains the measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream; performs signal sorting according to the measurement pulse description word data stream to obtain the sorting result, which can improve the simulation speed and simulation fidelity. BRIEF DESCRIPTION OF THE DRAWINGS

[0067] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0068] Figure 1 is a flow chart of a multi-resolution hybrid simulation method for a seeker provided by an embodiment of the present invention;

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

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

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

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

[0073] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of the present invention, and they are only examples of devices and methods consistent with some aspects of the embodiments of the present invention as detailed in the attached claims.

[0074] It should be noted that, although the functional modules are divided in the system schematic diagram and the logical order is shown in the flow chart, in some cases, the steps shown or described may be performed in a different order than the module division in the system or the flow chart. The terms "first / S100" and "second / S200" in the specification and claims and the above-mentioned drawings may be used to describe various concepts in this article, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, without departing from the scope of the embodiment of the present invention, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the words "if" and "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determination".

[0075] The terms "at least one", "multiple", "each", "any", etc. used in the present invention, at least one includes one, two or more, multiple includes two or more, 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 those commonly understood by those skilled in the art to which the present invention belongs. The terms used herein are only for the purpose of describing the embodiments of the present invention and are not intended to limit the present invention.

[0077] Radar seeker, also known as radio seeker, is a detection device. There are three main ways to verify and evaluate the working efficiency of the seeker system: formula derivation, field experiment and computer simulation. Using pure mathematical formulas to conduct static and qualitative analysis of the system is low-cost, but it is often difficult to consider all factors carefully, and the results obtained are not convincing enough; field experiments are the most realistic means to reflect the performance of the seeker system, but they require high costs and long preparation time, and are limited by the real environment. It is difficult to simulate the impact of different battlefield environments on system performance, and the flexibility and reusability are poor, so they are generally only used in the final stage of system development and testing; in recent years, the performance of computers has increased rapidly, and computer simulation capabilities have been greatly enhanced. Computer simulation has gradually become the most important means of system performance verification and evaluation. Computer simulation consumes less resources, is flexible and convenient to use, and has strong reusability. In addition, thanks to the continuous development of modern electronic information equipment towards digitalization, the fidelity of computer simulation can also be guaranteed. Therefore, research on simulation modeling of seeker systems can effectively promote the development of research in related fields of seeker systems.

[0078] At present, the research and development of the seeker simulation system is mainly divided into signal level and function level. Signal level simulation can truly reflect the flow process and processing steps of the internal signal of the system, and the simulation results are convincing, but there are outstanding problems such as slow simulation speed and large memory usage, which makes it difficult to adapt to the needs of large-scale system simulation or digital twins, especially when simulating complex broadband electromagnetic signal environments. The effect is very poor; Functional level simulation is a low-resolution simulation method. It aims to realize the system function and does not require the true reflection of signal changes. It usually only involves some spatial energy calculations and probability-based threshold detection calculations. Therefore, the simulation speed is high and the memory usage is small. However, in many links of functional level simulation, decisions need to be made based on human subjective judgment and some parameters need to be set manually. Therefore, some parameter settings are interfered by subjective factors, and the simulation results are easily affected. The persuasiveness is insufficient, and it cannot reflect the real processing and transformation of the signal. Therefore, it is often not used in system performance verification and evaluation. It can be seen that signal level simulation and functional level simulation have their own advantages and disadvantages, but neither can effectively simulate the seeker system in a broadband electromagnetic environment.

[0079] In view of this, if Figure 1 As shown, an embodiment of the present invention provides a multi-resolution hybrid simulation method for a seeker, which may include but is not limited to steps S100 to S700:

[0080] Step S100, obtaining initial pulse description word data of the radar;

[0081] Step S200, obtaining the relative distance between the radar and the seeker and the pulse arrival time according to the initial pulse description word data;

[0082] Step S300, judging whether the radar meets the radar interception condition according to the relative distance and the pulse arrival time, and if so, obtaining the transmitting antenna gain, and if not, returning to the step of obtaining the initial pulse description word data of the radar;

[0083] Step S400, obtaining a space link loss and a receiving antenna gain, and generating a seeker receiving pulse description word data according to the transmitting antenna gain, the space link loss and the receiving antenna gain;

[0084] Step S500, performing fusion processing on a plurality of signal pulses according to the pulse description word data received by the seeker to obtain a mixed intermediate frequency sampling signal stream;

[0085] Step S600, obtaining a measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream;

[0086] Step S700, performing signal sorting according to the measurement pulse description word data stream to obtain a sorting result.

[0087] In step S100 of some embodiments, initial pulse description word data of the radar is generated according to the radar parameter information, and the initial pulse description word data may include but is not limited to antenna scanning type, scanning parameters, signal pulse modulation type, modulation parameters, radar signal transmission power, radar pulse width, repetition type, repetition parameters, longitude and latitude of the radar, altitude, radar transmitting antenna gain, radar half-power beam width, etc. Optionally, the radar parameter information may be preset by the user, loaded from a local file, or loaded from a server, which is not limited here.

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

[0089] Step S210, obtaining the relative distance between the radar and the seeker according to the latitude and longitude of the radar in the initial pulse description word data;

[0090] Step S220, obtaining the pulse arrival time according to the relative distance and the speed of light.

[0091] In step S210 of some embodiments, the longitude and latitude of the radar of the initial pulse description word data are obtained, and the relative distance R between the radar and the seeker can be calculated from the longitude and latitude. Then, the relative azimuth of the radar and the seeker can be calculated. Before calculating the relative azimuth, the coordinate system needs to be converted. Optionally, the coordinate system conversion process is: geodetic coordinate system (L, B, H) -> geocentric rectangular coordinate system (X, Y, Z) -> north celestial east coordinate system (X N ,YN ,Z N ). Then the relative azimuth from the radiation source radar to the seeker can be calculated in the North Celestial East coordinate system. The expression of the relative azimuth is:

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

[0093] Wherein, direction indicates the relative azimuth from the radar to the seeker.

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

[0095] TimeDelay=R / c

[0096] The pulse arrival time is calculated as:

[0097] TOA=TOE+TimeDelay

[0098] Wherein, TimeDelay represents the time delay from radar to seeker; R represents the relative distance between radar and 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 may include but is not limited to steps S310 to S350:

[0100] Step S310, constructing airspace interception conditions according to the relative distance;

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

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

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

[0104] Step S350, determining whether the radar satisfies the air domain interception condition, the time domain interception condition, and the frequency domain interception condition; if the radar satisfies the air domain interception condition, the time domain interception condition, and the frequency domain interception condition, obtaining the transmit antenna gain; if the radar does not satisfy the air domain interception condition, the time domain interception condition, or the frequency domain interception condition, returning to the step of obtaining the initial pulse description word data of the radar;

[0105] Among them, the radar interception condition includes the air 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 radiation source radar antenna are obtained, and the maximum line of sight of the seeker to the radiation source radar signal is calculated by the height of the seeker and the height of the radiation source radar antenna. By comparing the relative distance and the maximum line of sight, it can be determined whether the seeker is likely to receive the pulse signal emitted by the radar. Among them, the maximum line of sight refers to the maximum direct line of sight distance that the radar can detect the target due to the influence of the curvature of the earth's surface. When the relative distance from the radar to the seeker exceeds the maximum line of sight, the pulse signal emitted by the radar will be blocked by the ground. Exemplarily, the calculation formula for the maximum line of sight is:

[0107]

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

[0109] According to the relative distance and the maximum line of sight, the airspace interception condition can be constructed, which is as follows:

[0110]

[0111] When the relative distance is less than or equal to the maximum line of sight, it means that the airspace screening has passed, that is, the seeker may receive the pulse signal emitted by the radar; when the relative distance is greater than the maximum line of sight, it means that the airspace screening has not passed, that is, the seeker may not receive the pulse signal emitted by the radar.

[0112] In step S320 of some embodiments, the working start time and working range of the seeker can be preset, and the time domain capture condition can be constructed after the pulse arrival time is calculated. The pulse can be screened in the time domain according to the time domain capture condition to determine whether the pulse arrival time is within the working time range of the seeker. Exemplarily, the time domain capture condition that passes the time domain screening is:

[0113] t begin<TOA<t begin +duration

[0114] Among them, duration represents the working time of the seeker; t begin Indicates the start time of the seeker's operation.

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

[0116]

[0117] In the formula,

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

[0119] In step S350 of some embodiments, the radar is intercepted in airspace, time domain and frequency domain, i.e., it is determined whether the radar meets the airspace interception condition, time domain interception condition and frequency domain interception condition. If the radar does not meet one of the airspace interception condition, time domain interception condition or frequency domain interception condition, the step of obtaining the initial pulse description word data of the radar is returned, i.e., the step of generating the initial pulse description word data of the next radar is entered. If the radar meets the airspace interception condition, time domain interception condition and frequency domain interception condition, the next step is entered to obtain the transmitting antenna gain. Based on the signal transmission power stored in the initial pulse description word data, the transmitting antenna gain can be determined according to the antenna scanning type and specific parameters of the radar. The radar antenna rotates according to a certain rule during operation. The angle between the radar antenna pointing and the radar to the receiver determines the transmitting antenna gain of the radar signal received by the seeker. The transmitting antenna gain can be calculated by the antenna pattern function.

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

[0121]

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

[0123] The mathematical expression of the directional pattern outside the main lobe is:

[0124]

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

[0126] After obtaining the above pattern function, the antenna gain can be calculated from the angle between the radar antenna pointing and the radar to the receiver. The antenna gain is the maximum gain of the antenna multiplied by the antenna pattern function, that is:

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

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

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

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

[0131] Step S420, obtaining a 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;

[0134] Step S450, converting the relative orientation to obtain the incident direction of the radar signal;

[0135] Step S460, presetting the antenna direction;

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

[0137] Step S480: Obtain the receiving antenna gain according to the angle through an antenna pattern function.

[0138] In steps S410 to S430 of some embodiments, the space link loss can be calculated by pulse carrier frequency, signal propagation distance, meteorological conditions and specific parameters. The space link loss is obtained by adding the distance loss caused by signal propagation and the signal power attenuation caused by meteorological conditions. During the propagation of electromagnetic signals in the atmosphere, losses will occur due to the influence of meteorology, and the losses in various meteorological conditions can be obtained by corresponding empirical formulas. Exemplarily, the calculation formula for distance loss is:

[0139]

[0140] Where λ is the wavelength of the electromagnetic wave signal; R is the relative distance between the radar and the seeker.

[0141] In some embodiments, in step S440 to step S480, the receiving antenna gain of the corresponding antenna channel is calculated. The incident direction of the radar signal can be obtained by converting the relative orientation of the radar to the seeker, and the conversion formula is:

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

[0143] Where, θ is the incident direction of the radar signal; θ0 is the relative position from the radar to the seeker; % is the remainder operation.

[0144] In the simulation system, it can be assumed that the four antennas point at 0°, 90°, 180°, and 270°, and that they are all Gaussian antennas. The receiving antenna gain can be determined based on the angle between the incident direction of the radar signal and the antenna. Assume that the antenna points at Then the angle between the antenna pointing and the incident direction of the radar signal is:

[0145]

[0146] After the angle is determined, the receiving antenna gain can be determined using the same method as that used to obtain the transmitting antenna gain.

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

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

[0149] In the formula, PDW represents the pulse description word data received by the seeker; TOA represents the pulse arrival time; PW represents the pulse width; Pr represents the pulse receiving power; RF represents the carrier frequency corresponding to the pulse; mod type represents the pulse intra-pulse modulation type; modpara represents the pulse intra-pulse modulation parameters.

[0150] In step S500 of some embodiments, a function-level simulation-signal-level simulation conversion is performed, and a simulation acceleration method for broadband electromagnetic environments is proposed. In view of the problem that computer memory is unable to withstand large-bandwidth, full-time intermediate frequency sampling signal flows, it is proposed to segment the full-time intermediate frequency sampling signal flow based on a time window and only simulate the intermediate frequency sampling signal flow near the pulse signal.

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

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

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

[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 flow, the total signal flow and the time relationship, the next signal pulse is fused with the current signal pulse to obtain a fused pulse, the fused pulse is used as the current signal pulse, and the step of returning to the step of using the current pulse sampling time flow of the current signal pulse as the total time flow and the current pulse sampling signal flow of the current signal pulse as the total signal flow is returned, until all the signal pulses are fused to obtain the mixed intermediate frequency sampling signal flow.

[0156] In step S510 of some embodiments, according to the pulse description word data received by each seeker, 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 around 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] Among them, 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 individual pulse sampling time streams; [SampleSignal-1], [SampleSignal-2], [SampleSignal-3]... are individual pulse sampling signal streams.

[0160] In step S520 of some embodiments, the basic idea of ​​fusion of intermediate frequency sampling signal streams is to use the sampling time stream and sampling signal stream of the first signal pulse as the total time stream and total signal stream, and fuse the subsequent signal pulses one by one. Optionally, the current pulse sampling time stream of the current signal pulse is used as the total time stream, and the current pulse sampling signal stream of the current signal pulse is used as the total signal stream, so as to lay a good foundation for the subsequent fusion of the signal pulses one by one.

[0161] In step S530 of some embodiments, since the pulse description word information is sorted by pulse arrival time in the generation of the pulse description word data received by the seeker by default, it can be ensured that the start time of the current signal pulse is before the start time of the subsequent signal pulse. Then the time relationship between the current signal pulse and the next signal pulse has the following three types: the current signal pulse and the next signal pulse are not overlapped, the current signal pulse and the next signal pulse are partially overlapped, and the next signal pulse is completely overlapped in the current signal pulse. Figure 2 The adjacent pulse situations shown include non-aliasing of adjacent pulses, partial aliasing of adjacent pulses, and complete aliasing of adjacent pulses.

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

[0163] Step S541, when the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling time stream to the total time stream;

[0164] Step S542, when the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling signal stream to the total signal stream;

[0165] Step S543, when the time relationship is partial aliasing of adjacent pulses, adding the next pulse sampling time stream to the total time stream, and removing repeated time values;

[0166] Step S544, when the time relationship is partial aliasing of adjacent pulses, pad the end of the total signal stream with 0 until the total signal stream is equal to the length of the total time stream, perform a traversal operation on the total signal stream from back to front, and add the next pulse sampling signal stream to the total signal stream until the length of the next pulse sampling signal stream exceeds;

[0167] Step S545, when the time relationship is that adjacent pulses are completely overlapped, 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 moment.

[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 overlapped, 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 some embodiments, in steps S543 to S544, 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, and adjacent pulses are partially overlapped, the total time stream is added to the next pulse sampling time stream and the repeated time value is removed, so that 0 is added to the end of the total signal stream to make it equal to the total time stream. Finally, starting from the end of the total signal stream, the sampling signal stream of the next signal pulse is traversed from back to front and added to the total signal stream until it exceeds the length of the next pulse sampling signal stream.

[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 overlapped at this time, then the total time flow is kept unchanged, the sampled signal flow of the next signal pulse is traversed from back to front, and the signal value of the sampled signal flow of the next signal pulse is added to the signal value of the total signal flow at the corresponding moment. The corresponding moment here refers to the corresponding part of the sampled time flow of the next signal pulse and the total time flow, so that the signal value of the next signal pulse can be superimposed with the signal value of the total signal flow corresponding to the corresponding total time flow.

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

[0172] In some embodiments, step S600 may 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 measured pulse description word data stream.

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

[0176] In step S620 of some embodiments, the time domain, frequency domain, and spatial domain parameter measurements are performed on the mixed intermediate frequency sampling signal stream of the radar radiation source in the composite pulse signal through a signal detection and measurement model. By digital signal processing of the mixed intermediate frequency sampling signal of the radar radiation source, parameters such as signal arrival time (TOA), pulse width (PW), pulse amplitude (PA), signal carrier frequency (RF), and arrival direction angle (DOA) are measured to form a measured pulse description word data stream, laying a solid foundation for subsequent sorting and identification processing.

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

[0178] Step S621, performing threshold detection on the composite pulse signal by envelope detection method to obtain a detection result and a detection threshold;

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

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

[0181] Step S624, obtaining a 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 moment as the pulse amplitude;

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

[0185] Step S628, measuring the arrival azimuth by using a four-antenna amplitude comparison direction finding method;

[0186] The measured 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 some embodiments, in step S621 to step S622, the signal arrival time is measured, referring to Figure 3 , Figure 3 The schematic diagram of the seeker measuring the signal arrival time. Exemplarily, the input composite pulse signal is threshold detected by the envelope detection method to obtain the detection result and the detection threshold. The detection result is compared with the detection threshold. If the detection result is greater than the detection threshold, it means that the composite pulse signal is detected, and the corresponding moment of the rising edge 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 previous pulse arrival time to form a new pulse arrival time and output as a parameter. Then the time counter is cleared and the counting starts again until the next pulse rising edge arrives.

[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 use a counter for timing, but the timing of the falling edge PFE (Pulse Falling Edge) is added. The calculation formula is:

[0189]

[0190] Where 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 quantized value of the time between the falling edge time PFE and the pulse arrival time TOA in units of Δt.

[0191] In steps S625 to S626 of some embodiments, the pulse amplitude is measured. The measurement of the pulse amplitude PA utilizes the measured pulse arrival time and pulse width. Exemplarily, the center moment of the composite pulse signal is calculated based on the measured pulse arrival time and pulse width, and the signal amplitude value at the center moment is used as the measured value of the pulse amplitude. Optionally, the calculation formula of the center moment is:

[0192]

[0193] Where, T center 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 in the form of a frequency domain threshold, and the radar signal arrival azimuth DOA is measured using a four-antenna amplitude comparison direction finding method.

[0195] In step S700 of some embodiments, radar signals can be sorted according to the generated measured pulse description word data stream to obtain sorting results. Exemplarily, 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 perform preliminary sorting of the radar signal. The purpose of pre-sorting is to dilute the pulse stream so that 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 pulse by processing the pulse arrival time TOA. By setting the detection threshold, the PRI transformation value exceeding the threshold can be identified, thereby estimating 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, radar signals can be effectively classified and identified.

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

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

[0198] like Figure 4 As shown, taking the composite simulation system equipped with the multi-resolution hybrid simulation method for the face guide head according to the embodiment of the present invention as an example, hybrid simulation is performed, and the following steps are involved:

[0199] Step 1: Generate the pulse description word data of the radar end according to the radar parameter information obtained by user settings or file loading, including antenna scanning type, scanning parameters, signal pulse modulation type, modulation parameters, radar signal transmission power, radar pulse width, repetition type, repetition parameters, radar latitude and longitude, altitude, radar transmitting antenna gain, radar half-power beam width and other information.

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

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

[0202] Step 4: Intercept the radar in air domain, time domain, and frequency domain. The specific steps are as follows:

[0203] Step 1: Perform airspace capture: After calculating the relative distance from the radiation source radar to the seeker, airspace screening can be performed. By comparing the relative distance and the maximum line of sight of the seeker to the radiation source radar signal, it is determined whether the seeker can receive the pulse signal emitted by the radar. If not, the subsequent process of generating the radar pulse description word is skipped directly, and the next radar pulse description word generation step is entered.

[0204] Step 2: Time domain capture: In the simulation system, the user sets the start time and working range of the seeker. After calculating the pulse arrival time, the pulse can be screened in the time domain to determine whether the pulse arrival time is within the seeker working time range. If the pulse arrival time is not within the seeker working time range, the subsequent process of generating the radar pulse description word is skipped directly, and the next radar pulse description word generation step is entered.

[0205] Step 3: Perform frequency domain capture: obtain the frequency scanning range, intermediate frequency bandwidth, and switching speed of the seeker receiver. According to the frequency scanning range, intermediate frequency bandwidth, switching speed, and pulse arrival time, construct the frequency domain capture conditions to determine whether the frequency domain capture conditions are met. If not, directly skip the subsequent process of generating the radar pulse description word and enter the next radar pulse description word generation step.

[0206] Step 5: Based on the signal transmission power stored in the pulse descriptor, 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 through the antenna pattern function.

[0207] Step 6: Add the distance loss caused by signal propagation and the signal power attenuation caused by meteorological conditions to obtain the spatial link loss.

[0208] Step 7: Calculate the receiving antenna gain of the corresponding antenna channel. In the simulation system, it is assumed that the four antennas point at 0°, 90°, 180° and 270°, and all are Gaussian antennas. Determine the receiving antenna gain based on the angle between the incident direction of the radar signal and the antenna. After determining the angle, determine the receiving antenna gain using the same method as the transmitting antenna gain.

[0209] Step 8: Generate the seeker receiving pulse descriptor data to obtain the signal pulse descriptor structure arriving at a certain reconnaissance antenna channel.

[0210] Step 9: Perform function-level simulation-signal-level simulation conversion. Generate the corresponding intermediate frequency sampling signal stream near the pulse according to each pulse description word information. After generating multiple sampling signal streams near the signal pulse according to the pulse description word information, perform fusion processing on the sampling signal streams. By continuously fusing the sampling time streams and sampling signal streams of adjacent pulses, a total time stream and signal stream can eventually be formed as a radiation source mixed intermediate frequency sampling signal stream.

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

[0212] Step 11: Measure parameters in the time domain, frequency domain, and spatial domain. 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. By digital signal processing of the mixed intermediate frequency signal of the radar radiation source, the signal arrival time, pulse width, pulse amplitude, signal carrier frequency, and arrival direction angle and other parameters are measured to form a measurement pulse description word data stream for subsequent sorting and identification processing. The specific steps are as follows:

[0213] Step 1: Measure the signal arrival time: perform envelope detection on the input signal, compare the detection result with the detection threshold, if it is greater than the threshold, it is considered that a pulse is detected, and the corresponding moment of the pulse rising edge is taken as the pulse arrival time. Each time the pulse rising edge arrives, the counting result of the time counter is accumulated with the previous pulse arrival time to form a new pulse arrival time and output as a parameter. Then clear the time counter and start counting again until the next pulse rising edge arrives.

[0214] Step 2: Measure the pulse width: The estimation of pulse width is similar to the estimation of pulse arrival time. Both use a counter for timing, but the timing of the falling edge is added. The pulse width is a quantized value of the time interval from the falling edge time to the rising edge time of the pulse, with the unit being the minimum quantized value.

[0215] Step 3: Measure the pulse amplitude: Measure the pulse amplitude according to the measured pulse arrival time and pulse width, calculate the center moment of the pulse, and use the signal amplitude at the center moment as the measured value of the pulse amplitude.

[0216] Step 4: Use the frequency domain threshold to measure the carrier frequency of the frequency domain parameter signal.

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

[0218] Step 12: Sorting and identification. Sorting and identification are divided into two steps, namely pre-sorting and main sorting. The formed measurement pulse description word data stream is pre-sorted by the OPTICS clustering algorithm, and the pre-sorted measurement pulse description word data stream is mainly sorted by the PRI transformation method, and the sorting results are output.

[0219] The embodiment of the present invention further provides a multi-resolution hybrid simulation device for a directional guide head, which can implement the multi-resolution hybrid simulation method for a directional guide head. The device includes:

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

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

[0222] The third module is used to determine whether the radar meets the radar interception condition according to the relative distance and the pulse arrival time, and if so, obtain the transmitting antenna gain; if not, return to the step of obtaining the initial pulse description word data of the radar;

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

[0224] The fifth module is used to perform fusion processing on a plurality of signal pulses according to the pulse description word data received by the seeker to obtain a mixed intermediate frequency sampling signal stream;

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

[0226] The seventh module is used to perform signal sorting according to the measurement pulse description word data stream to obtain a sorting result.

[0227] It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0228] The embodiment of the present invention further provides an electronic device, which includes a processor and a memory, wherein the memory stores a computer program, and the processor implements the above-mentioned multi-resolution hybrid simulation method for a directional guide head when executing the computer program. The electronic device can be any intelligent terminal including a tablet computer, a vehicle-mounted computer, etc.

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

[0230] refer to Figure 5 , Figure 5 The hardware structure of an electronic device of another embodiment is illustrated, and the electronic device includes:

[0231] The processor 801 may be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of the present invention.

[0232] The memory 802 can be implemented in the form of a read-only memory (ROM), a static storage device, a dynamic storage device or a random access memory (RAM). The memory 802 can store an operating system and other application programs. When the technical solution provided in the embodiment of this specification is implemented by software or firmware, the relevant program code is stored in the memory 802, and the processor 801 calls and executes a multi-resolution hybrid simulation method for a seeker according to an embodiment of the present invention;

[0233] Input / output interface 803, used to implement information input and output;

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

[0235] A bus 805 that transmits information between the various components of the device (e.g., the processor 801, the memory 802, the input / output interface 803, and the communication interface 804);

[0236] The processor 801 , the memory 802 , the input / output interface 803 and the communication interface 804 are connected to each other in communication within the device via a bus 805 .

[0237] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the multi-resolution hybrid simulation method for a directional guide head is implemented.

[0238] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0239] The embodiment of the present invention also provides a computer program product or a computer program, which includes 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 directional guide head.

[0240] A multi-resolution hybrid simulation method and device for a seeker according to an embodiment of the present invention has the following advantages:

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

[0242] 2. The traditional functional-level modeling of complex electromagnetic signal environments directly generates pulse description word data at the receiving end of the reconnaissance receiver, which cannot reflect the impact of signal propagation in space. The embodiment of the present invention improves the method for generating pulse description word data for functional-level modeling of complex electromagnetic signal environments. The construction of pulse description word data starts from the radar transmitting end, and the possible impacts on the signal transmission and propagation process such as signal type, intra-pulse modulation mode, carrier frequency information, transmitting antenna gain, and space link loss are taken into account, so that the pulse description word data that finally reaches the seeker receiver is more realistic.

[0243] 3. The embodiment of the present invention divides the full-time intermediate frequency sampling signal stream based on the time window, only simulates the intermediate frequency sampling signal stream near the pulse signal, and then fuses the aliased pulses. This effectively reduces the simulation memory usage without affecting the simulation results, and also speeds up the subsequent signal processing links.

[0244] In summary, in the process of modeling a complex signal environment, the embodiment of the present invention uses function-level simulation to improve the method for generating pulse descriptor data, starts the construction of pulse descriptors from the radar transmitter, and reflects the changes of signals during the transmission, propagation, and reception processes in the pulse descriptors through calculation, so that the pulse descriptor data arriving at the receiving end is more realistic; in the process of modeling the seeker, signal-level simulation is used to study the core modules of the seeker: signal detection and parameter measurement module, signal sorting module, etc., and selects an algorithm suitable for a broadband electromagnetic signal environment for simulation modeling; in response to the interface problem of simulation models of different granularities, the embodiment of the present invention uses a signal flow generation method based on a pulse descriptor-sampled signal mapping relationship to solve the problem of effective information transmission between modules of different granularities.

[0245] In some selectable embodiments, the function / operation mentioned in the block diagram may not occur in the order mentioned in the operation diagram. For example, depending on the function / operation involved, the two boxes shown in succession can actually be executed substantially simultaneously or the boxes can sometimes be executed in reverse order. In addition, the embodiment presented and described in the flow chart of the present invention is provided by way of example, for the purpose of providing a more comprehensive understanding of technology. The disclosed method is not limited to the operation and logic flow presented herein. Selectable embodiments are expected, wherein the order of various operations is changed and the sub-operation of a part for which is described as a larger operation is performed independently.

[0246] In addition, although the present invention is described in the context of functional modules, it should be understood that, unless otherwise specified, one or more of the functions and / or features described may be integrated into a single physical device and / or software module, or one or more functions and / or features may be implemented in separate physical devices or software modules. It is also understood that a detailed discussion of the actual implementation of each module is unnecessary for understanding the present invention. More specifically, in view of the properties, functions, and internal relationships of the various functional modules in the device disclosed herein, the actual implementation of the module will be understood within the conventional skills of the engineer. Therefore, those skilled in the art can implement the present invention set forth in the claims without excessive experimentation using ordinary techniques. It is also understood that the specific concepts disclosed are merely illustrative and are not intended to limit the scope of the present invention, which is determined by the full scope of the appended claims and their equivalents.

[0247] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc., which can store program codes.

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

[0249] More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk case (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be a paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, deciphering or, if necessary, processing in another suitable manner, and then stored in a computer memory.

[0250] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0251] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0252] Although the embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

[0253] The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the described embodiments. Those skilled in the art may make various equivalent modifications or substitutions without violating the spirit of the present invention. These equivalent modifications or substitutions are all included in the scope defined by the claims of the present invention.

Claims

1. A multi-resolution hybrid simulation method for a seeker, characterized in that: The following steps are involved: Get the initial pulse description word data of the radar; 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; According to the relative distance and the pulse arrival time, judging whether the radar meets the radar interception condition, if so, obtaining the transmitting antenna gain, if not, returning to the step of obtaining the initial pulse description word data of the radar; Acquire a space link loss and a receiving antenna gain, and generate a seeker receiving pulse description word data according to the transmitting antenna gain, the space link loss and the receiving antenna gain; According to the pulse description word data received by the seeker, 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, obtaining a measurement pulse description word data stream; Signal sorting is performed according to the measurement pulse description word data stream to obtain a sorting result.

2. The multi-resolution hybrid simulation method for a directional guide head according to claim 1, characterized in that: The method of obtaining the relative distance between the radar and the seeker and the pulse arrival time according to the initial pulse description 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 description word data; The pulse arrival time is obtained according to the relative distance and the speed of light.

3. The multi-resolution hybrid simulation method for a directional guide head according to claim 1, characterized in that: The step of judging whether the radar meets the radar interception condition according to the relative distance and the pulse arrival time, obtaining the transmitting antenna gain if the condition is met, and returning to the step of obtaining the initial pulse description word data of the radar if the condition is not met, comprises the following steps: According to the relative distance, constructing airspace interception conditions; According to the pulse arrival time, constructing a time domain interception condition; Obtaining the frequency scanning range, intermediate frequency bandwidth and switching frequency band speed of the seeker; Constructing frequency domain interception conditions according to the frequency scanning range, the intermediate frequency bandwidth, the speed of switching frequency bands, and the pulse arrival time; Determine whether the radar satisfies the air domain interception condition, the time domain interception condition and the frequency domain interception condition; if the radar satisfies the air domain interception condition, the time domain interception condition and the frequency domain interception condition, obtain the transmit antenna gain; if the radar does not satisfy the air domain interception condition, the time domain interception condition or the frequency domain interception condition, return to the step of obtaining the initial pulse description word data of the radar; Among them, the radar interception condition includes the air domain interception condition, the time domain interception condition and the frequency domain interception condition.

4. The multi-resolution hybrid simulation method for a directional guide head according to claim 1, characterized in that: The obtaining of the spatial link loss and the receiving antenna gain comprises the following steps: Obtain the wavelength of electromagnetic wave signals and the signal power attenuation caused by meteorological conditions; Obtaining a path loss according to the relative distance and the wavelength of the electromagnetic wave signal; Obtaining the spatial link loss according to the signal power attenuation and the path loss; Acquiring the relative position of the radar to the seeker; Converting the relative orientation to obtain the incident direction of the radar signal; Preset antenna pointing direction; Obtaining the angle between the antenna pointing direction and the incident direction of the radar signal; According to the angle, the receiving antenna gain is obtained through an antenna pattern function.

5. The multi-resolution hybrid simulation method for a directional guide head according to claim 1, characterized in that: The method of receiving pulse description word data from the seeker and fusing a plurality of signal pulses to obtain a mixed intermediate frequency sampling signal stream comprises the following steps: According to the pulse description word data received by the seeker, the adjacent time of a plurality of the signal pulses is sampled to obtain a plurality of pulse sampling signal streams and corresponding pulse sampling time streams; Taking the current pulse sampling time flow of the current signal pulse as the total time flow, and taking the current pulse sampling signal flow of the current signal pulse as the total signal flow; Obtaining the time relationship between the current signal pulse and the next signal pulse; According to the total time flow, the total signal flow and the time relationship, the next signal pulse is fused with the current signal pulse to obtain a fused pulse, and the fused pulse is used as the current signal pulse. The step of returning to the step of using the current pulse sampling time flow of the current signal pulse as the total time flow and the current pulse sampling signal flow of the current signal pulse as the total signal flow is performed until all the signal pulses are fused to obtain the mixed intermediate frequency sampling signal flow.

6. The multi-resolution hybrid simulation method for a directional guide head according to claim 5, characterized in that: The fusing process of the next signal pulse with the current signal pulse according to the total time flow, the total signal flow and the time relationship comprises the following steps: When the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling time stream to the total time stream; When the time relationship is that adjacent pulses are not aliased, adding the next pulse sampling signal stream to the total signal stream; When the time relationship is a partial aliasing of adjacent pulses, adding the next pulse sampling time stream to the total time stream, and removing repeated time values; When the time relationship is partial aliasing of adjacent pulses, the end of the total signal stream is padded with 0 until the total signal stream is 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 exceeds; When the time relationship is that adjacent pulses are completely overlapped, 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 moment.

7. The multi-resolution hybrid simulation method for a directional guide head according to claim 1, characterized in that: The step of obtaining a measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream comprises the following steps: Superimposing the mixed intermediate frequency sampling signal stream and the thermal noise signal stream to generate a composite pulse signal; Parameter measurement is performed on the composite pulse signal to generate the measured pulse description word data stream.

8. The multi-resolution hybrid simulation method for a directional guide head according to claim 7, characterized in that: The step of measuring parameters of the composite pulse signal and generating the measured pulse description word data stream comprises the following steps: Performing threshold detection on the composite pulse signal by envelope detection method to obtain detection result and detection threshold; Comparing the detection result with the detection threshold, if the detection result is greater than the detection threshold, taking the time corresponding to the rising edge of the composite pulse signal as the pulse arrival time; Obtaining the falling edge time of the composite pulse signal; Obtaining a pulse width according to the falling edge time and the pulse arrival time; According to the pulse width and the pulse arrival time, obtaining the center time of the composite pulse signal; Taking the signal amplitude value at the central moment as the pulse amplitude; Measure the signal carrier frequency through frequency domain threshold; The arrival azimuth is measured by the four-antenna amplitude comparison direction finding method; The measured 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 multi-resolution hybrid simulation device for a guide head, characterized in that: include: The first module is used to obtain the initial pulse description word data of the radar; The second module is used to obtain the relative distance between the radar and the seeker and the pulse arrival time according to the initial pulse description word data; The third module is used to determine whether the radar meets the radar interception condition according to the relative distance and the pulse arrival time, and if so, obtain the transmitting antenna gain; if not, return to the step of obtaining the initial pulse description word data of the radar; The fourth module is used to obtain the space link loss and the receiving antenna gain, and generate the seeker receiving pulse description word data according to the transmitting antenna gain, the space link loss and the receiving antenna gain; The fifth module is used to perform fusion processing on a plurality of signal pulses according to the pulse description word data received by the seeker to obtain a mixed intermediate frequency sampling signal stream; A sixth module is used to obtain a measurement pulse description word data stream according to the mixed intermediate frequency sampling signal stream; The seventh module is used to perform signal sorting according to the measurement pulse description word data stream to obtain a sorting result.

10. An electronic device, characterized in that: including a processor and a memory; The memory is used to store programs; The processor executes the program to implement the method according to any one of claims 1 to 8.

Images