A method and system for generating interference light field distribution characteristics of an arbitrary reflective surface velocity diagnostic system

Abstract

This invention relates to a method and system for generating the interference light field distribution characteristics of an arbitrary reflective surface velocity diagnostic system, belonging to the field of inertial confinement fusion diagnostic technology. Addressing the problem that existing simulation methods rely on a single input velocity field and cannot simulate complex dynamic processes, this invention constructs a customizable velocity field input system, supporting the invocation of typical velocity field functions or the import of external three-dimensional velocity field data; it simulates the temporal blurring effect of the slit in a fringe camera through convolution operations; and it rigorously models the entire physical process from Doppler frequency shift to interference fringe generation based on optical and wave theories. This invention effectively improves the diversity and realism of the simulation model, providing reliable data support for the verification of shock wave velocity field reconstruction algorithms.

Description

Technical Field

[0001] This invention belongs to the field of inertial confinement fusion diagnostic technology, and relates to a method and system for generating the interference optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system. Background Technology

[0002] With the global energy structure transformation and the continuous increase in demand for clean energy, nuclear energy, as an important energy form with high energy density and low carbon emissions, has received widespread research and attention, with controlled nuclear fusion being a key research focus. Inertial Confinement Fusion (ICF) is one method of achieving controlled nuclear fusion. Its principle involves injecting a large amount of energy into a target containing fusion fuel within an extremely short time. Under extremely high temperature, pressure, and density, the fuel within the target undergoes a compression implosion, leading to a thermonuclear fusion reaction. During the fusion process, factors such as cavity structure and uneven laser intensity can easily cause asymmetry in the compression implosion, resulting in fusion failure. To overcome the asymmetry in the implosion compression process, it is necessary to measure the velocity distribution of the shock wave to predict the compression state of the target, providing reliable reference data for optimizing the ICF compression implosion process.

[0003] The Velocity Interferometer System for AnyReflector (VISAR) is a high-precision dynamic velocity measurement and diagnostic technique based on the laser Doppler effect, and is a key technology for measuring shock wave velocity during ICF (Inductively Coupled Fluid Dynamics). This technique measures the frequency change of reflected light caused by the high-speed motion of the target surface and uses the principle of interference to convert the velocity information into phase changes of interference fringes, thereby achieving precise acquisition of the target's velocity over time.

[0004] Traditional VISAR systems use streak cameras for imaging and recording, which can only record the velocity field distribution along a single line over time. Because the image recorded when the streak camera slit is fully open suffers from aliasing during scanning, the slit opening must be reduced to obtain one-dimensional spatial information. While charge-coupled device (CCD) cameras or framing cameras can record two-dimensional spatial velocity field information in VISAR systems, their temporal resolution is limited. The lack of imaging and recording equipment with both high temporal and spatial resolution has hindered the research progress of ICF (Integrated Circuits for Imaging and Recording).

[0005] Compressed sensing technology can reconstruct aliased images recorded by a streak camera with the slit fully open. Combining streak cameras with compressed sensing technology to form compressed ultrafast photography (CUP) has ultra-high spatiotemporal resolution and can achieve two-dimensional spatial imaging. The CUP-VISAR system, which combines CUP technology with VISAR, can achieve two-dimensional velocity field recording at high temporal resolution, thus promoting further research in ICF.

[0006] Despite the significant advantages of CUP-VISAR technology in measuring two-dimensional shock wave velocity fields, problems and technical challenges remain in practical ICF research and engineering applications. Firstly, the optical system is complex, with numerous sources of system error. A VISAR system comprises core components such as a beam splitter, delay arm, interferometer cavity, imaging lens group, and fringe camera. After beam splitting, the laser beams enter different optical paths and are re-superimposed after a certain optical path difference to form interference fringes. Positional offsets, tilt angle errors, or optical path mismatches of optical elements in the system can all lead to interference phase shifts, thus affecting fringe stability and visibility. Vibrations, electromagnetic interference, and optical alignment errors generated in the shock experimental environment can all introduce additional phase noise, reducing fringe contrast and affecting velocity accuracy. In the CUP structure, compression coding and other methods further amplify the noise propagation effect, making the reconstruction results highly sensitive to system stability. Because the actual experimental operation of CUP-VISAR is prone to introducing errors and noise, the experiment is quite difficult; therefore, designing simulation methods is of great significance for research.

