CN112672423A - Low-orbit satellite multi-beam coverage area dividing method - Google Patents

Low-orbit satellite multi-beam coverage area dividing method Download PDF

Info

Publication number
CN112672423A
CN112672423A CN202011553560.1A CN202011553560A CN112672423A CN 112672423 A CN112672423 A CN 112672423A CN 202011553560 A CN202011553560 A CN 202011553560A CN 112672423 A CN112672423 A CN 112672423A
Authority
CN
China
Prior art keywords
circle
layer
coverage
beams
covering
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011553560.1A
Other languages
Chinese (zh)
Other versions
CN112672423B (en
Inventor
程钰间
钟开源
侯慧颖
樊勇
杨海宁
何宗锐
李廷军
赵明华
王洪斌
吴亚飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202011553560.1A priority Critical patent/CN112672423B/en
Publication of CN112672423A publication Critical patent/CN112672423A/en
Application granted granted Critical
Publication of CN112672423B publication Critical patent/CN112672423B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention provides a method for dividing a multi-beam coverage area of a low-orbit satellite, which adopts an outside-in area division strategy: covering the first edge scanning circle layer, determining the scanning angle and the number of covered beams of the current scanning circle layer to complete the beam covering of the current scanning circle layer, updating the edge parameters, and covering the scanning circle layer inwards layer by layer until the covering is completed. Under the condition of meeting the channel consistency, the optimal coverage of the ground equal-flux large-range multi-beam is realized, and the overall efficiency of the existing beam forming method is greatly improved. On the basis, a clustering color separation algorithm and a depth-first search-based frequency multiplexing color separation algorithm are provided, and the problem of low mutual interference color separation among non-uniform multi-beams is solved. The work has important significance for solving the division and the beam forming of the multi-beam coverage area of the low-orbit satellite.

Description

Low-orbit satellite multi-beam coverage area dividing method
Technical Field
The invention relates to the technical field of satellite communication, in particular to a multi-beam region division method for a low-orbit satellite.
Background
Satellite communication has the advantages of wide coverage area, high transmission quality, flexible deployment and the like, can provide seamless connection in a larger geographical area and even a global range, is convenient for rapid networking, and becomes an important component of civil and military communication systems. Compared with a geosynchronous orbit satellite, the low-orbit satellite has the characteristics of small communication delay and the like, more than 3 satellites can be seen to avoid signal shielding, the miniaturization of a user terminal is facilitated, and the low-orbit satellite becomes an important direction for the development of the current satellite.
The low earth orbit satellite adopts multi-beam subareas to cover the ground, namely, the whole coverage area is divided into a plurality of cells, and a plurality of independent spot beams are used for covering each cell. On the basis, the frequency resources of the whole system are divided into a plurality of frequency bands, and adjacent beams avoid mutual interference by adopting frequency division or code division. Compared with the satellite communication system covered by single beam, the satellite communication system with multi-beam and frequency reuse capability can obtain higher system capacity and flexibility. Of course, the small number of cells can effectively improve the utilization rate of frequency spectrum resources, reduce the interference between beams, and reduce the complexity and hardware cost of digital beam forming.
While terrestrial communications typically employ cellular sectorization, geostationary satellites also use a multi-spot beam sectorization approach similar to terrestrial cellular communications. Different from ground communication, the satellite-borne multi-beam antenna needs to meet the requirement of equal-flux coverage, namely the signal strength of the satellite reaching each beam cell on the ground is equal. The path loss for signal propagation varies due to the different distances between each cell and the satellite. The low-orbit satellite has large communication elevation angle and large path loss difference. Meanwhile, the antenna beam is also distorted during large-angle scanning, so that equal spot beams are difficult to realize. Under the conditions of meeting the equal flux coverage and the like, the shapes of the regions are efficiently and reasonably divided, and the technologies such as frequency reuse and the like are combined, so that the utilization rate of frequency spectrum resources is improved, and the complexity of beam forming is reduced.
Disclosure of Invention
The invention aims to solve the technical problem of providing a whole set of method for dividing a low-orbit satellite multi-beam coverage area, forming a multi-beam directional diagram and dividing frequency of multi-beams.
The invention adopts the technical scheme that the method for dividing the multi-beam coverage area of the low-orbit satellite comprises the following steps:
adopting an outside-in region division strategy: covering the first edge scanning circle layer, determining the scanning angle and the number of covered beams of the current scanning circle layer to complete the beam covering of the current scanning circle layer, updating the edge parameters, and covering the scanning circle layer inwards layer by layer until the covering is completed.
Specifically, a multi-beam covering method for completing full coverage when an antenna array is a full array is provided;
specifically, a multi-beam covering method for completing non-full covering when an antenna array is a full array is provided;
furthermore, in order to reduce the number of layers of the circle covered by the wave beam, a subarray type multi-beam shaping scheme is provided; furthermore, in order to solve the color separation problem caused by the multi-beam forming of the sub-array, a frequency multiplexing color separation method based on a clustering color separation algorithm and depth-first search is provided.
The method has the advantages that full-coverage and non-full-coverage multi-beam coverage models are established according to coverage scenes, an outside-in region division strategy and a sub-array based multi-beam forming method are provided on the basis, the optimal region distribution and directional diagram forming scheme meeting the coverage requirements can be quickly obtained on the premise of equal flux coverage and consistent channels, and the effect and the efficiency of the existing method are improved.
Drawings
FIG. 1 is a hierarchical overlay flow diagram;
FIG. 2 is a diagram of a layered full coverage model;
FIG. 3 is a diagram of a layered non-full coverage model;
FIG. 4 is a schematic diagram of a second circled subarray division;
FIG. 5 is a schematic diagram of a third circle of subarray division;
FIG. 6 shows the results of edge coverage in a 50 ° coverage area;
FIG. 7 is a diagram illustrating coverage of different parameters Δ;
figure 8 shows the result of subarray multi-beam coverage.
Detailed description of the invention
And establishing full-coverage and non-full-coverage multi-beam coverage models according to a coverage scene, and providing an outside-in region division strategy. The full coverage is that the equivalent isotropic radiated power EIRP at any position in the coverage area meets the coverage requirement. The non-full coverage is a more engineered application requirement, as is the case for achieving 80%, 90% EIRP coverage.
By adopting a multi-beam forming method based on a sub-array, the optimal area distribution and directional diagram forming scheme meeting the coverage requirement can be quickly obtained on the premise of equal flux coverage and consistent channels;
outside-in zone partitioning strategy: covering edge beams, completing one-layer beam division, updating edge parameters, and iterating layer by layer until the covering is completed, as shown in fig. 1.
Embodiments 1 and 2 are that the antenna array completes beam coverage for a full array; in order to reduce the number of turns of the beam coverage, embodiment 3 further proposes a manner of completing the beam coverage by partial array operation (a subarray-type multi-beam forming scheme), and then further proposes a color separation scheme in order to solve the color separation problem caused by the proposed subarray-type multi-beam forming.
Example 1
As shown in FIG. 2, for the full coverage model, the maximum in the graphThe radius of the circle represents the maximum angle required to be covered currently and corresponds to the edge parameter thetan. Determining the first layer beam coverage as + -theta1(n ═ 1), a gain flatness p is defined, i.e. the beam width of the coverage cell is pdB coverage. From the antenna end input, one half pdB beamwidth θ can be obtained when the beam is pointed normal0pAnd half pdB beamwidth θ swept to angle ααp
Under the condition of a full array, the wave beams are uniform, and the multi-beam area division process of full coverage from outside to inside comprises the following steps:
1) initializing the current scanning circle layer n to be 1, and inputting a parameter thetan=1、θ0pSetting a redundancy factor delta, presetting a stepping parameter delta of a scanning angle, wherein the delta is related to scanning precision;
2) firstly, the scanning angle alpha of the nth circle layer is determinednFinding a scanning angle alpha satisfying the following conditionn
Figure BDA0002858652890000031
And is
Figure BDA0002858652890000032
Wherein alpha isn=θnαp
Figure BDA0002858652890000036
Is alphanHalf pdB beamwidth;
Figure BDA0002858652890000033
is the scan angle (alpha)nHalf of-delta)
pdB beam width; 3) determining the number of covering beams N required by the nth layern
Figure BDA0002858652890000034
Wherein β is half of the adjacent beam spacing angle;
4) starting to prepare the coverage of the next circle layer, and updating the edge parameters
Figure BDA0002858652890000035
Judging whether theta is satisfiedn+1+Δ≤2θ0pIf yes, entering step 5); otherwise, updating n to n +1, and returning to the step 2);
5) further judge thetan+1The following three cases are classified into:
5-1) when 0<θn+1≤θ0pWhen the current area is divided into a whole area and a whole area, the last circle layer is covered by a beam positioned at the center of a circle, and the division of the whole area beam coverage area is finished;
5-2) when
Figure BDA0002858652890000041
Then, the last circle layer is covered by the beams with 4 beam tangent points at the circle center, and the division of the full-area beam coverage area is finished;
5-3) when
Figure BDA0002858652890000046
And then, covering the last two circle layers by using 7 beams, wherein the 2 nd from last circle layer uses 6 beams, the last circle layer is a beam positioned at the center of a circle, and the division of the full-area beam coverage area is finished.
That is, the multi-beam region dividing process is completed as long as the scanning angle α for each circle layer and the number N of beams for each layer are determined.
The redundancy factor delta is an important prior parameter, represents the redundancy meeting the full coverage, has important influence on the coverage result, and has the value range of delta
Figure BDA0002858652890000042
The quality of the coverage result can be observed by taking different delta values during specific operation.
Example 2
As shown in fig. 3, for the non-full coverage model, the coverage of the edge beams is also performed first, but each layer is not full coverage, but is planned as the beam tangency of each layer. The non-full coverage model is a special case of the full coverage model, and when Δ is 0, the full coverage is changed into the non-full coverage, and the two have uniformity.
Under the condition of full array, the process of dividing the multi-beam region from outside to inside by non-full coverage is as follows:
1) initializing the current scanning circle layer n to be 1, and inputting a parameter theta1、θ0pSetting a stepping parameter delta;
2) firstly, the scanning angle alpha of the nth circle layer is determinednFinding a scanning angle alpha satisfying the following conditionn
Figure BDA0002858652890000043
And is
Figure BDA0002858652890000044
3) Determining the number of covering beams N required by the nth layern
Figure BDA0002858652890000045
Wherein β is half of the adjacent beam spacing angle;
4) starting to prepare the coverage of the next circle layer, and updating the edge parameter thetan+1=αn-(θnn)=2αnnJudging whether theta is satisfiedn+1≤2θ0pIf yes, entering step 5); otherwise, updating n to n +1, and returning to the step 2);
5) further judge thetan+1The following two cases are classified into:
5-1) when 0<θn+1≤θ0pWhen the current area is divided into a whole area and a whole area, the last circle layer is covered by a beam positioned at the center of a circle, and the division of the whole area beam coverage area is finished;
5-2) when
Figure BDA0002858652890000051
And (4) returning to the step 2, and performing the last iteration.
Although the non-full coverage is less than 100%, when θ is equal to θ as shown in FIG. 3αpBeam theta changed to a dotted lineαxThen, 100% coverage, θ, can be achievedαxThe following relationship is satisfied:
Figure BDA0002858652890000052
the full coverage of the needed xdB beam width can be obtained.
Example 3
The multi-beam shaping method of the sub-array type can realize the shaping of the beam shapes of different areas under the condition of consistent channels. At this time, the number of times of use of the beam channel is required to be the same, and the subarrays are divided based on this precondition.
The subarray division rule is as follows:
layer 1 (outermost layer) is arranged according to full array N1Designing as N;
number of subarrays N of 2 nd circle layer2Is divided into two halves, as shown in FIG. 4, when N is an even number, N is2N/2; when N is an odd number, dividing the N into (N +1)/2 and (N-1)/2;
number of subarrays N of 3 rd circle layer3With N2Is divided into halves, as shown in FIG. 5, when N2When it is even, N3N 22; when N is present2Is divided into (N) when it is odd2+1)/2 and (N)2-1)/2。
Usually three hoop layers will suffice for the covering.
Because the same circle of beams has subarrays with different arrangement forms for covering, the value of the equivalent edge parameter α 'of each circle layer and the number N' of the covering beams of each circle layer need to be recalculated, and the specific process is as follows:
1) first, the coverage parameter [ alpha ] of the full array is calculated in the manner of example 1 or example 2123,…],[N1,N2,N3,…]And obtaining the equivalent half beam width theta 'of half pdB when the subarray beam points to the normal direction according to the antenna array topology as prior information before forming'0p
2) Initializing the current scanning circle layer n to be 1, and setting a stepping parameter deltaAnd a redundancy factor Δ; the coverage strategy of the circle 1 layer is not changed, the current circle 1 layer is assigned by utilizing the parameters of the circle 1 layer obtained in the mode of the embodiment 1 or the embodiment 2, and the equivalent edge parameter theta of the circle 1 layer1′=θ1Equivalent scan angle α'1=α1Number of equivalent coverage beams N'1=N1
3) The beam width of the subarray of the layer of the circle 2 is not changed much in one dimension compared with the beam width of the subarray of the circle 2, the beam width of the subarray of the other dimension is about twice as large as the beam width of the subarray of the circle 2, and the beam shape is an ellipse with the ratio of the long axis to the short axis being 2; the coverage area of the 2 nd circle is the positions of the 2 nd circle layer and the 3 rd circle layer of the full array strategy, then
Figure BDA0002858652890000061
And because the four subarrays are divided, then
Figure BDA0002858652890000062
Determining a scan angle α'2Half pdB beamwidth
Figure BDA0002858652890000063
Then calculating edge parameters
Figure BDA0002858652890000064
4) Preparing the subarray coverage of the layer of the 3 rd circle, which is equivalent to the 4 th to 6 th circles of the original coverage, and updating the equivalent edge parameters
Figure BDA0002858652890000065
Prepared from theta'3And theta'0pAnd (3) comparison:
when 0 is present<θ′3≤θ′0pThen, the last 1 central beam is needed to cover;
when in use
Figure BDA0002858652890000066
Then, the last 4 beams are needed to cover;
when in use
Figure BDA0002858652890000067
When, the last 7 beams are needed to cover.
Since the embodiment 3 uses the non-uniform beams, in order to solve the problem of color separation between the non-uniform beams, the embodiment uses a clustering color separation algorithm and a frequency multiplexing color separation algorithm based on depth-first search.
The clustering color separation algorithm has the following three characteristics: defining a distance measurement between beams as a similarity measurement between multiple beams; determining a criterion function for evaluating the quality of the clustering result; given a certain initial classification, the best clustering result of extreme values of the criterion function is obtained according to an iterative algorithm. The specific process is as follows:
distance matrix A defining the carrier-to-interference ratio C/I of a beamn×mWherein
Figure BDA0002858652890000068
Wherein n and m are the size of the array, PiFor a given beam power, i.e. the antenna gain, P, of the beamjFor a given co-channel beam power of a beam, a is expressed in dB.
Setting a classification target, wherein each color interval threshold value meets the condition that C/I is more than or equal to xi;
numbering all beams x1,x2,…,xnSelecting ω1={x1I.e. the beam # 1 falls under the first color ω1
According to the distance A between two beams21If A is21When x is more than xi, x is2Fall under omega1Otherwise, establishing new class omega2={x2};
When a step turns to xlWhen this is true, suppose that class k has already formed at that time, i.e., ω12,…,ωkCalculating the distance to each class, and taking the minimum distance min in a certain classkThen, the maximum max (min) of each class is takenk) If max (min)k) ≧ ζ, xlFall into the kth class, otherwise, x islEstablishing a new k + 1-th class, i.e., ωk+1={xk+1};
Based on the classification result, calculating
Figure BDA0002858652890000071
And obtaining final zeta 'to ensure that the relation C/I ≧ zeta' is satisfied among all the beams.
In one embodiment, the method further comprises the following steps: the color separation algorithm based on depth search converts the distance matrix A into an adjacency matrix G of an undirected graph, wherein
Figure BDA0002858652890000072
The physical significance of the method is that if C/I is larger than or equal to zeta, two points are not communicated, otherwise, the two points are communicated. That is, two beams smaller than the threshold are regarded as two points communicating with each other. Obviously, the two spots cannot be colored the same color, i.e., cannot be divided into co-frequency beams.
The problem is changed into a point coloring problem of a simple connected graph, the same frequency wave beams are colored in the same point, and the minimum point coloring number is the minimum color separation number. And (3) traversing the connected graph by adopting a depth-first search algorithm to obtain the optimal point coloring of the graph.
Simulation experiment
Take a rectangular grid planar array with an array size of 8 x 8 and a half-wavelength cell pitch as an example (all the shaped overlays mentioned below are developed based on this model).
Full coverage and non-full coverage multi-beam coverage models are established according to coverage scenes, and as shown in fig. 2 and 3, an outside-in region division strategy is provided. As shown in fig. 6, the 50 ° coverage area edge coverage results are presented.
The redundancy factor delta is an important prior parameter and characterizes the redundancy satisfying the full coverage, and different deltas have important influence on the coverage result. The value of the redundancy is not more than half of the width of the edge beam, otherwise, the redundancy is too large. Taking the example of covering 50 °, the result of covering with different Δ is shown in fig. 7. Δ ═ 0.5, NGeneral assembly=92;Δ=0.5,NGeneral assembly=92;Δ=1,NGeneral assembly=82;Δ=1.5,NGeneral assembly=78;Δ=2,NGeneral assembly=76;Δ=2.5,NGeneral assembly=77;Δ=3,NGeneral assembly=78;Δ=0,NGeneral assembly45. As Δ increases, NGeneral assemblyBecomes smaller but the number of layers required increases, and when Δ increases to a certain value, an inflection point appears, and when Δ continues to increase, NGeneral assemblyBut rather becomes larger. At this time, there is a critical optimum Δ, i.e., Δopt,Δopt=2。
In the ground mobile communication and high orbit satellite multi-beam, cellular networks, namely equal spot beam width models, are mostly adopted in the zoning strategy. The cells are arranged in regular hexagonal (cellular) order until the whole area is covered. The number of cells required for full coverage of 50 ° is 91, and 15 more beams are obtained compared with the number of beams obtained when the full coverage model Δ is 2 in this embodiment. The outside-in layered full coverage model reduces the number of multi-beam coverage.
For the multi-beam shaping method of the sub-array type, under the condition of consistent channels, the shaping of the beam shapes in different areas is realized. Taking the example of covering 50 °, the second turn is covered by a combination of forward and forward diagonal divisions as shown in fig. 4, and the result of the covering is shown in fig. 8. By combining the oblique division and the forward division, the beam coverage area has no obvious advantage, but the C/I problem is more advantageous, and the beam shape is closer to the target area. When a 50-degree area is covered, the coverage rate of the two areas is 100 percent; when covering a 55 ° area, the positive slope coverage is 87%, and the positive division coverage is 88%. Therefore, this method is not suitable for full coverage, which greatly increases the number of beams, contrary to the original intention of decreasing the number of beams in the sub-array.
Taking the subarray type multi-beam coverage as an example, according to the clustering and color separation algorithm, ζ is set to be 10dB, and 12 colors are required in total for the color separation result when the coverage is 50 °. When the coverage angles are 50 ° and 55 °, as the threshold ζ increases, the number of color separations required becomes larger; the same number of beams covers a larger range because the beams are further apart and require fewer separations.
If the minimum color separation number is obtained by directly adopting a color separation algorithm based on depth, the algorithm takes too long time, and an initial solution obtained by a clustering color separation algorithm is used as prior information to pre-classify all points (beams). The prior information is as visualized, all wave beams are divided into three groups, namely an outer ring, a middle ring and an inner ring, and then the optimal color separation number is obtained through a depth-first search algorithm.

