Radiation Source Simulation Method, Device, Storage Medium, and Program Product in a Dynamic Scenario

Abstract

Embodiments of the present disclosure disclose a method, apparatus, storage medium, and program product for simulating radiation sources in a dynamic scenario. Among them, the method includes: constructing a dynamic simulation scenario, setting a plurality of radiation sources in the dynamic simulation scenario and configuring radiation source parameters, where the radiation source parameters include pulse repetition interval and simulation step size; obtaining the frame period of the radiation source based on the pulse repetition interval; sampling and discretizing the simulation time of the radiation source based on the simulation step size to obtain a time discrete sequence; advancing the simulation according to the time discrete sequence, and for each advancement of a simulation step size, sampling and discretizing the current simulation step size based on the frame period to obtain the total number of loop calculations; obtaining the remaining time of the current simulation step size based on the total number of loop calculations; compensating the remaining time of the current simulation step size into the next simulation step size until the simulation ends. This method can solve the problem that the simulation time cannot be evenly divided by the sampling discretization of the pulse repetition interval through the remaining time compensation mechanism.

Description

Technical Field

[0001] The present disclosure relates to the technical field of electromagnetic environment simulation, and particularly to a method, device, storage medium, and program product for simulating radiation sources in a dynamic scenario. Background Art

[0002] In digital signal processing and simulation technologies, dealing with the problem that the simulation time is not divisible by the sampling discretization of the pulse repetition interval (PRI) is a common challenge. The traditional method, namely the truncation method, is widely adopted due to its simplicity. When dealing with non - divisible simulation time, this method directly discards the remaining part after division. Although this method is easy to implement, it causes a small period of time (less than the PRI part) to be missing in each simulation step. As the simulation steps accumulate, more and more time is discarded, resulting in a large number of pulses being lost, and then causing the disorder of the pulse sequence and the mismatch between the pulse amplitude and the device scanning period, ultimately causing errors to the measured device. Summary of the Invention

[0003] In view of this, embodiments of the present disclosure provide a method, device, storage medium, and program product for simulating radiation sources in a dynamic scenario, which can solve the problem that the simulation time is not divisible by the sampling discretization of the pulse repetition interval through a remaining time compensation mechanism.

[0004] In a first aspect, embodiments of the present disclosure provide a method for simulating radiation sources in a dynamic scenario, adopting the following technical solutions:

[0005] Construct a dynamic simulation scenario, set a number of radiation sources in the dynamic simulation scenario and configure radiation source parameters, where the radiation source parameters include the pulse repetition interval and the simulation step size;

[0006] Based on the pulse repetition interval, obtain the frame period of the radiation source;

[0007] Based on the simulation step size, sample and discretize the simulation time of the radiation source to obtain a time - discrete sequence;

[0008] Advance the simulation according to the time - discrete sequence. Each time a simulation step size is advanced, sample and discretize the current simulation step size based on the frame period to obtain the total number of loop calculations;

[0009] Based on the total number of loop calculations, obtain the remaining time of the current simulation step;

[0010] Compensate the remaining time of the current simulation step into the next simulation step until the simulation ends.

[0011] Optionally, the obtaining the frame period of the radiation source based on the pulse repetition interval includes:

[0012] If the radiation source parameters only include one pulse repetition interval, the frame period of the radiation source is equal to the pulse repetition interval;

[0013] If the radiation source parameters include multiple different pulse repetition intervals, the frame period of the radiation source is equal to the sum of the multiple different pulse repetition intervals added together.

[0014] Optionally, the calculation formula for the total number of loop calculations is:

[0015]

[0016] where i is the sequence number of the radiation source; x is the sequence number of the current simulation step; x - 1 is the sequence number of the previous simulation step; N ix is the total number of loop calculations of the i-th radiation source in the current simulation step; floor(·) is the floor function; Ts ix is the current simulation step of the i-th radiation source; PRI i(x-1)Left is the remaining time of the i-th radiation source in the previous simulation step.

[0017] Optionally, the calculation formula for the remaining time of the current simulation step is:

[0018] PRI ixLeft =(Ts ix +PRI i(x-1)Left ) - N ix ×PRI iFrame ;

[0019] where PRI ixLeft is the remaining time of the i-th radiation source in the current simulation step.

[0020] Optionally, the radiation source simulation method in the dynamic scenario further includes:

[0021] In each loop calculation, obtain the current PDW parameters;

[0022] where the current PDW parameters include the index of the current pulse repetition interval, the current pulse repetition interval, the pulse width, the signal modulation type, the pulse arrival time, the received carrier frequency, and the pulse amplitude.

[0023] Optionally, if the radiation source parameters only include one pulse repetition interval, the index of the current pulse repetition interval is 0;

[0024] If the radiation source parameters include multiple different pulse repetition intervals, obtain the total number of pulse repetition interval stagger or pulse repetition interval slip;

[0025] Obtain the index of the current pulse repetition interval based on the total number, where the calculation formula for the index of the current pulse repetition interval is:

[0026]

[0027] where i is the sequence number of the radiation source; x is the sequence number of the current simulation step; x - 1 is the sequence number of the previous simulation step; is the index of the current pulse repetition interval of the i-th radiation source at the current simulation step and the current loop calculation times; Index i(x-1) is the index value of the i-th radiation source after the end of the previous simulation step; n ix is the current loop calculation times of the i-th radiation source at the current simulation step, n ix ≤N ix ; N ix is the total number of loop calculations of the i-th radiation source at the current simulation step; % is the remainder function; M i is the total number of pulse repetition interval stagger or pulse repetition interval sliding of the i-th radiation source.

[0028] Optionally, the radiation source simulation method in the dynamic scenario further includes:

[0029] Construct the PDW parameter in each loop calculation into a PDW six-tuple matrix;

[0030] Combine the PDW six-tuple matrices of each radiation source into a PDW comprehensive matrix;

[0031] Output the PDW comprehensive matrix of each simulation step to a csv file for saving.

[0032] In a second aspect, the embodiments of the present disclosure further provide a radiation source simulation system in a dynamic scenario, adopting the following technical solutions:

[0033] A simulation scenario construction module, configured to construct a dynamic simulation scenario, set a plurality of radiation sources in the dynamic simulation scenario, and configure radiation source parameters, where the radiation source parameters include a pulse repetition interval and a simulation step;

[0034] A frame period acquisition module, configured to obtain the frame period of the radiation source based on the pulse repetition interval;

[0035] A discrete sequence acquisition module, configured to sample and discretize the simulation time of the radiation source based on the simulation step to obtain a time discrete sequence;

[0036] A total number acquisition module, configured to advance the simulation according to the time discrete sequence. For each advancement of a simulation step, sample and discretize the current simulation step based on the frame period to obtain the total number of loop calculations;

[0037] A remaining time acquisition module, configured to acquire the remaining time of the current simulation step based on the total number of loop calculations;

[0038] A simulation time compensation module, configured to compensate the remaining time of the current simulation step into the next simulation step until the simulation ends.

[0039] In a third aspect, an embodiment of the present disclosure further provides a computer device, adopting the following technical solution:

[0040] The computer device includes:

[0041] At least one processor; and,

[0042] A memory communicatively connected to the at least one processor; wherein,

[0043] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute any one of the above-mentioned radiation source simulation methods in a dynamic scenario.

[0044] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores computer instructions for causing a computer to execute any one of the above-mentioned radiation source simulation methods in a dynamic scenario.

[0045] In a fifth aspect, an embodiment of the present disclosure further provides a computer program product, including a computer program / instructions, and when the computer program / instructions are executed by a processor, the steps of any one of the above-mentioned methods are implemented.

[0046] The radiation source simulation method in a dynamic scenario provided by the embodiments of the present disclosure sets a plurality of radiation sources in a dynamic simulation scenario and configures corresponding parameters for each radiation source. This flexibility enables the simulation to more comprehensively simulate complex situations in the actual environment, including the interaction and influence between different radiation sources. By obtaining the frame period of the radiation source based on the pulse repetition interval and combining with the simulation step for time discretization, the dynamic changes of the radiation source at different time points are simulated, making the dynamic simulation scenario more real and closer to the actual situation, providing a basic condition for subsequent calculation of the remaining time. By advancing the simulation according to the time discrete sequence and performing sampling discretization again based on the frame period within each simulation step, it is convenient to compare the simulation step and the frame period, obtain the total number of loop calculations, and then obtain the remaining time of the current simulation step, and compensate this remaining time into the next simulation step. This compensation mechanism reduces the simulation error caused by the non-divisibility of the simulation step and the PRI, helps to maintain the continuity and stability of the simulation, and

[0047] Improve the simulation efficiency further, reduce the simulation error caused by time discretization, and enhance the reliability of the simulation.

[0048] The above description is only an overview of the technical solution of the present disclosure. In order to understand the technical means of the present disclosure more clearly, it can be implemented according to the content of the specification. And in order to make the above and other purposes, features and advantages of the present disclosure more obvious and understandable, the following preferred embodiments are specifically given and described in detail in conjunction with the accompanying drawings as follows. Brief Description of the Drawings

[0049] In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the accompanying drawings required for use in the embodiments will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some embodiments of the present disclosure. For those of ordinary skill in the art, other accompanying drawings can be obtained based on these drawings without creative efforts.

[0050] Figure 1 It is a schematic flowchart of the radiation source simulation method in a dynamic scenario provided by an embodiment of the present disclosure;

[0051] Figure 2 It is a schematic diagram of the time discrete sequence provided by an embodiment of the present disclosure;

[0052] Figure 3 It is a schematic flowchart of the time discrete sequence after resampling and discretization provided by an embodiment of the present disclosure;

[0053] Figure 4 It is another schematic flowchart of the time discrete sequence after resampling and discretization provided by an embodiment of the present disclosure;

[0054] Figure 5 It is yet another schematic flowchart of the time discrete sequence after resampling and discretization provided by an embodiment of the present disclosure;

[0055] Figure 6 It is a schematic diagram of the parameter values in the sgb file provided by an embodiment of the present disclosure;

[0056] Figure 7 It is a schematic diagram of the parameter values in the pulse file provided by an embodiment of the present disclosure;

[0057] Figure 8 It is a schematic diagram of the parameter values in the csv file provided by an embodiment of the present disclosure;

[0058] Figure 9 It is a schematic block diagram of the radiation source simulation system in a dynamic scenario provided by an embodiment of the present disclosure;