[0007] Secondly, the compressed sampling reconstruction algorithm itself has high computational complexity. CUP-VISAR typically requires iterative calculations on multiple frames of high-resolution stripe data. Especially after introducing prior regularization constraints, the solution process involves large-scale matrix operations and multiple convergence checks. When the data size reaches hundreds of rows and thousands of columns of pixels, the reconstruction time increases significantly. For time-consuming subsequent data processing, if the results are inaccurate due to experimental errors, it will seriously affect the research efficiency. Using simulation methods to customize the input field can eliminate error interference, obtain theoretical verification, and then corroborate the experimental results, which is of great value for further research.

[0008] Furthermore, existing research largely focuses on reconstruction algorithm optimization and experimental system improvement, while paying relatively little attention to the diversity and controllability of the input shock wave velocity field. In practical shock physics research, the velocity field may exhibit complex characteristics such as non-ideal symmetric distribution. Without a flexible and adjustable input field model, it is not conducive to system performance evaluation and algorithm applicability verification. Based on this, it is necessary to propose a custom shock wave input velocity field simulation method for CUP-VISAR systems to meet the simulation requirements of complex shock processes. Summary of the Invention

[0009] In view of this, the purpose of the present invention is to provide a method and system for generating the interference optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system.

[0010] To achieve the above objectives, the present invention provides the following technical solution: A method for generating the interferometric optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system includes the following steps: S1: Define the motion state distribution of the shock wave to be measured. The motion state distribution is obtained by calling a preset typical velocity field function or importing externally defined three-dimensional velocity field data. S2: Configure the parameters of laser wavelength, etalon refractive index, etalon initial thickness and stripe velocity constant, calculate the Doppler frequency shift according to the motion state distribution, and generate modulated signal light with velocity field information; S3: The modulated signal light is split into two paths, one of which is optically delayed and then superimposed with the other to form interference fringes, thus obtaining an ideal interference light intensity matrix; S4: Based on the physical response time of the slit of the fringe camera, perform a one-dimensional convolution operation on the ideal interference intensity matrix in the time axis direction to generate a simulated interference fringe image.

[0011] Furthermore, the typical velocity field function mentioned in S1 includes one or more of the following: static velocity field function, smooth velocity field function, step velocity field function, and sinusoidal velocity field function. The externally customized three-dimensional velocity field data is imported after preprocessing by numerical calculation software.

[0012] Furthermore, the fringe velocity constant described in S2 VPF The calculation formula is:

[0013] in, The wavelength of the laser. The optical delay time generated by the etalon.

[0014] Furthermore, S2 also includes: according to the set VPF and The optical delay time of the etalon is calculated using the formula. and according to Correct the physical thickness of the etalon.

[0015] Furthermore, the formula for calculating the intensity distribution of the interference fringes described in S3 is as follows:

[0016] in, , These represent the amplitudes of the two signal beams, , These are the angular frequencies of the two signal beams, respectively. , These represent the initial phases of the two signal beams. t It is a time variable.

[0017] Furthermore, S4 specifically includes: calculating the average simulation step size of the time axis array. dt Based on the slit physical response time of the streak camera Calculate the number of pixels occupied by the slit width tpx The tpx The calculation formula is:

[0018] in, This indicates the floor function.

[0019] Furthermore, S4 also includes: constructing a length of tpx All-1 window function w Along the time axis, the ideal interference intensity matrix Perform a one-dimensional convolution operation to obtain a simulated interference fringe image. The convolution operation formula is:

[0020] in, y Represents spatial location coordinates, t Represents the time coordinate. m This is the kernel index variable.

[0021] A system for generating the interferometric optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system includes: The velocity field definition module is used to define the motion state distribution of the shock wave to be measured. The motion state distribution is obtained by calling a preset typical velocity field function or importing externally defined three-dimensional velocity field data. The Doppler modulation signal generation module is used to configure parameters such as laser wavelength, etalon refractive index, etalon initial thickness and fringe velocity constant, calculate Doppler frequency shift according to the motion state distribution, and generate modulation signal light with velocity field information. An interference detection module is used to split the modulated signal light into two paths, so that one path generates an optical delay and is superimposed with the other path to form interference fringes, thereby obtaining an ideal interference light intensity matrix. The temporal blur processing module is used to perform a one-dimensional convolution operation on the ideal interference intensity matrix in the time axis direction based on the physical response time of the slit of the fringe camera, so as to generate a simulated interference fringe image.

[0022] Furthermore, the velocity field definition module includes a preset velocity field calling unit and an external data import unit. The preset velocity field calling unit is used to store static velocity field functions, smooth velocity field functions, step velocity field functions, and sinusoidal velocity field functions. The external data import unit is used to receive and parse externally customized three-dimensional velocity field data files.

[0023] Furthermore, the temporal blurring processing module includes a convolution kernel construction unit and a convolution operation unit. The convolution kernel construction unit is used to calculate the number of pixels occupied by the slit width based on the slit physical response time of the streak camera. tpx and construct a length of tpx All-1 window function w The convolution operation unit is used to perform a one-dimensional convolution operation on the ideal interference light intensity matrix along the time axis.

[0024] The beneficial effects of this invention are as follows: (1) A simulation system with customizable input velocity field was constructed, breaking the limitation of traditional simulation methods that can only generate interference fringes using fixed velocity field forms. Users can directly call preset typical velocity field functions or import externally customized three-dimensional velocity field data, realizing the simulation of various complex velocity distributions that may occur in actual experiments. This provides diverse test data for the verification of velocity field reconstruction algorithms and effectively improves the adaptability of the simulation model to different impact processes.

[0025] (2) A simulation mechanism for the diagnostic characteristics of a fringe camera based on convolution calculation was introduced to restore the temporal ambiguity features under real detection conditions. By performing convolution processing on the ideal interference light intensity along the time axis according to the physical response time of the fringe camera, the generated simulated interference fringes are highly consistent with the images recorded by the real camera in terms of visual effect and physical characteristics. This solves the problem of the disconnect between the idealized model and the real experimental environment in traditional simulation, and provides more reliable data support for the stability verification of subsequent reconstruction algorithms in complex backgrounds.

[0026] (3) From laser emission, Doppler frequency shift generated by reflection from the surface of the moving target, optical path difference time delay generated by the etalon, and superposition of two-way light interference, to the final slit integration of the striated camera, each step is strictly modeled mathematically based on the basic formulas of optics and wave theory, avoiding errors caused by simplified empirical models or approximations, and ensuring that the simulation results are highly consistent with the actual physical process in terms of striated morphology and spatiotemporal distribution.

[0027] (4) Through high-precision physical simulation, researchers can theoretically verify the system before the experiment, eliminate error interference, reduce repetitive work caused by experimental conditions, improve research efficiency, and provide strong technical support for the accurate diagnosis of shock wave velocity field in cutting-edge physics research such as inertial confinement fusion.

[0028] Other advantages, objectives, and features of the invention will be set forth in part in the description which follows, and in part will be apparent to those skilled in the art from the following examination, or may be learned from practice of the invention. The objectives and other advantages of the invention can be realized and obtained through the following description. Attached Figure Description

[0029] To make the objectives, technical solutions, and advantages of the present invention clearer, the preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings, wherein: Figure 1 Here is a flowchart of the simulation method; Figure 2 Interference optical fields were generated for sinusoidal perturbation velocity fields, and experimental test results were obtained. Among them, (a) is a three-dimensional sinusoidal velocity field image; (b) is the VISAR system output interference fringe image obtained by simulation; and (c) is the experimental test results. Figure 3 The simulation results are for a custom complex velocity field; (a) is a three-dimensional custom velocity field image; (b) is the VISAR system output interferometric fringe image obtained from the simulation; and (c) is the actual VISAR experimental image. Detailed Implementation

[0030] The following specific examples illustrate the implementation of the present invention. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of the present invention. Unless otherwise specified, the following embodiments and features can be combined with each other.

