Accelerated simulation method of gradient period polarization holographic grating and related device

Abstract

The invention discloses an accelerated simulation method of a gradient period polarization holographic grating and a related device, and the method comprises the steps: carrying out the D (kin) calculation of the direction of predicted incident light after passing through a grating layer, multiplying a calculation result with the amplitude of the corresponding incident light, calculating the diffraction efficiency, and carrying out the acceleration simulation of the gradient period polarization holographic grating. Compared with a traditional simulation method using complete strict coupling wave calculation, the method has the advantages that the process operation of multiple integration, matrix multiplication and matrix inversion is reduced, and the simulation speed can be greatly accelerated.

Description

Technical Field

[0001] The invention relates to an accelerated simulation method of a gradient period polarization body holographic grating and a related device, belonging to the field of holographic optical material simulation. Background Art

[0002] As the next generation of intelligent display technology, augmented reality (AR) has attracted widespread attention because it can provide an immersive and practical experience, enhance interactivity and visual communication, and improve efficiency, quality and user engagement. In AR display technology, polarization volume grating (PVG) has unique advantages in the fields of optics and optoelectronics, and is particularly suitable for augmented reality display systems.

[0003] In order to meet the growing field of view requirements, a gradient period polarization holographic grating has been proposed in recent years. Compared with the traditional fixed period PVG, the gradient period PVG enhances its diffraction characteristics and wavelength bandwidth capabilities. Simulating light passing through a gradient period polarization holographic grating has become a research direction that has attracted much attention. It involves using computer programming to simulate the behavior of light, that is, simulating a gradient period polarization holographic grating.

[0004] Traditional simulation methods using rigorous coupled wave analysis, such as simulation in Zemax software, require multiple integrations, matrix multiplications, and matrix inversions. In addition, since the gradually varying period PVG usually needs to be divided into multiple layers for simulation, the simulation efficiency is greatly reduced compared to the fixed period PVG, which makes it difficult to optimize the parameters of the optical waveguide design. Summary of the invention

[0005] The present invention provides an accelerated simulation method for a gradient period polarization holographic grating and a related device, which solves the problem of low simulation efficiency of the traditional method.

[0006] According to one aspect of the present disclosure, there is provided a method for accelerating simulation of a gradient period polarization volume holographic grating, comprising: The gradient period polarization holographic grating to be simulated is divided, and a mixing matrix of the gradient period polarization holographic grating is constructed according to each divided grating layer; wherein the mixing matrix includes amplitude phase factors in an electric field forward propagation mode, an electric field reverse propagation mode, a magnetic field forward propagation mode, and a magnetic field reverse propagation mode; The mixed matrix of the gradient period polarization volume holographic grating and the preset wave vector direction k in and amplitude E in The electromagnetic field boundary conditions are matched with the incident light, and the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated according to the matched incident light information and the mixing matrix of the gradient period polarization holographic grating; The direction of the incident light passing through the grating layer is deduced according to the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer and the preset grating parameters; D(k in ) Accelerate computing; D(k in ) Acceleration calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

[0007] In some embodiments of the present disclosure, a hybrid matrix of a graded period polarization volume holographic grating is constructed according to each grating layer after slicing, including: Each grating layer after slicing is regarded as a single-layer polarization volume holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarization volume holographic grating is constructed; According to the mixing matrix of the single-layer polarization volume holographic grating, the mixing matrix of the gradient period polarization volume holographic grating is obtained.

[0008] In some embodiments of the present disclosure, the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated using the formula: ; ; In the formula, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, e > 、e < 、h > and h < The W matrix is ​​split into sub-blocks in the electric field forward propagation mode, sub-blocks in the electric field reverse propagation mode, sub-blocks in the magnetic field forward propagation mode, and sub-blocks in the magnetic field reverse propagation mode. W is the eigenvector matrix derived from the single-layer polarizer holographic grating using the rigorous coupled wave algorithm.

[0009] In some embodiments of the present disclosure, the diffraction efficiency of the grating layer to the incident light is calculated using the formula: ; In the formula, E R is the amplitude of the reflected light, E Tis the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting, H is the mixing matrix of the graded period polarization holographic grating.

[0010] According to another aspect of the present disclosure, there is provided an accelerated simulation device for a gradient period polarization volume holographic grating, comprising: The global mixing matrix module divides the gradient period polarization holographic grating to be simulated, and constructs a mixing matrix of the gradient period polarization holographic grating according to each divided grating layer; wherein the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode and the magnetic field reverse propagation mode; The coupling relationship module combines the mixed matrix of the gradient period polarization volume holographic grating and the preset wave vector direction k in and amplitude E in The electromagnetic field boundary conditions are matched with the incident light, and the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated according to the matched incident light information and the mixing matrix of the gradient period polarization holographic grating; A direction module, which deduce the direction of the incident light passing through the grating layer according to the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer and the preset grating parameters; D(k in ) acceleration module, which predicts the direction of the incident light after passing through the grating layer. in ) Accelerate computing; Diffraction efficiency module, D(k in ) Acceleration calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

[0011] In some embodiments of the present disclosure, in the global mixing matrix module, a mixing matrix of a gradient period polarization volume holographic grating is constructed according to each grating layer after segmentation, including: Each grating layer after slicing is regarded as a single-layer polarization volume holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarization volume holographic grating is constructed; According to the mixing matrix of the single-layer polarization volume holographic grating, the mixing matrix of the gradient period polarization volume holographic grating is obtained.

[0012] In some embodiments of the present disclosure, in the coupling relationship module, the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated, and the formula is: ; ; In the formula, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, e > 、e < 、h > and h < The W matrix is ​​split into sub-blocks in the electric field forward propagation mode, sub-blocks in the electric field reverse propagation mode, sub-blocks in the magnetic field forward propagation mode, and sub-blocks in the magnetic field reverse propagation mode. W is the eigenvector matrix derived from the single-layer polarizer holographic grating using the rigorous coupled wave algorithm.

[0013] In some embodiments of the present disclosure, in the diffraction efficiency module, the diffraction efficiency of the grating layer to the incident light is calculated using the formula: ; In the formula, E R is the amplitude of the reflected light, E T is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting, H is the mixing matrix of the graded period polarization holographic grating.

[0014] According to another aspect of the present disclosure, a computer-readable storage medium is provided, which stores one or more programs, and the one or more programs include instructions. When the instructions are executed by a computing device, the computing device performs an accelerated simulation method of a gradient period polarization holographic grating.

[0015] According to another aspect of the present disclosure, a computer device is provided, comprising one or more processors and one or more memories, wherein one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs include instructions for executing an accelerated simulation method of a gradient period polarization holographic grating.

[0016] The beneficial effect achieved by the present invention is that the present invention predicts the direction of the incident light after passing through the grating layer D(k in ) is calculated, and then the result is multiplied with the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method using complete rigorous coupled wave calculation, it reduces the process operations of multiple integrations, matrix multiplications and matrix inversions, and can significantly accelerate the simulation speed. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 The flowchart of the accelerated simulation method of the gradient period polarization volume holographic grating; Figure 2 Schematic diagram of a common shape of a gradient periodic polarization holographic grating; Figure 3 It is a cross-sectional view of the interior of the holographic grating with a gradient periodic polarization; Figure 4 It is the relationship between the speedup ratio and the number of grating layers; Figure 5 The block diagram of the accelerated simulation device for the gradient period polarization holographic grating. DETAILED DESCRIPTION

[0018] The following will be combined with the drawings in the embodiments of the present disclosure to clearly and completely describe the technical solutions in the embodiments of the present disclosure. It is obvious that the described embodiments are only part of the embodiments of the present disclosure, rather than all the embodiments. The following description of at least one exemplary embodiment is actually only illustrative and is by no means intended to limit the present disclosure and its application or use. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present disclosure.

[0019] Unless otherwise specified, the relative arrangement of components and steps, the numerical expressions and numerical values ​​set forth in these embodiments do not limit the scope of the present disclosure.

[0020] At the same time, it should be understood that for the convenience of description, the sizes of the various parts shown in the drawings are not drawn according to the actual proportional relationship.

[0021] Technologies, methods, and equipment known to ordinary technicians in the relevant art may not be discussed in detail, but where appropriate, the technologies, methods, and equipment should be considered part of the specification.

[0022] In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limiting. Therefore, other examples of the exemplary embodiments may have different values.

[0023] It should be noted that similar symbols and letters refer to similar items in the following figures, and therefore, once an item is defined in one figure, it does not require further discussion in subsequent figures.

[0024] In order to solve the problem of low simulation efficiency of traditional methods, the present disclosure proposes an accelerated simulation method for a gradient period polarization holographic grating.

[0025] Figure 1 It is a schematic diagram of an embodiment of the accelerated simulation method of the gradient period polarization volume holographic grating disclosed in the present invention, Figure 1 The embodiments may be executed by an emulated terminal (such as a computer).

[0026] like Figure 1 As shown, step 1 of the embodiment is to divide the gradient period polarization holographic grating to be simulated, and construct a mixing matrix of the gradient period polarization holographic grating according to each grating layer after division; wherein the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode and the magnetic field reverse propagation mode.

[0027] It should be noted that before simulation, a series of definitions will be performed on the gradient period polarization holographic grating to be simulated, such as defining its parameters, shape, etc.

[0028] The defined parameters may include vertex, length, width, thickness above the grating, thickness below the grating, refractive index of the waveguide and the outside air, determination of the state of the light on the waveguide, position and amplitude change of the light after it is emitted from the exit surface, etc.; wherein, the position change may include the position change of the light reaching the upper wall and the lower wall of the waveguide after total reflection occurs, and the amplitude change may include calculating the amplitude change of the TE and TM components in the light.

[0029] See Figure 2 and 3 The defined shape is mainly to determine whether the gradient period polarizer holographic grating to be simulated is a trapezoidal grating, a circular grating, a triangular grating, etc.; the specific shape defined here can be selected in the established grating library. The positions of the four vertices of the trapezoidal grating, the center position and radius of the circular grating, the initialization of the circular grating parameters, and the result of the light emitting from the exit surface through the gradient period polarizer holographic grating, etc.

[0030] In some embodiments, constructing a mixed matrix of a gradient period polarization volume holographic grating according to each grating layer after segmentation may specifically include: 1) Each grating layer after slicing is regarded as a single-layer polarization volume holographic grating with a constant longitudinal period, and a mixing matrix of the single-layer polarization volume holographic grating is constructed.

[0031] It should be noted that the graded period polarization volume holographic grating can be divided into a single-layer polarization volume holographic grating with a fixed longitudinal period, and the mixing matrix of the single-layer polarization volume holographic grating is constructed by the formula of the coupling relationship module.

[0032] The hybrid matrix includes the matrix of the amplitude phase factor in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode, the magnetic field reverse propagation mode, the dielectric tensor Fourier series, etc. In order to describe the relative conversion between the forward and reverse propagation modes, the W matrix can be split into four sub-blocks of two rows and two columns, specifically the sub-block in the electric field forward propagation mode, the sub-block in the electric field reverse propagation mode, the sub-block in the magnetic field forward propagation mode, and the sub-block in the magnetic field reverse propagation mode, which can be represented by the letter e > 、e < 、h > 、h < , W is the eigenvector matrix derived from the single-layer polarization volume holographic grating using the rigorous coupled wave algorithm, which is used to represent the amplitudes corresponding to various propagation modes in the electromagnetic field.

[0033] 2) According to the mixing matrix of the single-layer polarization holographic grating, the mixing matrix of the gradient period polarization holographic grating is obtained; wherein, the mixing matrix of the gradient period polarization holographic grating is similar to the matrix of the single layer, and also includes the matrix of the amplitude phase factor in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode, and the magnetic field reverse propagation mode, the matrix of the dielectric tensor Fourier series, etc. It should be noted that the mixing matrix of the gradient period polarization holographic grating is initialized first, and then the elements in the mixing matrix of each single-layer polarization holographic grating are used to replace the elements in the initialized mixing matrix, and finally the final mixing matrix of the gradient period polarization holographic grating is obtained.

[0034] return Figure 1 , step 2 of the embodiment, the mixed matrix of the gradient period polarization holographic grating and the preset wave vector direction k in and amplitude E in The incident light is matched with the electromagnetic field boundary conditions, that is, the mixed matrix of the gradient period polarization holographic grating and the preset wave vector direction k in and amplitude E inSubstitute the incident light into the electromagnetic field boundary conditions for calculation, and calculate the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer according to the matching incident light information and the mixing matrix of the gradient period polarization holographic grating, that is, calculate the mixing matrix according to the formula of the coupling relationship module H .

[0035] It should be noted that the incident light information is predefined, and specifically the structure, amplitude, direction (projection in each direction), position, energy, wavelength independent variables, etc. of the light may be defined.

[0036] The coupling relationship between the electromagnetic fields on the two surfaces of the grating layer can be expressed as: ; ; In the formula, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, respectively.

[0037] return Figure 1 , step 3 of the embodiment, predicts the direction of the incident light after passing through the grating layer based on the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer.

[0038] It should be noted that in a waveguide system, different light rays emitted from the same pixel in the image source have the same incident direction on the input coupling grating after passing through the collimating lens. Therefore, the incident direction of such light rays is limited and predictable.

[0039] return Figure 1 , step 4 of the embodiment, predicting the direction of the incident light after passing through the grating layer D(k in ) calculation, namely the multiplication and inverse matrix operations in the diffraction efficiency calculation module. The computationally intensive matrix D depends only on k in .

[0040] return Figure 1 , step 5 of the embodiment, D(k in ) Calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

[0041] It should be noted that, specifically, D(k in ) Calculation results and corresponding amplitude E inMultiplying, the formula can be expressed as: ; In the formula, E R is the amplitude of the reflected light, E T is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting is H Split into a matrix of two rows and two columns.

[0042] The present invention first predicts the direction of the incident light after passing through the grating layer D(k in ) is calculated, and then the calculation result is multiplied with the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method using complete rigorous coupled wave calculation, the process operations of multiple integrations, matrix multiplications and matrix inversions are reduced. Therefore, the acceleration method proposed in the present invention can greatly accelerate the simulation speed.

[0043] In order to verify the above method, the above method was compared with the traditional strict coupled wave method. The verification experiment was to simulate and calculate 10,000 beams of light with different amplitudes in the same direction on a gradient periodic polarization holographic grating with different layers. The grating layers were simulated as 3, 5, 10, 15, 20, 25, and 30 respectively. When the number of simulated grating layers was 3, the above method took 12.15 seconds, while the traditional method took 649.49 seconds. The relationship between the acceleration ratio and the number of grating layers obtained by dividing the time consumed by the traditional method by the time consumed by the above method is shown in Figure 4 From the above comparison, it can be seen that the simulation efficiency of the above method far exceeds that of the traditional method.

[0044] Figure 5 It is a schematic diagram of an embodiment of an accelerated simulation device for a gradient period polarization volume holographic grating disclosed in the present invention, Figure 5 The embodiment is a virtual device that can be loaded and executed by a simulation terminal (such as a computer), including a global mixing matrix module, a coupling relationship module, a direction module, D(k in ) module and the diffraction efficiency module.

[0045] The global mixing matrix module of the embodiment is configured to divide the gradient period polarization holographic grating to be simulated, and construct a mixing matrix of the gradient period polarization holographic grating according to each grating layer after division; wherein the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode and the magnetic field reverse propagation mode.

[0046] It should be noted that in the global mixing matrix module, a mixing matrix of a gradually changing period polarization volume holographic grating is constructed according to each grating layer after slicing, including: treating each grating layer after slicing as a single-layer polarization volume holographic grating with a constant longitudinal period, constructing a mixing matrix of a single-layer polarization volume holographic grating; and obtaining a mixing matrix of a gradually changing period polarization volume holographic grating according to the mixing matrix of the single-layer polarization volume holographic grating.

[0047] The coupling relationship module of the embodiment is configured to combine the mixed matrix of the gradient period polarization volume holographic grating and the preset wave vector direction k in and amplitude E in The electromagnetic field boundary conditions are matched with the incident light, and the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated based on the matched incident light information and the mixing matrix of the gradient period polarization holographic grating.

[0048] It should be noted that in the coupling relationship module, the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated using the formula: ; ; In the formula, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, e > 、e < 、h > and h < The W matrix is ​​split into sub-blocks in the electric field forward propagation mode, sub-blocks in the electric field reverse propagation mode, sub-blocks in the magnetic field forward propagation mode, and sub-blocks in the magnetic field reverse propagation mode. W is the eigenvector matrix derived from the single-layer polarizer holographic grating using the rigorous coupled wave algorithm.

[0049] The direction module of the embodiment is configured to predict the direction of the incident light after passing through the grating layer according to the coupling relationship between the electromagnetic fields on two surfaces of the grating layer.

[0050] D(k) of the embodiment in ) module is configured to predict the direction of the incident light after passing through the grating layer D(k in )calculate.

[0051] The diffraction efficiency module of the embodiment is configured to convert D(k in ) Calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

[0052] It should be noted that in the diffraction efficiency module, the diffraction efficiency of the grating layer to the incident light is calculated using the formula: ; In the formula, E R is the amplitude of the reflected light, E T is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting, H is the mixing matrix of the graded period polarization holographic grating.

[0053] Similar to the above method, the above device first predicts the direction of the incident light after passing through the grating layer D(k in ) is calculated, and then the calculation result is multiplied with the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method using complete rigorous coupled wave calculation, the process operations of multiple integrations, matrix multiplications and matrix inversions are reduced. Therefore, the acceleration method proposed in the present invention can greatly accelerate the simulation speed.

[0054] Based on the same technical solution, the present disclosure also relates to a computer-readable storage medium, which stores one or more programs, and the one or more programs include instructions. When the instructions are executed by a computing device, the computing device performs an accelerated simulation method for a gradient period polarization holographic grating.

[0055] Based on the same technical solution, the present disclosure also relates to a computer device, including one or more processors and one or more memories, one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs include instructions for executing an accelerated simulation method of a gradient period polarization holographic grating.

[0056] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0057] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0058] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0059] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0060] The above are merely embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention are included in the scope of the claims of the present invention to be approved.

Claims

1. An accelerated simulation method for a gradient period polarization holographic grating, characterized in that: include: The gradient period polarization holographic grating to be simulated is divided, and a mixing matrix of the gradient period polarization holographic grating is constructed according to each divided grating layer; wherein the mixing matrix includes amplitude phase factors in an electric field forward propagation mode, an electric field reverse propagation mode, a magnetic field forward propagation mode, and a magnetic field reverse propagation mode; The mixed matrix of the gradient period polarization volume holographic grating and the preset wave vector direction k in and amplitude E in The incident light is matched using the electromagnetic field boundary conditions, and the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated based on the matched incident light information matrix and the mixed matrix of the gradient period polarization holographic grating; According to the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer and the preset grating parameters, the direction of the incident light after passing through the grating layer is deduced; D(k in )calculate; D(k in ) Calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

2. The accelerated simulation method of the gradient period polarization holographic grating according to claim 1 is characterized in that: According to each grating layer after slicing, a hybrid matrix of a gradient period polarization volume holographic grating is constructed, including: Each grating layer after slicing is regarded as a single-layer polarization volume holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarization volume holographic grating is constructed; According to the mixing matrix of the single-layer polarization volume holographic grating, the mixing matrix of the gradient period polarization volume holographic grating is obtained.

3. The accelerated simulation method of the gradient period polarization holographic grating according to claim 1, characterized in that: Calculate the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer, the formula is: ； ； Where, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, e > 、e < 、h > and h < The W matrix is ​​split into sub-blocks in the electric field forward propagation mode, sub-blocks in the electric field reverse propagation mode, sub-blocks in the magnetic field forward propagation mode, and sub-blocks in the magnetic field reverse propagation mode. W is the eigenvector matrix derived from the single-layer polarizer holographic grating using the rigorous coupled wave algorithm.

4. The accelerated simulation method of a gradient period polarization holographic grating according to claim 1, characterized in that: Calculate the diffraction efficiency of the grating layer for the incident light, the formula is: ； In the formula, E R is the amplitude of the reflected light, E T is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting, H is the mixing matrix of the graded period polarization holographic grating.

5. An accelerated simulation device for a gradient period polarization holographic grating, characterized in that: include: The global mixing matrix module divides the gradient period polarization holographic grating to be simulated, and constructs a mixing matrix of the gradient period polarization holographic grating according to each divided grating layer; wherein the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field reverse propagation mode, the magnetic field forward propagation mode and the magnetic field reverse propagation mode; The coupling relationship module combines the mixed matrix of the gradient period polarization volume holographic grating and the preset wave vector direction k in and amplitude E in The electromagnetic field boundary conditions are matched with the incident light, and the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated according to the matched incident light information matrix and the mixed matrix of the gradient period polarization holographic grating; A direction module, which deduce the direction of the incident light after passing through the grating layer according to the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer and the preset grating parameters; D(k in ) acceleration module, which predicts the direction of the incident light after passing through the grating layer. in ) Accelerate computing; Diffraction efficiency module, D(k in ) Calculation results and corresponding amplitude E in Multiply them together to calculate the diffraction efficiency of the grating layer for the incident light.

6. The accelerated simulation device for a gradient period polarization holographic grating according to claim 5, characterized in that: In the global mixing matrix module, a mixing matrix of a gradient period polarization volume holographic grating is constructed according to each grating layer after segmentation, including: Each grating layer after slicing is regarded as a single-layer polarization volume holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarization volume holographic grating is constructed; According to the mixing matrix of the single-layer polarization volume holographic grating, the mixing matrix of the gradient period polarization volume holographic grating is obtained.

7. The accelerated simulation device for a gradient period polarization holographic grating according to claim 5, characterized in that: In the coupling relationship module, the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer is calculated. The formula is: ； ； Where, E > and E < The electric field in the forward and reverse propagation modes, H > and H < The magnetic field in forward and reverse propagation modes, respectively, H is the mixing matrix of the graded period polarization holographic grating, p > (d) and p < (-d) < are the amplitude phase factors in the forward and reverse propagation modes, e > 、e < 、h > and h < The W matrix is ​​split into sub-blocks in the electric field forward propagation mode, sub-blocks in the electric field reverse propagation mode, sub-blocks in the magnetic field forward propagation mode, and sub-blocks in the magnetic field reverse propagation mode. W is the eigenvector matrix derived from the single-layer polarizer holographic grating using the rigorous coupled wave algorithm.

8. The accelerated simulation device for a gradient period polarization holographic grating according to claim 5, characterized in that: In the diffraction efficiency module, the diffraction efficiency of the grating layer to the incident light is calculated using the formula: ； In the formula, E R is the amplitude of the reflected light, E T is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the unit matrix, H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint cutting, H is the mixing matrix of the graded period polarization holographic grating.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores one or more programs, and the one or more programs include instructions. When the instructions are executed by a computing device, the computing device executes any one of the methods of claims 1 to 4.

10. A computer device, characterized in that: include: One or more processors and one or more memories, one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs include instructions for executing any of the methods described in claims 1 to 4.

Images