Method for constructing electromagnetic interference simulation model for shared iron tower

Abstract

The application relates to the technical field of data processing, and discloses a shared tower-oriented electromagnetic interference simulation model construction method, which comprises the following steps: establishing a three-dimensional terrain model of a region to be constructed, collecting electromagnetic data of all sampling positions of the three-dimensional terrain model; determining the geographical information complexity and sampling selection degree of the sampling positions, and screening and optimizing the sampling positions; obtaining equipotential lines and electromagnetic fluctuation degrees of the optimized sampling positions; and dividing the region to be constructed into different grid regions according to the electromagnetic data, the geographical information complexity and the electromagnetic fluctuation degrees of all the optimized sampling positions, and obtaining an electromagnetic interference simulation model of a shared tower. The application aims to adaptively divide the grid regions and improve the quality of the electromagnetic interference simulation model.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, in particular to a method for constructing an electromagnetic interference simulation model for a shared tower. BACKGROUND

[0002] A shared tower refers to installing communication equipment such as optical cables, communication base stations, and mobile antennas on a power tower to realize the co-construction and sharing of power and communication infrastructure. The various communication equipment installed on the power tower is susceptible to mutual influence and may produce electromagnetic interference, leading to problems such as signal attenuation and increased bit error rate. Meanwhile, the terrain affects the propagation path of electromagnetic waves and also changes the intensity, reflection, refraction, and other behaviors of the waves. Therefore, when establishing an electromagnetic interference simulation model for a shared tower, the terrain and electromagnetic data at different locations need to be considered simultaneously. Furthermore, to improve the accuracy of the electromagnetic interference simulation model, the electric field region needs to be divided into different grid regions during electromagnetic signal simulation.

[0003] When dividing the grid regions, the size of the grid regions affects the quality and simulation speed of the electromagnetic interference simulation model. A larger grid region can improve computational efficiency, but it may affect the authenticity of the simulation. A smaller grid region improves the simulation effect, but the dramatic increase in data volume may affect the normal operation of the equipment. Therefore, the specific division method of the grid regions needs to be determined based on the specific conditions of the terrain and electromagnetic data at different locations within the electric field region during electromagnetic signal simulation. SUMMARY

[0004] The present application provides a method for constructing an electromagnetic interference simulation model for a shared tower to solve the problem of unreasonable division of grid regions affecting the quality and simulation speed of the electromagnetic interference simulation model. The technical solution adopted is as follows:

[0005] One embodiment of the present application provides a method for constructing an electromagnetic interference simulation model for a shared tower, which includes the following steps:

[0006] Establish a three-dimensional terrain model of the to-be-constructed region of the electromagnetic interference simulation model, and collect electromagnetic data at all sampling locations of the three-dimensional terrain model multiple times;

[0007] Determine the geographical information complexity of each sampling location based on the difference between the heights of adjacent sampling locations and the height difference between each sampling location and the shared tower, determine the sampling selection preference of each sampling location based on the geographical information complexity of all sampling locations between the sampling location and the shared tower and the distance between different sampling locations, and select all preferred sampling locations based on the sampling selection preference;

[0008] According to all the electromagnetic data collected at all the preferred sampling positions in the three-dimensional terrain model, equipotential lines of each preferred sampling position are fitted respectively, and the electromagnetic fluctuation of each preferred sampling position is calculated according to the differences between the equipotential lines of different preferred sampling positions, the distances between different preferred sampling positions, and the differences between all the electromagnetic data collected at the same preferred sampling position.

[0009] According to the electromagnetic data, the geographical information complexity and the electromagnetic fluctuation of all the preferred sampling positions, the region to be constructed is divided into different grid regions, and the electromagnetic interference simulation model of the shared tower is obtained by processing all the grid regions divided from the region to be constructed, the electromagnetic data of all the sampling positions in all the grid regions, and the three-dimensional terrain model of the region to be constructed using electromagnetic wave propagation simulation software.

