Electromagnetic scattering simulation method, device, equipment, medium and computer program product

Abstract

The invention relates to an electromagnetic scattering simulation method, device and equipment, a medium and a computer program product. The method comprises the following steps: acquiring a target base station of an area where target intelligent sensing equipment is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station; dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing a scatterer electromagnetic scattering field simulation model according to the plurality of triangular large surface elements; according to the scatterer electromagnetic scattering field simulation model, determining the induction current density of the surface of the scatterer through a large-surface-element physical optical method; and determining an electromagnetic scattering field of the scatterer to the target intelligent sensing equipment according to the induced current density. The method can give consideration to both the calculation speed and the calculation amount, and improves the calculation efficiency.

Description

Technical Field

[0001] The present application relates to the field of electromagnetic simulation technology, and in particular to an electromagnetic scattering simulation method, apparatus, device, medium and computer program product. Background Art

[0002] Substations are the core of the power grid. Fully automated control and management of substation equipment and other power equipment, including online real-time monitoring, dispatching, video surveillance of on-site operations, and outdoor facility status, as well as timely response and rapid resolution of potential abnormal disturbances, rely on the widespread deployment of intelligent sensing devices within the substation. Deploying 5G mobile communication technology networks in substations complicates the electromagnetic environment and wireless channels of intelligent devices due to the high transmit power of 5G mobile communication base station antennas and the shielding and scattering of 5G mobile communication intermediate frequency signals by metal scatterers within the substation. This makes the intelligent sensing devices highly susceptible to external electromagnetic interference. Achieving electromagnetic compatibility requires specific determination of electromagnetic radiation field strength at different locations.

[0003] Currently, electromagnetic scattering field calculations are generally performed using the method of moments or the physical optics method. The method of moments offers excellent computational stability and high accuracy, but its high accuracy presupposes a highly refined mesh in the simulation model. The finer the mesh, the higher the accuracy. However, this places a significant burden on computer memory and computational resources, making it impossible to calculate electromagnetic scattering fields within substations using current general-purpose computers alone. Traditional physical optics methods ignore the induced currents in shadowed areas when calculating the target surface's induced currents. While this greatly simplifies the calculation process, speeds up the computation, and reduces computational memory requirements, they make numerous assumptions, resulting in low accuracy and, consequently, only an approximate representation of the scattered field. Further improving the accuracy and efficiency of the physical optics method in ultra-high frequency electromagnetic scattering studies requires a finer meshing of the target model, with higher frequencies requiring finer meshing. However, this also leads to excessive computational effort and slow computational speed.

[0004] Therefore, how to provide an electromagnetic scattering simulation method that takes into account both computing speed and computational complexity is a technical problem that needs to be solved urgently. Summary of the Invention

[0005] Based on this, it is necessary to provide an electromagnetic scattering simulation method, device, equipment, medium and computer program product that can take into account both calculation speed and calculation amount to address the above technical problems.

[0006] In a first aspect, the present application provides an electromagnetic scattering simulation method. The method comprises:

[0007] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0008] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0009] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0010] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0011] In one embodiment, the determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer includes:

[0012] For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established;

[0013] The induced current density on the surface of the scatterer is determined according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

[0014] In one embodiment, determining the induced current density on the surface of the scatterer according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements includes:

[0015] The induced current density on the surface of the scatterer is determined based on the adjacent triangular large-surface elements, the total number of common edges of all the triangular large-surface elements, the areas of the adjacent triangular large-surface elements, the lengths of the common edges of the adjacent triangular large-surface elements, the first position vectors of the adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edges, using the triangular large-surface element physical optics basis function model.

[0016] In one embodiment, constructing a scatterer electromagnetic scattering field simulation model based on the plurality of triangular large surface elements includes:

[0017] According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

[0018] In one embodiment, the preset simulation conditions include:

[0019] The scatterer is an ideal pure conductor, and the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0020] In one embodiment, the method further comprises:

[0021] Acquire radiation field data of the target base station at the location of the target intelligent sensing device;

[0022] Vector superposition is performed based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0023] In a second aspect, the present application further provides an electromagnetic scattering simulation device. The device comprises:

[0024] An acquisition module is used to acquire a target base station in an area where a target intelligent sensing device is located and scatterers in a radiation field formed by electromagnetic waves emitted by the target base station;

[0025] A simulation module, configured to divide the surface of the scatterer into a plurality of triangular large surface elements, and construct an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0026] A first calculation module is used to determine the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0027] The second calculation module is used to determine the electromagnetic scattering field of the scatterer to the target intelligent sensing device according to the induced current density.

[0028] In a third aspect, the present application further provides a computer device. The computer device includes a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are performed:

[0029] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0030] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0031] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0032] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0033] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the following steps:

[0034] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0035] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0036] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0037] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0038] In a fifth aspect, the present application further provides a computer program product. The computer program product includes a computer program that, when executed by a processor, implements the following steps:

[0039] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0040] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0041] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0042] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0043] The embodiments of the present application have the following beneficial effects:

[0044] The electromagnetic scattering simulation method, apparatus, equipment, medium, and computer program product provided in the embodiments of the present application can reduce the number of computational grids and reduce computational complexity by triangulating large surface elements. At the same time, they can more accurately describe the propagation characteristics of high-frequency electromagnetic waves by combining large surface element physical optics methods, taking into account both computational speed and computational accuracy, simplifying computational complexity while ensuring computational accuracy, and improving computational efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 1 is a flow chart of an electromagnetic scattering simulation method according to an embodiment;

[0046] Figure 2 Schematic diagram of the structure of a scatterer electromagnetic scattering field simulation model in one embodiment;

[0047] Figure 3 Schematic diagram of the division of a triangular large surface element in one embodiment;

[0048] Figure 4 Schematic diagram showing a comparison between a traditional physical optics method and a large-surface physical optics method in one embodiment;

[0049] Figure 5 FIG. 4 is a structural block diagram of an electromagnetic scattering simulation device in one embodiment. DETAILED DESCRIPTION

[0050] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.

[0051] The electromagnetic scattering simulation method provided in the embodiment of the present application can be applied to a terminal or a server. The terminal communicates with the server through a network. The data storage system can store data that the server needs to process. The data storage system can be integrated on the server, or it can be placed on the cloud or other network servers. The terminal can be, but is not limited to, various personal computers, laptops, smart phones, tablets, Internet of Things devices and portable wearable devices. The Internet of Things devices can be smart speakers, smart TVs, smart air conditioners, smart car-mounted devices, etc. Portable wearable devices can be smart watches, smart bracelets, head-mounted devices, etc. The server can be implemented as an independent server or a server cluster consisting of multiple servers.

[0052] Example 1

[0053] In one embodiment, Figure 1 As shown, an electromagnetic scattering simulation method is provided, comprising:

[0054] S1. Acquire a target base station in the area where the target intelligent sensing device is located and a scatterer in the radiation field formed by the electromagnetic waves emitted by the target base station;

[0055] S2. Divide the surface of the scatterer into a number of triangular large surface elements, and construct an electromagnetic scattering field simulation model of the scatterer based on the triangular large surface elements;

[0056] S3. Determine the induced current density on the surface of the scatterer by using the large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0057] S4. Determine the electromagnetic scattering field of the scatterer to the target intelligent sensing device based on the induced current density.

[0058] Specifically, the target base station includes a fourth-generation mobile communication technology base station, a fifth-generation mobile communication technology base station or other base stations that emit directional radio electromagnetic waves to form a radiation field. The target intelligent sensing device includes an intelligent sensing device that is affected by the scattering of scatterers in the radiation field formed by the electromagnetic waves emitted by the target base station, and may include secondary equipment. The secondary equipment includes electrical devices used to monitor, control, and protect primary equipment (such as transformers, circuit breakers, etc.) in the power system. Scatterers include electrically large-sized metal equipment. Among them, when the linear dimensions of the equipment (such as length, width, and height) are much larger than the wavelength of the electromagnetic wave, the object is called an electrically large-sized device. In this case, the electromagnetic wave will undergo multiple reflections and interferences on its surface, resulting in complex electromagnetic phenomena.

