CN113162676B - GSO rail position efficiency evaluation method based on rail position multistage joint risk - Google Patents
GSO rail position efficiency evaluation method based on rail position multistage joint risk Download PDFInfo
- Publication number
- CN113162676B CN113162676B CN202110324683.6A CN202110324683A CN113162676B CN 113162676 B CN113162676 B CN 113162676B CN 202110324683 A CN202110324683 A CN 202110324683A CN 113162676 B CN113162676 B CN 113162676B
- Authority
- CN
- China
- Prior art keywords
- resource
- frequency
- orbit
- value
- track
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 116
- 239000011159 matrix material Substances 0.000 claims abstract description 87
- 239000013598 vector Substances 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007123 defense Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18519—Operations control, administration or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/02—Computing arrangements based on specific mathematical models using fuzzy logic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/60—Business processes related to postal services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Software Systems (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Networks & Wireless Communication (AREA)
- Astronomy & Astrophysics (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Automation & Control Theory (AREA)
- Biomedical Technology (AREA)
- Fuzzy Systems (AREA)
- Operations Research (AREA)
- Molecular Biology (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
Abstract
The invention discloses a GSO rail position efficiency evaluation method based on rail position multilevel joint risk, which comprises the following steps: establishing a risk arc segment table to record evaluation intervals of different frequency bands; aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an actual in-track satellite table in an evaluation range, an occupied space resource table in the evaluation range and a resource table to be coordinated; setting the risk degree of the on-orbit satellite, the risk degree of the occupied space resource and the risk degree of the resource to be coordinated as primary indexes, and setting the primary indexes and a decision matrix corresponding to the secondary indexes; respectively calculating weights for the decision matrixes of the first-level index and the second-level index; calculating the fuzzy satisfaction degree of each secondary index, setting grades and giving boundary values of the grades by experts; generating an evaluation matrix; and calculating the ambiguity vectors weighted by the secondary indexes and the primary indexes in sequence, and selecting the maximum value in the ambiguity vectors as an efficiency evaluation value. The method and the device can efficiently and specifically evaluate the use risk of the space frequency track resource.
Description
Technical Field
The invention relates to the technical field of space frequency track resources, in particular to a GSO (ground satellite based) track performance evaluation method based on track level multistage joint risk.
Background
Through development of more than half a century, the space frequency orbit resource has very important strategic significance for politics, economy, national defense construction and the like of the country and also becomes a valuable strategic resource which must be struggled by countries in the world, so that the countries in the world pay more attention to the development of satellites and the demand on the frequency orbit resource is increased. But because the method has the advantages of being finite and declaring that the cost of a space frequency track resource is high, a reasonable efficiency evaluation method is provided, which is helpful for people to know and master the possibility of putting a frequency track resource into use, so that the use value is clear, and a reliable basis is provided for the use decision of the frequency track resource. Therefore, the efficiency evaluation has very important significance for improving the effectiveness of national strategic resources.
The existing frequency orbit resource efficiency evaluation completely depends on manpower, a perfect method is not used for assisting, simple judgment is carried out by experts according to the existing experience, and then the orbit position can be clearly evaluated through actual declaration, coordination and use.
Disclosure of Invention
The invention aims to provide a GSO rail position efficiency evaluation method based on rail position multistage joint risk, so that the use risk of space frequency rail resources is effectively evaluated, and each stage of indexes under the use risk are quantitatively calculated.
The technical solution for realizing the purpose of the invention is as follows: a GSO rail position efficiency evaluation method based on rail position multilevel joint risk comprises the following steps:
step 1, establishing a risk arc segment table, and recording evaluation intervals corresponding to different frequency bands;
Step 3, aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an occupied space resource table in an evaluation range, wherein the field in the table comprises the actual number A of the satellites in orbit which accord with the occupied space resource condition 21 Minimum rail spacing A 22 Maximum interference value of uplink A 23 And a downlink maximum interference value A 24 ；
Step 4, needlingEstablishing a resource table to be coordinated according to the track position and frequency band information of the frequency orbit resource to be evaluated, wherein the field in the table comprises the actual number A of in-orbit satellites meeting the resource condition to be coordinated 31 Minimum rail position interval A 32 Maximum interference value of uplink A 33 And a downlink maximum interference value A 34 ；
Step 5, grading the indexes to be evaluated to ensure the risk degree A of the on-orbit satellite 1 Risk degree of occupied space resource A 2 Risk degree of resource to be coordinated A 3 Giving a decision matrix of the first-level index, giving a decision matrix corresponding to the second-level index for each first-level index, and recording the decision matrix in a database;
step 6, respectively calculating weights for the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index;
step 7, aiming at the first-level index A 1 、A 2 And A 3 Calculating fuzzy satisfaction of each secondary index, setting 5 grades, namely A, B, C, D and E, and giving boundary values of the grades by experts;
step 8, calculating ambiguity vectors aiming at the values obtained in the steps 2 to 4 respectively to generate an evaluation matrix;
step 9, calculating a second-level index weighted ambiguity vector based on the second-level index weight;
and step 10, calculating a primary index weighted ambiguity vector based on the primary index weight, and selecting the maximum value in the ambiguity vector as an efficiency evaluation value.
Compared with the prior art, the invention has the following remarkable advantages: (1) reasonable evaluation can be accurately carried out on a space frequency track resource, and the evaluation system flow is complete; (2) the operation is simple, the realization is convenient, and the requirement on the professional knowledge of a user is low; (3) the method is suitable for large-scale use, has strong applicability, and greatly improves the accuracy of frequency resource declaration.
Drawings
FIG. 1 is a general operational framework diagram of the GSO rail level performance evaluation method based on rail level multi-level joint risk according to the present invention.
Fig. 2 is an evaluation index grading diagram of the GSO rail performance evaluation method based on rail level multi-level joint risk according to the present invention.
Fig. 3 is a schematic diagram of an index function of the GSO rail performance evaluation method based on rail level multi-level joint risk according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and the embodiments.
With reference to fig. 1, the GSO rail position performance evaluation method based on rail position multilevel joint risk of the present invention includes the steps of:
step 1, establishing a risk arc segment table, and recording evaluation intervals corresponding to different frequency bands;
Step 3, aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an occupied space resource table in an evaluation range, wherein the field in the table comprises the actual number A of the satellites in orbit which accord with the occupied space resource condition 21 Minimum rail spacing A 22 Maximum interference value of uplink A 23 And a downlink maximum interference value A 24 ；
Step 4, aiming at the track position and frequency band information of the frequency orbit resource to be evaluated, establishing a resource table to be coordinated, wherein the field in the table comprises the actual number A of in-orbit satellites meeting the resource condition to be coordinated 31 Minimum rail spacing A 32 Maximum interference value of uplink A 33 And a downlink maximum interference value A 34 ；
Step 5, grading the indexes to be evaluated to ensure the risk degree A of the on-orbit satellite 1 Risk degree of occupied space resource A 2 Risk degree of resource to be coordinated A 3 Giving a decision matrix of the first-level index, giving a decision matrix corresponding to the second-level index for each first-level index, and recording the decision matrix in a database;
step 6, respectively calculating weights for the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index;
step 7, aiming at the first-level index A 1 、A 2 And A 3 Calculating fuzzy satisfaction of each secondary index, setting 5 grades, namely A, B, C, D and E, and giving boundary values of the grades by experts;
step 8, calculating ambiguity vectors aiming at the values obtained in the steps 2 to 4 respectively to generate an evaluation matrix;
step 9, calculating a second-level index weighted ambiguity vector based on the second-level index weight;
And step 10, calculating a primary index weighted ambiguity vector based on the primary index weight, and selecting the maximum value in the ambiguity vector as an efficiency evaluation value.
Further, as a specific example, step 1 establishes a risk arc segment table, and records evaluation intervals corresponding to different frequency bands, as specifically shown in table 1:
TABLE 1 Risk arc table
The fields in the table sequentially comprise a sequence number id, a frequency band frequency, a frequency lowest value freq _ from, a frequency highest value freq _ to and an evaluation arc segment interval around _ eval, wherein the units of the freq _ from and the freq _ to are MHz, and the units of the around _ eval are degrees.
Further, as a specific example, in step 2, for the orbit position and frequency band information of the frequency orbit resource to be evaluated, an actual orbit satellite table within the evaluation range is established, and a field in the table includes the number a of actual orbit satellites 11 Minimum orbital spacing A of actual in-orbit satellites 12 The construction process is as follows:
(1) searching for on-orbit satellites within evaluation range from UCS database
The UCS database is a global on-orbit satellite database periodically provided by an American scientist website www.ucsusa.org, frequency band information is extracted from a comment field of the UCS database, and when a frequency band comprises a target frequency band, whether the track position is in an evaluation range is continuously judged, wherein the evaluation range is [ track position-6 of a frequency track resource to be evaluated, and track position +6 of the frequency track resource to be evaluated ]. When the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
(2) Finding in-orbit satellites within evaluation range from satbeam database
The website www.satbeams.com/satellites provides a global specialty database, satbeam, that includes satellite overlays, recording information for approximately over 380 geostationary satellites.
And extracting frequency band information from comments fields of the satpeak database, and when the frequency band comprises a target frequency band, continuously judging whether the track position is in an evaluation range, wherein the evaluation range is [ track position-6 of the frequency track resource to be evaluated, and track position +6 of the frequency track resource to be evaluated ]. When the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
after the actual in-orbit satellite table is established, order A 11 Equal to the number of satellites actually in orbit, A 12 Actual on-orbit satellites minimum orbital spacing.
Further, as a specific example, the step 3 establishes an occupied space resource table within an evaluation range for the track position and frequency band information of the frequency track resource to be evaluated, specifically as follows:
the fields in the table include the satellite name sat _ name and the track interval orb _ sep, the uplink interference up _ interference and the downlink interference down _ interference. Looking up the long _ nom field of srs database geo table provided by ITU (International Telecommunication Union), when the long _ nom value is within the evaluation range, recording the corresponding ntc _ id, and substituting ntc _ id into the notice table to look up ntf _ rsn value corresponding to ntc _ id, the following two cases will occur:
In the first case: when ntf _ rsn is N, it indicates that the track N data is enabled;
in the second case: when ntf _ rsn is C, ntc _ id at the moment is substituted into a com _ el table, the adm field and sat _ name field of the com _ el are recorded and substituted into a biu table in an ITU website, a status value with the same adm and sat _ name is searched, and when the status is C, the track N data is enabled;
the srs database is a frequency track resource database, geo is a geostationary space station table, long _ nom represents the longitude of a space station, notice is a general data table for notification, ntc _ id represents a unique identifier of notice, and ntf _ rsn represents a notification reason; the N data is the notification registration information which is provided by the satellite using main part to the ITU and is used in the satellite network in brief and actually; c represents satellite network coordination data;
when the N data is started, ntc _ id is substituted into an s _ beam table, all freq _ min and freq _ max corresponding to ntc _ id are searched, freq _ min-freq _ max are substituted into a risk arc segment table, and whether the corresponding frequency is the same as the target frequency band or not is checked; wherein s _ beam is a satellite antenna beam table, freq _ min is the lowest frequency of a beam, freq _ max is the highest frequency of the beam, and frequency is a frequency band;
When the frequency ranges are the same, recording the satellite name and the orbit space of the orbit resource to be evaluated into an occupied resource table, and recording A 21 Number of actual orbiting satellites, A, to meet the conditions of occupied space resources 22 Is the minimum rail position interval; then, aiming at the track resource with the minimum track position interval value and the track resource to be evaluated, respectively selecting the beam with the maximum interference value from the uplink and the downlink by calculation, respectively calculating the interference values between the two beams of the uplink and the downlink, and storing the interference values as the uplink interference value and the downlink interference value of the occupied space resource, wherein the calculation process specifically comprises the following steps:
(1.1) respectively substituting the beam _ names into a grp table, searching all grp _ ids under the beam _ names, then searching freq _ min and freq _ max corresponding to the grp _ ids, substituting freq _ min-freq _ max into a risk arc segment table, and judging whether the frequency band is the same as that of the track resource to be evaluated; substituting grp _ id with the same frequency band into an e _ as _ stn table, searching for gain under the grp _ id, comparing the gains of all grp _ id, and recording the maximum gain and the corresponding bmwdth; wherein, the grp table is an attribute table of the repeater group under the beam, grp _ id is a unique identifier of the repeater group, freq _ min is a lowest frequency value, freq _ max is a highest frequency value, e _ as _ stn table is an attribute table of the earth station corresponding to the repeater group, gain is gain of the repeater group, and bmwdth is the beam width of the earth station antenna;
(1.2) searching pwr _ ds _ max corresponding to grp _ id with the maximum gain in an emiss table, wherein the emiss table is a carrier attribute table, and pwr _ ds _ max is the maximum power density;
(1.3) judging whether the value of the track spacing 1.15 times is in the range of [1,100 λ/D ], λ represents the wavelength and D represents the radius, and calculating the gain recorded in the (1.1) th step and the pwr _ ds _ max recorded in the (1.2) th step as follows:
gain+pwr_ds_max+29-25log(A 22 *1.15)
wherein 100 λ/D — bmwdth 25/9, bmwdth being recorded in step (1.1);
comparing the values calculated by the formula in the step (1.3), selecting the beam with the maximum value to calculate the interference, and enabling A 23 Represents the uplink maximum interference value:
A 23 ＝max{I/N up ,I/N′ up }、
I/N up ＝G(θ) Rxsat +P TxES +G(φ) TxES -FSL up -T sat -B carrier -k
I/N′ up ＝G′(θ) Rxsat +P′ TxES +G′(φ) TxES -FSL′ up -T′ sat -B′ carrier -k
where I represents the total interference power and N up Representing the noise power, N ', produced by the rail resource with the minimum rail bit interval value' up Representing the noise power generated by the frequency track resource to be evaluated;
G(θ) Rxsat representing the reception gain, P, of the frequency-orbital resources to be evaluated in the direction of the earth station corresponding to the frequency-orbital resources having the smallest value of the orbital spacing TxES Representing the transmission power of the frequency-orbit resource corresponding to the earth station with the minimum orbit spacing value, G (phi) TxES Indicating the side lobe gain, T, of the frequency-orbital resource corresponding to the earth station with the minimum orbital spacing value sat 、B carrier And FSL up Respectively representing the noise temperature, the carrier bandwidth and the free space loss of an uplink of a frequency-track resource system to be evaluated, wherein K represents a boltzmann constant, and K is-228.6 dBJ/K; g' (theta) Rxsat Representing the reception gain, P, of the frequency-orbital resource having the minimum value of the orbital spacing in the direction of the earth station corresponding to the frequency-orbital resource to be evaluated T ′ xES Represents the transmission power, G' (φ), of the earth station corresponding to the frequency-orbit resource to be evaluated TxES Representing the side lobe gain, T 'of the corresponding earth station with the frequency track resource to be evaluated' sat 、B′ carrier And FSL' up Respectively representing the noise temperature, the carrier bandwidth and the free space loss of an uplink of the frequency-rail resource system with the minimum rail position interval value;
2. for the downlink, substituting ntc _ id of the track resource into the s _ beam table, looking up beam _ name of E at emi _ rcp under ntc _ id, where E indicates that the beam belongs to the downlink, and then performing the following operation on each beam _ name respectively:
(2.1) recording gain corresponding to the beam _ name;
(2.2) substituting the beam _ name into the grp table, searching all corresponding grp _ ids, substituting all grp _ ids into the emiss table, and searching the corresponding maximum pwr _ ds _ max;
(2.3) adding the gain recorded in the step (2.1) and the pwr _ ds _ max recorded in the step (2.2);
comparing the values calculated in step (2.3), selecting the beam with the maximum value to calculate interference, and enabling A 24 Represents the downlink maximum interference value:
A 24 ＝max{I/N down ,I/N′ down }
I/N down ＝G(α) RxES +P sat +G(φ) Txsat -FSL down -T sat -B carrier -k
I/N′ down ＝G′(α) RxES +P′ sat +G′(φ) Txsat -FSL′ down -T′ sat -B′ carrier -k
where I represents the total interference power and N up Representing the noise power, N ', produced by the rail resource with the minimum rail bit interval value' up Representing the noise power generated by the frequency track resource to be evaluated;
G(α) RxES represents the side lobe gain, P, of the earth station corresponding to the frequency orbit resource to be evaluated in the frequency orbit resource direction with the minimum orbit position interval value sat Representing the transmitted power of the frequency-rail resource with the minimum rail position interval value, G (phi) Txsat Representing antenna gain, T, of frequency-rail resource with minimum rail-bit spacing value sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of a downlink of a frequency-track resource system to be evaluated, wherein K represents a boltzmann constant, and K is-228.6 dBJ/K and has the unit of dBJ/K; g' (alpha) RxES Representing the side lobe gain, P 'of the earth station corresponding to the minimum orbit bit spacing value in the direction of the frequency-orbit resource to be evaluated' sat Representing the transmission power of the frequency-rail resource to be evaluated, G' (phi) Txsat Antenna gain, T, representing the frequency-track resource to be evaluated sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of the downlink of the frequency-track resource with the minimum track bit interval value.
Further, as a specific example, in step 4, a to-be-coordinated resource table is established for the track position and frequency band information of the frequency track resource to be evaluated, specifically as follows:
looking up srs a long _ nom field of a geo table of a database, and recording ntc _ id corresponding to the long _ nom when the long _ nom is in an evaluation range; substituting ntc _ id into a node table to search a ntf _ rsn value corresponding to the ntc _ id, wherein when ntf _ rsn is C, C indicates that the frequency-track resource has applied for coordination data, substituting ntc _ id of the C data frequency band into a com _ el table, wherein the com _ el table is a public content table, searching and recording sat _ name corresponding to the ntc _ id, searching the sat _ name again in the com _ el table, and when sat _ name is the same and ntf _ rsn is N, skipping the record and searching other frequency-track resources in an evaluation range; when there is no investigation When the sat _ name is the same and the value ntf _ rsn is N, substituting ntc _ id into an s _ beam table, searching all freq _ min-freq _ max corresponding to ntc _ id, substituting the freq _ min-freq _ max into a risk arc segment table, and checking whether the corresponding frequecy is the same as the target frequency band; when the frequency ranges are the same, recording the satellite name and the track position interval of the orbit resource to be evaluated into an occupied resource table, and enabling A to be 31 Representing the number of actual in-orbit satellites that meet the conditions of the resource to be coordinated, A 32 Representing a minimum rail bit interval;
then, aiming at the track resource with the minimum track position interval value and the track resource to be evaluated, selecting the beam with the maximum interference value from the uplink and the downlink respectively, calculating the interference value between the two beams of the uplink and the downlink, storing the interference value as the uplink interference value and the downlink interference value of the resource to be coordinated, wherein the calculation process is the same as the process of calculating the interference value in the step 3, and enabling A to be 33 Represents the maximum interference value of the uplink, A 34 Representing the downlink maximum interference value.
Further, as a specific example, in step 5, the to-be-evaluated indexes are classified, and as shown in fig. 2, the risk degree a of the on-orbit satellite is ordered 1 Risk degree of occupied space resource A 2 Risk degree of resource to be coordinated A 3 Giving a decision matrix of the first-level index, giving a decision matrix corresponding to the second-level index for each first-level index, and recording the decision matrix into a database, wherein the specific steps are as follows:
let the risk degree A of the on-orbit satellite 1 The number A of the actual in-orbit satellites is a first-level index 11 Minimum orbital spacing of actual on-orbit satellites A 12 Is A 1 The secondary index of (1);
let the risk degree A of occupied space resource 2 The number A of the actual orbit satellites meeting the occupied space resource condition is a first-level index 21 Minimum rail spacing A 22 Maximum interference value of uplink A 23 And a downlink maximum interference value A 24 Is A 2 The secondary index of (1);
let the resource risk degree A to be coordinated 3 The number of actual in-orbit satellites which are first-level indexes and meet the condition of the resource to be coordinatedQuantity A 31 Minimum rail spacing A 32 Maximum interference value of uplink A 33 And a downlink maximum interference value A 34 Is A 2 The secondary index of (1);
the decision matrix of the first-level index is as shown in table 2:
TABLE 2
| A 1 | A 2 | A 3 | |
| A 1 | 1 | 5 | 7 |
| A 2 | 1/5 | 1 | 5 |
| A 3 | 1/7 | 1/5 | 1 |
A 1 The decision matrix of the next two-level index is shown in table 3:
TABLE 3
| A 1 | A 11 | A 12 |
| A 11 | 1 | 1/3 |
| A 12 | 3 | 1 |
A 2 The decision matrix of the next two-level index is shown in table 4:
TABLE 4
| A 2 | A 21 | A 22 | A 23 | A 24 |
| A 21 | 1 | 1/3 | 1/3 | 1/3 |
| A 22 | 3 | 1 | 1 | 1 |
| A 23 | 3 | 1 | 1 | 1 |
| A 24 | 3 | 1 | 1 | 1 |
A 3 The decision matrix of the next two-level index is shown in table 5:
TABLE 5
| A 3 | A 31 | A 32 | A 33 | A 34 |
| A 31 | 1 | 1/3 | 1/3 | 1/3 |
| A 32 | 3 | 1 | 1 | 1 |
| A 33 | 3 | 1 | 1 | 1 |
| A 34 | 3 | 1 | 1 | 1 |
Further, as a specific example, the weights of the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index in step 6 are respectively calculated as follows:
For A 1 The decision matrix of the next two-level index firstly normalizes the value of each column of the matrix to omega ij ，
Sum by row, ω 1 ＝ω 11 +ω 12 ，ω 2 ＝ω 21 +ω 22 Will ω i Normalized to obtain
To obtain A 1 Weight vector of the next two-level index
By the same token, respectively obtain A 2 Weight vector of the next two-level index
A 3 Weight vector of the next two-level index
Weight vector of primary index decision matrix
Further, as a specific example, step 7 describes for the primary index A 1 、A 2 And A 3 Calculating fuzzy satisfaction degree of each secondary index, setting 5 grades, namely A, B, C, D and E, and giving boundary values of the grades by experts, namely a in figure 3 1 ,b 2 ,a 2 ,b 3 ,a 3 ,b 4 ,a 4 ,b 5 The method comprises the following steps:
A 11 and A 21 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 31 Has a boundary value of [1,4,6,8,10,12,14,18 ]]，A 12 、A 22 And A 32 Has a boundary value of [0.5,1,1.5,2,2.5,3,3.5,4 ]]，A 23 And A 24 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 33 And A 34 Has a boundary value of [1,4,6,8,10,12,14,18 ]]。
Further, as a specific example, in step 8, the ambiguity vectors are respectively calculated for the values obtained in steps 2 to 4, and an evaluation matrix is generated, specifically as follows:
1. for A 11 、A 21 、A 23 、A 24 、A 31 、A 33 And A 34 Respectively setting i as 1,2,3,4 and 5, and sequentially substituting the following formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, wherein the evaluation matrix is r 11 、r 21 、r 23 、r 24 、r 31 、r 33 、r 34 ；
2. For A 12 、A 22 、A 32 Respectively setting i as 1,2,3,4 and 5, sequentially substituting formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, and then respectively subtracting the values in the matrix from 1 to obtain a new evaluation matrix, wherein the evaluation matrix is r 12 、r 22 And r 32 ；
Further, as a specific example, step 9 calculates a secondary-index-weighted ambiguity vector based on the secondary index weight, specifically as follows:
summarizing the evaluation matrixes in the step 8, and respectively ordering the matrixes R 1 、R 2 And R 3 Is an index A 1 、A 2 And A 3 In which R is 1 ＝[r 11 ,r 12 ]，R 2 ＝[r 21 ,r 22 ,r 23 ,r 24 ]，R 3 ＝[r 31 ,r 32 ,r 33 ,r 34 ](ii) a Let beta 1 、β 2 And beta 3 In the form of a vector of degrees of ambiguity,
moment recording deviceArray B ═ beta 1 ,β 2 ,β 3 ) Wherein
And
are all obtained in step 6.
Further, as a specific example, in step 10, the first-level index weight-based ambiguity vector is calculated, and a maximum value in the ambiguity vector is selected as the performance evaluation value, specifically:
based on first-level index weight
Calculating a first-level index weighted ambiguity vector; order matrix
And selecting the maximum value in the matrix C as the performance evaluation value.
The invention is described in further detail below with reference to the figures and the embodiments.
Examples
As shown in fig. 1, the overall operational framework of the method comprises the following steps:
Step 1, establishing a risk arc segment table, and recording evaluation intervals corresponding to different frequency bands as shown in table 1. The fields in the table sequentially comprise a sequence number id, a frequency band frequency, a frequency lowest value freq _ from, a frequency highest value freq _ to and an evaluation arc segment interval around _ eval, wherein the units of the freq _ from and the freq _ to are MHz, and the units of the around _ eval are degrees.
And 2, aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an actual in-track satellite table in the evaluation range, wherein fields in the table comprise a satellite name sat _ name and a track position interval orb _ sep, and the unit of the track position interval is degree. Taking the network data ID of the frequency orbit resource to be evaluated as 90500037 as an example, the name of the satellite network is PALAPA-B1, the orbit position is 108 degrees E, and the frequency band to which the satellite network belongs is the C frequency band, an actual in-orbit satellite table is constructed, and the construction process is as follows:
(1) searching for on-orbit satellites within evaluation range from UCS database
The UCS database is a global orbiting satellite database periodically provided by an American scientist website www.ucsusa.org, frequency band information is extracted from a comment field of the UCS database, and when a frequency band comprises a C frequency band, whether the track position is in an evaluation range is continuously judged, wherein the evaluation range is [102 DEG E, 114 DEG E ]. When the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
(2) Finding in-orbit satellites within evaluation range from satbeam database
The website www.satbeams.com/satellites provides a global specialty database, satbeam, that includes satellite overlays, recording information for approximately over 380 geostationary satellites. And extracting frequency band information from the comments field of the satpeak database, and when the frequency band comprises a C frequency band, continuously judging whether the track is in an evaluation range, wherein the evaluation range is [102 DEG E, 114 DEG E ]. When the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
after the search, the number of the qualified actual orbit satellites is 7, and the minimum orbit position interval is 0.21 degrees, so that A 11 ＝7，A 12 ＝0.21。
And 3, aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an occupied space resource table in the evaluation range, wherein fields in the table comprise a satellite name sat _ name, a track position interval orb _ sep, an uplink interference up _ interference and a downlink interference down _ interference. Looking up the long _ nom field of srs database geo table provided by ITU (International telecommunications Union), when the long _ nom value is within the evaluation range, i.e., [102 ° E, 114 ° E ], recording the corresponding ntc _ id, and substituting ntc _ id into the notice table to look up ntf _ rsn value corresponding to ntc _ id, the following two cases will occur:
1. When ntf _ rsn is N, it indicates that the track N data is enabled;
2. when ntf _ rsn is C, ntc _ id at the moment is substituted into a com _ el table, the adm field and sat _ name field of the com _ el are recorded and substituted into a biu table in an ITU website, a status value with the same adm and sat _ name is searched, and when the status is C, the track N data is enabled;
when the N data is enabled, ntc _ id is substituted into the s _ beam table, all freq _ min and freq _ max corresponding to ntc _ id are searched, freq _ min-freq _ max are substituted into the risk arc segment table, and whether the corresponding frequecy is the C frequency band or not is checked.
And when the frequency band is the C frequency band, recording the name of the satellite and the track bit interval between the satellite and the track resource to be evaluated into an occupied resource table. Through the search, the number of the actual orbit satellites meeting the conditions is 23, and the minimum orbit position interval is 0.5 degrees, so that A 21 ＝23，A 22 ＝0.5。
Then, for the track resource with the minimum rail position interval value and the track resource to be evaluated, selecting the beam with the maximum interference value from the uplink and the downlink respectively, and calculating the interference value between the two beams as follows:
1. substituting ntc _ id of the track resource into an s _ beam table, searching beam _ name of ntc _ id lower emi _ rcp ═ R, and then respectively performing the following operations on each beam _ name:
(1) And respectively substituting the beam _ name into the grp table, searching all grp _ ids under the beam _ name, then searching freq _ min and freq _ max corresponding to the grp _ ids, substituting the freq _ min-freq _ max into the risk arc segment table, and judging whether the frequency band is the C frequency band or not. When the frequency band is the C frequency band, substituting grp _ id into the e _ as _ stn table, searching for gain under the grp _ id, comparing the gains of all grp _ id, and recording the maximum gain and the corresponding bmwdth;
(2) searching pwr _ ds _ max corresponding to grp _ id with the maximum gain in an emiss table;
(3) and judging whether the track bit interval value is within the range of [1, 100 lambda/D ] or not, and calculating gain recorded in the step (1) and pwr _ ds _ max recorded in the step (2) when the track bit interval value is within the range, namely gain + pwr _ ds _ max +29-25log (A22 is 1.15). Wherein 100 λ/D — bmwdth × 25/9, bmwdth being recorded in step (1);
comparing the values calculated in step (3) and selecting the beam having the maximum valueCalculating interference, let A 23 Represents an uplink interference value, A 23 ＝max{I/N up ,I/N′ up }，I/N up ＝G(θ) Rxsat +P TxES +G(φ) TxES -FSL up -T sat -B carrier -k，I/N′ up ＝G′(θ) Rxsat +P′ TxES +G′(φ) TxES -FSL′ up -T s ′ at -B′ carrier -k, wherein G (θ) Rxsat Representing the reception gain, P, of the frequency-orbital resources to be evaluated in the direction of the earth station corresponding to the frequency-orbital resources having the smallest value of the orbital spacing TxES Representing the transmission power of the frequency-orbit resource corresponding to the earth station with the minimum orbit spacing value, G (phi) TxES Indicating the side lobe gain, FSL, of the frequency-orbital resource corresponding to the earth station with the minimum orbital spacing value up 、T sat And B carrier Respectively representing the free space loss, the noise temperature and the carrier bandwidth of an uplink of the frequency-track resource system to be evaluated, K represents a boltzmann constant, K is-228.6 (dBJ/K), G' (theta) Rxsat Representing the receiving gain, P ', of the frequency-orbit resource with the minimum orbit bit interval value in the direction of the corresponding earth station of the frequency-orbit resource to be evaluated' TxES Indicating the transmission power, G' (phi), of the frequency-orbital resource to be evaluated relative to the earth station TxES Representing the side lobe gain, T 'of the earth station corresponding to the frequency track resource to be evaluated' sat 、B′ carrier And FSL' up Respectively representing the noise temperature, the carrier bandwidth and the free space loss of an uplink of the frequency-rail resource system with the minimum rail bit interval value. Calculated to obtain A 23 ＝45.542；
2. For the downlink, substituting ntc _ id of the track resource into an s _ beam table, looking up beam _ name of emi _ rcp ═ E under ntc _ id, and then performing the following operation on each beam _ name respectively:
(1) recording gain corresponding to the beam _ name;
(2) substituting the beam _ name into a grp table, searching for all corresponding grp _ ids, substituting all grp _ ids into an emiss table, and searching for the corresponding maximum pwr _ ds _ max;
(3) adding the gain recorded in the step (1) and pwr _ ds _ max recorded in the step (2);
comparing the values calculated in step (3), selecting the beam with the maximum value to calculate interference, and enabling A 24 Represents a downlink interference value, A 24 ＝max{I/N down ,I/N′ down }，I/N down ＝G(α) RxES +P sat +G(φ) Txsat -FSL down -T sat -B carrier -k，I/N′ down ＝G′(α) RxES +P′ sat +G′(φ) Txsat -FSL′ down -T′ sat -B′ carrier -k, wherein G (α) RxES Represents the side lobe gain, P, of the earth station corresponding to the frequency-orbit resource to be evaluated in the direction of the frequency-orbit resource with the minimum orbit position interval value sat Representing the transmitted power of the frequency-rail resource with the minimum rail position interval value, G (phi) Txsat Representing antenna gain, T, of frequency-rail resource with minimum rail-bit spacing value sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of a downlink of a frequency-track resource system to be evaluated, K represents a boltzmann constant, K is-228.6 (dBJ/K), G' (alpha) RxES Representing the side lobe gain, P 'of the earth station corresponding to the minimum orbit bit spacing value in the direction of the frequency-orbit resource to be evaluated' sat Representing the transmission power of the frequency-rail resource to be evaluated, G' (phi) Txsat Antenna gain, T, representing the frequency-track resource to be evaluated sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of the downlink of the frequency-track resource with the minimum track bit interval value. Calculated to obtain A 24 ＝0。
And 4, establishing a resource table to be coordinated according to the track position and frequency band information of the frequency track resource to be evaluated. The fields in the table include the satellite name sat _ name and the track interval orb _ sep, the uplink interference up _ interference and the downlink interference down _ interference. Looking up srs the long _ nom field of the database geo table, and recording the ntc _ id corresponding to the time when the long _ nom is in the evaluation range. Substituting ntc _ id into the node table to search ntf _ rsn value corresponding to ntc _ id, when ntf _ rsn is C, substituting ntc _ id of the N data frequency band into the s _ beam table to search ntc _ id corresponding to all freq _ m And in-frequency _ max, substituting the frequency _ min-frequency _ max into a risk arc segment table, and checking whether the corresponding frequency is in the C frequency band or not. And when the frequency band is the C frequency band, recording the satellite name and the track bit interval with the orbit resource to be evaluated into an occupied resource table. Through the search, the number of the resources to be coordinated which meet the conditions is 29, and the minimum track spacing is 0.2 degrees, so that A 31 ＝29，A 22 ＝0.2。
And then, aiming at the track resource with the minimum track position interval value and the track resource to be evaluated, selecting the beam with the maximum interference value from the uplink and the downlink respectively, and calculating the interference value between the two beams, wherein the calculation process is the same as the process of calculating the interference value in the step 3. Is calculated to obtain A 33 ＝0，A 34 ＝0。
And 5, grading the indexes to be evaluated as shown in figure 2.
Giving a decision matrix of the primary indexes, giving a decision matrix of the secondary indexes of each primary index, and recording the decision matrix into a database.
Decision matrix of primary index is shown in Table 2, A 1 Decision matrix of next two-level index is shown in Table 3, A 2 Decision matrix of next two-level index is shown in Table 4, A 3 The decision matrix of the next two-level index is shown in table 5.
And 6, respectively calculating the weight of the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index. The values of each column of the matrix are first normalized to ω ij ，
Will omega ij Summing by rows, e.g. ω 1 ＝ω 11 +ω 12 After that, ω will be i Normalized to obtain
By calculation, the results are as follows:
the decision matrix and weights of the first-level indexes are as shown in table 6:
TABLE 6
The decision matrix and weights of the following two-level indexes A1 are shown in Table 7:
TABLE 7
The decision matrix and weights of the secondary indexes under A2 are shown in Table 8:
TABLE 8
The decision matrix and weights of the next two-level indexes of A3 are shown in Table 9:
TABLE 9
Step 7, for A 1 、A 2 And A 3 Calculating each secondary index A ij The blur satisfaction degree of (c). 5 levels, A, B, C, D and E, are set, and the boundary value of each level, a in FIG. 3, is given by the expert 1 ,b 2 ,a 2 ,b 3 ,a 3 ,b 4 ,a 4 ,b 5 . The specific values are as follows: a. the 11 And A 21 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 31 Has a boundary value of [1,4,6,8,10,12,14,18 ]]，A 12 、A 22 And A 32 Has a boundary value of [0.5,1,1.5,2,2.5,3,3.5,4 ]]，A 23 And A 24 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 33 And A 34 Has a boundary value of [1,4,6,8,10,12,14,18 ]]。
Step 8, calculating ambiguity vectors respectively according to the values obtained in the steps 2 to 4, generating an evaluation matrix, and calculating as follows:
1) for A 11 、A 21 、A 23 、A 24 、A 31 、A 33 And A 34 Respectively setting i as 1,2,3,4 and 5, sequentially substituting formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, and recording the evaluation matrix as r 11 、r 21 …r 34 . Through calculation, the following results are obtained: r is 11 ＝[0,0,0,0.5,0.5]，r 21 ＝[0,0,0,0,1]，r 23 ＝[0,0,0,0,1]，r 24 ＝[1,0,0,0,0]，r 31 ＝[0,0,0,0,1]，r 33 ＝[1,0,0,0,0]，r 34 ＝[1,0,0,0,0]；
2) For A 12 、A 22 、A 32 Respectively setting i as 1,2,3,4 and 5, sequentially substituting formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, then respectively subtracting the values in the matrix from 1 to obtain a new evaluation matrix, and recording the evaluation matrix as r 12 、r 22 And r 32 . Through calculation, the following results are obtained: r is 12 ＝[0,1,1,1,1]，r 22 ＝[0,1,1,1,1]，r 32 ＝[0,1,1,1,1]。
And 9, calculating a weighted ambiguity vector according to the secondary index weight. Summarizing the evaluation matrixes in the step 8, and respectively ordering the matrixes R 1 、R 2 And R 3 Is an index A 1 、A 2 And A 3 The evaluation matrix of (1). Wherein R is 1 ＝[r 11 ,r 12 ]，R 2 ＝[r 21 ,r 22 ,r 23 ,r 24 ]，R 3 ＝[r 31 ,r 32 ,r 33 ,r 34 ]. Let beta 1 、β 2 And beta 3 In the form of a vector of degrees of ambiguity,
let matrix B ═ β 1 ,β 2 ,β 3 ). Wherein,
and
are all obtained in step 6. Calculated, matrix B is as in table 10:
watch 10
| β 1 | 0 | 0.75 | 0.75 | 0.88 |
| β 2 | 0.3 | 0.3 | 0.3 | 0.3 |
| β 3 | 0.6 | 0.3 | 0.3 | 0.3 |
Step 10, order matrix
By calculation, the matrix C is [0.12,0.61,0.61,0.69,0.8 ═ C]. The maximum value in the matrix C is selected as the performance evaluation value, so the final performance evaluation value of the matrix is 0.8.
The invention can accurately make reasonable evaluation on the space frequency orbit resource, and the evaluation system has perfect flow; the operation is simple, the realization is convenient, and the requirement on the professional knowledge of a user is low; the method is suitable for large-scale use, has strong applicability, and greatly improves the accuracy of frequency resource declaration.
Claims (10)
1. A GSO rail position efficiency evaluation method based on rail position multilevel joint risk is characterized by comprising the following steps:
Step 1, establishing a risk arc segment table, and recording evaluation intervals corresponding to different frequency bands;
step 2, aiming at the orbit position and frequency band information of the frequency orbit resource to be evaluated, establishing an actual orbit satellite table in an evaluation range, wherein the field in the table comprises the actual orbit satellite quantity A 11 Minimum orbital spacing of actual on-orbit satellites A 12 ；
Step 3, aiming at the track position and frequency band information of the frequency track resource to be evaluated, establishing an occupied space resource table in an evaluation range, wherein the field in the table comprises the actual number A of the satellites in orbit which accord with the occupied space resource condition 21 Minimum rail spacing A 22 Maximum interference value of uplink A 23 And a downlink maximum interference value A 24 ；
Step 4, aiming at the track position and frequency band information of the frequency orbit resource to be evaluated, establishing a resource table to be coordinated, wherein the field in the table comprises the actual number A of in-orbit satellites meeting the resource condition to be coordinated 31 Minimum, isTrack spacing A 32 Maximum interference value of uplink A 33 And a downlink maximum interference value A 34 ；
Step 5, grading the indexes to be evaluated to ensure the risk degree A of the on-orbit satellite 1 Risk degree of occupied space resource A 2 Risk degree of resource to be coordinated A 3 Giving a decision matrix of the first-level index, giving a decision matrix corresponding to the second-level index for each first-level index, and recording the decision matrix in a database;
Step 6, respectively calculating weights for the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index;
step 7, aiming at the first-level index A 1 、A 2 And A 3 Calculating fuzzy satisfaction of each secondary index, setting 5 grades, namely A, B, C, D and E, and giving boundary values of the grades by experts;
step 8, calculating ambiguity vectors aiming at the values obtained in the steps 2 to 4 respectively to generate an evaluation matrix;
step 9, calculating a second-level index weighted ambiguity vector based on the second-level index weight;
and step 10, calculating a primary index weighted ambiguity vector based on the primary index weight, and selecting the maximum value in the ambiguity vector as an efficiency evaluation value.
2. The GSO rail position effectiveness evaluation method based on rail position multistage joint risk according to claim 1, wherein the step 1 of establishing a risk arc segment table records evaluation intervals corresponding to different frequency bands, as shown in table 1:
TABLE 1 Risk arc Table
The fields in the table sequentially comprise a sequence number id, a frequency band frequency, a frequency lowest value freq _ from, a frequency highest value freq _ to and an evaluation arc segment interval around _ eval, wherein the units of the freq _ from and the freq _ to are MHz, and the units of the around _ eval are degrees.
3. The GSO orbit performance evaluation method based on orbit multistage joint risk according to claim 1, wherein step 2 establishes an actual orbit satellite table within an evaluation range for the orbit position and frequency band information of the frequency orbit resource to be evaluated, specifically as follows:
1. searching for on-orbit satellites in evaluation range from UCS database
Extracting frequency band information from comments fields of a UCS database, and when the frequency band comprises a target frequency band, continuously judging whether the track position is in an evaluation range, wherein the evaluation range is [ track position-6 of a frequency track resource to be evaluated, and track position +6 of the frequency track resource to be evaluated ]; when the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
2. finding in-orbit satellites within evaluation range from satbeam database
Extracting frequency band information from comments fields of the satpeak database, and when the frequency band comprises a target frequency band, continuously judging whether the track position is in an evaluation range, wherein the evaluation range is [ track position-6 of the frequency track resource to be evaluated, and track position +6 of the frequency track resource to be evaluated ]; when the orbit position is in the evaluation range, recording the name of the satellite and the orbit position interval between the satellite and the orbit resource to be evaluated to an actual orbit satellite table;
After the actual in-orbit satellite table is established, order A 11 Is equal to the number of actual in-orbit satellites, A 12 Actual on-orbit satellites minimum orbital spacing.
4. The GSO rail position performance evaluation method based on rail position multilevel joint risk according to claim 1, wherein step 3 establishes an occupied space resource table within an evaluation range for the rail position and frequency band information of the frequency rail resource to be evaluated, specifically as follows:
fields in the occupied space resource table comprise a satellite name sat _ name, a track interval orb _ sep, an uplink interference up _ interference and a downlink interference down _ interference; looking up the rail longitude long _ nom field of an srs database geo table provided by the ITU, recording the corresponding ntc _ id when the long _ nom value is within the evaluation range, and substituting ntc _ id into the notice table to look up the ntf _ rsn value corresponding to ntc _ id, the following two conditions occur:
in the first case: when ntf _ rsn is N, it indicates that the track N data is enabled;
in the second case: when ntf _ rsn is C, ntc _ id at the moment is substituted into a com _ el table, the adm field and sat _ name field of the com _ el are recorded and substituted into a biu table in an ITU website, a status value with the same adm and sat _ name is searched, and when the status is C, the track N data is enabled;
The srs database is a frequency track resource database, geo is a geostationary space radio station table, long _ nom represents longitude of a space radio station, notice is a general data table for notification, ntc _ id represents a unique identifier of notice, and ntf _ rsn represents a notification reason; the N data is notification registration information which is provided by the satellite main part to the ITU and is used for the satellite network in brief and actually; c represents satellite network coordination data;
when the N data is started, ntc _ id is substituted into an s _ beam table, all freq _ min and freq _ max corresponding to ntc _ id are searched, freq _ min-freq _ max are substituted into a risk arc segment table, and whether the corresponding frequency is the same as the target frequency band or not is checked; wherein s _ beam is a satellite antenna beam table, freq _ min is the lowest frequency of a beam, freq _ max is the highest frequency of the beam, and frequency is a frequency band;
when the frequency ranges are the same, recording the satellite name and the orbit space of the orbit resource to be evaluated into an occupied resource table, and recording A 21 Number of actual orbiting satellites, A, to meet the conditions of occupied space resources 22 Is the minimum rail position interval; then, aiming at the orbit resource with the minimum orbit position interval value and the orbit resource to be evaluated, respectively selecting the wave beam with the maximum interference value from the uplink and the downlink by calculation, and respectively counting Calculating an interference value between an uplink beam and a downlink beam, and storing the interference value as an uplink interference value and a downlink interference value of occupied space resources, wherein the calculation process specifically comprises the following steps:
1. for the uplink, respectively substituting ntc _ id of the orbit resource into an s _ beam table, looking up beam _ name of ntc _ id with emi _ rcp being R, where R represents that the beam belongs to the uplink, and beam _ name is the satellite antenna beam identifier, and then respectively doing the following operation for each beam _ name:
(1.1) respectively substituting the beam _ names into a grp table, searching all grp _ ids under the beam _ names, then searching freq _ min and freq _ max corresponding to the grp _ ids, substituting freq _ min-freq _ max into a risk arc segment table, and judging whether the frequency band is the same as that of the track resource to be evaluated; substituting grp _ id with the same frequency band into an e _ as _ stn table, searching for gain under the grp _ id, comparing the gains of all grp _ id, and recording the maximum gain and the corresponding bmwdth; wherein, the grp table is an attribute table of the repeater group under the beam, grp _ id is a unique identifier of the repeater group, freq _ min is a lowest frequency value, freq _ max is a highest frequency value, e _ as _ stn table is an attribute table of the earth station corresponding to the repeater group, gain is gain of the repeater group, and bmwdth is the beam width of the earth station antenna;
(1.2) searching pwr _ ds _ max corresponding to grp _ id with the maximum gain in an emiss table, wherein the emiss table is a carrier attribute table, and pwr _ ds _ max is the maximum power density;
(1.3) judging whether the value of the track spacing 1.15 times is in the range of [1,100 λ/D ], λ represents the wavelength and D represents the radius, and calculating the gain recorded in the (1.1) th step and the pwr _ ds _ max recorded in the (1.2) th step as follows:
gain+pwr_ds_max+29-25log(A 22 *1.15)
wherein 100 λ/D — bmwdth 25/9, bmwdth being recorded in step (1.1);
comparing the values calculated by the formula in the step (1.3), selecting the beam with the maximum value to calculate the interference, and enabling A 23 Represents the uplink maximum interference value:
A 23 ＝max{I/N up ,I/N′ up }、
I/N up ＝G(θ) Rxsat +P TxES +G(φ) TxES -FSL up -T sat -B carrier -k
I/N′ up ＝G′(θ) Rxsat +P′ TxES +G′(φ) TxES -FSL′ up -T′ sat -B′ carrier -k
where I represents the total interference power and N up Representing the noise power, N ', produced by the rail resource with the minimum rail bit interval value' up Representing the noise power generated by the frequency track resource to be evaluated;
G(θ) Rxsat representing the reception gain, P, of the frequency-orbital resources to be evaluated in the direction of the earth station corresponding to the frequency-orbital resources having the smallest value of the orbital spacing TxES Indicating the transmission power of the frequency-orbit resource corresponding to the earth station with the minimum orbit spacing value, G (phi) TxES Indicating the side lobe gain, T, of the frequency-orbital resource corresponding to the earth station with the minimum orbital spacing value sat 、B carrier And FSL up Respectively representing the noise temperature, the carrier bandwidth and the free space loss of an uplink of a frequency-track resource system to be evaluated, wherein K represents a Boltzmann constant, and K is-228.6 dBJ/K; g' (theta) Rxsat Representing the receiving gain, P ', of the frequency-orbit resource with the minimum orbit bit interval value in the direction of the corresponding earth station of the frequency-orbit resource to be evaluated' TxES Indicating the transmission power, G' (phi), of the frequency-orbital resource to be evaluated relative to the earth station TxES Representing the side lobe gain, T 'of the earth station corresponding to the frequency track resource to be evaluated' sat 、B′ carrier And FSL' up Respectively representing the noise temperature, the carrier bandwidth and the free space loss of an uplink of the frequency-rail resource system with the minimum rail position interval value;
2. for the downlink, substituting ntc _ id of the track resource into the s _ beam table, looking up beam _ name of E at emi _ rcp under ntc _ id, where E indicates that the beam belongs to the downlink, and then performing the following operation on each beam _ name respectively:
(2.1) recording gain corresponding to the beam _ name;
(2.2) substituting the beam _ name into the grp table, searching all corresponding grp _ ids, substituting all grp _ ids into the emiss table, and searching the corresponding maximum pwr _ ds _ max;
(2.3) adding the gain recorded in the step (2.1) and the pwr _ ds _ max recorded in the step (2.2);
comparing the values calculated in step (2.3), selecting the beam with the maximum value to calculate interference, and enabling A 24 Represents the downlink maximum interference value:
A 24 ＝max{I/N down ,I/N′ down }
I/N down ＝G(α) RxES +P sat +G(φ) Txsat -FSL down -T sat -B carrier -k
I/N′ down ＝G′(α) RxES +P′ sat +G′(φ) Txsat -FSL′ down -T′ sat -B′ carrier -k
where I represents the total interference power and N up Representing the noise power, N ', produced by the rail resource with the minimum rail bit interval value' up Representing the noise power generated by the frequency track resource to be evaluated;
G(α) RxES represents the side lobe gain, P, of the earth station corresponding to the frequency-orbit resource to be evaluated in the direction of the frequency-orbit resource with the minimum orbit position interval value sat Representing the transmitted power of the frequency-rail resource with the minimum rail position interval value, G (phi) Txsat Representing antenna gain, T, of frequency-rail resource with minimum rail-bit spacing value sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of a downlink of a frequency-track resource system to be evaluated, wherein K represents a boltzmann constant, and K is-228.6 dBJ/K and has the unit of dBJ/K; g' (alpha) RxES Representing the side lobe gain, P 'of the earth station corresponding to the minimum orbit bit spacing value in the direction of the frequency-orbit resource to be evaluated' sat Representing the transmission power of the frequency-rail resource to be evaluated, G' (phi) Txsat Antenna gain, T, representing the frequency-track resource to be evaluated sat 、B carrier And FSL down Respectively representing the noise temperature, the carrier bandwidth and the free space loss of the downlink of the frequency-track resource with the minimum track bit interval value.
5. The GSO rail position effectiveness evaluation method according to claim 4, wherein the step 4 establishes a resource table to be coordinated for the rail position and frequency band information of the frequency rail resource to be evaluated, specifically as follows:
Looking up srs the long _ nom field of the geo table of the database, and when the long _ nom is in the evaluation range, recording the ntc _ id corresponding to the moment; substituting ntc _ id into a node table to search a ntf _ rsn value corresponding to the ntc _ id, wherein when ntf _ rsn is C, C indicates that the frequency-track resource has applied for coordination data, substituting ntc _ id of the C data frequency band into a com _ el table, wherein the com _ el table is a public content table, searching and recording sat _ name corresponding to the ntc _ id, searching the sat _ name again in the com _ el table, and when sat _ name is the same and ntf _ rsn is N, skipping the record and searching other frequency-track resources in an evaluation range; when the fact that the frequency bands have the same sat _ name and ntf _ rsn are N is not found, ntc _ id is substituted into an s _ beam table, all frequency _ min-frequency _ max corresponding to ntc _ id are found, frequency _ min-frequency _ max are substituted into a risk arc segment table, and whether the corresponding frequency is the same as the target frequency band or not is checked; when the frequency ranges are the same, recording the satellite name and the rail position interval of the orbit resource to be evaluated into an occupied resource table, and enabling A to be 31 Representing the number of actual on-orbit satellites eligible for the resource to be coordinated, A 32 Representing a minimum rail spacing;
then, aiming at the track resource with the minimum track position interval value and the track resource to be evaluated, selecting the beam with the maximum interference value from the uplink and the downlink respectively, calculating the interference value between the two beams of the uplink and the downlink, storing the interference value as the uplink interference value and the downlink interference value of the resource to be coordinated, wherein the calculation process is the same as the process of calculating the interference value in the step 3, and enabling A to be 33 Represents the maximum interference value of the uplink, A 34 Representing the downlink maximum interference value.
6. The GSO (ground satellite based orbit performance) performance evaluation method based on the multi-level joint risk of the orbit position as claimed in any one of claims 1 to 5, wherein the step 5 grades the index to be evaluated to make the risk degree A of the on-orbit satellite 1 Risk degree of occupied space resource A 2 Risk degree of resource to be coordinated A 3 Giving a decision matrix of the first-level index, giving a decision matrix corresponding to the second-level index for each first-level index, and recording the decision matrix into a database, wherein the specific steps are as follows:
let the risk degree A of the on-orbit satellite 1 The number A of the actual in-orbit satellites is a first-level index 11 Minimum orbital spacing of actual on-orbit satellites A 12 Is A 1 The secondary index of (1);
let the risk degree A of occupied space resource 2 The number A of the actual orbit satellites meeting the occupied space resource condition is a first-level index 21 Minimum rail spacing A 22 Maximum interference value of uplink A 23 And a downlink maximum interference value A 24 Is A 2 The secondary index of (1);
let the resource risk degree A to be coordinated 3 The actual number A of the satellites in orbit which meet the conditions of the resource to be coordinated is a first-level index 31 Minimum rail spacing A 32 Maximum interference value of uplink A 33 And a downlink maximum interference value A 34 Is A 2 The secondary index of (1);
the decision matrix of the first-level index is as shown in table 2:
TABLE 2
A 1 The decision matrix of the next two-level index is shown in table 3:
TABLE 3
A 2 The decision matrix of the next two-level index is shown in table 4:
TABLE 4
A 3 The decision matrix of the next two-level index is shown in table 5:
TABLE 5
7. The method for evaluating GSO rail position effectiveness based on rail position multi-level joint risk according to any one of claim 6, wherein the weights are calculated for the decision matrix of the first-level index and the decision matrix of the second-level index under the first-level index in step 6 respectively as follows:
for A 1 The decision matrix of the next two-level index, the value of each column is normalized to omega ij ，
Sum by row, ω 1 ＝ω 11 +ω 12 ，ω 2 ＝ω 21 +ω 22 Will be ω i Normalized to obtain
To obtain A 1 Weight vector of the next two-level index
By the same token, respectively obtain A 2 Weight vector of the next two-level index
A 3 Weight vector of the next two-level index
Weight vector of primary index decision matrix
8. The GSO rail level performance assessment method according to claim 7, wherein said step 7 is performed for a first level index A 1 、A 2 And A 3 Calculating the fuzzy satisfaction degree of each secondary index, setting 5 grades, namely A, B, C, D and E, and giving boundary values of the grades by experts, wherein the boundary values are as follows:
A 11 And A 21 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 31 Has a boundary value of [1,4,6,8,10,12,14,18 ]]，A 12 、A 22 And A 32 Has a boundary value of [0.5,1,1.5,2,2.5,3,3.5,4 ]]，A 23 And A 24 Has a boundary value of [1,2,3,4,5,6,7,9 ]]，A 33 And A 34 Has a boundary value of [1,4,6,8,10,12,14,18 ]]。
9. The rail position multistage joint risk based GSO rail position performance evaluation method of claim 8, wherein in step 8, ambiguity vectors are calculated for the values obtained in steps 2 to 4, respectively, to generate evaluation matrices, specifically as follows:
1. for A 11 、A 21 、A 23 、A 24 、A 31 、A 33 And A 34 Respectively setting i as 1,2,3,4 and 5, and sequentially substituting the following formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, wherein the evaluation matrix is r 11 、r 21 、r 23 、r 24 、r 31 、r 33 、r 34 ；
2. For A 12 、A 22 、A 32 Respectively setting i as 1,2,3,4 and 5, sequentially substituting formula (1), formula (2) and formula (3) to obtain an evaluation matrix with one row and five columns, and then respectively subtracting the values in the matrix from 1 to obtain a new evaluation matrix, wherein the evaluation matrix is r 12 、r 22 And r 32 ；
10. The GSO rail performance evaluation method based on rail level multistage joint risk according to claim 9, wherein the step 9 of calculating the ambiguity vector weighted by the secondary index based on the secondary index weight comprises the following steps:
summarizing the evaluation matrixes in the step 8, and respectively ordering the matrixes R 1 、R 2 And R 3 Is an index A 1 、A 2 And A 3 In which R is 1 ＝[r 11 ,r 12 ]，R 2 ＝[r 21 ,r 22 ,r 23 ,r 24 ]，R 3 ＝[r 31 ,r 32 ,r 33 ,r 34 ](ii) a Let beta 1 、β 2 And beta 3 In the form of a vector of degrees of ambiguity,
let matrix B ═ β 1 ,β 2 ,β 3 ) Wherein
And
are all obtained in step 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110324683.6A CN113162676B (en) | 2021-03-26 | 2021-03-26 | GSO rail position efficiency evaluation method based on rail position multistage joint risk |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110324683.6A CN113162676B (en) | 2021-03-26 | 2021-03-26 | GSO rail position efficiency evaluation method based on rail position multistage joint risk |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113162676A CN113162676A (en) | 2021-07-23 |
| CN113162676B true CN113162676B (en) | 2022-07-29 |
Family
ID=76884905
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110324683.6A Active CN113162676B (en) | 2021-03-26 | 2021-03-26 | GSO rail position efficiency evaluation method based on rail position multistage joint risk |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113162676B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113627799B (en) * | 2021-08-12 | 2023-09-15 | 天津(滨海)人工智能军民融合创新中心 | GSO frequency rail resource efficiency evaluation method based on military benefit |
| CN113988650A (en) * | 2021-10-30 | 2022-01-28 | 天津(滨海)人工智能军民融合创新中心 | Space frequency track resource comprehensive benefit evaluation method |
| CN114002717B (en) * | 2021-10-30 | 2024-08-09 | 天津(滨海)人工智能军民融合创新中心 | GSO rail position selection method based on operational risk |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010021656A2 (en) * | 2008-08-19 | 2010-02-25 | Trimble Navigation Limited | Gnss signal processing methods and apparatus with ambiguity selection |
| CN104504482A (en) * | 2014-10-24 | 2015-04-08 | 南京理工大学 | Fuzzy comprehensive assessment method for railway snow disaster risk |
| CN108629511A (en) * | 2018-05-04 | 2018-10-09 | 上海微小卫星工程中心 | A kind of satellite synthetic effectiveness evaluation method based on multifactor fuzzy theory reasoning and Analytic hierarchy process |
| CN109858693A (en) * | 2019-01-25 | 2019-06-07 | 中国科学院国家空间科学中心 | A kind of prediction technique for declaring situation towards satellite network data |
| CN111222776A (en) * | 2019-12-31 | 2020-06-02 | 中国科学院国家空间科学中心 | Satellite network coordination situation assessment method and system based on convolutional neural network |
| CN112149958A (en) * | 2020-08-24 | 2020-12-29 | 成都天奥集团有限公司 | Method for calculating multi-dimensional efficiency evaluation index of space-ground base station network resources |
-
2021
- 2021-03-26 CN CN202110324683.6A patent/CN113162676B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010021656A2 (en) * | 2008-08-19 | 2010-02-25 | Trimble Navigation Limited | Gnss signal processing methods and apparatus with ambiguity selection |
| CN104504482A (en) * | 2014-10-24 | 2015-04-08 | 南京理工大学 | Fuzzy comprehensive assessment method for railway snow disaster risk |
| CN108629511A (en) * | 2018-05-04 | 2018-10-09 | 上海微小卫星工程中心 | A kind of satellite synthetic effectiveness evaluation method based on multifactor fuzzy theory reasoning and Analytic hierarchy process |
| CN109858693A (en) * | 2019-01-25 | 2019-06-07 | 中国科学院国家空间科学中心 | A kind of prediction technique for declaring situation towards satellite network data |
| CN111222776A (en) * | 2019-12-31 | 2020-06-02 | 中国科学院国家空间科学中心 | Satellite network coordination situation assessment method and system based on convolutional neural network |
| CN112149958A (en) * | 2020-08-24 | 2020-12-29 | 成都天奥集团有限公司 | Method for calculating multi-dimensional efficiency evaluation index of space-ground base station network resources |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113162676A (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113162676B (en) | GSO rail position efficiency evaluation method based on rail position multistage joint risk | |
| CN107330589B (en) | Quantitative evaluation method and system for satellite network coordination risk | |
| CN102523032B (en) | Method for determining orbit position of geosynchronous earth orbit satellite | |
| CN110147870B (en) | Optimization method for task allocation of multiple unmanned aerial vehicles after earthquake disaster | |
| CN106534392A (en) | Positioning information acquiring method, positioning method and apparatus | |
| CN107770783B (en) | Base station capacity expansion transformation scheme design method and related equipment | |
| He et al. | Random forests based path loss prediction in mobile communication systems | |
| CN109995440A (en) | Aerial angle method for adjusting combined, device, equipment and medium | |
| CN112149949B (en) | Satellite communication system efficiency evaluation method based on hesitation blurring and analytic hierarchy process | |
| CN113627799B (en) | GSO frequency rail resource efficiency evaluation method based on military benefit | |
| CN114399652B (en) | Multi-index system GEO satellite orbit selection method based on radar chart intelligent matching | |
| CN113610363B (en) | GSO frequency rail resource efficiency evaluation method based on rail position input use risk | |
| CN113610360B (en) | GSO frequency rail resource efficiency evaluation method based on strategic value | |
| CN113973336B (en) | Method, device, equipment and storage medium for determining interference cells in network | |
| Siderius | Proceedings of the Workshop on Land Evaluation for Land-use Planning and Conservation in Sloping Areas, Enschede, 17-21 December 1984 | |
| CN114002717B (en) | GSO rail position selection method based on operational risk | |
| CN115390113A (en) | BD3 high-precision positioning method and system for power system | |
| CN114358460A (en) | GSO frequency-rail resource efficiency evaluation method based on demand constraint | |
| KR101593876B1 (en) | Method and system for assessing environmental impacts capable of establishing weighting factors for environmental impact categories using internet search statistics | |
| CN114826370B (en) | Rail position selection method of GSO satellite | |
| Sun et al. | An efficient local search algorithm for minimum positive influence dominating set problem | |
| CN105224687A (en) | A kind of low consumed Apriori method | |
| Lomakin et al. | Modeling Single-Frequency DVB-T2 Network | |
| CN114665950B (en) | Method for calculating support degree of space spectrum resource to satellite communication system uplink | |
| CN113610361A (en) | GSO frequency-rail resource efficiency evaluation method based on coordination difficulty |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |