A method of NGSO inter-constellation interference avoidance

Abstract

The application discloses a kind of NGSO constellation interferences avoidance methods, the method includes: the satellite ephemeris of target NGSO constellation and interference NGSO constellation is updated, satellite ephemeris information is obtained;Satellite ephemeris information is processed, and NGSO constellation interferences information is obtained;Use the NGSO constellation information, construct useful signal model and interference signal model;Use the useful signal model and the interference signal model, calculate the expected value of target NGSO ground station in each available communication satellite remaining communication duration;The expected value is the ratio of the mathematical expectation of useful signal and the sum of interference signal and noise power;The expected value is processed, and NGSO constellation interferences avoidance model is obtained.The method of the application can comprehensively, accurately analyze and select suitable satellite to establish communication link to avoid NGSO constellation interferences.

Description

Technical Field

[0001] This invention relates to the field of satellite communication technology, and in particular to a method for avoiding interference between NGSO constellations. Background Technology

[0002] In recent years, with the development and innovation of technologies such as "multiple satellites launched on a single rocket," "rocket recovery," and "satellite miniaturization," commercial spaceflight has once again emerged globally. Satellite communication technology is no longer limited to single geostationary orbit (GSO) satellites. Faced with urgent needs for lower communication latency, lower communication link loss, faster data transmission speeds, larger system capacity, and lower launch costs, non-geostationary orbit (NGSO) mega-satellite communication systems have become the preferred choice for major domestic and international telecommunications operators and commercial technology companies. They have launched the construction of mega-satellite communication constellation systems composed of hundreds or thousands of NGSO satellites. With the implementation of these NGSO mega-satellite communication constellation system plans, space frequency and orbit resources will become even more scarce. Co-frequency interference (CFI) has become an unavoidable problem restricting the development of mega-NGSO constellations. Existing interference avoidance methods for mega-constellations require spectrum sensing capabilities and are difficult to cope with complex and time-varying interference environments. Summary of the Invention

[0003] The technical problem this invention aims to solve is to address the complexity of downlink interference scenarios between NGSO constellations and the difficulty of interference avoidance. It provides an interference avoidance method for NGSO constellations, proposing an intelligent satellite selection strategy that maximizes the expected signal-to-interference-plus-noise ratio (MESINR). The interfered earth station establishes a communication link by selecting the satellite with the highest expected MESINR within its line of sight, thus improving the communication quality for users of the target constellation. This invention considers the ephemeris of both the target and interfering NGSO constellations, accurately modeling the useful and interfering signals of the NGSO satellite constellations, and can mitigate downlink interference between different NGSO constellations to a certain extent.

[0004] To address the aforementioned technical problems, this invention discloses an inter-constellation interference avoidance method for NGSO systems, the method comprising:

[0005] S1, update the ephemeris of the target NGSO constellation and the interfering NGSO constellation to obtain satellite ephemeris information;

[0006] The satellite ephemeris information includes the satellite orbital inclination i, the right ascension of the ascending node Ω0 at the time of TLE data update, the satellite orbital eccentricity e, the satellite orbital perigee argument E0 at the time of TLE data update, and the average number of orbits N of the satellite per day. s ;

[0007] S2, process the satellite ephemeris information to obtain NGSO constellation inter-constellation interference information;

[0008] The inter-NGSO constellation interference information includes the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation, the remaining visibility time of each communication satellite, and the number of interfering satellites visible to the target NGSO constellation ground station from the interfering NGSO constellation.

[0009] S3, using the NGSO constellation information, construct a useful signal model and an interference signal model;

[0010] The useful signal model is a model of the useful signals received by the target NGSO constellation ground station from the target NGSO satellite;

[0011] The interference signal model is a model of the ground station of the target NGSO constellation receiving interference signals from interfering NGSO satellites;

[0012] S4. Using the useful signal model and the interference signal model, calculate the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite;

[0013] The expected value is the mathematical expectation of the ratio of the useful signal to the sum of the interference signal and noise power;

[0014] S5, process the expected value to obtain the NGSO constellation interference avoidance model.

[0015] As an optional implementation, in this embodiment of the invention, processing the satellite ephemeris information to obtain NGSO inter-constellation interference information includes:

[0016] S21, using the satellite orbital period calculation model, calculate the average number of times N orbits the Earth per day for the satellite. s The satellite's orbital period is obtained through processing.

[0017] The expression for the satellite orbital period calculation model is as follows:

[0018]

[0019] Among them, T s Let be the satellite's orbital period, and 'a' be the satellite's semi-major axis. T eThe value represents the Earth's rotation period, taken as 86164 s, and μ is the Earth's gravitational constant, which is 3.986005 × 10⁻⁶. 14 m 3 / s 2 ;

[0020] S22, Process the satellite's orbital period to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system;

[0021] S23, The coordinates of the NGSO satellite in the geocentric fixed coordinate system are processed to obtain the interference information between the NGSO constellations.

[0022] As an optional implementation, in this embodiment of the invention, processing the satellite's orbital period to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system includes:

[0023] S221, Process the orbital period of the satellite to obtain the average angular velocity of the NGSO satellite;

[0024] The average angular velocity of the NGSO satellite is:

[0025]

[0026] The average angular velocity of the NGSO satellite is ω;

[0027] S222, The average angular velocity of the NGSO satellite is processed to obtain the latitude parameter of the NGSO satellite at the calculation time t;

[0028] The latitude parameter E of the NGSO satellite at time t is calculated. t The expression is:

[0029] E t =E0+ω×t

[0030] Where E0 is the perigee argument of the satellite orbit at the time of TLE data update, and t is the time elapsed between the calculation time and the time of TLE data update;

[0031] S223, Process the satellite orbital inclination i to obtain the NGSO satellite orbital precession rate;

[0032] The expression for the orbital precession rate of the NGSO satellite is:

[0033]

[0034] Among them, Ω r J is the orbital precession rate of the NGSO satellite, and J2 is the second-order gravitational potential coefficient, with a value of 1082.6 × 10⁻⁶. -6 Re denoted as Earth's radius, taken as 6,378,000 m; denoted as a, the semi-major axis of the satellite; and denoted as μ, the Earth's gravitational constant.

[0035] S224, Process the orbital precession rate of the NGSO satellite to obtain the ascending right ascension of the NGSO satellite at the calculation time t;

[0036] The expression for the ascending right ascension of the NGSO satellite at calculation time t is:

[0037] Ω t =Ω0+Ω r ×t

[0038] Among them, Ω t The ascending right ascension of the NGSO satellite at time t;

