Modeling method and device for propagation coefficient of orbital angular momentum in multilayer board

Abstract

The invention belongs to the technical field of wireless communication, and provides a propagation coefficient modeling method and device of orbital angular momentum in a multilayer board. The method comprises the following steps: constructing an ITU multilayer board iterative propagation coefficient model; determining a vertical polarization reflection coefficient, a horizontal polarization reflection coefficient, a vertical polarization transmission coefficient and a horizontal polarization transmission coefficient according to an ITU multilayer board iterative propagation coefficient model; according to the vertical polarization reflection coefficient, constructing a vertical polarization multilayer board Laguerre reflection model; constructing a Laguerre reflection model of the horizontal polarization multilayer board according to the horizontal polarization reflection coefficient; constructing a vertical polarization multilayer board Laguerre transmission model according to the vertical polarization transmission coefficient; and constructing a Laguerre transmission model of the horizontal polarization multilayer board according to the horizontal polarization transmission coefficient. According to the method, a clear OAM propagation model can be constructed, and more accurate OAM propagation coefficients, including reflection coefficients and transmission coefficients, of vortex electromagnetic waves can be obtained.

Description

Technical Field

[0001] The present invention relates to the field of wireless communication technology, and in particular to a method and device for modeling a propagation coefficient of orbital angular momentum in a multilayer plate. Background Art

[0002] Orbital Angular Momentum (OAM) is an inherent physical quantity of electromagnetic waves, and also a new dimension of wireless transmission, and one of the potential key technologies of 6G. Vortex electromagnetic waves carrying orbital angular momentum are a new physical dimension and transmission carrier that is different from traditional plane electromagnetic waves. They have a large number of orthogonal OAM eigenstates, and can achieve high-quality and large-capacity information transmission on limited channel resources in a new way of information transmission, which revolutionizes the bottleneck problems of 6G spectrum efficiency, transmission capacity and anti-interference ability.

[0003] Indoor communication environment is one of the hot spots for 6G deployment. Indoor environments usually contain various forms of scatterers, which cause electromagnetic waves to propagate through reflection, scattering, transmission or diffraction. Since deterministic modeling relies on the propagation characterization of the interaction between electromagnetic waves and materials, the propagation mechanism of electromagnetic waves and scatterers is crucial for indoor communication deployment.

[0004] Since OAM waves and plane electromagnetic waves have different propagation characteristics, the propagation models of plane waves such as reflection, transmission, scattering and diffraction are vaguely applicable in OAM channel modeling. Directly using plane wave models and parameters for modeling will affect the accuracy of the OAM channel model. A clear understanding of the OAM propagation model is a prerequisite for accurate OAM channel modeling, and channel modeling is an important basis for the design and deployment of OAM indoor communications. Therefore, the basis for the smooth deployment of OAM technology in indoor environments is to clarify the interaction mechanism between OAM waves and building materials, so as to build a clear OAM propagation model. Summary of the invention

[0005] In view of the defects in the prior art, the present invention provides a method and device for modeling the propagation coefficient of orbital angular momentum in a multilayer plate, so as to construct a clear OAM propagation model and obtain a more accurate OAM propagation coefficient of vortex electromagnetic waves.

[0006] A method for modeling propagation coefficients of orbital angular momentum in a multilayer plate, characterized by comprising:

[0007] Construct an iterative propagation coefficient model for ITU multilayer boards;

[0008] Determine the vertical polarization reflection coefficient, horizontal polarization reflection coefficient, vertical polarization transmission coefficient and horizontal polarization transmission coefficient according to the ITU multilayer iterative propagation coefficient model;

[0009] The vertically polarized multi-layer plate Laguerre reflection model is constructed based on the vertically polarized reflection coefficient.

[0010] The Laguerre reflection model of horizontally polarized multilayer plates is constructed based on the horizontally polarized reflection coefficient.

[0011] The Laguerre transmission model of vertically polarized multilayer plates is constructed based on the vertically polarized transmission coefficient.

[0012] The Laguerre transmission model of horizontally polarized multilayer plates is constructed based on the horizontally polarized transmission coefficient.

[0013] Furthermore, the ITU multilayer iterative propagation coefficient model is characterized by the following formula:

[0014] A n =0.5exp(jk n d n cosθ n )[A n+1 (1+Y n+1 )+B n+1 (1-Y n+1 )];

[0015] B n =0.5exp(-jk n d n cosθ n )[A n+1 (1-Y n+1 )+B n+1 (1+Y n+1 )];

[0016] F n =0.5exp(jk n d n cosθ n )[F n+1 (1+W n+1 )+G n+t (1-W n+1 )];

[0017] G n =0.5exp(-jk n d n cosθ n )[F n+1 (1-W n+1 )+G n+1 (1+W n+1 )];

[0018] Among them, A n , B n , F n , G n , W n+1 , Yn+1 , k n is the constant term expression of the recursive algorithm, j is a complex unit, θ n is the propagation direction of the nth layer, d n is the width of the multilayer board.

[0019] Furthermore, the vertical polarization reflection coefficient, the horizontal polarization reflection coefficient, the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient are characterized by the following formulas respectively:

[0020]

[0021] Among them, R ⊥ is the vertical polarization reflection coefficient, R P is the horizontal polarization reflection coefficient, T ⊥ is the vertical polarization transmission coefficient, T P is the horizontal polarization transmission coefficient, and A0, B0, G0, and F0 represent the four recursive parameters of the recursive algorithm.

[0022] Furthermore, a vertically polarized multilayer plate Laguerre reflection model is constructed according to the vertically polarized reflection coefficient, including:

[0023] The predicted reflection coefficient of the vertically polarized multilayer plate Laguerre reflection model is calculated based on the vertically polarized reflection coefficient using the method shown in the following formula:

[0024]

[0025] Where Γ1 is the predicted reflection coefficient of the vertically polarized multilayer Laguerre reflection model, E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0026] Furthermore, a horizontally polarized multilayer plate Laguerre reflection model is constructed according to the horizontally polarized reflection coefficient, including:

[0027] The predicted reflection coefficient of the horizontally polarized multilayer plate Laguerre reflection model is calculated based on the horizontally polarized reflection coefficient using the method shown in the following formula:

[0028]

[0029] Where Γ2 is the predicted reflection coefficient of the Laguerre reflection model for horizontally polarized multilayer plates, and E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0030] Furthermore, a vertically polarized multilayer plate Laguerre transmission model is constructed according to the vertically polarized transmission coefficient, including:

[0031] The predicted transmission coefficient of the vertically polarized multilayer Laguerre transmission model is calculated based on the vertically polarized transmission coefficient using the method shown in the following formula:

[0032]

[0033] Where T1 is the predicted transmission coefficient of the vertically polarized multilayer Laguerre reflection model, E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ t is the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0034] Furthermore, a horizontally polarized multilayer plate Laguerre transmission model is constructed according to the horizontally polarized transmission coefficient, including:

[0035] The predicted transmission coefficient of the horizontally polarized multilayer plate Laguerre transmission model is calculated based on the horizontally polarized transmission coefficient using the method shown in the following formula:

[0036]

[0037] Where T2 is the predicted transmission coefficient of the Laguerre transmission model for horizontally polarized multilayer plates, and E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ t is the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0038] A propagation coefficient modeling device of orbital angular momentum in a multilayer plate comprises a first construction unit, a determination unit, a second construction unit, a third construction unit, a fourth construction unit, and a fifth construction unit, wherein:

[0039] The first construction unit is used to construct an ITU multilayer board iterative propagation coefficient model;

[0040] A determination unit, used for determining a vertical polarization reflection coefficient, a horizontal polarization reflection coefficient, a vertical polarization transmission coefficient and a horizontal polarization transmission coefficient according to an ITU multilayer iterative propagation coefficient model;

[0041] A second construction unit is used to construct a vertically polarized multi-layer plate Laguerre reflection model according to the vertically polarized reflection coefficient;

[0042] A third construction unit is used to construct a horizontally polarized multilayer plate Laguerre reflection model according to the horizontally polarized reflection coefficient;

[0043] A fourth construction unit is used to construct a vertically polarized multi-layer plate Laguerre transmission model according to the vertically polarized transmission coefficient;

[0044] The fifth construction unit is used to construct a horizontally polarized multilayer plate Laguerre transmission model according to the horizontally polarized transmission coefficient. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the specific embodiments of the present invention, the following will briefly introduce the drawings required for the specific embodiments or the prior art description. In all the drawings, each element or part is not necessarily drawn according to the actual scale.

[0046] Figure 1 A flow chart of a method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate provided by one embodiment of the present invention;

[0047] Figure 2 A schematic diagram of the structure of an OAM wavefront measurement system provided by one embodiment of the present invention;

[0048] Figure 3 A schematic diagram of the OAM propagation characteristic measurement principle provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0049] The following embodiments of the technical solution of the present invention are described in detail in conjunction with the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and are therefore only used as examples, and cannot be used to limit the protection scope of the present invention.

[0050] It should be noted that, unless otherwise specified, the technical terms or scientific terms used in this application should have the common meanings understood by those skilled in the art to which the present invention belongs.

[0051] like Figure 1 As shown, in one embodiment, a method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate is provided, comprising:

[0052] 1. Construct an ITU multilayer board iterative propagation coefficient model;

[0053] Preferably, the ITU multilayer iterative propagation coefficient model is characterized by the following formula:

[0054] A N+1 =1B N+1 =0 F N+1 =1 G N+1 =0

[0055] Then n=N,N-1,etc....0:

[0056] A n =0.5exp(jk n d n cosθ n )[A n+1 (1+Y n+1 )+B n+1 (1-Y n+1 )],

[0057] B n=0.5exp(-jk n d n cosθ n )[A n+1 (1-Y n+1 )+B n+1 (1+Y n+1 )],

[0058] F n =0.5exp(jk n dcosθ n )[F n+1 (1+W n+1 )+G n+t (1-W n+1 )],

[0059] G n =0.5exp(-jk n d n cosθ n )[F n+1 (1-W n+1 )+G n+1 (1+W n+1 )],

[0060]

[0061] Among them, A n , B n , F n , G n , W n+1 , Y n+1 , k n etc. are the constant term expressions of the recursive algorithm. When n=0, we get A0, B0, G0 and F0; j is a complex unit; θ n is the propagation direction of the nth layer, d n is the width of the multilayer board, ε rn is the relative dielectric constant of the nth layer, and λ is the wavelength.

[0062] 2. Determine the vertical polarization reflection coefficient, horizontal polarization reflection coefficient, vertical polarization transmission coefficient and horizontal polarization transmission coefficient based on the constructed ITU multilayer iterative propagation coefficient model;

[0063] The ITU multilayer iterative propagation coefficient model will determine the A n , B n , F n , G n When A n , B n , F n , G nWhen n=0, the values ​​of A0, B0, G0 and F0 are calculated, and then the vertical polarization reflection coefficient, horizontal polarization reflection coefficient, vertical polarization transmission coefficient and horizontal polarization transmission coefficient are determined according to the two polarization states of the electromagnetic wave (vertical polarization and horizontal polarization) and the constructed ITU multilayer iterative propagation coefficient model.

[0064] Optionally, the vertical polarization reflection coefficient, the horizontal polarization reflection coefficient, the vertical polarization transmission coefficient, and the horizontal polarization transmission coefficient are respectively characterized by the following formulas:

[0065]

[0066] Among them, R ⊥ is the vertical polarization reflection coefficient, R P is the horizontal polarization reflection coefficient, T ⊥ is the vertical polarization transmission coefficient, T P is the horizontal polarization transmission coefficient, and A0, B0, G0, and F0 represent the four recursive parameters of the recursive algorithm.

[0067] 3. Construct the Laguerre reflection model of vertically polarized multilayer plate based on the vertically polarized reflection coefficient;

[0068] The vertically polarized multilayer plate Laguerre reflection model is constructed based on the vertically polarized reflection coefficient, including:

[0069] The predicted reflection coefficient of the vertically polarized multilayer plate Laguerre reflection model is calculated based on the vertically polarized reflection coefficient using the method shown in the following formula:

[0070]

[0071] Where Γ1 is the predicted reflection coefficient of the vertically polarized multilayer Laguerre reflection model, E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0072] Preferably, the reflected electric field is characterized by the following formula:

[0073]

[0074] Among them, k r and Reflected wave vector in cylindrical coordinates, k zr The z-axis wave vector of the reflection surface, Z1 and Z2 represent the z-axis of the incident coordinate system and the reflection coordinate system, and d represents the horizontal propagation distance.

[0075] Preferably, the incident electric field is characterized by the following formula:

[0076]

[0077] Where ω0 is the initial beam waist radius, E0 represents the electric field at the initial beam waist radius ω0, l is the mode, is the Laguerre polynomial, and p is the radial exponent.

[0078] Preferably,

[0079] Among them, x i ,y i and z i It represents the propagation distance of the vortex electromagnetic wave in the x-axis, y-axis and z-axis of the incident electric field.

[0080] 4. Construct the Laguerre reflection model of horizontally polarized multilayer plate based on the horizontally polarized reflection coefficient;

[0081] The Laguerre reflection model of horizontally polarized multilayer plates is constructed based on the horizontally polarized reflection coefficient, including:

[0082] The predicted reflection coefficient of the horizontally polarized multilayer plate Laguerre reflection model is calculated based on the horizontally polarized reflection coefficient using the method shown in the following formula:

[0083]

[0084] Where Γ2 is the predicted reflection coefficient of the Laguerre reflection model for horizontally polarized multilayer plates, and E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z iis the z-axis propagation distance of the incident electric field.

[0085] Preferably, the reflected electric field is characterized by the following formula:

[0086]

[0087] Among them, k r and Reflected wave vector in cylindrical coordinates, k zr The z-axis wave vector of the reflection surface, Z1 and Z2 represent the z-axis of the incident coordinate system and the reflection coordinate system, and d represents the horizontal propagation distance.

[0088] Preferably, the incident electric field is characterized by the following formula:

[0089]

[0090] Where ω0 is the initial beam waist radius, E0 represents the electric field at the initial beam waist radius ω0, l is the mode, is the Laguerre polynomial, and p is the radial exponent.

[0091] Preferably,

[0092] Among them, x i ,y i and z i It represents the propagation distance of the vortex electromagnetic wave in the x-axis, y-axis and z-axis of the incident electric field.

[0093] 5. Construct the Laguerre transmission model of vertically polarized multilayer plates based on the vertically polarized transmission coefficient;

[0094] The vertically polarized multilayer plate Laguerre transmission model is constructed based on the vertically polarized transmission coefficient, including:

[0095] The predicted transmission coefficient of the vertically polarized multilayer Laguerre transmission model is calculated based on the vertically polarized transmission coefficient using the method shown in the following formula:

[0096]

[0097] Where T1 is the predicted transmission coefficient of the vertically polarized multilayer Laguerre reflection model, E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ t is the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0098] Preferably, the transmitted electric field is characterized by the following formula:

[0099]

[0100] Among them, k t and Transmitted wave vector in cylindrical coordinates, k zi is the z-axis wave vector of the incident plane, Z1 and Z3 represent the z-axis of the incident coordinate system and the transmission coordinate system.

[0101] Preferably, the incident electric field is characterized by the following formula:

[0102]

[0103] Where ω0 is the initial beam waist radius, E0 represents the electric field at the initial beam waist radius ω0, l is the mode, is the Laguerre polynomial, and p is the radial exponent.

[0104] Preferably,

[0105] Among them, x i ,y i and z i It represents the propagation distance of the vortex electromagnetic wave in the x-axis, y-axis and z-axis of the incident electric field.

[0106] 6. Construct the Laguerre transmission model of horizontally polarized multilayer plates based on the horizontally polarized transmission coefficient.

[0107] The Laguerre transmission model of horizontally polarized multilayer plates is constructed based on the horizontally polarized transmission coefficient, including:

[0108] The predicted transmission coefficient of the horizontally polarized multilayer plate Laguerre transmission model is calculated based on the horizontally polarized transmission coefficient using the method shown in the following formula:

[0109]

[0110] Where T2 is the predicted transmission coefficient of the Laguerre transmission model for horizontally polarized multilayer plates, and E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ tis the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

[0111] Preferably, the transmitted electric field is characterized by the following formula:

[0112]

[0113] Among them, k t and Transmitted wave vector in cylindrical coordinates, k zi is the z-axis wave vector of the incident plane, Z1 and Z3 represent the z-axis of the incident coordinate system and the transmission coordinate system.

[0114] Preferably, the incident electric field is characterized by the following formula:

[0115]