Claims (6)

1. A method for dividing a low-orbit satellite multi-beam coverage area is characterized in that an outside-in area division strategy is adopted: covering the first edge scanning circle layer, determining the scanning angle and the number of covered beams of the current scanning circle layer to complete the beam covering of the current scanning circle layer, updating the edge parameters, and covering the scanning circle layer inwards layer by layer until the covering is completed.
2. The method of claim 1, wherein when the antenna array is a full array, the step of performing beam coverage comprises:
1) initializing the current scanning circle layer n to be 1, and setting an edge parameter thetan=1Half pdB beamwidth θ with beam pointing normal0pP is a preset gain flatness; setting a redundancy factor delta and presetting a stepping parameter delta of a scanning angle;
2) firstly, the scanning angle alpha of the nth circle layer is determinednFinding a scanning angle alpha satisfying the following conditionn
Figure FDA0002858652880000011
And is
Figure FDA0002858652880000012
Wherein alpha isn=θnαp
Figure FDA0002858652880000013
To a scanning angle alphanHalf pdB beamwidth;
Figure FDA0002858652880000014
is the scan angle (alpha)n- δ) half pdB beamwidth;
3) determining nth turn layer requirementsNumber of covered beams Nn
Figure FDA0002858652880000015
Wherein β is half of the adjacent beam spacing angle;
4) starting to prepare the coverage of the next circle layer, and updating the edge parameters
Figure FDA0002858652880000016
Judging whether theta is satisfiedn+1+Δ≤2θ0pIf yes, entering step 5); otherwise, updating n to n +1, and returning to the step 2);
5) further judge thetan+1The following three cases are classified into:
5-1) when 0<θn+1≤θ0pWhen the current area is divided into a whole area and a whole area, the last circle layer is covered by a beam positioned at the center of a circle, and the division of the whole area beam coverage area is finished;
5-2) when
Figure FDA0002858652880000017
Then, the last circle layer is covered by the beams with 4 beam tangent points at the circle center, and the division of the full-area beam coverage area is finished;
5-3) when
Figure FDA0002858652880000018
And then, covering the last two circle layers by using 7 beams, wherein the 2 nd from last circle layer uses 6 beams, the last circle layer is a beam positioned at the center of a circle, and the division of the full-area beam coverage area is finished.
3. The method of claim 2, wherein when the antenna array is full and is not full coverage, Δ ═ 0;
finding a scanning angle alpha satisfying the following condition in step 2) of completing beam coveragenComprises the following steps:
Figure FDA0002858652880000021
and is
Figure FDA0002858652880000022
Determining the number N of the covering beams needed by the nth circle layer in the step 3)n
Figure FDA0002858652880000023
Wherein β is half of the adjacent beam spacing angle;
the method for updating the edge parameters in the step 4) comprises the following steps: thetan+1=2αnn(ii) a Judging whether the conditions of the step 5) are carried out or not according to the following two conditions:
when 0 is present<θn+1≤θ0pWhen the current area is divided into a whole area and a whole area, the last circle layer is covered by a beam positioned at the center of a circle, and the division of the whole area beam coverage area is finished;
when in use
Figure FDA0002858652880000024
And returning to the step 2) for the last iteration.
4. The method of claim 2, wherein the number of subarrays in a full array antenna array is N, and when multi-beam forming in a subarray is used, the specific steps for completing beam coverage are as follows:
1) firstly, carrying out subarray division, wherein the division rule is as follows:
the outermost 1 st circle layer uses a full matrix N1Covering beams by N;
the number N of subarrays for beam coverage of the 2 nd circle layer2Dividing into halves, when N is even number, N2N/2; when N is an odd number, dividing the N into (N +1)/2 and (N-1)/2;
the number N of subarrays for the 3 rd circle layer to carry out beam covering3With N2Is divided into halves when N2When it is even, N3=N22; when N is present2Is divided into (N) when it is odd2+1)/2 and (N)2-1)/2;
2) Calculating the scanning angle alpha for completing the 1 st, 2 nd and 3 rd circle wave beam coverage when the antenna array is a full array123]The number of required coverage beams [ N ]1,N2,N3]And corresponding edge parameter [ theta ]123];
3) Firstly, obtaining the equivalent half-beam width theta of half pdB when the subarray beam points to the normal direction according to the antenna array topology0p(ii) a And then recalculating the equivalent scanning angles of the 1 st, 2 nd and 3 rd circle layers and the equivalent covering beam number to complete the multi-beam forming of the sub-array:
3-1) initializing a current scanning ring layer n to be 1, and setting a stepping parameter delta and a redundancy factor delta; the coverage strategy of the 1 st circle layer is the same as that of a full array, and the equivalent edge parameter theta of the 1 st circle layer1′=θ1Equivalent scan angle α'1=α1Number of equivalent coverage beams N'1=N1
3-2) the beam shape of the subarray of the 2 nd circle layer is an ellipse with the ratio of the long axis to the short axis being 2; calculating to obtain the equivalent scanning angle of the 2 nd circle layer
Figure FDA0002858652880000031
Number of equivalent coverage beams
Figure FDA0002858652880000032
Updating edge parameters
Figure FDA0002858652880000033
3-3) preparing the subarray coverage of the 3 rd circle layer, and updating the equivalent edge parameters
Figure FDA0002858652880000034
Then theta 'is prepared'3And theta'0pAnd (3) comparison:
when 0 is present<θ′3≤θ′0pWhen the temperature of the water is higher than the set temperature,covering the last circle layer by using a beam positioned at the center of a circle, and finishing the division of the full-area beam coverage area;
when in use
Figure FDA0002858652880000035
Then, the last circle layer is covered by the beams with 4 beam tangent points at the circle center, and the division of the full-area beam coverage area is finished;
when in use
Figure FDA0002858652880000036
And then, covering the last two circle layers by using 7 beams, wherein the 2 nd from last circle layer uses 6 beams, the last circle layer is a beam positioned at the center of a circle, and the division of the full-area beam coverage area is finished.
5. The method of claim 4, wherein when the sub-array multi-beam forming is used for beam covering, the color separation between non-uniform beams is performed by using a clustering color separation algorithm or a frequency multiplexing color separation algorithm based on depth-first search.
6. A method as claimed in claim 2, 3 or 4, characterized in that the redundancy factor Δ is over a range of values
Figure FDA0002858652880000037
CN202011553560.1A 2020-12-24 2020-12-24 Low-orbit satellite multi-beam coverage area dividing method Active CN112672423B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011553560.1A CN112672423B (en) 2020-12-24 2020-12-24 Low-orbit satellite multi-beam coverage area dividing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011553560.1A CN112672423B (en) 2020-12-24 2020-12-24 Low-orbit satellite multi-beam coverage area dividing method

Publications (2)

Publication Number Publication Date
CN112672423A true CN112672423A (en) 2021-04-16
CN112672423B CN112672423B (en) 2022-05-03

Family

ID=75408527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011553560.1A Active CN112672423B (en) 2020-12-24 2020-12-24 Low-orbit satellite multi-beam coverage area dividing method

Country Status (1)

Country Link
CN (1) CN112672423B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113825238A (en) * 2021-09-03 2021-12-21 天地信息网络研究院(安徽)有限公司 Flexible signaling beam multi-scene application method in satellite communication
CN113872653A (en) * 2021-09-30 2021-12-31 重庆两江卫星移动通信有限公司 Beam forming method based on earth matching
CN114389675A (en) * 2021-12-24 2022-04-22 军事科学院系统工程研究院网络信息研究所 Multi-satellite multi-beam combined forming coverage area design method
CN115038095A (en) * 2022-06-15 2022-09-09 成都德杉科技有限公司 Multi-beam coverage planning method for low-earth orbit satellite
CN115242298A (en) * 2022-09-22 2022-10-25 鹏城实验室 Method and system for simulating satellite beam coverage area based on space geometric model

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107086888A (en) * 2017-03-02 2017-08-22 重庆邮电大学 A kind of two-layer hybrid satellite network optimization design and its covering performance appraisal procedure
CN107872348A (en) * 2017-09-30 2018-04-03 北京邮电大学 It is a kind of to divide domain hierarchy control method, device and system towards Incorporate network
CN109104236A (en) * 2018-07-18 2018-12-28 北京天链测控技术有限公司 A kind of method that low orbit satellite is communicated with grounded receiving station
US10419106B1 (en) * 2017-09-29 2019-09-17 Amazon Technologies, Inc. Satellite-based content delivery network (CDN) in an extraterrestrial environment
US20200077411A1 (en) * 2018-08-29 2020-03-05 Qualcomm Incorporated Elevation restriction beamforming in wireless systems
WO2020051508A1 (en) * 2018-09-06 2020-03-12 Lynk Global, Inc. Cellular core network and radio access network infrastructure and management in space
CN111083754A (en) * 2019-11-27 2020-04-28 成都天奥集团有限公司 Method for maintaining terminal tracking area irrelevant to cell switching under low-earth-orbit satellite communication system
CN111757267A (en) * 2019-03-29 2020-10-09 华为技术有限公司 Satellite tracking area updating method and related device

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107086888A (en) * 2017-03-02 2017-08-22 重庆邮电大学 A kind of two-layer hybrid satellite network optimization design and its covering performance appraisal procedure
US10419106B1 (en) * 2017-09-29 2019-09-17 Amazon Technologies, Inc. Satellite-based content delivery network (CDN) in an extraterrestrial environment
CN107872348A (en) * 2017-09-30 2018-04-03 北京邮电大学 It is a kind of to divide domain hierarchy control method, device and system towards Incorporate network
CN109104236A (en) * 2018-07-18 2018-12-28 北京天链测控技术有限公司 A kind of method that low orbit satellite is communicated with grounded receiving station
US20200077411A1 (en) * 2018-08-29 2020-03-05 Qualcomm Incorporated Elevation restriction beamforming in wireless systems
WO2020046484A1 (en) * 2018-08-29 2020-03-05 Qualcomm Incorporated Elevation restriction beamforming in wireless systems
WO2020051508A1 (en) * 2018-09-06 2020-03-12 Lynk Global, Inc. Cellular core network and radio access network infrastructure and management in space
CN111757267A (en) * 2019-03-29 2020-10-09 华为技术有限公司 Satellite tracking area updating method and related device
CN111083754A (en) * 2019-11-27 2020-04-28 成都天奥集团有限公司 Method for maintaining terminal tracking area irrelevant to cell switching under low-earth-orbit satellite communication system

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
KENNETH N. SHERMAN: "Phased Array Shaped Multi-Beam Optimization for LEO Satellite Communications Using a Genetic Algorithm", 《IEEE》 *
MENG SHAOFEI: "Analysis of detection capabilities of LEO reconnaissance satellite constellation based on coverage performance", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 *
SHKELZEN CAKAJ: "The Coverage Analysis for Low Earth Orbiting Satellites at Low Elevation", 《INTERNATIONAL JOURNAL OF ADVANCED COMPUTER SCIENCE AND APPLICATIONS》 *
ZHAOFENG WU: "A Graph-Based Satellite Handover Framework for LEO Satellite Communication Networks", 《IEEE COMMUNICATIONS LETTERS》 *
梁广: "星载多波束相控阵天线等通量覆盖设计", 《电波科学学报》 *
贾铂奇: "阵列天线多波束赋形技术研究", 《信息科技辑》 *
郝才勇: "卫星通信近期发展综述", 《电子技术应用》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113825238A (en) * 2021-09-03 2021-12-21 天地信息网络研究院(安徽)有限公司 Flexible signaling beam multi-scene application method in satellite communication
CN113825238B (en) * 2021-09-03 2024-02-20 天地信息网络研究院(安徽)有限公司 Flexible signaling beam multi-scene application method in satellite communication
CN113872653A (en) * 2021-09-30 2021-12-31 重庆两江卫星移动通信有限公司 Beam forming method based on earth matching
CN114389675A (en) * 2021-12-24 2022-04-22 军事科学院系统工程研究院网络信息研究所 Multi-satellite multi-beam combined forming coverage area design method
CN114389675B (en) * 2021-12-24 2022-09-20 军事科学院系统工程研究院网络信息研究所 Multi-satellite multi-beam combined forming coverage area design method
CN115038095A (en) * 2022-06-15 2022-09-09 成都德杉科技有限公司 Multi-beam coverage planning method for low-earth orbit satellite
CN115242298A (en) * 2022-09-22 2022-10-25 鹏城实验室 Method and system for simulating satellite beam coverage area based on space geometric model
CN115242298B (en) * 2022-09-22 2022-12-23 鹏城实验室 Method and system for simulating satellite beam coverage area based on space geometric model