[0039] S225, Process the ascending right ascension of the NGSO satellite at the calculation time t to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system;

[0040] The coordinate expression of the NGSO satellite in the geocentric fixed coordinate system is as follows:

[0041]

[0042] in, S represents the j-th satellite in the NGSO constellation. j The coordinate vector of time t at a fixed reference center at the Earth's center

[0043] x j,t y j,t z j,t For the j-th satellite S in the NGSO constellation j The coordinate values ​​of the coordinate vector at time t at a fixed reference center in the Earth's center in the x, y, and z directions, E j,t For the j-th satellite S j Latitude parameter at time t, Ω j,t For the j-th satellite S j At time t, the right ascension of the ascendant occurs.

[0044] As an optional implementation, in this embodiment of the invention, processing the coordinates of the NGSO satellite in the geocentric fixed coordinate system to obtain inter-NGSO constellation interference information includes:

[0045] S231, Process the latitude and longitude information of the target NGSO satellite ground station G to obtain the coordinate vector of the target NGSO ground station in a fixed reference system centered on the earth;

[0046] The expression for the coordinate vector of the target NGSO ground station in the Earth-centered fixed reference frame is:

[0047]

[0048] in, Let G be the coordinate vector of the target NGSO ground station in a fixed reference frame centered on Earth. The latitude and longitude information of the target NGSO satellite ground station G is (lat, lon), where lat represents the latitude of the target NGSO satellite ground station, lon represents the longitude of the target NGSO satellite ground station, and R... e Ω is the Earth's radius, t is the time elapsed between the calculation time and the TLE data update time, and Ω is the Earth's radius. e This represents the Earth's angular velocity of rotation at the equator, with a value of 7.29 × 10⁻⁶. -5 rad / s, x g,t y g,t z g,t These are the coordinates of the target NGSO satellite ground station G in the x, y, and z directions in the geocentric fixed reference frame;

[0049] S232, Process the coordinate vector of the target NGSO ground station in the fixed reference frame centered on the earth to obtain the antenna elevation angle of the target NGSO constellation ground station;

[0050] The antenna elevation angle expression for the target NGSO constellation ground station is as follows:

[0051]

[0052] in Represents vector The model, The distance between the target NGSO constellation ground station G and the NGSO constellation satellites. S represents the j-th satellite in the NGSO constellation. j The coordinate vector at time t is the fixed reference center at the Earth's center. ε is the antenna elevation angle of the target NGSO constellation ground station. The antenna elevation angle constraint of the target NGSO constellation ground station is ε≥ε0, where ε0 represents the minimum antenna elevation angle of the target NGSO constellation ground station.

[0053] S233, The coordinate vector of the target NGSO ground station in the fixed reference frame centered on the earth and the antenna elevation angle of the target NGSO constellation ground station are processed to obtain the NGSO constellation interference information;

[0054] The inter-NGSO constellation interference information includes the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation, the remaining visibility time of each communication satellite, and the number of interfering NGSO constellations visible to the target NGSO constellation ground station from the interfering NGSO constellation.

[0055] As an optional implementation, in this embodiment of the invention, the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation is the set of times t. The number N ob1,t ;

[0056] The set To meet A collection of satellites, among which, For the target NGSO constellation N c Of the N satellites, N is the only one capable of establishing a communication link with the target NGSO ground station. c The target NGSO constellation satellite total number, Represents vector The model, For the target NGSO constellation ground station G to the target NGSO constellation satellite The distance between them This indicates the target NGSO constellation's j-th satellite. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0057] As an optional implementation, in this embodiment of the invention, the remaining visible time of each communication satellite is expressed as follows:

[0058]

[0059] in, For time t, for N ob1,t There are 10 available communication satellites, and each satellite does not meet the requirements. At that time, the remaining visible time for each communication satellite was...

[0060] As an optional implementation, in this embodiment of the invention, the number N of the interfering NGSO constellation visible to the target NGSO constellation ground station is... ob2,t For time t set Quantity;

[0061] The set To meet A collection of interfering satellites To interfere with the NGSO constellation N interf N satellites that interfered with the target NGSO ground station interf To interfere with the total number of NGSO constellation satellites, Represents vector The model, To target the NGSO constellation ground station G to interfere with the NGSO constellation satellites The distance between them This indicates interference with the j-th satellite in the NGSO constellation. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0062] As an optional implementation, in this embodiment of the invention, the useful signal model expression is:

[0063]

[0064] in, For time t, the target NGSO constellation ground station G is within line of sight of available communication satellites. The expected signal power of the transmitter at the ground station G receiver. The transmit power of the j-th available satellite in the target NGSO constellation within the line of sight of the ground station G of the target NGSO constellation is expressed in W. This is the radiation pattern of the target NGSO constellation satellite antenna, with units of relative ratios; 0 indicates that the target NGSO constellation ground station G is beam-aligned with the satellite during communication. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; λ j The carrier wavelength of the j-th satellite in the NGSO constellation system within the line of sight of the ground station G of the target NGSO constellation is represented in meters (m). The distance between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t is expressed in meters.

[0065] As an optional implementation, in this embodiment of the invention, the useful interference signal model expression is:

[0066]

[0067] Among them, I j,t Let be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of ground station G of the target constellation system at time t, received at the target NGSO ground station G. The equivalent transmit power of the j-th potential interfering satellite of the interfering NGSO constellation within the line of sight of the target NGSO constellation ground station G in the overlapping frequency band is expressed in W. This represents the radiation pattern of the interfering NGSO constellation satellite antennas, expressed as a relative ratio. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; φ j,t θ represents the angle between the antenna of the j-th interfering NGSO constellation and the target NGSO constellation ground station G at time t. j,tλ represents the angle between the communication link of the target NGSO constellation ground station G and the j-th interfering satellite of the interfering NGSO constellation at time t; j The wavelength of the signal from the j-th satellite interfering with the NGSO constellation, within the line of sight of the target NGSO constellation ground station G, is expressed in meters (m). The distance, in meters (m), represents the distance between the target NGSO constellation ground station G and the j-th satellite of the interfering NGSO constellation within its line of sight at time t.

[0068] O represents the geocenter corresponding to the geocentric fixed coordinate system. D represents the distance, in meters, between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t; j,t This represents the distance between the j-th potential interfering satellite of the NGSO constellation and the Earth's center within the line of sight of ground station G.

[0069] in Interfering with NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame. The target NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame.

[0070] As an optional implementation, in this embodiment of the invention, the expected value expression of the target NGSO ground station within the remaining communication time of each available communication satellite is:

[0071]

[0072] Among them, E(SINR) j,t SINR is the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite. j,t This represents the signal-to-interference-plus-noise ratio (SIR) of the j-th satellite within the line of sight of the target NGSO constellation ground station at time t. t represents the signal power received by the target NGSO constellation ground station from the j-th available satellite within its line of sight at time t; t0 represents the time when the target NGSO constellation ground station establishes a communication link or switches to a communication satellite. This represents the remaining communication time of the j-th available communication satellite in the target NGSO constellation at time t; This represents the sum of interference signal power received by the target NGSO constellation ground station from potential interfering satellites within its line of sight at time t; N represents the equivalent noise of the interfered earth station receiver, which is Gaussian white noise.ob2,t To interfere with the number of visible jamming satellites of the NGSO constellation ground station, I j,t Let t be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of the ground station G of the target constellation system at time t, and the signal power of the receiver at the target NGSO ground station G.

[0073] Compared with the prior art, the embodiments of the present invention have the following beneficial effects:

[0074] (1) This invention proposes an analysis method for the number of target NGSO constellation communication satellites visible from the target NGSO constellation ground station and the remaining available communication time. The visible range of the target NGSO constellation ground station and the NGSO constellation orbit are modeled into the geocentric rectangular coordinate system. The visible range of communication satellites is represented by solid geometry, which can vividly and quickly determine the number of visible communication satellites and the remaining communication time.

[0075] (2) This invention proposes a model for the ground station of the target NGSO constellation to receive useful signals from the target NGSO satellite and a model for the ground station of the target NGSO constellation to receive interference signals from the interfering NGSO satellite. Under free space link attenuation, the antenna models of the ground station of the target NGSO constellation, the target NGSO constellation satellite and the interfering NGSO constellation satellite are considered simultaneously. The power of the useful signal of the target NGSO constellation and the interference signal of the interfering NGSO constellation at the receiver of the ground station of the target NGSO constellation can be calculated quickly.

[0076] (3) This invention proposes an inter-satellite interference avoidance method for NGSO constellations based on maximizing the expected signal-to-interference-plus-noise ratio. The remaining communication time of each available communication satellite within the visible range of the ground station of the target NGSO constellation is divided in time. The expected signal-to-interference-plus-noise ratio of each satellite is compared within the remaining communication time of each available satellite. The optimal satellite can be selected from the available communication satellites for communication, thereby mitigating the co-frequency interference of interfering NGSO constellation satellites on the downlink of the target NGSO constellation. Attached Figure Description

[0077] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0078] Figure 1 This is a flowchart illustrating an NGSO inter-constellation interference avoidance method disclosed in an embodiment of the present invention;

[0079] Figure 2This is a flowchart illustrating another NGSO constellation interference avoidance method disclosed in an embodiment of the present invention. Detailed Implementation

[0080] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0081] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this invention are used to distinguish different objects, not to describe a specific order. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, apparatus, product, or device that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or devices.

[0082] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of the invention. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0083] This invention discloses a method for avoiding inter-NGSO constellation interference. The method includes: updating the satellite ephemeris of the target NGSO constellation and the interfering NGSO constellation to obtain satellite ephemeris information; processing the satellite ephemeris information to obtain inter-NGSO constellation interference information; using the NGSO constellation information to construct a useful signal model and an interference signal model; the useful signal model is a model of the target NGSO constellation ground station receiving useful signals from the target NGSO satellite; the interference signal model is a model of the target NGSO constellation ground station receiving interference signals from the interfering NGSO satellite; using the useful signal model and the interference signal model, calculating the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite; the expected value is the mathematical expectation of the ratio of the useful signal to the sum of the interference signal and noise power; processing the expected value to obtain the inter-NGSO constellation interference avoidance model. This invention's method can comprehensively and accurately analyze and select suitable satellites to establish communication links to avoid co-channel interference between NGSO constellations. Detailed explanations follow.

[0084] Example 1

[0085] Please see Figure 1 , Figure 1 This is a flowchart illustrating an NGSO inter-constellation interference avoidance method disclosed in an embodiment of the present invention. Figure 1 The described NGSO inter-constellation interference avoidance method is applied in the field of satellite communication technology to achieve NGSO inter-constellation interference avoidance. The embodiments of this invention are not limited thereto. Figure 1 As shown, the NGSO inter-constellation interference avoidance method may include the following operations:

[0086] S1, update the ephemeris of the target NGSO constellation and the interfering NGSO constellation to obtain satellite ephemeris information;

[0087] The satellite ephemeris information includes the satellite orbital inclination i, the right ascension of the ascending node Ω0 at the time of TLE data update, the satellite orbital eccentricity e, the satellite orbital perigee argument E0 at the time of TLE data update, and the average number of orbits N of the satellite per day. s ;

[0088] S2, process the satellite ephemeris information to obtain NGSO constellation inter-constellation interference information;

[0089] The inter-NGSO constellation interference information includes the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation, the remaining visibility time of each communication satellite, and the number of interfering satellites visible to the target NGSO constellation ground station from the interfering NGSO constellation.

[0090] S3, using the NGSO constellation information, construct a useful signal model and an interference signal model;

[0091] The useful signal model is a model of the useful signals received by the target NGSO constellation ground station from the target NGSO satellite;

[0092] The interference signal model is a model of the ground station of the target NGSO constellation receiving interference signals from interfering NGSO satellites;

[0093] S4. Using the useful signal model and the interference signal model, calculate the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite;

[0094] The expected value is the mathematical expectation of the ratio of the useful signal to the sum of the interference signal and noise power;

[0095] S5, process the expected value to obtain the NGSO constellation interference avoidance model.

[0096] Optionally, the processing of the satellite ephemeris information to obtain NGSO inter-constellation interference information includes:

[0097] S21, using the satellite orbital period calculation model, calculate the average number of times N orbits the Earth per day for the satellite. s The satellite's orbital period is obtained through processing.

[0098] The expression for the satellite orbital period calculation model is as follows:

[0099]

[0100] Among them, T s Let be the satellite's orbital period, and 'a' be the satellite's semi-major axis. T e The value represents the Earth's rotation period, taken as 86164 s, and μ is the Earth's gravitational constant, which is 3.986005 × 10⁻⁶. 14 m 3 / s 2 ;

[0101] S22, Process the satellite's orbital period to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system;

[0102] S23, The coordinates of the NGSO satellite in the geocentric fixed coordinate system are processed to obtain the interference information between the NGSO constellations.

[0103] Optionally, processing the satellite's orbital period to obtain the NGSO satellite's coordinates in a geocentric fixed coordinate system includes:

[0104] S221, Process the orbital period of the satellite to obtain the average angular velocity of the NGSO satellite;

[0105] The average angular velocity of the NGSO satellite is:

[0106]

[0107] The average angular velocity of the NGSO satellite is ω;

[0108] S222, The average angular velocity of the NGSO satellite is processed to obtain the latitude parameter of the NGSO satellite at the calculation time t;

[0109] The latitude parameter E of the NGSO satellite at time t is calculated. t The expression is:

[0110] E t =E0+ω×t

[0111] Where E0 is the perigee argument of the satellite orbit at the time of TLE data update, and t is the time elapsed between the calculation time and the time of TLE data update;

[0112] S223, Process the satellite orbital inclination i to obtain the NGSO satellite orbital precession rate;

[0113] The expression for the orbital precession rate of the NGSO satellite is:

[0114]

[0115] Among them, Ω r J is the orbital precession rate of the NGSO satellite, and J2 is the second-order gravitational potential coefficient, with a value of 1082.6 × 10⁻⁶. -6 R e denoted as Earth's radius, taken as 6,378,000 m; denoted as a, the semi-major axis of the satellite; and denoted as μ, the Earth's gravitational constant.

[0116] S224, Process the orbital precession rate of the NGSO satellite to obtain the ascending right ascension of the NGSO satellite at the calculation time t;

[0117] The expression for the ascending right ascension of the NGSO satellite at calculation time t is:

[0118] Ω t =Ω0+Ω r ×t

[0119] Among them, Ω t The ascending right ascension of the NGSO satellite at time t;

[0120] S225, Process the ascending right ascension of the NGSO satellite at the calculation time t to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system;

[0121] The coordinate expression of the NGSO satellite in the geocentric fixed coordinate system is as follows:

[0122]

[0123] in, S represents the j-th satellite in the NGSO constellation. j The coordinate vector of time t at a fixed reference center at the Earth's center

[0124] x j,t y j,t z j,t For the j-th satellite S in the NGSO constellation j The coordinate values ​​of the coordinate vector at time t at a fixed reference center in the Earth's center in the x, y, and z directions, E j,t For the j-th satellite S jLatitude parameter at time t, Ω j,t For the j-th satellite S j At time t, the right ascension of the ascendant occurs.

[0125] Optionally, the process of processing the coordinates of the NGSO satellite in the geocentric fixed coordinate system to obtain inter-NGSO constellation interference information includes:

[0126] S231, Process the latitude and longitude information of the target NGSO satellite ground station G to obtain the coordinate vector of the target NGSO ground station in a fixed reference system centered on the earth;

[0127] The expression for the coordinate vector of the target NGSO ground station in the Earth-centered fixed reference frame is:

[0128]

[0129] in, Let G be the coordinate vector of the target NGSO ground station in a fixed reference frame centered on Earth. The latitude and longitude information of the target NGSO satellite ground station G is (lat, lon), where lat represents the latitude of the target NGSO satellite ground station, lon represents the longitude of the target NGSO satellite ground station, and R... e Ω is the Earth's radius, t is the time elapsed between the calculation time and the TLE data update time, and Ω is the Earth's radius. e This represents the Earth's angular velocity of rotation at the equator, with a value of 7.29 × 10⁻⁶. -5 rad / s, x g,t y g,t z g,t These are the coordinates of the target NGSO satellite ground station G in the x, y, and z directions in the geocentric fixed reference frame;

[0130] S232, Process the coordinate vector of the target NGSO ground station in the fixed reference frame centered on the earth to obtain the antenna elevation angle of the target NGSO constellation ground station;

[0131] The antenna elevation angle expression for the target NGSO constellation ground station is as follows:

[0132]

[0133] in Represents vector The model, The distance between the target NGSO constellation ground station G and the NGSO constellation satellites. S represents the j-th satellite in the NGSO constellation. jThe coordinate vector at time t is the fixed reference center at the Earth's center. ε is the antenna elevation angle of the target NGSO constellation ground station. The antenna elevation angle constraint of the target NGSO constellation ground station is ε≥ε0, where ε0 represents the minimum antenna elevation angle of the target NGSO constellation ground station.

[0134] S233, The coordinate vector of the target NGSO ground station in the fixed reference frame centered on the earth and the antenna elevation angle of the target NGSO constellation ground station are processed to obtain the NGSO constellation interference information;

[0135] The inter-NGSO constellation interference information includes the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation, the remaining visibility time of each communication satellite, and the number of interfering NGSO constellations visible to the target NGSO constellation ground station from the interfering NGSO constellation.

[0136] Optionally, the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation is the set at time t. The number N ob1,t ;

[0137] The set To meet A collection of satellites, among which, For the target NGSO constellation N c Of the N satellites, N is the only one capable of establishing a communication link with the target NGSO ground station. c The target NGSO constellation satellite total number, Represents vector The model, For the target NGSO constellation ground station G to the target NGSO constellation satellite The distance between them This indicates the target NGSO constellation's j-th satellite. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0138] Optionally, the remaining visible time expression for each communication satellite is:

[0139]

[0140] in, For time t, for N ob1,t There are 10 available communication satellites, and each satellite does not meet the requirements. At that time, the remaining visible time for each communication satellite was...

[0141] Optionally, the number N of the interfering NGSO constellation visible to the target NGSO constellation ground station is... ob2,t For time t set Quantity;

[0142] The set To meet A collection of interfering satellites To interfere with the NGSO constellation N interf N satellites that interfered with the target NGSO ground station interf To interfere with the total number of NGSO constellation satellites, Represents vector The model, To target the NGSO constellation ground station G to interfere with the NGSO constellation satellites The distance between them This indicates interference with the j-th satellite in the NGSO constellation. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0143] Optionally, the useful signal model expression is:

[0144]

[0145] in, For time t, the target NGSO constellation ground station G is within line of sight of available communication satellites. The expected signal power of the transmitter at the ground station G receiver. The transmit power of the j-th available satellite in the target NGSO constellation within the line of sight of the ground station G of the target NGSO constellation is expressed in W. This is the radiation pattern of the target NGSO constellation satellite antenna, with units of relative ratios; 0 indicates that the target NGSO constellation ground station G is beam-aligned with the satellite during communication. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; λ j The carrier wavelength of the j-th satellite in the NGSO constellation system within the line of sight of the ground station G of the target NGSO constellation is represented in meters (m). The distance between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t is expressed in meters.

[0146] Optionally, the useful interference signal model expression is:

[0147]

