Method and device for simulating antenna radiation characteristics of multilayer non-parallel medium

Abstract

The application provides a kind of multi-layer non-parallel medium antenna radiation characteristic simulation method and device, wherein the method comprises: obtaining the model information of the grid model of multi-layer non-parallel medium and the initialization information of incident ray;Model information includes the geometric information and material information of triangular grid;The initialization information of ray includes: incident origin, direction vector and vector field;Multi-order ray tracing is carried out on the split rays generated after the incident ray from incident origin to the interface of multi-layer non-parallel medium;The tracking content includes the split point, direction vector and vector field of each order split ray;Based on the tracking content obtained after multi-order ray tracing, the vector field of all split rays with split point located on the outer surface of multi-layer non-parallel medium is coherently superimposed to obtain the outer surface vector field;According to the outer surface vector field, the far-field pattern for simulating antenna radiation characteristics is calculated.The scheme can quickly calculate the antenna radiation characteristics of multi-layer non-parallel medium.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of modeling technology, and in particular to a method and apparatus for simulating the radiation characteristics of antennas in multilayer non-parallel media. Background Technology

[0002] Multilayer non-parallel media are commonly used in radar system structures such as radomes, which can affect the antenna radiation characteristics, causing distortions in the main beam shape, changes in beamwidth, reduced maximum gain, increased sidelobe levels, and increased wide-angle sidelobe levels. In practical engineering, full-wave simulation methods or instrumental measurements are mainly used to calculate the far-field radiation pattern of multilayer non-parallel media. However, full-wave simulation methods are inefficient and memory-intensive, while instrumental measurements are time-consuming and often limited by the complex radar system structure.

[0003] Therefore, there is an urgent need to provide a method that can quickly calculate the radiation characteristics of antennas in multilayer non-parallel media. Summary of the Invention

[0004] This invention provides a method and apparatus for simulating the antenna radiation characteristics of multilayer non-parallel media, which can quickly calculate the antenna radiation characteristics of multilayer non-parallel media.

[0005] In a first aspect, embodiments of the present invention provide a method for simulating the radiation characteristics of an antenna in a multilayer non-parallel medium, comprising:

[0006] The model information of the mesh model of the multilayer non-parallel medium and the initialization information of the incident ray are obtained; the model information includes the geometric information and material information of the triangular mesh; the initialization information of the ray includes: the incident origin, the direction vector, and the vector field;

[0007] Based on the model information of the mesh model and the initialization information of the ray, multi-order ray tracing is performed on the splitting ray generated after the ray enters from the incident origin to the interface of the multi-layer non-parallel medium; the tracing content includes the splitting point, direction vector and vector field of each order splitting ray;

[0008] Based on the tracing content obtained after multi-order ray tracing, the vector fields of all split rays with split points located on the outer surface of the multilayer non-parallel medium are coherently superimposed to obtain the outer surface vector field.

[0009] The far-field radiation pattern used to simulate the antenna radiation characteristics is calculated based on the outer surface vector field.

[0010] Secondly, embodiments of the present invention also provide an antenna radiation characteristic simulation device for multilayer non-parallel media, comprising:

[0011] The acquisition unit is used to acquire model information of a mesh model of a multi-layer non-parallel medium and initialization information of incident rays; the model information includes geometric information and material information of the triangular mesh; the initialization information of the rays includes: the incident origin, direction vector, and vector field;

[0012] The tracking unit is used to perform multi-order ray tracking on the split rays generated after the rays are incident from the origin to the interface of the multi-layer non-parallel medium, based on the model information of the mesh model and the initialization information of the rays; the tracking content includes the splitting point, direction vector and vector field of each order split ray;

[0013] The superposition unit is used to coherently superimpose the vector fields of all split rays whose split points are located on the outer surface of the multilayer non-parallel medium based on the tracking content obtained after multi-order ray tracing, so as to obtain the outer surface vector field.

[0014] The calculation unit is used to calculate the far-field radiation pattern for simulating the radiation characteristics of the antenna based on the vector field of the outer surface.

[0015] Thirdly, embodiments of the present invention also provide an electronic device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, it implements the method described in any embodiment of this specification.

[0016] Fourthly, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the methods described in any embodiment of this specification.

[0017] This invention provides a method and apparatus for simulating the radiation characteristics of antennas in multilayer non-parallel media. By acquiring the model information of the mesh model of the multilayer non-parallel media and the initialization information of the incident rays, multi-order ray tracing can be performed on the split rays generated after the rays are incident from the origin to the interface of the multilayer non-parallel media. Based on the tracing content obtained from the multi-order ray tracing, the vector fields of all split rays whose split points are located on the outer surface of the multilayer non-parallel media are coherently superimposed, thereby quickly obtaining the outer surface vector field. This outer surface vector field is then used to calculate the far-field radiation pattern used to simulate the antenna radiation characteristics. Therefore, this method can quickly calculate the antenna radiation characteristics of multilayer non-parallel media. Attached Figure Description

[0018] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0019] Figure 1 This is a flowchart of a simulation method for the radiation characteristics of an antenna in a multilayer non-parallel medium, provided by an embodiment of the present invention.

[0020] Figure 2 This is a hardware architecture diagram of an electronic device provided in an embodiment of the present invention;

[0021] Figure 3 This is a structural diagram of an antenna radiation characteristic simulation device for multilayer non-parallel media provided in an embodiment of the present invention. Detailed Implementation

[0022] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0023] Please refer to Figure 1 This invention provides a method for simulating the radiation characteristics of an antenna in a multilayer non-parallel medium. The method includes:

[0024] Step 100: Obtain the model information of the mesh model of the multilayer non-parallel medium and the initialization information of the incident ray; the model information includes the geometric information and material information of the triangular mesh; the initialization information of the ray includes: the incident origin, the direction vector, and the vector field;

[0025] Step 102: Based on the model information of the mesh model and the initialization information of the ray, perform multi-order ray tracing on the splitting ray generated after the ray enters from the incident origin to the interface of the multi-layer non-parallel medium; the tracing content includes the splitting point, direction vector and vector field of each order splitting ray;

[0026] Step 104: Based on the tracking content obtained after multi-order ray tracing, the vector fields of all split rays with split points located on the outer surface of the multilayer non-parallel medium are coherently superimposed to obtain the outer surface vector field.

[0027] Step 106: Calculate the far-field radiation pattern for simulating the antenna radiation characteristics based on the outer surface vector field.

[0028] In this embodiment of the invention, by acquiring the model information of the mesh model of the multilayer non-parallel medium and the initialization information of the incident rays, multi-order ray tracing can be performed on the split rays generated after the rays are incident from the origin to the interface of the multilayer non-parallel medium. Based on the tracing content obtained from the multi-order ray tracing, the vector fields of all split rays whose split points are located on the outer surface of the multilayer non-parallel medium are coherently superimposed, thereby quickly obtaining the outer surface vector field. This outer surface vector field is then used to calculate the far-field radiation pattern for simulating the antenna radiation characteristics. Therefore, this scheme can quickly calculate the antenna radiation characteristics of multilayer non-parallel media.

[0029] The following description Figure 1 The execution method for each step is shown.

[0030] First, for step 100, obtain the model information of the mesh model of the multi-layer non-parallel medium and the initialization information of the incident rays.

[0031] In this embodiment of the invention, a mesh model needs to be pre-constructed for a multi-layered non-parallel medium. This mesh model includes multiple mesh layers, and each layer is composed of multiple triangular meshes, with each triangular mesh being a face element. Therefore, the model information of the mesh model includes the geometric information of the triangular meshes, which can be described using triangular face elements. This geometric information mainly includes the vertex position information of the triangle, the indices of the three vertices corresponding to the face element, and the face element normal vector.

[0032] In addition, the medium is a multi-layered non-parallel medium, so the material information on both sides of the triangular mesh on each layer of the mesh is different. Therefore, the model information also includes the material information on both sides of the triangular mesh.

[0033] To determine the antenna radiation characteristics, initialization information of the ray is also required. This initialization information includes: the origin of incidence, the direction vector, and the vector field. The origin of incidence represents the initial position of the ray when it enters the non-parallel medium, the direction vector represents the incident direction of the ray, and the vector field represents the initial electric field along the incident direction of the ray.

[0034] It is understood that the initial incident rays into the non-parallel medium are multiple rays. Each ray has an incident origin, a direction vector, and a vector field. In this embodiment of the invention, the density and number of rays can be flexibly set according to the calculation accuracy and the size of the triangular mesh, so as to minimize the number of rays while ensuring high calculation accuracy.

[0035] Then, for step 102, based on the model information of the mesh model and the initialization information of the ray, multi-order ray tracing is performed on the splitting ray generated after the ray enters from the incident origin to the interface of the multi-layer non-parallel medium.

[0036] In this embodiment of the invention, since the medium is a multilayer non-parallel medium including multiple interfaces, when a ray is incident on an interface, it generates a splitting ray. This splitting ray includes splitting reflected rays and splitting refracted rays. As the splitting ray propagates based on its direction vector, it will also be incident on the next interface, thus generating a splitting ray at the next interface, and so on. In order to determine the vector field of the outer surface of the medium and then calculate the far-field radiation pattern used to simulate the antenna radiation characteristics, it is necessary to perform multi-order ray tracing on the splitting rays until the tracing ends.

[0037] Specifically, multi-level ray tracing can be performed in the following way:

[0038] S1: Based on the model information of the mesh model and the initialization information of the ray, determine the first-order splitting ray generated after the ray enters the multilayer non-parallel medium interface from the incident origin;

[0039] S2: Calculate the tracking content of the currently generated split ray of this order, and determine whether the current tracking end condition has been met. If so, end the tracking; otherwise, based on the tracking content of the split ray of this order, determine the next split ray generated after the split ray of this order is incident on the next interface, and continue to execute this step until the tracking ends.

[0040] In this embodiment of the invention, since the higher the order of the split ray generated by the ray, the more the vector field of the split ray attenuates, the tracking termination condition may include: the order of the currently generated split ray reaches a set tracking threshold, or the vector field of the currently generated split ray attenuates to a set field threshold.

[0041] The tracking threshold can be set based on experience, while the field threshold can be set based on the calculation accuracy.

[0042] In this embodiment of the invention, all rays incident from the origin to the interface of the multilayer non-parallel medium can be considered as the root node of a tree. The splitting ray generated from the splitting at this interface can be called the first-order splitting ray, the splitting ray generated from the further splitting of the first-order splitting ray can be called the second-order splitting ray, and so on, and the splitting ray generated from the further splitting of the (m-1)th-order splitting ray can be called the m-th-order splitting ray. Where m is a positive integer.

[0043] To improve the computational speed of split ray tracing, starting from the root node, split reflected rays generated from the same node are treated as a single reflection node, and split refracted rays generated from the same node are treated as a single refraction node. Rays under each node are computed in parallel. That is, when any ray from the root node is reflected, a reflection branch is generated at a reflection node of the same order, and all split reflected rays are under that reflection node. Similarly, when any ray from the root node is refracted, a refraction branch is generated at a refraction node of the same order, and all split refracted rays are under that refraction node. These split refracted or reflected rays serve as new emission sources for tracing, resulting in a new layer of nodes. This process continues until tracing ends, yielding a complete ray tracing tree of the ray cluster. This allows for parallel computation of all rays under each node, significantly improving the computational speed of ray splitting and tracing.

[0044] For example, if the initial number of rays is 100, these 100 rays serve as the root node of the tree. At the interface, they split, producing 100 first-order split reflection rays and 98 split refraction rays. These 100 first-order split reflection rays form a first-order reflection node, and the 98 first-order split refraction rays form a first-order refraction node. Tracing continues on these first-order reflection and refraction nodes. The 100 first-order split reflection rays further split, producing 100 second-order split reflection rays and 99 second-order split refraction rays. These 100 second-order split reflection rays form a reflection node, and the 99 second-order split refraction rays form a refraction node. Similarly, the 98 first-order split refraction rays further split, producing 98 second-order split reflection rays, 97 second-order split refraction rays, and so on. These 98 second-order split reflection rays form a reflection node, and the 97 second-order split refraction rays form a refraction node. In this way, a complete ray tracing tree for the ray cluster can be obtained.

