CN113014306A - NGEO and GEO satellite spectrum sharing method based on interference control - Google Patents

NGEO and GEO satellite spectrum sharing method based on interference control Download PDF

Info

Publication number
CN113014306A
CN113014306A CN202110201112.3A CN202110201112A CN113014306A CN 113014306 A CN113014306 A CN 113014306A CN 202110201112 A CN202110201112 A CN 202110201112A CN 113014306 A CN113014306 A CN 113014306A
Authority
CN
China
Prior art keywords
ngeo
geo
satellite
user
calculating
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.)
Pending
Application number
CN202110201112.3A
Other languages
Chinese (zh)
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110201112.3A priority Critical patent/CN113014306A/en
Publication of CN113014306A publication Critical patent/CN113014306A/en
Pending legal-status Critical Current

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
    • 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/18521—Systems of inter linked satellites, i.e. inter satellite service
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/14—Spectrum sharing arrangements between different networks
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/22—Traffic simulation tools or models
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/24—Cell structures
    • H04W16/28—Cell structures using beam steering

Landscapes

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

Abstract

The invention discloses an NGEO and GEO satellite spectrum sharing method based on interference control, which relates to the field of wireless communication, and aims to improve the spectrum utilization rate on the premise of ensuring the QoS of GEO, wherein under the double-satellite cognitive radio architecture, the specific method comprises the following steps: and calculating the distance between two satellites and users and the off-axis angle between a user receiving antenna and a satellite transmitting antenna according to the exchanged ephemeris, then planning NGEO beams and calculating the total transmission rate of all users in a beam coverage area, listing an optimization problem under the condition that the NGEO satellite has the maximum throughput while protecting the QoS of GEO users, solving the optimization problem, obtaining the transmitting power P of each NGEO satellite beam, and completing the power distribution of one-time spectrum sharing.

Description

NGEO and GEO satellite spectrum sharing method based on interference control
Technical Field
The invention relates to the field of wireless communication, in particular to a multi-beam power control technology for spectrum sharing of a non-geostationary satellite system (NGEO) and a geostationary orbit satellite system (GEO).
Background
With the increasing demand for global communication, the construction of spatial information networks is increased in various countries, and satellites as important components of the spatial information networks have the advantages that other communication forms do not have: the method has the advantages of long communication distance, no limitation of any complex geographic conditions between two communication points, no influence of natural disasters and artificial events between the two communication points, and the like, and in recent years, more and more satellites with huge numbers are transmitted as follows: SpaceX, One Web, LeoSat and the like, which are mostly composed of dozens to hundreds of satellites, provide challenges for satellite spectrum resource management, and the shortage of spectrum resources will become a bottleneck limiting the development of satellite communication, so that how to increase spectrum efficiency and improve spectrum utilization rate become problems to be solved urgently.
Next generation satellite communication requires a satellite to have higher spectrum efficiency to solve the current problem of lower spectrum utilization, and in order to achieve this goal, different satellite systems need to be able to reasonably coexist in the same frequency band by some method to increase spectrum efficiency. In order to solve the problem of low space spectrum utilization rate, the NGEO system and the GEO system can share the same spectrum, with the continuous increase of service requirements, the traditional low-orbit communication frequency band cannot meet the development trend of NGEO satellites, the Ka frequency band has rich frequency band resources, the Ka frequency band is the development trend of the NGEO satellites, however, when the NGEO works in the Ka frequency band, co-channel interference is generated between the NGEO and the GEO satellites which are deployed in the frequency band, how to coexist the NGEO and the GEO, further and further research is needed, how to reasonably share the spectrum, and reducing the co-channel interference brought by spectrum sharing is the key point of the current research. In order to ensure good satellite transmission performance, corresponding measures are required to reduce co-channel interference.
The most typical technology adopted for spectrum sharing is Cognitive Radio (CR), a father Mitola of software Radio in 1999 proposes a concept of Cognitive Radio, and provides a solution for alleviating the contradiction between spectrum resource shortage and low actual utilization rate, a Secondary User (SU) can share the same frequency band with a Primary User (PU), a Cognitive satellite network can be generally divided into a satellite-ground hybrid Cognitive network and a dual-satellite Cognitive network, and the application of the CR technology in the satellite network faces the problems of long path, large loss, long delay and the like.
Compared with a single-beam satellite, the multi-beam satellite can obtain optimal coverage by adjusting the number of beams, the networking scheme of spot beams and the transmission power of each spot beam, and the multi-beam satellite generally adopts a multiplexing technology to improve the spectrum utilization rate. The power control method is introduced into a dual-satellite cognitive network of NGEO and GEO, co-channel interference generated by spectrum sharing between NGEO and GEO can be reduced by reducing the transmitting power of an interference beam of an NGEO satellite system, and the Quality of Service (QoS) of the GEO system and a ground user is ensured.
Disclosure of Invention
The invention provides an NGEO and GEO satellite spectrum sharing method based on interference control, aiming at improving the spectrum utilization rate on the premise of ensuring the QoS of GEO.
A double-satellite cognitive radio architecture is established and comprises a GEO satellite and a plurality of NGEO satellites, wherein the GEO satellite and the NGEO satellites are multi-beam satellites, gateway stations of the two satellites are connected through high-speed lossless optical fibers, ephemeris, a frequency configuration scheme and an antenna directional diagram of each satellite are exchanged, a GEO user serves as a main user, the NGEO user serves as a secondary user, and the NGEO system reuses authorized frequency of the GEO system.
NGEO and GEO frequency spectrum sharing multi-beam power control algorithm based on frequency planning is characterized by comprising the following steps:
step A: and calculating the distance between the two satellites and the user according to the exchanged ephemeris, and the off-axis angle between the receiving antenna of the user and the transmitting antenna of the satellite. The method specifically comprises the following steps:
step A1: calculating the position coordinates (x) of the NGEO satellite in the Earth-Centered-Fixed (ECEF) coordinate system according to the shared ephemeris of the two satellitessl,ysl,zsl) GEO satellite (x)sg,ysg,zsg) And NGEO terrestrial users (x)ul,yul,zul) GEO ground user (x)ug,yug,zug)。
Step A11: six satellite orbits can be obtained according to the ephemeris of the satellite: semi-major axis a of the track, eccentricity e, inclination angle i of the track, ascension omega of the intersection point0Angular distance of near point ω, true angle of near point M0And calculating the coordinate mode of the satellite in the ECEF according to the six orbital elements as follows:
by the formula
Figure BDA0002948953380000021
Calculating the average angular velocity n0Wherein mu is 398600.5 × 109m3/s2From the formula tk＝t-tpCalculating the time t from the ephemeris epoch StartkFrom the formula Mk＝M0+n0tkCalculating mean and near point angle MkFrom the formula Ek＝Mk+esinEkTo solve the off-angle EkFrom the formula
Figure BDA0002948953380000022
And
Figure BDA0002948953380000023
solving the true near point angle fkFrom the formula uk＝fk+ omega calculation of the angle of ascent ukFrom the formula Ωk＝Ω0-Ωe(tk-tp) Calculating the corrected intersection point Huangjing omegakFrom the formula rk＝a(1-ecosEk) Calculating the radius r of the trackkFrom xk＝rk cosukAnd yk＝rksinukCalculating the x-axis and y-axis coordinates on the orbital plane, so that the satellite is inThe coordinates in ECEF are:
Figure BDA0002948953380000031
step A12: by latitude and longitude coordinates of the location of the user
Figure BDA0002948953380000032
And the altitude h of the satellite0And calculating the coordinates of the user in the ECEF according to the following calculation formula:
Figure BDA0002948953380000033
step A2: calculating a vector V of the NGEO satellite pointing to the NGEO user according to the position coordinates obtained in the step A1ngeo_ngeouserVector V of NGEO satellite pointing to GEO usersngeo_geouserVector V of GEO satellite pointing to NGEO usergeo_ngeouserVector V of GEO satellite pointing to GEO usergeo_geouserAnd beam pointing R of NGEO satelliteblBeam pointing R of GEObg。
Step A3: calculating the distance d from the NGEO satellite to the NGEO user from the coordinates and vectors in the step A2 and the step A1ngeo_ngeouserDistance d from NGEO satellite to GEO userngeo_geouserDistance d from GEO satellite to NGEO usergeo_ngeouserDistance d from GEO satellite to GEO usergeo_geouserOff-axis angle θ between NGEO transmit antennas to NGEO usersngeo_ngeouserOff-axis angle θ between NGEO users to GEO usersngeo_geouserOff-axis angle θ between GEO transmit antenna to NGEO usergeo_ngeouserOff-axis angle θ between GEO users to GEO usersgeo_geouserOff-axis angle θ for NGEO users to receive GEO interference signalsngeouser_geoOff-axis angle theta of GEO user receiving NGEO satellite interference signalgeouser_ngeo。
And B: the method comprises the following steps of reasonably planning the NGEO beam:
step B1: when a satellite spectrum sharing architecture is established, the NGEO satellite system and the GEO satellite system share respective frequency configuration schemes through optical fiber connection of a gateway station, and the NGEO avoids changing the self spot beam frequency allocation based on interference after the frequency configuration scheme of the GEO satellite beam is known.
Step B2: when the GEO satellite and the NGEO satellite adopt seven-color multiplexing, beam frequency allocation is carried out on the principle that the minimum times of superposition among co-frequency beams is generated in order to reduce the number of co-frequency interference spot beams. When the seven-color multiplexing of the GEO wave beam frequency is performed from left to right, the following steps are performed from top to bottom in sequence: f. of3,f4,f5,f1,f2,f6,f7When planning the frequency distribution of NGEO as f7,f6,f2,f1,f5,f4,f3。
And C: calculating the signal to interference and noise ratio (SINR) of the NGEO user and the GEO user in the sharing process of the NGEO spectrum and the GEO spectrum according to the related angle and distance obtained in the step A
Step C1: calculating SINR of user u in NGEO beam i according to the angle of the distance between the NGEO satellite and the user obtained in the step Angeo,i,uThe calculation formula is as follows:
Figure BDA0002948953380000041
wherein N isgeoNumber of beams, P, for NGEO satellitesngeo,iTransmit power, P, for NGEO satellite to beam igeo,jTransmission power, h, of beam j in a GEO satelliteDngeo,i,uAnd hIgeo,i,uRepresenting respectively the channel parameters, N, of the signals expected to be received by user u in NGEO beam i and the signals from the GEO satellite expected to be received by terrestrial user u in beam j within the coverage of the GEOLThe power of white gaussian noise received for an NGEO terrestrial user. F is the interference coefficient between the sending end and the receiving end of the interference signal, and the value of the F is superposed with the beam planning result in the step B, the visibility and the bandwidth of the sending end and the receiving end of the interference signalThe value of Γ is in a closed interval from 0 to 1, is 0 if two nodes are not visible, and is 1 if two nodes are visible and the bandwidths are completely overlapped. The channel parameters are calculated as follows:
Figure BDA0002948953380000042
Figure BDA0002948953380000043
where c is the speed of light, f is the communication frequency at that time, GT，GRAntenna gains for the transmit antenna and the receive antenna, respectively. The power of gaussian white noise received by an NGEO user is:
NL＝KTnBl (6)
where K is the Boltzmann constant, TnFor the noise temperature at which the ground user is located, BlIs the bandwidth of NGEO.
Step C2: calculating the SINR of the user u in the beam j in the GEO beam range according to the angle of the distance between the GEO satellite and the user obtained in the step Ageo,j,uThe calculation formula is as follows:
Figure BDA0002948953380000051
wherein N isleoNumber of beams for NGEO satellite, NGPower of white Gaussian noise received for GEO terrestrial users, where NGThe calculation formula of (a) is as follows:
NG＝KTnBg (8)
wherein B isgIs the bandwidth of the GEO.
Step D: finding the user with the worst communication quality from the SINRs of all GEO users in the GEO beam range obtained in the step C, and if the SINR of the user is lower than the protection threshold gamma of the GEOthStep E is performed.
Step E: and C, calculating the total transmission rate of all the users in the NGEO coverage area according to the SINR of the NGEO and GEO users obtained in the step C
Step E1: calculating the transmission rate of the single user in one beam according to the SINR of the single user located in the NGEO beam i obtained in step C1, wherein the calculation formula is as follows:
Figure BDA0002948953380000052
step E2: from R obtained in step E1ngeo,i,uAnd calculating the transmission rate of all users in the NGEO beam i:
Figure BDA0002948953380000053
step E3: from R obtained in step E2ngeo,iCalculating the total transmission rate R of all users under the coverage of NGEO full beamngeo：
Figure BDA0002948953380000054
And F, listing an optimization problem based on the thought that the NGEO system has the maximum throughput while the QoS of the GEO user is protected, and rewriting the optimization problem into a convex optimization problem which can be solved.
Step F1: the SINR of the NGEO and GEO users obtained in the step C and the transmission rate R of the NGEO and GEO obtained in the step EngeoAnd the power, bandwidth and transmission rate limits of the satellite, and the problem is converted into the following convex optimization problem:
Figure BDA0002948953380000061
Figure BDA0002948953380000062
Figure BDA0002948953380000063
Figure BDA0002948953380000064
Figure BDA0002948953380000065
0≤Pngeo,i≤Pngeo,i,max (17)
Pgeo≥0 (18)
the above optimization problem is explained as follows: equation (12) maximizes the communication capacity of the NGEO satellite system during spectrum sharing for the optimization problem objective function, and equation (13) is a limiting condition, ensuring that any authorized non-geostationary orbit satellite (NGEO) specified by ITU must not cause unacceptable interference to GEO satellites, and ensuring that any user of GEO satellites meets the requirements of quality of service, where γthFor the preset protection threshold, the formula (14) is a limiting condition, and the maximum communication capacity under each NGEO beam is limited not to exceed Cbeam_maxThe limit of the transmitting power of the satellite beam is given by the equations (15) to (18), the on-board power is limited, and the total power of the satellite beam cannot exceed the limited maximum power Pngeo_totalAnd Pgeo_totalAnd the power of the beam is positive and cannot be less than 0.
Step F2: in step F1, the optimization problem is a Fractional Programming (FP) problem with a multi-ratio, which is difficult to directly solve, so that the original objective function is rewritten into the following form by using quadratic transformation:
Figure BDA0002948953380000066
in the formula (ii)i,uThe update method of (1) is as follows:
Figure BDA0002948953380000067
the limitation of the optimization problem in step F2 is formula (13) to formula (18) in F1.
Step G: solving the optimization problem obtained in the step F to obtain the transmitting power P of each NGEO satellite beam ii。
Step G1: initialization vector PngeoIs a suitable value.
Step G2: the vector y is updated according to equation (20).
Step G3: p is obtained by solving the convex optimization problem after y is fixed in step F2ngeoAnd update PngeoAnd (5) vector quantity.
Step G4: repeating steps G2 and G3 until the objective function (19)
Figure BDA0002948953380000071
And (6) converging. And obtaining the proper transmitting power of each NGEO beam after power control.
The method has the advantages that a large number of simulation experiments show that the function rule control algorithm is fast in convergence, the QoS of the GEO user is guaranteed, meanwhile, the communication capacity of the NGEO satellite is improved to the maximum extent, the NGEO satellite and the GEO satellite can coexist to a certain extent, and the spatial spectrum efficiency is remarkably improved.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a schematic diagram of NGEO satellite and GEO satellite spectrum sharing of the present invention;
FIG. 3 is a schematic diagram of the beam frequency planning of the present invention;
FIG. 4 is a geometric diagram of an interference model used in co-channel interference analysis according to the present invention;
fig. 5 is a simulation diagram of signal quality for LEO satellite spectrum sharing different users without beam frequency planning;
FIG. 6 is a graph of signal quality simulations for different users sharing the GEO satellite spectrum without beam frequency planning;
fig. 7 is a simulation diagram of signal quality for LEO satellite spectrum sharing different users with beam frequency planning;
FIG. 8 is a graph of signal quality simulations of GEO satellite spectrum sharing different users with beam frequency planning;
FIG. 9 is a schematic diagram of LEO satellite beam power versus iteration number for beamless frequency planning;
FIG. 10 is a graphical illustration of LEO satellite beam power versus iteration number for beam frequency planning;
fig. 11 is a schematic diagram of a simulation of the total velocity of a LEO satellite at different interference thresholds after power control without beam frequency planning;
fig. 12 is a schematic diagram of a simulation of the total velocity of a LEO satellite at different interference thresholds after power control during beam frequency planning;
Detailed Description
First embodiment, the first embodiment will be described with reference to fig. 1 to 11
A schematic diagram of the NGEO satellite and GEO satellite spectrum sharing adopted in the present embodiment is shown in fig. 2.
The spectrum sharing network architecture shown in fig. 2 is composed of GEO satellites and NGEO satellites, the NGEO satellite system adopts a Walker constellation architecture, the Walker constellation is a constellation configuration with the characteristics of the same orbit height, uniformly distributed inclined circular orbit planes and the like, and the orbit height of the GEO satellites is 35786 Km. The GEO is used as a main user, the NGEO is used as a secondary user, the GEO satellite and the NGEO satellite share the same frequency band to provide service for users in a service area, and the NGEO system is used as the secondary user to communicate on the basis of not influencing the communication quality of the main user. The gateway stations of the two satellite systems exchange respective ephemeris, frequency allocation schemes and antenna patterns through high-speed lossless optical fiber connections. Since the gateways of the NGEO and GEO satellites are connected through the optical fiber and share the channel state information and the ephemeris information with each other, one party can predict when and where the interference occurs. All satellites use multi-beam antennas, frequency multiplexing is performed between beams, and it is assumed that their corresponding user antennas are always tracking the satellites, i.e., the elevation angles of the user antennas are dynamically changed.
Fig. 4 shows a geometric model of interference analysis used in spectrum-shared downlink co-channel interference analysis, in which a solid line represents a signal path to be received and a dotted line represents an interference signal path, in which angle and distance parameters to be calculated for calculating SINR are indicated.
NGEO and GEO frequency spectrum sharing multi-beam power control algorithm based on frequency planning is characterized by comprising the following steps:
step A: and calculating the distance between the two satellites and the user according to the exchanged ephemeris, and the off-axis angle between the receiving antenna of the user and the transmitting antenna of the satellite. The method specifically comprises the following steps:
step A1: calculating the position coordinates (x) of the NGEO satellite in the Earth-Centered-Fixed (ECEF) coordinate system according to the shared ephemeris of the two satellitessl,ysl,zsl) GEO satellite (x)sg,ysg,zsg) And NGEO terrestrial users (x)ul,yul,zul) GEO ground user (x)ug,yug,zug)。
Step A11: six satellite orbits can be obtained according to the ephemeris of the satellite: semi-major axis a of the track, eccentricity e, inclination angle i of the track, ascension omega of the intersection point0Angular distance of near point ω, true angle of near point M0And calculating the coordinate mode of the satellite in the ECEF according to the six orbital elements as follows:
by the formula
Figure BDA0002948953380000081
Calculating the average angular velocity n0Wherein mu is 398600.5 × 109m3/s2From the formula tk＝t-tpCalculating the time t from the ephemeris epoch StartkFrom the formula Mk＝M0+n0tkCalculating mean and near point angle MkFrom the formula Ek＝Mk+esinEkTo solve the off-angle EkFrom the formula
Figure BDA0002948953380000082
And
Figure BDA0002948953380000083
solving the true near point angle fkFrom the formula uk＝fk+ omega calculation of the angle of ascent ukFrom the formula Ωk＝Ω0-Ωe(tk-tp) Calculating the corrected intersection point Huangjing omegakFrom the formula rk＝a(1-ecosEk) Calculating the radius r of the trackkFrom xk＝rk cosukAnd yk＝rksinukAnd calculating x-axis coordinates and y-axis coordinates on the orbital plane, wherein the coordinates of the satellite in the ECEF are as follows:
Figure BDA0002948953380000091
step A12: by latitude and longitude coordinates of the location of the user
Figure BDA0002948953380000092
And the altitude h of the satellite0And calculating the coordinates of the user in the ECEF according to the following calculation formula:
Figure BDA0002948953380000093
step A2: calculating a vector V of the NGEO satellite pointing to the NGEO user according to the position coordinates obtained in the step A1ngeo_ngeouserVector V of NGEO satellite pointing to GEO usersngeo_geouserVector V of GEO satellite pointing to NGEO usergeo_ngeouserVector V of GEO satellite pointing to GEO usergeo_geouserAnd beam pointing R of NGEO satelliteblBeam pointing R of GEObg。
Step A3: calculating the distance d from the NGEO satellite to the NGEO user from the coordinates and vectors in the step A2 and the step A1ngeo_ngeouserDistance d from NGEO satellite to GEO userngeo_geouserDistance d from GEO satellite to NGEO usergeo_ngeouserDistance d from GEO satellite to GEO usergeo_geouserOff-axis angle θ between NGEO transmit antennas to NGEO usersngeo_ngeouserOff-axis angle θ between NGEO users to GEO usersngeo_geouserOff-axis angle θ between GEO transmit antenna to NGEO usergeo_ngeouserOff-axis angle θ between GEO users to GEO usersgeo_geouserOff-axis angle θ for NGEO users to receive GEO interference signalsngeouser_geoOff-axis angle theta of GEO user receiving NGEO satellite interference signalgeouser_ngeo。
And B: the method comprises the following steps of reasonably planning the NGEO beam:
step B1: when a satellite spectrum sharing architecture is established, the NGEO satellite system and the GEO satellite system share respective frequency configuration schemes through optical fiber connection of a gateway station, and the NGEO avoids changing the self spot beam frequency allocation based on interference after the frequency configuration scheme of the GEO satellite beam is known.
Step B2: when the GEO satellite and the NGEO satellite adopt seven-color multiplexing, beam frequency allocation is carried out on the principle that the minimum times of superposition among co-frequency beams is generated in order to reduce the number of co-frequency interference spot beams. When the seven-color multiplexing of the GEO wave beam frequency is performed from left to right, the following steps are performed from top to bottom in sequence: f. of3,f4,f5,f1,f2,f6,f7When planning the frequency distribution of NGEO as f7,f6,f2,f1,f5,f4,f3。
As illustrated in figure 3 of the specification.
And C: calculating the signal to interference and noise ratio (SINR) of the NGEO user and the GEO user in the sharing process of the NGEO spectrum and the GEO spectrum according to the related angle and distance obtained in the step A
Step C1: calculating SINR of user u in NGEO beam i according to the angle of the distance between the NGEO satellite and the user obtained in the step Angeo,i,uThe calculation formula is as follows:
Figure BDA0002948953380000101
wherein: n is a radical ofgeoNumber of beams, P, for NGEO satellitesngeo,iTransmit power, P, for NGEO satellite to beam igeo,jTransmission power, h, of beam j in a GEO satelliteDngeo,i,uAnd hIgeo,i,uRepresenting respectively the channel parameters, N, of the signals expected to be received by user u in NGEO beam i and the signals from the GEO satellite expected to be received by terrestrial user u in beam j within the coverage of the GEOLThe power of white gaussian noise received for an NGEO terrestrial user. And F is an interference coefficient between the transmitting end and the receiving end of the interference signal, the value of the F is related to the beam planning result in the step B, the visibility of the transmitting end and the receiving end of the interference signal and the overlapping size of the bandwidth, the value range of the F is in a closed interval from 0 to 1, if the two nodes are invisible, the F is 0, and if the two nodes are visible and the bandwidths are completely overlapped, the value of the F is 1. The channel parameters are calculated as follows:
Figure BDA0002948953380000102
Figure BDA0002948953380000103
where c is the speed of light, f is the communication frequency at that time, GT，GRAntenna gains for the transmit antenna and the receive antenna, respectively. The power of gaussian white noise received by an NGEO user is:
NL＝KTnBl (6)
where K is the Boltzmann constant, TnFor the noise temperature at which the ground user is located, BlIs the bandwidth of NGEO.
Step C2: calculating the SINR of the user u in the beam j in the GEO beam range according to the angle of the distance between the GEO satellite and the user obtained in the step Ageo,j,uThe calculation formula is as follows:
Figure BDA0002948953380000111
wherein N isleoNumber of beams for NGEO satellite, NGPower of white Gaussian noise received for GEO terrestrial users, where NGThe calculation formula of (a) is as follows:
NG＝KTnBg (8)
wherein B isgIs the bandwidth of the GEO.
Step D: finding the user with the worst communication quality from the SINRs of all GEO users in the GEO beam range obtained in the step C, and if the SINR of the user is lower than the protection threshold gamma of the GEOthStep E is performed.
Step E: and C, calculating the total transmission rate of all the users in the NGEO coverage area according to the SINR of the NGEO and GEO users obtained in the step C
Step E1: calculating the transmission rate of the single user in one beam according to the SINR of the single user located in the NGEO beam i obtained in step C1, wherein the calculation formula is as follows:
Figure BDA0002948953380000112
step E2: from R obtained in step E1ngeo,i,uAnd calculating the transmission rate of all users in the NGEO beam i:
Figure BDA0002948953380000113
step E3: from R obtained in step E2ngeo,iCalculating the total transmission rate R of all users under the coverage of NGEO full beamngeo：
Figure BDA0002948953380000114
And F, listing an optimization problem based on the thought that the NGEO system has the maximum throughput while the QoS of the GEO user is protected, and rewriting the optimization problem into a convex optimization problem which can be solved.
Step F1: the SINR of the NGEO and GEO users obtained in the step C and the transmission rate R of the NGEO and GEO obtained in the step EngeoAnd the power, bandwidth and transmission rate limits of the satellite, and the problem is converted into the following convex optimization problem:
Figure BDA0002948953380000121
Figure BDA0002948953380000122
Figure BDA0002948953380000123
Figure BDA0002948953380000124
Figure BDA0002948953380000125
0≤Pngeo,i≤Pngeo,i,max (17)
Pgeo≥0 (18)
the above optimization problem is explained as follows: equation (12) maximizes the communication capacity of the NGEO satellite system during spectrum sharing for the optimization problem objective function, and equation (13) is a limiting condition, ensuring that any authorized non-geostationary orbit satellite (NGEO) specified by ITU must not cause unacceptable interference to GEO satellites, and ensuring that any user of GEO satellites meets the requirements of quality of service, where γthFor the preset protection threshold, the formula (14) is a limiting condition, and the maximum communication capacity under each NGEO beam is limited not to exceed Cbeam_maxFormula (15) -formula (18)) For limiting the transmitting power of the satellite beam, the on-board power is limited, and the total power of the beam of the satellite cannot exceed the limited maximum power Pngeo_totalAnd Pgeo_totalAnd the power of the beam is positive and cannot be less than 0.
Step F2: in step F1, the optimization problem is a Fractional Programming (FP) problem with a multi-ratio, which is difficult to directly solve, so that the original objective function is rewritten into the following form by using quadratic transformation:
Figure BDA0002948953380000126
in the formula (ii)i,uThe update method of (1) is as follows:
Figure BDA0002948953380000131
the limitation of the optimization problem in step F2 is formula (13) to formula (18) in F1.
Step G: solving the optimization problem obtained in the step F to obtain the transmitting power P of each NGEO satellite beam ii。
Step G1: initialization vector PngeoIs a suitable value.
Step G2: the vector y is updated according to equation (20).
Step G3: p is obtained by solving the convex optimization problem after y is fixed in step F2ngeoAnd update PngeoAnd (5) vector quantity.
Step G4: repeating steps G2 and G3 until the objective function (19)
Figure BDA0002948953380000132
And (6) converging. And obtaining the proper transmitting power of each NGEO beam after power control.
The invention is subjected to simulation verification, and the parameters are as follows:
TABLE 1 orbital parameters of LEO and GEO satellites
Figure BDA0002948953380000133
TABLE 2 simulation parameters
Figure BDA0002948953380000134
Figure BDA0002948953380000141
Simulation description: for the convenience of simulation explanation, it is assumed that there is one LEO user and one GEO user under each GEO beam, both of which are located at the center of the beam, the simulation time is a fraction of the period during which the LEO beam passes the top, and the simulation time interval is set to 1 s.
Firstly, simulating the SINRs received by the LEO and GEO users before and after beam planning, wherein the simulation results are shown in FIG. 5 and FIG. 6, figure 5 is a simulation diagram of the signal quality of users with LEO and GEO satellite spectrum sharing on different beams without beam frequency planning, fig. 6 is a simulation diagram of the signal quality of users with different beams for LEO and GEO satellite spectrum sharing in the presence of beam frequency planning, and it can be seen from fig. 5 and 6 that in the worst case without frequency planning, it is possible that all beams will generate co-channel interference, resulting in that the SINR of all users cannot satisfy the communication quality of the users, when beam frequency planning is performed, as shown in fig. 7 and 8, the beams of the LEO satellite and the GEO satellite do not all generate co-channel interference at the same time, under the parameters set by simulation, only one co-frequency beam is superposed in a period of time to generate co-frequency interference. The user suffering interference cannot guarantee the communication quality.
Fig. 9 and 10 show the convergence of the power of each beam of the LEO satellite with the whole algorithm when the interference is the most serious in the simulation process, and it can be seen from the simulation diagram that the convergence of the algorithm is good, and fig. 9 shows the simulation situation when the beam frequency planning is not performed, that all beams reduce the transmission power in order to ensure the service quality of GEO users because there are more overlapped co-frequency beams compared to the case of planning by the beam frequency in fig. 10, and after the beam frequency planning, only beam 1 shown in fig. 10 reduces the transmission power in order to ensure that the communication of GEO is not affected.
Fig. 11 and 12 are simulation diagrams of the total velocity of LEO satellites under different reference thresholds after power control, fig. 11 is a simulation diagram of a beam without frequency planning, fig. 12 is a simulation diagram of a beam after frequency planning, and it can be seen from comparison between the two diagrams that the total velocity of LEO satellites is 10 when beam planning is not performed8In bps, the total rate is 10 after frequency planning9In bps order, the effect is improved significantly, specifically, when the protection threshold of the GEO user is 13dB, the worst case rate of fig. 11 is 1.8 × 108bps, and the worst rate under this threshold in FIG. 12 is 1.08 × 109bps。
The simulation shows that the function rule control algorithm is fast in convergence, the QoS of the GEO user is guaranteed, meanwhile, the communication capacity of the LEO satellite is improved to the maximum extent, the LEO satellite and the GEO satellite can coexist to a certain extent, and the spatial spectrum efficiency is improved.
While the present invention has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present invention and that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims, it is to be understood that various dependent claims may be combined with features described herein in other combinations than those described in the original claims, and it is to be understood that features described in connection with the individual embodiments may be used in other described embodiments.

Claims (7)

1. Firstly, establishing a double-satellite cognitive radio architecture, wherein the double-satellite cognitive radio architecture comprises a GEO satellite and N NGEO satellites, and the GEO satellite and each NGEO satellite are multi-beam satellites, wherein the two satellites are as follows: the gateway stations of the GEO satellite and one of the NGEO satellites are connected through high-speed lossless optical fibers, respective ephemeris, frequency configuration schemes and antenna directional diagrams are exchanged, under the double-satellite cognitive radio architecture, a GEO user serves as a primary user, an NGEO user serves as a secondary user, and the NGEO satellite reuses the authorized frequency of the GEO satellite;
the method is characterized in that: under the dual-satellite cognitive radio architecture, the NGEO and GEO satellite spectrum sharing method based on interference control comprises the following steps:
step one, calculating the distance between two satellites and a user according to the exchanged ephemeris, and calculating the off-axis angle between a user receiving antenna and a satellite transmitting antenna;
step two, planning the NGEO wave beam;
thirdly, obtaining the SINR of the NGEO user and the GEO user in the sharing process of the NGEO and GEO frequency spectrums according to the off-axis angle and the distance obtained in the first step;
step four, the user with the worst communication quality is obtained from the signal to interference plus noise ratio SINR of all GEO users in the GEO satellite beam range obtained in the step three, and if the SINR of the user is lower than the protection threshold gamma of the GEOthExecuting the step five;
step five, calculating the total transmission rate of all users in the NGEO coverage area according to the SINR of the NGEO and GEO users obtained in the step C;
step six, under the condition that the NGEO satellite has the maximum throughput while the QoS of the GEO user is protected, listing optimization problems, and rewriting the listed optimization problems into a convex optimization problem capable of being solved;
step seven, solving the optimization problem listed in the step six to obtain the transmitting power P of each NGEO satellite beam iiI is a positive integer; and completing power allocation of NGEO and GEO satellite spectrum sharing based on interference control.
2. The method for NGEO and GEO satellite spectrum sharing based on interference control as claimed in claim 1, wherein in step one, the distance between the two satellites and the user is calculated according to the exchanged ephemeris, and the off-axis angle between the user receiving antenna and the satellite transmitting antenna is specifically:
step one, calculating the position coordinate (x) of the NGEO satellite in an Earth-Centered Earth-Fixed coordinate system (Earth-Centered, Earth-Fixed, ECEF) according to the ephemeris shared by the two satellitessl,ysl,zsl) GEO satellite (x)sg,ysg,zsg) And NGEO terrestrial users (x)ul,yul,zul) GEO ground user (x)ug,yug,zug) The method specifically comprises the following steps:
step A1, acquiring six numbers of satellite orbits according to ephemeris of the satellites: semi-major axis a of the track, eccentricity e, track inclination angle i, elevation intersection right ascension omega0Angular distance of near point omega and true angle of near point M0The method for calculating the coordinate of the satellite in the ECEF according to the six orbital parameters comprises the following steps:
according to the formula
Figure FDA0002948953370000021
Calculating the average angular velocity n0(ii) a Wherein mu is 398600.5 × 109m3/s2According to the formula tk＝t-tpCalculating the time t from the ephemeris epoch Startk(ii) a According to formula Mk＝M0+n0tkCalculating mean and near point angle Mk(ii) a According to formula Ek＝Mk+esinEkTo solve the off-angle EkAccording to the formula
Figure FDA0002948953370000022
And
Figure FDA0002948953370000023
solving the true near point angle fk(ii) a According to the formula uk＝fk+ omega calculation of the angle of ascent uk(ii) a According to the formula omegak＝Ω0-Ωe(tk-tp) Calculating the corrected intersection point Huangjing omegak(ii) a According to the formula rk＝a(1-ecosEk) Calculating the radius r of the trackk(ii) a According to the formula xk＝rkcosukAnd yk＝rksinukCalculating x-axis and y-axis coordinates on a track plane; the coordinates of the satellite in ECEF are then:
Figure FDA0002948953370000024
step A2: according to the longitude and latitude coordinates of the position of the user
Figure FDA0002948953370000025
And the altitude h of the satellite0And calculating the coordinates of the user in the ECEF according to the following calculation formula:
Figure FDA0002948953370000026
the first step is: calculating a vector V of the NGEO satellite pointing to the NGEO user according to the position coordinates obtained in the step one by onengeo_ngeouserVector V of NGEO satellite pointing to GEO usersngeo_geouserVector V of GEO satellite pointing to NGEO usergeo_ngeouserVector V of GEO satellite pointing to GEO usergeo_geouserAnd beam pointing R of NGEO satelliteblAnd beam pointing R of GEObg；
Step one is three: calculating the distance d from the NGEO satellite to the NGEO user according to the coordinates and the vectors in the step onengeo_ngeouserDistance d from NGEO satellite to GEO userngeo_geouserDistance d from GEO satellite to NGEO usergeo_ngeouserDistance d from GEO satellite to GEO usergeo_geouserOff-axis angle theta between NGEO transmit antenna to NGEO userngeo_ngeouserOff-axis angle θ between NGEO user to GEO userngeo_geouserOff-axis angle θ between GEO transmit antenna to NGEO usergeo_ngeouserOff-axis angle θ between GEO user to GEO usergeo_geouserOff-axis angle theta of NGEO user receiving GEO interference signalngeouser_geoAnd off-axis angle theta of GEO user receiving NGEO satellite interference signalgeouser_ngeo。
3. The interference control based NGEO and GEO satellite spectrum sharing method according to claim 2, wherein in step two, NGEO beam is planned; the method specifically comprises the following steps:
step two, according to the fact that the known NGEO satellite and the GEO satellite share respective frequency configuration schemes through optical fiber connection of a gateway station when a satellite spectrum sharing architecture is established, the NGEO satellite avoids changing the self spot beam frequency distribution based on interference after the frequency configuration scheme of the GEO satellite beam is known;
step two: when the GEO satellite and the NGEO satellite are multiplexed, in order to reduce the number of the co-channel interference spot beams, the beam frequency is allocated on the principle that the least times of superposition between the co-channel beams is generated, and when the GEO satellite beam frequency seven-color multiplexing is performed from left to right, the following steps are performed from top to bottom: f. of3,f4,f5,f1,f2,f6,f7When planning a frequency assignment of NGEO satellites as f7,f6,f2,f1,f5,f4,f3。
4. The method for sharing NGEO and GEO satellite spectrum based on interference control according to claim 3, characterized in that in step three, the SINR of NGEO users and GEO users in the process of sharing NGEO and GEO spectrum is obtained according to the off-axis angle and the distance obtained in step one; the specific method comprises the following steps:
step three, calculating the SINR of the user u in the NGEO beam i according to the angle of the distance between the NGEO satellite and the user obtained in the step onengeo,i,uThe calculation formula is as follows:
Figure FDA0002948953370000031
wherein N isgeoNumber of beams, P, for NGEO satellitesngeo,iTransmit power, P, for NGEO satellite to beam igeo,jTransmission power, h, of beam j in a GEO satelliteDngeo,i,uAnd hIgeo,i,uRepresenting respectively the channel parameters, N, of the signals expected to be received by user u in NGEO beam i and the signals from the GEO satellite expected to be received by terrestrial user u in beam j within the coverage of the GEOLThe power of white gaussian noise received for an NGEO terrestrial user. F is an interference coefficient between the sending end and the receiving end of the interference signal, the value of the F is related to the beam planning result in the step B, the visibility of the sending end and the receiving end of the interference signal and the overlapping size of the bandwidth, the value range of the F is in a closed interval from 0 to 1, if the two nodes are invisible, the F is 0, and if the two nodes are visible and the bandwidths are completely overlapped, the value of the F is 1;
the channel parameters are calculated as follows:
Figure FDA0002948953370000041
Figure FDA0002948953370000042
where c is the speed of light, f is the communication frequency at that time, GT，GRAntenna gains for the transmit antenna and the receive antenna, respectively.
The power of gaussian white noise received by an NGEO user is:
NL＝KTnBl (12)
where K is the Boltzmann constant, TnFor the noise temperature at which the ground user is located, BlIs the bandwidth of NGEO;
step three, calculating the SINR of the user u in the beam j in the GEO beam range according to the angle of the distance between the GEO satellite and the user obtained in the step onegeo,j,uThe calculation formula is as follows:
Figure FDA0002948953370000043
wherein N isleoNumber of beams for NGEO satellite, NGPower of white Gaussian noise received for GEO terrestrial users, where NGThe calculation formula of (2) is as follows:
NG＝KTnBg (14)
wherein B isgIs the bandwidth of the GEO satellite.
5. The interference control-based NGEO and GEO satellite spectrum sharing method according to claim 4, wherein in the fifth step, the specific method for calculating the total transmission rate of all users located in the NGEO coverage area according to the SINR of the NGEO and GEO users obtained in the third step is as follows:
step five, calculating the transmission rate of a single user in one beam according to the SINR of the single user in the NGEO beam i obtained in the step three, wherein the calculation formula is as follows:
Figure FDA0002948953370000044
step five two: r obtained in the fifth stepleo,i,uAnd calculating the transmission rate of all users in the NGEO beam i:
Figure FDA0002948953370000051
step five and step three: from R obtained in step five or twongeo,iCalculating the total transmission rate R of all users under the coverage of NGEO full beamngeo：
Figure FDA0002948953370000052
6. The interference control-based NGEO and GEO satellite spectrum sharing method according to claim 5, characterized in that in step six, under the condition of ensuring the NGEO satellite to have the maximum throughput while protecting the QoS of the GEO users, the specific method for listing the optimization problem and rewriting the listed optimization problem into the convex optimization problem capable of being solved is as follows:
step six: according to the SINR of the NGEO users and GEO users obtained in the third step and the transmission rate R of the NGEO satellites and GEO obtained in the fifth stepleoAnd the power, bandwidth and transmission rate limits of the satellite, and the problem is converted into the following convex optimization problem:
Figure FDA0002948953370000053
Figure FDA0002948953370000054
Figure FDA0002948953370000055
Figure FDA0002948953370000056
Figure FDA0002948953370000057
0≤Pngeo,i≤Pngeo,i,max (23)
Pgeo≥0 (24)
the above optimization problem is explained as follows: equation (18) maximizes the communications capacity of the NGEO satellite system during spectrum sharing for the optimization problem objective function, and equation (13) is a constraint that ensures that any authorized non-geostationary orbit satellite (NGEO) specified by the ITU must not cause unacceptable interference with GEO satellitesAny user of the GEO satellite meets the quality of service requirements, where γthFor the preset protection threshold, the formula (20) is a limiting condition, and the maximum communication capacity under each NGEO beam is limited not to exceed Cbeam_maxThe equations (21) to (24) are the limits of the satellite beam transmission power, the on-board power is limited, the total power of the satellite beam cannot exceed the limited maximum power Pngeo_totalAnd Pgeo_totalAnd the power of the wave beam is positive and can not be less than 0;
step six two, step six in the optimization problem is the Fractional Programming (FP) problem of the multi-ratio, this kind of problem is difficult to solve directly, so adopt the quadratic transformation to change the original objective function over into the following form:
Figure FDA0002948953370000061
in the formula (ii)i,uThe update method of (1) is as follows:
Figure FDA0002948953370000062
the constraints of the optimization problem in step sixty-two are formula (19) to formula 24) in F1.
7. The interference control-based NGEO and GEO satellite spectrum sharing method according to claim 6, wherein in the seventh step, the optimization problem listed in the sixth step is solved in the seventh step, and the transmission power P of each NGEO satellite beam i is obtainediThe specific method for completing power allocation of NGEO and GEO satellite spectrum sharing based on interference control for one time comprises the following steps:
step seven one, initializing vector Pngeo；
Step seven and two: according to equation 26) update vector y;
step seven and three: p is obtained by solving the convex optimization problem after y is fixed in the step six or twongeoAnd update PngeoAnd (5) vector quantity.
CN202110201112.3A 2021-02-23 2021-02-23 NGEO and GEO satellite spectrum sharing method based on interference control Pending CN113014306A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110201112.3A CN113014306A (en) 2021-02-23 2021-02-23 NGEO and GEO satellite spectrum sharing method based on interference control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110201112.3A CN113014306A (en) 2021-02-23 2021-02-23 NGEO and GEO satellite spectrum sharing method based on interference control

Publications (1)

Publication Number Publication Date
CN113014306A true CN113014306A (en) 2021-06-22

Family

ID=76407408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110201112.3A Pending CN113014306A (en) 2021-02-23 2021-02-23 NGEO and GEO satellite spectrum sharing method based on interference control

Country Status (1)

Country Link
CN (1) CN113014306A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114916051A (en) * 2022-05-24 2022-08-16 桂林电子科技大学 LEO satellite power control method based on BP neural network

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108882245A (en) * 2018-07-03 2018-11-23 中国人民解放军陆军工程大学 A kind of GEO and LEO cognition satellite network and its dynamic frequency allocation method
US20180343055A1 (en) * 2016-05-03 2018-11-29 Theia Group, Incorporated Low earth orbit satellite constellation system for communications with re-use of geostationary satellite spectrum
CN109743098A (en) * 2018-12-30 2019-05-10 清华大学 Frequency spectrum sharing method, device and the electronic equipment of spectra sharing system
WO2020058503A1 (en) * 2018-09-21 2020-03-26 Methera Global Communications Limited A method of, and apparatus for, improved satellite communications

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180343055A1 (en) * 2016-05-03 2018-11-29 Theia Group, Incorporated Low earth orbit satellite constellation system for communications with re-use of geostationary satellite spectrum
CN108882245A (en) * 2018-07-03 2018-11-23 中国人民解放军陆军工程大学 A kind of GEO and LEO cognition satellite network and its dynamic frequency allocation method
WO2020058503A1 (en) * 2018-09-21 2020-03-26 Methera Global Communications Limited A method of, and apparatus for, improved satellite communications
CN109743098A (en) * 2018-12-30 2019-05-10 清华大学 Frequency spectrum sharing method, device and the electronic equipment of spectra sharing system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
H. WANG ET AL.: ""Coexistence Downlink Interference Analysis Between LEO System and GEO System in Ka Band"", 《2018 IEEE/CIC INTERNATIONAL CONFERENCE ON COMMUNICATIONS IN CHINA (ICCC)》 *
张伟忠: ""多波束卫星***频率复用与干扰避免算法"", 《中国优秀硕士学位论文全文数据库》 *
王闯等: ""空间信息网络中面向双卫星的频谱共享方法"", 《航空学报》 *
胡晓月等: ""基于认知无线电的GEO与LEO卫星频谱共存"", 《中国空间科学技术》 *
魏文秋: ""GEO/LEO双层卫星通信网络同频干扰避免技术研究"", 《中国优秀硕士学位论文全文数据库》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114916051A (en) * 2022-05-24 2022-08-16 桂林电子科技大学 LEO satellite power control method based on BP neural network

Similar Documents

Publication Publication Date Title
CN108882245B (en) GEO and LEO cognitive satellite network and dynamic frequency allocation method thereof
Wang et al. Ultra-dense LEO satellite-based communication systems: A novel modeling technique
Werner et al. Analysis of system parameters for LEO/ICO-satellite communication networks
Vatalaro et al. Analysis of LEO, MEO, and GEO global mobile satellite systems in the presence of interference and fading
US5669062A (en) Methods of demand-based adaptive channel reuse for telecommunications systems
CN111865398B (en) Satellite-ground transmission method under large-scale LEO satellite deployment
CN110518965A (en) A kind of Non GEO satellite is to geosynchronous satellite uplink interference bypassing method
Wang et al. Coexistence downlink interference analysis between LEO system and GEO system in Ka band
CN113131989A (en) Design method for frequency spectrum sharing simulation time parameter of NGSO constellation system
CN112399429B (en) Communication scene modeling method and system for satellite communication system
CN112532298A (en) Wireless resource scheduling method and architecture system based on low-earth-orbit constellation satellite communication
Li et al. Optimal beam power control for co-existing multibeam GEO and LEO satellite system
Lagunas et al. Carrier allocation for hybrid satellite-terrestrial backhaul networks
Wang et al. Multi-layer LEO satellite constellation design for seamless global coverage
Wang et al. Spectrum optimization for cognitive satellite communications with cournot game model
Jia et al. An analytic approach for modeling uplink performance of mega constellations
CN116156421A (en) Differentiated service transmission method based on double-layer satellite heterogeneous network
Jia et al. Joint multi-beam power control for LEO and GEO spectrum-sharing networks
CN113014306A (en) NGEO and GEO satellite spectrum sharing method based on interference control
US11178550B2 (en) Coordination of spectrum allocation and interference avoidance among high-altitude networks
Chen et al. Coverage analysis of sagin with sectorized beam pattern under shadowed-rician fading channels
Lee et al. Interference mitigation for reverse spectrum sharing in B5G/6G satellite-terrestrial networks
Li et al. Satellite communication on the non-geostationary system and the geostationary system in the Fixed-satellite service
Yan et al. Interference Analysis of NGSO Constellation to GEO Satellite Communication System Based on Spatio–Temporal Slices
Huangt et al. A GSO protected area calculation model based on controllable NGSO system parameters

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210622

RJ01 Rejection of invention patent application after publication