[0148] Among them, I j,t Let be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of ground station G of the target constellation system at time t, received at the target NGSO ground station G. The equivalent transmit power of the j-th potential interfering satellite of the interfering NGSO constellation within the line of sight of the target NGSO constellation ground station G in the overlapping frequency band is expressed in W. This represents the radiation pattern of the interfering NGSO constellation satellite antennas, expressed as a relative ratio. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; φ j,t θ represents the angle between the antenna of the j-th interfering NGSO constellation and the target NGSO constellation ground station G at time t. j,t λ represents the angle between the communication link of the target NGSO constellation ground station G and the j-th interfering satellite of the interfering NGSO constellation at time t; j The wavelength of the signal from the j-th satellite interfering with the NGSO constellation, within the line of sight of the target NGSO constellation ground station G, is expressed in meters (m). The distance, in meters (m), represents the distance between the target NGSO constellation ground station G and the j-th satellite of the interfering NGSO constellation within its line of sight at time t.

[0149] O represents the geocenter corresponding to the geocentric fixed coordinate system. D represents the distance, in meters, between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t; j,t This represents the distance between the j-th potential interfering satellite of the NGSO constellation and the Earth's center within the line of sight of ground station G.

[0150] in Interfering with NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame. The target NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame.

[0151] Optionally, the expected value expression for the target NGSO ground station within the remaining communication time of each available communication satellite is:

[0152]

[0153] Among them, E(SINR) j,t SINR is the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite. j,t This represents the signal-to-interference-plus-noise ratio (SIR) of the j-th satellite within the line of sight of the target NGSO constellation ground station at time t. t represents the signal power received by the target NGSO constellation ground station from the j-th available satellite within its line of sight at time t; t0 represents the time when the target NGSO constellation ground station establishes a communication link or switches to a communication satellite. This represents the remaining communication time of the j-th available communication satellite in the target NGSO constellation at time t; This represents the sum of interference signal power received by the target NGSO constellation ground station from potential interfering satellites within its line of sight at time t; N represents the equivalent noise of the interfered earth station receiver, which is Gaussian white noise. ob2,t To interfere with the number of visible jamming satellites of the NGSO constellation ground station, I j,t Let t be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of the ground station G of the target constellation system at time t, and the signal power of the receiver at the target NGSO ground station G.

[0154] As can be seen, this invention proposes an analysis method for the number of visible target NGSO constellation communication satellites and their remaining usable communication time from a target NGSO constellation ground station. It models the visibility range of the target NGSO constellation ground station and the NGSO constellation orbit into a geocentric rectangular coordinate system, representing the visible communication satellite range through solid geometry. This allows for a clear and quick determination of the number of visible communication satellites and the remaining communication time. Furthermore, this invention proposes a model for the target NGSO constellation ground station receiving useful signals from target NGSO satellites and a model for the target NGSO constellation ground station receiving interference signals from interfering NGSO satellites. Under free-space link attenuation, it simultaneously considers the target NGSO satellites... The invention employs antenna models of a ground station, target NGSO constellation satellites, and interfering NGSO constellation satellites to quickly calculate the power of the useful signal of the target NGSO constellation and the interfering signal of the interfering NGSO constellation at the receiver of the target NGSO constellation ground station. The invention divides time by the remaining communication time of each available communication satellite within the line of sight of the target NGSO constellation ground station. By comparing the expected signal-to-interference-plus-noise ratio (SNR) of each satellite within its remaining communication time, the optimal satellite can be selected for communication, thereby mitigating co-channel interference of the interfering NGSO constellation satellites on the downlink of the target NGSO constellation.

[0155] Example 2

[0156] Please see Figure 2 , Figure 2 This is a flowchart illustrating another NGSO inter-constellation interference avoidance method disclosed in an embodiment of the present invention. Figure 2 The described NGSO inter-constellation interference avoidance method is applied in the field of satellite communication technology to achieve NGSO inter-constellation interference avoidance. The embodiments of the present invention are not limited.

[0157] like Figure 2 As shown, the NGSO inter-constellation interference avoidance method may include the following operations:

[0158] Step S1: Update the ephemeris of the target NGSO constellation and the interfering NGSO constellation satellites;

[0159] The target NGSO constellation is defined as the NGSO constellation that is being interfered with, the target NGSO constellation ground station is defined as the NGSO constellation ground station that is being interfered with, and the interfering NGSO constellation is defined as the NGSO constellation that is different from the target NGSO constellation and interferes with the target NGSO constellation ground station.

[0160] Satellite ephemeris, also known as two-line orbital elements (TLE), is used to accurately calculate, predict, depict, and track the time, position, velocity, and other operational status of satellites, and is publicly available. Ground stations of the target NGSO constellation can accurately track the positions of satellites in both the target and interfering NGSO constellations based on updated TLE data. TLE data includes the satellite's orbital inclination i, the right ascension of the ascending node Ω0 at the time of the TLE update, the satellite's orbital eccentricity e, the argument of perigee E0 at the time of the TLE update, and the average number of orbits N per day. s ;

[0161] Based on the average number of times N satellites orbit the Earth per day s Calculate the satellite orbital period T s (in seconds, s) and the satellite's semi-major axis a (in meters, m).

[0162]

[0163] therefore: Where T e The value represents the Earth's rotation period, taken as 86164 s; μ represents the Earth's gravitational constant, 3.986005 × 10⁻⁶. 14 m 3 / s 2 ;

[0164] Calculate the coordinates of the NGSO satellite in the geocentric fixed coordinate system;

[0165] The average angular velocity of the NGSO satellite is The unit is rad / s;

[0166] The latitude parameter E of the NGSO satellite at time t. t =E0+ω×t, where t represents the time elapsed between the calculation time and the TLE data update time;

[0167] NGSO satellite orbital precession rate Where J2 represents the second-order gravitational potential coefficient, with a value of 1082.6 × 10⁻⁶. -6 , i represents the satellite orbital inclination, R e The Earth's radius is 6,378,000 m, and 'a' represents the satellite's semi-major axis.

[0168] The ascending right ascension Ω of the NGSO satellite at time t. t =Ω0+Ω r ×t, where t represents the time elapsed between the calculation time and the TLE data update time;

[0169] NGSO constellation i-th satellite S i The coordinate vector of time t at a fixed reference center at the Earth's center

[0170]

[0171] in, S represents the j-th satellite in the NGSO constellation. j The coordinate vector of time t at a fixed reference center at the Earth's center

[0172] x j,t y j,t z j,t For the j-th satellite S in the NGSO constellation j The coordinate values ​​of the coordinate vector at time t at a fixed reference center in the Earth's center in the x, y, and z directions, E j,t For the j-th satellite S j Latitude parameter at time t, Ω j,t For the j-th satellite S j At time t, the right ascension of the ascendant occurs.