[0045] Since there are many rays under each node, the ray content under each node can be calculated in parallel to improve the calculation speed.

[0046] In this embodiment of the invention, the tracking content includes: the splitting point, direction vector, and vector field of each order splitting ray. The method for calculating the tracking content of the currently generated splitting ray of that order includes: for each splitting ray of that order, performing the following steps: determining the splitting point on the interface that generated the splitting ray based on the direction vector of the incident ray that generated the splitting ray and the model information of the mesh model; calculating the direction vector of the splitting ray based on the material information on both sides of the interface, the direction vector of the incident ray that generated the splitting ray, and the splitting point; wherein, the splitting point of the first-order splitting ray is the origin of incidence; and calculating the vector field of the splitting ray based on the vector field of the incident ray that generated the splitting ray.

[0047] Specifically, the direction vectors of the split reflected ray and the split refracted ray are calculated according to the following formula:

[0048] r = i - 2n·(i·n)

[0049]

[0050] Where r is the direction vector of the split reflected ray, t is the direction vector of the split refracted ray, i is the direction vector of the incident ray that produces the split ray, n is the outward normal unit vector of the interface that produces the split ray, and k1 and k2 are the propagation constants of the space where the incident ray and the space where the split refracted ray are produced, respectively. If the incident ray undergoes total internal reflection at the interface and there is no splitting refracted ray, it indicates that the incident ray was not refracted into the next medium at the interface.

[0051] The methods for calculating the vector field of split rays include:

[0052] The incident ray that generates the split ray is projected onto the interface in the directions of the vertical polarization vector and the parallel polarization vector, respectively, so as to decompose the vector field of the incident ray into the vertical incident wave component and the horizontal incident wave component; wherein, the formulas for calculating the vertical polarization vector s and the parallel polarization vector p are: s = n × k, p = k × s; k is the propagation constant of the space in which the incident ray is located;

[0053] The vector fields of split reflected rays and split refracted rays are calculated according to the following formula:

[0054]

[0055]

[0056] in, These are the vector fields of the split reflected ray and the split refracted ray, respectively. These are the vertical polarization vector and the parallel polarization vector, θ. i =cos -1 (n·k) is the incident angle of the incident ray, r ∥ r ⊥ These are the Fresnel coefficients and complex amplitude reflection coefficients for parallel and perpendicular polarization, respectively, t ∥ t ⊥ E represents the transmission coefficients for parallel polarization and vertical polarization, respectively. ∥ E ⊥ These are the horizontal incident wave component and the vertical incident wave component, respectively.

[0057] Finally, steps 104 and 106 will be explained simultaneously.

[0058] Since each layer of the multilayer non-parallel medium is composed of triangular surface elements, after multi-order ray tracing, the vector field corresponding to each triangular surface element on the outer surface of the multilayer non-parallel medium can be calculated. Then, the outer surface vector field is used as the excitation source, the far field is calculated according to the surface integral formula, and the far field pattern is calculated based on the far field.

[0059] Based on the equivalence principle, the vector field of the outer surface is calculated, and the far-field radiation field is calculated by integrals covering the outer surface.

[0060]

[0061] In the above formula, J = n × H T M = n × E T These are the equivalent surface current and the equivalent surface magnetic current, respectively. After mathematical calculations and simplification, the output far-field is:

[0062]

[0063] Where l is the free space wavenumber, Z0 is the free space wave impedance, u is the unit vector of dS towards the observation point, r' is the observation point P(x', y', z'), r is the distance from the observation point P to (x, y, z), and n is the unit vector of the outward normal of the medium interface; E T H T Let be the total electric field and total magnetic field on the surface.

[0064] After obtaining the above far field, the far field radiation pattern used to simulate the antenna radiation characteristics can be calculated using the calculation formula of the far field radiation pattern.

[0065] like Figure 2 , Figure 3 As shown, this embodiment of the invention provides a simulation device for the radiation characteristics of an antenna in a multilayer non-parallel medium. The device embodiment can be implemented through software, hardware, or a combination of both. From a hardware perspective, as... Figure 2 The diagram shown is a hardware architecture diagram of an electronic device used in an embodiment of the present invention to simulate the antenna radiation characteristics of a multilayer non-parallel medium. (Except for...) Figure 2 In addition to the processor, memory, network interface, and non-volatile memory shown, the electronic device in the embodiment may also include other hardware, such as a forwarding chip responsible for processing packets. Taking software implementation as an example, such as... Figure 3 As shown, a device in a logical sense is formed by the CPU of its host electronic device reading the corresponding computer program from the non-volatile memory into the main memory for execution. This embodiment provides a simulation device for the antenna radiation characteristics of a multilayer non-parallel medium, comprising:

[0066] The acquisition unit 301 is used to acquire model information of the mesh model of the multilayer non-parallel medium and initialization information of the incident ray; the model information includes the geometric information and material information of the triangular mesh; the initialization information of the ray includes: the incident origin, the direction vector, and the vector field;

[0067] The tracking unit 302 is used to perform multi-order ray tracking on the split rays generated after the rays are incident from the origin to the interface of the multi-layer non-parallel medium, based on the model information of the mesh model and the initialization information of the rays; the tracking content includes the splitting point, direction vector and vector field of each order split ray;

[0068] The superposition unit 303 is used to coherently superimpose the vector fields of all split rays whose split points are located on the outer surface of the multilayer non-parallel medium based on the tracking content obtained after multi-order ray tracing, to obtain the outer surface vector field.

[0069] The calculation unit 304 is used to calculate the far-field radiation pattern for simulating the radiation characteristics of the antenna based on the vector field of the outer surface.

[0070] In one embodiment of the present invention, the tracking unit specifically includes: determining, based on the model information of the mesh model and the initialization information of the ray, the first-order splitting ray generated after the ray enters the multilayer non-parallel medium interface from the incident origin; calculating the tracking content of the currently generated splitting ray and determining whether the current tracking termination condition has been met; if so, ending the tracking; otherwise, based on the tracking content of the splitting ray, determining the next-order splitting ray generated after the splitting ray enters the next interface, and continuing to execute this step until the tracking ends.

[0071] In one embodiment of the present invention, the splitting ray includes splitting reflected rays and splitting refracted rays;

[0072] All rays incident from the origin to the interface of the multilayer non-parallel medium are taken as the root node of the tree. Starting from the root node, the split reflected rays generated by the same node are taken as a reflection node, and the split refracted rays generated by the same node are taken as a refraction node. The rays under each node are calculated and tracked in parallel.

[0073] In one embodiment of the present invention, the tracking termination condition includes: the order of the currently generated split ray reaches a set tracking threshold, or the vector field of the currently generated split ray decays to a set field threshold.

[0074] In one embodiment of the present invention, when the tracking unit calculates the tracking content of the currently generated split ray of that order, it specifically includes:

[0075] For each split ray of this order, the following steps are performed: Based on the direction vector of the incident ray that generated the split ray and the model information of the mesh model, determine the split point on the interface that generated the split ray; based on the material information on both sides of the interface, the direction vector of the incident ray that generated the split ray, and the split point, calculate the direction vector of the split ray; wherein, the split point of the first-order split ray is the origin of incidence; calculate the vector field of the split ray based on the vector field of the incident ray that generated the split ray.

[0076] In one embodiment of the present invention, the splitting ray includes splitting reflected rays and splitting refracted rays;

[0077] The tracking unit, when calculating the direction vector of the split ray, specifically includes:

[0078] The direction vectors of split reflected rays and split refracted rays are calculated according to the following formula:

[0079] r = i - 2n·(i·n)

[0080]

[0081] Where r is the direction vector of the split reflected ray, t is the direction vector of the split refracted ray, i is the direction vector of the incident ray that produces the split ray, n is the outward normal unit vector of the interface that produces the split ray, and k1 and k2 are the propagation constants of the space where the incident ray and the space where the split refracted ray are produced, respectively. The incident ray will undergo total internal reflection at the interface, and there will be no split refracted ray.

[0082] When the tracking unit calculates the vector field of the split ray based on the vector field of the incident ray that generated the split ray, it specifically includes:

[0083] The incident ray that produces the split ray is projected onto the interface in the directions of the vertical polarization vector and the parallel polarization vector, respectively, so as to decompose the vector field of the incident ray into the vertical incident wave component and the horizontal incident wave component.

