A laser jamming full-link digital parallel simulation system and method

Abstract

The present application belongs to the field of photoelectric countermeasure simulation and digital twin technology, and relates to a laser jamming full-link digital parallel simulation system and method. The system comprises a scene initialization unit, a laser parameter configuration unit, a target characteristic definition unit, an atmospheric transmission calculation unit, a load imaging simulation unit, a jamming effect evaluation unit and a visualization terminal connected in sequence. The method realizes digital simulation of the whole process of laser jamming through the whole process coupling of scene configuration, atmospheric transmission calculation, load imaging simulation, jamming effect evaluation and visualization output. The present application solves the problems of lack of full-link integration, insufficient dynamic process simulation, weak semi-physical verification capability and poor scalability in the prior art by constructing a full-link integrated simulation architecture covering ground radiation, atmospheric transmission, load imaging and image quality evaluation, and can be widely applied to photoelectric countermeasure equipment effectiveness evaluation, countermeasure strategy optimization and equipment performance verification.

Description

Technical Field

[0001] This invention belongs to the field of optoelectronic countermeasures simulation and digital twin technology, specifically relating to a laser interference full-link digital parallel simulation system and method. Background Technology

[0002] Laser jamming technology uses high-energy lasers to irradiate optical imaging systems (such as visible light CCD / CMOS cameras, infrared thermal imagers, and laser seekers), causing detector saturation, crosstalk, thermal damage, or even permanent failure, thereby effectively suppressing optoelectronic reconnaissance and guidance equipment.

[0003] Laser interference involves multiple stages, including ground-based radiation, atmospheric transmission, detector photoelectric response, and image processing, exhibiting typical multi-physics coupling characteristics. Existing simulation systems and related technologies suffer from the following main shortcomings: (1) It focuses on a single link, such as atmospheric transmission or detector damage, and lacks an integrated simulation platform that covers the entire link of ground radiation, atmospheric transmission, detector response and image quality assessment. The models of each link are independent and cannot truly reflect the complete physical process of laser interference. (2) The simulation results are static and difficult to use for dynamic countermeasure analysis because the temporal randomness of atmospheric turbulence, the transient characteristics of laser pulses and the time-varying characteristics of detector thermal accumulation are not considered. (3) The lack of hardware interfaces with physical equipment such as lasers, detector arrays and imaging systems makes it impossible to achieve closed-loop verification between simulation and actual measurement, and the accuracy of the model is difficult to correct through actual measurement. (4) The model parameters are fixed and it is difficult to adapt to the flexible configuration of different laser wavelengths, different detector types and different meteorological conditions, resulting in insufficient versatility.

[0004] Existing research on laser interference simulation falls into three categories: First, laser spot simulation based on atmospheric transmission, which uses atmospheric radiative transmission models such as MODTRAN and LOWTRAN to calculate laser energy attenuation and diffusion, but does not couple with detector response models, thus failing to generate realistic photoelectric imaging interference images; second, detector damage simulation based on physical models, which simulates the saturation and damage effects of lasers on CCD / CMOS, but does not consider atmospheric turbulence and background scene radiation, resulting in significant deviations between simulation results and real-world scenarios; and third, flare superposition methods based on image processing, which directly superimpose Gaussian spot templates onto the original image. While simple to operate, these methods lack physical mechanism support, cannot reflect the dynamic impact of parameter changes on interference effects, and are difficult to simulate nonlinear effects.

[0005] Therefore, there is an urgent need to build a full-link digital parallel simulation platform that can dynamically, accurately, and repeatably simulate the entire process of laser interference, so as to systematically support optoelectronic countermeasures simulation and performance evaluation and solve the above-mentioned defects of existing technologies. Summary of the Invention

[0006] To address the aforementioned technical problems, the present invention provides a solution. The present invention offers a laser interference end-to-end digital parallel simulation system and method.

[0007] In a first aspect, the present invention provides a laser interference full-link digital parallel simulation system, including a scene initialization unit, a laser parameter configuration unit, a target characteristic definition unit, an atmospheric transmission calculation unit, a load imaging simulation unit, an interference effect evaluation unit, and a visualization terminal. The scene initialization unit is used to import geographic information data, bind material reflection properties, and set the corresponding meteorological conditions for the simulation. The laser parameter configuration unit is used to select the laser interference mode, configure the laser wavelength and power parameters, and define the geometric parameters of laser emission. The target characteristic definition unit is used to specify the motion trajectory of the simulated target, set the thermal radiation parameters of the simulated target, and load the corresponding countermeasures for the simulated target. The atmospheric transmission calculation unit is used to calculate laser energy attenuation, generate turbulent phase screen and iterative angular spectrum propagation based on input parameters, and calculate the spot energy distribution of laser after atmospheric transmission. The payload imaging simulation unit is used to complete the modeling of the modulation transfer function of the optical system, the simulation of the photoelectric response of the detector (detector) and the injection of imaging noise based on the energy distribution of the light spot after atmospheric transmission, and to generate photoelectric imaging results under laser interference. The interference effect evaluation unit is used to perform saturation region detection, spot geometric feature extraction, and thermal damage prediction based on the generated photoelectric imaging results, and outputs a quantitative evaluation result of the laser interference effect. The visualization terminal is used to visualize the quantitative evaluation results of spot energy distribution, detector status and laser interference effect in real time and generate standardized simulation evaluation reports. The scene initialization unit, laser parameter configuration unit, and target characteristic definition unit are connected to the input of the atmospheric transmission calculation unit. The output of the atmospheric transmission calculation unit is connected to the input of the load imaging simulation unit. The output of the load imaging simulation unit is connected to the input of the interference effect evaluation unit. The output of the interference effect evaluation unit is connected to the input of the visualization terminal.