Also Published As

Publication number Publication date
CN112672423B (en) 2022-05-03

Similar Documents

Publication Publication Date Title
CN112672423B (en) Low-orbit satellite multi-beam coverage area dividing method
Thornton et al. Optimizing an array of antennas for cellular coverage from a high altitude platform
Fonseca et al. Multi-beam reflector antenna system combining beam hopping and size reduction of effectively used spots
CN103022728B (en) Method for designing partial-feed paraboloid multi-beam antenna feed array
CN112636804B (en) Multi-unmanned aerial vehicle base station cooperative transmission method based on millimeter wave array
US9544029B2 (en) Method for antenna array with steerable spotlight beams
CN107732464B (en) Design method, system and medium of multivariable shaped beam antenna
CN104103910A (en) Optimization design method of single-aperture and multi-beam antenna
CN106252901A (en) Wideband three beam array antenna
US8145123B2 (en) Satellite communication system and method for dividing the coverage area thereof
CN109004366A (en) A kind of antenna assembly and the method for wave beam adjustment
CN108511914A (en) Coated by dielectric realizes the convergent device and method of wide bandwidth base station antenna wave beam
Park et al. Technology development for wireless communications system using stratospheric platform in Korea
CN116633404A (en) Combined beam forming and deployment method for double intelligent reflecting surface auxiliary communication
Chang et al. Quantization effects of phase shifters on 5G mmWave antenna arrays
CN109103610B (en) Multi-beam antenna with non-uniform sub-beam coverage and design method
CN108668352B (en) Ocean communication QoS guarantee method based on multi-beam power distribution
Lai et al. A multipart 5G base-station antenna using series-fed patch antenna sub-arrays
Dicandia et al. Improving the 5G massive MIMO performance for aerial base stations by exploiting triangular lattice arrays
Xu et al. Optimisation of beam pattern of high-altitude platform antenna using conventional beamforming
Wu et al. Towards ubiquitous coverage of high altitude platforms aided 5G+ for massive Internet of Things: A cell perspective
Tao et al. Shape optimization design method of high performance Multi-beam antenna for mixed-payload communication satellite
Zhang et al. Generalized cell planning technique applied for 4x12 frequency-reuse
CN115588851B (en) Validity verification method for array element channel phase calibration of spherical phased array antenna
KR102189953B1 (en) Design method for micro-strip antenna with unit metal patches

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant