Electromagnetic situation information determination method and device, satellite and storage medium

Abstract

The invention relates to an electromagnetic situation information determination method and device, a satellite and a storage medium. The method comprises the following steps: executing multiple interference power detection in a satellite coverage range of an own satellite, and determining electromagnetic situation information corresponding to the satellite coverage range according to an interference power detection result; wherein the ith interference power detection comprises the following steps: acquiring an initial sensing point set which comprises a plurality of first sensing points, and acquiring an interference power value at each first sensing point; a plurality of second sensing points are generated in the satellite coverage range according to the interference power values corresponding to the first sensing points, a target sensing point set is generated according to the second sensing points, and distribution of the second sensing points is determined according to the interference power values corresponding to the first sensing points; the target sensing point set is an initial sensing point set of the (i + 1) th interference power detection. By adopting the method, space electromagnetic situation information sensing resources can be saved.

Description

Technical Field

[0001] The present application relates to the field of satellite communication technology, and in particular to a method, device, satellite and storage medium for determining electromagnetic situation information. Background Art

[0002] Due to the rapid increase in the number of low-orbit satellites, there is a large amount of overlap in the spectrum resources declared by different operators, especially in the Ku and Ka bands, which easily leads to co-channel interference. In the uplink co-channel interference scenario where low-orbit constellations coexist, the satellite system of one's own uses the sensing beam to actively perceive the spatial spectrum situation of the other constellation, thereby using the known spectrum distribution of the other system to avoid potential co-channel interference. This is a key technical means to achieve efficient coexistence and sharing of spectrum resources between one's own system and other systems.

[0003] For an unknown co-channel interference environment of another system, the ideal situation is that the satellite senses all possible uplink wave directions. However, when the number of spatial sampling points is too large, the sensing process consumes a large amount of time slots and beam resources. Summary of the invention

[0004] Based on this, it is necessary to provide an electromagnetic situation information determination method, device, satellite and storage medium that can save space electromagnetic situation information perception resources and shorten the space electromagnetic situation information generation cycle in response to the above technical problems.

[0005] In a first aspect, the present application provides a method for determining electromagnetic situation information, comprising:

[0006] Perform multiple interference power detections within the satellite coverage of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage according to the results of the interference power detection. The results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to the own satellite;

[0007] The i-th interference power detection in the multiple interference power detections includes:

[0008] Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each first sensing point;

[0009] According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set of the i+1th interference power detection.

[0010] In one embodiment, generating a plurality of second sensing points within a satellite coverage range according to interference power values ​​corresponding to the first sensing points includes:

[0011] According to the interference power value corresponding to each first perception point, multiple first perception points are clustered to obtain multiple perception point clusters; according to the clustering power value corresponding to each perception point cluster, the number of newly added perception points corresponding to each perception point cluster is determined; according to the number of newly added perception points corresponding to each perception point cluster, second perception points are generated in the perception area corresponding to each perception point cluster, wherein the perception area corresponding to the perception point cluster is the area where the perception points included in the perception point cluster are located within the satellite coverage range.

[0012] In one of the embodiments, clustering is performed on a plurality of first perception points according to interference power values ​​corresponding to the first perception points to obtain a plurality of perception point clusters, including:

[0013] According to the interference power value corresponding to each first perception point, a plurality of first perception points are clustered to obtain an initial perception point cluster; for each initial perception point cluster, the initial perception point cluster is clustered according to the position information of each first perception point in the initial perception point cluster to obtain a plurality of sub-clusters included in the initial perception point cluster; the sub-clusters included in each initial perception point cluster are used as perception point clusters.

[0014] In one of the embodiments, the number of newly added sensing points corresponding to the sensing point cluster is positively correlated with the clustering power value corresponding to the sensing point cluster.

[0015] In one embodiment, determining the number of newly added sensing points corresponding to each sensing point cluster according to the clustering power value corresponding to each sensing point cluster includes:

[0016] The total number of newly added sensing points in a single round is obtained, and the allocation ratio corresponding to each sensing point cluster is determined according to the clustering power value corresponding to each sensing point cluster; according to the allocation ratio corresponding to each sensing point cluster, the number of newly added sensing points is allocated to each sensing point cluster from the total number, wherein the sum of the number of newly added sensing points allocated to each sensing point cluster is equal to the total number.

[0017] In one embodiment, determining the allocation ratio corresponding to each sensing point cluster according to the clustering power value corresponding to each sensing point cluster includes:

[0018] The cluster position corresponding to each sensing point cluster is obtained; for each sensing point cluster, the allocation ratio is determined according to the cluster power value corresponding to the sensing point cluster, the cluster position and the number of first sensing points included in the sensing point cluster.

[0019] In one of the embodiments, generating a target perception point set according to a plurality of second perception points comprises one of the following:

[0020] taking the plurality of second perception points as a target perception point set;

[0021] A target number of perception points are selected from the plurality of second perception points and the plurality of first perception points to serve as a target perception point set.

[0022] In one embodiment, generating a target perception point set according to a plurality of second perception points includes:

[0023] Obtain the maximum number of sensing points and the target number of sensing points within the satellite coverage, wherein the target number includes the number of sensing points generated during previous interference power detection processes and the number of sensing points in the initial state; determine whether a loop stop condition is met based on the maximum number of sensing points and the target number; if the loop stop condition is not met, execute the step of generating a target sensing point set based on multiple second sensing points.

[0024] In one embodiment, determining whether a loop stop condition is satisfied based on the maximum number of sensing points and the target number includes:

[0025] Determine the difference between the maximum number of perception points and the target number; if the difference is greater than or equal to the total number of newly added perception points in a single round, it is determined that the loop stop condition is not met; if the difference is less than the total number of newly added perception points in a single round, it is determined that the loop stop condition is met.

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

[0027] If the loop stop condition is met, the number of additional sensing points is determined based on the difference between the maximum number of sensing points and the target number; multiple third sensing points are generated within the satellite coverage range based on the number of additional sensing points; the interference power at each third sensing point is detected to obtain the interference power value corresponding to each third sensing point.

[0028] In one embodiment, determining electromagnetic situation information corresponding to a satellite coverage range according to a result of interference power detection includes:

[0029] The electromagnetic situation information is determined based on the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

[0030] In one embodiment, the electromagnetic situation information is determined according to the interference power value corresponding to each sensing point within the satellite coverage range determined by each interference power detection, including:

[0031] Obtain the location information of each sensing point within the satellite coverage area; perform interpolation fitting on the location information and interference power value of each sensing point based on a preset interpolation method to determine the electromagnetic situation information.

[0032] In a second aspect, the present application also provides an electromagnetic situation information determination device, comprising:

[0033] A determination module is used to perform multiple interference power detections within the satellite coverage of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage according to the results of the interference power detection. The results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to the own satellite;

[0034] The i-th interference power detection in the multiple interference power detections includes:

[0035] Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each first sensing point;

[0036] According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set of the i+1th interference power detection.

[0037] In one of the embodiments, the determination module is specifically used to cluster multiple first perception points according to the interference power value corresponding to each first perception point to obtain multiple perception point clusters; determine the number of newly added perception points corresponding to each perception point cluster according to the cluster power value corresponding to each perception point cluster; and generate second perception points in the perception area corresponding to each perception point cluster according to the number of newly added perception points corresponding to each perception point cluster, wherein the perception area corresponding to the perception point cluster is the area where the perception points included in the perception point cluster are located within the satellite coverage range.

[0038] In one of the embodiments, the determination module is specifically used to cluster multiple first perception points according to the interference power value corresponding to each first perception point to obtain an initial perception point cluster; for each initial perception point cluster, the initial perception point cluster is clustered according to the position information of each first perception point in the initial perception point cluster to obtain multiple sub-clusters included in the initial perception point cluster; and the sub-clusters included in each initial perception point cluster are used as perception point clusters.

[0039] In one of the embodiments, the number of newly added sensing points corresponding to the sensing point cluster is positively correlated with the clustering power value corresponding to the sensing point cluster.

[0040] In one of the embodiments, the determination module is specifically used to obtain the total number of newly added perception points in a single round, and determine the allocation ratio corresponding to each perception point cluster according to the clustering power value corresponding to each perception point cluster; according to the allocation ratio corresponding to each perception point cluster, allocate the number of newly added perception points to each perception point cluster from the total number, wherein the sum of the number of newly added perception points allocated to each perception point cluster is equal to the total number.

[0041] In one of the embodiments, the determination module is specifically used to obtain the cluster position corresponding to each perception point cluster; for each perception point cluster, the allocation ratio is determined according to the cluster power value corresponding to the perception point cluster, the cluster position and the number of first perception points included in the perception point cluster.

[0042] In one embodiment, the determination module is specifically configured to select the plurality of second perception points as a target perception point set; or select a target number of perception points from the plurality of second perception points and the plurality of first perception points as a target perception point set.

[0043] In one of the embodiments, the determination module is specifically used to obtain the maximum number of perception points and the target number of perception points included in the satellite coverage, the target number including the number of perception points generated in previous interference power detection processes and the number of perception points in the initial state; based on the maximum number of perception points and the target number, determine whether the loop stop condition is met; if the loop stop condition is not met, execute the step of generating a target perception point set based on multiple second perception points.

[0044] In one embodiment, the determination module is specifically used to determine the difference between the maximum number of perception points and the target number; if the difference is greater than or equal to the total number of newly added perception points in a single round, it is determined that the loop stop condition is not met; if the difference is less than the total number of newly added perception points in a single round, it is determined that the loop stop condition is met.

[0045] In one of the embodiments, the determination module is also used to determine the number of additional sensing points based on the difference between the maximum number of sensing points and the target number if the loop stop condition is met; generate multiple third sensing points within the satellite coverage range based on the number of additional sensing points; and detect the interference power at each third sensing point to obtain the interference power value corresponding to each third sensing point.

[0046] In one of the embodiments, the determination module is specifically used to determine the electromagnetic situation information based on the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

[0047] In one of the embodiments, the determination module is specifically used to obtain the position information of each sensing point within the satellite coverage area; interpolate and fit the position information and interference power value of each sensing point based on a preset interpolation method to determine the electromagnetic situation information.

[0048] In a third aspect, the present application further provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in any one of the first aspects is implemented.

[0049] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements any of the methods described in the first aspect above.

[0050] In a fifth aspect, the present application further provides a computer program product, including a computer program, which, when executed by a processor, implements any of the methods described in the first aspect above.

[0051] The above-mentioned electromagnetic situation information determination method, device, satellite and storage medium perform multiple interference power detections within the satellite coverage range of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage range based on the results of the interference power detection. The results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to the own satellite; wherein the i-th interference power detection in the multiple interference power detections includes: obtaining an initial perception point set, the initial perception point set includes multiple first perception points within the satellite coverage range, and obtaining the interference power value at each first perception point; generating multiple second perception points within the satellite coverage range according to the interference power value corresponding to each first perception point, and generating a target perception point set based on the multiple second perception points, the distribution of the multiple second perception points is determined according to the interference power value corresponding to each first perception point, and the target perception point set is the initial perception point set of the i+1-th interference power detection. In this way, the second sensing point is adaptively generated multiple times according to the interference power value corresponding to the first sensing point, and the uplink space spectrum situation of the other party's low-orbit constellation is sensed based on each sensing point to generate real-time space electromagnetic situation information. This can reduce the number of sensing points used under the same sensing accuracy, thereby saving the beam resources involved in sensing as much as possible and shortening the space electromagnetic situation information generation cycle. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the drawings required for use in the embodiments of the present application or related technical descriptions will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0053] Figure 1 It is a flowchart of a method for determining electromagnetic situation information in one embodiment;

[0054] Figure 2 A schematic diagram of the pitch angle and azimuth angle of a sensing point in one embodiment;

[0055] Figure 3 A schematic flow chart of a step of generating a plurality of second sensing points within a satellite coverage range according to interference power values ​​corresponding to each first sensing point in an embodiment;

[0056] Figure 4 A schematic diagram of a flow chart of steps for obtaining multiple sensing point clusters in one embodiment;

[0057] Figure 5 A schematic diagram of a flow chart of a step of determining the number of newly added sensing points corresponding to each sensing point cluster in one embodiment;

[0058] Figure 6 A schematic diagram of a flow chart of a step of determining an allocation ratio corresponding to each sensing point cluster in an embodiment;

[0059] Figure 7 A schematic diagram of a flow chart of a step of generating a target perception point set according to a plurality of second perception points in one embodiment;

[0060] Figure 8 A schematic diagram of a flow chart of a step of determining whether a cycle stop condition is met in one embodiment;

[0061] Fig. 9 It is a flowchart of a method for determining electromagnetic situation information in another embodiment;

[0062] Fig.10 It is a flowchart of a method for determining electromagnetic situation information in another embodiment;

[0063] Fig.11 It is a flowchart of a method for determining electromagnetic situation information in another embodiment;

[0064] Fig.12 A schematic diagram of a constellation configuration in one embodiment;

[0065] Fig.13 Schematic diagram of the elevation angle and azimuth angle coordinates of the first sensing point when i=1 in one embodiment;

[0066] Fig.14 This is a schematic diagram of actual co-frequency interference spatial electromagnetic situation information in an embodiment;

[0067] Fig.15 A schematic diagram of electromagnetic situation information determined according to a spherical random point selection method in one embodiment;

[0068] Fig.16 A schematic diagram of electromagnetic situation information determined according to a spherical uniform point selection method in one embodiment;

[0069] Fig.17 A schematic diagram of electromagnetic situation information determined according to the method in an embodiment of the present application;

[0070] Fig.18It is a curve showing the change of the normalized WMSE of the space electromagnetic situation information obtained by different schemes in one embodiment as the number of sensing points increases;

[0071] Fig.19 is a curve showing a change in the ratio of the number of perception points saved compared to the spherical uniform point selection method with the normalized WMSE in one embodiment;

[0072] Fig. 20 It is a structural block diagram of an electromagnetic situation information determination device in one embodiment;

[0073] Fig.21 FIG. 4 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0074] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0075] A low-orbit communication constellation is a network composed of a large number of low-Earth orbit satellites. Its core features are low latency, global coverage, high bandwidth, and high redundancy of the system. With the advancement of technology, especially the miniaturization of satellites, the reduction of launch costs, and the improvement of network management capabilities, the construction of low-orbit constellations has become more feasible and efficient. Early low-orbit satellite communication systems were mostly small-scale experimental, but now, low-orbit constellations have rapidly expanded in number, coverage, and service capabilities, forming a communication network with global service capabilities. These changes not only improve global Internet connectivity, especially in remote areas that are difficult to cover with traditional ground facilities, but also provide solutions for high-bandwidth and high-speed communication needs. Typical low-orbit satellite communication systems include: Iridium system, OneWeb system, StarLink system, etc.

[0076] Due to the rapid rise of low-orbit satellite communication systems, existing spectrum resources are facing an increasingly severe tension, which is mainly due to the interference risk when multiple low-orbit satellite systems operate in the same frequency band. Due to the surge in the number of low-orbit satellites, there is a lot of overlap in the spectrum resources declared by different operators, especially in the Ku and Ka bands, which are prone to co-channel interference. At the same time, the lack of international unified standards for spectrum allocation and management, especially the lack of supervision and coordination between low-orbit constellation systems, has exacerbated this conflict problem. Therefore, in the uplink co-channel interference scenario where low-orbit constellations coexist, the satellite system of the party dynamically uses the frequency according to local conditions and time, that is, using the sensing beam to actively perceive the spatial spectrum situation of the other party's constellation, thereby using the known spectrum distribution of the other party's system to avoid potential co-channel interference, is a key technical means to achieve efficient coexistence and sharing of spectrum resources between the own system and the other party's system.

[0077] In the process of dynamically sensing the uplink space spectrum situation of the other party's system, the own satellite selects the target frequency band that is more likely to produce co-frequency interference, and uses the on-board sensing beam to sense the co-frequency interference power of the target frequency band in different uplink wave directions. According to the limited sensing information, the predicted co-frequency interference power in all uplink wave directions is obtained by the interpolation fitting method, that is, the current real-time space electromagnetic situation information. The real-time space electromagnetic situation information of the co-frequency interference of the other party's system in the target frequency band can be generated by repeating the above operation. The real-time space electromagnetic situation information can reflect the spatial distribution of the uplink co-frequency interference power of the other party's system in the target frequency band under the current constellation configuration and the beam snapshot of the other party's system, and is a good spectrum situation perception method.

[0078] For the unknown co-frequency interference environment of other systems, the most ideal situation is that our own satellite senses in all possible uplink wave directions. However, when there are too many space sensing points, the sensing process needs to consume a large amount of time slots and beam resources, which will inevitably lead to a decrease in the updating frequency of space electromagnetic situation information and communication efficiency. Therefore, it is urgent to propose a method to save the beam resources involved in sensing as much as possible and shorten the generation cycle of space electromagnetic situation information while ensuring the accuracy of space electromagnetic situation information perception.

[0079] In view of this, the present application provides an electromagnetic situation information determination method that can save space electromagnetic situation information perception resources and shorten the space electromagnetic situation information generation cycle. The electromagnetic situation information determination method provided in the embodiment of the present application, its execution subject can be an electromagnetic situation information determination device, and the electromagnetic situation information determination device can be implemented by software, hardware, or a combination of software and hardware. It can be embedded in or independent of the processor in the computer device in the form of hardware, or it can be stored in the memory of the computer device in the form of software. In the following method embodiments, the execution subject is a computer device as an example for explanation, wherein the computer device can be a server, and the embodiment of the present application does not limit the specific type of the computer device.

[0080] In an exemplary embodiment, Figure 1 As shown, a method for determining electromagnetic situation is provided, comprising the following steps. Among them:

[0081] Multiple interference power detections are performed within the satellite coverage of one's own satellite, and the electromagnetic situation information corresponding to the satellite coverage is determined based on the results of the interference power detection. The results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to one's own satellite.

[0082] The i-th interference power detection in the multiple interference power detections includes:

[0083] Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each first sensing point;

[0084] According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set of the i+1th interference power detection.

[0085] Optionally, a target frequency band is determined based on the spectrum utilization information of one's own satellite system and the other party's satellite system, wherein the target frequency band may be a user frequency band where the spectrum utilization information of one's own satellite system and the transmitting satellite system overlaps most seriously. Interference power detection will be performed based on the target frequency band to construct space electromagnetic situation information.

[0086] Optionally, in the initial state, that is, when i is equal to 1, the multiple first perception points included in the initial perception point set can be selected within the satellite coverage according to a spherical random sampling method or a spherical uniform sampling method.

[0087] Optional, such as Figure 2 As shown, the satellite coverage range can be determined based on the maximum elevation angle of the target satellite. It is used to describe the maximum pointing angle of the target satellite antenna in the vertical direction. It can be expressed by the spatial rectangular coordinate system. Indicates the location information of the first perception point, and can also be converted into elevation angle and azimuth The two-dimensional angle coordinates ( Indicates the position information of the first sensing point, where the pitch angle is used to describe the vertical position of the first sensing point relative to the target satellite, and the azimuth angle is used to describe the horizontal position of the first sensing point relative to the target satellite. Figure 2 middle is the minimum elevation angle of the user equipment, is the radius of the Earth, is the satellite orbit altitude of one's own satellite.

[0088] Exemplarily, in the initial state, the spherical uniform sampling method is used at a maximum pitch angle of Determining multiple first sensing points within the satellite coverage area can be done by the following methods:

[0089] First, get the preset initial pitch angle interval number , and then The maximum pitch angle is evenly divided by the interval of Pitch angle, at a pitch angle of The first sensing point is selected by evenly dividing the target angle interval in the corresponding azimuth interval sampling. The target angle interval can be selected according to Determined, understandably, through right By weighting, multiple first perception points evenly distributed within the satellite coverage can be obtained.

[0090] Optionally, the interference power value at each first sensing point may be a co-channel interference power value in the uplink incoming wave direction corresponding to the two-dimensional angular coordinates of all first sensing points sensed by the satellite sensing beam. The specific interference power value may be calculated by the following formula:

[0091]

[0092] in, is the co-channel interference power (in dBw / Hz), is the user equipment transmit power (in dBw / Hz), is the antenna gain of the user equipment of the other party’s satellite system, is the propagation path loss, is the antenna gain of the satellite receiver.

[0093] Among them, when the user equipment adopts an omnidirectional antenna, is a constant of 0, The size of is determined by the angle between the sensing beam vector of the satellite and the vector pointing from the satellite to the user equipment of the other satellite system and the antenna pattern. The propagation path calculation method is:

[0094]

[0095] in, The uplink frequency band of the user equipment, i.e. the target frequency band, in MHz. is the distance between the user equipment of the other satellite system and the own satellite, in km, is the polarization loss, which can be taken as a constant of 3.

[0096] Optionally, the following describes generating multiple second sensing points within the satellite coverage according to the interference power values ​​corresponding to the first sensing points:

[0097] In one possible implementation, the first perception points may be sorted according to the size of the interference power value, and then more second perception points may be generated around the first perception points with larger interference power values, and fewer second perception points may be generated around the first perception points with smaller interference power values. The specific method for generating the second perception points around a first perception point may be to determine the generation range of the corresponding second perception points with the first perception point as the center, and then randomly or uniformly generate multiple second perception points within the generation range.

[0098] In another possible implementation manner, each first perception point may be clustered to obtain a plurality of perception point clusters, and then a corresponding second perception point may be generated in each perception point cluster.

[0099] Optionally, it can be understood that generating the second perception point is to determine the location information of the second perception point, and then detecting the interference power at each second perception point to obtain the interference power value of each second perception point. The process of determining the interference power value at each perception point in the embodiment of the present application will not be repeated.

[0100] Optionally, generating a target perception point set based on multiple second perception points may be to use multiple first perception points and multiple second perception points as the target perception point set; or, to use multiple second perception points as the target perception point set; or, to select a target number of perception points from multiple second perception points and multiple first perception points as the target perception point set.

[0101] It can be understood that multiple second perception points are used as target perception point sets or a target number of perception points are selected from multiple second perception points and multiple first perception points as target perception point sets, and the target perception point set is used as the initial perception point set for the next interference power detection. In this way, except for the first time, the number of first perception points used each time is the same, thereby ensuring that the computational complexity does not gradually increase with each round of iteration.

[0102] Optionally, when i>1, when the target perception point set is used as the initial perception point set for the next interference power detection, the interference power values ​​of each first perception point in the initial perception point set have been detected, and the interference power value of each perception point can be directly obtained, or the interference power of each first perception point can be detected again to obtain the corresponding interference power value. The embodiments of the present application are not limited to this.

[0103] Optionally, a stop condition for cyclically executing interference power detection is described:

[0104] In a possible implementation manner, a preset number of executions may be obtained, and when the number of executions of the interference power detection is equal to the preset number of executions, the interference power detection is no longer executed.

[0105] In another possible implementation method, the maximum number of perception points can be obtained, and the target number of perception points can be determined based on the number of perception points generated in previous interference power detection processes and the number of perception points in the initial state. When the target number of perception points is greater than or equal to the maximum number of perception points, or when the target number of perception points and the maximum number of perception points are less than a preset threshold, interference power detection is no longer performed.

[0106] The above-mentioned electromagnetic situation information determination method performs multiple interference power detections within the satellite coverage range of its own satellite, and determines the electromagnetic situation information corresponding to the satellite coverage range based on the results of the interference power detection. The results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to its own satellite. Among them, the i-th interference power detection in the multiple interference power detections includes: obtaining an initial perception point set, the initial perception point set includes multiple first perception points within the satellite coverage range, and obtaining the interference power value at each first perception point; generating multiple second perception points within the satellite coverage range based on the interference power value corresponding to each first perception point, and generating a target perception point set based on the multiple second perception points, the distribution of the multiple second perception points is determined based on the interference power value corresponding to each first perception point, and the target perception point set is the initial perception point set of the i+1-th interference power detection. In this way, the second sensing point is adaptively generated multiple times according to the interference power value corresponding to the first sensing point, and the uplink space spectrum situation of the other party's low-orbit constellation is sensed based on each sensing point to generate real-time space electromagnetic situation information. This can reduce the number of sensing points used under the same sensing accuracy, thereby saving the beam resources involved in sensing as much as possible and shortening the space electromagnetic situation information generation cycle.

[0107] In an exemplary embodiment, optionally, Figure 3 As shown, generating multiple second sensing points within the satellite coverage according to the interference power values ​​corresponding to the first sensing points includes the following steps 301 to 303. Among them:

[0108] S301: Clustering a plurality of first perception points according to an interference power value corresponding to each first perception point to obtain a plurality of perception point clusters.

[0109] In one possible implementation, clustering is performed on multiple first perception points according to the interference power value corresponding to each first perception point to obtain multiple perception point clusters, which may include: obtaining a preset number of clusters, and then clustering the multiple first perception points into perception point clusters of a corresponding number based on a clustering algorithm according to the interference power value corresponding to each first perception. It can be understood that, for a perception point cluster, the interference power values ​​of each perception point in the cluster are highly similar.

[0110] The clustering algorithm may be a K-means algorithm, an optimization algorithm of the K-means algorithm (eg, a K-means++ algorithm), a hierarchical clustering algorithm, or a density-based clustering algorithm.

[0111] In another possible implementation, Figure 4 As shown, according to the interference power value corresponding to each first perception point, clustering is performed on multiple first perception points to obtain multiple perception point clusters, and the following steps 401 to 403 may also be included. Among them:

[0112] S401: Clustering multiple first perception points according to interference power values ​​corresponding to the first perception points to obtain an initial perception point cluster.

[0113] Optionally, a first preset cluster number K may be obtained, and clustering may be performed based on a clustering algorithm according to interference power values ​​corresponding to the first perception points to obtain a plurality of initial perception point clusters, where the number of the initial perception point clusters is K.

[0114] S402: For each initial perception point cluster, clustering processing is performed on the initial perception point cluster according to the position information of each first perception point in the initial perception point cluster to obtain a plurality of sub-clusters included in the initial perception point cluster.

[0115] Optionally, for each initial perception point cluster, clustering is performed again according to the position information of each first perception point included in the cluster, and the initial perception point cluster is clustered into multiple sub-clusters.

[0116] Optionally, a second preset number of clusters L may be obtained, and each first perception point in an initial perception point cluster may be clustered into L subclusters. It is understandable that if the number of first perception points in an initial perception point cluster is less than the second preset number of clusters L, each first perception point is a separate subcluster.

[0117] Optionally, the clustering process may be performed on each initial perception point cluster in parallel, or the clustering process may be performed on each initial perception point cluster in sequence according to a preset order.

[0118] Optionally, in the clustering process, the distance calculation method between any two first perception points can be Euclidean distance, Manhattan distance, or the like. When the position information of each first perception point in the initial perception point cluster is calculated using the two-dimensional angular coordinates of each first perception point When expressed, the distance between any two first perception points can be determined by the distance formula between two points in polar coordinates. The calculation formula is as follows:

[0119]

[0120] in, is the distance between any two first perception points, is the two-dimensional angle coordinate of a first perception point, is the two-dimensional angular coordinate of another first perception point.

[0121] S403: taking the subclusters included in each initial perception point cluster as a perception point cluster.

[0122] Optionally, through the above clustering process, at most KL sub-clusters can be obtained, that is, at most KL perception point clusters can be obtained.

[0123] In the above, according to the interference power value corresponding to each first perception point, a clustering process is performed on the multiple first perception points to obtain an initial perception point cluster. For each initial perception point cluster, the initial perception point cluster is clustered according to the position information of each first perception point in the initial perception point cluster to obtain a plurality of sub-clusters included in the initial perception point cluster. The sub-clusters included in each initial perception point cluster are used as perception point clusters. In this way, the multiple first perception points can be finely divided to obtain a plurality of perception point clusters. Subsequently, second perception points are generated based on the plurality of perception point clusters to more accurately perceive the change in interference power.

[0124] S302: Determine the number of newly added sensing points corresponding to each sensing point cluster according to the clustering power value corresponding to each sensing point cluster.

[0125] Optionally, the cluster power value corresponding to each sensing point cluster may be the interference power value corresponding to the cluster center of each sensing point cluster. The cluster center may not necessarily be the actual power value of a sensing point in the sensing point cluster, but is calculated based on certain rules.

[0126] Optionally, for subclusters in the same initial perception point cluster, the clustering power value corresponding to the initial perception point cluster may be used as the clustering power value of each subcluster.

[0127] Exemplarily, the mean of the interference power values ​​of all the perception points in the perception point cluster may be used as the interference power value of the perception point at the center of the cluster, or the median of the interference power values ​​of all the perception points in the perception point cluster may be used as the interference power value of the perception point at the center of the cluster.

[0128] Optionally, the number of newly added perception points corresponding to the perception point cluster is positively correlated with the clustering power value corresponding to the perception point cluster. The larger the clustering power value corresponding to the perception point cluster, the greater the number of newly added perception points in the corresponding perception point cluster. In this way, perception can be performed through multiple perception points in the perception area corresponding to the perception point cluster, which can improve the perception accuracy of the perception area.

[0129] S303: Generate second perception points in the perception areas corresponding to the perception point clusters according to the number of newly added perception points corresponding to the perception point clusters.

[0130] The sensing area corresponding to the sensing point cluster is the area where the sensing points included in the sensing point cluster are located within the satellite coverage area.

[0131] Optionally, the perception area corresponding to the perception point cluster can be determined by the cluster position corresponding to the perception point cluster and the position information of each first perception point included in the perception point cluster. The cluster position corresponding to the perception point cluster can be the position information corresponding to the cluster center of the perception point cluster.

[0132] Optionally, a sensing region can be ,in, and They are respectively the minimum and maximum values ​​of the Euclidean distance between the position information corresponding to the first perception point in the perception point cluster and the cluster center of the perception point cluster in the rectangular coordinate system.

[0133] Optionally, generating the second perception points in the perception area corresponding to each perception point cluster may be performed by randomly generating a corresponding number of newly added perception points in the corresponding perception area, or by uniformly generating a corresponding number of newly added perception points in the corresponding perception area.

[0134] For example, the number of newly added sensing points corresponding to the sensing point cluster is For example, in Randomly select r in the range of Randomly select from the range , and As the two-dimensional angle coordinates of the second perception point, a total of A second perception point.

[0135] In the above, according to the interference power value corresponding to each first perception point, multiple first perception points are clustered to obtain multiple perception point clusters; according to the cluster power value corresponding to each perception point cluster, the number of newly added perception points corresponding to each perception point cluster is determined; according to the number of newly added perception points corresponding to each perception point cluster, second perception points are generated in the perception area corresponding to each perception point cluster, and the second perception points can be adaptively generated according to the interference power value and perception can be performed, which can save space electromagnetic situation information perception resources as much as possible and shorten the space electromagnetic situation information generation cycle while meeting the perception accuracy.

[0136] In an exemplary embodiment, optionally, Figure 5 As shown, according to the clustering power value corresponding to each sensing point cluster, determining the number of newly added sensing points corresponding to each sensing point cluster includes the following steps 501 to 502. Among them:

[0137] S501, obtaining the total number of newly added sensing points in a single round, and determining the allocation ratio corresponding to each sensing point cluster according to the clustering power value corresponding to each sensing point cluster.

[0138] Optionally, the total number of newly added sensing points in a single round can be preset. It can be understood that each time multiple second sensing points are generated within the satellite coverage range according to the interference power values ​​corresponding to each first sensing point, the number of second sensing points generated is fixed and known.

[0139] In a possible implementation manner, the allocation ratio corresponding to each sensing point cluster may be determined according to the sum of the cluster power values ​​corresponding to each sensing point cluster and the cluster power values ​​corresponding to each sensing point cluster.

[0140] In another possible implementation method, the allocation ratio corresponding to each perception point cluster can be determined according to the cluster power value corresponding to each perception point cluster and the number of first perception points included in each perception point cluster, wherein, for a perception point cluster, its allocation ratio is positively correlated with the cluster power value, and negatively correlated with the number of included first perception points, that is to say, the higher the cluster power value, the higher the allocation ratio, and the more the number of included first perception points, the lower the allocation ratio, and the second perception point can be adaptively generated according to the interference power value, and it can avoid generating too many perception points in a perception area, resulting in duplication of perception points.

[0141] In another possible implementation, considering the uniformity of the generation of the second perception points, for each perception point cluster, the number of included first perception points may be weighted based on the clustering position of the perception point cluster, and then the allocation ratio corresponding to the perception point cluster is determined, such as Figure 6 As shown, according to the clustering power value corresponding to each sensing point cluster, determining the allocation ratio corresponding to each sensing point cluster includes the following steps 601 to 602. Among them:

[0142] S601, obtaining the cluster position corresponding to each perception point cluster.

[0143] Optionally, the cluster position corresponding to each perception point cluster may be represented by a two-dimensional angular coordinate.

[0144] S602: For each sensing point cluster, determine an allocation ratio according to a cluster power value corresponding to the sensing point cluster, a cluster position, and the number of first sensing points included in the sensing point cluster.

[0145] Optionally, for each sensing point cluster, the corresponding allocation ratio is determined according to the cluster power value corresponding to the sensing point cluster, the cluster position and the number of first sensing points included in the sensing point cluster, and the pitch angle in the cluster position can be determined. The sine value of weights the number of included first perception points.

[0146] Optionally, taking the case where the sensing point cluster is a sub-cluster included in each initial sensing point cluster, the allocation ratio may be calculated by the following formula:

[0147]

[0148] in, For the The initial perception point cluster includes the The allocation ratio corresponding to the sub-clusters is is the clustering power value corresponding to the subcluster, is the number of first perception points included in the sub-cluster, is the pitch angle in the cluster position of this sub-cluster.

[0149] S502: Allocate the number of newly added sensing points to each sensing point cluster from the total number according to the allocation ratio corresponding to each sensing point cluster.

[0150] The sum of the number of newly added sensing points allocated to each sensing point cluster is equal to the total number.

[0151] Optionally, the number of newly added sensing points allocated to each sensing point cluster may be determined based on the product of the allocation ratio corresponding to each sensing point cluster and the total number.

[0152] Optionally, it can be expressed by the following formula:

[0153]

[0154] in, For the The initial perception point cluster includes the The number of newly added sensing points corresponding to each sub-cluster.

[0155] In the above, the total number of newly added perception points in a single round is obtained, and the allocation ratio corresponding to each perception point cluster is determined according to the clustering power value corresponding to each perception point cluster. According to the allocation ratio corresponding to each perception point cluster, the number of newly added perception points is allocated to each perception point cluster from the total number. New perception points can be adaptively generated according to the clustering power value of each perception point cluster, thereby saving space electromagnetic situation information perception resources as much as possible and shortening the space electromagnetic situation information generation cycle under the premise of the same perception accuracy.

[0156] In an exemplary embodiment, optionally, Figure 7 As shown, generating a target perception point set according to a plurality of second perception points includes the following steps 701 to 703. Among them:

[0157] S701, obtaining a maximum number of sensing points, and obtaining a target number of sensing points within a satellite coverage area.

[0158] The target number includes the number of sensing points generated during each interference power detection process and the number of sensing points in the initial state.

[0159] Optionally, the maximum number of perception points may be the number of perception points that meets the perception accuracy requirement, and may be determined based on historical perception.

[0160] Optionally, the target number can be determined by the sum of the number of first perception points in the initial perception point set in the initial state and the number of second perception points generated in previous interference power detection processes.

[0161] Optionally, there may be duplicate points between the second perception points generated in the interference power detection process, or there may be duplicate points between the generated second perception points and the first perception points in the initial perception point set. When determining the target number, deduplication processing is required.

[0162] S702: Determine whether the loop stop condition is met based on the maximum number of perception points and the target number.

[0163] In one possible implementation, the maximum number of perception points and the target number are compared. If the target number is greater than or equal to the maximum number of perception points, it is determined that the loop stop condition is met; if the target number is less than the maximum number of perception points, it is confirmed that the loop stop condition is not met.

[0164] In another possible implementation, Figure 8 As shown, according to the maximum number of sensing points and the target number, determining whether the loop stop condition is met may include the following steps 801 to 803. Among them:

[0165] S801, determining the difference between the maximum number of perception points and the target number.

[0166] Optionally, in an embodiment of the present application, the target number is always less than or equal to the maximum number of perception points. Therefore, the difference between the maximum number of perception points and the target number can be determined by subtracting the target number from the maximum number of perception points.

[0167] S802: If the difference is greater than or equal to the total number of newly added perception points in a single round, it is determined that the loop stop condition is not met.

[0168] Optionally, since the number of second perception points generated each time is fixed, when the difference between the target number and the maximum number of perception points is greater than or equal to the total number of newly added perception points in a single round, it can be determined that the loop stop condition is not met, and a target perception point set can be generated based on multiple second perception points, and the next interference power detection step can be performed to increase the target number.

[0169] S803: If the difference is less than the total number of newly added sensing points in a single round, it is determined that the loop stop condition is met.

[0170] Optionally, if the difference is less than the total number of newly added perception points in a single round, if the loop is not stopped, then the target number obtained after the second perception point is generated next time exceeds the maximum value of the perception points. This does not meet the requirement of reducing the number of perception points used as much as possible under the same perception accuracy. Therefore, when the difference is less than the total number of newly added perception points in a single round, it is determined that the loop stop condition is met and the loop is stopped.

[0171] S703: If the loop stop condition is not met, the step of generating a target perception point set according to multiple second perception points is executed.

[0172] The above-mentioned method obtains the maximum value of the number of perception points and the target number of perception points included in the satellite coverage area, and determines whether the loop stop condition is met based on the maximum value of the number of perception points and the target number. If the loop stop condition is not met, the step of generating a target perception point set based on multiple second perception points is executed. Multiple interference power detections can be performed by determining the maximum value of the number of perception points and the target number, thereby obtaining multiple perception points to accurately determine the degree of interference distribution within the satellite coverage area.

[0173] In an exemplary embodiment, optionally, Fig. 9 As shown, the method further includes the following steps 901 to 903, wherein:

[0174] S901: If the loop stop condition is met, the number of additional perception points is determined according to the difference between the maximum number of perception points and the target number.

[0175] Optionally, if the loop stop condition is met, the maximum number of perception points is greater than or equal to the target number. When the maximum number of perception points is greater than the target number, the number of additional perception points can be determined based on the difference.

[0176] S902: Generate multiple third sensing points within the satellite coverage area according to the number of additional sensing points.

[0177] Optionally, a plurality of third perception points may be generated within the entire satellite coverage area by using a spherical uniform sampling method or a spherical random sampling method, and the number of the third perception points is the same as the number of the supplementary perception points.

[0178] Optionally, the maximum pitch angle can be discretized into N pitch angles, and the target pitch angle is determined according to the sine value of each pitch angle. Within the satellite coverage range corresponding to the target pitch angle, the azimuth angle is The pitch angles with the highest probability are selected as the target pitch angle.

[0179]

[0180] in, is the nth pitch angle, is the selection probability corresponding to the nth pitch angle.

[0181] S903: Detect the interference power at each third perception point to obtain the interference power value corresponding to each third perception point.

[0182] Optionally, when detecting the interference power at each third sensing point, the satellite sensing beam is also used to sense the co-frequency interference power value in the uplink direction corresponding to the two-dimensional angular coordinates of all third sensing points. This is the same as the interference power detection method for the first sensing point, and the embodiments of the present application will not elaborate on this.

[0183] If the above-mentioned loop stop condition is met, the number of additional sensing points is determined according to the difference between the maximum number of sensing points and the target number, and multiple third sensing points are generated within the satellite coverage range according to the number of additional sensing points. The interference power at each third sensing point is detected to obtain the interference power value corresponding to each third sensing point. When the loop stop condition is met and the target number is less than the maximum number of sensing points, the third sensing point is generated and the interference power value at the third sensing point is determined, so that the electromagnetic situation information can be determined according to the interference power value of each sensing point within the satellite coverage range.

[0184] In an exemplary embodiment, optionally, determining electromagnetic situation information corresponding to a satellite coverage range according to a result of interference power detection includes:

[0185] The electromagnetic situation information is determined based on the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

[0186] In one possible implementation, Fig.10 As shown, Fig.10 As shown, according to the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection, the electromagnetic situation information is determined, including the following steps 1001 to 1002. Among them:

[0187] S1001, obtaining the location information of each sensing point within the satellite coverage area.

[0188] Optionally, the location information and interference power value corresponding to each generated sensing point are stored, and when the number of sensing points within the satellite coverage area is equal to the maximum number of sensing points, the location information and interference power value of each sensing point are obtained.

[0189] S1002: Perform interpolation fitting on the position information and interference power value of each sensing point based on a preset interpolation method to determine electromagnetic situation information.

[0190] Optionally, the preset interpolation method may include polynomial interpolation, spline interpolation, Kriging method, etc., and the spline interpolation may include cubic spline interpolation or quadratic spline interpolation. The embodiments of the present application do not limit the preset interpolation method.

[0191] Optionally, through a preset interpolation method, perception can be performed through a limited number of sensing points to determine the electromagnetic situation information corresponding to the satellite coverage range.

[0192] In another possible implementation, the electromagnetic situation information within the entire satellite coverage area may be determined based on a preset statistical method according to the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

[0193] Optionally, the preset statistical method may include linear regression, multivariate linear regression, nonlinear regression, etc., to determine the interference power value at the sensing point not selected within the satellite coverage.

[0194] By determining the electromagnetic situation information based on the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection, the same-frequency interference of the other party's constellation system to one's own satellite system can be avoided based on the electromagnetic situation information, thereby achieving same-frequency coexistence.

[0195] As an optional implementation, Fig.11 As shown, the electromagnetic situation information determination method provided in the embodiment of the present application may include the following specific steps:

[0196] S1101, obtaining an initial perception point set.

[0197] The initial sensing point set includes multiple first sensing points within the satellite coverage area.

[0198] S1102: Obtain interference power values ​​at each first sensing point.

[0199] S1103: Clustering the multiple first perception points according to the interference power value corresponding to each first perception point to obtain an initial perception point cluster.

[0200] S1104: For each initial perception point cluster, clustering processing is performed on the initial perception point cluster according to the position information of each first perception point in the initial perception point cluster to obtain a plurality of sub-clusters included in the initial perception point cluster.

[0201] S1105: taking the subclusters included in each initial perception point cluster as a perception point cluster.

[0202] S1006: Obtain the cluster power value and cluster position corresponding to each sensing point cluster.

[0203] S1107: For each sensing point cluster, determine an allocation ratio according to a cluster power value corresponding to the sensing point cluster, a cluster position, and the number of first sensing points included in the sensing point cluster.

[0204] S1108, obtaining the total number of newly added sensing points in a single round.

[0205] S1109: Allocate the number of newly added sensing points to each sensing point cluster from the total number according to the allocation ratio corresponding to each sensing point cluster.

[0206] The sum of the number of newly added sensing points allocated to each sensing point cluster is equal to the total number.

[0207] S1110: Generate second perception points in the perception areas corresponding to the perception point clusters according to the number of newly added perception points corresponding to the perception point clusters.

[0208] The sensing area corresponding to the sensing point cluster is the area where the sensing points included in the sensing point cluster are located within the satellite coverage area.

[0209] S1111, generating a target perception point set according to multiple second perception points.

[0210] S1112, obtaining the maximum number of sensing points, and obtaining the target number of sensing points included in the satellite coverage.

[0211] The target number includes the number of sensing points generated during each interference power detection process and the number of sensing points in the initial state.

[0212] S1113, determining the difference between the maximum number of perception points and the target number.

[0213] S1114: If the difference is greater than or equal to the total number of newly added perception points in a single round, the target perception point set is used as the initial perception point set, and steps 1101 to 1113 are executed.

[0214] S1115: If the difference is less than the total number of newly added perception points in a single round, the number of additional perception points is determined according to the difference between the maximum number of perception points and the target number.

[0215] S1116, generating multiple third sensing points within the satellite coverage according to the number of the additional sensing points.

[0216] S1117, detecting the interference power at each third perception point to obtain the interference power value corresponding to each third perception point.

[0217] S1118, obtaining the position information of each sensing point within the satellite coverage area determined by each interference power detection.

[0218] S1119, interpolating and fitting the position information and interference power value of each sensing point based on a preset interpolation method to determine the electromagnetic situation information.

[0219] Exemplarily, for a low-orbit non-geostationary orbit satellite constellation system of another party, the electromagnetic situation information is determined according to the method provided in the embodiment of the present application, wherein the typical constellation system Set-1 specified by 3GPP 5G NTN is selected, and its basic configuration parameters are: orbital altitude 600km, orbital inclination 87.5 degrees, eccentricity 0, constellation configuration Star constellation, number of orbital planes is 17, there are 30 satellites on each orbital plane, phase factor is 8, and the own system is a satellite that is subject to co-frequency interference from the other party's constellation system. Since the algorithm execution process is consistent for any co-frequency interference environment, the own system is located near one of the satellites in the Set-1 constellation and is subject to uplink co-frequency interference from the nearest 19 Set-1 satellites. The constellation configuration is as follows Fig.12 As shown in the figure, the dot in the center of the figure represents its own satellite, the cross represents the other satellite that interferes with its own satellite, and the dots on the edge represent the other satellite that does not interfere with its own satellite. The communication parameters of the constellation system and user equipment (UE) also refer to the technical specifications in 3GPP 5G NTN. The user equipment can be various user terminals that can access the communication network. The communication parameters of the satellite system and the communication parameters of the user terminal are shown in Table 1 and Table 2 respectively.

[0220] Table 1

[0221]

[0222] Table 2

[0223]

[0224] First, you can select any UE communication frequency band (1MHz) as the target frequency band, according to the UE minimum elevation angle Use the following formula to calculate the maximum pitch angle of your satellite ,exist When the angle is 30°, It is 52.3316°.

[0225]

[0226] in, is the radius of the Earth, is the satellite orbit altitude.

[0227] Then, at the maximum pitch angle The initial sensing point is generated by spherical uniform sampling method within the satellite coverage area, that is, the first sensing point when i=1, and the selection As the initial uniform perception pitch angle interval number, The pitch angle interval is divided equally, and the pitch angle is The corresponding azimuth interval is The angle interval is The number of the first sensing points is determined to be 705. The elevation angle and azimuth angle coordinates of the first sensing point are as follows: Fig.13 As shown, the first sensing points are evenly distributed within the satellite coverage area, and the maximum number of sensing points is selected. , (where the spherical uniform sampling method is used to determine The number of pitch angle intervals when there are 40 sensing points in a single round is selected, the total number of newly added sensing points in a single round is selected as B=500, the first preset number of clusters K=10, the second preset number of clusters L=30, and the co-frequency interference power value of the target user frequency band is sensed according to the uplink beam sensing direction corresponding to the two-dimensional angular coordinates of the initial sensing point. In this example, the interference power value of the first sensing point is calculated according to the parameters shown in Tables 1 and 2.

[0228] Then, the K-means++ clustering algorithm is used to cluster the first perception point (the initial perception point set in the first iteration) into K=10 initial perception point clusters according to the interference power value. The power value corresponding to the cluster center of each cluster is recorded as , the number of perception points contained in each cluster is recorded as For example, after clustering, the historical perception points corresponding to the cluster center P1=-219.87dBW have a total of M1=57 points, and the co-channel interference power values ​​corresponding to these first perception points are all close to -219.87dBW.

[0229] Then, the K-means++ clustering algorithm is used to divide the first perception points in each initial perception point cluster into L=30 sub-clusters according to the two-dimensional angular coordinates corresponding to the first perception points in K=10 initial perception point clusters. If the number of perception points in the cluster is less than L, each first perception point is divided into a sub-cluster separately. When clustering, the distance between any two first perception points adopts the distance between the two points in polar coordinates. At this time, there are a maximum of KL=300 perception point clusters, and the two-dimensional angle coordinates corresponding to the cluster center of each perception point cluster are marked as , the number of perception points contained in each perception point cluster is recorded as For example, the perception point cluster <2,15> corresponds to the 15th subcluster in the second initial perception point cluster, which has a total of M 2,15 =33 first perception points, the corresponding cluster power value is the power value of the parent cluster center, which is P2=-223.24dBW, and the corresponding cluster position is the position of the sub-cluster center. is (32.15°,145.31°).

[0230] Then, the number of second perception points generated is allocated to each perception point cluster from the total number of newly added perception points B = 500 in a single round, and the perception point cluster<k,l> The corresponding number of newly added perception points is calculated by the following formula:

[0231]

[0232] Among them, the higher the co-frequency interference power value of the perception point cluster, the higher the allocation ratio, the higher the number of first perception points contained in the perception point cluster, the lower the allocation ratio. At the same time, considering the uniformity of spherical sampling, the number of first perception points is weighted by the sine value of the pitch angle of the two-dimensional angular coordinate of the center of the perception point cluster.

[0233] Then, the perception area of ​​each perception point cluster is delineated and the second perception point is randomly generated according to the number of newly added perception points of each perception point cluster.<k,l> ,exist Randomly select r in the range of Randomly select from the range , and As the two-dimensional angular coordinates of the second perception point, each perception point cluster generates The two-dimensional angular coordinates corresponding to the second perception points of all perception point clusters are taken as the target perception point set, and also as the initial perception point set for the next round of iteration. Except for the first round, the number of initial perception points used in each round of iteration is the same, thereby ensuring that the algorithm complexity does not gradually increase with each round of iteration.

[0234] Then, the sensing beam sequentially senses the co-frequency interference power of the other system in the direction of the second sensing point according to the target sensing point set, thereby completing a round of iteration. The number of target sensing points increases with the number of iterations. After that, or when the difference between the maximum number of perception points and the total number of newly added perception points in a single round is less than 500, the iteration is stopped, and for the third perception points that are the difference between the duplicate points generated in the second perception point generation process and the difference after the iteration, multiple third perception points are determined according to the spherical uniform sampling method.

[0235] Then, the co-channel interference power values ​​and two-dimensional angular coordinate information of all sensing points within the satellite coverage are interpolated by cubic spline with an elevation angle interval of 100 and an azimuth angle interval of 543 (corresponding to the spherical uniform sampling method with an elevation angle interval of 100) to generate real-time spatial electromagnetic situation information of the co-channel interference of other systems.

[0236] Finally, the method in the embodiment of the present application is compared with the spherical random point method and the spherical uniform point method to verify the effectiveness of the method in the embodiment of the present application. The spherical random point method selects the point according to the pitch angle. is the pitch angle of the sensing point, and the azimuth is The spherical uniform point selection method is the same as the method for selecting the first perception point in the initial perception point set when i=1 in the embodiment of the present application. The total number of perception points is increased by increasing the number of pitch angle intervals. Under a random other party system beam snapshot, the frequency bands of target users in different schemes are The contour map of the interpolation fitting result of the point space electromagnetic situation information and the corresponding two-dimensional angle coordinate distribution of the sensing point are shown in the figure below: Figures 14 to 17 As shown, Fig.14 It is the actual co-frequency interference space electromagnetic situation information. Fig.15 The electromagnetic situation information is determined by the spherical random point method. Fig.16 The electromagnetic situation information is determined by the spherical uniform point selection method. Fig.17 It is the electromagnetic situation information determined according to the method in the embodiment of the present application. According to the simulation results, it can be seen that, unlike the spherical random point selection method and the spherical uniform point selection method, the two-dimensional angular coordinates of the sensing points obtained in the embodiment of the present application are more concentrated in the uplink wave direction where the co-frequency interference power changes dramatically.

[0237] Optionally, in order to quantitatively compare the perceptual accuracy performance differences of different schemes, a weighted mean square error (WMSE) is defined, and its specific value is calculated using the following formula:

[0238]

[0239] Where N is the number of target sensing points, is the actual co-channel interference power in the uplink direction of the nth sensing point, It is the co-channel interference power predicted by air electromagnetic situation information in the uplink direction of the n-th sensing point.

[0240] Optionally, simulation verifies the change curve of the normalized WMSE of the spatial electromagnetic situation information obtained by three different schemes under 100 other party system beam snapshots as the number of sensing points increases, and verifies that the sensing point set selects a spherical uniform sampling sensing point set when the number of pitch angle intervals is 100. The simulation results are as follows: Fig.18 As shown, 1801 is the change curve corresponding to the method in the embodiment of the present application, 1802 is the change curve corresponding to the spherical uniform point selection method, and 1803 is the change curve corresponding to the spherical random point selection method.

[0241] Depend on Fig.18 It can be seen from the simulation results that when the WMSE performance of the spatial electromagnetic situation information is the same, the number of perception points used by the method in the embodiment of the present application is less than that of the other two comparison schemes. This shows that the method in the embodiment of the present application uses fewer perception points under the premise of the same perception accuracy, thereby achieving the purpose of saving spatial electromagnetic situation information perception resources as much as possible and shortening the spatial electromagnetic situation information generation cycle.

[0242] Optionally, when the perception accuracy is the same, the curve of the ratio of perception points saved by the method in the embodiment of the present application compared with the spherical uniform point selection method as the normalized WMSE is as follows: Fig.19 As shown, it can be seen that when the perception accuracy is the same, the method in the embodiment of the present application can save up to 55.2% of the perception points compared to the spherical uniform point selection method. Therefore, with the method in the embodiment of the present application, it is possible to save space electromagnetic situation information perception resources as much as possible and shorten the space electromagnetic situation information generation cycle under the premise of ensuring perception accuracy, thereby constructing effective satellite uplink space electromagnetic situation information through interpolation fitting to avoid the co-frequency interference of the other party's constellation system to the own satellite system and achieve co-frequency coexistence.

[0243] It should be understood that, although the various steps in the flowcharts involved in the above-mentioned embodiments are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear explanation in this article, the execution of these steps does not have a strict order restriction, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-mentioned embodiments can include multiple steps or multiple stages, and these steps or stages are not necessarily executed 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 part of the steps or stages in other steps.

[0244] Based on the same inventive concept, the embodiment of the present application also provides an electromagnetic situation information determination device for implementing the electromagnetic situation information determination method involved above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations in one or more electromagnetic situation information determination device embodiments provided below can refer to the limitations of the electromagnetic situation information determination method above, and will not be repeated here.

[0245] In an exemplary embodiment, Fig. 20 As shown, an electromagnetic situation information determination device 2000 is provided, comprising: a determination module 2001, wherein:

[0246] The determination module 2001 is used to perform multiple interference power detections within the satellite coverage of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage according to the results of the interference power detection, and the results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to the own satellite;

[0247] The i-th interference power detection in the multiple interference power detections includes:

[0248] Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each first sensing point;

[0249] According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set of the i+1th interference power detection.

[0250] In an exemplary embodiment, the determination module 2001 is specifically used to cluster multiple first perception points according to the interference power value corresponding to each first perception point to obtain multiple perception point clusters; determine the number of newly added perception points corresponding to each perception point cluster according to the clustering power value corresponding to each perception point cluster; and generate second perception points in the perception area corresponding to each perception point cluster according to the number of newly added perception points corresponding to each perception point cluster, wherein the perception area corresponding to the perception point cluster is the area where the perception points included in the perception point cluster are located within the satellite coverage range.

[0251] In an exemplary embodiment, the determination module 2001 is specifically used to cluster multiple first perception points according to the interference power value corresponding to each first perception point to obtain an initial perception point cluster; for each initial perception point cluster, the initial perception point cluster is clustered according to the position information of each first perception point in the initial perception point cluster to obtain multiple sub-clusters included in the initial perception point cluster; and the sub-clusters included in each initial perception point cluster are used as perception point clusters.

[0252] In an exemplary embodiment, the number of newly added sensing points corresponding to the sensing point cluster is positively correlated with the clustering power value corresponding to the sensing point cluster.

[0253] In an exemplary embodiment, the determination module 2001 is specifically used to obtain the total number of newly added perception points in a single round, and determine the allocation ratio corresponding to each perception point cluster according to the clustering power value corresponding to each perception point cluster; according to the allocation ratio corresponding to each perception point cluster, allocate the number of newly added perception points to each perception point cluster from the total number, wherein the sum of the number of newly added perception points allocated to each perception point cluster is equal to the total number.

[0254] In an exemplary embodiment, the determination module 2001 is specifically used to obtain the cluster position corresponding to each perception point cluster; for each perception point cluster, the allocation ratio is determined according to the cluster power value corresponding to the perception point cluster, the cluster position and the number of first perception points included in the perception point cluster.

[0255] In an exemplary embodiment, the determination module 2001 is specifically configured to use the plurality of second perception points as a target perception point set; or to select a target number of perception points from the plurality of second perception points and the plurality of first perception points as a target perception point set.

[0256] In an exemplary embodiment, the determination module 2001 is specifically used to obtain the maximum number of sensing points and the target number of sensing points included in the satellite coverage, the target number including the number of sensing points generated in previous interference power detection processes and the number of sensing points in the initial state; based on the maximum number of sensing points and the target number, determine whether the loop stop condition is met; if the loop stop condition is not met, execute the step of generating a target sensing point set based on multiple second sensing points.

[0257] In an exemplary embodiment, the determination module 2001 is specifically used to determine the difference between the maximum number of perception points and the target number; if the difference is greater than or equal to the total number of newly added perception points in a single round, it is determined that the loop stop condition is not met; if the difference is less than the total number of newly added perception points in a single round, it is determined that the loop stop condition is met.

[0258] In an exemplary embodiment, the determination module 2001 is also used to determine the number of additional sensing points based on the difference between the maximum number of sensing points and the target number if the loop stop condition is met; generate multiple third sensing points within the satellite coverage range based on the number of additional sensing points; and detect the interference power at each third sensing point to obtain the interference power value corresponding to each third sensing point.

[0259] In an exemplary embodiment, the determination module 2001 is specifically used to determine the electromagnetic situation information according to the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

[0260] In an exemplary embodiment, the determination module 2001 is specifically used to obtain the position information of each sensing point within the satellite coverage area; interpolate and fit the position information and interference power value of each sensing point based on a preset interpolation method to determine the electromagnetic situation information.

[0261] Each module in the electromagnetic situation information determination device of the above device can be implemented in whole or in part by software, hardware and a combination thereof. Each module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to each module above.

[0262] In an exemplary embodiment, a computer device is provided. The computer device may be a service, and its internal structure diagram may be as follows: Fig.21As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a method for determining electromagnetic situation information is implemented.

[0263] Those skilled in the art will understand that Fig.21 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0264] In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps described in any of the above method embodiments when executing the computer program.

[0265] 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 steps described in any of the above method embodiments are implemented.

[0266] In one embodiment, a computer program product is provided, including a computer program, which implements the steps described in any of the above method embodiments when executed by a processor.

[0267] A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment method can be completed 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 the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can 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 can include random access memory (RAM) or external cache memory, etc. As an 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 each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, an artificial intelligence (AI) processor, etc., but are not limited to this.

[0268] The technical features of the above embodiments may 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 application.

[0269] The above-described embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the present application. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the attached claims.

Claims

1. A method for determining electromagnetic situation information, characterized in that: The method comprises: Perform multiple interference power detections within the satellite coverage of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage according to the results of the interference power detection, wherein the results of the interference power detection are used to characterize the degree of uplink interference of the other satellite system to the own satellite; The i-th interference power detection in the multiple interference power detections includes: Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each of the first sensing points; According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set for the i+1th interference power detection.

2. The method according to claim 1, characterized in that: Generating a plurality of second sensing points within the satellite coverage range according to the interference power values ​​corresponding to the first sensing points includes: Clustering the plurality of first perception points according to the interference power value corresponding to each of the first perception points to obtain a plurality of perception point clusters; Determining the number of newly added sensing points corresponding to each of the sensing point clusters according to the clustering power values ​​corresponding to each of the sensing point clusters; According to the number of newly added perception points corresponding to each of the perception point clusters, the second perception point is generated in the perception area corresponding to each of the perception point clusters, wherein the perception area corresponding to the perception point cluster is the area where the perception points included in the perception point cluster are located within the satellite coverage area.

3. The method according to claim 2, characterized in that The clustering of the plurality of first perception points according to the interference power value corresponding to each of the first perception points to obtain a plurality of perception point clusters includes: Clustering the plurality of first perception points according to the interference power value corresponding to each of the first perception points to obtain an initial perception point cluster; For each of the initial perception point clusters, clustering processing is performed on the initial perception point cluster according to the position information of each of the first perception points in the initial perception point cluster to obtain a plurality of subclusters included in the initial perception point cluster; The subclusters included in each of the initial perception point clusters are used as the perception point clusters.

4. The method according to claim 2, characterized in that: The number of newly added sensing points corresponding to the sensing point cluster is positively correlated with the clustering power value corresponding to the sensing point cluster.

5. The method according to claim 3, characterized in that: The determining, according to the clustering power value corresponding to each of the sensing point clusters, the number of newly added sensing points corresponding to each of the sensing point clusters includes: Obtaining the total number of newly added sensing points in a single round, and determining the allocation ratio corresponding to each of the sensing point clusters according to the clustering power values ​​corresponding to each of the sensing point clusters; According to the allocation ratio corresponding to each of the perception point clusters, a number of newly added perception points is allocated to each of the perception point clusters from the total number, wherein the sum of the numbers of newly added perception points allocated to each of the perception point clusters is equal to the total number.

6. The method according to claim 5, characterized in that The determining, according to the clustering power values ​​corresponding to the sensing point clusters, the allocation ratios corresponding to the sensing point clusters includes: Obtain the cluster position corresponding to each sensing point cluster; For each of the sensing point clusters, an allocation ratio is determined according to a cluster power value corresponding to the sensing point cluster, a cluster position, and the number of first sensing points included in the sensing point cluster.

7. The method according to claim 1, characterized in that The generating a target perception point set according to the plurality of second perception points comprises one of the following: Taking the plurality of second perception points as the target perception point set; A target number of perception points are selected from the plurality of second perception points and the plurality of first perception points to serve as the target perception point set.

8. The method according to claim 1, characterized in that The generating a target perception point set according to the plurality of second perception points comprises: Obtaining a maximum number of sensing points, and obtaining a target number of sensing points included in the satellite coverage, wherein the target number includes the number of sensing points generated in previous interference power detection processes and the number of sensing points in an initial state; Determining whether a loop stop condition is met according to the maximum number of sensing points and the target number; If the loop stop condition is not met, the step of generating the target perception point set according to the multiple second perception points is executed.

9. The method according to claim 8, characterized in that The determining whether a loop stop condition is satisfied according to the maximum number of sensing points and the target number includes: Determining a difference between the maximum number of perception points and the target number; If the difference is greater than or equal to the total number of newly added sensing points in a single round, it is determined that the loop stop condition is not met; If the difference is less than the total number of newly added perception points in a single round, it is determined that the loop stop condition is met.

10. The method according to claim 8, characterized in that The method further comprises: If the loop stop condition is met, the number of additional perception points is determined according to the difference between the maximum number of perception points and the target number; Generating a plurality of third sensing points within the satellite coverage area according to the number of the additional sensing points; The interference power at each of the third perception points is detected to obtain the interference power value corresponding to each of the third perception points.

11. The method according to any one of claims 1 to 10, characterized in that: The step of determining the electromagnetic situation information corresponding to the satellite coverage range according to the result of the interference power detection includes: The electromagnetic situation information is determined according to the interference power value corresponding to each sensing point within the satellite coverage area determined by each interference power detection.

12. The method according to claim 11, characterized in that Determining the electromagnetic situation information according to the interference power value corresponding to each sensing point within the satellite coverage range determined by each interference power detection includes: Acquire the location information of each of the sensing points within the coverage area of ​​the satellite; Based on a preset interpolation method, interpolation fitting is performed on the position information and interference power value of each sensing point to determine the electromagnetic situation information.

13. An electromagnetic situation information determination device, characterized in that: The device comprises: A determination module, configured to perform multiple interference power detections within the satellite coverage of the own satellite, and determine the electromagnetic situation information corresponding to the satellite coverage according to the result of the interference power detection, wherein the result of the interference power detection is used to characterize the degree of uplink interference of the other satellite system to the own satellite; The i-th interference power detection in the multiple interference power detections includes: Acquire an initial sensing point set, where the initial sensing point set includes a plurality of first sensing points within the satellite coverage area, and acquire an interference power value at each of the first sensing points; According to the interference power values ​​corresponding to the first perception points, multiple second perception points are generated within the satellite coverage range, and a target perception point set is generated based on the multiple second perception points. The distribution of the multiple second perception points is determined according to the interference power values ​​corresponding to the first perception points. The target perception point set is the initial perception point set for the i+1th interference power detection.

14. 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 12 are implemented.

15. 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 12 are implemented.

16. 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 12 are implemented.

Images