[0008] Secondly, the present invention provides a laser interference full-link digital parallel simulation method based on a laser interference full-link digital parallel simulation system, comprising: The scene initialization unit imports geographic information data, binds material reflection properties, and sets the meteorological conditions corresponding to the simulation. The laser parameter configuration unit configures the laser interference mode, laser wavelength and power parameters, and laser emission geometric parameters. The target characteristic definition unit specifies the motion trajectory of the simulation target, sets the thermal radiation parameters of the simulation target, and loads the corresponding countermeasures for the simulation target. The atmospheric transmission calculation unit calculates laser energy attenuation, generates a turbulent phase screen, and iterates angular spectrum propagation to obtain the laser spot energy distribution after atmospheric transmission. The payload imaging simulation unit, based on the energy distribution of the light spot after atmospheric transmission, completes the modeling of the modulation transfer function of the optical system, the simulation of the photoelectric response of the detector, and the injection of imaging noise, generating photoelectric imaging results under laser interference. Based on the generated photoelectric imaging results, the interference effect evaluation unit completes saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect. The visualization terminal provides real-time visualization of the quantitative evaluation results of the spot energy distribution, detector status, and laser interference effect, and generates standardized simulation evaluation reports.

[0009] In some optional embodiments, the geographic information data includes DEM elevation maps and surface classification maps; the bidirectional reflectance distribution function parameters are Schlick model parameters called from a pre-built BRDF library; and meteorological conditions include atmospheric visibility, turbulence intensity, and temperature and humidity profile parameters.

[0010] In some optional embodiments, the laser parameter configuration unit integrates a laser hardware control interface and a parameter storage unit; the laser interference modes include pulse mode, continuous mode and composite modulation mode; the geometric parameters of laser emission include elevation angle and azimuth angle.

[0011] In some optional embodiments, the target characteristic definition unit supports CSV format motion path import and manual trajectory drawing; thermal radiation parameters include emissivity in the range of 0.1-1.0 and target temperature in the range of 200K-1500K; countermeasures include smoke screen release timing parameters and decoy launch strategy parameters.

[0012] In some optional embodiments, laser energy attenuation is calculated based on input parameters using the path integral method and based on the Beer-Lambert law, including: set up For transmission distance, For transmission distance Laser power density at that location This represents the initial power density of the laser emitter. The geometric distance for laser transmission. Let be the integral variable along the transmission path. Indicates the path position measured from the transmitting end. Atmospheric attenuation coefficient, It is the path location. The function of molecular absorption coefficient is The molecular scattering coefficient is The aerosol absorption coefficient is The aerosol scattering coefficient is Atmospheric attenuation coefficient is the molecular absorption coefficient. Molecular scattering coefficient Aerosol absorption coefficient Aerosol scattering coefficient The sum of them is: ; Power density after laser transmission Represented as: .

[0013] In some optional embodiments, based on Kolmogorov's turbulence theory, a power spectrum inversion method is used to generate a turbulent phase screen, including: set up Here, is the phase structure function, which describes the mean square value of the phase difference between two points in space. The distance between two points in space. Let be the atmospheric coherence length, then: ; set up The wavelength of the laser. For wave number, , The zenith angle is the angle between the laser transmission direction and the zenith direction. For transmission distance, The atmospheric refractive index structure constant is located at the path position. , The wind speed is in the vertical direction. The surface refractive index structure constant represents the intensity of near-surface atmospheric turbulence. It is the path location. The function is calculated using the Hufnagel-Valley model:

[0014] .

[0015] In some optional embodiments, iterative angular spectrum propagation is employed, using an angular spectrum propagation method to solve for beam diffusion with an accuracy proportional to the laser wavelength, including: set up for Spatial frequency of direction, for Spatial frequency of direction, Let be the angular spectrum transfer function. For transmission distance, The imaginary unit, The wavelength of the laser. For wave number, Let the step size be the propagation step size, then: ; The transmission distance is obtained through iterative calculation using Fourier transform. The complex amplitude distribution of the light spot at that location.

