Terahertz wave regulation and control method and system for inhomogeneous medium discrete data and storage medium

Abstract

The invention discloses a terahertz wave regulation and control method and system for inhomogeneous medium discrete data and a storage medium. The method comprises the following steps: establishing an electromagnetic model based on a multi-dimensional discrete data set; calculating the propagation trajectory of the terahertz wave in the medium by adopting a non-uniform medium ray tracing method; calculating transmission change characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium; the transmission change characteristics comprise an attenuation characteristic, a phase shift characteristic and a polarization change characteristic of the terahertz wave in a transmission process; and performing terahertz wave regulation and control according to the attenuation characteristic and / or the phase shift characteristic and / or the polarization change characteristic of the terahertz wave. According to the method, the problem that the current modeling and electromagnetic calculation scheme is difficult to be suitable for electromagnetic transmission process calculation and terahertz wave regulation and control of the high-dimensional complex inhomogeneous medium is solved.

Description

Technical Field

[0001] The present invention belongs to the field of communication technology, and in particular relates to a terahertz wave control method, system and storage medium for discrete data in an inhomogeneous medium. Background Art

[0002] Calculating the electromagnetic transmission characteristics of inhomogeneous media (gases and fluids) is of vital importance in fields such as communications, remote sensing, and navigation. Related issues include electromagnetic transmission in rainy and foggy weather, the Earth's atmosphere, particularly the ionosphere, plasma, and near-space sheaths. Commonly used methods for modeling inhomogeneous media include analytical models, semi-empirical and semi-analytical models, and fluid dynamics. Fluid dynamics not only constructs a detailed three-dimensional model of the medium but also calculates parameter values ​​at any point, meaning the model exhibits good continuity and is therefore widely used. However, in practice, experimental data is typically presented as discrete points, requiring the parameters of any point to be determined before simulation can be used. When modeling massive amounts of discrete data, this calculation method can lead to severe model distortion.

[0003] Currently, one solution for modeling such massive discrete data is to use a layered modeling approach, dividing the one-dimensional data into several layers, replacing gradual variations with uniformity within the layers. While this solution facilitates computation, it struggles to simulate the electromagnetic wave transmission characteristics of real inhomogeneous media. Furthermore, considering that terahertz wave calculations in inhomogeneous media are primarily used to achieve high-energy, highly directional, and long-distance terahertz transmission, targets in the terahertz frequency band often have ultra-large dimensions, which makes the use of classical electromagnetic calculation methods challenging. To address these computational challenges in ultra-large terahertz wave scenarios, the current approach is to approximate the inhomogeneous medium using a certain number of thin layers. The electromagnetic wave behavior of each layer is then solved and accumulated based on reflection laws or boundary conditions to examine the medium's transmission characteristics. However, as the complexity of the medium increases, the computational accuracy of this algorithm decreases significantly, making it difficult to adapt to the analysis of detailed electromagnetic transmission processes.

[0004] In order to solve the problem that current modeling and electromagnetic calculation schemes are difficult to apply to the calculation of electromagnetic transmission process and terahertz wave control in high-dimensional complex inhomogeneous media, a terahertz wave control method, system and storage medium for discrete data of inhomogeneous media are proposed. Summary of the Invention

[0005] The embodiments of the present invention provide a terahertz wave control method, system and storage medium for discrete data in inhomogeneous media, to at least solve the problem that current modeling and electromagnetic calculation schemes are difficult to apply to the calculation of electromagnetic transmission processes and terahertz wave control in high-dimensional complex inhomogeneous media.

[0006] According to one embodiment of the present invention, a terahertz wave control method for discrete data in an inhomogeneous medium is provided, comprising:

[0007] Building electromagnetic models based on multidimensional discrete data sets;

[0008] The propagation trajectory of terahertz waves in the medium is calculated using the inhomogeneous medium ray tracing method;

[0009] Calculating the transmission change characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium; the transmission change characteristics include the attenuation characteristics, phase shift characteristics and polarization change characteristics of the terahertz wave during the transmission process;

[0010] The terahertz wave is regulated according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave.

[0011] In an exemplary embodiment, the step of establishing an electromagnetic model based on a multi-dimensional discrete data set comprises the following steps:

[0012] Adopt the adjacent data query method to obtain the data near the observation point;

[0013] Calculate the required direction parameters based on the data interpolation near the observation point;

[0014] Determine the inner and outer boundaries of the model calculation domain based on the spatial characteristics of the calculation environment;

[0015] Calculate the required medium parameters, refractive index; the dielectric constant requires different calculation methods depending on the type of inhomogeneous medium.

[0016] In an exemplary embodiment, the method of obtaining data near a calculation observation point by using a nearby data query method includes the following steps:

[0017] If the observation point and the computational domain data point basically coincide, the dataset multiple partitioning query method or direct traversal of the data file is used to find one or more data points closest to the observation point as the data near the observation point;

[0018] If the observation point does not coincide with the computational domain data point, a multiple subdivision query method is used to traverse and find the nearest point. Multiple points around the nearest point are selected based on distance and orientation. Based on the positions of multiple points around the nearest point and the judgment of adjacent surfaces, multiple data points around the observation point are obtained as data near the observation point.

[0019] In an exemplary embodiment, the interpolation calculation of the required directional parameters based on the data near the observation point is to perform an interpolation operation on the data near the observation point based on the directional parameter difference function and the dimensional data of the observation point to obtain the required directional parameters; the dimensional data includes one-dimensional data, two-dimensional data and three-dimensional data.

[0020] In an exemplary embodiment, determining the inner boundary and the outer boundary of the model computing domain according to the spatial characteristics of the computing environment comprises the steps of:

[0021] The outer boundary of the model calculation domain is obtained based on the geometric spatial characteristics of the data near the observation point; the outer boundary represents the boundary where the electromagnetic wave leaves the calculation domain;

[0022] The inner boundary of the model calculation domain is determined according to the boundary where transmission and reflection occur in the calculation environment.

[0023] In an exemplary embodiment, the method of calculating the propagation trajectory of a terahertz wave in a medium using an inhomogeneous medium ray tracing method includes:

[0024] Acquire initial conditions; the initial conditions include the spatial position, azimuth and elevation angle of a given ray;

[0025] Determine the range of ray tracing based on domain boundary conditions;

[0026] Calculate the medium refractive index based on the dielectric constant of the current calculation point;

[0027] The Runge-Kutta method is used to calculate the ray equation in the inhomogeneous medium to obtain the electromagnetic trajectory of the terahertz wave, which is the propagation trajectory of the terahertz wave in the medium.

[0028] In an exemplary embodiment, the method of calculating the transmission variation characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium includes the steps of:

[0029] The phase change of the electric wave along the trajectory is calculated according to the terahertz wave transmission trajectory, which is the phase shift characteristic of the terahertz wave during the transmission process;

[0030] The energy change of the radio wave along the trajectory is calculated based on the terahertz wave transmission trajectory, which is the attenuation characteristic of the terahertz wave during the transmission process;

[0031] The polarization change of the radio wave along the trajectory is calculated based on the terahertz wave transmission trajectory, which is the planned change characteristic of the terahertz wave during the transmission process.

[0032] In an exemplary embodiment, the terahertz wave control according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave comprises the steps of:

[0033] Calculating a functional relationship between a control parameter and a change in a characteristic of the terahertz wave based on the attenuation characteristic and / or the phase shift characteristic and / or the polarization change characteristic of the terahertz wave;

[0034] Calculating the required control parameters based on the functional relationship between the control parameters and the changes in the characteristics of the terahertz wave and the target characteristics of the terahertz wave;

[0035] Configure according to the required control parameters to achieve control of terahertz waves.

[0036] According to yet another embodiment of the present invention, a computer-readable storage medium is provided, which stores a computer program for electronic data exchange, wherein the computer program enables a computer to execute the above-mentioned terahertz wave control method for discrete data in an inhomogeneous medium.

[0037] According to another embodiment of the present invention, a terahertz wave control system for discrete data in an inhomogeneous medium is provided, comprising:

[0038] Terahertz wave control processor;

[0039] and one or more programs, wherein the one or more programs are stored in a memory and configured to be executed by the terahertz wave control, the programs causing a computer to execute the above-mentioned terahertz wave control method for discrete data in an inhomogeneous medium.

[0040] The terahertz wave control method, system and storage medium for discrete data in inhomogeneous media of the present invention have the following advantages:

[0041] (1) Compared with the traditional inhomogeneous medium modeling method, the electromagnetic model is established based on multidimensional discrete data sets. It can effectively handle the electromagnetic modeling problems of a large number of discrete data sets, reduce the deviation of the model, improve the model accuracy, and facilitate the subsequent electromagnetic calculation of engineering data.

[0042] (2) The propagation trajectory of terahertz waves in the medium is calculated by using the inhomogeneous medium ray tracing method, and the transmission change characteristics of terahertz waves are calculated based on the propagation trajectory of terahertz waves in the medium. Compared with the traditional calculation method of using thin layers to approximate inhomogeneous media and solving the electromagnetic wave behavior of each layer according to the reflection law or boundary conditions and accumulating them, it can effectively realize the electromagnetic characteristics analysis of the inhomogeneous electromagnetic environment of large and large electromagnetic fields. At the same time, it clarifies the calculation method of the trajectory, phase, polarization and energy changes of electromagnetic waves during transmission, which is convenient for subsequent engineering control.

[0043] (3) Terahertz waves are controlled based on their attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics. Compared with traditional medium transmission characteristic analysis schemes, this method can effectively control the transmission of terahertz waves in large and large non-uniform electromagnetic environments, which is conducive to expanding the multi-scenario application and intelligent application of devices. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 This is a flow chart of a terahertz wave control method for discrete data in an inhomogeneous medium according to an embodiment of the present invention;

[0045] Figure 2 is a flow chart of step S01 of the terahertz wave control method according to an embodiment of the present invention;

[0046] Figure 3This is a flowchart of sub-step S011 of the terahertz wave control method according to an embodiment of the present invention;

[0047] Figure 4 is a schematic diagram of an embodiment of the present invention using a terahertz plane wave as an incident source;

[0048] Figure 5 is a flowchart of sub-step S013 of the terahertz wave control method according to an embodiment of the present invention;

[0049] Figure 6 is a flow chart of step S02 of the terahertz wave control method according to an embodiment of the present invention;

[0050] Figure 7 is a flow chart of step S03 of the terahertz wave control method according to an embodiment of the present invention;

[0051] Figure 8 is a flow chart of step S04 of the terahertz wave control method according to an embodiment of the present invention;

[0052] Figure 9 It is a structural schematic diagram of a terahertz wave control system for discrete data in an inhomogeneous medium according to an embodiment of the present invention. DETAILED DESCRIPTION

[0053] The present invention will be described in detail below with reference to specific embodiments. The following examples will help those skilled in the art further understand the invention, but are not intended to limit the present invention in any form. It should be noted that those skilled in the art may make several changes and improvements without departing from the scope of the present invention. These all fall within the scope of protection of the present invention. The embodiments of the present invention will be described in further detail below with reference to the accompanying drawings. No specific limitations are imposed in this application based on the specific settings of the actual application environment.

[0054] The present invention adopts a ray tracing method in the terahertz and near-terahertz bands, and utilizes high- and low-order interpolation and nearby point query to process discrete data sets and solve ray equations in inhomogeneous media environments. It can efficiently realize complete transmission characteristic analysis such as path, energy, phase and polarization changes during terahertz wave propagation, providing a feasible new path for the electromagnetic characteristic analysis of terahertz waves in complex inhomogeneous media.

[0055] According to one embodiment of the present invention, a terahertz wave control method for discrete data in an inhomogeneous medium is provided. The flow chart is as follows: Figure 1 Shown, including:

[0056] Step S01: establishing an electromagnetic model based on a multidimensional discrete data set;

[0057] Step S02: Calculate the propagation trajectory of the terahertz wave in the medium using the inhomogeneous medium ray tracing method;

[0058] Step S03: calculating the transmission change characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium;

[0059] Step S04: regulating the terahertz wave according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave.

[0060] The transmission change characteristics include attenuation characteristics, phase shift characteristics and polarization change characteristics of the terahertz wave during the transmission process.

[0061] In an exemplary embodiment, the step S01 is to establish an electromagnetic model based on a multi-dimensional discrete data set, as shown in the flowchart. Figure 2 As shown, the steps include:

[0062] Step S011: using a nearby data query method to obtain data near the observation point;

[0063] Step S012: Calculate the required direction parameters based on the data near the observation point through interpolation;

[0064] Step S013: determining the inner and outer boundaries of the model calculation domain according to the spatial characteristics of the calculation environment;

[0065] Step S014: Calculate the required medium parameters: refractive index; the dielectric constant needs to be calculated in different ways depending on the type of inhomogeneous medium.

[0066] In this embodiment, the cylindrical area data to be used for electromagnetic calculation is given, and the cylinder radius is a; the dielectric area data points are represented by the dielectric constant Determined by, where; in the data domain, the dielectric constant of the medium is determined by the formula ε r ( ρ i ) = ( 1 + ε r 0 ) 2 ) [ 1 − ( 1 − ε r 0 1 + ε r 0 ) ) cos( πρ i a ) ] Calculation, ρ is the cylindrical coordinate radius vector in the cylindrical coordinate system (ρ, φ, z), ρ i is the cylindrical coordinate radius vector value of data point i, a is the radius of the medium cylinder, which is 50 times the wavelength; the relevant parameters of the observation point are calculated using data query and data interpolation methods; and the electromagnetic model of the data is established in combination with the electromagnetic parameters of the medium.

[0067] In an exemplary embodiment, the sub-step S011 uses a nearby data query method to obtain data near the observation point, as shown in the flowchart. Figure 3 As shown, the steps include:

[0068] Step S0111: If the observation point and the computational domain data point are substantially coincident, then a data set multi-partition query method or a direct traversal of the data file is used to find one or more data points closest to the observation point as the data near the observation point;

[0069] Step S0112: If the observation point does not coincide with the computational domain data point, a multiple subdivision query method is used to traverse and find the nearest point. Multiple points around the nearest point are selected based on distance and orientation. Based on the positions of the multiple points around the nearest point and the judgment of the adjacent faces, multiple data points around the observation point are obtained as data near the observation point.

[0070] In this embodiment, data near the calculated observation point is obtained through adjacent data query, wherein one-dimensional data modeling searches for data near at least 3 observation points, two-dimensional data modeling searches for data near at least 9 observation points, and three-dimensional data modeling needs to combine adjacent surface judgment to search for data near at least 27 observation points.

[0071] The data near the calculation point can be obtained through nearby data query. This query is divided into two situations.

[0072] Case 1: The observation point P basically coincides with the computational domain data point.

[0073] Use the multiple partitioning method of the data set or directly traverse the data file to quickly find the data point Q closest to the observation point P. Data point Q is used as the data near the observation point, and then point Q is used directly for subsequent parameter calculations.

[0074] like Figure 4 Figure 2 shows a schematic diagram of using a terahertz plane wave as the incident source. In this case, the observation point P does not coincide with the data point.

[0075] That is, case 2: the observation point P does not coincide with the data point.

[0076] In this case, a multiple-partition query method is used to traverse and find the nearest point Q; 8 points around the nearest point Q (i.e., the 8 points around point Q that are closest to point Q) are selected based on distance and orientation; the nearest point Q and its 8 surrounding points constitute the 9 observation point vicinity data required for two-dimensional data modeling; 27 data points around the observation point P are obtained by judging the adjacent faces, which serve as the 27 observation point vicinity data required for three-dimensional data modeling.

[0077] In an exemplary embodiment, the interpolation calculation of the required directional parameters based on the data near the observation point is to perform an interpolation operation on the data near the observation point based on the directional parameter difference function and the dimensional data of the observation point to obtain the required directional parameters; the dimensional data includes one-dimensional data, two-dimensional data and three-dimensional data.

[0078] In this embodiment, based on the data found near the observation point, one-dimensional data modeling directly interpolates the found data points to calculate the required parameters, and two-dimensional data and three-dimensional data modeling are based on the found data, and the required parameters are interpolated in the corresponding two and three dimensions respectively.

[0079] The parameter interpolation function in the x-axis direction is , where x is the data point to be calculated, x i or x j is the queried data point, and f is the direction parameter value of the discrete data coordinate point; the y-axis direction parameter value can be calculated in the same way.

[0080] In an exemplary embodiment, the sub-step S013 is to determine the inner and outer boundaries of the model calculation domain according to the spatial characteristics of the calculation environment. The flowchart is as follows: Figure 5 As shown, the steps include:

[0081] Step S0131: obtaining the outer boundary of the model calculation domain based on the geometric spatial characteristics of the data near the observation point; the outer boundary represents the boundary where the electromagnetic wave leaves the calculation domain;

[0082] Step S0132: Determine the inner boundary of the model calculation domain according to the occurrence boundary of transmission and reflection in the calculation environment.

[0083] In this embodiment, the inner and outer boundaries of the model's computational domain are defined based on the spatial characteristics of the computational environment. The outer boundary is defined as an open domain, indicating the boundary where electromagnetic waves leave the computational domain. In a cylindrical coordinate system (ρ, φ, z), the outer boundary is determined by ρ = a, where a is the boundary constant (i.e., the radius of the dielectric cylinder). The ray tracing calculation ends when the ρ component of the electromagnetic trajectory value r exceeds a.

[0084] The inner boundary is defined as a boundary where transmission and reflection occur, and is usually determined in a preset manner. In the calculation domain of this embodiment, there is no inner boundary.

[0085] In step S014, the required medium parameter refractive index is calculated. The required medium parameter refractive index is given by Decision, among which It should be noted that the dielectric constant requires different calculation methods depending on the type of inhomogeneous medium; the Debye model can be used for cloud-related models, and the Drew model can be used for plasma-related models.

[0086] In an exemplary embodiment, the step S02 uses the inhomogeneous medium ray tracing method to calculate the propagation trajectory of the terahertz wave in the medium, as shown in the flowchart. Figure 6 Shown, including:

[0087] Step S021, obtaining initial conditions; the initial conditions include the spatial position, azimuth angle and pitch angle of a given ray;

[0088] Step S022: determining the range of ray tracing according to the domain boundary conditions;

[0089] Step S023, calculating the medium refractive index according to the dielectric constant of the current calculation point;

[0090] Step S024: Use the Runge-Kutta method to calculate the ray equation in the inhomogeneous medium to obtain the electromagnetic trajectory of the terahertz wave, that is, the propagation trajectory of the terahertz wave in the medium.

[0091] In this embodiment, the initial conditions are obtained, which include the given ray spatial position r0, azimuth angle and pitch angle , where the initial condition is the given ray azimuth =0° and =90°, x=-50, which means a plane wave is incident on a cylinder with a radius of a;

[0092] The range of ray tracing is determined according to the domain boundary conditions described in the above embodiment. In this embodiment, the domain boundary condition is ρ=a;

[0093] The medium refractive index is calculated based on the dielectric constant of the current calculation point. Here, the medium refractive index is given by Interpolation calculation;

[0094] The Runge-Kutta method is used to numerically calculate the ray equation in the inhomogeneous medium and obtain the electromagnetic trajectory of the terahertz wave. The ray equation is:

[0095] (1)

[0096] The terahertz wave vector is (k r , k θ , k φ ), p' represents the integral path of electromagnetic rays, r is the electromagnetic trajectory, , c is the speed of light, ω is the frequency of the radio wave; the Runge-Kutta method is used to directly calculate the ray equation to obtain the propagation trajectory of the terahertz wave in the medium.

[0097] In an exemplary embodiment, the step S03 is to calculate the transmission variation characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium. The flow chart is as follows: Figure 7 As shown, the steps include:

[0098] Step S031: Calculate the phase change of the electric wave along the trajectory according to the terahertz wave transmission trajectory, that is, the phase shift characteristics of the terahertz wave during the transmission process;

[0099] Step S032: Calculate the energy change of the radio wave along the trajectory according to the terahertz wave transmission trajectory, that is, the attenuation characteristics of the terahertz wave during the transmission process;

[0100] Step S033: Calculate the polarization change of the electric wave along the trajectory according to the terahertz wave transmission trajectory, that is, the planned change characteristics of the terahertz wave during the transmission process.

[0101] In this embodiment, a THz wave transmission trajectory is determined, and based on the determined transmission trajectory, the phase change of the radio wave along the trajectory is calculated. The calculation formula is:

[0102] (2)

[0103] Among them, P0 and P1 are the starting and ending points of the ray, n is the refractive index along the trajectory, k0 is the wave vector, is the frequency of radio waves, and c is the speed of light.

[0104] Based on the transmission trajectory, the energy change of the radio wave along the trajectory is calculated using the following formula:

[0105] (3)

[0106] Among them, E0 is the electric field value at the initial point, and E is the electric field value at the end point.

[0107] Based on the transmission trajectory, the polarization change of the radio wave along the trajectory is calculated using the following formula:

[0108] (4)

[0109] in, is the normal vector along the trajectory, is the secondary normal vector along the trajectory, q0 is the polarization intensity at the starting point P0, and and The ratio of the directional components, q is the torque integral along the trajectory.

[0110] In an exemplary embodiment, the step S04 is to control the terahertz wave according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave, and the flow chart is as follows: Figure 8 As shown, the steps include:

[0111] Step S041: calculating a functional relationship between a control parameter and a change in a characteristic of the terahertz wave according to the attenuation characteristic and / or the phase shift characteristic and / or the polarization change characteristic of the terahertz wave;

[0112] Step S042: calculating the required control parameters according to the functional relationship between the control parameters and the change of the terahertz wave characteristics and the target characteristics of the terahertz wave;

[0113] Step S043: configure according to the required control parameters to achieve control of the terahertz wave.

[0114] In this embodiment, the functional relationship between the control parameters and the change of the terahertz wave characteristics is obtained by training the change relationship between different control parameters and the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave, that is, F( and / or Att and / or e) = g(m), where m represents a control parameter;

[0115] Obtain the terahertz wave target characteristics required for the current application scenario, and obtain F( and / or Att and / or e), according to the function F( and / or Att and / or e)=g(m) to calculate the control parameter m, which is the required control parameter;

[0116] Configure according to the required control parameters to achieve control of terahertz waves.

[0117] According to another embodiment of the present invention, a computer-readable storage medium is provided, which stores a computer program for electronic data exchange, wherein the computer program enables a computer to execute the terahertz wave control method for discrete data in an inhomogeneous medium according to any one of the above embodiments.

[0118] According to another embodiment of the present invention, a terahertz wave control system for discrete data in an inhomogeneous medium is provided. The structural diagram is shown in FIG. Figure 9 Shown, including:

[0119] Terahertz wave control processor;

[0120] Memory;

[0121] and one or more programs, wherein the one or more programs are stored in a memory and configured to be executed by the terahertz wave control, and the programs enable a computer to execute the terahertz wave control method for discrete data in an inhomogeneous medium according to any one of the above embodiments.

[0122] Furthermore, in addition to the above-mentioned embodiments, the algorithm in this application is applicable to any non-uniform medium with macroscopic local continuity, such as clouds, rain, fog and the atmosphere; they are not listed one by one in the application. It should be noted that the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method and environment including a series of elements include not only those elements, but also other elements that are not explicitly listed, or also include elements that are inherent to such process, method and environment. In the absence of further restrictions, the statement only lists typical elements and does not exclude the existence of other identical elements in the process, method and environment including the elements.

[0123] Of course, those skilled in the art should realize that the above embodiments are only used to illustrate the present invention and are not intended to limit the present invention. As long as they are within the scope of the present invention, any changes or modifications to the above embodiments will fall within the scope of protection of the present invention.

Claims

1. A terahertz wave control method for discrete data in an inhomogeneous medium, characterized in that: include: Building electromagnetic models based on multidimensional discrete data sets; The propagation trajectory of terahertz waves in the medium is calculated using the inhomogeneous medium ray tracing method; Calculating the transmission change characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium; the transmission change characteristics include the attenuation characteristics, phase shift characteristics and polarization change characteristics of the terahertz wave during the transmission process; The terahertz wave is regulated according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave.

2. The terahertz wave control method for discrete data in inhomogeneous media according to claim 1, characterized in that: The method of establishing an electromagnetic model based on a multidimensional discrete data set comprises the following steps: Adopt the adjacent data query method to obtain the data near the observation point; Calculate the required direction parameters based on the data interpolation near the observation point; Determine the inner and outer boundaries of the model calculation domain based on the spatial characteristics of the calculation environment; Calculate the required medium parameters, refractive index; the dielectric constant requires different calculation methods depending on the type of inhomogeneous medium.

3. The terahertz wave control method for discrete data in inhomogeneous media according to claim 2, characterized in that: The method of using the adjacent data query method to obtain the data near the observation point includes the following steps: If the observation point and the computational domain data point basically coincide, the dataset multiple partitioning query method or direct traversal of the data file is used to find one or more data points closest to the observation point as the data near the observation point; If the observation point does not coincide with the computational domain data point, a multiple subdivision query method is used to traverse and find the nearest point. Multiple points around the nearest point are selected based on distance and orientation. Based on the positions of multiple points around the nearest point and the judgment of adjacent surfaces, multiple data points around the observation point are obtained as data near the observation point.

4. The terahertz wave control method for discrete data in inhomogeneous media according to claim 2, characterized in that: The interpolation calculation of the required directional parameters based on the data near the observation point is to perform an interpolation operation on the data near the observation point based on the directional parameter difference function and the dimensional data of the observation point to obtain the required directional parameters; the dimensional data includes one-dimensional data, two-dimensional data and three-dimensional data.

5. The terahertz wave control method for discrete data in inhomogeneous media according to claim 2, characterized in that: The method of determining the inner boundary and the outer boundary of the model calculation domain according to the spatial characteristics of the calculation environment includes the following steps: The outer boundary of the model calculation domain is obtained based on the geometric spatial characteristics of the data near the observation point; the outer boundary represents the boundary where the electromagnetic wave leaves the calculation domain; The inner boundary of the model calculation domain is determined according to the boundary where transmission and reflection occur in the calculation environment.

6. The terahertz wave control method for discrete data in inhomogeneous media according to claim 1, characterized in that: The method of calculating the propagation trajectory of the terahertz wave in the medium using the inhomogeneous medium ray tracing method includes: Acquire initial conditions; the initial conditions include the spatial position, azimuth and elevation angle of a given ray; Determine the range of ray tracing based on domain boundary conditions; Calculate the medium refractive index based on the dielectric constant of the current calculation point; The Runge-Kutta method is used to calculate the ray equation in the inhomogeneous medium to obtain the electromagnetic trajectory of the terahertz wave, which is the propagation trajectory of the terahertz wave in the medium.

7. The terahertz wave control method for discrete data in inhomogeneous media according to claim 1, characterized in that: The method of calculating the transmission change characteristics of the terahertz wave according to the propagation trajectory of the terahertz wave in the medium includes the following steps: The phase change of the electric wave along the trajectory is calculated according to the terahertz wave transmission trajectory, which is the phase shift characteristic of the terahertz wave during the transmission process; The energy change of the radio wave along the trajectory is calculated based on the terahertz wave transmission trajectory, which is the attenuation characteristic of the terahertz wave during the transmission process; The polarization change of the radio wave along the trajectory is calculated based on the terahertz wave transmission trajectory, which is the planned change characteristic of the terahertz wave during the transmission process.

8. The terahertz wave control method for discrete data in inhomogeneous media according to claim 1, characterized in that: The method of regulating the terahertz wave according to the attenuation characteristics and / or phase shift characteristics and / or polarization change characteristics of the terahertz wave comprises the following steps: Calculating a functional relationship between a control parameter and a change in a characteristic of the terahertz wave based on the attenuation characteristic and / or the phase shift characteristic and / or the polarization change characteristic of the terahertz wave; Calculating the required control parameters based on the functional relationship between the control parameters and the changes in the characteristics of the terahertz wave and the target characteristics of the terahertz wave; Configure according to the required control parameters to achieve control of terahertz waves.

9. A computer-readable storage medium storing a computer program for electronic data exchange, wherein: The computer program enables a computer to execute the method according to any one of claims 1 to 8.

10. A terahertz wave control system for discrete data in an inhomogeneous medium, characterized in that include: Terahertz wave control processor; Memory; as well as One or more programs, wherein the one or more programs are stored in a memory and configured to be executed by the terahertz wave control, and the programs enable a computer to execute the method according to any one of claims 1 to 8.

Images