CN110278024B - System capacity optimization method and device for satellite communication constellation - Google Patents

System capacity optimization method and device for satellite communication constellation Download PDF

Info

Publication number
CN110278024B
CN110278024B CN201910727794.4A CN201910727794A CN110278024B CN 110278024 B CN110278024 B CN 110278024B CN 201910727794 A CN201910727794 A CN 201910727794A CN 110278024 B CN110278024 B CN 110278024B
Authority
CN
China
Prior art keywords
optimal solution
interference
transmitting power
mapping relation
satellite communication
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.)
Active
Application number
CN201910727794.4A
Other languages
Chinese (zh)
Other versions
CN110278024A (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.)
Shanghai Qingshen Technology Development Co ltd
Tsinghua University
Original Assignee
Shanghai Qingshen Technology Development Co ltd
Tsinghua University
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 Shanghai Qingshen Technology Development Co ltd, Tsinghua University filed Critical Shanghai Qingshen Technology Development Co ltd
Priority to CN201910727794.4A priority Critical patent/CN110278024B/en
Publication of CN110278024A publication Critical patent/CN110278024A/en
Application granted granted Critical
Publication of CN110278024B publication Critical patent/CN110278024B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay
    • H04B7/18513—Transmission in a satellite or space-based system
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay
    • H04B7/18519—Operations control, administration or maintenance
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Relay Systems (AREA)

Abstract

The invention provides a method and a device for optimizing system capacity of a satellite communication constellation, which relate to the technical field of satellite communication, wherein the satellite communication constellation comprises a plurality of non-geosynchronous orbit satellites, and the method comprises the following steps: obtaining constraint conditions of system capacity performance in a satellite communication constellation, wherein the constraint conditions comprise constraint conditions of a mapping relation between beams and an earth station, constraint conditions of transmitting power and constraint conditions of interference values, and the interference values are represented by the mapping relation and the transmitting power; calculating according to the constraint conditions to obtain an optimal solution of the transmitting power and an optimal solution of the mapping relation; and carrying out communication according to the optimal solution of the transmitting power and the optimal solution of the mapping relation, and ensuring that the maximum system capacity is achieved under the condition of not interfering the same-frequency geosynchronous orbit satellite communication system.

Description

System capacity optimization method and device for satellite communication constellation
Technical Field
The invention relates to the technical field of satellite communication, in particular to a method and a device for optimizing system capacity of a satellite communication constellation.
Background
The Satellite communication system includes a geosynchronous orbit (GSO) Satellite communication system and a non-geosynchronous orbit (NGSO) Satellite communication system, classified according to Satellite orbits. Compared with GSO satellite, the NGSO satellite system positioned in middle and low orbit has the advantages of short propagation delay, small link loss and the like. Because the coverage range of a single NGSO satellite is very limited, a constellation can be formed through the orbit design to realize global and full-time coverage.
With the increase of NGSO satellites, there are more and more satellites in the same frequency and adjacent orbits, which makes the problem of co-channel interference between the NGSO constellation and the GSO satellite system increasingly serious. In addition, in practical applications, the situation of resource maldistribution often occurs because of the need to consider more variables and more complex modulation modes, so that the NGSO constellation cannot achieve the optimal system capacity.
Disclosure of Invention
The invention aims to provide a method and a device for optimizing system capacity of a satellite communication constellation, which ensure that the maximum system capacity is achieved under the condition of not interfering a same-frequency geosynchronous orbit communication system.
In a first aspect, an embodiment of the present invention provides a method for optimizing system capacity of a satellite communication constellation, where the satellite communication constellation includes a plurality of non-geosynchronous orbit satellites, and the method includes:
obtaining constraint conditions of system capacity performance in the satellite communication constellation, wherein the constraint conditions comprise constraint conditions of a mapping relation between beams and earth stations, constraint conditions of transmitting power and constraint conditions of interference values, and the interference values are characterized by the mapping relation and the transmitting power;
calculating according to the constraint condition to obtain an optimal solution of the transmitting power and an optimal solution of the mapping relation;
and communicating according to the optimal solution of the transmitting power and the optimal solution of the mapping relation.
In an optional embodiment, the constraint condition includes that a preset mapping relationship that the beam and the earth station need to satisfy, the interference value needs to be smaller than an interference threshold value, the transmission power needs to be within a power threshold range, and a time-space characteristic variation range of each satellite in the satellite communication constellation corresponds to a preset orbit configuration.
In an optional implementation manner, the calculating the optimal solution of the transmission power and the optimal solution of the mapping relationship according to the constraint condition includes:
carrying out solving calculation according to the interference characteristic and the constraint characteristic to obtain the transmitting power and the mapping relation on each link;
and carrying out iterative calculation on the transmitting power and the mapping relation on each link to respectively obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation.
In an optional embodiment, performing solution calculation according to the interference characteristic and the constraint characteristic to obtain the transmission power and the mapping relationship on each link includes:
and calculating the interference value and the interference threshold value by accumulating and representing the interference value by the product of the transmitting power and the mapping relation of the plurality of non-geosynchronous earth orbit satellites to obtain the transmitting power and the mapping relation on each link which meet the requirement of communication service.
In an alternative embodiment, the method further comprises:
and calculating to obtain the maximum throughput according to the optimal solution of the transmitting power and the optimal solution of the mapping relation.
In an optional implementation manner, calculating the maximum throughput according to the optimal solution of the transmission power and the optimal solution of the mapping relationship, further includes:
the maximum throughput is calculated according to the following formula:
Figure BDA0002158223730000031
wherein, ai,jRepresenting the mapping of the ith beam to the jth earth station, pi,jRepresenting the magnitude of the transmitted power of the ith beam to the jth earth station, ci,jRepresenting the maximum communication rate per link, qi,jIs the product coefficient, IthIs the interference threshold.
In an alternative embodiment, the method further comprises:
and adjusting the modulation mode of the satellite communication constellation.
In an optional embodiment, the modulation scheme includes one or more of an adaptive modulation and coding scheme, a pulse code modulation scheme, a phase keying modulation scheme, a quadrature amplitude modulation scheme, and a minimum frequency shift keying modulation scheme.
In an optional implementation manner, before the calculating the optimal solution of the transmission power and the optimal solution of the mapping relation according to the constraint condition, the method further includes:
and (4) performing convex optimization problem transformation operation to enable the transformed convex optimization problem to meet the standard form of the convex optimization problem.
In a second aspect, an embodiment of the present invention provides an apparatus for optimizing system capacity of a satellite communication constellation, where the satellite communication constellation includes a plurality of non-geosynchronous orbit satellites, and the apparatus includes:
an obtaining module, configured to obtain constraint conditions of system capacity performance in the satellite communication constellation, where the constraint conditions include a constraint condition of a mapping relationship between a beam and an earth station, a constraint condition of transmit power, and a constraint condition of an interference value, and the interference value is characterized by the mapping relationship and the transmit power;
the calculation module is used for calculating to obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation according to the constraint condition;
and the communication module is used for carrying out communication according to the optimal solution of the transmitting power and the optimal solution of the mapping relation.
The embodiment of the invention provides a method and a device for optimizing system capacity of a satellite communication constellation, wherein the satellite communication constellation comprises a plurality of non-geosynchronous orbit satellites, the optimal solution of transmitting power and the optimal solution of the mapping relation with an earth station are calculated according to the constraint condition of the system capacity performance in the satellite communication constellation, the joint optimization of the transmitting power and the mapping relation is realized, and the maximum system capacity can be achieved under the condition that the same-frequency geosynchronous orbit communication system is not interfered.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of a method for optimizing system capacity of a satellite communication constellation according to an embodiment of the present invention;
fig. 2 is a flowchart of a system capacity optimization method for a satellite communication constellation according to another embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating a system throughput curve versus variation of a non-geosynchronous orbit satellite communication constellation according to an embodiment of the present invention;
fig. 4 is a schematic diagram illustrating comparison and variation of a receiving-end lumped interference curve of a non-geosynchronous orbit satellite communication constellation according to an embodiment of the present invention;
FIG. 5 is a second schematic diagram illustrating a system throughput curve versus variation of a non-geosynchronous orbit satellite communication constellation according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a comparison variation of a receiving-end lumped interference curve of a non-geosynchronous orbit satellite communication constellation according to an embodiment of the present invention;
fig. 7 is a functional block diagram of an apparatus for optimizing system capacity of a satellite communication constellation according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
With the increase of NGSO satellites, there are more and more satellites in the same frequency and adjacent orbits, which makes the problem of co-channel interference between the NGSO constellation and the GSO satellite system increasingly serious. According to the ITU regulations, the medium and low orbit satellites must take measures to avoid harmful interference on the GSO satellite during communication, and the transmitting power is controlled on the premise of meeting the communication requirements of the medium and low orbit satellites, so that potential co-frequency interference can be effectively avoided. Currently, there is little research on power control, and the main control schemes are: adaptive power control technology, downlink adaptive FEC coding technology, feedback loop power control technology, power distribution technology and the like. These schemes basically monitor the link quality at the receiving end of the NGSO satellite system, such as the signal-to-noise ratio, the equivalent power flux density, etc., and adjust the transmission power according to a certain rule or method to avoid potential interference when co-channel interference may occur.
Since these schemes have certain difficulty in engineering implementation, in practice, it is common practice to turn off the transmitter thereof to avoid interference, and for other non-interfering links, the corresponding transmitter is usually in the maximum power transmission state. Some of the links currently will not continue to provide communication services to earth stations because the transmitter is turned off. Therefore, to ensure continuity of communication, the earth station is usually under multi-satellite coverage, and when a certain link cannot continue communication, the earth station will access a new NGSO satellite link to ensure continuity of communication. At this time, the system performance such as call blocking rate, switching overhead, etc. is directly affected by what rule to access the NGSO satellite system. In consideration of the characteristics of the NGSO satellite system communication, the currently common multi-satellite coverage access strategies include a distance priority strategy, a coverage time priority strategy, a maximum elevation angle priority strategy and the like. In addition, due to the problems of fluctuation of the earth surface topography, change of the atmospheric environment, ground object shielding and the like, the modulation mode needs to be adjusted according to the actual situation during communication so as to ensure the communication quality under different environments.
Therefore, in order to ensure that the NGSO communication system is continuous and stable and does not interfere with other GSO communication systems, many variables are often considered in design, and meanwhile, various control strategies such as power control, multi-satellite access, modulation modes and the like are complex, so that on-satellite resource distribution is likely to be uneven, and the performance of the whole NGSO constellation system cannot reach the best.
Based on this, the method and the device for optimizing the system capacity of the satellite communication constellation provided by the embodiment of the invention can ensure that the maximum system capacity is achieved under the condition that the same-frequency geosynchronous orbit communication system is not interfered.
To facilitate understanding of the present embodiment, first, a method for optimizing system capacity of a satellite communication constellation disclosed in the embodiment of the present invention is described in detail.
Fig. 1 is a flowchart of a method for optimizing system capacity of a satellite communication constellation according to an embodiment of the present invention.
Referring to fig. 1, a method for optimizing system capacity of a satellite communication constellation is applied to the technical field of satellite communication systems in mobile communication, the satellite communication constellation comprises a plurality of non-geosynchronous orbit satellites, and the method comprises the following steps:
step S102, obtaining constraint conditions of system capacity performance in a satellite communication constellation, wherein the constraint conditions comprise constraint conditions of a mapping relation between beams and an earth station, constraint conditions of transmitting power and constraint conditions of interference values, and the interference values are represented through the mapping relation and the transmitting power.
And step S104, calculating according to the constraint conditions to obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation.
And S106, communicating according to the optimal solution of the transmitting power and the optimal solution of the mapping relation.
In the practical application preferred embodiment, the optimal solution of the transmitting power and the optimal solution of the mapping relation with the earth station are calculated according to the constraint condition of the system capacity performance in the satellite communication constellation, so that the joint optimization of the transmitting power and the mapping relation is realized, and the maximum system capacity can be achieved under the condition that the same-frequency geosynchronous orbit communication system is not interfered.
The constraint conditions comprise a preset mapping relation which is required to be met by a beam and an earth station, an interference value which is required to be smaller than an interference threshold value, transmission power which is required to be within a power threshold value range, and a space-time characteristic change range of each satellite in a satellite communication constellation which corresponds to a preset orbit configuration.
It should be noted that, the interference avoidance constraint of the NGSO satellite system is specified according to the existing protocol framework of the international telecommunications union, the NGSO satellite system must not generate harmful interference to the GSO satellite system, the interference evaluation index (interference threshold value) is usually the lumped interference-to-noise ratio or the equivalent power flux density at the receiving end of the GSO satellite system, and different interference protection standards are set for different frequency bands and different interference evaluation indexes, which is not limited herein.
The transmit power constraint of the NGSO satellite system, in order to avoid interference to other systems in the space, and simultaneously meet the communication quality requirement of the system itself, the transmit power P of the NGSO system needs to be limited in a certain interval: wherein the upper limit PmaxIs the maximum value, lower limit P, that the transmitter can transmitminThe minimum value of the transmitting power of the carrier-to-noise ratio threshold of the receiving end is met so as to ensure the quality of the link of the receiving end; when the transmitting power P is equal to zero, the transmitter is turned off and stops transmitting.
And (3) the orbital configuration of the NGSO satellite system is restricted, namely the NGSO satellite system needs to run on a set orbit, and the space-time characteristic change range of each NGSO satellite needs to correspond to the orbital configuration preset by the NGSO satellite.
In order to obtain the optimal solution calculation method of the transmission power and the mapping relationship more clearly, step S104 in the above embodiment includes:
1. and solving and calculating according to the interference characteristic and the constraint characteristic to obtain the transmitting power and the mapping relation on each link.
The interference value and the interference threshold value are calculated by accumulating and representing the interference value through the product of the transmitting power and the mapping relation of a plurality of non-geosynchronous orbit satellites, and the transmitting power and the mapping relation on each link which meet the communication service requirement are obtained.
2. And carrying out iterative calculation on the transmitting power and the mapping relation on each link to respectively obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation.
Further, the method further comprises: and calculating to obtain the maximum throughput according to the optimal solution of the transmitting power and the optimal solution of the mapping relation, and optimizing the system capacity according to the maximum throughput.
The maximum throughput is calculated according to the following formula:
Figure BDA0002158223730000081
wherein, ai,jRepresenting the mapping of the ith beam to the jth earth station, pi,jRepresenting the magnitude of the transmitted power of the ith beam to the jth earth station, ci,jRepresenting the maximum communication rate per link, qi,jIs the product coefficient, IthIs the interference threshold.
The embodiment of the invention needs to simultaneously meet the constraint conditions, and on the basis, the system resources are distributed to improve the system capacity through the joint optimization of the transmitting power and the mapping relation.
Considering a scene of co-frequency coexistence of an NGSO satellite base system and a GSO satellite system, the NGSO system consists of S NGSO satellites and N earth stations, and each NGSO satellite is provided with K electric scanning phased array beams.
Defining a transmit power matrix P ═ P for an NGSO constellation systemi,j](S·K)×NMapping relation matrix A of beams and earth stations on NGSO satellite star ═ ai,j](S·K)×NWherein p isi,jRepresenting the magnitude of the transmitted power of the ith beam to the jth earth station, each pi,jSatisfy Pmin≤pi,j≤PmaxOr pi,j＝0；ai,jRepresenting the mapping relation between the ith beam and the jth earth station, if the ith beam of the NGSO constellation establishes a link with the jth NGSO earth station, corresponding to ai,jThe value is 1, otherwise the value is 0. The specific implementation method is shown in fig. 2, and comprises the following steps:
step S201, determining an optimization target, that is, a certain system performance of the NGSO constellation system, which is not limited herein, and in the embodiment of the present invention, the system capacity performance may be taken as an example, setting a corresponding constraint condition according to the system capacity performance, and constructing a constraint equation satisfying the constraint condition by using the constraint condition, where the constraint equation includes an interference constraint equation, a constraint equation for the power matrix P, and a constraint equation for the beam mapping relationship matrix a.
Here, in order to calculate the optimal solution of the transmission power and the mapping relationship more accurately, before calculating the optimal solution of the transmission power and the optimal solution of the mapping relationship according to the constraint condition, the method generally further includes:
and (4) performing convex optimization problem transformation operation to enable the transformed convex optimization problem to meet the standard form of the convex optimization problem.
Step S202, properly converting the system performance and a constraint equation to enable the converted system performance and the constraint equation to meet the standard form of a convex optimization problem;
step S203, solving the convex optimization problem, analyzing the interference characteristics and constraint characteristics of all potential links between the satellite and the earth station, and respectively obtaining the transmitting power p on each potential linki,jAnd mapping relation a of beams on NGSO satellite and earth stationi,j。
Step S204, iterative computation pi,jAnd ai,jOf (2) an optimal solution
Figure BDA0002158223730000091
And
Figure BDA0002158223730000092
obtaining an optimal solution matrix PoptAnd Aopt。
In the solving process, each potential link transmitting power pi,jIn the matrix P0Obtaining a transmitting power matrix P after optimizing the interference characteristic and the constraint characteristic1Wherein the matrix P1Determining a mapping relation matrix A1By mapping the relationship matrix A1To optimize the transmission power matrix P1Obtaining a new transmit matrix P2，P2And determining to obtain a mapping relation matrix A2By analogy, continuously iterating each other until convergence to obtain the final
Figure BDA0002158223730000101
And
Figure BDA0002158223730000102
step S205, adjusting the modulation mode according to the propagation state and link margin on each link, and maximizing the link rate under the condition of meeting the margin required by engineering, so as to enable the system capacity CoptAnd the optimization is achieved.
As a possible embodiment, the system capacity index is considered as the throughput T in this examplem. Simulation parameters of the NGSO and GSO constellation satellite systems are respectively shown in the following tables 1 and 2:
TABLE 1
Figure BDA0002158223730000104
TABLE 2
Figure BDA0002158223730000103
Figure BDA0002158223730000111
Defining a transmit power matrix P ═ Pi,j](S·K)×NWherein S, K, N is given by the above table, each pi,jSatisfies 0 ≤ pi,j≤Pmax。
Defining a beam link allocation matrix a ═ ai,j](S·K)×NIs provided with ai,j∈ {0,1} represents the mapping of the ith beam of the NGSO constellation to the jth NGSO earth station, assuming that a beam can only establish a communication link with at most one earth station at a time, and vice versa, and thus the beam allocation element ai,jSatisfies the following conditions:
Figure BDA0002158223730000112
by the above definition, we can convert the throughput TmRepresented by the formula:
Figure BDA0002158223730000113
wherein, ci,jRepresenting the maximum communication rate of each link.
Because the same communication frequency point is used by the NGSO constellation system and the GSO satellite system, when the ith wave beam of the NGSO constellation and the jth earth station are establishedDuring communication link, the ith wave beam of NGSO constellation will generate co-frequency interference to the receiving end of GSO earth station, and the corresponding interference value Ii,jCan be expressed as:
Figure BDA0002158223730000114
wherein G isngt(θ1) For the transmit gain of the NGSO on-satellite beam in the direction of the GSO earth station, Ggr(θ2) For the GSO earth station in the direction of the NGSO satellite in which the ith beam is located, di,jIs the distance, p ', from the NGSO satellite where the ith beam is located to the GSO earth station'i,jTransmitting power p for NGSO beamsi,jConverting the equivalent transmitting power to the overlapped frequency band, wherein the relation is as follows:
Figure BDA0002158223730000121
wherein, WngFor the downstream communication bandwidth, W, of the NGSO constellation systemgIs the downstream communication bandwidth of GSO constellation system, therefore, the interference value Ii,jAnd a transmission power pi,jCan be expressed in the form:
Figure BDA0002158223730000122
wherein q isi,jFor coefficients, the lumped interference I received by the GSO earth station from the NGSO constellation is calculated as follows:
Figure BDA0002158223730000123
in summary, the NGSO constellation system should satisfy the following constraints:
system transmission power pi,jSatisfies 0 ≤ pi,j≤Pmax。
A beam of the system can only establish a communication link with one earth station and vice versa.
The system should satisfy the interference avoidance condition and have low interference valueIn the interference discrimination criterion Ith(the interference to noise ratio threshold is-12.2 dB).
Applying a joint optimization method to obtain the system performance throughput TmThe specific steps are as follows:
firstly, an optimization problem and a constraint equation of the problem are constructed, namely:
Figure BDA0002158223730000131
the problem is converted into a convex optimization problem, namely corresponding conditions of a model are converted, so that the convex optimization problem meets the standard form of the convex optimization problem. The transformation process is as follows: the objective function being to find the throughput TmTo convert this maximization problem into a minimization problem; optimizing the variable as the transmit power pi,jAnd beam link ai,jWherein p isi,jIs a continuous variable, ai,jFor non-continuous variables, 0-1 non-continuous variable ai,jRelaxed to a value range of [0, 1]The end result is:
Figure BDA0002158223730000132
solving the problem: lagrangian functions are listed:
Figure BDA0002158223730000133
wherein λ, μ, η are lagrangian multipliers, the lagrangian multipliers are initialized, and the problem is solved by using a lagrangian dual decomposition method, wherein the dual function of the original problem is as follows:
Figure BDA0002158223730000149
decomposing the dual function into a main problem and (S.K) × N sub-problems, wherein the optimal solution of each sub-problem is solved
Figure BDA0002158223730000141
And
Figure BDA0002158223730000142
satisfying the optimization condition of Karush-Kuhn-Tucker (KKT), and iteratively calculating pi,jAnd ai,jOf (2) an optimal solution
Figure BDA0002158223730000143
And
Figure BDA0002158223730000144
in the iterative calculation
Figure BDA0002158223730000145
In the process of (3), because the constraint conditions C2 and C3 have the mutual coupling and restriction relationship, a dichotomy model in a graph theory is introduced, the problem is solved and converted into a dichotomy graph minimum weight matching problem, the Kuhn-Munkres algorithm is utilized to complete minimum weight matching, and the calculation is updated to obtain
Figure BDA0002158223730000146
Repeating iteration until convergence, and calculating the optimal solution of the convex optimization problem
Figure BDA0002158223730000147
And
Figure BDA0002158223730000148
obtaining an optimal solution matrix PoptAnd AoptAnd calculating the throughput TmAn optimal value;
based on the calculation of the optimal value of the throughput, the modulation mode of the satellite communication constellation can be readjusted to ensure that the throughput T ismA better value is reached.
The modulation method includes, but is not limited to, one or more of an adaptive modulation and coding method (ACM), a pulse code modulation method (PCM), a phase shift keying modulation method (PSK/DPSK/QPSK), a quadrature amplitude modulation method (QAM), and a minimum shift keying modulation Method (MSK).
Aiming at the NGSO constellation system under the GSO same-frequency coexistence scene, the transmission power in the satellite which avoids the interference and the distribution of the optimized beam link are realized through a joint optimization mode, and the system resources such as a modulation mode and the like are adjusted to improve the whole capacity of the system.
As an optional embodiment, if the NGSO constellation system is configured to ensure that the normal communication service of the NGSO constellation system is performed because the interference value exceeds the interference threshold, the corresponding link is switched on and off, so that the mapping relationship or the allocated power transmission changes, and at this time, the performance of the system is optimized again by the joint optimization method of the transmission power and the mapping relationship provided by the embodiment of the present invention.
To highlight the optimization performance of the real-time optimization method of the present invention, the embodiment of the present invention is compared and analyzed with the conventional allocation strategies 1 and 2, respectively, wherein the allocation strategies 1 and 2 are described as follows:
allocation policy 1: the current NGSO constellation system adopts a space-domain isolation mode to avoid harmful co-channel interference on GSO, a link separation angle is set to be 6 degrees, and satellite beam resources are distributed according to the maximum elevation angle rule from the angle of an earth station on the basis.
Allocation policy 2: adopting an interference avoidance strategy which is the same as the distribution strategy 1, and setting a space domain isolation with a link separation angle of 6 degrees; under this constraint, at the link communication rate ci,jAnd as a weight value, solving the maximum weight distribution problem by adopting a KM algorithm, and completing corresponding satellite beam resource distribution according to the solving result.
As can be seen from fig. 4 and 6, the harmful interference of the NGSO constellation system to the GSO satellite system can be effectively avoided by all of the three methods. As shown in fig. 3, when the NGSO constellation system uses the method of the embodiment of the present invention, the corresponding system throughput is entirely higher than the system throughput when the allocation policy 1 is used, which is higher by 6.1% on average; as shown in fig. 5 and 6, the difference between the corresponding NGSO constellation system throughputs is small in the method and the allocation strategy 2 according to the embodiment of the present invention under the condition that the lumped interference constraint is weak; but when the orbit condition is worse, the lumped interference constraint received by the GSO earth station becomes stronger, and compared with the distribution strategy 2, the worst value of the throughput is improved by 11.6 percent.
From the above analysis, it can be seen that, by adjusting the transmission power of the transmitter in the NGSO constellation system, the mapping relationship between the beams on the NGSO satellite and the earth station, and adjusting the modulation mode, the embodiments of the present invention can ensure that resource allocation is optimized without interfering with the same-frequency geosynchronous orbit communication system, so as to optimize the system capacity.
As shown in fig. 7, an embodiment of the present invention further provides an apparatus for optimizing system capacity of a satellite communication constellation, where the satellite communication constellation includes a plurality of non-geosynchronous orbit satellites, and the apparatus includes:
the system comprises an acquisition module, a receiving module and a processing module, wherein the acquisition module is used for acquiring constraint conditions of system capacity performance in a satellite communication constellation, the constraint conditions comprise constraint conditions of a mapping relation between a beam and an earth station, constraint conditions of transmitting power and constraint conditions of an interference value, and the interference value is represented by the mapping relation and the transmitting power;
the calculation module is used for calculating according to the constraint condition to obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation;
and the communication module is used for carrying out communication according to the optimal solution of the transmitting power and the optimal solution of the mapping relation.
The system capacity optimization device for the satellite communication constellation provided by the embodiment of the invention has the same technical characteristics as the system capacity optimization method for the satellite communication constellation provided by the embodiment, so that the same technical problems can be solved, and the same technical effect can be achieved.
The computer program product of the method and apparatus for optimizing system capacity of a satellite communication constellation provided in the embodiments of the present invention includes a computer-readable storage medium storing a program code, where instructions included in the program code may be used to execute the method described in the foregoing method embodiments, and specific implementation may refer to the method embodiments, and will not be described herein again.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In addition, in the description of the embodiments of the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
The embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the steps of the method for optimizing the system capacity of the satellite communication constellation provided in the foregoing embodiment are implemented.
The embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method for optimizing the system capacity of the satellite communication constellation according to the embodiment are executed.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein.