[0084] The vector fields of split reflected rays and split refracted rays are calculated according to the following formula:

[0085]

[0086]

[0087] in, These are the vector fields of the split reflected ray and the split refracted ray, respectively. These are the vertical polarization vector and the parallel polarization vector, θ. i =cos -1 (n·k) is the incident angle of the incident ray, r ∥ r ⊥ These are the Fresnel coefficients and complex amplitude reflection coefficients for parallel and perpendicular polarization, respectively, t ∥ t ⊥ E represents the transmission coefficients for parallel polarization and vertical polarization, respectively. / / E ⊥ These are the horizontal incident wave component and the vertical incident wave component, respectively.

[0088] It is understood that the structures illustrated in the embodiments of the present invention do not constitute a specific limitation on an antenna radiation characteristic simulation device for multilayer non-parallel media. In other embodiments of the present invention, an antenna radiation characteristic simulation device for multilayer non-parallel media may include more or fewer components than illustrated, or combine some components, split some components, or arrange different components. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.

[0089] The information interaction and execution process between the modules in the above-mentioned device are based on the same concept as the method embodiment of the present invention, and the specific details can be found in the description in the method embodiment of the present invention, and will not be repeated here.

[0090] This invention also provides an electronic device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, it implements a simulation method for antenna radiation characteristics of a multilayer non-parallel medium according to any embodiment of this invention.

[0091] This invention also provides a computer-readable storage medium storing a computer program. When executed by a processor, the computer program causes the processor to perform a simulation method for antenna radiation characteristics of a multilayer non-parallel medium according to any embodiment of this invention.

[0092] Specifically, a system or apparatus equipped with a storage medium may be provided, on which software program code implementing the functions of any of the embodiments described above is stored, and the computer (or CPU or MPU) of the system or apparatus may read and execute the program code stored in the storage medium.

[0093] In this case, the program code read from the storage medium can itself implement the function of any of the above embodiments, and therefore the program code and the storage medium storing the program code constitute part of the present invention.

[0094] Examples of storage media used to provide program code include floppy disks, hard disks, magneto-optical disks, optical disks (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), magnetic tapes, non-volatile memory cards, and ROMs. Alternatively, program code can be downloaded from a server computer via a communication network.

[0095] Furthermore, it should be clear that not only can the program code read by the computer be executed, but also the operating system or other components operating on the computer can be instructed based on the program code to perform some or all of the actual operations, thereby realizing the function of any of the embodiments described above.

[0096] Furthermore, it is understood that the program code read from the storage medium is written to the memory set in the expansion board inserted into the computer or to the memory set in the expansion module connected to the computer. Then, based on the instructions of the program code, the CPU or other components installed on the expansion board or expansion module execute some and all of the actual operations, thereby realizing the function of any of the above embodiments.

[0097] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0098] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as ROM, RAM, magnetic disk, or optical disk.

[0099] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for simulating the radiation characteristics of an antenna in a multilayer non-parallel medium, characterized in that, include: Obtain model information and incident ray initialization information for a mesh model of a multilayer non-parallel medium; The model information includes the geometric and material information of the triangular mesh; The initialization information of the ray includes: the origin of incidence, the direction vector, and the vector field; Based on the model information of the mesh model and the initialization information of the ray, multi-order ray tracing is performed on the splitting ray generated after the ray enters from the incident origin to the interface of the multi-layer non-parallel medium; the tracing content includes the splitting point, direction vector and vector field of each order splitting ray; Based on the tracing content obtained after multi-order ray tracing, the vector fields of all split rays with split points located on the outer surface of the multilayer non-parallel medium are coherently superimposed to obtain the outer surface vector field. The far-field radiation pattern used to simulate the antenna radiation characteristics is calculated based on the outer surface vector field. The step of performing multi-order ray tracing on the splitting ray generated after the ray enters the multi-layer non-parallel medium interface from the origin of incidence includes: determining the first-order splitting ray generated after the ray enters the multi-layer non-parallel medium interface from the origin of incidence based on the model information of the mesh model and the initialization information of the ray; calculating the tracing content of the currently generated splitting ray and determining whether the tracing termination condition has been met. If so, the tracing is terminated; otherwise, based on the tracing content of the splitting ray, the next-order splitting ray generated after the splitting ray enters the next interface is determined, and this step is continued until the tracing is terminated. The split rays include split reflected rays and split refracted rays; all rays incident from the origin of incidence to the interface of the multilayer non-parallel medium are taken as the root node of the tree. Starting from the root node, the split reflected rays generated by the same node are taken as a reflection node, and the split refracted rays generated by the same node are taken as a refraction node. The rays under each node are calculated and tracked in parallel. The tracking termination conditions include: the order of the currently generated split ray reaches a set tracking threshold, or the vector field of the currently generated split ray decays to a set field threshold.