[0031] The accompanying drawings are for illustrative purposes only and are schematic diagrams, not actual pictures. They should not be construed as limiting the invention. To better illustrate the embodiments of the invention, some parts in the drawings may be omitted, enlarged, or reduced, and do not represent the actual product dimensions. It is understandable to those skilled in the art that some well-known structures and their descriptions may be omitted in the drawings.

[0032] In the accompanying drawings of the embodiments of the present invention, the same or similar reference numerals correspond to the same or similar components. In the description of the present invention, it should be understood that if terms such as "upper," "lower," "left," "right," "front," and "rear" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, they are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, the terms used to describe positional relationships in the drawings are only for illustrative purposes and should not be construed as limiting the present invention. For those skilled in the art, the specific meaning of the above terms can be understood according to the specific circumstances.

[0033] like Figure 1 As shown, Figure 1 This is a flowchart of the simulation method. The diagram visually illustrates the complete logical flow from defining the velocity field in S1, through Doppler modulation in S2, interferometric detection in S3, and time fuzzing processing in S4, ultimately outputting a simulated interference fringe image.

[0034] like Figure 2 As shown, Figure 2 Interference light fields were generated for sinusoidal perturbation velocity fields, and experimental test results were presented. (a) shows a three-dimensional sinusoidal velocity field image, with velocity values ​​fluctuating periodically in space. (b) shows the simulated VISAR system output interference fringe image, where the fringes exhibit a regular curvature corresponding to the velocity field. (c) shows the experimental test results; the comparison shows that the morphology of the simulated fringes highly matches the measured results, verifying the simulation capability of this method for typical input velocity fields.

[0035] like Figure 3 As shown, Figure 3 The simulation results are for a custom complex velocity field. (a) shows a 3D custom velocity field image with irregular spatiotemporal variation characteristics. (b) shows the simulated VISAR system output interferometric fringe image, exhibiting obvious fringe curvature and density variations in the corresponding region. (c) shows the actual VISAR experimental image. The simulation results are consistent with the real experimental data in fringe morphology and distribution characteristics, demonstrating the effectiveness and versatility of this method for complex input velocity fields.

[0036] This invention provides a method for generating the interferometric optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system. This method is implemented through software and runs on a high-performance computing device. The program pre-stores four typical velocity field function models: stationary, smooth, step, and sinusoidal step.

[0037] In practice, S1 is executed first to define the motion state distribution of the shock wave to be measured. The user can select the input mode through the human-computer interaction interface. If the standard mode is selected, the system calls a preset velocity field function, such as a sinusoidal velocity field function. ,in A For amplitude, f For frequency, users need to set specific amplitude and frequency parameters. If the custom mode is selected, the system receives an external three-dimensional velocity field data file imported by the user. This file contains the correspondence between spatiotemporal coordinates and velocity values. The system uses a trilinear interpolation algorithm to map the external data onto the current simulation mesh, ensuring the continuity of the velocity field in space and time.

[0038] Next, execute S2 to configure the relevant physical parameters and generate a Doppler modulation signal. The system sets the laser wavelength. The refractive index of the etalon n and initial thickness d To meet measurement accuracy requirements, a fringe velocity constant (VPF) needs to be set. VPF The velocity change corresponding to one cycle of stripe movement is represented by the following formula:

[0039] in, The wavelength of the laser. The optical delay time generated for the etalon. The system determines this delay time based on user input. VPF and known The required optical delay time can be derived by using the above formula. The system then corrects the physical thickness of the etalon accordingly. Subsequently, the system simulates the reflection process of the probe beam on the surface of the moving target. Based on the Doppler effect, it calculates the frequency shift of the probe beam according to the transient velocity of the target surface and generates a modulated signal light carrying velocity information.

[0040] Then, step S3 is executed to simulate the interferometric detection process. The system inputs the modulated signal light into the interferometer model and splits it into two paths. The etalon delay determined in step S2 is then used. This causes an optical delay in one of the signal beams. Simultaneously, a pre-defined linear phase shift is introduced along the spatial axis to simulate the tilted interference effect of the initial reference fringes in the experiment. The two signal beams are then superimposed, and the total electric field after superposition is calculated. E Represented as:

[0041] in, , These represent the amplitudes of the two signal beams, , These are the angular frequencies of the two signal beams, respectively. , These represent the initial phases of the two signal beams. According to the principle of wave superposition, the instantaneous energy received by the detector is equal to the square of the total electric field. Proportional. After expanding and ignoring high-frequency terms, the effective light intensity distribution is obtained. Its formula is:

[0042] This process transforms unobservable frequency variations into recordable light intensity distributions, forming ideal interference fringes.

[0043] Finally, S4 is executed to perform temporal blurring of the interferometric image and generate the result. Since the slit of a real fringe camera has a physical width, a time integration effect occurs when the light beam passes through. To simulate this process, the system first calculates the average difference between adjacent time points on the time axis array to obtain the average simulation step size. dt The physical response time of the slit in the streak camera is defined as... The number of pixels occupied by the slit width is obtained by dividing the average step size downwards. tpx The calculation formula is:

[0044] The system build length is tpx All-1 window function w The convolution operator is used to perform a one-dimensional convolution operation on the ideal interference intensity matrix generated by S3 along the time axis using a window function. The convolution formula is as follows:

[0045] in, For ideal interference light intensity, y For spatial location, t For time, m The convolution kernel is indexed. Finally, the system generates and displays three result images: an interference fringe pattern with time on the horizontal axis and space on the vertical axis; a three-dimensional velocity field pattern with time and space as the base; and a two-dimensional pseudo-color image reflecting the velocity distribution.

[0046] This invention also provides a system for generating interferometric optical field distribution characteristics for an arbitrary reflecting surface velocity diagnostic system. This system is used to implement the steps in the above method embodiments, and its structure corresponds to... Figure 1 The simulation process logic is shown below.

[0047] The system includes a velocity field definition module, a Doppler modulation signal generation module, an interferometric detection module, and a time ambiguity processing module.

[0048] The velocity field definition module is used to define the motion state distribution of the shock wave under test. Specifically, it receives a mode selection command from the user. If the standard mode is selected, it calls preset typical velocity field functions, including static velocity field functions, smooth velocity field functions, step velocity field functions, and sinusoidal velocity field functions, supporting user settings for parameters such as amplitude and frequency. If the custom mode is selected, it imports an external three-dimensional velocity field data file and maps the spatiotemporal coordinates and velocity values ​​from the external data to the simulation mesh using a trilinear interpolation algorithm. The velocity field definition module further includes a preset velocity field calling unit and an external data import unit: the preset velocity field calling unit stores the aforementioned four typical velocity field functions; the external data import unit parses the external file and extracts the spatiotemporal distribution information of the velocity field.

[0049] The Doppler modulation signal generation module configures parameters such as laser wavelength, etalon refractive index, initial etalon thickness, and fringe velocity constant (VPF). It calculates the Doppler frequency shift based on the motion state distribution output by the velocity field definition module, generating a modulated signal light with velocity field information. This module has a built-in VPF calculation unit, using the formula... Optical delay time generated by reverse etalon The physical thickness of the etalon is then adjusted accordingly.

[0050] The interferometric detection module splits the modulated signal light into two paths, causing one path to experience optical delay, which then superimposes with the other path to form interference fringes, thus obtaining an ideal interference intensity matrix. This module is based on the formula... Calculate the interference light intensity, where , These represent the amplitudes of the two signal beams, , These are the angular frequencies of the two signal beams, respectively. , These represent the initial phases of the two signal beams. t It is a time variable.

[0051] The temporal blur processing module is used to simulate the slit physical response characteristics of a streak camera. It includes a convolution kernel building unit and a convolution operation unit: the convolution kernel building unit calculates the average simulation stride of the time axis array. dt Based on the slit physical response time of the streak camera Through formula Get the number of pixels occupied by the slit width tpx and construct a length of tpx All-1 window function w The convolution operation unit operates along the time axis on the ideal interference intensity matrix. To perform a one-dimensional convolution operation, the formula is: ,in, y Represents spatial location coordinates,t Represents the time coordinate. m The kernel index variable is used to output a simulated interference fringe image.

[0052] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A method for generating the interference optical field distribution characteristics of an arbitrary reflective surface velocity diagnostic system, characterized in that: Includes the following steps: S1: Define the motion state distribution of the shock wave to be measured. The motion state distribution is obtained by calling a preset typical velocity field function or importing externally defined three-dimensional velocity field data. S2: Configure the parameters of laser wavelength, etalon refractive index, etalon initial thickness and stripe velocity constant, calculate the Doppler frequency shift according to the motion state distribution, and generate modulated signal light with velocity field information; S3: The modulated signal light is split into two paths, one of which is optically delayed and then superimposed with the other to form interference fringes, thus obtaining an ideal interference light intensity matrix; S4: Based on the physical response time of the slit of the fringe camera, perform a one-dimensional convolution operation on the ideal interference intensity matrix in the time axis direction to generate a simulated interference fringe image.

2. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 1, characterized in that: The typical velocity field function described in S1 includes one or more of the following: static velocity field function, smooth velocity field function, step velocity field function, and sinusoidal velocity field function. The externally defined three-dimensional velocity field data is imported after preprocessing by numerical calculation software.

3. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 1, characterized in that: The stripe velocity constant described in S2 VPF The calculation formula is: in, The wavelength of the laser. The optical delay time generated by the etalon.

4. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 3, characterized in that: S2 also includes: based on the settings VPF and The optical delay time of the etalon is calculated using the formula. and according to Correct the physical thickness of the etalon.

5. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 1, characterized in that: The formula for calculating the light intensity distribution of the interference fringes described in S3 is: in, , These represent the amplitudes of the two signal beams, , These are the angular frequencies of the two signal beams, respectively. , These represent the initial phases of the two signal beams. t It is a time variable.

6. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 1, characterized in that: S4 specifically includes: calculating the average simulation step size of the time axis array. dt Based on the slit physical response time of the streak camera Calculate the number of pixels occupied by the slit width tpx The tpx The calculation formula is: in, This indicates the floor function.

7. The method for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 6, characterized in that: S4 also includes: a construction length of tpx All-1 window function w Along the time axis, the ideal interference intensity matrix Perform a one-dimensional convolution operation to obtain a simulated interference fringe image. The convolution operation formula is: in, y Represents spatial location coordinates, t Represents the time coordinate. m This is the kernel index variable.

8. A system for generating the interference optical field distribution characteristics of an arbitrary reflecting surface velocity diagnostic system, characterized in that: include: The velocity field definition module is used to define the motion state distribution of the shock wave to be measured. The motion state distribution is obtained by calling a preset typical velocity field function or importing externally defined three-dimensional velocity field data. The Doppler modulation signal generation module is used to configure parameters such as laser wavelength, etalon refractive index, etalon initial thickness and fringe velocity constant, calculate Doppler frequency shift according to the motion state distribution, and generate modulation signal light with velocity field information. An interference detection module is used to split the modulated signal light into two paths, so that one path generates an optical delay and is superimposed with the other path to form interference fringes, thereby obtaining an ideal interference light intensity matrix. The temporal blur processing module is used to perform a one-dimensional convolution operation on the ideal interference intensity matrix in the time axis direction based on the physical response time of the slit of the fringe camera, so as to generate a simulated interference fringe image.

9. The system for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 8, characterized in that: The velocity field definition module includes a preset velocity field calling unit and an external data import unit. The preset velocity field calling unit is used to store static velocity field functions, smooth velocity field functions, step velocity field functions, and sinusoidal velocity field functions. The external data import unit is used to receive and parse externally customized three-dimensional velocity field data files.

10. The system for generating the interference optical field distribution characteristics of the arbitrary reflecting surface velocity diagnostic system according to claim 8, characterized in that: The temporal blurring module includes a convolution kernel construction unit and a convolution operation unit. The convolution kernel construction unit is used to calculate the number of pixels occupied by the slit width based on the slit physical response time of the streak camera. tpx and construct a length of tpx All-1 window function w The convolution operation unit is used to perform a one-dimensional convolution operation on the ideal interference light intensity matrix along the time axis.

Images