[0016] In some optional embodiments, the payload imaging simulation unit completes the modulation transfer function modeling of the optical system based on the spot energy distribution after atmospheric transmission, including: Construct the chain product of the modulation transfer function of the optical system; let... Let be the overall modulation transfer function of the optical system. For spatial frequency, Let be the diffraction-limited modulation transfer function, then: ; set up The cutoff frequency, , The entrance pupil diameter of the optical system. For the focal length of the optical system, For aberration modulation transfer function, The detector modulation transfer function. The modulation transfer function for a circuit is expressed as follows: For a circular aperture diffraction limit, the modulation transfer function is: ; The detector photoelectric response simulation includes: solving the charge accumulation equation based on the CCD / CMOS photoelectric conversion principle; and setting... The amount of signal charge accumulated in a single pixel. η Quantum efficiency represents the proportion of incident photons that are converted into electrons. Let be the electron charge constant, and take the value of . , Let be Planck's constant, and take the value of . , Let be the laser frequency, and the relationship between laser frequency and wavelength is: , At the speed of light, The laser power density incident on the pixel, The area of ​​a single pixel. For the integration time, the signal charge is: ; Imaging noise injection includes adding Gaussian noise and Poisson noise according to the signal-to-noise ratio requirements; set up The standard deviation of Gaussian noise. For signal strength, Given the signal-to-noise ratio (SNR), the relationship between the standard deviation of Gaussian noise and the SNR is as follows: ; set up Let Variance be the variance of Poisson noise. For signal strength, the variance of Poisson noise is equal to the signal strength, expressed as: .

[0017] In some optional embodiments, the interference effect evaluation unit, based on the generated photoelectric imaging results, performs saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect, including: Saturated regions are detected and saturated areas are calculated based on a grayscale threshold segmentation algorithm. The geometric features of the light spot are extracted, and the turbulence effect is considered to correct the far-field light spot radius; the geometric features of the light spot include the light spot radius, the light spot area, and the light spot shape factor; The laser power density distribution is calculated based on the heat conduction equation, and the thermal damage prediction results of the detector are obtained based on the laser power density distribution and the critical damage power density. The thermal damage prediction results include thermal damage probability and thermal damage level. The quantitative evaluation results of the output laser interference effect include at least one of the following: saturation area, spot radius, spot area, spot shape factor, thermal damage probability, and damage level.

[0018] The beneficial effects of this invention are as follows: This invention is the first to construct a full-link digital parallel simulation architecture for laser interference, encompassing ground scenarios, atmospheric transmission, payload imaging, and image quality assessment. The tightly coupled components can realistically reflect the complete physical process of laser interference, overcoming the shortcomings of existing technologies where each component is independent and unable to achieve full-process closed-loop simulation. It also addresses the problems of lack of full-link integration, insufficient dynamic process simulation, weak semi-physical verification capabilities, and poor scalability in existing technologies. This invention can be widely applied to the effectiveness evaluation of optoelectronic countermeasures equipment, optimization of countermeasure strategies, and verification of equipment performance. Furthermore, by calculating laser energy attenuation, generating turbulent phase screens, and iterative angular spectrum propagation, this invention achieves dynamic simulation of beam drift, expansion, and intensity flicker caused by atmospheric turbulence, enabling the reproduction of the dynamic countermeasures process of laser interference, significantly improving simulation accuracy and realism. Attached Figure Description Figure 1This is a system schematic diagram of a laser interference end-link digital parallel simulation system provided in Embodiment 1 of the present invention; Figure 2 This is a schematic diagram of a laser interference full-link digital parallel simulation method provided in Embodiment 2 of the present invention. Detailed Implementation

[0019] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations.

[0020] Example 1 As an example, to address the problems existing in the prior art, this embodiment provides a laser interference end-to-end digital parallel simulation system.

[0021] The following is a detailed description of the implementation details of a laser interference full-link digital parallel simulation system according to this embodiment. The following content is only for the convenience of understanding and is not necessary for implementing this solution.

[0022] This embodiment of a laser interference end-to-end digital parallel simulation system can be applied to electronic devices with communication, computing, and data storage capabilities. (See attached...) Figure 1 As shown in the figure, this embodiment provides a laser interference full-link digital parallel simulation system, including a scene initialization unit, a laser parameter configuration unit, a target characteristic definition unit, an atmospheric transmission calculation unit, a load imaging simulation unit, an interference effect evaluation unit, and a visualization terminal. The scene initialization unit is used to import geographic information data, bind material reflection properties, and set the corresponding meteorological conditions for the simulation. The laser parameter configuration unit is used to select the laser interference mode, configure the laser wavelength and power parameters, and define the geometric parameters of laser emission. The target characteristic definition unit is used to specify the motion trajectory of the simulated target, set the thermal radiation parameters of the simulated target, and load the corresponding countermeasures for the simulated target. The atmospheric transmission calculation unit is used to calculate laser energy attenuation, generate turbulent phase screen and iterative angular spectrum propagation based on input parameters, and calculate the spot energy distribution of laser after atmospheric transmission. The payload imaging simulation unit is used to complete the modeling of the modulation transfer function of the optical system, the simulation of the photoelectric response of the detector, and the injection of imaging noise based on the energy distribution of the light spot after atmospheric transmission, and to generate photoelectric imaging results under laser interference. The interference effect evaluation unit is used to perform saturation region detection, spot geometric feature extraction, and thermal damage prediction based on the generated photoelectric imaging results, and outputs a quantitative evaluation result of the laser interference effect. The visualization terminal is used to visualize the quantitative evaluation results of spot energy distribution, detector status and laser interference effect in real time and generate standardized simulation evaluation reports. The scene initialization unit, laser parameter configuration unit, and target characteristic definition unit are connected to the input of the atmospheric transmission calculation unit. The output of the atmospheric transmission calculation unit is connected to the input of the load imaging simulation unit. The output of the load imaging simulation unit is connected to the input of the interference effect evaluation unit. The output of the interference effect evaluation unit is connected to the input of the visualization terminal.

[0023] In some optional embodiments, the geographic information data includes DEM elevation maps and surface classification maps; the bidirectional reflectance distribution function parameters are Schlick model parameters called from a pre-built BRDF library; and meteorological conditions include atmospheric visibility, turbulence intensity, and temperature and humidity profile parameters.

[0024] In some optional embodiments, the laser parameter configuration unit integrates a laser hardware control interface and a parameter storage unit; the laser interference modes include pulse mode, continuous mode and composite modulation mode; the geometric parameters of laser emission include elevation angle and azimuth angle.

[0025] In some optional embodiments, the target characteristic definition unit supports CSV format motion path import and manual trajectory drawing; thermal radiation parameters include emissivity in the range of 0.1-1.0 and target temperature in the range of 200K-1500K; countermeasures include smoke screen release timing parameters and decoy launch strategy parameters.

[0026] In some optional embodiments, laser energy attenuation is calculated based on input parameters using the path integral method and based on the Beer-Lambert law, including: set up For transmission distance, For transmission distance Laser power density at that location This represents the initial power density of the laser emitter. The geometric distance for laser transmission. Let be the integral variable along the transmission path. Indicates the path position measured from the transmitting end. Atmospheric attenuation coefficient, It is the path location. The function of molecular absorption coefficient is The molecular scattering coefficient is The aerosol absorption coefficient is The aerosol scattering coefficient is Atmospheric attenuation coefficient is the molecular absorption coefficient. Molecular scattering coefficient Aerosol absorption coefficient Aerosol scattering coefficient The sum of them is: ; Power density after laser transmission Represented as: .

[0027] In some optional embodiments, based on Kolmogorov's turbulence theory, a power spectrum inversion method is used to generate a turbulent phase screen, including: set up Here, is the phase structure function, which describes the mean square value of the phase difference between two points in space. The distance between two points in space. Let be the atmospheric coherence length, then: ; set up The wavelength of the laser. For wave number, , The zenith angle is the angle between the laser transmission direction and the zenith direction. For transmission distance, The atmospheric refractive index structure constant is located at the path position. , The wind speed is in the vertical direction. The surface refractive index structure constant represents the intensity of near-surface atmospheric turbulence. It is the path location. The function is calculated using the Hufnagel-Valley model:

[0028] .

[0029] In some optional embodiments, iterative angular spectrum propagation is employed, using an angular spectrum propagation method to solve for beam diffusion with an accuracy proportional to the laser wavelength, including: set up for Spatial frequency of direction, for Spatial frequency of direction, Let be the angular spectrum transfer function. For transmission distance, The imaginary unit, The wavelength of the laser. For wave number, Let the step size be the propagation step size, then: ; In some optional embodiments, the payload imaging simulation unit completes the modulation transfer function modeling of the optical system based on the spot energy distribution after atmospheric transmission, including: Construct the chain product of the modulation transfer function of the optical system; let... Let be the overall modulation transfer function of the optical system. For spatial frequency, Let be the diffraction-limited modulation transfer function, then: ; set up The cutoff frequency, , The entrance pupil diameter of the optical system. For the focal length of the optical system, For aberration modulation transfer function, The detector modulation transfer function. The modulation transfer function for a circuit is expressed as follows: For a circular aperture diffraction limit, the modulation transfer function is: ; The detector photoelectric response simulation includes: solving the charge accumulation equation based on the CCD / CMOS photoelectric conversion principle; and setting... The amount of signal charge accumulated in a single pixel. η Quantum efficiency represents the proportion of incident photons that are converted into electrons. Let be the electron charge constant, and take the value of . , Let be Planck's constant, and take the value of . , Let be the laser frequency, and the relationship between laser frequency and wavelength is: , At the speed of light, The laser power density incident on the pixel, The area of ​​a single pixel. For the integration time, the signal charge is: ; Imaging noise injection includes adding Gaussian noise and Poisson noise according to the signal-to-noise ratio requirements; set up The standard deviation of Gaussian noise. For signal strength, Given the signal-to-noise ratio (SNR), the relationship between the standard deviation of Gaussian noise and the SNR is as follows: ; set up Let Variance be the variance of Poisson noise. For signal strength, the variance of Poisson noise is equal to the signal strength, expressed as: .