[0059] Figure 10 It is a schematic structural diagram of a computer device provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0060] The embodiments of the present disclosure are described in detail below with reference to the accompanying drawings.

[0061] It should be clear that the following embodiments of the present disclosure are described by specific specific examples, and those skilled in the art can easily understand other advantages and effects of the present disclosure from the contents disclosed in this specification. Obviously, the described embodiments are only a part of the embodiments of the present disclosure, rather than all of the embodiments. The present disclosure can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present disclosure. It should be noted that the following embodiments and features in the embodiments can be combined with each other in the absence of conflict. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of the present disclosure.

[0062] It should be noted that various aspects of the embodiments within the scope of the appended claims are described below. It should be apparent that the aspects described herein may be embodied in a wide variety of forms, and any specific structures and / or functions described herein are merely illustrative. Based on this disclosure, those skilled in the art should understand that an aspect described herein may be implemented independently of any other aspect, and two or more of these aspects may be combined in various ways. For example, any number of the aspects described herein may be used to implement the device and / or practice the method. In addition, one or more of the aspects described herein may be used.

[0063] Other structure and / or functionality besides those described above may be used to implement this apparatus and / or practice this method.

[0064] It should also be noted that the illustrations provided in the following embodiments are only schematic illustrations of the basic concept of the present disclosure. The drawings only show components related to the present disclosure rather than being drawn according to the number, shape and size of components in actual implementation. In actual implementation, the type, quantity and proportion of each component may be changed arbitrarily, and the component layout may also be more complicated.

[0065] Additionally, in the following description, specific details are provided to facilitate a thorough understanding of the examples. However, it will be understood by those skilled in the art that the aspects described may be practiced without these specific details.

[0066] Reference Figure 1 The present disclosure provides a radiation source simulation method in a dynamic scene, comprising the following steps:

[0067] S1: Construct a dynamic simulation scenario, set a number of radiation sources in the dynamic simulation scenario, and configure the radiation source parameters, where the radiation source parameters include pulse repetition interval and simulation step size;

[0068] S2: Based on the pulse repetition interval, obtain the frame period of the radiation source;

[0069] S3: Sample and discretize the simulation time of the radiation source based on the simulation step size to obtain a time discrete sequence;

[0070] S4: Advance the simulation according to the time discrete sequence. For each advancement of a simulation step size, sample and discretize the current simulation step size based on the frame period to obtain the total number of loop calculations;

[0071] S5: Based on the total number of loop calculations, obtain the remaining time of the current simulation step size;

[0072] S6: Compensate the remaining time of the current simulation step size into the next simulation step size until the simulation ends.

[0073] The radiation source simulation method provided by the present disclosure sets a number of radiation sources in the dynamic simulation scenario and configures corresponding parameters for each radiation source. This flexibility enables

[0074] the simulation to more comprehensively simulate complex situations in the actual environment, including the interactions and influences between different radiation sources. By obtaining the frame period of the radiation source based on the pulse repetition interval and combining the simulation step size for time discretization, the dynamic changes of the radiation source at different time points are simulated, making the dynamic simulation scenario more realistic and closer to the actual situation, providing a basic condition for subsequent calculation of the remaining time. By advancing the simulation according to the time discrete sequence and performing sampling and discretization again based on the frame period within each simulation step size, it is convenient to compare the simulation step size and the frame period, obtain the total number of loop calculations, and then obtain the remaining time of the current simulation step size. Compensating this remaining time into the next simulation step size, this compensation mechanism reduces the simulation error caused by the non - divisibility of the simulation step size and the PRI, helps to maintain the continuity and stability of the simulation, further improves the simulation efficiency, reduces the simulation error caused by time discretization, and improves the reliability of the simulation.

[0075] In summary, the radiation source simulation method in this dynamic scenario realizes high - precision and high - efficiency simulation of the behavior of radiation sources in a complex dynamic environment through refined time processing, flexible parameter configuration, optimized computing resource allocation, and strong scene adaptability.

[0076] In S1, use radar simulation tools (such as WRSim (Waveform Radar Simulator), MATLAB / Simulink, STK (Systems Tool Kit), etc.) to construct a dynamic scene for radar simulation (referred to as a dynamic simulation scene). Set several (for example, 3) radiation sources at different three-dimensional coordinate positions, and configure a motion trajectory and radiation source parameters for each radiation source. The radiation source parameters include the Pulse Repetition Interval (abbreviated as PRI). The Pulse Repetition Interval refers to the time interval between two consecutive pulses emitted by the radar.

[0077] In S2, calculate the frame period based on the Pulse Repetition Interval. The frame period refers to the time interval for continuously emitting a group of pulses (referred to as a frame). Macroscopically, the frame period is the smallest and repeatedly occurring Pulse Repetition Interval (PRI) period.

[0078] If the Pulse Repetition Interval of the radiation source is a fixed value (abbreviated as fixed PRI), that is, the radiation source parameters only include one Pulse Repetition Interval, then the frame period is equal to the Pulse Repetition Interval, that is, PRI iFrame =PRI i =PRI 固定 where PRI iFrame is the frame period of the i-th radiation source, PRI i is the Pulse Repetition Interval of the i-th radiation source, i is the sequence number of the radiation source, 1 ≤ i ≤ I, I is the total number of radiation sources, and PRI 固定 is a fixed value.