Claims (8)

1. A method for optimizing system capacity of a satellite communication constellation, wherein the satellite communication constellation includes a plurality of non-geosynchronous orbit satellites, the method comprising:
obtaining constraint conditions of system capacity performance in the satellite communication constellation, wherein the constraint conditions comprise constraint conditions of a mapping relation between beams and earth stations, constraint conditions of transmitting power and constraint conditions of interference values, and the interference values are characterized by the mapping relation and the transmitting power;
calculating according to the constraint condition to obtain an optimal solution of the transmitting power and an optimal solution of the mapping relation;
communicating according to the optimal solution of the transmitting power and the optimal solution of the mapping relation;
the method further comprises the following steps:
calculating to obtain the maximum throughput according to the optimal solution of the transmitting power and the optimal solution of the mapping relation;
calculating to obtain the maximum throughput according to the optimal solution of the transmitting power and the optimal solution of the mapping relation, and further comprising:
the maximum throughput is calculated according to the following formula:
Figure FDA0002489625780000011
Figure FDA0002489625780000012
Figure FDA0002489625780000013
Figure FDA0002489625780000014
Figure FDA0002489625780000015
Figure FDA0002489625780000016
the satellite communication constellation comprises S non-geosynchronous orbit satellites and N earth stations, wherein each non-geosynchronous orbit satellite is provided with K electric scanning phased array wave beams, ai,jRepresenting the mapping of the ith beam to the jth earth station, pi,jRepresenting the magnitude of the transmitted power of the ith beam to the jth earth station, ci,jRepresenting the maximum communication rate per link, qi,jIs the product coefficient, IthIn order to achieve an interference threshold that is not exceeded by lumped interference for co-frequency geosynchronous orbit communication systems,
Figure FDA0002489625780000021
is the interference value.
2. The method according to claim 1, wherein the constraint conditions include a predetermined mapping relationship that the beam and the earth station need to satisfy, an interference value that needs to be smaller than an interference threshold value, a transmission power that needs to be within a power threshold range, and a variation range of the space-time characteristics of each satellite in the satellite communication constellation corresponding to a predetermined orbital configuration.
3. The method of claim 2, wherein calculating the optimal solution of the transmit power and the optimal solution of the mapping relationship according to the constraint condition comprises:
carrying out solving calculation according to the interference characteristic and the constraint characteristic to obtain the transmitting power and the mapping relation on each link;
and carrying out iterative calculation on the transmitting power and the mapping relation on each link to respectively obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation.
4. The method of claim 3, wherein the obtaining the transmission power and the mapping relationship on each link by performing solution calculation according to the interference characteristic and the constraint characteristic comprises:
and calculating the interference value and the interference threshold value by accumulating and representing the interference value by the product of the transmitting power and the mapping relation of the plurality of non-geosynchronous earth orbit satellites to obtain the transmitting power and the mapping relation on each link which meet the requirement of communication service.
5. The system capacity optimization method of claim 1, further comprising:
and adjusting the modulation mode of the satellite communication constellation.
6. The method according to claim 5, wherein the modulation scheme comprises one or more of an adaptive modulation and coding scheme, a pulse code modulation scheme, a phase keying modulation scheme, a quadrature amplitude modulation scheme, and a minimum frequency shift keying modulation scheme.
7. The method of claim 1, wherein before calculating the optimal solution of the transmit power and the optimal solution of the mapping relationship according to the constraint condition, the method further comprises:
and (4) performing convex optimization problem transformation operation to enable the transformed convex optimization problem to meet the standard form of the convex optimization problem.
8. An apparatus for system capacity optimization of a satellite communication constellation, the satellite communication constellation including a plurality of non-geosynchronous orbit satellites, the apparatus comprising:
an obtaining module, configured to obtain constraint conditions of system capacity performance in the satellite communication constellation, where the constraint conditions include a constraint condition of a mapping relationship between a beam and an earth station, a constraint condition of transmit power, and a constraint condition of an interference value, and the interference value is characterized by the mapping relationship and the transmit power;
the calculation module is used for calculating to obtain the optimal solution of the transmitting power and the optimal solution of the mapping relation according to the constraint condition;
the communication module is used for carrying out communication according to the optimal solution of the transmitting power and the optimal solution of the mapping relation;
the calculation module is further configured to calculate a maximum throughput according to the optimal solution of the transmission power and the optimal solution of the mapping relationship;
the calculation module is further configured to calculate a maximum throughput according to the following formula:
Figure FDA0002489625780000041
Figure FDA0002489625780000042
Figure FDA0002489625780000043
Figure FDA0002489625780000044
Figure FDA0002489625780000045
Figure FDA0002489625780000046
the satellite communication constellation comprises S non-geosynchronous orbit satellites and N earth stations, wherein each non-geosynchronous orbit satellite is provided with K electric scanning phased array wave beams, ai,jRepresenting the mapping of the ith beam to the jth earth station, pi,jRepresenting the magnitude of the transmitted power of the ith beam to the jth earth station, ci,jRepresenting the maximum communication rate per link, qi,jIs the product coefficient, IthIn order to achieve an interference threshold that is not exceeded by lumped interference for co-frequency geosynchronous orbit communication systems,
Figure FDA0002489625780000047
is the interference value.
CN201910727794.4A 2019-08-07 2019-08-07 System capacity optimization method and device for satellite communication constellation Active CN110278024B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910727794.4A CN110278024B (en) 2019-08-07 2019-08-07 System capacity optimization method and device for satellite communication constellation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910727794.4A CN110278024B (en) 2019-08-07 2019-08-07 System capacity optimization method and device for satellite communication constellation