[0030] In some optional embodiments, the interference effect evaluation unit, based on the generated photoelectric imaging results, performs saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect, including: Saturated regions are detected and saturated areas are calculated based on a grayscale threshold segmentation algorithm. The geometric features of the light spot are extracted, and the turbulence effect is considered to correct the far-field light spot radius; the geometric features of the light spot include the light spot radius, the light spot area, and the light spot shape factor; The laser power density distribution is calculated based on the heat conduction equation, and the thermal damage prediction results of the detector are obtained based on the laser power density distribution and the critical damage power density. The thermal damage prediction results include thermal damage probability and thermal damage level. The quantitative evaluation results of the output laser interference effect include at least one of the following: saturation area, spot radius, spot area, spot shape factor, thermal damage probability, and damage level.

[0031] Example 2 Based on the same principles as the system shown in Embodiment 1 of the present invention, as illustrated in the appendix... Figure 2 As shown, the embodiments of the present invention also provide a laser interference full-link digital parallel simulation method based on a laser interference full-link digital parallel simulation system, including the following steps 110-170.

[0032] Step 110: The scene initialization unit imports geographic information data, binds material reflection attributes, and sets the meteorological conditions corresponding to the simulation. The laser parameter configuration unit configures the laser interference mode, laser wavelength and power parameters, and laser emission geometric parameters. The target characteristic definition unit specifies the motion trajectory of the simulation target, sets the thermal radiation parameters of the simulation target, and loads the corresponding countermeasures for the simulation target.

[0033] In some optional embodiments, the geographic information data includes DEM elevation maps and surface classification maps; the bidirectional reflectance distribution function parameters are Schlick model parameters called from a pre-built BRDF library; and meteorological conditions include atmospheric visibility, turbulence intensity, and temperature and humidity profile parameters.

[0034] In some optional embodiments, the laser parameter configuration unit integrates a laser hardware control interface and a parameter storage unit; the laser interference modes include pulse mode, continuous mode and composite modulation mode; the geometric parameters of laser emission include elevation angle and azimuth angle, for example, the elevation angle range is (-30°~+90°) and the azimuth angle range is 0°~360°.

[0035] Hardware interfaces include the SCPI instruction set, RS422 protocol, and Camera Link interface. These standardized hardware interfaces support hardware-in-the-loop simulation with physical lasers and detector arrays, enabling closed-loop verification between simulation and actual measurements. In practical applications, a loosely coupled modular design supports flexible configuration for different laser wavelengths, detector types, and weather conditions, exhibiting excellent scalability and reusability.

[0036] In some optional embodiments, the target characteristic definition unit supports CSV format motion path import and manual trajectory drawing; thermal radiation parameters include emissivity in the range of 0.1-1.0 and target temperature in the range of 200K-1500K; countermeasures include smoke screen release timing parameters and decoy launch strategy parameters.

[0037] In some optional embodiments, a laser interference end-to-end digital parallel simulation method is performed by a simulation server cluster.

[0038] Step 120: The atmospheric transmission solution unit calculates the laser energy attenuation, generates the turbulent phase screen, and iterates the angular spectrum propagation to obtain the laser spot energy distribution after atmospheric transmission.

[0039] In some optional embodiments, laser energy attenuation is calculated based on input parameters using the path integral method and based on the Beer-Lambert law, including: set up For transmission distance, For transmission distance Laser power density at that location This represents the initial power density of the laser emitter. The geometric distance for laser transmission. Let be the integral variable along the transmission path. Indicates the path position measured from the transmitting end. Atmospheric attenuation coefficient, It is the path location. The function of molecular absorption coefficient is The molecular scattering coefficient is The aerosol absorption coefficient is The aerosol scattering coefficient is Atmospheric attenuation coefficient is the molecular absorption coefficient. Molecular scattering coefficient Aerosol absorption coefficient Aerosol scattering coefficient The sum of them is: ; Power density after laser transmission Represented as: .

[0040] In some optional embodiments, based on Kolmogorov's turbulence theory, a power spectrum inversion method is used to generate a turbulent phase screen, including: set up Here, is the phase structure function, which describes the mean square value of the phase difference between two points in space. The distance between two points in space. Let be the atmospheric coherence length, then: ; set up The wavelength of the laser. For wave number, , The zenith angle is the angle between the laser transmission direction and the zenith direction. For transmission distance, The atmospheric refractive index structure constant is located at the path position. , The wind speed is in the vertical direction. The surface refractive index structure constant represents the intensity of near-surface atmospheric turbulence. It is the path location. The function is calculated using the Hufnagel-Valley model:

[0041] .

[0042] By simulating the phase modulation of the beam wavefront by atmospheric turbulence, phase perturbation is provided for subsequent angular spectrum propagation.

[0043] In some optional embodiments, iterative angular spectrum propagation is employed, using an angular spectrum propagation method to solve for beam diffusion with an accuracy proportional to the laser wavelength, including: set up for Spatial frequency of direction, for Spatial frequency of direction, Let be the angular spectrum transfer function. For transmission distance, The imaginary unit, The wavelength of the laser. For wave number, Let the step size be the propagation step size, then: ; The transmission distance is obtained through iterative calculation using Fourier transform. The complex amplitude distribution of the light spot at that location.