[0116] Where ω0 is the initial beam waist radius, E0 represents the electric field at the initial beam waist radius ω0, l is the mode, is the Laguerre polynomial, and p is the radial exponent.

[0117] Preferably,

[0118] Among them, x i ,y i and z i It represents the propagation distance of the vortex electromagnetic wave in the x-axis, y-axis and z-axis of the incident electric field.

[0119] In one of the embodiments, the accuracy of the predicted reflection coefficient (including the predicted reflection coefficient of the vertically polarized multilayer plate Laguerre reflection model and the predicted reflection coefficient of the horizontally polarized multilayer plate Laguerre reflection model) is calculated based on the actual reflection coefficient (including the actual vertical polarization reflection coefficient and the actual horizontal polarization reflection coefficient) obtained by actual measurement to obtain the fit between the actual reflection coefficient and the predicted reflection coefficient.

[0120] It also includes calculating the accuracy of the predicted transmission coefficient (including the predicted transmission coefficient of the vertically polarized multilayer plate Laguerre transmission model and the predicted transmission coefficient of the horizontally polarized multilayer plate Laguerre transmission model) based on the actual transmission coefficient (including the actual vertical polarization transmission coefficient and the actual horizontal polarization transmission coefficient) obtained by actual measurement, and obtaining the fit between the actual transmission coefficient and the predicted transmission coefficient.

[0121] The actual reflection coefficient and the actual transmission coefficient are obtained according to the actual measurement, including:

[0122] Orbital angular momentum (OAM) wavefront measurement;

[0123] Calculates wavefront measurement data purity.

[0124] Optionally, a millimeter wave OAM propagation characteristic measurement scheme based on the free space method is used to perform OAM wavefront measurement. Figure 2 As shown in the figure, a frequency domain measurement system platform for millimeter wave OAM propagation characteristics is built according to the measurement principle and scheme. The entire OAM wavefront measurement system mainly includes a vector network analyzer (VNA), a virtual array scanning frame, a software control program, and a probe antenna. The workflow of the automated wavefront measurement platform is as follows: first, set the measurement parameters of the OAM wavefront measurement system and complete the connection of the device; then, align the center positions of the transmitter (Transmitter, Tx) and the receiver (Receiver, Rx), and deploy Rx to the first measurement point; then, use the software control program to coordinate the movement of Rx and the data storage of VNA, and measure the amplitude and phase data of the OAM wave radiated in a certain spatial cross section point by point; finally, use the software to converge the points into a surface to obtain the OAM wavefront data.

[0125] like Figure 3 The OAM propagation characteristic measurement principle process shown is as follows: First, the signal source transmitter (Transmitter, hereinafter referred to as Tx) generates a plane wave signal, and then the closely attached spiral phase plate (Spiral Phase Plate, hereinafter referred to as SPP) modulates the plane wave signal into a vortex wave signal, and then the lens suppresses the divergence of the vortex wave. After the OAM wave propagates through free space, it is reflected or transmitted under the action of the building material with an angle of θi, and finally demodulated and received by the combination of SPP and horn antenna at the receiver (hereinafter referred to as Rx). Figure 3 The metal plate in (a) is the reflection measurement control group. Figure 3 (b) is the material being tested through transmission, Figure 3 (c) is the transmission measurement control group.

[0126] Preferably, as an improvement, after obtaining the wavefront measurement data, the modal purity of the OAM wavefront is calculated. Specifically, the modal purity of the OAM wavefront is calculated using the method shown in the following formula:

[0127]

[0128]

[0129] in,, is Laguerre-Gaussian phase sampling, A l Indicates that the mode classification is a positive real number, l is the mode, is the phase, P l A is the modal purity of the OAM wavefront, where k is the number of modes. k Represents the power of mode k. To calculate the purity, divide the target mode Al by the sum of all modes.

[0130] The calculated OAM wavefront measurement modal purity is to verify the reliability of the OAM wavefront measurement data, thereby ensuring an accurate measurement environment. Then, based on this accurate measurement environment, the accurate incident electric field, reflected electric field, and transmitted electric field (i.e., the actual incident electric field, the actual reflected electric field, and the actual transmitted electric field) can be measured, and then the actual propagation coefficient (the actual propagation coefficient includes the actual reflection coefficient and the actual transmission coefficient) can be calculated based on the accurate reflected electric field and the transmitted electric field. The actual reflection coefficient is the ratio of the actual reflected electric field to the actual incident electric field, and the actual transmission coefficient is the ratio of the actual transmitted electric field to the actual incident electric field.

[0131] The theoretical values ​​(predicted reflection coefficient and predicted transmission coefficient) and measured values ​​(actual reflection coefficient and actual transmission coefficient) are plotted in the form of angle on the horizontal axis and propagation coefficient on the vertical axis, and then the MAE and SD between them are calculated. The MAE and SD evaluation indicators are used to evaluate the degree of fit of the system to the measured values. Some coefficient parameters are optimized, such as inaccurate calculation of the initial beam waist radius and insufficient integral calculation range, and some abnormal values ​​are eliminated and repeated measurements are performed for multiple times to improve the fit.

[0132] The multilayer plate Laguerre reflection model and multilayer plate Laguerre transmission model of the present invention can improve the two key indicators of MAE and SD on the basis of the existing models. Among them, the MAE for OAM mode 1 and OAM mode 2 are all lower than 0.1, the SD for OAM mode 1 can be reduced to 0.05, and the SD for the mode can be reduced to 0.06. Compared with the existing reflection and transmission models, the accuracy can be increased to 50% at most.

[0133] In one embodiment, a propagation coefficient modeling device of orbital angular momentum in a multilayer plate is provided, characterized in that it includes a first construction unit, a determination unit, a second construction unit, a third construction unit, a fourth construction unit, and a fifth construction unit, wherein:

[0134] The first construction unit is used to construct an ITU multilayer board iterative propagation coefficient model;

[0135] A determination unit, used for determining a vertical polarization reflection coefficient, a horizontal polarization reflection coefficient, a vertical polarization transmission coefficient and a horizontal polarization transmission coefficient according to an ITU multilayer iterative propagation coefficient model;

[0136] A second construction unit is used to construct a vertically polarized multi-layer plate Laguerre reflection model according to the vertically polarized reflection coefficient;

[0137] A third construction unit is used to construct a horizontally polarized multilayer plate Laguerre reflection model according to the horizontally polarized reflection coefficient;

[0138] A fourth construction unit is used to construct a vertically polarized multi-layer plate Laguerre transmission model according to the vertically polarized transmission coefficient;

[0139] The fifth construction unit is used to construct a horizontally polarized multilayer plate Laguerre transmission model according to the horizontally polarized transmission coefficient.

[0140] In the description of the present invention, a large number of specific details are described. However, it is understood that embodiments of the present invention can be practiced without these specific details. In some instances, well-known methods, structures and techniques are not shown in detail so as not to obscure the understanding of this description.

[0141] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0142] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and specification of the present invention.

Claims

1. A method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate, characterized in that: include: Construct an iterative propagation coefficient model for ITU multilayer boards; Determining a vertical polarization reflection coefficient, a horizontal polarization reflection coefficient, a vertical polarization transmission coefficient, and a horizontal polarization transmission coefficient according to the ITU multilayer iterative propagation coefficient model; Constructing a vertically polarized multi-layer plate Laguerre reflection model according to the vertically polarized reflection coefficient; Constructing a horizontally polarized multilayer plate Laguerre reflection model according to the horizontally polarized reflection coefficient; Constructing a vertically polarized multilayer plate Laguerre transmission model according to the vertically polarized transmission coefficient; A Laguerre transmission model of a horizontally polarized multilayer plate is constructed according to the horizontally polarized transmission coefficient.

2. A method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate as claimed in claim 1, characterized in that: The ITU multilayer iterative propagation coefficient model is characterized by the following formula: A n =0.5exp(jk n d n cosθ n )[A n+1 (1+Y n+1 )+B n+1 (1-Y n+1 )]; B n =0.5exp(-jk n d n cosθ n )[A n+1 (1-Y n+1 )+B n+1 (1+Y n+1 )]; F n =0.5exp(jk n d n cosθ n )[F n+1 (1+W n+1 )+G n+t (1-W n+1 )]; G n =0.5exp(-jk n d n cosθ n )[F n+1 (1-W n+1 )+G n+1 (1+W n+1 )]; Among them, A n , B n , F n , G n , W n+1 , Y n+1 , k n is the constant term expression of the recursive algorithm, j is a complex unit, θ n is the propagation direction of the nth layer, d n is the width of the multilayer board.