[0079] If the Pulse Repetition Interval of the radiation source is not a fixed value (abbreviated as non-fixed PRI), that is, the radiation source adopts the PRI staggering or PRI sliding mode, and the set radiation source parameters will include multiple different Pulse Repetition Intervals, then the frame period of the radiation source is equal to the sum of multiple different Pulse Repetition Intervals, that is, where PRI ij is the j-th Pulse Repetition Interval of the i-th radiation source, j is the sequence number of the Pulse Repetition Interval, 1 ≤ j ≤ J, and J is the total number of different Pulse Repetition Intervals of the i-th radiation source. In the above calculation process, the unit of the Pulse Repetition Interval is unified to microseconds (μs). For example, a certain radiation source adopts the PRI staggering mode. In this mode, a total of 3 different Pulse Repetition Intervals are set, with values of 200 μs, 211 μs, and 217 μs respectively. Add the 3 values to get a frame period of 628 μs.

[0080] In S3, the radiation source parameters also include the simulation step, which refers to the minimum time unit for the simulation system to advance in time, that is, the time interval for the simulation system to update the state, perform calculations, and process events each time.

[0081] Each radiation source emits signals according to its own motion trajectory and radiation source parameters for radar simulation. During the simulation process, the simulation step samples and discretizes the simulation time of the radiation source to obtain a time-discrete sequence. Refer to Figure 2 the schematic diagram of the time-discrete sequence shown. The simulation step divides the overall simulation time into multiple time slices, where the simulation step is fixed at Ts.

[0082] In S4, the time-discrete sequence is a set of time points divided at fixed time intervals during the simulation process. Each time point represents a simulation step, and the time-discrete sequence is used to control the way of time advancement during the simulation to ensure the consistency and accuracy of the simulation step. According to the time-discrete sequence, the simulation time point is initialized, and the simulation time is advanced according to the time-discrete sequence. Each time a simulation step is advanced, the current simulation step is sampled and discretized based on the frame period. Furthermore, the time-discrete sequence is sampled and discretized again by the frame period. At this time, three situations will occur, as follows:

[0083] In the first situation, the current simulation step is greater than or equal to the frame period, and the current simulation step can be divided evenly by the frame period. Refer to Figure 3 the schematic diagram of the process of the time-discrete sequence after resampling and discretization shown. When advancing to the first and second simulation steps of the i-th radiation source, after being intercepted by multiple frame periods, there is no remaining time.

[0084] In the second situation, the current simulation step is greater than the frame period, but the current simulation step cannot be divided evenly by the frame period, and there will be some remaining time after the current simulation step is intercepted by one or more frame periods. Refer to Figure 4 another schematic diagram of the process of the time-discrete sequence after resampling and discretization shown. After the first simulation step of the i-th radiation source is intercepted by multiple frame periods, there is still some remaining time, and after the second simulation step of the i-th radiation source is intercepted by multiple frame periods, there is also some remaining time. Among them, the remaining time is Figure 4 the gray part in.

[0085] In the third situation, the current simulation step is less than the frame period, and the current simulation step cannot be divided evenly by the frame period. Refer to Figure 5 the schematic diagram of the process of the time-discrete sequence after resampling and discretization shown. Figure 5 The gray part in is the time when the frame period exceeds the current simulation step.

[0086] The total number of loop calculations for the current simulation step is calculated using the following formula in all of the above three cases:

[0087]

[0088] In Formula 1, x is the sequence number of the current simulation step; x - 1 is the sequence number of the previous simulation step; N ix is the total number of loop calculations for the i-th radiation source at the current simulation step; floor(·) is the floor function, which is used to round a real number down to the nearest integer that is not greater than the original number; Ts ix is the current simulation step of the i-th radiation source; PRI i(x-1)Left is the remaining time of the i-th radiation source at the previous simulation step. When x = 1, PRI i(x-1)Left = 0, that is, the initial value of PRI i(x-1)Left is 0.

[0089] In S5, for the total number of loop calculations for the current simulation step, the remaining time of the current simulation step is calculated using the following formula:

[0090] PRI ixLeft = (Ts ix + PRI i(x-1)Left ) - N ix ×PRI iFrame ; (Formula 2)

[0091] In Formula 2, PRI ixLeft is the remaining time of the i-th radiation source at the current simulation step.

[0092] Specifically, in the above first case (i.e., the current simulation step is greater than or equal to the frame period and the current simulation step is divisible by the frame period), there is no remaining time for each simulation step. Therefore, N ix > 0 and PRI ixLeft = 0. In the above second case (i.e., the current simulation step is greater than the frame period but the current simulation step is not divisible by the frame period), there is remaining time for each simulation step, and there may be a divisible situation. Therefore, N ix > 0 and 0 ≤ PRI ixLeft ≤ PRI iFrame ≤ Ts ix . In the above third case (i.e., the current simulation step is less than the frame period and the current simulation step is not divisible by the frame period), there is remaining time for each simulation step, and there may be a divisible situation. Therefore, N ix ≥ 0 and 0 ≤ PRI ixLeft ≤ Ts ix < PRI iFrame .

[0093] In summary, the above formula 1 and formula 2 can be used for loop calculation in all these three cases,

[0094] except that the total number of loop calculations and the values of the remaining time are different.

[0095] In S6, update the remaining time of the current simulation step to the next simulation step, and the next simulation step is Ts i(x+1) , Ts i(x+1) = Ts ix + PRI ixLeft , when advancing to the next simulation step according to the time discrete sequence, x = x + 1, update the next simulation step to the current simulation step and continue the iteration until the simulation ends.

[0096] Among them, the conditions for the end of the simulation include: First, the simulation time of all radiation sources reaches the corresponding predetermined time threshold; Second, the simulation of all radiation sources reaches the corresponding preset target or performance index; Third, the simulation of all radiation sources is manually interrupted by the user (such as through keys, commands, etc.); Fourth, all radiation sources trigger the corresponding preset events during the simulation. When any of the above conditions is met, it is determined that the simulation ends.

[0097] When performing radiation source simulation in a dynamic scenario, the remaining time is generated in each simulation step due to the change of PRI. If these remaining times are not processed, it may lead to simulation inaccuracy because they may affect the emission timing and periodicity of the radiation source signal. To solve this problem, this application adopts a compensation mechanism, that is, accumulates the time that cannot be fully consumed in the current simulation step (i.e., the remaining time) and compensates it to the next simulation step. This method solves the problem that the remaining time is discarded when the simulation step and PRI cannot be divided evenly, thus maintaining the continuity and accuracy of the simulation, ensuring that the signal emission of the radiation source is consistent with the actual change of PRI. This method is applicable not only to the simulation mode with fixed PRI, but also to the simulation modes with PRI stagger and PRI slip.

[0098] Furthermore, the radiation source parameters also include pulse width (Pulse Width, abbreviated as PW) and signal modulation type (Modulation Operation Parameter, abbreviated as MOP). Among them, the pulse width refers to the duration of the radar transmitted pulse; the signal modulation type refers to the method

[0099] used to control and modulate the pulse signal.

[0100] When N ix ≠ 0, perform N operations on the current simulation step ixThe calculation of the nth cycle. In each cycle calculation, the Pulse Description Words (abbreviated as PDW parameters, also known as PDW streams) of each radiation source are obtained respectively. For any one of the radiation sources, its current PDW parameters are obtained. The current PDW parameters include the index of the current Pulse Repetition Interval (Index), the current Pulse Repetition Interval, the Pulse Width, the Signal Modulation Type, the Time of Arrival (abbreviated as TOA), the Received Carrier Frequency (abbreviated as RCF), and the Pulse Amplitude (abbreviated as PA).

[0101] Among them, the current Pulse Repetition Interval, the Pulse Width, and the Signal Modulation Type can be directly obtained according to the radiation source parameters. During the simulation process, each radiation source will cycle through the Pulse Repetition Intervals. If the Pulse Repetition Interval of the radiation source is a fixed value, the value of the current Pulse Repetition Interval will always be this fixed value. If the Pulse Repetition Interval of the radiation source is not a fixed value, the radiation source will cycle through multiple different Pulse Repetition Intervals, and the value of the current Pulse Repetition Interval will change according to a certain rule. The index refers to the index of the current PRI in the entire PRI sequence. This index provides a unique identifier for each different PRI, which is crucial for distinguishing and managing the PRI change patterns from different radiation sources; the Time of Arrival of the pulse refers to the exact time when the pulse signal reaches the receiver after being emitted from the radiation source; the Received Carrier Frequency refers to the carrier frequency of the signal detected by the receiver; the Pulse Amplitude refers to the peak intensity of the pulse signal.

[0102] If it is a fixed PRI, the index of the current Pulse Repetition Interval is always 0; if it is a non-fixed PRI, the total number of multiple different PRI values within the frame period is obtained, that is, the total number of different PRI values within a PRI modulation period, specifically referring to the total number of PRI stagger or PRI slip. Based on obtaining the index of the current Pulse Repetition Interval, where the formula for the index of the current Pulse Repetition Interval is as follows:

[0103] In the formula 3,

[0104]

[0105] In formula 3, is the index of the current Pulse Repetition Interval of the ith radiation source at the current simulation step and the current cycle calculation times; Index i(x-1) is the index value of the ith radiation source after the end of the previous simulation step; n ix is the current cycle calculation times of the ith radiation source at the current simulation step, n ix ≤N ix; % is the remainder function (also known as modulo operation), used to determine the remainder after one number is divided by another; M i is the total number of pulse repetition interval stagger or pulse repetition interval sliding for the i-th radiation source.

[0106] Collect the pulse transmission time (Time of Transmission, abbreviated as TOT, unit: second), emission-reception delay (Emission-Reception Delay, abbreviated as Delay, unit: second), Doppler frequency caused by emission-reception movement (Doppler Frequency Shift, abbreviated as DopFre, unit: hertz), and received signal strength (Received Signal Strength, abbreviated as RSS, unit: decibel-milliwatt) of each radiation source. Refer to Figure 6 the schematic diagram of parameter values in the sgb file shown, and output the pulse transmission time, emission-reception delay, Doppler frequency caused by emission-reception movement, received signal strength, and index of each radiation source to the corresponding sgb file respectively.

[0107] Collect the signal carrier frequency (Carrier Frequency, abbreviated as CF, unit: hertz), bandwidth (Bandwidth, abbreviated as BW, unit: hertz), and initial phase (Initial Phase, abbreviated as PhI, unit: rad) of each radiation source. Refer to Figure 7 the schematic diagram of parameter values in the pulse file shown, and output the signal modulation type, signal carrier frequency, bandwidth, initial phase, pulse repetition interval, and pulse width of each radiation source to the corresponding pulse file respectively.

[0108] Read the sgb file and pulse file output by each radiation source respectively, and obtain the pulse arrival time, received carrier frequency, and pulse amplitude based on the parameters in the two files. Among them, the calculation formula for the pulse arrival time is as follows:

[0109]

[0110] In formula 4, is the pulse arrival time of the i-th radiation source at the current simulation step and the current loop calculation times; is the pulse transmission time of the i-th radiation source at the current simulation step and the current loop calculation times; is the signal propagation delay of the i-th radiation source at the current simulation step and the current loop calculation times, is the distance between the i-th radiation source and the receiving end at the current simulation step and the current loop calculation times, with the unit of meter; c is the speed of light, taking 299792458.0.

[0111] The calculation formula for the received carrier frequency is as follows:

[0112]

[0113] In Equation (5), is the carrier frequency received by the i-th radiation source at the current simulation step and the current loop calculation times; is the signal carrier frequency of the i-th radiation source at the current simulation step and the current loop calculation times; is the Doppler frequency caused by the transceiver movement of the i-th radiation source at the current simulation step and the current loop calculation times, is the Doppler velocity of the i-th radiation source at the current simulation step and the current loop calculation times.

[0114] The calculation formula for the pulse amplitude is as follows:

[0115]

[0116] In Equation (6), is the pulse amplitude reaching the receiving end of the i-th radiation source at the current simulation step and the current loop calculation times; is the equivalent radiated power of the i-th radiation source at the current simulation step and the current loop calculation times; is the space propagation loss of the electromagnetic wave from the transmitting end to the receiving end of the i-th radiation source at the current simulation step and the current loop calculation times, is the receiving antenna gain of the i-th radiation source at the current simulation step and the current loop calculation times; is the system loss of the i-th radiation source at the current simulation step and the current loop calculation times.

[0117] Obtain the PDW parameters for each loop calculation in the above manner, sort them in ascending order of TOA, and construct the PDW parameters in each loop calculation into a PDW six-tuple matrix (i.e., the PDW stream of a single radiation source at the current simulation step). The expression of the PDW six-tuple matrix is as follows:

[0118]

[0119] In Equation (7), is the pulse arrival time of the i-th radiation source at the last loop calculation of the current simulation step, is the carrier frequency received by the i-th radiation source in the last cycle calculation of the current simulation step; is the pulse repetition interval of the i-th radiation source at the current simulation step and the current cycle calculation times; is the pulse repetition interval of the i-th radiation source in the last cycle calculation of the current simulation step; is the pulse width of the i-th radiation source at the current simulation step and the current cycle calculation times; is the pulse width of the i-th radiation source in the last cycle calculation of the current simulation step; is the pulse amplitude of the i-th radiation source in the last cycle calculation of the current simulation step; is the signal modulation type of the i-th radiation source at the current simulation step and the current cycle calculation times; The signal modulation type of the i-th radiation source in the last cycle calculation of the current simulation step.

[0120] Sorted by TOA from small to large, the PDW six-tuple matrix of each radiation source is combined into a PDW comprehensive matrix (i.e., the PDW stream of dynamic multi-radiation sources at the current simulation step). The expression of the PDW comprehensive matrix is as follows:

[0121]

[0122] In formula (8), y is the current sorting value after sorting the TOAs obtained by all radiation sources at the current simulation step from small to large. y ≤ Y, where Y is the maximum sorting value at the current simulation step. According to y, other parameters of the corresponding radiation source can be determined in the same batch of cycle calculations.

[0123] Refer to Figure 8 the schematic diagram of parameter values in the shown csv file. When the simulation ends, the PDW streams of dynamic multi-radiation sources at each simulation step are output to the csv file for storage.

[0124] In summary, this method can not only solve the problem of simulation time loss in fixed PRI, PRI stagger, PRI sliding and other modes, ensure the integrity of pulses in the whole simulation process and the matching of the change law of pulse amplitude with the timing of the device scanning period, greatly improve the readiness rate of the measured device, but also systematically record and save the detailed PDW stream data in the simulation process of dynamic multi-radiation sources. By outputting these data to the CSV file, it is not only convenient for subsequent data analysis and processing, but also improves the traceability and verifiability of the simulation results, thus providing an efficient, flexible and reliable data management method for electromagnetic environment simulation.

[0125] Refer to Figure 9 , the present disclosure provides a radiation source simulation system in a dynamic scenario, including:

[0126] The simulation scenario construction module 101 is used to construct a dynamic simulation scenario, set a plurality of radiation sources in the dynamic simulation scenario and configure the radiation source parameters, where the radiation source parameters include the pulse repetition interval and the simulation step size;

[0127] The frame period acquisition module 102 is used to acquire the frame period of the radiation source based on the pulse repetition interval;

[0128] The discrete sequence acquisition module 103 is used to sample and discretize the simulation time of the radiation source based on the simulation step size to obtain a time discrete sequence;

[0129] The total number of times acquisition module 104 is used to advance the simulation according to the time discrete sequence. For each advancement of one simulation step size, sample and discretize the current simulation step size based on the frame period to obtain the total number of loop calculations;

[0130] The remaining time acquisition module 105 is used to acquire the remaining time of the current simulation step size based on the total number of loop calculations;

[0131] The simulation time compensation module 106 is used to compensate the remaining time of the current simulation step size into the next simulation step size until the simulation ends.

[0132] The various change methods and specific examples in the above-provided radiation source simulation method under the dynamic scenario are equally applicable to the radiation source simulation system under the dynamic scenario provided by the present disclosure. Through the foregoing detailed description of the radiation source simulation method under the dynamic scenario, those skilled in the art can clearly know the implementation method of the radiation source simulation system under the dynamic scenario. For the sake of brevity of the specification, it will not be elaborated herein.

[0133] The computer device according to an embodiment of the present disclosure includes a memory and a processor. The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0134] The processor may be a central processing unit (CPU) or other forms of processing units with data processing capabilities and / or instruction execution capabilities, and may control other components in the computer device to execute the desired functions. In an embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory, so that the computer device executes all or part of the steps of the radiation source simulation method under the dynamic scenario of the foregoing embodiments of the present disclosure.

[0135] Those skilled in the art should be able to understand that, in order to solve the technical problem of how to obtain good user experience effects, the embodiments may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included in the protection scope of the present disclosure.

[0136] Such as Figure 10 FIG. is a schematic structural diagram of a computer device provided by an embodiment of the present disclosure. It shows a schematic structural diagram of a computer device suitable for implementing the computer device in the embodiments of the present disclosure. Figure 10 The shown computer device is only an example and should not impose any limitation on the functions and usage scope of the embodiments of the present disclosure.

[0137] Such as Figure 10 As shown, the computer device may include a processor (such as a central processing unit, a graphics processing unit, etc.), which may perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) or a program loaded from a storage device into a random access memory (RAM). In the RAM, various programs and data required for the operation of the computer device are also stored. The processor, ROM, and RAM are connected to each other through a bus. An input / output (I / O) interface is also connected to the bus.

[0138] Generally, the following devices may be connected to the I / O interface: an input device including, for example, a sensor or a visual information acquisition device; an output device including, for example, a display screen; including, for example

[0139] a storage device such as a magnetic tape or a hard disk; and a communication device. The communication device may allow the computer device to communicate wirelessly or wirelessly with other devices (such as edge computing devices) to exchange data. Although Figure 10 a computer device with various devices is shown, it should be understood that it is not required to implement or have all the shown devices. Instead, more or fewer devices may be implemented or had.

[0140] Specifically, according to the embodiments of the present disclosure, the process described above with reference to the flowchart may be implemented as a computer software program. For example, the embodiments of the present disclosure include a computer program product, which includes a computer program carried on a non-transitory computer-readable medium, and the computer program includes program codes for executing the method shown in the flowchart. In such an embodiment, the computer program may be downloaded and installed from the network through the communication device, or installed from the storage device, or installed from the ROM. When the computer program is executed by the processor, all or part of the steps of the radiation source simulation method in the dynamic scenario of the embodiments of the present disclosure are executed.

[0141] For a detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which will not be repeated herein.

[0142] A computer-readable storage medium according to an embodiment of the present disclosure stores non-transitory computer-readable instructions thereon. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the method for simulating a radiation source in a dynamic scenario of each embodiment of the present disclosure described above are executed.

[0143] The above computer-readable storage media include, but are not limited to: optical storage media (such as CD-ROMs and DVDs), magneto-optical storage media (such as MOs), magnetic storage media (such as magnetic tapes or external hard drives), media with built-in rewritable non-volatile memories (such as

[0144] memory cards) and media with built-in ROMs (such as ROM cartridges).

[0145] For a detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which will not be repeated herein.

[0146] The basic principles of the present disclosure have been described above in conjunction with specific embodiments. However, it should be noted that the advantages, benefits, effects, etc. mentioned in the present disclosure are only examples and not limitations, and it cannot be considered that these advantages, benefits, effects, etc. are essential for each embodiment of the present disclosure. In addition, the above-disclosed specific details are only for the purposes of illustration and facilitation of understanding, and not for limitation. The above details do not limit the present disclosure to necessarily adopt the above specific details for implementation.

[0147] In the present disclosure, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. The block diagrams of the devices, apparatuses, equipment, and systems involved in the present disclosure are only illustrative examples and do not intend to require or imply that they must be connected, arranged, and configured in the manner shown in the block diagrams. As those skilled in the art will recognize, these devices, apparatuses, equipment, and systems can be connected, arranged, and configured in any manner. Words such as "including", "comprising", "having", etc. are open-ended words, meaning "including but not limited to", and can be used interchangeably with them. The words "or" and "and" used herein refer to the word "and / or", and can be used interchangeably with it, unless the context clearly indicates otherwise. The word "such as" used herein refers to the phrase "such as but not limited to", and can be used interchangeably with it.

[0148] In addition, as used herein, the "or" used in the listing of items starting with "at least one" indicates a separate listing, so that for example, the listing of "at least one of A, B, or C" means A

[0149] or B or C, or AB or AC or BC, or ABC (i.e., A and B and C). In addition, the term "exemplary" does not mean that the examples described are preferred or better than other examples.

