A Multipoint Beam Area Planning Method and System for Satellite Communication

Abstract

The present invention provides a method for single-point beam area planning for satellite communication, including: calculating the distance between the grid center of the sub-satellite point and the grid center of the specified single-point beam center, the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered, and the distance between the grid center of the single-point beam center and the random grid center; calculating the distance between the satellite and the grid center of the single-point beam center and the distance between the satellite and the random grid center based on the calculated distances; calculating the included angle between the direction where the single-point beam center is located and the direction from the satellite to the random grid center based on the distances; determining whether the included angle satisfies a preset included angle condition, and forming the coverage area of the single-point beam by all the grid centers that satisfy the preset included angle condition.

Description

Technical Field

[0001] The present invention relates to the field of wireless communication, specifically to the field of satellite communication, and more specifically to a method for multi-point beam area planning for a high-throughput satellite (HTS) communication system. Background Art

[0002] Satellite communication has the inherent advantages of wide-area coverage and being unrestricted by distance and geographical conditions. It is an important part of 6G communication in the future space-air-ground integration. As non-terrestrial communication complements terrestrial cellular network communication, good satellite planning is conducive to promoting the realization of the future 6G communication network.

[0003] The purpose of satellite planning is mainly to plan the satellite's beam to achieve full coverage of the specified communication area. Among them, in order to achieve full coverage of the specified communication area to be covered, it is necessary to perform multi-point beam area planning on the specified communication area to be covered.

[0004] Most of the existing research technologies assume that the coverage area corresponding to the point beam is an ideal area with a regular shape. However, in an actual HTS communication system, the coverage area corresponding to the point beam is often an irregular shape. Therefore, when performing multi-point beam area planning on the specified communication area to be covered based on the existing research technologies, problems such as serious overlapping coverage between point beams, incomplete coverage of the communication area to be covered, waste of point beam resources, serious interference between point beams, and reduction of the performance of the HTS communication system will occur. Therefore, when facing an actual HTS communication system, designing a reasonable multi-point beam area planning method is the key to improving the performance of the HTS communication system.

[0005] Based on the ideal situation where the coverage area corresponding to a spot beam is a regular shape, domestic and foreign researchers have conducted a series of studies on the solution of the coverage range of spot beams and the regional planning of spot beams. For example, in reference [1], it is assumed that the coverage area shape of a single spot beam is circular or elliptical. After determining the center pointing of the spot beam, the process of projecting the spot beam onto the ground area coverage is regarded as a cone. Using the relevant knowledge of spatial points, lines, and planes, the boundary of the coverage area of a single spot beam is obtained by calculating the intersection of the cone projected by the spot beam and the Earth, and thus the coverage area of the spot beam is determined. In reference [2], it is assumed that the coverage surface of the spot beam is elliptical, and the elliptical parameters including the major semi-axis, minor semi-axis, and center point are calculated using plane geometry knowledge to solve the coverage area of the spot beam. In reference [3], a single-wave pointing coverage model is established, the beam cone is divided into several ray clusters, and the calculation of the coverage area of the spot beam is realized by solving the beam pointing point. In reference [4], under the ideal situation where the coverage area corresponding to the spot beam is a regular shape, a hybrid spot beam regional planning method applicable to HTS is proposed. This method can achieve the optimization goal of seamless coverage with the least number of spot beams through the particle swarm intelligent optimization algorithm. In reference [5], a multi-beam regional planning strategy is proposed. This strategy adopts a beam division strategy from the outside to the inside, scans and beam covers layer by layer from the outermost layer of the communication area to the inside until optimal coverage is achieved, and the beam interference problem is alleviated through frequency reuse. In reference [6], the Cellular Beam Footprint Planning (CBFP) method is proposed, but this method has a serious problem of overlapping coverage between spot beams, which affects the system performance.

[0006] However, due to the arc characteristics of the Earth's surface, the coverage area corresponding to the spot beam will deform to varying degrees with the change of longitude and latitude and the different attitude offsets of HTS. Therefore, when the existing research technologies are directly applied to the HTS communication system, problems such as serious overlapping coverage between spot beams, unreasonable waste of spot beam resources, and reduction of the performance of the HTS communication system will occur.

[0007] List of references:

[0008] Reference [1] Li Dezhi, Dou Zhaohui, Lü Bo, Wei Bo. Research on the Algorithm of Satellite Spot Beam Coverage Area [J]. Manned Spaceflight, 2009, 15(04): 53 - 55. (Chinese)

[0009] Reference [2] Dong Yanlei, Wang Chunting, Sun Wei. Research and Application of the GEO Satellite Spot Beam Coverage Algorithm [J]. Radio Communications Technology, 2016, 42(04): 83 - 86.

[0010] Reference [3] Wei Cheng, Wei Yongshang. Calculation method for multi-beam ground communication characteristics of satellite antennas based on GPU [P]. Heilongjiang Province: CN112559192A, 2021-03-26.

[0011] Reference [4] Zhao Yang, Wang Hailong, Huang Kaibo, Bai Jing. Design of high-throughput satellite coverage scheme based on population density distribution [J]. Digital Communication World, 2018(10): 28-29.

[0012] Reference [5] Cheng Yujian, Zhong Kaiyuan, Hou Huiying, Fan Yong, Yang Haining, He Zongrui, Li Tingjun, Zhao Minghua, Wang Hongbin, Wu Yafei. Method for dividing multi-beam coverage area of low-earth orbit satellites: Sichuan Province, CN112672423A [P]. 2021-04-16.

[0013] Reference [6] H. Cao, Y. Su, Y. Zhou, and J. Hu, “QoS guaranteed load balancing in broadband multi-beam satellite networks,” IEEE International Conference on Communications (ICC), pp. 1–6, 2019. Summary of the Invention

[0014] Therefore, the object of the present invention is to overcome the defects of the above-mentioned prior art and provide a method capable of accurately and reasonably planning multi-point beam regions.

[0015] The object of the present invention is achieved by the following technical solutions:

[0016] According to a first aspect of the present invention, there is provided a method for single - point beam area planning for satellite communication, wherein the earth's surface is equally divided into multiple parts according to a preset grid division accuracy. The grid center represents the corresponding grid area and the area to be covered, and the area to be covered is mapped into a non - empty matrix. The mapping of the grid center in the area to be covered on the non - empty matrix corresponds to an element in the non - empty matrix. The method includes the steps of: S1. Calculate the distance between the grid center of the sub - satellite point and the grid center of the specified single - point beam center, the distance between the grid center of the sub - satellite point and a random grid center in the area to be covered, and the distance between the grid center of the single - point beam center and the random grid center; S2. Calculate the distance between the satellite and the grid center of the single - point beam center and the distance between the satellite and the random grid center based on the distances calculated in step S1; S3. Calculate the included angle between the direction of the single - point beam center and the direction from the satellite to the random grid center based on the distances calculated in step S2; S4. Determine whether the included angle satisfies a preset included - angle condition, and form the coverage area of the single - point beam with all the grid centers that satisfy the preset included - angle condition.

[0017] In some embodiments of the present invention, the value of the grid division accuracy is 1 / 10. The grid division accuracy of 1 / 10 means that the unit longitude and unit latitude are respectively equally spaced and divided into 10 equal parts, so that the area corresponding to the unit longitude and unit latitude is divided into 100 grids in total.

[0018] In some embodiments of the present invention, in step S1, the distance from the grid center of the sub - satellite point to the specified single - point beam center is calculated in the following manner:

[0019]

[0020] where Z represents the grid center of the sub - satellite point, P k represents the specified single - point beam center, k represents the specified single - point beam, represents the distance from the grid center of the sub - satellite point to the single - point beam center, R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub - satellite point, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam k;

[0021] The distance between the grid center of the sub - satellite point and a random grid center in the area to be covered is calculated in the following manner:

[0022]

[0023] Among them, Z represents the grid center of the sub-satellite point, and A represents a random grid center. represents the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered, and R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the said random grid center;

[0024] Calculate the distance between the grid center of the single-point beam center and the random grid center in the following manner:

[0025]

[0026] Among them, A represents a random grid center, and P k represents the single-point beam center, represents the distance between the grid center of the single-point beam center and the random grid center, k represents the specified single-point beam, and R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the said random grid center, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam k.

[0027] In some embodiments of the present invention, in the step S2, calculate the distance between the satellite and the grid center of the single-point beam center in the following manner:

[0028]

[0029] Among them, represents the distance between the satellite and the grid center of the single-point beam center, and R e represents the radius of the earth, and H s represents the distance between the satellite and its sub-satellite point, represents the distance between the grid center of the sub-satellite point and the grid center of the specified single-point beam center;

[0030] Calculate the distance between the satellite and the said random grid center in the following manner:

[0031]

[0032] Among them, L SADenote the distance between the satellite and the center of the random grid as R e Denote the radius of the Earth as H s Denote the distance between the satellite and its sub-satellite point Denote the distance between the grid center of the sub-satellite point and the center of a random grid within the area to be covered

[0033] In some embodiments of the present invention, in step S3, calculate the included angle between the direction of the single-point beam center and the direction from the satellite to the center of the random grid in the following manner:

[0034]

[0035] Wherein Denote the distance between the satellite and the grid center of this single-point beam center as L SA Denote the distance between the satellite and the center of the random grid Denote the distance between the grid center of this single-point beam center and the center of the random grid

[0036] In some embodiments of the present invention, the preset included angle condition is: Wherein Denote the included angle between the direction of the single-point beam center and the direction from the satellite to the center of the random grid as θ 3dB Denote the included angle between the direction where the radiation gain of the satellite transmitting antenna is less than the maximum gain value by 3 dB and the direction of the point beam center

[0037] According to the second aspect of the present invention, there is provided a multi-point beam area planning method for satellite communication. The method includes the following steps: T1. Obtain the model parameters of the HTS communication system to be planned, the relevant parameters of the Earth, the longitude and latitude of the specified communication area to be covered, and the configuration parameters of the satellite; T2. Based on the parameters obtained in step T1 and the preset maximum overlap coverage rate, perform multiple rounds of iterative calculations for a preset number of rounds to obtain the minimum number of point beams that meet the maximum overlap coverage rate constraint. Wherein step T2 includes using the simulated annealing algorithm to perform the following steps in each round according to the preset number of rounds: T21. Calculate the coverage range corresponding to the set of point beam centers composed of the current number of point beams according to the method described in the first aspect of the present invention, and determine whether it can fully cover the area to be covered; T22. Calculate the overlap coverage rate corresponding to the set of point beam centers that can fully cover the area to be covered, determine whether the overlap coverage rate corresponding to this set of point beam centers is greater than the overlap coverage rate corresponding to the historical set of point beam centers, and use the set of point beam centers with the larger overlap coverage rate among the two as the minimum number of point beams corresponding to the current round that meet the maximum overlap coverage rate constraint

[0038] In some embodiments of the present invention, the step T1 further includes setting a maximum overlap coverage rate, an initial number of point beams, and a maximum number of point beams.

[0039] In some embodiments of the present invention, the maximum overlap coverage rate is 30%, and the initial number of point beams and the maximum number of point beams are randomly set.

[0040] In some embodiments of the present invention, the constraint of the maximum overlap coverage rate is as follows:

[0041] Constraint conditions for the area to be covered:

[0042] Constraint conditions for specifying the coverage area of a single point beam:

[0043] Constraint conditions for the overlap coverage rate between any two point beams:

[0044] Overlap coverage constraint condition for the area to be covered: η Ω ≤η max ;

[0045] Constraint condition for the payload of HTS point beam resources: K≤K max ;

[0046] where K represents the number of point beams, k represents a point beam, Ψ Ω represents a non-empty matrix corresponding to the area to be covered by the specified communication, Ψ k represents a non-empty matrix corresponding to the coverage area of the point beam k, N + represents the set of positive integers, h represents a point beam, g represents a point beam, η h,g represents the overlap coverage rate between the point beam h and the point beam g, η max represents the maximum overlap coverage rate, η Ω represents the overlap coverage rate of the area to be covered by the specified communication.

[0047] According to the third aspect of the present invention, a satellite system is provided, and the system includes: a high-throughput satellite, satellite terminal users, and a network control center. The system performs multi-point beam area planning by using the steps of the methods described in the first aspect and the second aspect of the present invention.

[0048] Compared with the prior art, the advantages of the present invention are as follows:

[0049] 1. In view of the fact that the coverage area corresponding to a spot beam will undergo varying degrees of deformation with changes in longitude and latitude and different HTS attitude offsets, the present invention designs a calculation method for the coverage area of a single spot beam based on a specified gridded communication area to solve the coverage range of any shape of the spot beam within the specified communication area, improving the accuracy of solving the coverage range of the spot beam.

[0050] 2. In view of the inevitable overlapping coverage between spot beams, the present invention defines the maximum overlapping coverage rate. Under the constraint of the maximum overlapping coverage rate, a reasonable multi-spot beam area planning for the area to be covered is realized to obtain the minimum number of spot beams capable of fully covering the area to be covered, effectively reducing the number of spot beams capable of achieving global coverage.

[0051] In addition, compared with the ideal situation where the coverage area corresponding to the spot beam in the existing method is a regular shape, the present invention has stronger practicability and accuracy, and can improve the throughput of the spot beam of the HTS communication system while reducing the interference between spot beams. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] The following further describes the embodiments of the present invention with reference to the drawings, where:

[0053] Figure 1 is a satellite system structure diagram according to an embodiment of the present invention;

[0054] Figure 2 is a structure diagram of the coverage area of a single spot beam according to an embodiment of the present invention;

[0055] Figure 3 is a multi-spot beam area planning diagram using the multi-spot beam area planning method applied to the present invention according to an embodiment of the present invention;

[0056] Figure 4 is a multi-spot beam area planning diagram using the CBFP method under ideal conditions according to an embodiment of the present invention;

[0057] Figure 5 is a multi-spot beam area planning diagram using the CBFP method under actual conditions according to an embodiment of the present invention;

[0058] Figure 6 is a comparison diagram of the SINR performance cumulative distribution density function on a satellite terminal user according to an embodiment of the present invention;

[0059] Figure 7 is a comparison diagram of the CDF performance index of the spot beam throughput according to an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0060] In order to make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below through specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0061] As mentioned in the background art, traditional methods in the prior art will cause problems such as serious overlapping coverage between point beams, incomplete coverage of the communication area to be covered, waste of point beam resources, serious interference between point beams, and degradation of the performance of the HTS communication system when performing multi-point beam area planning for the communication area to be covered. In view of the defects of the existing methods, the present invention proposes a method capable of accurately and reasonably performing multi-point beam area planning, which grids the communication area to be covered, calculates the range of the coverage area corresponding to a single point beam, and calculates the minimum number of point beams that can completely cover the communication area to be covered and meet the maximum overlapping coverage rate constraint based on the coverage area of the single point beam.

[0062] Before introducing the detailed steps of each part of the present invention, the main technical idea of the present invention will be introduced first:

[0063] The method for accurately and reasonably performing multi-point beam area planning proposed by the present invention for the HTS communication system mainly includes two aspects: the calculation method for the non-regular coverage area of a single point beam and the planning method for the non-regular coverage area of multi-point beams. First, the present invention divides the specified communication area into grids, and then proposes a calculation method for the coverage area of a single point beam based on the gridded communication area (Fine Grid Division based Single-Spot Beam Coverage Calculation, FGD-SBCC) according to the definition of the single point beam width. The calculation method for the coverage area of the single point beam can be used to solve the coverage area of a single point beam of any shape. Finally, the present invention defines the overlapping coverage rate based on the overlapping coverage generated between point beams and the finiteness of HTS point beam resources, and proposes a multi-point beam area planning method based on the coverage area of a single point beam (Fine Grid Division based Multi-Spot Beam Footprint Planning, FGD-MBFP) to solve the multi-point beam area planning method with the optimization goal of achieving seamless coverage of the entire specified communication area with the minimum number of point beams under the maximum overlapping coverage rate constraint.

[0064] To better understand the present invention, the present invention will be described in detail from four aspects: the satellite system model, the gridding of the specified communication area, the calculation method for the coverage area of a single point beam based on the gridded communication area, and the multi-point beam area planning method based on the coverage area of a single point beam through specific embodiments below.

[0065] 1. Satellite system model

[0066] The multi - point beam area planning method proposed by the present invention can be applied to any satellite communication system containing multiple HTSs. However, for a clearer description of the relevant technical features of the present invention, the present invention is introduced by taking a satellite system model containing one HTS as an example, as Figure 1 shown. The satellite system model contains one HTS, a Network Control Center (NCC), and multiple satellite terminal users (UTs). By adjusting the number, direction, and width of the beams, the satellite system can achieve full coverage of the entire communication area. In the satellite system, the HTS provides K point beams through the NCC. The direction of the point beam is determined by the position of the point beam center, and the point beam center refers to the position where the radiation gain of the satellite transmitting antenna is the largest. The point beam width is represented by the half - power beam angle θ 3dB of the point beam, where θ 3dB is the angle between the direction where the radiation gain of the satellite transmitting antenna is less than the maximum gain value by 3 dB and the direction of the maximum gain (point beam center). For a better understanding and demonstration of the present invention, the present invention defines C Ω to represent the specified communication area to be covered, C k to represent the coverage area of point beam k, P k to represent the beam center of point beam k, and Π = {P k |k ∈ [1, K]} to represent the set of all point beam centers, where K refers to the number of point beams provided by the HTS. It should be noted that C k is both the coverage area of point beam k and the area corresponding to the radiation gain of the transmitting antenna of point beam k being reduced by 3 dB compared to the maximum gain, that is, the area corresponding to θ 3dB . For the interference problem between point beams, the satellite system alleviates the interference between point beams through frequency reuse, where the frequency reuse factor is denoted as R.

[0067] 2. Grid - based specified communication area

[0068] The present invention grids the specified communication area and represents the coverage range corresponding to each grid with the grid center. The specified communication area also refers to the area to be covered, and all the grid centers within the area to be covered represent the area to be covered. According to an embodiment of the present invention, the present invention equally divides the unit longitude and unit latitude into d parts (d ∈ N + , N +(where \(d\) represents a positive integer), and define the grid division precision \(\delta\), where \(\delta = 1 / d\), \(\delta\in(0,1]\). The grid division precision value adopted in the present invention is \(1 / 10\). It should be noted that the grid division precision value adopted in the present invention is not unique, and the grid division precision value may be different according to different application scenarios. Based on the grid division precision, the present invention divides the longitude range and latitude range of the area to be covered into \(M\) parts and \(N\) parts respectively and generates \(M\times N\) grids, and when the grid division precision \(\delta\) is small enough, the coverage range corresponding to each grid is represented by the grid center. Among them, the area to be covered is represented by \(C\). Ω is represented, and the longitude range and latitude range are \([Lon min , Lon max and \([Lat min , Lat max respectively. The calculation formulas for the \(M\) parts and \(N\) parts are and The \(Lon min and \(Lon max respectively refer to the minimum longitude and maximum longitude of the area to be covered \(C Ω and \(Lon min , \(Lon max \in[0^{\circ}, 180^{\circ}E]\). The \(Lat min and \(Lat max respectively refer to the minimum latitude and maximum latitude of the area to be covered \(C Ω and \(Lat min , \(Lat max \in[0^{\circ}, 90^{\circ}N]\).

[0069] After meshing the area to be covered, the present invention obtains all the grid centers in the meshed area to be covered and maps the area to be covered into a non-empty matrix, and represents the area to be covered by the non-empty matrix. Among them, the mapping of the grid center in the area to be covered on the non-empty matrix corresponds to an element in the non-empty matrix. According to an embodiment of the present invention, the present invention represents the area to be covered \(C by the non-empty matrix Ω . The matrix element Ω in the non-empty matrix \(\varPsi is the mapping of the grid center \(g Ω in the area to be covered \(C ij on the non-empty matrix \(\varPsi Ω . When the grid center \(g ij is in the area to be covered \(C Ω , the matrix element Otherwise, the matrix element Therefore, when all the grid centers are in the area to be covered \(C Ω , for we can all get Among them, the matrix element with subscript ij represents the matrix element in the non-empty matrix Ψ Ω at the row and column numbers of its position. The longitude and latitude of the grid center g ij are The subscript ij of the grid center g ij represents the index of the longitude and latitude where the grid center g ij is located. Since the longitude range [Lon min , Lon max and the latitude range [Lat min , Lat max can be equally divided into and parts respectively, the range of the subscript ij of the grid center g ij is i ∈ [1, M] and j ∈ [1, N] respectively, and the range of the non-empty matrix Ψ Ω is M × N.

[0070] Similar to the method for obtaining the range of the area to be covered, the present invention obtains all the grid centers within the coverage area of a single-point beam and maps the coverage area of the single-point beam into a non-empty matrix, and represents the coverage area of the single-point beam with the non-empty matrix. Among them, the mapping of the grid center within the coverage area of the single-point beam on the non-empty matrix corresponds to an element in the non-empty matrix. According to an embodiment of the present invention, the present invention uses the non-empty matrix to represent the coverage area of the specified point beam k. It is easy to know that the coverage area C k of the point beam k is included in the area to be covered C Ω . Then the non-empty matrix Ψ k corresponding to the coverage area C k of the point beam k is also included in the non-empty matrix Ψ Ω corresponding to the area to be covered C Ω , that is And the matrix element k of the non-empty matrix Ψ is the mapping of the grid center g k within the coverage area C xy of the point beam k on the non-empty matrix Ψ k . When the grid center g xy is within the coverage area C k of the point beam k, the matrix element Otherwise, the matrix element Therefore, when the grid center is within the coverage area C k of the point beam k, for we can all get Among them, the matrix element with subscript xy represents the matrix element The row and column numbers of the position in the non-empty matrix Ψ k The matrix elements The superscript k of represents the coverage area C of the point beam k k The grid center g xy The longitude and latitude And the subscript xy of the grid center g xy Represents the index of the longitude and latitude where the grid center g xy Is located. It should be noted that when the grid division accuracy δ is small enough, the grid center within the coverage area C of the point beam k k Can represent the beam center P of the point beam k The longitude and latitude coordinates of the grid center k Is Its subscript Represents the index of the longitude and latitude where the grid center Is located, and Since the beam center P of the point beam k Is located within the coverage area of the point beam k, the beam center of the point beam k is mapped to the matrix element k In the non-empty matrix Ψ k The value of is 1, that is

[0071] In the process of multi-point beam planning for the to-be-covered area C Ω Since the overlapping coverage between point beams is inevitable, the present invention needs to control the overlapping coverage rate within a specified range while achieving full coverage of the to-be-covered area C Ω Therefore, in order to better understand and demonstrate the present invention, the present invention defines the overlapping coverage rate of the to-be-covered area C Ω As η Ω The overlapping coverage rate η Ω Refers to the ratio of the number of overlapping grid centers to the total number of grid centers within the to-be-covered area C Ω And defines the non-empty matrix corresponding to all point beam coverage areas as U, ab Wherein, u k Represents the matrix element where the grid center within the point beam coverage area is mapped to the non-empty matrix U, K represents the number of point beams, Ψ k Represents the non-empty matrix corresponding to the coverage area C of the point beam k Ω M×N represents the range of the non-empty matrix U, and the value of M×N is the total number of grid centers within the to-be-covered area C Ω When the coverage areas corresponding to all point beams fully cover the to-be-covered area C Ω And there are no overlapping grid centers within the coverage areas corresponding to each point beam, the matrix element u in the non-empty matrix U at this timeab = 1 when the coverage areas corresponding to all point beams completely cover the area to be covered C Ω and there are overlapping grids within the coverage areas corresponding to each point beam. At this time, for the matrix elements u in the non-empty matrix U ab ≥ 2, while in the areas outside the coverage areas corresponding to all point beams, for the matrix elements u of the non-empty matrix U ab = 0. Therefore, the number of grid centers with overlapping coverage within the area to be covered C Ω is equal to the number of matrix elements u in the non-empty matrix U where ab u ≥ 2, where the number of grid centers with overlapping coverage within the area to be covered C Ω is represented by , and the total number of all grid centers within the area to be covered C Ω is represented by . Thus, the overlapping coverage rate η Ω of the area to be covered C Ω can be expressed by the following calculation formula: To better understand and demonstrate the present invention, when calculating the overlapping coverage rate of the coverage areas corresponding to two point beams, define as the number of grid centers with overlapping coverage of point beam h and point beam g within the area to be covered, that is Define as the total number of grid centers within the coverage areas of point beam h and point beam g, that is Therefore, the calculation formula for the overlapping coverage rate between two point beams such as point beam h and point beam g can be expressed as:

[0072] III. Calculation method for the coverage area of a single point beam based on the gridified communication area

[0073] When facing an actual HTS communication system, based on the gridified area to be covered, the present invention calculates the range of the coverage area of a single point beam through the calculation method for the coverage area of a single point beam, and the range of the coverage area of a single point beam obtained through calculation and solution is the basis for realizing the multi-point beam area planning for the area to be covered. To better understand and demonstrate the present invention, before elaborating on the calculation method for the coverage area of a single point beam, the present invention first introduces the relevant variables for calculating the range of the coverage area of a single point beam, such as Figure 2 as shown, define the orbital height of the HTS (marked as S in the figure) as H s , the radius of the earth as R e , the sub-satellite point of the HTS as Z, the earth's geocenter as O, the longitude and latitude coordinates of the HTS S as (Lon s , Lat s ), and point A as the random grid center point g xy, where \(x\in[1,M]\), \(y\in[1,N]\), and the longitude and latitude coordinates of point \(A\) are \((Lon a ,Lat a ). Given the beam center \(P k of the grid center where the beam center \(P k of the point beam \(k\) has longitude and latitude coordinates of is the angle \(\angle ASP k between the direction from HTS S to \(A\) and the direction from HTS S to the beam center \(P k of the point beam \(k\).

[0074] Next, the calculation method of the coverage area of a single point beam proposed by the present invention is introduced. According to an embodiment of the present invention, the steps for calculating the range of the coverage area of a single point beam of the present invention include steps T1 to T7:

[0075] Step T1: Use the calculation formula for the distance between any two points on the sphere to calculate the distance between the grid center of the sub-satellite point and the grid center of the specified single point beam center

[0076]

[0077] where \(Z\) represents the grid center of the sub-satellite point, \(P k represents the specified single point beam center, \(k\) represents the specified single point beam, \(R e represents the radius of the earth, \(Lon s and \(Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam \(k\);

[0078] Step T2: Based on the distance calculated in Step T1, use the cosine theorem of a triangle to calculate the distance between the satellite and the grid center of the single point beam center

[0079]

[0080] where \(R e represents the radius of the earth, \(H s represents the distance between the satellite \(S\) and its sub-satellite point, represents the distance between the grid center of the sub-satellite point and the grid center of the specified single point beam center;

[0081] Step T3: Use the calculation formula for the distance between any two points on the sphere to calculate the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered

[0082]

[0083] Among them, Z represents the grid center of the sub-satellite point, A represents a random grid center, and R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the said random grid center;

[0084] Step T4: Based on the distance calculated in step T3, use the cosine theorem of a triangle to calculate the distance L between the satellite and the said random grid center SA :

[0085]

[0086] Among them, R e represents the radius of the earth, H s represents the distance between the satellite S and its sub-satellite point, represents the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered;

[0087] Step T5: Use the formula for calculating the distance between any two points on a sphere to calculate the distance between the grid center of the single-point beam center and the said random grid center

[0088]

[0089] Among them, A represents a random grid center, P k represents the single-point beam center, R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the said random grid center, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam k;

[0090] Step T6: Based on the distances calculated in steps T2 to T5, use the cosine theorem of a triangle to calculate the included angle between the direction where the single-point beam center is located and the direction from the satellite to the said random grid center

[0091]

[0092] Among them denotes the distance between the satellite and the grid center of the single-point beam center, L SA denotes the distance between the satellite and the random grid center denotes the distance between the grid center of the single-point beam center and the random grid center

[0093] Step T7, determine whether the included angle described in step T6 satisfies a preset included angle condition, and the preset included angle condition is where θ 3dB denotes the included angle between the direction where the radiation gain of the satellite transmitting antenna is less than the maximum gain value by 3 dB and the direction of the point beam center, and all grid centers that satisfy the preset included angle condition form the coverage area of the single-point beam. According to the description in the above "II. Grid the specified communication area", "When then otherwise ", the coverage area C k corresponding to the non-empty matrix Ψ k of the point beam k can be updated according to the matrix element update formula, and the matrix element update formula is as follows:

[0094]

[0095] IV. Multi-point beam area planning method based on single-point beam coverage area

[0096] Before elaborating on the multi-point beam area planning method, the present invention first introduces multiple constraint conditions applied to the multi-point beam area planning method. Since in an actual HTS communication system, overlapping coverage between point beams is inevitable, the present invention constructs multiple constraint conditions based on the grid-shaped area to be covered to obtain the minimum number of point beams that can fully cover the area to be covered and satisfy the constraint conditions. The multiple constraint conditions respectively include the constraint conditions of the area to be covered: Constraint conditions for specifying the coverage area of a single-point beam: Constraint conditions for the overlapping coverage rate between any two point beams: Overlapping coverage constraint conditions for the area to be covered: and the constraint condition of the HTS point beam resource load: K ≤ K max . Where K represents the number of point beams, k represents a point beam, Ψ Ω represents the non-empty matrix corresponding to the area to be covered by the specified communication, Ψ k represents the non-empty matrix corresponding to the coverage area of the point beam k, N + represents the set of positive integers, h represents a point beam, g represents a point beam, η h,g represents the overlapping coverage rate between the point beam h and the point beam g, η maxDenote the maximum overlap coverage rate as η Ω Denote the overlap coverage rate of the specified communication area to be covered. Among the multiple constraints, the constraint of the area to be covered ensures that the multi-point beams can fully cover the area to be covered. The constraint of specifying the coverage area of a single-point beam ensures that the matrix corresponding to the coverage area of the single-point beam is non-empty. The overlap coverage rate constraint between any two point beams ensures that the overlap coverage rate between the two point beams is not greater than the maximum overlap coverage rate. The overlap coverage constraint of the area to be covered ensures that the overlap coverage rate of the area to be covered is not greater than the maximum overlap coverage rate. The constraint of the HTS point beam resource load ensures that the number of point beams in the area to be covered is not greater than the maximum number of point beams that the satellite can carry.

[0097] Next, a multi-point beam area planning method based on the coverage area of a single-point beam will be introduced. It should be noted that when the point beam width θ 3dB is given, the HTS multi-point beam area planning is determined by a set of point beam centers, that is, the multi-point beam center set Π={P k |k∈[1,K]}. The overall idea of the multi-point beam area planning method is based on the gridified area to be covered. Under the constraint of the HTS point beam resource load constraint, the number of point beams is gradually increased in steps of 1. During the process of increasing the number of point beams, the given number of point beams K is calculated according to the formula of the overlap coverage rate between two point beams and the overlap coverage constraint of the area to be covered to generate a set of point beam centers, and obtain a set of point beam centers that can fully cover the area to be covered with the least number of point beams and satisfy the maximum overlap coverage rate. However, since the optimization problem of achieving full coverage of the area to be covered with the least number of point beams proposed in the present invention is NP-hard and has a high complexity when solving this optimization problem, preferably, the present invention uses an improved simulated annealing algorithm (Modified Simulated Annealing, MSA) to solve this optimization problem. According to an embodiment of the present invention, the steps of the multi-point beam area planning for the area to be covered by the present invention include P1 to P4:

[0098] Step P1: Obtain the HTS communication system model parameters, earth-related parameters, the longitude and latitude of the specified communication area to be covered, and the satellite configuration parameters. Set the maximum overlap coverage rate, the initial number of point beams, and the maximum number of point beams in the simulated annealing algorithm, where the maximum overlap coverage rate is 30%, and the initial number of point beams and the maximum number of point beams are random values;

[0099] Step P2: Iterate the simulated annealing algorithm according to the preset number of rounds. During the operation of each round of the simulated annealing algorithm, gradually increase the number of point beams in steps of 1 based on the constraint conditions of the HTS spot beam resource load, and perform the following steps during the operation of each round of the simulated annealing algorithm:

[0100] Step P21: Generate a set of point beam centers based on the current number of point beams according to the formula for the overlapping coverage rate between two point beams and the overlapping coverage constraint conditions of the area to be covered;

[0101] Step P22: Calculate the coverage area corresponding to each single point beam center in the set of point beam centers described in Step P21 based on the calculation method of the coverage area of a single point beam, and the coverage areas corresponding to each single point beam center form the coverage area corresponding to the set of point beam centers;

[0102] Step P23: Determine whether the coverage range corresponding to the set of point beam centers composed of the current number of point beams can fully cover the area to be covered, and obtain the set of point beam centers that can fully cover the area to be covered;

[0103] Step P3: Obtain a set of point beam centers that can fully cover the area to be covered based on the constraint conditions of the area to be covered. Use the overlapping coverage rate η Ω of the area to be covered C Ω formula to calculate the overlapping coverage rate of this set of point beam centers, and determine whether the overlapping coverage rate corresponding to this set of point beam centers is greater than the overlapping coverage rate corresponding to the historical set of point beam centers. Take the set of point beam centers with the larger overlapping coverage rate among the two as the minimum number of point beams corresponding to the current round that satisfies the maximum overlapping coverage rate constraint. It should be noted that when the overlapping coverage rate corresponding to this set of point beam centers is greater than the maximum overlapping coverage rate, in order to avoid the simulated annealing algorithm falling into a locally optimal state, the present invention calculates a reference value to avoid falling into a locally optimal solution where represents the overlapping coverage rate of the current set of point beam centers, η Ω represents the overlapping coverage rate of the area to be covered C Ω and γ represents the annealing rate. When the reference value λ≥random(0,1), still take this set of point beam centers as the minimum number of point beams corresponding to the current round that satisfies the maximum overlapping coverage rate constraint;

[0104] Step P4: Calculate the overlapping coverage rate of the coverage areas corresponding to the set of beam center points of each number of point beams in each round of each simulation algorithm constructed based on the present invention until the preset number of rounds is reached, and retain the minimum number of point beams that satisfy the maximum overlapping coverage rate constraint, where the preset number of rounds refers to the maximum number of rounds set artificially that the simulated annealing algorithm can execute.

[0105] To illustrate the performance advantages of the present invention in the actual communication scenario of HTS, the present invention uses a simulation experiment to compare the multi-point beam area planning method applied to the present invention with the Cellular Beam Footprint Planning (CBFP) method. The experimental parameters of the simulation experiment are shown in Table 1. The working frequency band of the HTS communication system is the Ka band, the point beam width θ 3dB = 0.3°, the point beam bandwidth is 500 MHz, the frequency reuse factor R = 7, and the radiation pattern of the multi-point beam transmitting antenna follows the rec.ITU-S.672 standard. The formula for calculating the antenna gain in the rec.ITU-S.672 standard is:

[0106]

[0107] where the parameter The parameter a = 2.88 and the parameter b = 6.32.

[0108] Table 1 FGD-MBFP simulation parameter settings

[0109]

[0110]

[0111] Based on each simulation parameter in Table 1, as Figure 3 shown, obtain the multi-point beam area planning map for fully covering the specified communication area to be covered with longitude and latitude ranges of [120°E, 150°E] and [30°N, 60°N] under the constraint condition that the given maximum overlapping coverage rate η max = 30% for the multi-point beam area planning method applied to the present invention. And Figure 4 is the multi-point beam area planning map for covering the specified communication area to be covered with longitude and latitude ranges of [120°E, 150°E] and [30°N, 60°N] using the CBFP method in the ideal state, Figure 5 while

[0112] is the multi-point beam area planning map for fully covering the specified communication area to be covered with longitude and latitude ranges of [120°E, 150°E] and [30°N, 60°N] using the CBFP method in the actual state. Figure 5As shown, when the CBFP method is used to fully cover the specified communication area to be covered in the actual state, it can be clearly found that the shape of the single-point beam coverage area in the actual scenario is irregular, and the shape of the single-point beam coverage area will deform to varying degrees with the change of longitude and latitude. Figure 4 and Figure 5 By comparison, when the CBFP method is used to plan the multi-point beam area for the communication area to be covered in the actual state, excessive overlapping coverage will occur between the point beams, resulting in serious interference between the point beams. Figure 3 and Figure 5 By comparison, it can be known that the multi-point beam area planning method applied to the present invention can achieve full coverage of the communication area to be covered with a smaller number of point beams.

[0113] As shown in Table 2, when fully covering the same communication area to be covered, 132 point beams are required when using the CBFP method in the actual state, while only 50 point beams are required for the multi-point beam area planning method applied to the present invention. The reduced number of point beams exceeds 60%, significantly saving the HTS point beam resources. Moreover, the multi-point beam area planning method applied to the present invention reduces the overlapping coverage rate generated by the CBFP method from 98% to 25%, achieving full coverage of the communication area to be covered with the least number of point beams under the constraint of the maximum overlapping coverage rate, and alleviating the interference between the point beams to a certain extent, improving the performance of the HTS communication system.

[0114] Table 2 Comparison of simulation results between CBFP method and FGD-MBFP method

[0115]

[0116] Figure 6 Shows the comparison relationship between the cumulative distribution density function (CDF) of the Signal to Interference plus Noise Ratio (SINR) performance of the multi-point beam area planning method (FGD-MBFP) applied to the present invention and the CBFP method on satellite terminal users. It can be clearly found in Figure 6 that the SINR of all users in the HTS communication system using the multi-point beam area planning method applied to the present invention is above 0 dB, while about 10% of the users in the HTS communication system using the CBFP method have a SINR below 0 dB.

[0117] As shown in Table 3, the average SINR of users adopting the multi-point beam area planning method (FGD-MBFP) applied to the present invention is 8.15 dB, while the average SINR of users adopting the CBFP method is 5.17. Comparing the two average SINRs of users, it is found that the average SINR of users adopting the multi-point beam area planning method applied to the present invention has increased by 58.6%. Thus, it can be seen that the HTS communication system adopting the multi-point beam area planning method applied to the present invention can provide signals of better quality, which benefits from the fact that the multi-point beam area planning method applied to the present invention can effectively reduce the overlapping coverage rate between point beams, alleviate the interference between point beams, and improve the SINR performance of users in the HTS communication system.

[0118] Table 3 Comparison of User SINR between CBFP Method and FGD-MBFP Method

[0119] Method Minimum SINR [dB] Maximum SINR [dB] Average SINR [dB] FGD-MBFP 1.10 13.61 8.15 CBFP -12.27 12.84 5.17

[0120] Figure 7 Shows the comparison relationship between the multi-point beam area planning method (FGD-MBFP) applied to the present invention and the CBFP method in terms of the CDF performance index of point beam throughput. In Figure 7 it can be clearly found that the throughput of all point beams in the HTS communication system adopting the multi-point beam area planning method applied to the present invention can reach 0.8 Gbps and above, while only 60% of the point beams in the HTS communication system adopting the CBFP method can reach a throughput similar to that of the multi-point beam area planning method applied to the present invention.

[0121] As shown in Table 4, the average throughput of point beams in the HTS communication system adopting the multi-point beam area planning method (FGD-MBFP) applied to the present invention is 1.03 Gbps, while the average throughput of point beams in the HTS communication system adopting the CBFP method is 0.79 Gbps. Comparing the two average throughputs of point beams, it is found that the average throughput adopting the multi-point beam area planning method applied to the present invention has increased by 30.4%, which benefits from the fact that the multi-point beam area planning method applied to the present invention can effectively improve the SINR performance of end users.

[0122] Table 4 CDF Comparison of Point Beam Throughput between CBFP Method and FGD-MBFP Method

[0123]

[0124] It should be noted that although the above steps are described in a specific order, it does not mean that the steps must be executed in the above specific order. In fact, some of these steps can be executed concurrently or even in a different order, as long as the required functions can be achieved.

[0125] The present invention may be a system, a method, and / or a computer program product. The computer program product may include a computer-readable storage medium having thereon computer-readable program instructions for causing a processor to implement various aspects of the present invention.

[0126] A computer-readable storage medium may be a tangible device that retains and stores instructions for use by an instruction execution device. A computer-readable storage medium may include, for example, but is not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium would include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disc (DVD), a memory stick, a floppy disk, a mechanically encoded device such as a punch card or raised structures in a groove having instructions stored thereon, and any suitable combination of the foregoing.

[0127] The embodiments of the present invention have been described above. The above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The choice of terms used herein is intended to best explain the principles of the embodiments, the practical application, or improvements made to the technology in the market, or to enable other ordinary skilled persons in the art to understand the embodiments disclosed herein.

Claims

1. A single-point beam area planning method for satellite communication, wherein, The Earth's surface is equally divided into multiple parts according to a preset grid division accuracy. The grid center represents the corresponding grid area and the area to be covered, and the area to be covered is mapped to a non-empty matrix. The mapping of the grid center in the area to be covered on the non-empty matrix corresponds to an element in the non-empty matrix. The method includes the steps: S1. Calculate the distance between the grid center of the sub-satellite point and the grid center of the specified single-point beam center, the distance between the grid center of the sub-satellite point and a random grid center in the area to be covered, and the distance between the grid center of this single-point beam center and this random grid center; S2. Calculate the distance between the satellite and the grid center of this single-point beam center and the distance between the satellite and the random grid center based on the distances calculated in step S1; S3. Calculate the included angle between the direction where this single-point beam center is located and the direction from the satellite to the random grid center based on the distances calculated in step S2; S4. Determine whether the included angle satisfies a preset included angle condition, and form the coverage area of this single-point beam with all the grid centers that satisfy the preset included angle condition.

2. The method according to claim 1, characterized in that, The value of the grid division accuracy is 1 / 10, where the grid division accuracy 1 / 10 means that the unit longitude and unit latitude are respectively equally divided into 10 equal parts to divide the area corresponding to the unit longitude and unit latitude into 100 grids.

3. The method according to claim 2, characterized in that, In step S1, Calculate the distance from the grid center of the sub-satellite point to the specified single-point beam center in the following manner: Among them, Z represents the grid center of the sub-satellite point, and P k represents the specified single-point beam center, k represents the specified single-point beam, represents the distance from the grid center of the sub-satellite point to the center of this single-point beam, R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam k; Calculate the distance between the grid center of the sub-satellite point and a random grid center in the area to be covered in the following manner: Among them, Z represents the grid center of the sub-satellite point, and A represents a random grid center. represents the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered, and R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the random grid center. Calculate the distance between the grid center of this single-point beam center and this random grid center in the following manner: Among them, A represents a random grid center, P k represents the center of the single-point beam, represents the distance between the grid center of the single-point beam center and the random grid center, k represents the specified single-point beam, R e represents the radius of the earth, Lon s and Lat s respectively represent the longitude and latitude of the grid center corresponding to the sub-satellite point, Lon a and Lat a respectively represent the longitude and latitude of the random grid center, and respectively represent the longitude and latitude of the grid center corresponding to the beam center of the point beam k.

4. The method according to claim 3, characterized in that, In step S2, Calculate the distance between the satellite and the grid center of this single-point beam center in the following manner: Among them, represents the distance between the satellite and the grid center of the center of the single-point beam, R e represents the radius of the earth, H s represents the distance between the satellite and its sub-satellite point, represents the distance between the grid center of the sub-satellite point and the grid center of the specified single-point beam center; Calculate the distance between the satellite and the random grid center in the following manner: Among them, L SA represents the distance between the satellite and the center of the random grid, R e represents the radius of the Earth, H s represents the distance between the satellite and its sub-satellite point, represents the distance between the grid center of the sub-satellite point and a random grid center within the area to be covered.

5. The method according to claim 4, wherein In step S3, calculate the included angle between the direction where this single-point beam center is located and the direction from the satellite to the random grid center in the following manner: Among them, represents the distance between the satellite and the grid center of the single-point beam center, L SA represents the distance between the satellite and the random grid center, represents the distance between the grid center of the single-point beam center and the random grid center.

6. The method according to claim 1, characterized in that The preset included angle condition is: Wherein represents the included angle between the direction where the center of the single-point beam is located and the direction from the satellite to the center of the random grid, θ 3dB represents the included angle between the direction where the radiation gain of the satellite transmitting antenna is less than the maximum gain value by 3 dB and the direction where the center of the point beam is located.

7. A method for multi - point beam area planning for satellite communication, characterized in that, The method includes the following steps: T1. Obtain the model parameters of the HTS communication system to be planned, the relevant parameters of the Earth, the longitude and latitude of the specified communication area to be covered, and the configuration parameters of the satellite; T2. Based on the parameters obtained in step T1 and the preset maximum overlap coverage rate, perform multiple rounds of iterative calculations for a preset number of rounds to obtain the minimum number of point beams that satisfy the maximum overlap coverage rate constraint, where step T2 includes using the simulated annealing algorithm to perform the following steps in each round according to the preset number of rounds: T21. Calculate the set of point beam centers composed of the current number of point beams according to the method of any one of claims 1-6, and determine whether the coverage range corresponding to the set of point beam centers can achieve full coverage of the area to be covered; among them, the set of point beams is obtained in the following manner: Based on the to-be-covered area after grid division, under the constraints of the HTS spot beam resource payload, the number of spot beams is gradually increased in steps of 1. During the process of increasing the number of spot beams, the given number of spot beams is used to generate a set of spot beam centers according to the formula for the overlapping coverage rate between two spot beams and the overlapping coverage constraint conditions of the to-be-covered area, where K is the given number of spot beams, is the number of grid centers covered by spot beam h and spot beam g within the to-be-covered area, is the total number of grid centers within the coverage areas of spot beam h and spot beam g, and η h,g represents the overlapping coverage rate between spot beam h and spot beam g; T22. Calculate the overlapping coverage rate corresponding to the set of point beam centers that can achieve full coverage of the area to be covered, and determine whether the overlapping coverage rate corresponding to the set of point beam centers is greater than the overlapping coverage rate corresponding to the historical set of point beam centers. Use the set of point beam centers with the larger overlapping coverage rate among the two as the minimum number of point beams corresponding to the current round that satisfies the maximum overlapping coverage rate constraint. Among them, the maximum overlapping coverage rate constraint is: Constraints for the area to be covered: Constraints for specifying the coverage area of a single-point beam: Overlap coverage rate constraint condition between any two point beams: η h,g ≤η max For g ∈ [1, K], h ≠ g; Overlap coverage constraint condition for the area to be covered: η Ω ≤ η max ; Constraint conditions for the HTS spot beam resource payload: K ≤ K max ; Among them, K represents the number of point beams, k represents a point beam, and Ψ Ω represents a non-empty matrix corresponding to the area to be covered by the specified communication, and Ψ k represents a non-empty matrix corresponding to the area covered by the point beam k, and N + represents the set of positive integers, h represents a point beam, g represents a point beam, and η h,g represents the overlapping coverage rate between the point beam h and the point beam g, and η max represents the maximum overlapping coverage rate, and η Ω represents the overlapping coverage rate of the area to be covered by the specified communication.

8. The method according to claim 7, characterized in that The step T1 further includes setting the maximum overlapping coverage rate, the initial number of point beams, and the maximum number of point beams.

9. The method according to claim 8, wherein The maximum overlapping coverage rate is 30%, and the initial number of point beams and the maximum number of point beams are randomly set.

10. A satellite system, the system comprising: A high-throughput satellite, satellite terminal users, and a network control center, characterized in that the system uses the method described in any one of claims 7-9 for multi-point beam area planning.

11. A computer-readable storage medium, characterized in that, It stores a computer program, and the computer program can be executed by a processor to implement the steps of the method described in any one of claims 1-6, 7-9.

12. An electronic device, characterized in that, Including: One or more processors; A storage device for storing one or more programs, and when the one or more programs are executed by the one or more processors, the electronic device is caused to implement the steps of the method described in any one of claims 1-6, 7-9.

Images