A satellite internet of things uplink interference analysis method, device and storage medium

Abstract

The present invention discloses a satellite Internet of Things uplink interference analysis method, device and storage medium, the method comprising: modeling using Poisson point distribution ground terminals according to terminal distribution density to obtain a terminal distribution model; calculating the regionalized transmission field strength so that the regionalized uplink aggregate interference is only related to the distribution parameters and the transmission parameters, and has nothing to do with the exact location and number of the terminal; processing the superposition field at the receiving end and analyzing the interference signal. The method has good fitting performance for different numbers of terminals, can effectively reduce the calculation complexity during long-term analysis, and can be used for efficient interference analysis and calculation.

Description

Technical Field

[0001] The present invention relates to a satellite Internet of Things uplink interference analysis method, device and storage medium, and belongs to the technical field of satellite communications. Background Art

[0002] With the development of satellite communication, low-orbit satellite communication has become a hot spot in the field of space-based information systems. At the same time, due to the continuous maturity of micro-satellite mass production technology and the rapid development of launch technologies such as reusable launch vehicles and multiple satellites in one rocket, a new round of low-orbit satellite IoT constellation development boom has been triggered. Foreign technology companies such as OneWeb, SpaceX, and Amazon have launched their own space IoT plans. Domestic Aerospace Science and Technology, Aerospace Science and Industry, Chinese Academy of Sciences, China Electronics Technology Group, and many private enterprises have also announced the development plans of low-orbit satellite constellations, including Hongyan, Xingyun, Tianxiang and other systems. However, the negative effect of the explosive development of satellite constellations is the tension of frequency and orbit resources. Before the satellite system is designed and launched, it is necessary to declare the frequency and orbit according to the relevant ITU regulations and frequency division rules, and to coordinate the necessary frequency and orbit with the existing system to ensure the coexistence and sharing of various systems. Due to the large number of satellites in the low-orbit satellite constellation system and the global coverage of the NGSO constellation, there are multiple interfering satellites within the visible range of users or gateway stations, resulting in more serious aggregate interference, and the probability of collinear extreme interference will also increase. In addition, compared with the GSO system, the NGSO satellite constellation system has the characteristics of a large number of satellites, the constantly changing positions and beam pointing of satellites, and the earth stations often accessing multiple satellites. At the same time, the rapid development of ground wireless communications has led to a sharp increase in the amount of calculation when calculating the interference of NGSO system satellites, and the system has a large amount of calculation. Different from the relatively static interference scenario of the GEO system, the relative position of the NGSO satellite and the ground station is changing, which leads to the time-varying spatial geometric relationship between satellites and between satellites and ground stations, increasing the complexity of interference calculation. There have been many studies on the interference analysis between NGSO constellations. Zhou Aosong et al. studied the interference calculation method between NGSO constellations, but did not provide relevant simulation proofs. Li Rui et al. calculated and simulated the instantaneous interference at different sampling points based on the time-varying characteristics of the NGSO constellation, but only considered the interference between the NGSO constellation and the GSO satellite. For interference between NGSO systems, Jin Jin et al. proposed a method of using the angle between satellite links as an equivalent interference indicator and gave the availability probability between systems on a global scale using the OneWeb system and O3b system as examples. However, due to the differences in transmission power and link loss between different systems, the threshold value of the angle is also different. Therefore, the interference signal cannot be described by a unified angle, so it is only suitable for interference analysis between two systems. Summary of the invention

[0003] Purpose of the invention: In order to overcome the deficiencies in the prior art, the present invention provides a satellite Internet of Things uplink interference analysis method. Starting from the mechanism of co-frequency interference generation, the regionalized transmission field strength is used as the modeling core of the uplink aggregate interference of dense Internet of Things terminals, and the Poisson point distribution model is introduced to perform equivalent modeling of the ground terminal, so that the regionalized uplink aggregate interference is only related to the distribution parameters and the transmission parameters, but has nothing to do with the exact location and number of the terminal. It can reduce the complexity of interference calculation and effectively reduce the analysis time, thereby realizing efficient calculation of regionalized uplink aggregate interference.

[0004] Technical solution: To achieve the above purpose, the technical solution adopted by the present invention is:

[0005] In a first aspect, a satellite Internet of Things uplink interference analysis method is provided, comprising:

[0006] Step 1: According to the terminal distribution density λ b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model;

[0007] Step 2: According to the terminal distribution model obtained in step 1, the uplink regionalized transmission field strength of the dense IoT terminals is calculated and superimposed to obtain the combined field strength E;

[0008] Step 3: According to the combined field strength E obtained in step 2, the interference intensity of the receiving end is calculated using link calculation to obtain the satellite Internet of Things uplink interference analysis result.

[0009] In some embodiments, in step 1, according to the terminal distribution density λ b The model is modeled using Poisson point distribution ground terminals, including:

[0010] 1a) Satellite IoT uplink interference scenario is set as: the terminal transmission power is P t , an omnidirectional antenna is used for transmission, the receiving end uses the same antenna model, and the satellite orbit height is h; in some specific embodiments, the receiving end antenna model adopts ITU-R S.1528-0;

[0011] 1b) Assume that the distribution density of ground terminals with potential interference in a circle with radius R is λ b , the location information is L j (r j ,θ j ), where θ1, θ2, θ3,…,θ j independent and identically distributed in uniform distribution U[0,2π] and r1, r2, r3, ..., r j Independent, j represents the jth terminal, r j represents the distance between the jth terminal and the center O, θ j represents the polar angle of the jth terminal;

[0012] 1c) Assume that the number of terminals in the circle is a random variable M b , then we get M b =j's probability model P:

[0013]

[0014] In some embodiments, in step 2, calculating the uplink regionalized transmission field strength of dense IoT terminals includes:

[0015] 2a) Take n points in the area above the terminal distribution, denoted as S, then S = {S1, S2, S3, ..., S n}, define r i→j For any point S in the space i To ground terminal L j The horizontal distance of the link is calculated as:

[0016]

[0017] where d i→j Point S i The spatial distance to the jth terminal, h is the satellite orbit height;

[0018] 2b) Calculate r i→j The probability density function f(r i→j ):

[0019]

[0020] For each r i→j , whose angles all obey the uniform distribution of [0,2π];

[0021] 2c) Calculate the jth terminal at point S i The field strength E j→i :

[0022]

[0023] Where P j is the transmission power of the jth terminal, G t is the terminal transmitting antenna gain. In this scenario, all terminals have the same antenna model, which is omnidirectional transmitting antenna.

[0024] In some embodiments, in step 3, calculating the interference strength of the receiving end includes:

[0025] 3a) In the scenario where the uplink of the low-orbit satellite system is interfered, calculate the interference source signal power received by the receiver at time t:

[0026]

[0027] Where I(t) is the interference signal power of a single beam received at time t, P j is the jth terminal transmission power, λ j is the signal frequency of the jth terminal, G t (θ j (t)) is the antenna gain of the interference signal transmitting antenna at a certain angle away from its main axis, is the antenna gain of the receiving antenna of the interfered system at a certain angle away from its main axis, θ j (t) and They represent the angle between the main axis of the transmitting antenna of the jth interference terminal and the interference link and the angle between the main axis of the receiving antenna and the interference link at time t, respectively. j (t) represents the link distance;

[0028] 3b) Calculate the aggregate interference power I received by the satellite system using the interference power caused by the interference source within the communication angle range total (t):

[0029]

[0030] Where N i (t) represents the number of terminals within the beam coverage at the current moment, and N1 represents the interference beam;

[0031] 3c) Use the equivalent noise temperature to construct the system noise power:

[0032] N=K·T·W

[0033] Where K is the Boltzmann constant, T is the system equivalent noise temperature, and W is the bandwidth;

[0034] 3d) Calculate the interference-to-noise ratio:

[0035]

[0036] Where E j Represents the field strength generated by the jth terminal at time t.

[0037] In a second aspect, a satellite Internet of Things uplink interference analysis device is provided, comprising:

[0038] The terminal modeling module is configured to b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model;

[0039] The interference field strength calculation module is configured to calculate and superimpose the uplink regionalized transmission field strength of the dense IoT terminals according to the obtained terminal distribution model to obtain the combined field strength E;

[0040] The interference analysis module is configured to calculate the interference intensity of the receiving end according to the obtained combined field strength E by using link calculation to obtain the satellite Internet of Things uplink interference analysis result.

[0041] In a third aspect, a satellite Internet of Things uplink interference analysis device is provided, comprising: a processor and a storage medium;

[0042] The storage medium is used to store instructions;

[0043] The processor is used to operate according to the instructions to execute the steps of the method according to the first aspect.

[0044] According to a fourth aspect, a storage medium is provided, on which a computer program is stored, characterized in that when the computer program is executed by a processor, the steps of the method according to the first aspect are implemented.

[0045] Beneficial effect: The present invention provides a satellite Internet of Things uplink interference analysis method, which takes the regionalized transmission field strength as the modeling core of the dense Internet of Things terminal uplink aggregate interference, and introduces the Poisson point distribution model to perform equivalent modeling of the ground terminal, so that the regionalized uplink aggregate interference is only related to the distribution parameters and transmission parameters, but has nothing to do with the exact location and number of the terminal. It can reduce the complexity of interference calculation and effectively reduce the analysis time. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 It is a flowchart of the implementation process of the method of the embodiment of the present invention;

[0047] Figure 2 This is a schematic diagram of field intensity superposition in an embodiment of the present invention;

[0048] Figure 3 The absolute error of the interference-to-noise ratio achieved by the method of the embodiment of the present invention;

[0049] Figure 4 The relative error of the interference-to-noise ratio achieved by the method of the embodiment of the present invention. DETAILED DESCRIPTION

[0050] The present invention will be further described below in conjunction with the accompanying drawings and embodiments. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and cannot be used to limit the protection scope of the present invention.

[0051] In the description of the present invention, "several" means more than one, "many" means more than two, "greater than", "less than", "exceed", etc. are understood to exclude the number itself, and "above", "below", "within", etc. are understood to include the number itself. If there is a description of "first" or "second", it is only used for the purpose of distinguishing the technical features, and cannot be understood as indicating or implying the relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features.

[0052] In the description of the present invention, the description with reference to the terms "one embodiment", "some embodiments", "illustrative embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0053] Example 1

[0054] A satellite Internet of Things uplink interference analysis method, comprising:

[0055] Step 1: According to the terminal distribution density λ b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model;

[0056] Step 2: According to the terminal distribution model obtained in step 1, the uplink regionalized transmission field strength of the dense IoT terminals is calculated and superimposed to obtain the combined field strength E;

[0057] Step 3: According to the combined field strength E obtained in step 2, the interference intensity of the receiving end is calculated using link calculation to obtain the satellite Internet of Things uplink interference analysis result.

[0058] See also Figure 1 The specific implementation steps of the present invention are as follows:

[0059] Step 1: Terminal modeling

[0060] 1a) Satellite IoT uplink interference scenario is set as: the terminal transmission power is P t , using an omnidirectional antenna for transmission, the receiving antenna model adopts ITU-R S.1528-0, and the satellite orbit height is h;

[0061] 1b) Assume that the distribution density of ground terminals with potential interference in a circle with radius R is λ b , the location information is L j (r j ,θ j ), where θ1, θ2, θ3,…,θj independent and identically distributed in uniform distribution U[0,2π] and r1, r2, r3, ..., r j Independent, j represents the jth terminal, r j represents the distance between the jth terminal and the center O, θ j represents the polar angle of the jth terminal;

[0062] 1c) Assume that the number of terminals in the circle is a random variable M b , then we get M b =j's probability model P:

[0063]

[0064] Step 2: Calculate the interference field strength. Figure 2 This is a schematic diagram of field intensity superposition in an embodiment of the present invention;

[0065] 2a) Take n points in the area above the terminal distribution, denoted as S, then S = {S1, S2, S3, ..., S n}, define r i→j For any point S in the space i To ground terminal L j The horizontal distance of the link is calculated as:

[0066]

[0067] where d i→j Point S i The spatial distance to the jth terminal, h is the satellite orbit height;

[0068] 2b) Calculate r i→j The probability density function f(r i→j ):

[0069]

[0070] For each r i→j , whose angles all obey the uniform distribution of [0,2π];

[0071] 2c) Calculate the jth terminal at point S i The field strength E j→i :

[0072]

[0073] Where P j is the transmission power of the jth terminal, G t is the terminal transmitting antenna gain. In this scenario, all terminals have the same antenna model, which is omnidirectional transmitting antenna.

[0074] Step 3: Interference analysis

[0075] 3a) In the scenario where the uplink of the low-orbit satellite system is interfered, calculate the interference source signal power received by the receiver at time t:

[0076]

[0077] Where I(t) is the interference signal power of a single beam received at time t, P j is the jth terminal transmission power, λ j is the signal frequency of the jth terminal, G t (θ j (t)) is the antenna gain of the interference signal transmitting antenna at a certain angle away from its main axis, is the antenna gain of the receiving antenna of the interfered system at a certain angle away from its main axis, θ j (t) and They represent the angle between the main axis of the transmitting antenna of the jth interference terminal and the interference link and the angle between the main axis of the receiving antenna and the interference link at time t, respectively. j (t) represents the link distance;

[0078] 3b) Calculate the aggregate interference power I received by the satellite system using the interference power caused by the interference source within the communication angle range total (t):

[0079]

[0080] Where N i (t) represents the number of terminals within the beam coverage at the current moment, and N1 represents the interference beam;

[0081] 3c) Use the equivalent noise temperature to construct the system noise power:

[0082] N=K·T·W

[0083] Where K is the Boltzmann constant, T is the system equivalent noise temperature, and W is the bandwidth;

[0084] 3d) Calculate the interference-to-noise ratio:

[0085]

[0086] Where E j Represents the field strength generated by the jth terminal at time t.

[0087] Figure 3 The absolute error of the interference-to-noise ratio achieved by the method of the embodiment of the present invention;

[0088] Figure 4 The relative error of the interference-to-noise ratio achieved by the method of the embodiment of the present invention.

[0089] Example 2

[0090] In a second aspect, this embodiment provides a satellite Internet of Things uplink interference analysis device, including:

[0091] The terminal modeling module is configured to b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model;

[0092] The interference field strength calculation module is configured to calculate and superimpose the uplink regionalized transmission field strength of the dense IoT terminals according to the obtained terminal distribution model to obtain the combined field strength E;

[0093] The interference analysis module is configured to calculate the interference intensity of the receiving end according to the obtained combined field strength E by using link calculation to obtain the satellite Internet of Things uplink interference analysis result.

[0094] Example 3

[0095] In a third aspect, this embodiment provides a satellite Internet of Things uplink interference analysis device, including a processor and a storage medium;

[0096] The storage medium is used to store instructions;

[0097] The processor is used to operate according to the instructions to execute the steps of the method according to embodiment 1.

[0098] Example 4

[0099] In a fourth aspect, this embodiment provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of the method described in Example 1 are implemented.

[0100] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0101] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0102] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0103] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0104] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A satellite Internet of Things uplink interference analysis method, characterized in that: include: Step 1: According to the terminal distribution density λ b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model; Step 2: According to the terminal distribution model obtained in step 1, the uplink regionalized transmission field strength of the dense IoT terminals is calculated and superimposed to obtain the combined field strength E; Step 3: According to the combined field strength E obtained in step 2, the interference intensity of the receiving end is calculated by link calculation to obtain the satellite Internet of Things uplink interference analysis result, which specifically includes: 3a) In the scenario where the uplink of the low-orbit satellite system is interfered, calculate the interference source signal power received by the receiver at time t: Where I(t) is the interference signal power of a single beam received at time t, P j is the jth terminal transmission power, λ j is the signal frequency of the jth terminal, G t (θ j (t)) is the antenna gain of the interference signal transmitting antenna at a certain angle away from its main axis, is the antenna gain of the receiving antenna of the interfered system at a certain angle away from its main axis, θ j (t) and They represent the angle between the main axis of the transmitting antenna of the jth interference terminal and the interference link and the angle between the main axis of the receiving antenna and the interference link at time t, respectively. j (t) represents the link distance; 3b) Calculate the aggregate interference power I received by the satellite system using the interference power caused by the interference source within the communication angle range total (t): Where N i (t) represents the number of terminals within the beam coverage at the current moment, and N1 represents the interference beam; 3c) Use the equivalent noise temperature to construct the system noise power: N=K·T·W Where K is the Boltzmann constant, T is the system equivalent noise temperature, and W is the bandwidth; 3d) Calculate the interference-to-noise ratio: Where E j Represents the field strength generated by the jth terminal at time t.

2. The satellite Internet of Things uplink interference analysis method according to claim 1, characterized in that: According to the terminal distribution density λ b The model is modeled using Poisson point distribution ground terminals, including: 1a) Satellite IoT uplink interference scenario is set as: the terminal transmission power is P t , using an omnidirectional antenna for transmission, the receiving end uses the same antenna model, and the satellite orbit height is h; 1b) Assume that the distribution density of ground terminals with potential interference in a circle with radius R is λ b , the location information is L j (r j ,θ j ), where θ1, θ2, θ3,…,θ j independent and identically distributed in uniform distribution U[0,2π] and r1, r2, r3, ..., r j Independent, j represents the jth terminal, r j represents the distance between the jth terminal and the center O, θ j represents the polar angle of the jth terminal; 1c) Assume that the number of terminals in the circle is a random variable M b , then we get M b =j's probability model P:

3. The satellite Internet of Things uplink interference analysis method according to claim 1, characterized in that: Calculate the regionalized uplink transmission field strength of dense IoT terminals, including: Calculate the jth terminal at point S i The field strength E j→i : Among them, take n points in the area above the terminal distribution, denoted as S, then S={S1,S2,S3,......,S n }, define r i→j For any point S in the space i To ground terminal L j The horizontal distance, d i→j Point S i The spatial distance to the jth terminal, f(r i→j ) is r i→j The probability density function, P j is the transmission power of the jth terminal, G t is the terminal transmitting antenna gain.

4. The satellite Internet of Things uplink interference analysis method according to claim 3 is characterized in that: d i→j The calculation methods include: h is the satellite orbit altitude.

5. The satellite Internet of Things uplink interference analysis method according to claim 3 is characterized in that: f(r i→j ) include: For each r i→j , and their angles all obey the uniform distribution of [0,2π].

6. The satellite Internet of Things uplink interference analysis method according to claim 2 is characterized in that: The receiving end antenna model adopts ITU-R S.1528-0.

7. A satellite Internet of Things uplink interference analysis device, characterized in that: include: The terminal modeling module is configured to b The Poisson point distribution ground terminal is used for modeling to obtain the terminal distribution model; The interference field strength calculation module is configured to calculate and superimpose the uplink regionalized transmission field strength of the dense IoT terminals according to the obtained terminal distribution model to obtain the combined field strength E; The interference analysis module is configured to calculate the interference intensity of the receiving end according to the obtained combined field strength E by using link calculation to obtain the satellite Internet of Things uplink interference analysis result, which is specifically used for: 3a) In the scenario where the uplink of the low-orbit satellite system is interfered, calculate the interference source signal power received by the receiver at time t: Where I(t) is the interference signal power of a single beam received at time t, P j is the jth terminal transmission power, λ j is the signal frequency of the jth terminal, G t (θ j (t)) is the antenna gain of the interference signal transmitting antenna at a certain angle away from its main axis, is the antenna gain of the receiving antenna of the interfered system at a certain angle away from its main axis, θ j (t) and They represent the angle between the main axis of the transmitting antenna of the jth interference terminal and the interference link and the angle between the main axis of the receiving antenna and the interference link at time t, respectively. j (t) represents the link distance; 3b) Calculate the aggregate interference power I received by the satellite system using the interference power caused by the interference source within the communication angle range total (t): Where Ni(t) represents the number of terminals within the beam coverage at the current moment, and N1 represents the interfering beam; 3c) Use the equivalent noise temperature to construct the system noise power: N=K·T·W Where K is the Boltzmann constant, T is the system equivalent noise temperature, and W is the bandwidth; 3d) Calculate the interference-to-noise ratio: Where E j Represents the field strength generated by the jth terminal at time t.

8. A satellite Internet of Things uplink interference analysis device, characterized in that: include: including processor and storage medium; The storage medium is used to store instructions; The processor is configured to operate according to the instructions to execute the steps of the method according to any one of claims 1 to 6.

9. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images