GEO satellite frequency orbit resource availability evaluation method and system and storage medium

Abstract

The invention discloses a GEO satellite frequency orbit resource availability evaluation method and system and a storage medium, relates to the technical field of spatial frequency orbit resource management, and aims to solve the problem that a comprehensive and effective GEO satellite frequency orbit resource availability evaluation method is lacked in the prior art. According to the method, influence factors such as the satellite orbit, the declaration frequency, the occupied bandwidth, the link direction, the polarization type and the coverage area are comprehensively considered, the availability sequence of a use scheme is evaluated and given, a standby scheme is given under the condition that the availability is low, the use efficiency of GEO satellite frequency orbit resources can be greatly improved, and the current situation that the GEO satellite frequency orbit resources are insufficient is effectively relieved. By means of the effectiveness evaluation result of the GEO satellite frequency orbit resources, on one hand, a frequency orbit resource use scheme can be rapidly formed, labor cost and time cost are saved, on the other hand, the current GEO satellite frequency orbit resource use situation can be comprehensively reflected, and auxiliary decision support is provided for satellite network data declaration and utilization.

Description

Technical Field

[0001] The present invention relates to the technical field of space frequency and orbital resource management, and more particularly, to a method, system, and storage medium for evaluating the availability of GEO satellite frequency and orbital resources. Background Art

[0002] Satellite frequency and orbital position resources, especially GEO satellite frequency and orbital position resources, are limited resources. The orbital position resources of GEO satellites mainly refer to the 360° orbital plane at 36,000 kilometers above the equator. Here, the satellites are relatively stationary with respect to the Earth. On the premise of limited orbital position resources, the more GEO satellites are launched, the more tense the orbital position resources become, and the greater the difficulty of orbital position declaration, coordination, and use. The frequency resources of GEO satellites mainly refer to the frequency bands used by various services of GEO satellites, including UHF / VHF, L, S, C, X, Ku, Ka, EHF, etc. With the increase in the number of GEO satellites, frequency interference and coordination become more and more difficult. At present, there is no clear evaluation standard for the availability of GEO satellite frequency and orbital resources. The evaluation of the availability of GEO satellite frequency and orbital resources only relies on manual database checking and subjective judgment based on the existing situation of in-orbit satellites and satellite network data, resulting in many review problems when submitting satellite network data, and it is easy to miss effective resources, reducing resource utilization. Therefore, evaluating the rational use of GEO satellite frequency and orbital resources can better guide the application efficiency of GEO satellites. Summary of the Invention

[0003] The technical problem to be solved by the present invention is:

[0004] In the prior art, there is a lack of a comprehensive and effective method for evaluating the availability of GEO satellite frequency and orbital resources.

[0005] The technical solution adopted by the present invention to solve the above technical problem:

[0006] The present invention provides a method for evaluating the availability of GEO satellite frequency and orbital resources, including the following steps:

[0007] Step 1, Evaluation task interaction

[0008] Determine the evaluation requirement indicators, including: preliminary selected orbital position, allowable arc segment range, frequency band type, allowable frequency band range, occupied bandwidth, link direction, polarization type, coverage area;

[0009] Step 2, Evaluation index assignment

[0010] Determine the evaluation indicators, assign values to the evaluation indicators, and perform quantitative calculations on the evaluation indicators;

[0011] Step 3, Evaluation index weight assignment

[0012] Step 4, Input of preliminary selected frequency and orbital plan

[0013] Select the initial frequency and orbit plan through the evaluation task interaction input;

[0014] Step 5, Frequency overlap detection

[0015] Compare the allowable frequency band range of the input initial frequency and orbit with the satellite frequency bands in the basic database. If there is a frequency overlap area, proceed to the next step of calculating the orbit position interval; if there is no frequency overlap area, determine whether all satellites in the basic database have been traversed;

[0016] Step 6, Orbit position interval calculation

[0017] For different frequency bands, design a high-risk interval vector for the radian; according to the initial orbit positions of the input initial frequency and orbit and the satellite orbit positions read from the database, calculate the orbit position interval between satellites. If the orbit position is within the high-risk interval, conduct a coverage area overlap evaluation; if the orbit position is not within the high-risk interval, determine whether all satellites in the basic database have been traversed;

[0018] Step 7, Coverage area overlap analysis

[0019] Compare the coverage area of the input initial frequency and orbit with the coverage areas of the satellites in the database. If the coverage areas overlap, conduct an interference evaluation; if the coverage areas do not overlap, determine whether all satellites in the basic database have been traversed;

[0020] Step 8, Output of evaluation data and results.

[0021] Furthermore, in step 2, trapezoidal fuzzy numbers are used to quantitatively calculate the evaluation indicators, specifically:

[0022]

[0023] Among them, A i (x) represents the membership function, and α and b represent the satisfaction membership function parameters.

[0024] Furthermore, in step 3, the evaluation index weights are determined based on the second-order coefficient of variation, that is:

[0025]

[0026] Among them, v represents the coefficient of variation.

[0027] Furthermore, it also includes generating alternative plans. The alternative plans have different orbit positions from the initial frequency and orbit plan, and the other parameters remain unchanged; the generation of the alternative plans is based on the following method:

[0028] a) Select the frequency and orbit whose orbital parameters do not conflict with the initial frequency and orbit plan;

[0029] b) Select the frequency track with the least difficulty in conflict selection coordination.

[0030] The present invention provides a GEO satellite frequency orbit resource availability evaluation system, which has program modules corresponding to the steps of the method described in any one of the above technical solutions, and executes the steps in the above-mentioned GEO satellite frequency orbit resource availability evaluation method when running.

[0031] The present invention provides a computer-readable storage medium, which stores a computer program configured to implement the steps in the GEO satellite frequency orbit resource availability evaluation method described in any one of the above technical solutions when called by a processor.

[0032] Compared with the prior art, the beneficial effects of the present invention are:

[0033] The present invention provides a GEO satellite frequency orbit resource availability evaluation method, system and storage medium, which can, for the frequency orbit usage plan of GEO satellites, based on the global GEO satellite network data declaration data and the global GEO in-orbit satellite data, comprehensively consider influencing factors such as satellite orbit, declared frequency, occupied bandwidth, link direction, polarization type, coverage area, etc., evaluate and give the availability ranking of the usage plan, and give an alternative plan in the case of low availability, which can greatly improve the usage efficiency of GEO satellite frequency orbit resources and effectively alleviate the current situation of tense GEO satellite frequency orbit resources. Specifically manifested in:

[0034] First, the evaluation index system is perfect, comprehensively considering the global in-orbit satellites, satellite network data declaration data, and the advantageous orbit resources mastered by traditional satellite operators, the index system is more complete, and the evaluation results are more accurate;

[0035] Second, there are various calculation methods for evaluation indexes, which can adapt to the evaluation of different frequency orbit schemes;

[0036] Third, it can quickly form a frequency orbit resource usage plan, saving labor costs and time costs;

[0037] Fourth, it can comprehensively reflect the current situation of GEO satellite frequency orbit resource usage, providing auxiliary decision-making support for satellite network data declaration and utilization. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 It is a schematic diagram of the evaluation task interaction content in the embodiment of the present invention;

[0039] Figure 2 It is a diagram of the selected evaluation data and the generated alternative plan in the embodiment of the present invention;

[0040] Figure 3 It is a flowchart of the evaluation in the embodiment of the present invention;

[0041] Figure 4 Schematic diagram of the evaluation index system in the embodiments of the present invention;

[0042] Figure 5 Schematic diagram of the on-orbit satellite interference risk degree in the embodiments of the present invention;

[0043] Figure 6 Schematic diagram of the risk caused by the resources already in use in the embodiments of the present invention;

[0044] Figure 7 Schematic diagram of the risk caused by the coordination materials in the embodiments of the present invention. Detailed implementation manners

[0045] In order to enable those skilled in the art to better understand the solution of the present invention, the exemplary embodiments or examples of the present invention will be described below in conjunction with the accompanying drawings. Obviously, the described embodiments or examples are only a part of the embodiments or examples of the present invention, rather than all of them. Based on the embodiments or examples in the present invention, all other embodiments or examples obtained by those of ordinary skill in the art without creative efforts shall fall within the scope of protection of the present invention.

[0046] To make the above objects, features, and advantages of the present invention more obvious and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0047] The present invention provides a method for evaluating the availability of GEO satellite frequency and orbit resources, as Figure 3 shown, including the following steps:

[0048] Step 1, Evaluation task interaction

[0049] As Figure 1 shown, determine the relevant requirement indicators according to the requirements of the evaluation party: frequency band type, occupied bandwidth, link direction, polarization type, and coverage area as the input parameters of the primary selected frequency and orbit scheme. At the same time, calculate the arc segment required to cover this area according to the coverage area, that is, the allowable arc segment range. Calculate the available orbit positions according to the allowable arc segment range and the orbit position data in the basic database, and further select the primary selected orbit positions. Further determine the allowable frequency band range of this frequency and orbit scheme according to the link direction, occupied bandwidth, and frequency band type. Determine the interference limit of this system according to the link direction, occupied bandwidth, frequency band type, polarization mode, and in combination with the "Radio Regulations".

[0050] Step 2, Evaluation index assignment customization configuration

[0051] As Figure 4As shown in the figure, the evaluation indicators include: the actual number of satellites within the proximity range, the minimum actual satellite spacing, the worst link interference of adjacent satellites, the worst beam interference of adjacent satellites, the worst link interference of satellites with the minimum spacing, the worst beam interference of satellites with the minimum spacing, the worst link interference of satellites, the worst beam interference of satellites, the number of occupied satellite networks within the proximity range, the minimum spacing of occupied satellite networks, the worst link interference of occupied satellite networks, the worst beam interference of occupied satellite networks, the worst link interference of occupied satellite networks with the minimum spacing, the worst beam interference of occupied satellite networks with the minimum spacing, the worst link interference of occupied satellite networks, the worst beam interference of occupied satellite networks, the number of coordination materials of adjacent status-priority satellite networks, the minimum spacing of status-priority satellite network coordination materials, the worst link interference of adjacent status-priority satellite network coordination materials, the worst beam interference of adjacent status-priority satellite network coordination materials, the worst link interference of coordination materials with the minimum spacing of priority coordination, the worst beam interference of coordination materials with the minimum spacing of priority coordination, the worst link interference of status-priority satellite network coordination materials, and the worst beam interference of status-priority satellite network coordination materials. Non-boolean indicators are quantified through the calculation of fuzzy vectors. Fuzzy vectors process fuzzy evaluation indicators through precise digital means, and can complete a more scientific, reasonable, and practical quantitative evaluation of fuzzy indicators; the evaluation result is a vector, rather than a point value, containing relatively rich information, which can accurately depict the evaluated object and further process it to obtain reference information. The specific method is as follows:

[0052] Let the fuzzy satisfaction level be N, and use X 1 , X 2 , …, X N , X i = (i - 1) / (N - 1), where i = 1, 2, …, N respectively represent the satisfaction level from low to high, X 1 = 0 represents the lowest satisfaction level, and X N = 1 represents the highest satisfaction level. Domain experts give the fuzzy membership functions corresponding to each satisfaction level according to the data change range of single indicators. For example, the satellite revisit period is an attribute indicator where the smaller the better, and trapezoidal fuzzy numbers are used for quantitative calculation. The membership function is shown in the following formula:

[0053]

[0054] Among them, A i (x) represents the membership function, and a and b represent the parameters of the satisfaction membership function, which are given specific values by experts.

[0055] Similarly, for the type where the larger the better and the moderate index type, the same method is used to give the satisfaction levels and their membership functions.

[0056] For each index, after determining the satisfaction level and its membership function, substitute the index data into the above formula to calculate the fuzzy vector of the index. The corresponding parameter values of the membership function are given by experts.

[0057] Taking the worst-link interference of neighboring satellites as an example, assume there are N = 5 satisfaction levels, which are (0, 0.25, 0.5, 0.75, 1) respectively. The parameters of its satisfaction membership function are given by experts, which are a 1 = 10, a 2 = 21, a 3 = 40, a 4 = 80, b 2 = 15, b 3 = 36, b 4 = 60, b 5 = 100. Currently, the worst-link interference of neighboring satellites is 25 dB. Substitute it into the formula to calculate the fuzzy satisfaction membership degree of the revisit period. Its fuzzy vector is {0, 0.2, 0.8, 0, 0}, and its satisfaction degree is [10.75 0.5 0.25 0][0 0.2 0.8 0 0] = 0.55.

[0058] Step 3: Index weight customization and configuration

[0059] Use the Delphi subjective weighting method to achieve the subjective assignment of parameter weight values, use the coefficient of variation objective weighting method to achieve the objective assignment of parameter weight values, and be able to re-customize the index weights according to the needs of the evaluation party or select empirical learning values.

[0060] In order to eliminate the influence of different dimensions of each evaluation index, the coefficient of variation of each index is used to measure the degree of difference between indexes. The coefficient of variation objective weighting method has unique advantages and realizes the determination of multi-index weight factors with different dimensions.

[0061] The determination method of the conventional coefficient of variation is as follows. The coefficient of variation of each index is determined by the following formula.

[0062]

[0063] In the formula,

[0064] v i —— represents the coefficient of variation of the i-th index

[0065] σ i - represents the standard deviation of the i-th index

[0066] —— represents the arithmetic mean of the i-th index.

[0067] The weight w i of each index is defined as:

[0068]

[0069] The steps of using the coefficient of variation objective weighting method to solve the frequency-orbit demonstration problem are relatively simple. The most important part in the solution process is the calculation of the data statistics of the evaluation indicators. Specifically, using the coefficient of variation method to solve the index weights is divided into three steps in total:

[0070] 1) Extract the index data from the input frequency-orbit scheme and the basic database, and conduct statistical analysis on the evaluation indicators. According to the information contained in each indicator itself, solve the two statistics of the arithmetic mean and standard deviation of this indicator.

[0071] 2) According to the definition of the coefficient of variation, directly calculate the coefficient of variation of each level of indicators using the arithmetic mean and standard deviation. The process of solving the coefficient of variation is a dimensionless process. Therefore, the calculated coefficient of variation of the evaluation indicators does not contain the dimensional factors of the original indicator data.

[0072] 3) Sum up the coefficients of variation of the indicators within each level of the indicator layer, and respectively calculate the proportion of the coefficients of variation of different indicators within their respective subsystems, so as to obtain the proportion of each indicator within the corresponding evaluation system.

[0073] The above calculation process is the process of determining each indicator by the conventional coefficient of variation weighting method. However, there are many problems with the traditional coefficient of variation method. For example, the accuracy of this method is limited by the value range of the mean (especially when the mean is very small and tends to zero, the denominator of the formula approaches zero, and the formula becomes inaccurate), and the method itself lacks a reasonable physical explanation and meaning. On the basis of the above method, the second-order coefficient of variation method is used to solve the above problems. See the following formula:

[0074]

[0075] The above formula gives the second-order coefficient of variation of the j-th indicator. The second-order coefficient of variation takes into account the difference between the average value of the squares of each indicator and the square of its average value. The requirement for the value range of the mean by the second-order coefficient of variation is lower than that of the coefficient of variation in the conventional variation parameter weighting method, and the accuracy becomes higher.

[0076] On this basis, the weight value determined by the second-order coefficient of variation is expressed by the following formula, that is

[0077]

[0078] The method of determining the index weight value by the second-order coefficient of variation method is: calculate the average value of the second-order coefficient of variation corresponding to each indicator, and obtain the weight of the second-order coefficient of variation of each indicator accounting for the second-order coefficient of variation of all indicators.

[0079] Step 4. Input of the preliminary selected frequency-orbit scheme

[0080] The evaluation party inputs the preliminary selected frequency and orbit plan through the evaluation task interaction. The parameters of the input frequency and orbit plan mainly include: preliminary selected orbit position, allowable arc range, frequency band type, allowable frequency band range, occupied bandwidth, link direction, polarization type, and coverage area.

[0081] Step 5, Frequency overlap detection

[0082] Compare the allowable frequency band range input by the evaluation party with the satellite frequency bands in the basic database. If there is a frequency overlap area, proceed to the next step of calculating the orbit position interval; if there is no frequency overlap area, determine whether all satellites in the basic database have been traversed.

[0083] Step 6, Orbit position interval calculation

[0084] For different frequency bands, design a high-risk interval vector for the arc. For example, the tentative orbit position interval configuration is 2 degrees for the Ka band, 3 degrees for the Ku band, 2 degrees for the X band, and 4 degrees for the C band. Subsequently, further improve this interval through spectrum compatibility data information and optimization methods.

[0085] According to the preliminary selected orbit position input by the evaluation party and the satellite orbit positions read from the database, calculate the orbit position interval between satellites. If the orbit position is within the high-risk interval, conduct an analysis of coverage area overlap; if the orbit position is not within the high-risk interval, determine whether all satellites in the basic database have been traversed.

[0086] Step 7, Coverage area overlap analysis

[0087] Compare the coverage area in the preliminary selected frequency and orbit plan input by the evaluation party with the coverage area of the satellites in the database. If the coverage areas overlap, conduct an interference analysis; if the coverage areas do not overlap, determine whether all satellites in the basic database have been traversed.

[0088] Step 8, Generation of alternative plans

[0089] As Figure 2 shown, in addition to the frequency and orbit plan input by the evaluation party, according to a certain step of the orbit position, within the arc segment that meets the requirements of the evaluation party, select other orbital positions as alternative plans. During the process of selecting other orbital positions, the system provides a default step of the orbit position and supports the evaluation party to customize it. The default step of the orbit position is 0.2 degrees, and the 360-degree geostationary orbit arc segment in the equatorial plane corresponds to 1800 optional orbit positions. The alternative plans only differ from the plan input by the evaluation party in terms of the orbit position, and the other parameters remain unchanged. It mainly includes two types of recommended orbit positions:

[0090] a) Select orbits whose orbital parameters have no conflict with ours;

[0091] b) Select orbits with conflicts and the least coordination difficulty.

[0092] Step 9, Evaluation of Data and Result Output

[0093] As Figures 5 to 7 shown, in the evaluation process of the evaluation index system of the present invention, relevant weights and input values are recorded, and relevant parameter values can be displayed in the form of a mind map. For all the recommended frequency-orbit schemes evaluated, the details of the scheme scores can be viewed, the risk degree can be traced back, and the specific index items and basic data that cause the evaluation value to be too low can be located, helping the evaluator to master the current satellite or satellite network environment.

[0094] The evaluation data and results output by the present invention can meet the requirements of the evaluator to visualize the evaluation process data and the final results throughout the process. The satellite data and satellite network basic data currently participating in the evaluation are displayed in the form of a list; and the data participating in this evaluation can be exported according to the needs of the evaluator.

[0095] A method (algorithm) for evaluating the availability of GEO satellite frequency-orbit resources proposed by the present invention is the underlying technical core of the present invention, and various products can be derived based on the algorithm.

[0096] Based on the method proposed by the present invention, a system for evaluating the availability of GEO satellite frequency-orbit resources is developed using a programming language. The system has program modules corresponding to the steps of the above technical solution, and executes the steps in the above method for evaluating the availability of GEO satellite frequency-orbit resources when running.

[0097] The computer program of the developed system (software) is stored on a computer-readable storage medium, and the computer program is configured to implement the steps of the above method for evaluating the availability of GEO satellite frequency-orbit resources when called by a processor. That is, the present invention is materialized on a carrier to become a computer program product.

[0098] The various embodiments of the systems and technologies described herein can be implemented in digital electronic circuit systems, integrated circuit systems, dedicated ASICs (application-specific integrated circuits), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: implemented in one or more computer programs, the one or more computer programs can be executed and / or interpreted on a programmable system including at least one programmable processor, the programmable processor can be a dedicated or general programmable processor, and can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit the data and instructions to the storage system, the at least one input device, and the at least one output device.

[0099] The computing procedures (also referred to as programs, software, software applications, or code) in the present invention include machine instructions for a programmable processor, and these computing procedures can be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. As used herein, the terms "machine-readable medium" and "computer-readable medium" refer to any computer program product, device, and / or apparatus (e.g., magnetic disks, optical disks, memories, programmable logic device PLD) for providing machine instructions and / or data to a programmable processor, including a machine-readable medium that receives machine instructions as a machine-readable signal. The term "machine-readable signal" refers to any signal for providing machine instructions and / or data to a programmable processor.

[0100] Although the present invention is disclosed as above, the scope of protection of the present invention is not limited thereto. Those skilled in the art of the present invention can make various changes and modifications without departing from the spirit and scope of the present disclosure, and these changes and modifications will all fall within the scope of protection of the present invention.

Claims

1. A method for assessing the availability of GEO satellite frequency orbital resources, characterized in that: The steps include: Step 1: Evaluate task interactions Determine the evaluation requirements, including: initial orbital position, allowable arc range, frequency band type, allowable frequency band range, occupied bandwidth, link direction, polarization type, and coverage area; Step 2: Assign evaluation indicators Determine the evaluation indicators, assign values ​​to the evaluation indicators, and perform quantitative calculations on the evaluation indicators; Step 3: Assign weights to evaluation indicators Step 4: Input the preliminary frequency and track plan Interactively input the preliminary frequency and orbit plan through the evaluation task; Step 5: Frequency overlap detection Compare the allowable frequency band range of the input primary frequency orbit with the satellite frequency band in the basic database. If there is a frequency overlap area, proceed to the next orbital spacing calculation; if there is no frequency overlap area, determine whether all satellites in the basic database have been traversed. Step 6: Calculate the orbital interval Design arc high-risk interval vectors for different frequency bands; calculate the orbital interval between satellites based on the initial orbital position of the input primary frequency orbit and the satellite orbital position read from the database. If the orbital position is within the high-risk interval, perform coverage overlap assessment; if the orbital position is not within the high-risk interval, perform a judgment on "whether all satellites in the basic database have been traversed"; Step 7: Coverage overlap analysis Compare the coverage area of ​​the input primary frequency track with the coverage area of ​​the satellites in the database. If the coverage areas overlap, perform interference assessment; if the coverage areas do not overlap, determine whether all satellites in the basic database have been traversed. Step 8: Output evaluation data and results.

2. The GEO satellite frequency orbit resource availability assessment method according to claim 1, characterized in that: In step 2, trapezoidal fuzzy numbers are used to quantify the evaluation index, specifically: Among them, A i (x) represents the membership function, and a and b represent the satisfaction membership function parameters.

3. The GEO satellite frequency orbit resource availability assessment method according to claim 2, characterized in that: Step 3 determines the evaluation index weight based on the second-order coefficient of variation, namely: Here, v represents the coefficient of variation.

4. The GEO satellite frequency orbit resource availability assessment method according to claim 3, characterized in that: The method further includes generating an alternative plan, wherein the alternative plan has a different track position from the primary frequency and track plan, and the other parameters remain unchanged; the generation of the alternative plan is based on the following method: a) Select a frequency track whose orbital parameters do not conflict with the primary frequency track plan; b) Choose the conflicting frequency track with the least coordination difficulty.

5. A GEO satellite frequency orbit resource availability assessment system, characterized in that: The system has a program module corresponding to the steps of the method described in any one of claims 1 to 4, and executes the steps in the GEO satellite frequency orbit resource availability assessment method when running.

6. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and the computer program is configured to implement the steps in the GEO satellite frequency orbit resource availability assessment method according to any one of claims 1 to 4 when called by a processor.

Images