2. The method according to claim 1, characterized in that, The calculation of the tracking content of the currently generated split ray of this order includes: For each split ray of this order, the following steps are performed: Based on the direction vector of the incident ray that generated the split ray and the model information of the mesh model, determine the split point on the interface that generated the split ray; based on the material information on both sides of the interface, the direction vector of the incident ray that generated the split ray, and the split point, calculate the direction vector of the split ray; wherein, the split point of the first-order split ray is the origin of incidence; calculate the vector field of the split ray based on the vector field of the incident ray that generated the split ray.

3. The method according to claim 2, characterized in that, The splitting rays include splitting reflected rays and splitting refracted rays; The calculation of the direction vector of the splitting ray includes: The direction vectors of split reflected rays and split refracted rays are calculated according to the following formula: in, r The direction vector of the split reflected ray. t The direction vector of the split refracted ray. i To generate the direction vector of the incident ray of this split ray, n To generate the outward normal unit vector of the interface of this splitting ray, k 1. k 2 are the propagation constants of the space containing the incident ray that produces the split ray and the space containing the refracted ray that produces the split ray, respectively. If the incident ray undergoes total internal reflection at the interface, there will be no split refracted ray.

4. The method according to claim 2, characterized in that, The calculation of the vector field of the split ray based on the vector field of the incident ray that generates the split ray includes: The incident ray that produces the split ray is projected onto the interface in the directions of the vertical polarization vector and the parallel polarization vector, respectively, so as to decompose the vector field of the incident ray into the vertical incident wave component and the horizontal incident wave component. The vector fields of split reflected rays and split refracted rays are calculated according to the following formula: in, , These are the vector fields of the split reflected ray and the split refracted ray, respectively. , These are the vertical polarization vector and the parallel polarization vector, respectively. θ i = cos -1 ( n · k ( ) represents the angle of incidence of the incident ray. r ∥ , r ⊥ These are the Fresnel coefficients and complex amplitude reflection coefficients for parallel polarization and perpendicular polarization, respectively. t ∥ , t ⊥ These are the transmission coefficients for parallel polarization and vertical polarization, respectively. , These are the horizontal incident wave component and the vertical incident wave component, respectively.

5. A simulation device for antenna radiation characteristics in a multilayer non-parallel medium, characterized in that, A method for simulating the antenna radiation characteristics of a multilayer non-parallel medium as described in any of claims 1-4 includes: The acquisition unit is used to acquire model information of a mesh model of a multi-layer non-parallel medium and initialization information of incident rays; the model information includes geometric information and material information of the triangular mesh; the initialization information of the rays includes: the incident origin, direction vector, and vector field; The tracking unit is used to perform multi-order ray tracking on the split rays generated after the rays are incident from the origin to the interface of the multi-layer non-parallel medium, based on the model information of the mesh model and the initialization information of the rays; the tracking content includes the splitting point, direction vector and vector field of each order split ray; The superposition unit is used to coherently superimpose the vector fields of all split rays whose split points are located on the outer surface of the multilayer non-parallel medium based on the tracking content obtained after multi-order ray tracing, so as to obtain the outer surface vector field. The calculation unit is used to calculate the far-field radiation pattern for simulating the radiation characteristics of the antenna based on the vector field of the outer surface.

6. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, and the processor, when executing the computer program, implements the method as described in any one of claims 1-4.

7. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the method of any one of claims 1-4.

Images