[0150] It should also be noted that in the systems and methods of the present disclosure, each component or each step can be decomposed and / or recombined. These decompositions and / or recombinations should be regarded as equivalent solutions of the present disclosure.

[0151] Various changes, substitutions, and alterations to the technology described herein can be made without departing from the teachings defined by the appended claims. In addition, the scope of the claims of the present disclosure is not limited to the specific aspects of the processes, machines, manufactures, compositions of events, means, methods, and acts described above. Current or later-developed processes, machines, manufactures, compositions of events, means, methods, or acts that perform substantially the same function or achieve substantially the same result as the corresponding aspects described herein can be utilized. Thus, the appended claims include such processes, machines, manufactures, compositions of events, means, methods, or acts within their scope.

[0152] The above description of the disclosed aspects is provided to enable any person skilled in the art to make or use the present disclosure. Various modifications to these aspects will be readily apparent to those skilled in the art, and the general principles defined herein can be applied to other aspects without departing from the scope of the present disclosure. Therefore, the present disclosure is not intended to be limited to the aspects shown herein, but rather to the broadest scope consistent with the principles and novel features disclosed herein.

[0153] The above description has been presented for purposes of illustration and description. In addition, this description is not intended to limit the embodiments of the present disclosure to the forms disclosed herein. Although several example aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations thereof.

Claims

1. A radiation source simulation method in a dynamic scene, characterized in that: include: Constructing a dynamic simulation scene, setting a plurality of radiation sources in the dynamic simulation scene and configuring radiation source parameters, wherein the radiation source parameters include a pulse repetition interval and a simulation step size; Based on the pulse repetition interval, acquiring a frame period of the radiation source; Based on the simulation step size, the simulation time of the radiation source is sampled and discretized to obtain a time discrete sequence; The simulation is advanced according to the time discrete sequence, and each time a simulation step is advanced, the current simulation step is sampled and discretized based on the frame period to obtain the total number of cyclic calculations; The calculation formula for the total number of cycles is: Where i is the sequence number of the radiation source; x is the sequence number of the current simulation step; x-1 is the sequence number of the previous simulation step; N ix is the total number of loop calculations of the ith radiation source in the current simulation step; floor(·) is the rounding function; Ts ix is the current simulation step of the ith radiation source; PRI i(x-1)Left is the remaining time of the ith radiation source in the previous simulation step; PRI iFrame is the frame period of the i-th radiation source; Based on the total number of loop calculations, the remaining time of the current simulation step is obtained; The calculation formula for the remaining time of the current simulation step is: PRI ixLeft =(Ts ix +PRI i(x-1)Left )-N ix ×PRI iFrame ; Among them, PRI ixLeft is the remaining time of the i-th radiation source in the current simulation step; The remaining time of the current simulation step is compensated to the next simulation step until the simulation ends.

2. The radiation source simulation method in a dynamic scene according to claim 1, characterized in that: The acquiring the frame period of the radiation source based on the pulse repetition interval comprises: If the radiation source parameters include only one pulse repetition interval, the frame period of the radiation source is equal to the pulse repetition interval; If the radiation source parameters include a plurality of different pulse repetition intervals, the frame period of the radiation source is equal to the sum of the plurality of different pulse repetition intervals.

3. The radiation source simulation method in a dynamic scene according to claim 1, characterized in that: Also includes: In each cycle calculation, obtain the current PDW parameters; The current PDW parameters include the index of the current pulse repetition interval, the current pulse repetition interval, the pulse width, the signal modulation type, the pulse arrival time, the received carrier frequency, and the pulse amplitude.

4. The radiation source simulation method in a dynamic scene according to claim 3, characterized in that: If the radiation source parameters include only one pulse repetition interval, the index of the current pulse repetition interval is 0; If the radiation source parameters include a plurality of different pulse repetition intervals, obtaining a total number of pulse repetition interval variations or pulse repetition interval slips; The index of the current pulse repetition interval is obtained based on the total number, wherein the calculation formula of the index of the current pulse repetition interval is: Where i is the sequence number of the radiation source; x is the sequence number of the current simulation step; x-1 is the sequence number of the previous simulation step; Index is the index of the current pulse repetition interval of the i-th radiation source under the current simulation step and the current number of loop calculations; i(x-1) is the index value of the ith radiation source after the previous simulation step ends; n ix is the current cycle calculation number of the i-th radiation source at the current simulation step, n ix ≤N ix ; N ix is the total number of loop calculations for the ith radiation source at the current simulation step; % is the remainder function; M i is the total number of pulse repetition interval variations or pulse repetition interval slips of the i-th radiation source.

5. The radiation source simulation method in a dynamic scene according to claim 3, characterized in that: Also includes: The PDW parameters in each cycle calculation are constructed as a PDW six-tuple matrix; Combine the PDW six-tuple matrix of each radiation source into a PDW comprehensive matrix; The PDW comprehensive matrix of each simulation step is exported to a csv file for saving.

6. A computer device, characterized in that: The computer device comprises: at least one processor; and, a memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the radiation source simulation method in a dynamic scene described in any one of claims 1-5.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the radiation source simulation method in a dynamic scene as described in any one of claims 1-5.

8. A computer program product comprising computer instructions, characterized in that: When the computer instruction is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.

Images