[0173] Step S2: Analyze the number of communication satellites visible to the target NGSO constellation ground station from the target NGSO constellation, the remaining visibility time of each communication satellite, and the number of interfering satellites visible to the target NGSO constellation ground station from the interfering NGSO constellation.

[0174] The visible range of the target NGSO satellite ground station on the NGSO satellite orbital plane is expressed;

[0175] Let the latitude and longitude of the target NGSO satellite ground station G be (lat, lon), where lat represents the latitude of the target NGSO satellite ground station and lon represents the longitude of the target NGSO satellite ground station; the coordinate vector of the target NGSO ground station in a fixed reference frame centered on the Earth. for;

[0176]

[0177] in, Let G be the coordinate vector of the target NGSO ground station in a fixed reference frame centered on Earth. The latitude and longitude information of the target NGSO satellite ground station G is (lat, lon), where lat represents the latitude of the target NGSO satellite ground station, lon represents the longitude of the target NGSO satellite ground station, and R... e Ω is the Earth's radius, t is the time elapsed between the calculation time and the TLE data update time, and Ω is the Earth's radius. e This represents the Earth's angular velocity of rotation at the equator, with a value of 7.29 × 10⁻⁶. -5 rad / s, x g,t y g,t z g,t These are the coordinates of the target NGSO satellite ground station G in the x, y, and z directions in the geocentric fixed reference frame;

[0178] The antenna elevation angle constraint of the target NGSO constellation ground station is ε≥ε0. The target NGSO constellation ground station can establish a communication link with any target NGSO constellation satellite that meets the antenna elevation angle constraint. The target constellation ground station G and the NGSO constellation satellite S i The angle between the line connecting the antenna and the horizontal plane is the antenna elevation angle.

[0179]

[0180] in Represents vector The model, The distance between the target NGSO constellation ground station G and the NGSO constellation satellites. S represents the j-th satellite in the NGSO constellation. j The coordinate vector at time t is the fixed reference center at the Earth's center. ε is the antenna elevation angle of the target NGSO constellation ground station. The antenna elevation angle constraint of the target NGSO constellation ground station is ε≥ε0, where ε0 represents the minimum antenna elevation angle of the target NGSO constellation ground station.

[0181] Analyze the number of visible satellites and remaining communication time of the target NGSO satellite ground station from the target NGSO constellation;

[0182] Assume the total number of target NGSO constellation satellites is N. c The j-th satellite S j,t The coordinates at time t can be obtained from step S1, then the target NGSO constellation N c The satellites among the 3 satellites that can establish a communication link with the target NGSO ground station are denoted as . but Indicates satisfaction A set of satellites, denoted as N ob1,t Represents the set at time t The number, i.e., the total number of N target NGSO ground stations within the visible range at time t. ob1,t The target NGSO constellation satellites can establish communication links with it, and their coordinates are as follows: j∈[1,N ob1,t ]; Represents vector The model, For the target NGSO constellation ground station G to the target NGSO constellation satellite The distance between them This indicates the target NGSO constellation's j-th satellite. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0183] At time t for N ob1,t There are 10 available communication satellites, and each satellite is calculated to be just below the required standard. The moment The remaining communication time for each available satellite is then...

[0184] Analyze the number of visible satellites at the target NGSO satellite ground station for interference with the NGSO constellation;

[0185] Assume the total number of interfering NGSO constellation satellites is N. inter The j-th interfering satellite S j,t The coordinates at time t can be obtained from step S1, then the interfering NGSO constellation N inter The satellites that interfered with the target NGSO ground station are denoted as but Indicates satisfaction The set of interfering satellites, denoted as N. ob2,t Represents the set at time t The number, i.e., the total number of N target NGSO ground stations within the visible range at time t. ob2,t Several interfering NGSO constellation satellites caused interference, with their coordinates being respectively... i∈[1,N ob2,t ]; NGSO constellation N interf N satellites that interfered with the target NGSO ground station interf To interfere with the total number of NGSO constellation satellites, Represents vector The model, To target the NGSO constellation ground station G to interfere with the NGSO constellation satellites The distance between them This indicates interference with the j-th satellite in the NGSO constellation. The coordinate vector at time t at a fixed reference center in the Earth's center.

[0186] Step S3: Construct a mathematical model of the target NGSO constellation ground station receiving useful signals from the target NGSO satellite and a model of the target NGSO constellation ground station receiving interference signals from interfering NGSO satellites;

[0187] At time t, communication satellites are available within the line of sight of the target NGSO constellation ground station G. The desired signal power of the transmitter at the ground station G receiver is denoted as Considering only free-space link loss, it can be calculated using the following formula;

[0188]

[0189] in: The transmit power of the j-th available satellite in the target NGSO constellation within the line of sight of the ground station G of the target NGSO constellation is expressed in W. This is the radiation pattern of the target NGSO constellation satellite antenna, with units of relative ratios; 0 indicates that the target NGSO constellation ground station G is beam-aligned with the satellite during communication. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; λ j The carrier wavelength of the j-th satellite in the NGSO constellation system within the line of sight of the ground station G of the target NGSO constellation is represented in meters (m). The distance between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t is expressed in meters.

[0190] At time t, the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of the ground station G of the target constellation system at the receiver of the target NGSO ground station G is denoted as I. j,t Under the condition of only considering free space link loss, it can be calculated by the following formula;

[0191]

[0192] Among them, I j,t Let be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of ground station G of the target constellation system at time t, received at the target NGSO ground station G. The equivalent transmit power of the j-th potential interfering satellite of the interfering NGSO constellation within the line of sight of the target NGSO constellation ground station G in the overlapping frequency band is expressed in W. This represents the radiation pattern of the interfering NGSO constellation satellite antennas, expressed as a relative ratio. This represents the antenna radiation pattern of the target NGSO constellation ground station G, expressed as a relative ratio; φ j,tθ represents the angle between the antenna of the j-th interfering NGSO constellation and the target NGSO constellation ground station G at time t. j,t λ represents the angle between the communication link of the target NGSO constellation ground station G and the j-th interfering satellite of the interfering NGSO constellation at time t; j The wavelength of the signal from the j-th satellite interfering with the NGSO constellation, within the line of sight of the target NGSO constellation ground station G, is expressed in meters (m). The distance, in meters (m), represents the distance between the target NGSO constellation ground station G and the j-th satellite of the interfering NGSO constellation within its line of sight at time t.