[0059] For example, first obtain the target base station in the area where the target intelligent sensing device is located and the scatterers in the radiation field formed by the electromagnetic waves emitted by the target base station, and refer to Figure 2 A rectangular coordinate system is established with the center point of the target base station antenna reflector as the coordinate origin. The scatterer surface is divided into several large triangular elements. Based on these large triangular elements, a simulation model of the scatterer's electromagnetic scattering field is constructed. Based on this surface model, the induced current density on the scatterer surface is determined using the large-element physical optics method. Finally, the electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined based on the induced current density. The method of moments offers excellent computational stability and high accuracy, but its high accuracy requires a high degree of mesh refinement in the simulation model. The finer the model mesh, the higher the accuracy. However, this also places a significant burden on computer memory and computational resources. Current general-purpose computers alone cannot complete the calculation of the electromagnetic scattering field in the substation. Compared to the method of moments, by dividing the scatterer surface into several large triangular elements, large-element simulation can significantly reduce the number of meshes and computational effort. By adopting this technical solution, the number of computational grids can be reduced through triangulated large-surface division, reducing computational complexity. At the same time, the large-surface physical optics method can be combined to more accurately describe the propagation characteristics of high-frequency electromagnetic waves, taking into account both computational speed and computational accuracy, simplifying computational complexity while ensuring computational accuracy, and improving computational efficiency.

[0060] In one embodiment, S2 includes:

[0061] According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

[0062] For example, according to the vector wave equation, the electromagnetic scattering field generated by the scatterer in the target base station at the target intelligent sensing device point P can be expressed as an integral equation:

[0063]

[0064] Where A (r) is the magnetic vector potential function, is the magnetic calibration function, ω is the angular frequency of the electromagnetic wave emitted by the target base station antenna, and j represents the imaginary unit.

[0065] By combining Green's function and electromagnetic wave number with the induced current charge density of the scatterer, the magnetic vector potential function A at point P can be obtained: (r) By substituting the magnetic vector potential function, we can get the electric field dyadic Green's function G e (r, r′). Considering the skin effect of the scatterer housing, electromagnetic scattering is mainly generated by a very thin layer of conductor on the surface of the device housing. Therefore, only the induced current distribution on the housing surface needs to be considered. Based on the current continuity equation, a scatterer surface model composed of several triangular large surface elements is combined with preset simulation conditions to construct a scatterer electromagnetic scattering field simulation model, and the electric field integral equation is obtained:

[0066]

[0067] Where ω is the angular frequency of the electromagnetic wave emitted by the target base station antenna, j represents the imaginary unit, μ is the vacuum magnetic permeability, r and r′ represent the distance from point P and any unit surface element dS on the surface of the scatterer shell to the origin of the coordinate system O, respectively, g(r,r′) is the Green's function, J(r′) is the induced current density generated on the surface of the scatterer under the electromagnetic wave radiation of the target base station antenna, and ε is the dielectric constant.

[0068] By adopting this technical solution, the electric field dyadic Green's function and the current continuity equation can be combined with a scatterer surface model composed of several triangular large surface elements for calculation. The scatterer surface model can be converted into a model divided into several triangular large surface elements. In addition, the skin effect of the scatterer shell needs to be considered, and only the induced current distribution on the shell surface needs to be considered, which simplifies the calculation of the model and simplifies the calculation of the complex electromagnetic scattering field simulation model to the calculation of the electric field integral equation, thereby improving the calculation efficiency.

[0069] In one embodiment, the preset simulation conditions include:

[0070] The scatterer is an ideal pure conductor and its size is much larger than the wavelength of the electromagnetic wave emitted by the target base station. The surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0071] Specifically, based on the simulation conditions that the scatterer is an ideal pure conductor, the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave, the electric field integral equation can be further simplified as follows:

[0072]

[0073] Where ω is the angular frequency of the electromagnetic wave emitted by the target base station antenna, j represents the imaginary unit, μ0 is the vacuum permeability, r and r′ represent the distance from point P and any unit surface element dS on the surface of the scatterer shell to the origin of the coordinate system O, R′ represents the distance between the unit surface element dS and point P, n is the external normal vector of the unit surface element dS, and H i is the magnetic field of the incident electromagnetic wave emitted by the target base station antenna, k and k′ are the electromagnetic wave number and the corresponding propagation vector, respectively.