Publications (2)

Publication Number Publication Date
CN110278024A CN110278024A (en) 2019-09-24
CN110278024B true CN110278024B (en) 2020-07-24

Family

ID=67965868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910727794.4A Active CN110278024B (en) 2019-08-07 2019-08-07 System capacity optimization method and device for satellite communication constellation

Country Status (1)

Country Link
CN (1) CN110278024B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111865398B (en) * 2020-07-01 2022-06-03 哈尔滨工业大学（深圳） Satellite-ground transmission method under large-scale LEO satellite deployment
CN112994778B (en) * 2021-02-07 2022-07-29 哈尔滨工业大学 High-throughput satellite beam adaptive scheduling method based on service priority
CN113536550B (en) * 2021-06-29 2024-01-12 中国人民解放军63920部队 Situation map generation method and device, nonvolatile storage medium and processor
CN113595616B (en) * 2021-07-27 2022-05-17 北京邮电大学 NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization
CN113825238B (en) * 2021-09-03 2024-02-20 天地信息网络研究院(安徽)有限公司 Flexible signaling beam multi-scene application method in satellite communication
CN113691332B (en) * 2021-09-10 2022-05-27 北京科技大学 Co-channel interference characterization method and device of low-earth-orbit satellite communication system
CN114301512B (en) * 2021-12-10 2022-10-14 西安电子科技大学 Method and system for evaluating bearing capacity of multilayer low-orbit satellite network
CN114553300B (en) * 2022-02-18 2023-02-17 清华大学 Method and system for rapidly allocating wireless resources of multi-antenna earth station of satellite system
CN116865843B (en) * 2023-09-05 2023-12-26 武汉能钠智能装备技术股份有限公司 System and method for distributing links between satellites of satellite system
CN117544255A (en) * 2023-11-09 2024-02-09 中国人民解放军军事科学院系统工程研究院 NGSO satellite interference condition matrix manifold visualization method and device
CN117375706B (en) * 2023-12-04 2024-03-12 成都本原星通科技有限公司 Low-orbit inter-satellite interference optimization method and system for receiving end

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102355294A (en) * 2011-11-01 2012-02-15 东南大学 Multipoint coordinated beam forming and power allocation method for single base station power constraint
CN105744531A (en) * 2016-02-04 2016-07-06 中国空间技术研究院 Inline interference suppression based geostationary orbit (GEO) and non-geostationary orbit (NGEO) communication satellite spectrum sharing method
CN108419286A (en) * 2018-01-18 2018-08-17 北京邮电大学 A kind of allocation algorithm in face of 5G UAV Communications united beam and power
US10104641B2 (en) * 2013-03-01 2018-10-16 Intel IP Corporation Network selection in a heterogeneous network

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7286855B2 (en) * 1995-02-22 2007-10-23 The Board Of Trustees Of The Leland Stanford Jr. University Method and apparatus for adaptive transmission beam forming in a wireless communication system
US5822680A (en) * 1996-11-07 1998-10-13 Teledesic Llc Frequency sharing for satellite communication system
FR2782214B1 (en) * 1998-08-06 2006-11-10 Alsthom Cge Alkatel METHOD FOR ALLOCATING LINKS BETWEEN A SET OF ZONES AND A SET OF SATELLITES
CN107342804B (en) * 2016-04-28 2019-07-19 北京大学 Inhibit the method for adjacent star interference in satellite communication system
CN106027138B (en) * 2016-05-05 2018-07-13 清华大学 Evade the earth station system collinearly interfered with synchronous satellite and method
CN107276660B (en) * 2017-06-22 2018-08-31 清华大学 Resource allocation methods and device in non-orthogonal multiple air-ground coordination communication system
CN108882245B (en) * 2018-07-03 2023-02-28 中国人民解放军陆军工程大学 GEO and LEO cognitive satellite network and dynamic frequency allocation method thereof
CN108631840A (en) * 2018-07-05 2018-10-09 东南大学 Multi-beam satellite system robust pre-coding method based on outage probability constraint
CN109412682B (en) * 2018-12-19 2021-03-19 北京卫星信息工程研究所 Method for distributing frequency resources of medium-low orbit constellation satellite beams
CN109743098B (en) * 2018-12-30 2020-06-16 清华大学 Spectrum sharing method, related device, communication method and computer readable medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102355294A (en) * 2011-11-01 2012-02-15 东南大学 Multipoint coordinated beam forming and power allocation method for single base station power constraint
US10104641B2 (en) * 2013-03-01 2018-10-16 Intel IP Corporation Network selection in a heterogeneous network
CN105744531A (en) * 2016-02-04 2016-07-06 中国空间技术研究院 Inline interference suppression based geostationary orbit (GEO) and non-geostationary orbit (NGEO) communication satellite spectrum sharing method
CN108419286A (en) * 2018-01-18 2018-08-17 北京邮电大学 A kind of allocation algorithm in face of 5G UAV Communications united beam and power

Also Published As

Publication number Publication date
CN110278024A (en) 2019-09-24

Similar Documents

Publication Publication Date Title
CN110278024B (en) System capacity optimization method and device for satellite communication constellation
US9974081B2 (en) Method of power allocation and base station using the same
Kyrgiazos et al. A terabit/second satellite system for European broadband access: a feasibility study
Christopoulos et al. Linear and nonlinear techniques for multibeam joint processing in satellite communications
CN113595616B (en) NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization
Ge et al. Joint user pairing and power allocation for NOMA-based GEO and LEO satellite network
CA3163380C (en) Systems and methods for operating a multi-band satellite terminal
Wang et al. Spectrum optimization for cognitive satellite communications with cournot game model
Wang et al. Admission control and power allocation for NOMA-based satellite multi-beam network
CN116156421A (en) Differentiated service transmission method based on double-layer satellite heterogeneous network
Wang et al. On fairness optimization for NOMA-enabled multi-beam satellite systems
CN115882924A (en) NOMA-based beam hopping satellite communication system design method
CN113365288B (en) NB-IoT system uplink resource allocation method based on SWIPT
Lee et al. Interference Mitigation for Reverse Spectrum Sharing in B5G/6G Satellite-Terrestrial Networks
Lagunas et al. Resource allocation for cognitive satellite uplink and fixed-service terrestrial coexistence in Ka-band
Yu et al. Adaptive spatial modulation and thresholds optimization for MIMO systems in correlated Rayleigh channels
Adeane et al. Centralised and distributed power allocation algorithms in cooperative networks
Bhowal et al. Performance evaluation of RIS-assisted full-duplex MIMO bidirectional communications with a realistic channel model
KR102263034B1 (en) Method for transmitting down-link signal of wireless communication system
CN114337902A (en) IRS (inter-cell interference) assisted millimeter wave multi-cell interference suppression method
Stallo et al. Link performance analysis of multi-user detection techniques for W-band multi-beam satellites
Jo et al. Analysis on the neural network-aided satellite resource allocation schemes
KR20090108358A (en) Apparatus and method for interference cancellation using co-channel interferrence in broadband wireless access system
CN113824490B (en) Soft switching method based on satellite-to-ground link uplink non-orthogonal multiple access
Li et al. Max-Min fair resource allocation for min-rate guaranteed services in distributed antenna systems

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