[0010] Further, the method for obtaining the geographical information complexity of the sampling position is as follows:

[0011] Any one sampling position is recorded as a target sampling position, all the sampling positions in a window centered on the target sampling position are recorded as the adjacent sampling positions of the target sampling position, and the standard deviation of the heights of all the adjacent sampling positions of the target sampling position is recorded as the first standard deviation of the target sampling position. The absolute value of the difference between the height of the target sampling position and the height of the shared tower is recorded as the first absolute value of the target sampling position.

[0012] The geographical information complexity of the target sampling position is determined according to the first standard deviation and the first absolute value of the target sampling position.

[0013] Further, the method for determining the geographical information complexity of the target sampling position according to the first standard deviation and the first absolute value of the target sampling position includes the following specific method:

[0014] The sum of the first standard deviation and the first absolute value of the target sampling position is recorded as the geographical information complexity of the target sampling position.

[0015] Further, the method for determining the sampling selection degree of the sampling position is as follows:

[0016] The cumulative sum of the geographical information complexity of all the sampling positions contained in the line segment determined by the target sampling position and the shared tower as the end points is recorded as the first cumulative sum of the target sampling position.

[0017] The distance between the target sampling position and the nearest sampling position is recorded as the first distance of the target sampling position.

[0018]

[0019] ​The ratio of the product of the geographical information complexity of the target sampling position and the first distance to the first accumulated sum is recorded as the sampling selection priority of the target sampling position.

[0020] Further, the screening method of the preferred sampling position is:

[0021] The sampling position with the sampling selection priority greater than the preset preferred threshold is recorded as the preferred sampling position.

[0022] Further, the method for obtaining the equipotential line of the preferred sampling position is:

[0023] The average of all electromagnetic data collected at the same preferred sampling position of the three-dimensional terrain model is recorded as the average electromagnetic data of the same preferred sampling position.

[0024] Any one preferred sampling position is recorded as the target preferred sampling position. The average electromagnetic data collected at all preferred sampling positions of the three-dimensional terrain model is used to interpolate and calculate the electromagnetic data of different positions, to obtain the electromagnetic data interpolation of all positions equal to the average electromagnetic data of the target preferred sampling position. The average electromagnetic data of the target preferred sampling position and the electromagnetic data interpolation of the positions equal to the average electromagnetic data of the target preferred sampling position are used for curve fitting to obtain the equipotential line of the target preferred sampling position.

[0025] Further, the electromagnetic fluctuation degree of each preferred sampling position is calculated according to the difference between the equipotential lines of different preferred sampling positions, the distance between different preferred sampling positions, and the difference between all electromagnetic data collected at the same preferred sampling position, including the specific method:

[0026] The adjacent equipotential line difference of the target preferred sampling position is determined according to the difference between the equipotential lines of different preferred sampling positions. The distance between the two preferred sampling positions corresponding to the adjacent equipotential line difference of the target preferred sampling position is recorded as the adjacent equipotential line distance of the target preferred sampling position. The ratio of the adjacent equipotential line difference to the adjacent equipotential line distance of the target preferred sampling position is recorded as the planar affected degree of the target preferred sampling position.

[0027] The normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling position and the planar affected degree of the target preferred sampling position is recorded as the electromagnetic fluctuation degree of the target preferred sampling position.

[0028] Further, the method for determining the adjacent equipotential line difference of the target preferred sampling position according to the difference between the equipotential lines of different preferred sampling positions includes the specific method:

[0029] Calculate the Frechet distance between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position. The minimum Frechet distance between the target preferred sampling position and the equipotential line of the target preferred sampling position is recorded as the difference between adjacent equipotential lines of the target preferred sampling position.

[0030] Furthermore, the specific method for dividing the region to be constructed into different grid regions based on the electromagnetic data, geographic information complexity, and electromagnetic fluctuation of all preferred sampling locations includes:

[0031] The wavelength corresponding to the electromagnetic data at the same preferred sampling position is taken as the x-axis value, the frequency corresponding to the electromagnetic data at the same preferred sampling position is taken as the y-axis value, and the electromagnetic fluctuation at the preferred sampling position is taken as the z-axis value. Surface fitting is performed on the data points corresponding to different frequencies and wavelengths of the electromagnetic data at the same preferred sampling position to obtain the fitted surface of the same preferred sampling position.

[0032] Based on the electromagnetic fluctuation degree and geographical information complexity of all preferred sampling locations, calculate the stopping feature value corresponding to the target preferred sampling location when it extends to different numbers of preferred sampling locations in different directions, until the stopping feature value is greater than or equal to the preset stopping threshold, and then stop calculating the stopping feature value in the corresponding direction.

[0033] The preferred sampling positions corresponding to all stopping feature values ​​less than the stopping threshold corresponding to the target preferred sampling position are denoted as the sampling positions of the same sub-region; the convex hull of all preferred sampling positions of the same sub-region sampling positions is denoted as the grid region; and the region composed of different grid regions that overlap is considered as a single grid region.

[0034] All regions in the area to be constructed that have not been assigned to a grid region are denoted as grid regions.

[0035] Furthermore, the formula for calculating the stopping feature value is:

[0036]

[0037] in, Indicates the first The preferred sampling position is oriented in the direction extend The stopping feature value corresponding to the preferred sampling position; Indicates the first The preferred sampling position is oriented in the direction The extended first The geographical information complexity of each preferred sampling location; Indicates the first The preferred sampling position is oriented in the direction The extended first Electromagnetic fluctuations at preferred sampling locations.

[0038] The beneficial effects of the present application are:

[0039] The present application first considers that when electromagnetic information is interfered by electromagnetic interference and ground type, the electromagnetic information of adjacent positions is relatively complex, the electromagnetic information of different sampling positions is relatively dense when constructing an electromagnetic interference simulation model according to the electromagnetic information, the complexity of the electromagnetic information of the sampling position is evaluated, the geographical information complexity of the sampling position is obtained, and the sampling selection degree of the sampling position is obtained in combination with the effectiveness of the sampling position, all optimal sampling positions are selected according to the sampling selection degree, that is, reasonable sampling points for constructing an electromagnetic interference simulation model are selected; when the signal transmission on the shared tower is radioactive, and the signal is interfered by buildings, terrain and other reasons in the transmission process, the transmission signals with equal distance to the source in the same two-dimensional plane will be different, the equipotential lines of each optimal sampling position are constructed respectively, and the electromagnetic fluctuation degree of each optimal sampling position is calculated according to the equipotential lines; finally, according to the electromagnetic data, geographical information complexity and electromagnetic fluctuation degree of all optimal sampling positions, the region to be constructed is adaptively divided into grid regions, the problem that the division of the grid region is unreasonable, affecting the quality and simulation speed of the electromagnetic interference simulation model is solved, and the electromagnetic interference simulation model of the shared tower is obtained by using electromagnetic wave propagation simulation software according to the divided grid region, and the quality of the electromagnetic interference simulation model is improved. BRIEF DESCRIPTION OF DRAWINGS

[0040] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description only constitute some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor.

[0041] Figure 1 The flowchart of the method for constructing an electromagnetic interference simulation model for a shared tower provided by an embodiment of the present application is shown in the figure.

[0042] Figure 2 The schematic diagram of the equipotential line of the optimal sampling position provided by an embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0043] The technical solutions in the embodiments of the present application will be described clearly and completely in combination with the drawings in the embodiments of the present application. Obviously, the described embodiments only constitute some embodiments of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0044] Referring to Figure 1 , which shows a flow chart of a method for constructing an electromagnetic interference simulation model for a shared tower according to an embodiment of the present application. The method comprises the following steps:

[0045] Step S001, a three-dimensional terrain model of a region to be constructed of the electromagnetic interference simulation model is established, and electromagnetic data of all sampling positions of the three-dimensional terrain model is collected multiple times.

[0046] DEM data of the region to be constructed of the electromagnetic interference simulation model is collected using a remote sensing satellite, LiDAR point cloud data of the region to be constructed of the electromagnetic interference simulation model is collected using a LiDAR scanner carried by a drone, and a three-dimensional terrain model is constructed from the DEM data and the LiDAR point cloud data.

[0047] The DEM data is digital elevation model data, which is usually raster data, and each point in the LiDAR point cloud data has a corresponding three-dimensional coordinate. The construction of a three-dimensional terrain model from DEM data and LiDAR point cloud data is a known technique and will not be described in detail.

[0048] The test personnel obtain electromagnetic data at different sampling positions of the three-dimensional terrain model using a handheld signal receiving device, and obtain electromagnetic data at the same sampling position multiple times.

[0049] Preferably, in an embodiment of the present application, the number of times of obtaining electromagnetic data at the same sampling position is 5. In actual application, as other implementation manners, the implementer can determine the value of the number of times of obtaining electromagnetic data at the same sampling position according to actual conditions, and the present application does not make special limitations.

[0050] At this point, the three-dimensional terrain model and the electromagnetic data collected multiple times at different sampling positions of the three-dimensional terrain model are obtained.

[0051] Step S002, the geographical information complexity of each sampling position is determined according to the difference between the heights of adjacent sampling positions and the difference between the height of each sampling position and the height of the shared tower, the sampling selection priority of each sampling position is determined according to the geographical information complexity of all sampling positions between the sampling position and the shared tower and the distance between different sampling positions, and all preferred sampling positions are selected according to the sampling selection priority.

[0052] The various communication devices installed on the tower are susceptible to mutual interference, potentially causing electromagnetic interference that leads to signal attenuation and increased bit error rate. When the electromagnetic wave source is located at a low altitude, obstacles such as tall buildings and mountains can cause shadows to appear on the waves, while when the electromagnetic wave source is located at a high altitude, the electromagnetic waves can propagate further. Different types of ground surfaces, such as water bodies, forests, and cities, have different characteristics in terms of reflection, diffraction, and absorption. For example, water bodies reflect electromagnetic waves strongly, while forests attenuate signals significantly. Therefore, different types of ground surfaces also affect the propagation of electromagnetic waves. Thus, terrain affects the propagation path of electromagnetic waves and can also change the wave's intensity, reflection, and refraction. When establishing an electromagnetic interference simulation model for a shared tower, it is necessary to consider both the terrain and electromagnetic data from different locations simultaneously.

[0053] When electromagnetic information is affected by electromagnetic interference and ground-type interference, the electromagnetic information at adjacent locations is quite complex. Therefore, when constructing an electromagnetic interference simulation model based on this electromagnetic information, the electromagnetic information at different sampling locations should also be relatively dense. At the same time, it is also necessary to consider the effectiveness of each sampling location and select reasonable sampling points for constructing the electromagnetic interference simulation model.

[0054] FDTD (Finite-Difference Time-Domain) simulation is a numerical simulation method used to calculate the distribution and propagation of electromagnetic fields in the time domain. The sampling locations used in FDTD simulations are the appropriate sampling points for constructing electromagnetic interference simulation models. The setting of the sampling locations in FDTD simulations has a significant impact on the effectiveness of the established electromagnetic interference simulation model; therefore, it is necessary to set the sampling locations appropriately. Ideally, the sampling locations in FDTD simulations should be located at key locations in the region to be constructed, such as near the wave source, the center of the region of interest, or areas where the field distribution changes drastically.

[0055] The geographic information complexity of each sampling location is determined based on the height difference between adjacent sampling locations and the height difference between each sampling location and the shared tower.

[0056] Let any sampling location be designated as the target sampling location, and then define the sampling location as the center of the target sampling location. All sampling locations within the window are denoted as neighboring sampling locations of the target sampling location. The standard deviation of the height of all neighboring sampling locations of the target sampling location is denoted as the first standard deviation of the target sampling location. The absolute value of the difference between the height of the target sampling location and the height of the shared tower is denoted as the first absolute value of the target sampling location. The sum of the first standard deviation and the first absolute value of the target sampling location is denoted as the geographic information complexity of the target sampling location.

[0057] The same method can be used to obtain the geographic information complexity of each sampling location in the 3D terrain model.

[0058] The greater the geographical information complexity of the sampling position is, the more complex the electromagnetic information of the sampling position is, the more dense the electromagnetic information provided by the sampling position is, and the greater the possibility that the sampling position is located at a lower position is. When the sampling position is selected as a reasonable sampling point of the FDTD simulation, the possibility that a shadow zone of a wave appears is greater.

[0059] The sampling selection priority of each sampling position is determined according to the geographical information complexity of all sampling positions between the sampling position and the shared tower and the distance between different sampling positions.

[0060] An accumulated sum of the geographical information complexity of all sampling positions contained in a line segment with the target sampling position and the shared tower as end points is recorded as a first accumulated sum of the target sampling position. A distance between the target sampling position and the nearest sampling position is recorded as a first distance of the target sampling position. A ratio of a product of the geographical information complexity of the target sampling position and the first distance to the first accumulated sum is recorded as the sampling selection priority of the target sampling position.

[0061] The first accumulated sum of the target sampling position is used to evaluate the complexity of electromagnetic information of each sampling position in a path from the target sampling position to the shared tower. The greater the complexity of electromagnetic information of each sampling position in the path from the target sampling position to the shared tower is, the greater the sampling selection priority of the target sampling position is. When the target sampling position is selected as a reasonable sampling point of the FDTD simulation, more dense electromagnetic information can be extracted. In the process of calculating the sampling selection priority of the target sampling position, the first distance of the target sampling position is used to avoid that the set reasonable sampling points of the FDTD simulation are too dense.

[0062] The sampling selection priority of each sampling position of the three-dimensional terrain model can be obtained in the same way.

[0063] The sampling position with the sampling selection priority greater than a preset preferred threshold is recorded as a preferred sampling position.

[0064] It can be understood that the preferred sampling position is a reasonable sampling point of the FDTD simulation, that is, a reasonable sampling point for constructing an electromagnetic interference simulation model.

[0065] At this point, all the preferred sampling positions are screened out.

[0066] In step S003, the equipotential lines of each preferred sampling position are fitted according to all electromagnetic data collected by all the preferred sampling positions in the three-dimensional terrain model. The electromagnetic fluctuation degree of each preferred sampling position is calculated according to the differences between the equipotential lines of different preferred sampling positions, the distances between different preferred sampling positions, and the differences between all electromagnetic data collected by the same preferred sampling position.

[0067] The signal transmission on the shared iron tower is radioactive, the transmission signals with equal distance to the transmission source in the same two-dimensional plane should be consistent, but the transmission signals are interfered by buildings, terrain and other reasons in the transmission process, and the transmission signals with equal distance to the transmission source in the same two-dimensional plane will exist differences, therefore, the equipotential lines can be constructed, and the complexity of the interference of the electromagnetic signal is determined according to the equipotential lines.

[0068] According to all electromagnetic data collected by all preferred sampling positions in the three-dimensional terrain model, the equipotential line of each preferred sampling position is fitted respectively.

[0069] Preferably, as an embodiment of the present application, the mean value of all electromagnetic data collected by the same preferred sampling position in the three-dimensional terrain model is recorded as the average electromagnetic data of the same preferred sampling position. Any one preferred sampling position is recorded as a target preferred sampling position, the electromagnetic data of different positions is calculated by interpolation according to the average electromagnetic data collected by all preferred sampling positions in the three-dimensional terrain model, the electromagnetic data interpolation of all positions equal to the average electromagnetic data of the target preferred sampling position is obtained, and the equipotential line of the target preferred sampling position is obtained by curve fitting according to the average electromagnetic data of the target preferred sampling position and the electromagnetic data interpolation of the positions equal to the average electromagnetic data of the target preferred sampling position.

[0070] Wherein, the interpolation of the electromagnetic data and the curve fitting of the data are all known technologies, and will not be repeated; specifically, the cubic spline interpolation is used to interpolate the electromagnetic data, and the least square method is used for curve fitting.

[0071] When the electromagnetic data of each preferred sampling position is not interfered, the equipotential lines of all preferred sampling positions are circular, when the electromagnetic data is interfered, the position of the equipotential line of the preferred sampling position will be deformed. Due to the interference of various communication equipment installed on the power iron tower, the type of ground and other factors, the equipotential line of each preferred sampling position is a deformed circle, and the schematic diagram of the equipotential line of the preferred sampling position is shown in 2. Figure 2 In the figure, 1 represents the preferred sampling position, and 2 represents the equipotential line of the preferred sampling position.

[0072] The Frechet distance between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position is calculated respectively, the minimum value of the Frechet distance with the equipotential line of the target preferred sampling position is recorded as the adjacent equipotential line difference of the target preferred sampling position, the distance between the two preferred sampling positions corresponding to the adjacent equipotential line difference of the target preferred sampling position is recorded as the adjacent equipotential line distance of the target preferred sampling position, and the ratio of the adjacent equipotential line difference and the adjacent equipotential line distance of the target preferred sampling position is recorded as the plane affected degree of the target preferred sampling position.

[0073] The equipotential line is a curve, and the Frechet distance between different curves is a known technology and will not be described again.

[0074] When the plane affected degree of the target preferred sampling position is smaller, the grid area divided by the electric field area where the target preferred sampling position is located should be smaller, so as to improve the quality of the electromagnetic interference simulation model.

[0075] The normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling position and the plane affected degree of the target preferred sampling position is recorded as the electromagnetic fluctuation degree of the target preferred sampling position.

[0076] It should be noted that the Z-Score standard normalization method is used to calculate the normalized value in this embodiment, and other methods such as the maximum and minimum value normalization method, sigmoid function and other methods of prior art can be used to calculate the normalized value in actual application, which is not limited here.

[0077] The electromagnetic fluctuation degree of each preferred sampling position of the three-dimensional terrain model can be obtained in the same way.

[0078] Thus, the electromagnetic fluctuation degree of all preferred sampling positions of the three-dimensional terrain model is obtained.

[0079] Step S004, according to the electromagnetic data, geographical information complexity and electromagnetic fluctuation degree of all preferred sampling positions, the region to be constructed is divided into different grid areas, and the electromagnetic wave propagation simulation software is used to process all grid areas divided by the region to be constructed, electromagnetic data of all sampling positions in all grid areas and three-dimensional terrain model of the region to be constructed, and obtain the electromagnetic interference simulation model of the shared tower.

[0080] The wavelength corresponding to the electromagnetic data of the same preferred sampling position is taken as the x-axis value, the frequency corresponding to the electromagnetic data of the same preferred sampling position is taken as the y-axis value, and the electromagnetic fluctuation degree of the preferred sampling position is taken as the z-axis value. The least square method is used to fit the data points corresponding to different frequencies and wavelengths of the electromagnetic data of the same preferred sampling position to obtain the fitting surface of the same preferred sampling position.

[0081] Surface fitting for different data points to obtain the fitted surface is a well-known technique and will not be elaborated further.

[0082] Based on the electromagnetic fluctuation degree and geographical information complexity of all preferred sampling locations, the stopping characteristic value corresponding to the extension of the target preferred sampling location to different numbers of preferred sampling locations in different directions is calculated. The calculation formula is as follows:

[0083]

[0084] in, Indicates the first The preferred sampling position is oriented in the direction extend The stopping feature value corresponding to the preferred sampling position; Indicates the first The preferred sampling position is oriented in the direction The extended first The geographical information complexity of each preferred sampling location; Indicates the first The preferred sampling position is oriented in the direction The extended first Electromagnetic fluctuations at preferred sampling locations.

[0085] in, The value starts from 1 and increases gradually in increments of 1 until the stopping feature value is greater than or equal to a preset stopping threshold; direction It can be any direction within the plane.

[0086] The stopping threshold is set to 0.75. The preferred sampling positions corresponding to all stopping feature values ​​less than the stopping threshold at the target preferred sampling position are denoted as the same sub-region sampling positions.

[0087] The convex hull of all preferred sampling locations that are determined to be sampling locations in the same sub-region is denoted as the grid region.

[0088] It is understandable that when different grid regions overlap, the region consisting of the overlapping different grid regions is considered as a single grid region.

[0089] All regions in the area to be constructed that have not been assigned to a grid region are denoted as grid regions.

[0090] At this point, the area to be constructed is divided into different grid areas.

[0091] Using electromagnetic wave propagation simulation software, an electromagnetic interference simulation model of the shared tower is obtained based on the electromagnetic data of all grid areas divided into the area to be constructed, the electromagnetic data of all sampling locations within all grid areas, and the three-dimensional terrain model of the area to be constructed.

[0092] The electromagnetic wave propagation simulation software can be CST Studio Suite, COMSOL Multiphysics, etc., and the electromagnetic interference simulation model obtained by the electromagnetic wave propagation simulation software is a known technology and will not be described here.

[0093] At this point, the construction of the electromagnetic interference simulation model of the shared tower is completed.

[0094] The above only describes the preferred embodiments of the present application and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for constructing an electromagnetic interference simulation model for a shared tower, characterized in that, The method comprises the following steps: A three-dimensional terrain model of a to-be-constructed area of an electromagnetic interference simulation model is established, and electromagnetic data of all sampling positions of the three-dimensional terrain model is collected multiple times; According to the difference between the heights of adjacent sampling positions and the height difference between each sampling position and the shared tower, the geographical information complexity of each sampling position is determined respectively, any one sampling position is recorded as a target sampling position, the geographical information complexity of all sampling positions contained in the line segment determined by the target sampling position and the shared tower as end points is recorded as the first cumulative sum of the target sampling position, the distance between the target sampling position and the nearest sampling position is recorded as the first distance of the target sampling position, and the ratio of the product of the geographical information complexity of the target sampling position and the first distance to the first cumulative sum is recorded as the sampling selection priority of the target sampling position, and all preferred sampling positions are selected according to the sampling selection priority; According to all electromagnetic data collected at all preferred sampling positions in the three-dimensional terrain model, the equipotential lines of each preferred sampling position are fitted respectively, the adjacent equipotential line difference of the target preferred sampling position is determined according to the difference between the equipotential lines of different preferred sampling positions, the distance between the two preferred sampling positions corresponding to the adjacent equipotential line difference of the target preferred sampling position is recorded as the adjacent equipotential line distance of the target preferred sampling position, and the ratio of the adjacent equipotential line difference to the adjacent equipotential line distance of the target preferred sampling position is recorded as the plane affected degree of the target preferred sampling position; and the normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling position and the plane affected degree of the target preferred sampling position is recorded as the electromagnetic fluctuation degree of the target preferred sampling position. The wavelength corresponding to the electromagnetic data of the same preferred sampling position is taken as the x-axis value, the frequency corresponding to the electromagnetic data of the same preferred sampling position is taken as the y-axis value, the electromagnetic fluctuation degree of the preferred sampling position is taken as the z-axis value, the data points corresponding to different frequencies and wavelengths of the electromagnetic data of the same preferred sampling position are subjected to surface fitting to obtain the fitting surface of the same preferred sampling position; the stopping characteristic value corresponding to the extension of the target preferred sampling position to different directions by different numbers of preferred sampling positions is calculated according to the electromagnetic fluctuation degree and the geographical information complexity of all preferred sampling positions, until the stopping characteristic value is greater than or equal to the preset stopping threshold value, the stopping characteristic value corresponding to the direction is stopped; the preferred sampling positions corresponding to all stopping characteristic values less than the stopping threshold value of the target preferred sampling position are recorded as the same sub-area sampling positions; the convex hull of all preferred sampling positions of the same sub-area sampling positions is recorded as the grid area; the area composed of different grid areas with overlaps is taken as a grid area; all areas in the to-be-constructed area that are not divided into grid areas are recorded as grid areas; and the electromagnetic interference simulation model of the shared tower is obtained by processing the three-dimensional terrain model of the to-be-constructed area, the electromagnetic data of all sampling positions in all grid areas divided by the electromagnetic wave propagation simulation software and all grid areas.

2. The method of claim 1, wherein, The method for obtaining the geographical information complexity of the sampling position is as follows: Centered on the target sampling location All sampling positions within the window are recorded as neighboring sampling positions of the target sampling position, and the standard deviation of the height of all neighboring sampling positions of the target sampling position is recorded as the first standard deviation of the target sampling position. An absolute value of a difference between the target sampling position and the height of the shared tower is recorded as a first absolute value of the target sampling position; According to the first standard deviation and the first absolute value of the target sampling position, a geographical information complexity of the target sampling position is determined.

3. The method of claim 2, wherein the method further comprises: The method of determining the geographical information complexity of the target sampling position according to the first standard deviation and the first absolute value of the target sampling position includes the following specific method: A sum of the first standard deviation and the first absolute value of the target sampling position is recorded as the geographical information complexity of the target sampling position.

4. The method of claim 1, wherein, The method of screening the preferred sampling position includes the following method: Sampling positions with a sampling selection degree greater than a preset preferred threshold are all recorded as the preferred sampling position.

5. The method of claim 1, wherein, The method of obtaining the equipotential line of the preferred sampling position includes the following method: A mean value of all electromagnetic data collected by the same preferred sampling position from the three-dimensional terrain model is recorded as average electromagnetic data of the same preferred sampling position. Any one of the preferred sampling positions is recorded as a target preferred sampling position, average electromagnetic data of all the preferred sampling positions from the three-dimensional terrain model is used to perform interpolation calculation on electromagnetic data of different positions, electromagnetic data interpolation of all positions with the same average electromagnetic data as the target preferred sampling position is obtained, curve fitting is performed on the average electromagnetic data of the target preferred sampling position and the electromagnetic data interpolation of the positions with the same average electromagnetic data as the target preferred sampling position, and an equipotential line of the target preferred sampling position is obtained.

6. The method of claim 1, wherein, The method of determining the adjacent equipotential line difference of the target preferred sampling position according to the difference between the equipotential lines of different preferred sampling positions includes the following specific method: Frechet distances between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position are respectively calculated, and a minimum value of the Frechet distances between the equipotential line of the target preferred sampling position is recorded as the adjacent equipotential line difference of the target preferred sampling position.

7. The method of claim 1, wherein the method is a method of constructing an electromagnetic interference simulation model for a shared tower. The calculation formula of the stop feature value is wherein, denotes the th preferred sampling location extends in the direction denotes the th preferred sampling location when the corresponding stop feature value; denotes the th preferred sampling location extends in the direction denotes the geographical information complexity of the th preferred sampling location; denotes the th preferred sampling location extends in the direction denotes the electromagnetic fluctuation of the th preferred sampling location.

Images