[0074] For example, the rotation vector n×H i It can be replaced by the induced current J(r′), that is:

[0075] J (r′) =2η i ·n×H i (4)

[0076] Where J(r′) is the induced current density generated on the surface of the scatterer under the electromagnetic wave radiation of the target base station antenna, η i is the coefficient that takes into account the influence of the shielding effect. If point P is within the illumination area of the electromagnetic wave of the target base station antenna, then η i = ±1, the illumination area is the area that can be directly illuminated by the incident wave of the target base station antenna, and its positive or negative sign is determined by the relationship between the incident angle of the electromagnetic wave of the target base station antenna and the external normal vector n of the unit surface element dS; if point P is in the shadow area of the electromagnetic wave of the target base station antenna, then η i =0, the shadow area is the area that the incident wave of the target base station antenna cannot directly illuminate, which further yields:

[0077]

[0078] By adopting such a technical solution, simulation can be performed based on the preset assumptions that the scatterer is an ideal pure conductor, the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave. Assuming that the scatterer is an ideal pure conductor and the size of the scatterer is much larger than the wavelength of the incident electromagnetic wave can greatly reduce the number of grid divisions and thus greatly reduce the amount of calculation. In addition, since the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave, ignoring the current in the shadow area significantly reduces the amount of calculation. At the same time, the calculation accuracy is maintained by phase correction of the illuminated area. This can take into account the calculation accuracy while significantly reducing the amount of calculation, and optimize the construction of the scatterer electromagnetic scattering field simulation model.

[0079] In one embodiment, S3 includes:

[0080] S31. For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established;

[0081] S32. Determine the induced current density on the surface of the scatterer according to the triangular large-surface physical optics method basis function model and the total number of common edges of all triangular large-surface elements.

[0082] For example, refer to Figure 3 , we can calculate the value of A for any pair of adjacent triangular large surface elements n + and A n - , establish the basis function model of triangular large surface element physical optics method:

[0083]

[0084] Where, T n + 、T n - and l n Represent any pair of adjacent triangular large surface elements A n + 、A n - The area and length of the common side; ρ n + and ρ n - Respectively represent any pair of adjacent triangular large surface elements by non-common vertices O + , O - The position vector indicated by ρ nc + and ρ nc - Respectively represent any pair of adjacent triangular large surface elements by non-common vertices O + , O -The position vector to the midpoint of the common edge, k n is the phase correction factor.

[0085] The induced current density J(r′) is represented by the triangular large-surface physical optics basis function model and converted into the unknown γ n The solution:

[0086]

[0087] Where N is the total number of common edges between triangular macro-surfaces after the surface of the scatterer device shell is divided into many triangular macro-surfaces. Edges that cannot form triangular macro-surfaces and the junction of the illuminated area and the shadow area are not included; γ n are the unknown coefficients to be solved.

[0088] By adopting such a technical solution, the induced current density of any pair of adjacent triangular large-surface elements can be calculated, and the induced current density of all triangular large-surface elements can be determined based on the induced current density of any pair of adjacent triangular large-surface elements and the total number of triangular large-surface elements, thereby determining the induced current density on the surface of the scatterer. Taking into account the influence of adjacent triangular large-surface elements, the calculation accuracy of the induced current density on the surface of the scatterer is improved.

[0089] In one embodiment, S32 includes:

[0090] The induced current density on the surface of the scatterer is determined by the basis function model of the triangular large-surface element physical optics method based on the adjacent triangular large-surface elements, the total number of common edges of all triangular large-surface elements, the area of adjacent triangular large-surface elements, the length of the common edges of adjacent triangular large-surface elements, the first position vectors of adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edge.

[0091] For example, refer to Figure 3 , can be achieved by using the common edge l of the triangular large surface element n Introduce a unit vector t perpendicular to the midpoint of n + , t n - , thus obtaining the relationship between the unit vector and the triangular large-surface physical optics basis function model, and then obtaining the induced current density of the scatterer shell:

[0092]

[0093] Where N is the total number of common edges between triangular large-surface elements after the surface of the scatterer device shell is divided into many triangular large-surface elements, t n + , t n - The common edge l of the triangle large surface elementn The vertical unit vector introduced at the midpoint of n + and ρ n - Respectively represent any pair of adjacent triangular large surface elements by non-common vertices O + , O - The position vector indicated by ρ nc + and ρ nc - Respectively represent any pair of adjacent triangular large surface elements by non-common vertices O + , O - The position vector to the midpoint of the common edge, k n is the phase correction factor.

[0094] Thus, the electromagnetic scattering field of the scatterer to the target intelligent sensing device at this time is obtained:

[0095]

[0096] By adopting such a technical solution, the direction of the basis function is defined by the unit vector to ensure the continuity and physical rationality of the current. The unit vector is combined with the phase factor to optimize the propagation characteristics of high-frequency electromagnetic waves and the high-frequency calculation accuracy. It is possible to determine the induced current density on the surface of the scatterer through the triangular large-surface element physical optics basis function model based on the adjacent triangular large-surface elements, the total number of common edges of all triangular large-surface elements, the area of adjacent triangular large-surface elements, the length of the common edges of adjacent triangular large-surface elements, the first position vector of the adjacent triangular large-surface elements, and the second position vector from the adjacent triangular large-surface elements to the midpoint of the common edge. It is possible to combine the adjacent triangular large-surface elements, the total number of common edges of the triangular large-surface elements, and the unit vector from the adjacent triangular large-surface elements to the midpoint of the common edge to ensure the continuity of the current and optimize the high-frequency calculation accuracy. Reference Figure 4 When calculating the induced current on the target surface, the traditional physical optics method ignores the induced current in the shadow area. While this greatly simplifies the calculation process, speeds up the calculation, and reduces the computational memory, it makes numerous assumptions during the calculation, resulting in low computational accuracy and, therefore, only an approximate representation of the scattered field. The large-surface physical optics method considers that the scatterer's surface has both illuminated and shadowed areas formed by electromagnetic waves. Ignoring the current in the shadow area significantly reduces the computational effort, while maintaining computational accuracy through phase correction in the illuminated area. This method achieves both high computational accuracy and a significant reduction in computational effort, optimizing the construction of a simulation model of the scatterer's electromagnetic scattering field.

[0097] In one embodiment, the method further comprises:

[0098] 101. Acquire radiation field data of a target base station at a location of a target intelligent sensing device;

[0099] 102. Perform vector superposition based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0100] Specifically, the total electric field strength at the target intelligent sensing device is formed by the vector superposition of the radiation field generated by the target base station antenna at this location and the electromagnetic scattering field generated by the surrounding scatterers (such as the surrounding large-sized metal devices), that is:

[0101]

[0102] Where ω is the angular frequency of the electromagnetic wave emitted by the target base station antenna, j represents the imaginary unit, μ0 is the vacuum permeability, r and r′ represent the distance from point P and any unit surface element dS on the surface of the scatterer shell to the origin of the coordinate system O, R′ represents the distance between the unit surface element dS and point P, and H i is the magnetic field of the incident electromagnetic wave emitted by the target base station antenna, k and k′ are the electromagnetic wave number and the corresponding propagation vector, respectively, E T(r) is the electric field strength of the radiation field at point P. This can be calculated by calculating the radiation field generated by the target base station antenna at point P. This is not detailed here. N is the total number of common edges between the triangular large-surface elements after the surface of the scatterer device housing is divided into many triangular large-surface elements.

[0103] In this embodiment, the number of computational grids can be reduced by triangulating large surface elements, thereby reducing computational complexity. At the same time, the large surface element physical optics method can be combined to more accurately describe the propagation characteristics of high-frequency electromagnetic waves, taking into account both computational speed and computational accuracy. While ensuring computational accuracy, the computational complexity is simplified, thereby improving computational efficiency.

[0104] It should be understood that, although the various steps in the flowcharts involved in the various embodiments described above are displayed in sequence according to the instructions of the arrows, these steps are not necessarily executed in sequence in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order restriction on the execution of these steps, and these steps can be executed in other orders. Moreover, at least a portion of the steps in the flowcharts involved in the various embodiments described above can include multiple steps or multiple stages, and these steps or stages are not necessarily executed and completed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a portion of steps or stages in other steps.

[0105] Example 2

[0106] Based on the same inventive concept, embodiments of the present application also provide an electromagnetic scattering simulation device for implementing the aforementioned electromagnetic scattering simulation method. The solution provided by this device is similar to the solution described in the aforementioned method. Therefore, the specific limitations in one or more of the following electromagnetic scattering simulation device embodiments can be found in the above-described limitations on the electromagnetic scattering simulation method and are not further elaborated here.

[0107] In one embodiment, Figure 5 As shown, an electromagnetic scattering simulation device is provided, including: an acquisition module for acquiring a target base station in an area where a target intelligent sensing device is located and a scatterer within a radiation field formed by electromagnetic waves emitted by the target base station; a simulation module for dividing the surface of the scatterer into a plurality of triangular large-surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large-surface elements; a first calculation module for determining an induced current density on the surface of the scatterer using a large-surface physical optics method based on the scatterer electromagnetic scattering field simulation model; and a second calculation module for determining the electromagnetic scattering field of the scatterer to the target intelligent sensing device based on the induced current density.

[0108] Furthermore, the first calculation module is also used to establish a triangular large-surface element physical optics basis function model for any pair of adjacent triangular large-surface elements; and to determine the induced current density on the surface of the scatterer based on the triangular large-surface element physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

[0109] Furthermore, the first calculation module is also used to determine the induced current density on the surface of the scatterer through the triangular large surface element physical optics basis function model based on the adjacent triangular large surface elements, the total number of common edges of all the triangular large surface elements, the area of the adjacent triangular large surface elements, the length of the common edges of the adjacent triangular large surface elements, the first position vector of the adjacent triangular large surface elements, and the second position vector from the adjacent triangular large surface element to the midpoint of the common edge.

[0110] Furthermore, the simulation module is also used to construct a scatterer electromagnetic scattering field simulation model based on the electric field dyadic Green's function and the current continuity equation and a scatterer surface model composed of several triangular large surface elements in combination with preset simulation conditions.

[0111] Furthermore, the preset simulation conditions include: the scatterer is an ideal pure conductor, the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0112] Furthermore, the second calculation module is also used to obtain radiation field data of the target base station at the location of the target intelligent sensing device; and to perform vector superposition based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0113] Each module in the electromagnetic scattering simulation device described above can be implemented in whole or in part through software, hardware, or a combination thereof. Each module can be embedded in or independent of a processor in a computer device in hardware form, or can be stored in a computer device memory in software form, so that the processor can call and execute the corresponding operations of each module.

[0114] Example 3

[0115] In one embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the following steps are implemented:

[0116] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0117] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0118] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0119] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0120] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0121] For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established;

[0122] The induced current density on the surface of the scatterer is determined according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

[0123] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0124] The induced current density on the surface of the scatterer is determined based on the adjacent triangular large-surface elements, the total number of common edges of all the triangular large-surface elements, the areas of the adjacent triangular large-surface elements, the lengths of the common edges of the adjacent triangular large-surface elements, the first position vectors of the adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edges, using the triangular large-surface element physical optics basis function model.

[0125] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0126] According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

[0127] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0128] The scatterer is an ideal pure conductor, and the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0129] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0130] Acquire radiation field data of the target base station at the location of the target intelligent sensing device;

[0131] Vector superposition is performed based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0132] Example 4

[0133] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:

[0134] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0135] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0136] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0137] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0138] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0139] For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established;

[0140] The induced current density on the surface of the scatterer is determined according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

[0141] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0142] The induced current density on the surface of the scatterer is determined based on the adjacent triangular large-surface elements, the total number of common edges of all the triangular large-surface elements, the areas of the adjacent triangular large-surface elements, the lengths of the common edges of the adjacent triangular large-surface elements, the first position vectors of the adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edges, using the triangular large-surface element physical optics basis function model.

[0143] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0144] According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

[0145] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0146] The scatterer is an ideal pure conductor, and the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0147] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0148] Acquire radiation field data of the target base station at the location of the target intelligent sensing device;

[0149] Vector superposition is performed based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0150] Example 5

[0151] In one embodiment, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the following steps:

[0152] Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station;

[0153] Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements;

[0154] Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer;

[0155] An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

[0156] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0157] For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established;

[0158] The induced current density on the surface of the scatterer is determined according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

[0159] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0160] The induced current density on the surface of the scatterer is determined based on the adjacent triangular large-surface elements, the total number of common edges of all the triangular large-surface elements, the areas of the adjacent triangular large-surface elements, the lengths of the common edges of the adjacent triangular large-surface elements, the first position vectors of the adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edges, using the triangular large-surface element physical optics basis function model.

[0161] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0162] According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

[0163] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0164] The scatterer is an ideal pure conductor, and the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

[0165] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0166] Acquire radiation field data of the target base station at the location of the target intelligent sensing device;

[0167] Vector superposition is performed based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

[0168] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties.

[0169] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, database or other media used in the embodiments provided in this application may include at least one of non-volatile and volatile memory. Non-volatile memory may include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory may include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processor involved in the various embodiments provided herein may be, but are not limited to, a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, and the like.

[0170] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0171] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present application. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.

Claims

1. An electromagnetic scattering simulation method, characterized in that: The method comprises: Acquire a target base station in an area where a target intelligent sensing device is located and a scatterer in a radiation field formed by electromagnetic waves emitted by the target base station; Dividing the surface of the scatterer into a plurality of triangular large surface elements, and constructing an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements; Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer; An electromagnetic scattering field of the scatterer to the target intelligent sensing device is determined according to the induced current density.

2. The method according to claim 1, characterized in that Determining the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer includes: For any pair of adjacent triangular large-surface elements, a triangular large-surface element physical optics basis function model is established; The induced current density on the surface of the scatterer is determined according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements.

3. The method according to claim 2, characterized in that The determining of the induced current density on the surface of the scatterer according to the triangular large-surface physical optics basis function model and the total number of common edges of all the triangular large-surface elements includes: The induced current density on the surface of the scatterer is determined based on the adjacent triangular large-surface elements, the total number of common edges of all the triangular large-surface elements, the areas of the adjacent triangular large-surface elements, the lengths of the common edges of the adjacent triangular large-surface elements, the first position vectors of the adjacent triangular large-surface elements, and the second position vectors from the adjacent triangular large-surface elements to the midpoint of the common edges, using the triangular large-surface element physical optics basis function model.

4. The method according to claim 1, wherein The constructing of the electromagnetic scattering field simulation model of the scatterer according to the plurality of triangular large surface elements comprises: According to the electric field dyadic Green's function and the current continuity equation, a scatterer surface model composed of several triangular large surface elements is constructed in combination with preset simulation conditions to construct an electromagnetic scattering field simulation model of the scatterer.

5. The method according to claim 4, characterized in that The preset simulation conditions include: The scatterer is an ideal pure conductor, and the size of the scatterer is much larger than the wavelength of the electromagnetic wave emitted by the target base station, and the surface of the scatterer has an illuminated area and a shadow area formed by the electromagnetic wave.

6. The method according to claim 1, characterized in that The method further comprises: Acquire radiation field data of the target base station at the location of the target intelligent sensing device; Vector superposition is performed based on the radiation field data and the electromagnetic scattering field to determine the total electric field strength at the location of the target intelligent sensing device.

7. An electromagnetic scattering simulation device, characterized in that: The device comprises: An acquisition module is used to acquire a target base station in an area where a target intelligent sensing device is located and scatterers in a radiation field formed by electromagnetic waves emitted by the target base station; A simulation module, configured to divide the surface of the scatterer into a plurality of triangular large surface elements, and construct an electromagnetic scattering field simulation model of the scatterer based on the plurality of triangular large surface elements; A first calculation module is used to determine the induced current density on the surface of the scatterer by a large-surface physical optics method based on the electromagnetic scattering field simulation model of the scatterer; The second calculation module is used to determine the electromagnetic scattering field of the scatterer to the target intelligent sensing device according to the induced current density.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images