[0044] Based on Kolmogorov's turbulence theory and the angular spectrum propagation method, dynamic simulations of beam drift, spread, and intensity scintillation caused by atmospheric turbulence were achieved. In practical applications, the angular spectrum propagation iterations were used to solve for beam spread with an accuracy of λ / 10.

[0045] The far-field spot energy distribution of the laser beam after atmospheric attenuation and turbulence modulation is accurately calculated by iterative angular spectrum propagation.

[0046] Step 130: Based on the energy distribution of the light spot after atmospheric transmission, the payload imaging simulation unit completes the modeling of the modulation transfer function of the optical system, the simulation of the detector photoelectric response, and the injection of imaging noise, generating photoelectric imaging results under laser interference.

[0047] In some optional embodiments, the payload imaging simulation unit completes the modulation transfer function modeling of the optical system based on the spot energy distribution after atmospheric transmission, including: Construct the chain product of the modulation transfer function of the optical system; let... Let be the overall modulation transfer function of the optical system. For spatial frequency, Let be the diffraction-limited modulation transfer function, then: ; set up The cutoff frequency, , The entrance pupil diameter of the optical system. For the focal length of the optical system, For aberration modulation transfer function, The detector modulation transfer function. The modulation transfer function for a circuit is expressed as follows: For a circular aperture diffraction limit, the modulation transfer function is: ; By constructing a chain product of the modulation transfer function of an optical system, we can comprehensively characterize the transfer characteristics of optical systems, detectors, and circuits to spatial frequencies.

[0048] The detector photoelectric response simulation includes: solving the charge accumulation equation based on the CCD / CMOS photoelectric conversion principle; and setting... The amount of signal charge accumulated in a single pixel. ηQuantum efficiency represents the proportion of incident photons that are converted into electrons. Let be the electron charge constant, and take the value of . , Let be Planck's constant, and take the value of . , Let be the laser frequency, and the relationship between laser frequency and wavelength is: , At the speed of light, The laser power density incident on the pixel, The area of ​​a single pixel. For the integration time, the signal charge is: ; When the laser power density exceeds the saturation threshold, charge overflow occurs. Let... Crosstalk rate is the percentage of overflow charge relative to the total charge. This represents the amount of overflow charge. For the total amount of charge generated, The potential difference between the potential wells depends on the device structure. Here is the Boltzmann constant, which takes the value of , Given the detector's operating temperature, the crosstalk rate is expressed as: .

[0049] By simulating the photoelectric response of the detector, it is possible to simulate the generation, accumulation, and overflow of charge in the detector under laser irradiation.

[0050] Imaging noise injection includes adding Gaussian noise and Poisson noise according to the signal-to-noise ratio requirements; set up The standard deviation of Gaussian noise. For signal strength, Given the signal-to-noise ratio (SNR), the relationship between the standard deviation of Gaussian noise and the SNR is as follows: ; set up Let Variance be the variance of Poisson noise. For signal strength, the variance of Poisson noise is equal to the signal strength, expressed as: .

[0051] In practical applications, the atmospheric transport calculation module is a parallel computing cluster composed of multiple GPU computing nodes connected via a high-speed network. Employing this GPU parallel computing cluster architecture, it achieves high-frame-rate real-time rendering, supports interactive parameter adjustment and real-time effect feedback, and can quickly complete laser interference performance simulation and comparative analysis under various operating conditions, significantly improving simulation efficiency.

[0052] By injecting imaging noise into the actual noise characteristics of the simulated imaging system, the realism of the simulated image is improved.

[0053] Step 140: Based on the generated photoelectric imaging results, the interference effect evaluation unit completes saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect.

[0054] In some optional embodiments, the interference effect evaluation unit, based on the generated photoelectric imaging results, performs saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect, including: Saturated regions are detected and saturated areas are calculated based on a grayscale threshold segmentation algorithm. The geometric features of the light spot are extracted, and the turbulence effect is considered to correct the far-field light spot radius; the geometric features of the light spot include the light spot radius, the light spot area, and the light spot shape factor; The laser power density distribution is calculated based on the heat conduction equation, and the thermal damage prediction results of the detector are obtained based on the laser power density distribution and the critical damage power density. The thermal damage prediction results include thermal damage probability and thermal damage level. The quantitative evaluation results of the output laser interference effect include at least one of the following: saturation area, spot radius, spot area, spot shape factor, thermal damage probability, and damage level.

[0055] Specifically, saturated regions are detected based on a grayscale threshold segmentation algorithm, with the threshold set to 95% of the maximum grayscale value.

[0056] Calculate the geometric parameters of the saturated spot, including its radius, area, and shape factor. Let... Let be the radius of the saturated spot. The total laser power, Let be the detector saturation threshold, i.e., the incident power density at which a pixel just reaches saturation. Then, the relationship between the saturation spot radius and the laser irradiance is expressed as: set up To account for the effective spot radius after considering turbulence effects, The radius of the light spot when there is no turbulence. For beam drift variance, To account for the beam spread variance and considering turbulence effects, the far-field spot radius is corrected as follows: .

[0057] Specifically, the risk of device meltdown is assessed based on the heat conduction equation. Let... For material density, For specific heat capacity, For position ,time The temperature at that location Thermal conductivity, For the Laplace operator, The material absorption coefficient, Given the laser power density distribution, the heat conduction equation can be expressed as: ; set up This is the critical damage power density; exceeding this value may cause permanent damage. For material density, For specific heat capacity, The melting point of the material. The initial temperature, The material absorption coefficient, Let be the laser pulse width. For continuous lasers, taking the equivalent thermal action time, the formula for calculating the critical damage power density is: .

[0058] By assessing the risk of device meltdown using the heat conduction equation, it is possible to predict whether laser irradiation will cause permanent damage to the detector.

[0059] Based on the charge diffusion equation and the heat conduction equation, the time-varying characteristics of detector saturation crosstalk and heat accumulation were simulated.

[0060] Step 150: The visualization terminal displays the quantitative evaluation results of the spot energy distribution, detector status and laser interference effect in real time and generates a standardized simulation evaluation report.

[0061] Specifically, the control visualization terminal generates and displays a comprehensive simulation report. The visualization content includes: laser spot morphology display, pseudo-color mapping energy distribution (0-65535 gray levels); detector status including saturation and temperature status, with detector status display: red flashing warning in oversaturated areas; parameter curve plotting: real-time display of changes in indicators such as MTF (Modulation Transfer Function), SNR (Signal-to-Noise Ratio), and spot radius; and quantitative report output: generating a PDF format evaluation report.

[0062] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A laser interference end-to-end digital parallel simulation system, characterized in that, It includes a scene initialization unit, a laser parameter configuration unit, a target characteristic definition unit, an atmospheric transmission calculation unit, a payload imaging simulation unit, an interference effect evaluation unit, and a visualization terminal; The scene initialization unit is used to import geographic information data, bind material reflection properties, and set the corresponding meteorological conditions for the simulation. The laser parameter configuration unit is used to select the laser interference mode, configure the laser wavelength and power parameters, and define the geometric parameters of laser emission. The target characteristic definition unit is used to specify the motion trajectory of the simulated target, set the thermal radiation parameters of the simulated target, and load the corresponding countermeasures for the simulated target. The atmospheric transmission calculation unit is used to calculate laser energy attenuation, generate turbulent phase screen and iterative angular spectrum propagation based on input parameters, and calculate the spot energy distribution of laser after atmospheric transmission. The payload imaging simulation unit is used to complete the modeling of the modulation transfer function of the optical system, the simulation of the photoelectric response of the detector, and the injection of imaging noise based on the energy distribution of the light spot after atmospheric transmission, and to generate photoelectric imaging results under laser interference. The interference effect evaluation unit is used to perform saturation region detection, spot geometric feature extraction, and thermal damage prediction based on the generated photoelectric imaging results, and outputs a quantitative evaluation result of the laser interference effect. The visualization terminal is used to visualize the quantitative evaluation results of spot energy distribution, detector status and laser interference effect in real time and generate standardized simulation evaluation reports. The scene initialization unit, laser parameter configuration unit, and target characteristic definition unit are connected to the input of the atmospheric transmission calculation unit. The output of the atmospheric transmission calculation unit is connected to the input of the load imaging simulation unit. The output of the load imaging simulation unit is connected to the input of the interference effect evaluation unit. The output of the interference effect evaluation unit is connected to the input of the visualization terminal.

2. The laser interference end-to-end digital parallel simulation system according to claim 1, characterized in that, The geographic information data includes DEM elevation maps and surface classification maps; the bidirectional reflectance distribution function parameters are Schlick model parameters called from the pre-built BRDF library; meteorological conditions include atmospheric visibility, turbulence intensity, and temperature and humidity profile parameters.

3. The laser interference end-to-end digital parallel simulation system according to claim 1, characterized in that, The laser parameter configuration unit integrates a laser hardware control interface and a parameter storage unit; the laser interference modes include pulse mode, continuous mode and composite modulation mode; the geometric parameters of laser emission include elevation angle and azimuth angle.

4. The laser interference end-to-end digital parallel simulation system according to claim 1, characterized in that, The target characteristic definition unit supports importing motion paths in CSV format and manual trajectory drawing; Thermal radiation parameters include emissivity in the range of 0.1-1.0 and target temperature in the range of 200K-1500K; countermeasures include smoke screen release timing parameters and decoy launch strategy parameters.

5. A laser interference full-link digital parallel simulation method based on a laser interference full-link digital parallel simulation system according to any one of claims 1-4, characterized in that, include: The scene initialization unit imports geographic information data, binds material reflection properties, and sets the meteorological conditions corresponding to the simulation. The laser parameter configuration unit configures the laser interference mode, laser wavelength and power parameters, and laser emission geometric parameters. The target characteristic definition unit specifies the motion trajectory of the simulation target, sets the thermal radiation parameters of the simulation target, and loads the corresponding countermeasures for the simulation target. The atmospheric transmission calculation unit calculates laser energy attenuation, generates a turbulent phase screen, and iterates angular spectrum propagation to obtain the laser spot energy distribution after atmospheric transmission. The payload imaging simulation unit, based on the energy distribution of the light spot after atmospheric transmission, completes the modeling of the modulation transfer function of the optical system, the simulation of the photoelectric response of the detector, and the injection of imaging noise, generating photoelectric imaging results under laser interference. Based on the generated photoelectric imaging results, the interference effect evaluation unit completes saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect. The visualization terminal provides real-time visualization of the quantitative evaluation results of the spot energy distribution, detector status, and laser interference effect, and generates standardized simulation evaluation reports.

6. The laser interference end-to-end digital parallel simulation method according to claim 5, characterized in that, Laser energy attenuation is calculated based on input parameters using the path integral method and the Beer-Lambert law, including: set up For transmission distance, For transmission distance Laser power density at that location The initial power density of the laser emitter. The geometric distance for laser transmission. Let be the integral variable along the transmission path. Indicates the path position measured from the transmitting end. Atmospheric attenuation coefficient, It is the path location. The function of molecular absorption coefficient is The molecular scattering coefficient is The aerosol absorption coefficient is The aerosol scattering coefficient is Atmospheric attenuation coefficient is the molecular absorption coefficient. Molecular scattering coefficient Aerosol absorption coefficient Aerosol scattering coefficient The sum of them is: ； Power density after laser transmission Represented as: 。 7. The laser interference end-to-end digital parallel simulation method according to claim 5, characterized in that, Based on Kolmogorov's turbulence theory, a power spectrum inversion method is used to generate a turbulent phase screen, including: set up Here, is the phase structure function, which describes the mean square value of the phase difference between two points in space. The distance between two points in space. Let be the atmospheric coherence length, then: ； set up The wavelength of the laser. For wave number, , The zenith angle is the angle between the laser transmission direction and the zenith direction. For transmission distance, The atmospheric refractive index structure constant is located at the path position. , The wind speed is in the vertical direction. The surface refractive index structure constant represents the intensity of near-surface atmospheric turbulence. It is the path location. The function is calculated using the Hufnagel-Valley model: 。 8. The laser interference end-to-end digital parallel simulation method according to claim 5, characterized in that, Iterative angular spectrum propagation, using the angular spectrum propagation method, solves for beam diffusion with a precision proportional to the laser wavelength, including: set up for Spatial frequency of direction, for Spatial frequency of direction, Let be the angular spectrum transfer function. For transmission distance, The imaginary unit, The wavelength of the laser. For wave number, If the step size is a single propagation step, then: ； The transmission distance is obtained through iterative calculation using Fourier transform. The complex amplitude distribution of the light spot at that location.

9. The laser interference end-to-end digital parallel simulation method according to claim 5, characterized in that, The payload imaging simulation unit, based on the light spot energy distribution after atmospheric transmission, completes the modulation transfer function modeling of the optical system, including: Construct the chain product of the modulation transfer function of the optical system; let... Let be the overall modulation transfer function of the optical system. For spatial frequency, Let be the diffraction-limited modulation transfer function, then: ； set up The cutoff frequency, , The entrance pupil diameter of the optical system. For the focal length of the optical system, For aberration modulation transfer function, The detector modulation transfer function. The modulation transfer function for a circuit is expressed as follows: For a circular aperture diffraction limit, the modulation transfer function is: ； The detector photoelectric response simulation includes: solving the charge accumulation equation based on the CCD / CMOS photoelectric conversion principle; and setting... The amount of signal charge accumulated in a single pixel. η Quantum efficiency represents the proportion of incident photons that are converted into electrons. Let be the electron charge constant, and take the value of . , Let be Planck's constant, and take the value of . , Let be the laser frequency, and the relationship between laser frequency and wavelength is: , At the speed of light, The laser power density incident on the pixel, The area of ​​a single pixel. For the integration time, the signal charge is: ； Imaging noise injection includes adding Gaussian noise and Poisson noise according to the signal-to-noise ratio requirements; set up The standard deviation of Gaussian noise. For signal strength, Given the signal-to-noise ratio (SNR), the relationship between the standard deviation of Gaussian noise and the SNR is as follows: ； set up Let Variance be the variance of Poisson noise. Given signal strength, the variance of Poisson noise is equal to the signal strength, expressed as: .

10. The laser interference end-to-end digital parallel simulation method according to claim 5, characterized in that, Based on the generated photoelectric imaging results, the interference effect evaluation unit completes saturation region detection, spot geometric feature extraction, and thermal damage prediction, and outputs a quantitative evaluation result of the laser interference effect, including: Saturated regions are detected and saturated areas are calculated based on a grayscale threshold segmentation algorithm. The geometric features of the light spot are extracted, and the turbulence effect is considered to correct the far-field light spot radius; the geometric features of the light spot include the light spot radius, the light spot area, and the light spot shape factor; The laser power density distribution is calculated based on the heat conduction equation, and the thermal damage prediction results of the detector are obtained based on the laser power density distribution and the critical damage power density. The thermal damage prediction results include thermal damage probability and thermal damage level. The quantitative evaluation results of the output laser interference effect include at least one of the following: saturation area, spot radius, spot area, spot shape factor, thermal damage probability, and damage level.

Images