[0193] O represents the geocenter corresponding to the geocentric fixed coordinate system. D represents the distance, in meters, between the target NGSO constellation ground station G and the j-th target NGSO constellation satellite within its line of sight at time t; j,t This represents the distance between the j-th potential interfering satellite of the NGSO constellation and the Earth's center within the line of sight of ground station G.

[0194] in Interfering with NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame. The target NGSO satellites The coordinate values ​​in the x, y, and z directions in a geocentric fixed reference frame.

[0195] Step S4: Calculate the expected ratio of the useful signal to the sum of the interference signal and noise power of the target NGSO ground station during the remaining communication time of each available communication satellite;

[0196]

[0197] Among them, E(SINR) j,t SINR is the expected value of the target NGSO ground station within the remaining communication time of each available communication satellite. j,t This represents the signal-to-interference-plus-noise ratio (SIR) of the j-th satellite within the line of sight of the target NGSO constellation ground station at time t. t represents the signal power received by the target NGSO constellation ground station from the j-th available satellite within its line of sight at time t; t0 represents the time when the target NGSO constellation ground station establishes a communication link or switches to a communication satellite. This represents the remaining communication time of the j-th available communication satellite in the target NGSO constellation at time t; This represents the sum of interference signal power received by the target NGSO constellation ground station from potential interfering satellites within its line of sight at time t; N represents the equivalent noise of the interfered earth station receiver, which is Gaussian white noise. ob2,t To interfere with the number of visible jamming satellites of the NGSO constellation ground station, I j,t Let t be the interference signal power of the transmitter of the j-th potential interfering satellite within the line of sight of the ground station G of the target constellation system at time t, and the signal power of the receiver at the target NGSO ground station G.

[0198] Step S5: Construct an NGSO constellation interference avoidance model based on maximizing the expected signal-to-interference-plus-noise ratio;

[0199] j t =argmaxE(SINR) j,t )

[0200] j t It is the objective function, which represents the index of the target NGSO constellation satellite that is ultimately selected each time a decision needs to be made.

[0201] The device embodiments described above are merely illustrative. The modules described as separate components may or may not be physically separate, and the components shown as modules may or may not be physical modules; that is, they may be located in one place or distributed across multiple network modules. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0202] Through the detailed description of the above embodiments, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, including read-only memory (ROM), random access memory (RAM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), one-time programmable read-only memory (OTPROM), electrically-Erasable Programmable Read-Only Memory (EEPROM), compact disc read-only memory (CD-ROM) or other optical disc storage, disk storage, magnetic tape storage, or any other computer-readable medium that can be used to carry or store data.

[0203] Finally, it should be noted that the NGSO constellation interference avoidance method disclosed in the embodiments of the present invention is only a preferred embodiment of the present invention and is only used to illustrate the technical solution of the present invention, not to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method of NGSO constellation interferences avoidance, comprising: The method comprises: S1, updating satellite ephemeris of a target NGSO constellation and an interfering NGSO constellation to obtain satellite ephemeris information; The satellite ephemeris information includes a satellite orbit inclination i, a satellite ascending node right ascension Ω0 at a TLE data update time, a satellite orbit eccentricity e, a satellite orbit perigee argument E0 at the TLE data update time, and a satellite average number of revolutions around the earth per day N s ; S2, processing the satellite ephemeris information to obtain NGSO constellation-to-NGSO constellation interference information; The NGSO constellation-to-NGSO constellation interference information comprises a target NGSO constellation ground station visible communication satellite number, a residual visible time of each communication satellite, and an interfering NGSO constellation ground station visible interfering satellite number of the target NGSO constellation; S3, constructing a useful signal model and an interference signal model; The useful signal model is a model of a target NGSO constellation ground station receiving a useful signal from a target NGSO satellite; The interference signal model is a model of a target NGSO constellation ground station receiving an interference signal from an interfering NGSO satellite; S4, using the useful signal model and the interference signal model, calculating the expected value E(SINR j,t ) of the target NGSO constellation ground station within the remaining visible time of each available communication satellite j,t ); E(SINR j,t ) is the expected value of the target NGSO constellation ground station within the remaining visible time of each available communication satellite, SINR j,t represents the signal-to-interference-and-noise ratio of the jth satellite within the visible range of the target NGSO constellation ground station at time t; The expected value is a mathematical expectation of a ratio of the useful signal to a sum of the interference signal and noise power; S5, processing the expected value to obtain an NGSO constellation-to-NGSO constellation interference avoidance model; j t = argmax E (SINR j,t ) j t is the objective function, which represents the index of the final selected target NGSO constellation satellite for each decision-making.

2. The NGSO constellation interference avoidance method of claim 1, wherein, The processing of the satellite ephemeris information to obtain the NGSO constellation-to-NGSO constellation interference information comprises: S21, calculating the number of orbits N of the satellite around the Earth per day using a satellite orbit period calculation model s processing to obtain the satellite orbit period; The satellite orbit operation period calculation model expression is: wherein T s is the satellite orbit period, a is the semi-major axis of the satellite, T e represents the Earth rotation period, takes the value of 86164 s, μ is the Earth gravitational constant, is 3.986005 x 10 14 m 3 / s 2 ; S22, processing the satellite orbit operation period to obtain coordinates of an NGSO satellite in a geocentric fixed coordinate system; S23, processing the coordinates of the NGSO satellite in the geocentric fixed coordinate system to obtain the NGSO constellation-to-NGSO constellation interference information.

3. The NGSO constellation interference avoidance method of claim 2 wherein, The processing of the satellite orbit operation period to obtain the coordinates of the NGSO satellite in the geocentric fixed coordinate system comprises: S221, processing the satellite orbit operation period to obtain an average angular velocity of the NGSO satellite; The average angular velocity of the NGSO satellite is ω; S222, processing the average angular velocity of the NGSO satellite to obtain a latitude parameter of the NGSO satellite at a calculation time t; E0 is a satellite orbit perigee argument at a TLE data update time, and t is a time experienced by the calculation time relative to a TLE data update time; The NGSO satellite calculates a latitude parameter E at time instant t t The expression is: E t = E0+ ω x t S223, processing the satellite orbit inclination i to obtain an NGSO satellite orbit precession rate; The NGSO satellite orbit precession rate expression is: S224, processing the NGSO satellite orbit precession rate to obtain a right ascension of the NGSO satellite at the calculation time t; where Ω r is the orbital precession rate of the NGSO satellite, J2 is the second-order gravitational coefficient, which is 1082.6 x 10 -6 , R e is the radius of the Earth, which is 6378000 m, a is the semi-major axis of the satellite, and μ is the gravitational constant of the Earth; The right ascension of the NGSO satellite at the calculation time t expression is: S225, processing the right ascension of the NGSO satellite at the calculation time t to obtain coordinates of the NGSO satellite in the geocentric fixed coordinate system; Ω t = Ω0+ Ω r × t where Ω t is the right ascension of the ascending node of the NGSO satellite at the time t; The coordinates of the NGSO satellite in the geocentric fixed coordinate system expression is: The processing of the coordinates of the NGSO satellite in the geocentric fixed coordinate system to obtain the NGSO constellation-to-NGSO constellation interference information comprises: wherein, denotes the jth satellite S of the NGSO constellation j coordinate vector of the time instant t in the earth-fixed reference center x j,t , y j,t , z j,t is a coordinate vector of the jth satellite S j at time t in the earth-fixed reference center in the x, y, z directions, E j,t is a latitude parameter of the jth satellite S j at time t, Ω j,t is a right ascension parameter of the jth satellite S j at time t.

4. The NGSO constellation interference avoidance method of claim 2 wherein, S231, processing latitude and longitude information of a target NGSO constellation ground station G to obtain a coordinate vector of the target NGSO constellation ground station in a geocentric fixed reference system; ​ An expression of a coordinate vector of a target NGSO constellation ground station in the earth-centered fixed reference frame is: wherein, is the coordinate vector of the target NGSO constellation ground station in the earth-centered fixed reference frame, the latitude and longitude information of the target NGSO constellation ground station G is (lat, lon), lat represents the latitude of the target NGSO constellation ground station, and lon represents the longitude of the target NGSO constellation ground station, R e is the radius of the earth, t is the time elapsed since the update time of the TLE data, Ω e represents the rotation angular velocity of the earth on the equator, and is 7.29 x 10 -5 rad / s. S232, processing the coordinate vector of the target NGSO constellation ground station in the earth-centered fixed reference frame to obtain an antenna elevation angle of the target NGSO constellation ground station; An expression of the antenna elevation angle of the target NGSO constellation ground station is: wherein denotes the vector of the module, is a vector between a target NGSO constellation ground station G and a NGSO constellation satellite, denotes the jth satellite S j of the target NGSO constellation ground station at time t in the geocentric fixed reference center, and ε is the antenna elevation angle of the target NGSO constellation ground station, the antenna elevation angle constraint of the target NGSO constellation ground station is ε≥ε0, and ε0represents the minimum antenna elevation angle of the target NGSO constellation ground station. S233, processing the coordinate vector of the target NGSO constellation ground station in the earth-centered fixed reference frame and the antenna elevation angle of the target NGSO constellation ground station to obtain NGSO constellation interferences information.

5. The NGSO constellation interference avoidance method of claim 4, wherein, The target NGSO constellation ground station visible communication satellite number of the target NGSO constellation ground station is the set of time t The number N ob1,t ; The set To meet The set of satellites, wherein, The target NGSO constellation N c Satellites of the target NGSO constellation capable of establishing a communication link with a ground station of the target NGSO constellation, N c The total number of satellites of the target NGSO constellation, Denotes the vector The modulus of, The coordinate vector of the jth satellite of the target NGSO constellation At time t in the Earth-Centered Fixed reference frame.

6. The NGSO constellation interference avoidance method of claim 4, wherein, N the number of interfering NGSO constellation visible satellites to the target NGSO constellation ground station ob2,t for time t set the number of interfering NGSO constellation visible satellites to the target NGSO constellation ground station The set To meet The set of interfering satellites, The number of interfering NGSO satellites, N interf The number of satellites in the constellation that interfere with the target NGSO constellation ground stations, N interf The total number of interfering NGSO constellation satellites, Denotes the norm of a vector The norm of a vector 7. The NGSO constellation interference avoidance method of claim 4, wherein, An expression of the useful signal model is: wherein, PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; j PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; PtargetNGSOsatellitej(t) is the transmit power of the jth available satellite of the target NGSO constellation in the target NGSO constellation ground station G visible range at time t, in W; 8. The NGSO constellation interference avoidance method of claim 7, wherein, An expression of the interference signal model is: where I j,t is the interference signal power at the target NGSO constellation ground station G receiver from the jth potential interfering satellite transmitter within the target NGSO constellation ground station G's field of view at time t, is the equivalent transmit power of the jth potential interfering satellite of the interfering NGSO constellation within the coincident frequency band in the target NGSO constellation ground station G's field of view at time t, in W; is the interfering NGSO constellation satellite antenna radiation pattern in relative units; φ j,t is the angle between the jth interfering satellite antenna pointing of the interfering NGSO constellation and the target NGSO constellation ground station G at time t; θ j,t is the angle between the target NGSO constellation ground station G's communication link and the jth interfering satellite of the interfering NGSO constellation at time t; λ j is the signal wavelength of the jth satellite of the interfering NGSO constellation within the target NGSO constellation ground station G's field of view, in m; is the distance between the target NGSO constellation ground station G and the jth satellite of the interfering NGSO constellation within its field of view at time t, in m; O represents the Earth's center corresponding to the Earth-fixed coordinate system; D j,t denotes the distance between the jth potential interfering satellite of the NGSO constellation and the Earth's center within the line of sight range of the ground station G, wherein are interfering NGSO satellites are coordinate values in the x, y, z directions in the Earth-fixed reference frame, are target NGSO satellites are coordinate values in the x, y, z directions in the Earth-fixed reference frame, g,t g,t g,t are coordinate values in the x, y, z directions in the Earth-fixed reference frame of the target NGSO constellation ground station G; is the vector from the target NGSO constellation ground station G to the jth satellite of the target NGSO constellation at time t, is the vector from the target NGSO constellation ground station G to the jth satellite of the interfering NGSO constellation at time t, denotes the coordinate vector of the jth satellite of the interfering NGSO constellation at time t in the Earth-fixed reference center.​​ 9. The NGSO constellation interference avoidance method of claim 8, wherein, An expression of an expected value of the target NGSO constellation ground station in a remaining visible time of each available communication satellite is: Wherein, t0 represents the time when the target NGSO constellation ground station establishes a communication link or switches the communication satellite; represents the residual visible time of the jth available communication satellite of the target NGSO constellation at time t; represents the sum of the interference signal power of the potential interfering satellites in the visible range received by the target NGSO constellation ground station at time t; N represents the equivalent noise of the interfered earth station receiver, which is Gaussian white noise, N ob2,t is the number of visible interfering satellites of the interfering NGSO constellation to the target NGSO constellation ground station.

Images