3. A method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate as claimed in claim 2, characterized in that: The vertical polarization reflection coefficient, horizontal polarization reflection coefficient, vertical polarization transmission coefficient and horizontal polarization transmission coefficient are respectively characterized by the following formulas: Among them, R ⊥ is the vertical polarization reflection coefficient, R P is the horizontal polarization reflection coefficient, T ⊥ is the vertical polarization transmission coefficient, T P is the horizontal polarization transmission coefficient, and A0, B0, G0, and F0 represent the four recursive parameters of the recursive algorithm.

4. A method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate as claimed in claim 3, characterized in that: The vertically polarized multilayer plate Laguerre reflection model is constructed based on the vertically polarized reflection coefficient, including: The predicted reflection coefficient of the vertically polarized multilayer plate Laguerre reflection model is calculated based on the vertically polarized reflection coefficient using the method shown in the following formula: Where Γ1 is the predicted reflection coefficient of the vertically polarized multilayer Laguerre reflection model, E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

5. The method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate according to claim 3, characterized in that: The Laguerre reflection model of horizontally polarized multilayer plates is constructed based on the horizontally polarized reflection coefficient, including: The predicted reflection coefficient of the horizontally polarized multilayer plate Laguerre reflection model is calculated based on the horizontally polarized reflection coefficient using the method shown in the following formula: Where Γ2 is the predicted reflection coefficient of the Laguerre reflection model for horizontally polarized multilayer plates, and E r (r r ,φ r ,z r ) is the reflected electric field, r r is the reflection distance, φ r is the reflection angle of electromagnetic wave, z r is the z-axis propagation distance of the reflected electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

6. The method for modeling the propagation coefficient of orbital angular momentum in a multilayer plate according to claim 3, characterized in that: The vertically polarized multilayer plate Laguerre transmission model is constructed based on the vertically polarized transmission coefficient, including: The predicted transmission coefficient of the vertically polarized multilayer Laguerre transmission model is calculated based on the vertically polarized transmission coefficient using the method shown in the following formula: Where T1 is the predicted transmission coefficient of the vertically polarized multilayer Laguerre reflection model, E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ t is the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

7. A method for modeling propagation coefficient of orbital angular momentum in a multilayer plate as claimed in claim 5, characterized in that: The Laguerre transmission model of horizontally polarized multilayer plates is constructed based on the horizontally polarized transmission coefficient, including: The predicted transmission coefficient of the horizontally polarized multilayer plate Laguerre transmission model is calculated based on the horizontally polarized transmission coefficient using the method shown in the following formula: Where T2 is the predicted transmission coefficient of the Laguerre transmission model for horizontally polarized multilayer plates, and E t (r t ,φ t ,z t ) is the transmitted electric field, r t is the straight-line propagation distance of the transmitted electric field, φ t is the electromagnetic wave transmission angle, z t is the z-axis propagation distance of the transmitted electric field, E i (r i ,φ i ,z i ) is the incident electric field, r i is the straight-line propagation distance of the incident electric field, φ i is the incident angle of the electromagnetic wave, z i is the z-axis propagation distance of the incident electric field.

8. A device for modeling the propagation coefficient of orbital angular momentum in a multilayer plate, characterized in that: The method comprises a first building block, a determining block, a second building block, a third building block, a fourth building block, and a fifth building block, wherein: The first construction unit is used to construct an ITU multilayer board iterative propagation coefficient model; The determining unit is used to determine the vertical polarization reflection coefficient, the horizontal polarization reflection coefficient, the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient according to the ITU multilayer iterative propagation coefficient model; The second construction unit is used to construct a vertically polarized multi-layer plate Laguerre reflection model according to the vertically polarized reflection coefficient; The third construction unit is used to construct a Laguerre reflection model of a horizontally polarized multilayer plate according to the horizontally polarized reflection coefficient; The fourth construction unit is used to construct a vertically polarized multilayer plate Laguerre transmission model according to the vertically polarized transmission coefficient; The fifth construction unit is used to construct a horizontally polarized multilayer plate Laguerre transmission model according to the horizontally polarized transmission coefficient.

Images