CN116980030B - Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position - Google Patents
Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position Download PDFInfo
- Publication number
- CN116980030B CN116980030B CN202311229905.1A CN202311229905A CN116980030B CN 116980030 B CN116980030 B CN 116980030B CN 202311229905 A CN202311229905 A CN 202311229905A CN 116980030 B CN116980030 B CN 116980030B
- Authority
- CN
- China
- Prior art keywords
- disturbed
- constellation
- satellite
- seed
- satellites
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013213 extrapolation Methods 0.000 title claims abstract description 33
- 238000004364 calculation method Methods 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
-
- 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 application discloses a method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position, which comprises the steps of obtaining longitude and latitude position information of a disturbed seed satellite by using the orbit extrapolation; dividing the satellite lower points of the disturbed seeds to obtain a grid set; calculating each grid probability; acquiring the position of a satellite lower point of a disturbance seed satellite, the space positions of other satellites of a disturbed constellation and the space positions of other satellites of a Shi Rao constellation; calculating an interference value of the disturbed constellation receiver according to the position parameters of the disturbed constellation receiver, the position parameters of the disturbed constellation transmitter and the position parameters of all transmitters of the disturbed constellation; and traversing the grid set, and calculating probability density distribution of interference values of the disturbed constellation receiver. According to the application, the orbit extrapolation data of seed satellites in two NGSO constellations are used, the positions and probabilities of the two NGSO constellation satellites are calculated by using a space position probability method to obtain final interference value probability distribution, and compared with an orbit extrapolation method, the calculation time is reduced by at least one order of magnitude.
Description
Technical Field
The application relates to the technical field of communication interference calculation among NGSO constellations, in particular to a method for calculating communication interference among NGSO constellations by using orbit extrapolation and spatial position.
Background
Two methods exist for communication interference simulation calculation between satellite constellations. One is to predict the orbit position of the constellation satellite in a period of time by satellite orbit extrapolation, i.e. setting time step, calculate the interference value of each step point, and finally calculate the accumulated time of different interference values to obtain the final interference value probability distribution. The other is to calculate the probability of the position of the constellation satellite by a space position probability, namely setting a space step length, calculating the interference value of each step length point by a mathematical method, and finally counting the accumulated probabilities of different interference values to obtain the final interference value probability distribution.
The satellite orbit extrapolation method has wide application range, can calculate the interference of the NGSO constellation system to the GSO constellation system and also can calculate the interference of the NGSO constellation system to the NGSO constellation system, but has the problems of long simulation time and low calculation efficiency in large-scale and huge constellation scenes; the spatial position probability method is suitable for the scene of NGSO to GSO, because the position probability of one NGSO constellation can be calculated, and the probability of the relative positions of two NGSO constellations can not be calculated.
Patent document CN113595616 discloses an NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization. The method mainly reduces interference among NGSO systems and cannot acquire the probability distribution of interference values of NGSO constellation satellites.
Disclosure of Invention
The application aims to provide a method for calculating communication interference between NGSO constellations by using orbit extrapolation and spatial position, which can combine the orbit extrapolation method and the spatial position probability method to quickly and efficiently calculate the positions and probabilities of all satellites of two NGSO constellations so as to obtain final interference value probability distribution.
The aim of the application can be achieved by the following technical scheme: a method for calculating interference in communication between NGSO constellations using orbital extrapolation and spatial location, comprising the steps of:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, dividing longitude and latitude ranges of points under the disturbed seed satellites according to longitude and latitude step sizes to obtain a grid set;
s3, calculating the probability of each grid according to the number of the points under the disturbed seed satellites in each grid;
s4, acquiring the position of the satellite lower point of the disturbed seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same moment according to the satellite lower point of the disturbed seed satellite;
s5, selecting a transmitter for a receiver of the disturbed constellation system according to the transmitter parameters of the disturbed constellation system; selecting a transmitter for a receiver of the scrambling constellation system according to the transmitter parameters of the scrambling constellation system;
s6, calculating an interference value of the disturbed constellation receiver according to the position and the parameter of the disturbed constellation receiver, the position and the parameter of the disturbed constellation transmitter and the positions and the parameters of all transmitters of the disturbed constellation;
and S7, repeating the steps S4-S6, traversing the grid set, and calculating probability density distribution of interference values of the disturbed constellation receiver.
Further: and S2, acquiring a grid set, wherein the basis of longitude and latitude step division is a deflection angle boundary value of a disturbed constellation receiver, and the deflection angle boundary value is obtained by calculating antenna parameters of the disturbed constellation receiver.
Further: the grid set is divided into a large grid set and a small grid set.
Further: in the step S7, when the step S4-S6 is repeated to traverse the grid set, the step of selecting whether the large grid or the small grid set is traversed is as follows:
s71, calculating and obtaining an interference off-axis angle of the disturbed constellation receiver according to the position of the disturbed constellation receiver, the position of the disturbed constellation transmitter and the positions of all transmitters of the disturbed constellation;
s72, comparing the interference off-axis angle with the off-axis angle boundary value:
if the interference off-axis angle is larger than the off-axis angle boundary value, continuing to traverse by adopting a large grid;
if the interference off-axis angle is smaller than the off-axis angle boundary value, traversing the current large grid by adopting the small grid.
Further: the basis for acquiring the latitude and longitude position information of the satellite point in the S1 is six parameters of the disturbed seed satellite and the disturbed seed satellite.
Further: and S4, the satellite lower points of the disturbed seed satellites in the S4 are in a single grid, and the satellite lower points of the disturbed seed satellites closest to the geometric center of the grid in the satellite lower point set.
Further: in the step S4, the step of obtaining the position of the satellite point of the disturbing seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same time according to the satellite point of the disturbed seed satellite is as follows:
s41, acquiring the position of the satellite point of the disturbing seed satellite at the same moment according to the satellite point of the disturbed seed satellite;
s42, according to constellation configuration parameters of the disturbed constellation; calculating the space positions of other satellites in the disturbed constellation by combining the undersea points of the disturbed seed satellites;
s43, constellation configuration parameters according to a Shi Rao constellation; and calculating the space positions of other satellites of the disturbance constellation by combining the positions of the points under the satellites of the disturbance seed satellites.
The method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position comprises the following steps of:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, determining longitude and latitude ranges of satellite points of the disturbed satellites according to six numbers of the disturbed seed satellites;
s3, acquiring antenna parameters of a disturbed constellation receiver, calculating an off-axis angle boundary value of the disturbed constellation receiver, and determining longitude and latitude space step sizes of the size grids according to the off-axis angle boundary value;
s4, traversing the longitude and latitude range of the satellite lower point of the disturbed seed by taking the longitude and latitude step length of the large grid as the length to obtain a large grid set;
s5, calculating the number of satellite lower points of the disturbed seed satellites falling in each grid according to the satellite lower point position information of the disturbed seed satellites in a period of time, and calculating the probability of each grid according to the number of satellite lower points of the disturbed seed satellites in each grid;
s6, selecting the satellite lower point of the disturbed seed satellite closest to the geometric center of the grid according to the satellite lower point set of the disturbed seed satellite in each grid; according to the satellite point of the disturbed seed satellite, the satellite point position of the disturbed seed satellite at the same moment is obtained;
s7, inputting constellation configuration parameters of a disturbed constellation, and calculating the space positions of other satellites of the disturbed constellation according to the undersea points of the disturbed seed satellites;
s8, inputting transmitter parameters of a disturbed constellation system, and selecting a transmitter for a receiver of the disturbed constellation system;
s9, inputting constellation configuration parameters of a scrambling constellation; calculating the space positions of other satellites of the disturbance constellation according to the undersea points of the disturbance seed satellites;
s10, inputting transmitter parameters of a scrambling constellation system, and selecting a transmitter for a receiver of the scrambling constellation system;
s11, calculating an interference off-axis angle of the disturbed constellation receiver according to the position of the disturbed constellation receiver, the positions of the disturbed constellation transmitters and all transmitter positions of the disturbed constellation;
s12, judging whether the current grid is required to be divided into small grids according to the interference off-axis angle and the off-axis angle boundary value of the interfered constellation receiver;
s13, if the current grid is required to be divided into small grids, traversing the longitude and latitude range of the current grid of the satellite under the disturbed satellite by taking the longitude and latitude step length of the small grid as the length to obtain a small grid set; repeating the steps S5-S10 to traverse the small grid set;
s14, calculating an interference value of the receiver according to the position and the antenna parameter of the receiver of the disturbed constellation, the position of the transmitter of the disturbed constellation, and the positions and the antenna parameters of all transmitters of the disturbed constellation;
s15, traversing the large grid set, and calculating probability density distribution of the interference values according to probability values and the interference values of the large grid set and all the small grid sets.
The application has the beneficial effects that:
1. the application uses the orbit extrapolation data of seed satellites in two NGSO constellations, calculates the interference value and probability of the two NGSO constellation satellites by using a space position probability method to obtain the final interference value probability distribution; the method avoids huge complex orbit extrapolation calculation of all satellites of two NGSO constellations, improves calculation efficiency, and reduces calculation time by at least one order of magnitude compared with the orbit extrapolation method in a large-scale or huge constellation scene.
2. The space positions of other satellites in the disturbed constellation and the space positions of other satellites in the disturbed constellation are obtained by combining the constellation configuration parameters of the disturbed constellation and the constellation configuration parameters of the disturbed constellation, so that the position accuracy is high, and the numerical accuracy obtained by calculating the probability distribution of the interference values among the constellations is high.
3. When the positions and probabilities of two NGSO constellation satellites are calculated to obtain final interference value probability distribution, the large grid or the small grid is judged to be used for traversing according to the interference off-axis angle condition of the interfered constellation receiver in the grid, so that the accuracy of calculating the interference value probability distribution is ensured, and the speed of calculating the interference value probability distribution value is improved.
Drawings
Fig. 1 is a flow chart of a method for calculating interference of communication between NGSO constellations by using orbit extrapolation and spatial position according to the present application.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar symbols indicate like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary only for explaining the present application and are not to be construed as limiting the present application.
The application discloses a method for calculating communication interference between NGSO constellations by using orbit extrapolation and spatial position, which comprises the following steps:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, dividing longitude and latitude ranges of points under the disturbed seed satellites according to longitude and latitude step sizes to obtain a grid set;
s3, calculating the probability of each grid according to the number of the points under the disturbed seed satellites in each grid;
s4, acquiring the position of the satellite lower point of the disturbed seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same moment according to the satellite lower point of the disturbed seed satellite;
s5, selecting a transmitter for a receiver of the disturbed constellation system according to the transmitter parameters of the disturbed constellation system; selecting a transmitter for a receiver of the scrambling constellation system according to the transmitter parameters of the scrambling constellation system;
s6, calculating an interference value of the disturbed constellation receiver according to the position and the parameter of the disturbed constellation receiver, the position and the parameter of the disturbed constellation transmitter and the positions and the parameters of all transmitters of the disturbed constellation;
and S7, repeating the steps S4-S6, traversing the grid set, and calculating probability density distribution of interference values of the disturbed constellation receiver.
The NGSO constellation refers to a non-stationary satellite orbit constellation, the GSO constellation refers to a stationary satellite orbit constellation, the disturbed seed satellites refer to basic satellites in the disturbed NGSO constellation, all satellites in the disturbed NGSO constellation are mutually related, and the accurate positions of other satellites in the disturbed NGSO constellation can be calculated through the disturbed seed satellites; the same scrambling seed satellite refers to a basic satellite in a scrambling NGSO constellation, all satellites in the scrambling NGSO constellation are related to each other, and the accurate positions of other satellites in the Shi Rao NGSO constellation can be calculated through the scrambling seed satellite.
Six parameters of the two seed satellites can be obtained by inquiry, and after the six parameters of the two seed satellites are obtained, the longitude and latitude position information of the satellite points of the two seed satellites in a period of time can be independently calculated by using an orbit extrapolation method.
And determining the latitude and longitude range of the satellite point of the disturbed seed satellite according to the latitude and longitude position information of the disturbed seed satellite.
The antenna parameters of the disturbed constellation receiver can be obtained by inquiry, an antenna gain diagram of the disturbed constellation receiver can be calculated according to the antenna parameters of the disturbed constellation receiver, the antenna gain change condition can be obtained according to the antenna gain diagram, and the off-axis angle for distinguishing the severe antenna gain change and the gentle antenna gain change limit is determined, namely the off-axis angle boundary value.
According to the off-axis angle boundary value, the grid set can be divided into a large grid set and a small grid set, wherein the large grid is adopted in the area larger than the off-axis angle boundary value, and the small grid is adopted in the area smaller than the off-axis angle boundary value.
When performing grid set traversal calculation, carrying out traversal calculation by adopting a large grid by default, judging an interference off-axis angle acquired by the grid when carrying out the large grid traversal calculation, and if the interference off-axis angle is larger than an off-axis angle boundary value, continuing to carry out traversal by adopting the large grid; if the interference off-axis angle is smaller than the off-axis angle boundary value, traversing the current large grid by adopting the small grid.
The interference off-axis angle can be obtained by calculating the receiver position of the disturbed constellation, the transmitter position of the disturbed constellation and all transmitter positions of the disturbed constellation.
The disturbed constellation receiver position may be a disturbed constellation ground receiver station receiver position; the transmitter parameters of the disturbed constellation system can be obtained by inquiry, and a transmitter is selected for the receiver of the disturbed constellation system according to the transmitter parameters of the disturbed constellation system to build a chain for satellite communication, so as to obtain the position of the disturbed constellation transmitter.
Shi Rao constellation receiver positions may be scrambling constellation ground receiver station receiver positions; the transmitter parameters of the scrambling constellation system can be obtained by inquiry, a transmitter is selected for a receiver of the scrambling constellation system according to the transmitter parameters of the scrambling constellation system, satellite communication link establishment is carried out, and the position of the scrambling constellation transmitter is obtained; all transmitter positions of the Shi Rao constellation may be obtained from all scrambling satellite positions in the scrambling constellation.
Preferably, when the position of the satellite lower point of the disturbance seed satellite, the space positions of other satellites of the disturbance constellation and the space position information of other satellites of the disturbance constellation at the same moment are obtained according to the satellite lower point of the disturbance seed satellite, the satellite lower point of the disturbance seed satellite closest to the geometric center of the grid in the grid is preferentially used as the basis of calculation. The satellite lower point of the disturbed seed satellite closest to the geometric center of the grid is used, the accuracy of the satellite lower point position of the disturbed seed satellite, the space positions of other satellites of the disturbed constellation and the space position information of other satellites of the disturbed constellation can be improved, and the accuracy of probability density distribution of the disturbance value of the disturbed constellation receiver obtained integrally is higher.
The step of obtaining the position of the satellite point of the disturbing seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same time may be: firstly, acquiring the position of the satellite point of the disturbing seed satellite at the same moment according to the satellite point of the disturbed seed satellite; then according to constellation configuration parameters of the disturbed constellation; calculating the space positions of other satellites in the disturbed constellation by combining the undersea points of the disturbed seed satellites; then according to constellation configuration parameters of the scrambling constellation; and calculating the space positions of other satellites of the disturbance constellation by combining the positions of the points under the satellites of the disturbance seed satellites.
As shown in fig. 1, the large grid set and the small grid set are used for traversing the longitude and latitude range of the point under the disturbed seed satellite, and the probability density distribution of the disturbance value of the disturbed constellation receiver is calculated, which comprises the following steps:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, determining longitude and latitude ranges of satellite points of the disturbed satellites according to six numbers of the disturbed seed satellites;
s3, acquiring antenna parameters of a disturbed constellation receiver, calculating an off-axis angle boundary value of the disturbed constellation receiver, and determining longitude and latitude space step sizes of the size grids according to the off-axis angle boundary value;
s4, traversing the longitude and latitude range of the satellite lower point of the disturbed seed by taking the longitude and latitude step length of the large grid as the length to obtain a large grid set;
s5, calculating the number of satellite lower points of the disturbed seed satellites falling in each grid according to the satellite lower point position information of the disturbed seed satellites in a period of time, and calculating the probability of each grid according to the number of satellite lower points of the disturbed seed satellites in each grid;
s6, selecting the satellite lower point of the disturbed seed satellite closest to the geometric center of the grid according to the satellite lower point set of the disturbed seed satellite in each grid; according to the satellite point of the disturbed seed satellite, the satellite point position of the disturbed seed satellite at the same moment is obtained;
s7, inputting constellation configuration parameters of a disturbed constellation, and calculating the space positions of other satellites of the disturbed constellation according to the undersea points of the disturbed seed satellites;
s8, inputting transmitter parameters of a disturbed constellation system, and selecting a transmitter for a receiver of the disturbed constellation system;
s9, inputting constellation configuration parameters of a scrambling constellation; calculating the space positions of other satellites of the disturbance constellation according to the undersea points of the disturbance seed satellites;
s10, inputting transmitter parameters of a scrambling constellation system, and selecting a transmitter for a receiver of the scrambling constellation system;
s11, calculating an interference off-axis angle of the disturbed constellation receiver according to the position of the disturbed constellation receiver, the positions of the disturbed constellation transmitters and all transmitter positions of the disturbed constellation;
s12, judging whether the current grid is required to be divided into small grids according to the interference off-axis angle and the off-axis angle boundary value of the interfered constellation receiver;
s13, if the current grid is required to be divided into small grids, traversing the longitude and latitude range of the current grid of the satellite under the disturbed satellite by taking the longitude and latitude step length of the small grid as the length to obtain a small grid set; repeating the steps S5-S10 to traverse the small grid set;
s14, calculating an interference value of the receiver according to the position and the antenna parameter of the receiver of the disturbed constellation, the position of the transmitter of the disturbed constellation, and the positions and the antenna parameters of all transmitters of the disturbed constellation;
s15, traversing the large grid set, and calculating probability density distribution of the interference values according to probability values and the interference values of the large grid set and all the small grid sets.
The foregoing is only a preferred embodiment of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art, who is within the scope of the present application, should make equivalent substitutions or modifications according to the technical scheme of the present application and the inventive concept thereof, and should be covered by the scope of the present application.
It is to be understood that the terms "center," "longitudinal," "transverse," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counter-clockwise," "axial," "radial," "circumferential," and the like are directional or positional relationships as indicated based on the drawings, merely to facilitate describing the application and to simplify the description, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be configured and operated in a particular orientation, and therefore should not be construed as limiting the application.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
In the present application, unless explicitly specified and limited otherwise, the terms "mounted," "connected," "secured," and the like are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally formed; may be mechanically connected, may be electrically connected or may be in communication with each other; either directly or indirectly, through intermediaries, or both, may be in communication with each other or in interaction with each other, unless expressly defined otherwise. The specific meaning of the above terms in the present application can be understood by those of ordinary skill in the art according to the specific circumstances.
In the present application, unless expressly stated or limited otherwise, a first feature "up" or "down" a second feature may be the first and second features in direct contact, or the first and second features in indirect contact via an intervening medium. Moreover, a first feature being "above," "over" and "on" a second feature may be a first feature being directly above or obliquely above the second feature, or simply indicating that the first feature is level higher than the second feature. The first feature being "under", "below" and "beneath" the second feature may be the first feature being directly under or obliquely below the second feature, or simply indicating that the first feature is less level than the second feature.
Claims (8)
1. A method for calculating interference in communication between NGSO constellations using orbital extrapolation and spatial location, comprising the steps of:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, dividing longitude and latitude ranges of points under the disturbed seed satellites according to longitude and latitude step sizes to obtain a grid set;
s3, calculating the probability of each grid according to the number of the points under the disturbed seed satellites in each grid;
s4, acquiring the position of the satellite lower point of the disturbed seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same moment according to the satellite lower point of the disturbed seed satellite;
s5, selecting a transmitter for a receiver of the disturbed constellation system according to the transmitter parameters of the disturbed constellation system; selecting a transmitter for a receiver of the scrambling constellation system according to the transmitter parameters of the scrambling constellation system;
s6, calculating an interference value of the disturbed constellation receiver according to the position and the parameter of the disturbed constellation receiver, the position and the parameter of the disturbed constellation transmitter and the positions and the parameters of all transmitters of the disturbed constellation;
and S7, repeating the steps S4-S6, traversing the grid set, and calculating probability density distribution of interference values of the disturbed constellation receiver.
2. The method for calculating the communication interference between NGSO constellations using orbital extrapolation and spatial location as recited in claim 1, wherein: and S2, acquiring a grid set, wherein the basis of longitude and latitude step division is a deflection angle boundary value of a disturbed constellation receiver, and the deflection angle boundary value is obtained by calculating antenna parameters of the disturbed constellation receiver.
3. The method for calculating the communication interference between NGSO constellations using orbital extrapolation and spatial location as recited in claim 2, wherein: and dividing the grid set into a large grid set and a small grid set according to the off-axis angle boundary value.
4. A method for calculating interference in communication between NGSO constellations using orbital extrapolation and spatial location as claimed in claim 3 wherein: in the step S7, when the step S4-S6 is repeated to traverse the grid set, the step of selecting whether the large grid or the small grid set is traversed is as follows:
s71, calculating and obtaining an interference off-axis angle of the disturbed constellation receiver according to the position of the disturbed constellation receiver, the position of the disturbed constellation transmitter and the positions of all transmitters of the disturbed constellation;
s72, comparing the interference off-axis angle with the off-axis angle boundary value:
if the interference off-axis angle is larger than the off-axis angle boundary value, continuing to traverse by adopting a large grid;
if the interference off-axis angle is smaller than the off-axis angle boundary value, traversing the current large grid by adopting the small grid.
5. The method for calculating the communication interference between NGSO constellations using orbital extrapolation and spatial location as recited in claim 1, wherein: the basis for acquiring the latitude and longitude position information of the satellite point in the S1 is six parameters of the disturbed seed satellite and the disturbed seed satellite.
6. The method for calculating the communication interference between NGSO constellations using orbital extrapolation and spatial location as recited in claim 1, wherein: the satellite lower points of the disturbed seed satellites in the S4 are the satellite lower points of the disturbed seed satellites closest to the geometric center of the grid in the single intra-grid satellite lower point set.
7. The method for calculating the communication interference between NGSO constellations using orbital extrapolation and spatial location as recited in claim 6, wherein: in the step S4, the step of obtaining the position of the satellite point of the disturbing seed satellite, the spatial positions of other satellites of the disturbed constellation and the spatial positions of other satellites of the Shi Rao constellation at the same time according to the satellite point of the disturbed seed satellite is as follows:
s41, acquiring the position of the satellite point of the disturbing seed satellite at the same moment according to the satellite point of the disturbed seed satellite;
s42, according to constellation configuration parameters of the disturbed constellation; calculating the space positions of other satellites in the disturbed constellation by combining the undersea points of the disturbed seed satellites;
s43, constellation configuration parameters according to a Shi Rao constellation; and calculating the space positions of other satellites of the disturbance constellation by combining the positions of the points under the satellites of the disturbance seed satellites.
8. The method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position is characterized in that a large grid set and a small grid set are utilized to traverse the longitude and latitude range of the point under the disturbed seed satellite, and the probability density distribution of the disturbance value of the disturbed constellation receiver is calculated, and the method comprises the following steps:
s1, acquiring a disturbed seed satellite in a disturbed NGSO constellation and longitude and latitude position information of a satellite point of the disturbed seed satellite in the disturbed NGSO constellation within a period of time by using an orbit extrapolation method;
s2, determining longitude and latitude ranges of satellite points of the disturbed satellites according to six numbers of the disturbed seed satellites;
s3, acquiring antenna parameters of a disturbed constellation receiver, calculating an off-axis angle boundary value of the disturbed constellation receiver, and determining longitude and latitude space step sizes of the size grids according to the off-axis angle boundary value;
s4, traversing the longitude and latitude range of the satellite lower point of the disturbed seed by taking the longitude and latitude step length of the large grid as the length to obtain a large grid set;
s5, calculating the number of satellite lower points of the disturbed seed satellites falling in each grid according to the satellite lower point position information of the disturbed seed satellites in a period of time, and calculating the probability of each grid according to the number of satellite lower points of the disturbed seed satellites in each grid;
s6, selecting the satellite lower point of the disturbed seed satellite closest to the geometric center of the grid according to the satellite lower point set of the disturbed seed satellite in each grid; according to the satellite point of the disturbed seed satellite, the satellite point position of the disturbed seed satellite at the same moment is obtained;
s7, inputting constellation configuration parameters of a disturbed constellation, and calculating the space positions of other satellites of the disturbed constellation according to the undersea points of the disturbed seed satellites;
s8, inputting transmitter parameters of a disturbed constellation system, and selecting a transmitter for a receiver of the disturbed constellation system;
s9, inputting constellation configuration parameters of a scrambling constellation; calculating the space positions of other satellites of the disturbance constellation according to the undersea points of the disturbance seed satellites;
s10, inputting transmitter parameters of a scrambling constellation system, and selecting a transmitter for a receiver of the scrambling constellation system;
s11, calculating an interference off-axis angle of the disturbed constellation receiver according to the position of the disturbed constellation receiver, the positions of the disturbed constellation transmitters and all transmitter positions of the disturbed constellation;
s12, judging whether the current grid is required to be divided into small grids according to the interference off-axis angle and the off-axis angle boundary value of the interfered constellation receiver;
s13, if the current grid is required to be divided into small grids, traversing the longitude and latitude range of the current grid of the satellite under the disturbed satellite by taking the longitude and latitude step length of the small grid as the length to obtain a small grid set; repeating the steps S5-S10 to traverse the small grid set;
s14, calculating an interference value of the receiver according to the position and the antenna parameter of the receiver of the disturbed constellation, the position of the transmitter of the disturbed constellation, and the positions and the antenna parameters of all transmitters of the disturbed constellation;
s15, traversing the large grid set, and calculating probability density distribution of the interference values according to probability values and the interference values of the large grid set and all the small grid sets.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311229905.1A CN116980030B (en) | 2023-09-22 | 2023-09-22 | Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311229905.1A CN116980030B (en) | 2023-09-22 | 2023-09-22 | Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116980030A CN116980030A (en) | 2023-10-31 |
| CN116980030B true CN116980030B (en) | 2023-12-15 |
Family
ID=88479933
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311229905.1A Active CN116980030B (en) | 2023-09-22 | 2023-09-22 | Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116980030B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1282154A (en) * | 1999-07-08 | 2001-01-31 | 环球星有限合伙人公司 | Low Earth's orbit distributed net communication system |
| CN111447034A (en) * | 2020-04-16 | 2020-07-24 | 清华大学 | Method for relieving same frequency interference of non-stationary orbit satellite to stationary orbit satellite |
| CA3211802A1 (en) * | 2021-02-26 | 2022-09-01 | China Academy Of Space Technology | Method for reducing frequency interference, and communication satellite system |
| WO2022183477A1 (en) * | 2021-03-05 | 2022-09-09 | 北京小米移动软件有限公司 | Interference processing method and apapratus, communication device and storage medium |
| CN115173974A (en) * | 2022-03-16 | 2022-10-11 | 军事科学院***工程研究院网络信息研究所 | Downlink interference calculation method based on space resource database and contour line interpolation |
| WO2023071716A1 (en) * | 2021-10-29 | 2023-05-04 | 华为技术有限公司 | Satellite communication method and satellite communication device |
| CN116112056A (en) * | 2022-12-12 | 2023-05-12 | 中国空间技术研究院 | Frequency interference avoidance method based on ground station pitching azimuth angle division |
| CN116633424A (en) * | 2023-07-25 | 2023-08-22 | 中国人民解放军战略支援部队航天工程大学 | Walker constellation configuration design method for avoiding large-scale constellation downlink interference |
-
2023
- 2023-09-22 CN CN202311229905.1A patent/CN116980030B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1282154A (en) * | 1999-07-08 | 2001-01-31 | 环球星有限合伙人公司 | Low Earth's orbit distributed net communication system |
| CN111447034A (en) * | 2020-04-16 | 2020-07-24 | 清华大学 | Method for relieving same frequency interference of non-stationary orbit satellite to stationary orbit satellite |
| CA3211802A1 (en) * | 2021-02-26 | 2022-09-01 | China Academy Of Space Technology | Method for reducing frequency interference, and communication satellite system |
| WO2022183477A1 (en) * | 2021-03-05 | 2022-09-09 | 北京小米移动软件有限公司 | Interference processing method and apapratus, communication device and storage medium |
| CN115315903A (en) * | 2021-03-05 | 2022-11-08 | 北京小米移动软件有限公司 | Interference processing method and device, communication equipment and storage medium |
| WO2023071716A1 (en) * | 2021-10-29 | 2023-05-04 | 华为技术有限公司 | Satellite communication method and satellite communication device |
| CN115173974A (en) * | 2022-03-16 | 2022-10-11 | 军事科学院***工程研究院网络信息研究所 | Downlink interference calculation method based on space resource database and contour line interpolation |
| CN116112056A (en) * | 2022-12-12 | 2023-05-12 | 中国空间技术研究院 | Frequency interference avoidance method based on ground station pitching azimuth angle division |
| CN116633424A (en) * | 2023-07-25 | 2023-08-22 | 中国人民解放军战略支援部队航天工程大学 | Walker constellation configuration design method for avoiding large-scale constellation downlink interference |
Non-Patent Citations (1)
| Title |
|---|
| 基于空间位置概率的NGSO通信星座干扰仿真分析研究;李伟，魏文康，刘畅，刘珊杉，武秀广;《电波科学学报》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116980030A (en) | 2023-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101565351B1 (en) | Method and apparatus for geo-locating mobile station | |
| JP7108626B2 (en) | Method and system for locating a terminal in a wireless communication system | |
| JP6290539B2 (en) | An adaptive method for estimating the number of electrons in the ionosphere. | |
| CN104093205A (en) | Method for deploying anchor nodes of wireless positioning system based on received signal strength indication | |
| AU2007260332B2 (en) | Calculation method for network-specific factors in a network of reference stations for a satellite-based positioning system | |
| US6714782B1 (en) | Method for the localizing of terrestrial mobile radio from an aircraft | |
| CN107484240B (en) | Method and device for positioning based on fingerprint | |
| CN110798804B (en) | Indoor positioning method and device | |
| EP1506429B1 (en) | Altitude determination and distribution in cellular communication systems | |
| CN114189809A (en) | Indoor positioning method based on convolutional neural network and high-dimensional 5G observation characteristics | |
| CN104144495B (en) | A kind of fingerprint positioning method based on direction sensor and wlan network | |
| CN113630172B (en) | Uplink interference analysis method based on space resource database and contour line interpolation | |
| CN112558076A (en) | Volume scanning mode calculation method based on networking weather radar coverage area and application | |
| CN112272407A (en) | Indoor Bluetooth positioning method and device based on geographic information system | |
| CN107289973B (en) | A kind of gravitational field suitability judgment method in Gravity Matching navigation | |
| CN116980030B (en) | Method for calculating communication interference between NGSO constellations by using orbit extrapolation and space position | |
| CN113075704B (en) | Satellite signal shielding judgment method | |
| US20040166877A1 (en) | Database for locating wireless terminals based on combinations of signal-strength measurements and geometry-of-arrival measurements | |
| CN115718905A (en) | VTS system-oriented multi-sensor information fusion method | |
| CN103167528B (en) | A kind of frequency allocation assessment, Apparatus and system | |
| ES2795274T3 (en) | Procedure and distributed location system PHDOA, PWDOA, of emission sources | |
| CN116248163B (en) | Frequency compatibility analysis method for large-scale low-orbit constellation | |
| CN114466451B (en) | Positioning method, positioning device and storage medium based on drive test frequency spectrum data | |
| CN111447549B (en) | Non-uniform UWB positioning error set network construction method and positioning error modeling method | |
| CN116471547A (en) | Arrival time positioning non-line-of-sight anchor screening method based on channel knowledge map |
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 |