CN114630362A - Method for calculating support degree of space spectrum resources to satellite communication system downlink - Google Patents

Method for calculating support degree of space spectrum resources to satellite communication system downlink Download PDF

Info

Publication number
CN114630362A
CN114630362A CN202210262381.5A CN202210262381A CN114630362A CN 114630362 A CN114630362 A CN 114630362A CN 202210262381 A CN202210262381 A CN 202210262381A CN 114630362 A CN114630362 A CN 114630362A
Authority
CN
China
Prior art keywords
earth station
satellite
grid point
point earth
downlink
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.)
Granted
Application number
CN202210262381.5A
Other languages
Chinese (zh)
Other versions
CN114630362B (en
Inventor
王利利
张周
于思源
万钰涵
胡向晖
桑玮
董飞鸿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Network Engineering Institute of Systems Engineering Academy of Military Sciences
Original Assignee
Institute of Network Engineering Institute of Systems Engineering Academy of Military Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Network Engineering Institute of Systems Engineering Academy of Military Sciences filed Critical Institute of Network Engineering Institute of Systems Engineering Academy of Military Sciences
Priority to CN202210262381.5A priority Critical patent/CN114630362B/en
Publication of CN114630362A publication Critical patent/CN114630362A/en
Application granted granted Critical
Publication of CN114630362B publication Critical patent/CN114630362B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/08Testing, supervising or monitoring using real traffic
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W84/00Network topologies
    • H04W84/02Hierarchically pre-organised networks, e.g. paging networks, cellular networks, WLAN [Wireless Local Area Network] or WLL [Wireless Local Loop]
    • H04W84/04Large scale networks; Deep hierarchical networks
    • H04W84/06Airborne or Satellite Networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Radio Relay Systems (AREA)

Abstract

The invention discloses a method for calculating the support of space spectrum resources to a satellite communication system downlink, which comprises the following steps: acquiring related calculation parameter information; searching a wave beam which meets the target frequency band condition in the space frequency spectrum resource to be calculated; dividing an earth area grid needing to be subjected to support degree calculation into K grid points; respectively judging whether each grid point earth station meets the lowest communication requirement of a satellite system downlink by using the transmitting power spectral density and the demodulation threshold required by a user; and counting the number of the grid point earth stations meeting the minimum communication requirement of the downlink of the satellite communication system, and calculating the support degree of the downlink. The method decomposes the target earth region by using a grid method, the calculation of the support degree is more comprehensive and credible, and a user can select to use the transmitting power spectral density to calculate the support degree when needing rough support degree, and select to use the demodulation threshold to calculate if needing precise support degree. The method is simple and has better project realizability.

Description

Method for calculating support degree of space spectrum resources to satellite communication system downlink
Technical Field
The invention relates to the technical field of wireless communication, in particular to a method for calculating the support of space spectrum resources to a satellite communication system downlink.
Background
The spatial spectrum resource refers to satellite frequency and orbit resource and is an essential element for the development of satellite systems. With the development of satellite communication technology, the competition of world countries for space spectrum resources is more and more in a trend of becoming white. When the space spectrum resources are utilized, the value of the space spectrum resources to the satellite system needs to be evaluated to judge whether the resources can be effectively applied to the satellite system or not so as to achieve the purpose of efficiently utilizing the space spectrum resources. The downlink support degree is the ratio of the number of information which can be successfully transmitted to the satellite system by the earth station in the satellite communication system, and in the actual calculation, a certain region on the earth needs to be specified, the region needs to be divided into grid points, and the support degree between the grid points and the satellite system for space spectrum resource application needs to be calculated respectively. The existing method for calculating the support of the space spectrum resource downlink is relatively deficient, so that an existing implementation scheme is not available at present, and a complete method for calculating the support of the space spectrum resource to the satellite communication system downlink is necessary to be researched and formulated in order to enhance the reserve of frequency-orbit resources for a satellite system in China.
Disclosure of Invention
The invention discloses a method for calculating the support of space spectrum resources to a satellite communication system downlink, which aims to solve the problem that the conventional method for calculating the support of the space spectrum resources to the satellite communication system downlink is relatively deficient, can quantitatively calculate the support of the space spectrum resources to the satellite system downlink, and can calculate whether the space spectrum resources meet the minimum communication requirement of the satellite system downlink or not from longitude and latitude areas. Downlink support refers to the number of earth stations in a satellite communication system that can successfully transmit information to the satellite system as a proportion of the total number of earth stations.
The invention discloses a method for calculating the support of space spectrum resources to a satellite communication system downlink, which comprises the following steps:
s1, acquiring related calculation parameter information; the method specifically comprises a space spectrum resource identifier ntc _ id required to be calculated, an earth region required to be subjected to support degree calculation, a target frequency band freq _ need, a selection identifier method _ flag of a calculation method, a transmitting power spectral density psd required by a user, and a demodulation threshold required by the user;
s2, searching beams which meet the target frequency band condition in the space spectrum resources to be calculated;
and acquiring all beams under the space spectrum resource according to the space spectrum resource identifier ntc _ id, sequentially reading the frequency range of each beam, comparing the frequency range with the target frequency band freq _ need, if the frequency range of each beam is overlapped with the target frequency band freq _ need, storing the beam information into an array freq _ list, and taking the beam as the beam meeting the target frequency band condition in the space spectrum resource.
S3, dividing the region grid of the earth needing to be subjected to support degree calculation into K grid points; according to the longitude and latitude of the area region on the map, the area region is evenly divided into K grid points, each grid point is provided with an earth station of a satellite communication system, the earth station is called a grid point earth station, and each grid point earth station is identified by a longitude and latitude sequence; judging a selection identifier method _ flag of the calculation method, if the method _ flag is 0, judging the grid point earth station by using the transmitting power spectral density psd required by a user, and turning to step S4; if the method _ flag is 1, the demodulation threshold required by the user is used to determine the grid point earth station, and the process proceeds to step S5.
S4, judging whether each mesh point earth station meets the lowest communication requirement of the satellite system downlink by using the transmitting power spectrum density required by the user, and turning to the step S6;
according to the K grid point earth stations obtained in the step S3, sequentially testing whether each grid point earth station meets the lowest communication requirement of a downlink, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the identification ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
s41, judging whether the mesh point earth station meets the lowest communication requirement of the downlink, detecting the mark ok _ flag field, if ok _ flag is 1, indicating that the mesh point earth station meets the lowest communication requirement of the downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating the beam; if ok _ flag is 0, the process proceeds to step S42.
S42, calculating the path distance d between the satellite and the grid point earth station in the satellite communication system, wherein the calculation formula is as follows:
Figure BDA0003550554850000031
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. And acquiring the longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id.
And S43, calculating the beam center frequency f. According to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2;
s44, calculating path loss LF. Obtaining path loss L from the earth station of the grid point to the satellite according to the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S45, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity.
And S46, calculating the satellite transmission gain G. And acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as the transmission gain G of the satellite.
And S47, calculating the Effective Isotropic Radiated Power (EIRP) of the grid point earth station. According to the path loss LFThe effective omnidirectional radiation power EIRP is obtained by the satellite transmitting power P and the satellite gain G, and the calculation formula is as follows:
EIRP=P+G-LF
and S48, detecting whether the earth station of the grid point reaches the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the transmitted power spectral density psd required by a user, if the EIRP is more than or equal to the psd, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink under the beam, and setting ok _ flag to 1; if EIRP < psd, the earth station of the mesh point can not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S5, calculating whether each grid point earth station meets the lowest communication requirement of the satellite communication system downlink by using the demodulation threshold required by the user, and turning to the step S6;
sequentially calculating whether each grid point earth station meets the lowest communication requirement of a downlink of a satellite communication system according to the K grid point earth stations obtained in the step S3, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
and S51, judging whether the mesh point earth station meets the minimum communication requirement of the satellite communication system downlink. Detecting an identification ok _ flag field of the mesh point earth station, if the ok _ flag field is 1, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating a beam of the mesh point earth station; if ok _ flag is 0, proceeding to step S52;
and S52, calculating a path d between the satellite and the earth station of the grid point in the satellite system. The calculation formula is as follows:
Figure BDA0003550554850000041
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. Acquiring longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id;
and S53, calculating the beam center frequency f. The beam center frequency is calculated according to the frequency range { freq _ min, freq _ max } of the beam, which is calculated by the formula:
f=(freq_min+freq_max)/2,
s54, calculating path loss LF. Obtaining path loss L from the mesh point earth station to the satellite communication system from the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S55, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity.
And S56, calculating the grid point earth station receiving gain. And calculating the contour line coArea of the projection range of the beam according to the space frequency spectrum resource identifier ntc _ id and the projection range of the beam on the earth, and obtaining the receiving gain of the grid point earth station according to the distance between the central point of the projection range and the grid point earth station.
And S57, calculating the satellite receiving gain G. And acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as satellite receiving gain G.
S58, calculating the receiving power P of the grid point earth stationr. According to the path loss LFSatellite transmit power P and a grid point earth station receive gain,obtaining grid point earth station received power PrThe calculation formula is as follows:
Pr=G+P-LF+gain;
and S59, detecting whether the earth station of the grid point reaches the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the demodulation threshold required by the user, if P isrIf the wave beam is greater than or equal to threshold, the earth station of the grid point is explained to reach the lowest communication requirement of a satellite communication system downlink under the wave beam, and ok _ flag is set to be 1; if PrIf < threshold, it means that the earth station of the grid point does not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S6, counting the number of grid point earth stations meeting the lowest communication requirement of the satellite communication system downlink, and calculating the downlink support degree; and circularly detecting ok _ flag fields of the K grid point earth stations, and setting a num variable to record the number of the grid point earth stations which reach the minimum communication requirement of a satellite communication system downlink. If the ok _ flag field of the grid point earth station is 1, adding 1 to the num variable; num is unchanged when the ok _ flag field of the mesh point earth station is 0. Taking the ratio of the number of the grid point earth stations meeting the minimum communication requirement to the total number of the grid point earth stations as the support degree support of the space spectrum resource to the downlink of the satellite communication systemupThe calculation formula is as follows:
Figure BDA0003550554850000051
wherein K is the total number of earth stations of the grid point.
The invention has the beneficial effects that:
(1) according to the method, the target earth region is decomposed by using a grid method, and the calculation of the support degree is more comprehensive and credible; (2) the calculation method of the invention is optional, when a user needs rough support, the user can choose to use the transmitting power spectral density to calculate the support, and if the user needs precise support, the user chooses to use the demodulation threshold to calculate; (3) the method is simple, has strong operability and better project realizability.
Drawings
FIG. 1 is a flow chart of an implementation of the method of the present invention;
FIG. 2 is a schematic diagram of a user calibrating a target area according to the present invention;
FIG. 3 is a schematic diagram of a segmentation grid according to the present invention.
Detailed Description
For a better understanding of the present disclosure, three examples are given herein.
FIG. 1 is a flow chart of an implementation of the method of the present invention; FIG. 2 is a schematic diagram of a user calibrating a target area according to the present invention; FIG. 3 is a schematic diagram of a segmentation grid according to the present invention.
The first embodiment is as follows:
the invention discloses a method for calculating the support of space spectrum resources to a satellite communication system downlink, which comprises the following steps:
s1, acquiring related calculation parameter information; the method specifically comprises a space spectrum resource identifier ntc _ id required to be calculated, an earth region required to be subjected to support degree calculation, a target frequency band freq _ need, a selection identifier method _ flag of a calculation method, a transmitting power spectral density psd required by a user, and a demodulation threshold required by the user;
s2, searching beams which meet the target frequency band condition in the space spectrum resources to be calculated;
and acquiring all beams under the space spectrum resource according to the space spectrum resource identifier ntc _ id, sequentially reading the frequency range of each beam, comparing the frequency range with the target frequency band freq _ need, if the frequency range of each beam is overlapped with the target frequency band freq _ need, storing the beam information into an array freq _ list, and taking the beam as the beam meeting the target frequency band condition in the space spectrum resource.
S3, dividing the region grid of the earth needing to be subjected to support degree calculation into K grid points; according to the longitude and latitude of the area region on the map, the area region is evenly divided into K grid points, each grid point is provided with an earth station of a satellite communication system, the earth station is called a grid point earth station, and each grid point earth station is identified by a longitude and latitude sequence; judging a selection identifier method _ flag of the calculation method, if the method _ flag is 0, judging the grid point earth station by using the transmitting power spectral density psd required by a user, and turning to step S4; if the method _ flag is 1, the demodulation threshold required by the user is used to determine the grid point earth station, and the process proceeds to step S5.
S4, judging whether each mesh point earth station meets the lowest communication requirement of the satellite system downlink by using the transmitting power spectrum density required by the user, and turning to the step S6;
according to the K grid point earth stations obtained in the step S3, sequentially testing whether each grid point earth station meets the lowest communication requirement of a downlink, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the identification ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
s41, judging whether the mesh point earth station meets the lowest communication requirement of the downlink, detecting the mark ok _ flag field, if ok _ flag is 1, indicating that the mesh point earth station meets the lowest communication requirement of the downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating the beam; if ok _ flag is 0, the process proceeds to step S42.
S42, calculating the path distance d between the satellite and the grid point earth station in the satellite communication system, wherein the calculation formula is as follows:
Figure BDA0003550554850000071
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. And acquiring the longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id.
And S43, calculating the beam center frequency f. According to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2;
s44, calculating the path loss LF. Obtaining path loss L from the earth station of the grid point to the satellite according to the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S45, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity.
And S46, calculating the satellite gain G. And acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as the transmission gain G of the satellite.
And S47, calculating the Effective Isotropic Radiated Power (EIRP) of the grid point earth station. According to the path loss LFThe effective omnidirectional radiation power EIRP is obtained by the satellite transmitting power P and the satellite gain G, and the calculation formula is as follows:
EIRP=P+G-LF
and S48, detecting whether the earth station of the mesh point meets the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the transmitting power spectral density psd required by a user, if the EIRP is more than or equal to the psd, indicating that the earth station of the mesh point meets the lowest communication requirement of a satellite communication system downlink under the wave beam, and setting ok _ flag to be 1; if EIRP < psd, the mesh point earth station can not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S5, calculating whether each grid point earth station meets the lowest communication requirement of the satellite communication system downlink by using the demodulation threshold required by the user, and turning to the step S6;
sequentially calculating whether each grid point earth station meets the lowest communication requirement of a downlink of a satellite communication system according to the K grid point earth stations obtained in the step S3, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
and S51, judging whether the mesh point earth station meets the minimum communication requirement of the satellite communication system downlink. Detecting an identification ok _ flag field of the mesh point earth station, if the ok _ flag field is 1, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating a beam of the mesh point earth station; if ok _ flag is 0, proceeding to step S52;
and S52, calculating a path d between the satellite and the earth station of the grid point in the satellite system. The calculation formula is as follows:
Figure BDA0003550554850000091
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. Acquiring longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id;
and S53, calculating the beam center frequency f. The beam center frequency is calculated according to the frequency range { freq _ min, freq _ max } of the beam, which is calculated by the formula:
f=(freq_min+freq_max)/2,
s54, calculating path loss LF. Obtaining path loss L from the mesh point earth station to the satellite communication system from the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S55, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity.
And S56, calculating the grid point earth station receiving gain. And calculating the contour line coArea of the projection range of the beam according to the space frequency spectrum resource identifier ntc _ id and the projection range of the beam on the earth, and obtaining the receiving gain of the grid point earth station according to the distance between the central point of the projection range and the grid point earth station.
And S57, calculating the satellite receiving gain G. And acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as satellite receiving gain G.
S58, calculating the receiving power P of the grid point earth stationr. According to the path loss LFThe satellite transmitting power P and the grid point earth station receiving gain are obtained to obtain the grid point earth station receiving power PrThe calculation formula is as follows:
Pr=G+P-LF+gain;
and S59, detecting whether the earth station of the grid point reaches the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the demodulation threshold required by the user, if P isrIf the wave beam is greater than or equal to threshold, the earth station of the grid point meets the lowest communication requirement of a downlink of a satellite communication system under the wave beam, and ok _ flag is set to be 1; if PrIf < threshold, it means that the earth station of the grid point does not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S6, counting the number of grid point earth stations meeting the lowest communication requirement of the satellite communication system downlink, and calculating the downlink support degree; and circularly detecting ok _ flag fields of the K grid point earth stations, and setting a num variable to record the number of the grid point earth stations which reach the minimum communication requirement of a satellite communication system downlink. If the ok _ flag field of the grid point earth station is 1, adding 1 to the num variable; num is unchanged when the ok _ flag field of the grid point earth station is 0. Taking the ratio of the number of the grid point earth stations meeting the minimum communication requirement to the total number of the grid point earth stations as the support degree support of the space spectrum resource to the downlink of the satellite communication systemupThe calculation formula is as follows:
Figure BDA0003550554850000101
wherein K is the total number of earth stations of the grid point.
Example two:
the invention discloses a method for calculating the support of space spectrum resources to a satellite communication system downlink, which comprises the following steps:
s1, acquiring related calculation parameter information, wherein the related calculation parameter information specifically comprises a space spectrum resource identifier ntc _ id required to be calculated, an earth region required to be subjected to support degree calculation, a target frequency band freq _ need, a selection identifier method _ flag of a calculation method, a transmitting power spectrum density psd required by a user and a demodulation threshold required by the user;
s2, searching beams which meet the target frequency band condition in the space spectrum resources to be calculated;
and bringing the spatial frequency spectrum resource identifier ntc _ id into an s _ beam table in a srs database, acquiring all beam beams _ name under the resource according to the spatial frequency spectrum resource identifier ntc _ id, sequentially reading the frequency range { freq _ min, freq _ max } (MHz) of each beam, comparing the frequency range { freq _ min, freq _ max } (MHz) with a target frequency band freq _ need, if the frequency bands of the two overlap, storing the beam information into an array freq _ list, and taking the beam as a beam meeting the target frequency band condition in the spatial frequency spectrum resource. srs the database includes the satellite network performance parameters currently in operation.
S3, dividing the region grid of the earth which needs to be subjected to support calculation into K grid points, calculating the segmentation grid of the external matrix according to the longitude and latitude of the region on the map, uniformly dividing the region into K grid points, wherein each grid point is provided with an earth station of a satellite communication system, called a grid point earth station, and each grid point earth station is identified by a longitude and latitude sequence; judging a selection identifier method _ flag of the calculation method, if the method _ flag is 0, judging the grid point earth station by using the transmitting power spectral density psd required by a user, and turning to step S4; if the method _ flag is 1, the demodulation threshold required by the user is used to determine the grid point earth station, and the process proceeds to step S5.
S4, judging whether each grid point earth station meets the lowest communication requirement of the satellite system downlink by using the transmitting power spectral density required by the user, and turning to the step S6;
according to the K grid point earth stations obtained in the step S3, testing whether each grid point earth station meets the lowest communication requirement of a downlink in sequence, setting a mark ok _ flag field in each grid point earth station to show whether the grid point earth station meets the requirement, and initializing the mark ok _ flag field to be 0; for each mesh point earth station, performing the following operations on each beam meeting the target frequency band condition in each spatial frequency spectrum resource in freq _ list:
s41, judging whether the mesh point earth station meets the lowest communication requirement of the downlink, detecting the mark ok _ flag field, if the ok _ flag is 1, indicating that the mesh point earth station meets the lowest communication requirement of the downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating the beam; if ok _ flag is 0, the flow proceeds to step S42.
S42, calculating the path distance d between the satellite and the grid point earth station in the satellite communication system, wherein the calculation formula is as follows:
Figure BDA0003550554850000111
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. And reading the bore _ long and bore _ lat fields of the s _ beam table by using the space spectrum resource identifier ntc _ id to acquire the longitude and latitude coordinates of the satellite. The unit of d is km.
And S43, calculating the beam center frequency f. According to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2;
s44, calculating path loss LF. Obtaining path loss L from the earth station of the grid point to the satellite according to the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S45, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id by using an emiss table according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity. The emiss table records the maximum transmit power of each satellite beam of the operating satellite network.
And S46, calculating the satellite transmission gain G. And acquiring a gain corresponding to the beam group identifier grp _ id by using the e _ as _ stn table according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as the transmission gain G of the satellite. The e _ as _ stn table records the transmission gain of each satellite of the operating satellite network.
And S47, calculating the Effective Isotropic Radiated Power (EIRP) of the grid point earth station. According to the path loss LFThe effective omnidirectional radiation power EIRP is obtained by the satellite transmission power P and the satellite transmission gain G, and the calculation formula is as follows:
EIRP=P+G-LF
and S48, detecting whether the earth station of the grid point reaches the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the transmitted power spectral density psd required by a user, if the EIRP is more than or equal to the psd, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink under the beam, and setting ok _ flag to 1; if EIRP < psd, the mesh point earth station can not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S5, calculating whether each grid point earth station reaches the lowest communication requirement of the satellite communication system downlink by using the demodulation threshold required by the user, and turning to the step S6;
sequentially calculating whether each grid point earth station meets the lowest communication requirement of a downlink of a satellite communication system according to the K grid point earth stations obtained in the step S3, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
and S51, judging whether the mesh point earth station meets the minimum communication requirement of the satellite communication system downlink. Detecting an identification ok _ flag field of the mesh point earth station, if the ok _ flag field is 1, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating a beam of the mesh point earth station; if ok _ flag is 0, proceeding to step S52;
and S52, calculating a path d between the satellite and the earth station of the grid point in the satellite system. The calculation formula is as follows:
Figure BDA0003550554850000131
where θ is the difference between the longitude of the satellite and the grid point earth station, and α is the grid point earth station latitude. Reading a bore _ long field and a bore _ lat field of an s _ beam table by using a space spectrum resource identifier ntc _ id to obtain a satellite longitude and latitude coordinate;
and S53, calculating the beam center frequency f. According to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2,
s54, calculating path loss LF. Obtaining path loss L from the mesh point earth station to the satellite communication system from the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
and S55, calculating the satellite transmission power P. And acquiring the maximum transmitting power pwr _ ds _ max of the beam group through the grp _ id by using an emiss table according to the beam group identification grp _ id corresponding to the beam, and taking the maximum transmitting power pwr _ ds _ max as the transmitting power P of the satellite. The beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity. The emiss table records the maximum transmit power of each satellite beam of the operating satellite network.
And S56, calculating the grid point earth station receiving gain. And calculating the contour line coArea of the projection range of the beam according to the space frequency spectrum resource identifier ntc _ id and the projection range of the beam on the earth, and obtaining the receiving gain of the grid point earth station according to the distance between the central point of the projection range and the grid point earth station.
And S57, calculating the satellite receiving gain G. All beam group identities grp _ id are obtained by bringing the spatial spectrum resource identities ntc _ id into the grp table in the srs database. And according to the beam group identification grp _ id corresponding to the beam, bringing the grp _ id into an e _ as _ stn table in an srs database, acquiring the corresponding gain of the grp _ id, and taking the gain as satellite receiving gain G. The e _ as _ stn table records the gain of each satellite of the network of satellites in operation.
S58, calculating the receiving power P of the grid point earth stationr. According to the path loss LFThe satellite transmitting power P and the grid point earth station gain are obtained to obtain the grid point earth station receiving power PrThe calculation formula is as follows:
Pr=G+P-LF+gain;
and S59, detecting whether the earth station of the grid point reaches the minimum communication requirement of the downlink of the satellite communication system under the beam. Judging according to the demodulation threshold required by the user, if P isrIf the wave beam is greater than or equal to threshold, the earth station of the grid point meets the lowest communication requirement of a downlink of a satellite communication system under the wave beam, and ok _ flag is set to be 1; if PrIf < threshold, it means that the earth station of the grid point does not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
S6, counting the number of grid point earth stations meeting the lowest communication requirement of the satellite communication system downlink, and calculating the downlink support degree; and circularly detecting ok _ flag fields of the K grid point earth stations, and setting a num variable to record the number of the grid point earth stations which reach the minimum communication requirement of a satellite communication system downlink. If the ok _ flag field of the grid point earth station is 1, the num variable is added with 1(ii) a Num is unchanged when the ok _ flag field of the grid point earth station is 0. Taking the ratio of the number of the grid point earth stations meeting the minimum communication requirement to the total number of the grid point earth stations as the support degree support of the space spectrum resource to the downlink of the satellite communication systemupThe calculation formula is as follows:
Figure BDA0003550554850000151
wherein K is the total number of grid points earth stations.
Example three:
network profile ntc _ id of track resource at spatial frequency: 111500204, the satellite network name is CHNSAT-130E.
Step 1, inputting parameters required by calculation.
1. Inputting the spatial spectrum resource identifier ntc _ id to be calculated: 111500204. the specific content of the target frequency band freq _ need is shown in table 1. The selection flag of the calculation method is 1, the transmitted power spectral density psd is 10, and the demodulation threshold is 18.
TABLE 1 target frequency band Table
freq_min freq_max
1980 2010
2025 2110
2170 2290
2500 2535
2655 2690
10700 11450
11700 13250
13750 14450
14500 14800
17300 21200
27000 31000
2. The earth region for which the support calculation is performed is input, and a target region of an arbitrary shape is marked on the map by the user and is referred to as a region, as shown in fig. 2.
And 2, searching beams meeting the target frequency band condition in the space spectrum resources to be calculated.
The spatial spectrum resource identifier ntc _ id is brought into an s _ beam table in a srs database, all beam names beam _ names under the resource are obtained, and { freq _ min, freq _ max } (MHz) of each beam _ name is read in sequence, as shown in table 2. And then comparing the frequency band with a target frequency band freq _ need input by a user, if an overlapped frequency band exists, storing the frequency band into an array freq _ list, and performing the next operation as a beam meeting the target frequency band condition in the spatial frequency spectrum resource, as shown in table 3.
Table 2 found beam sequence
id beam_name freq_min freq_max
1 CD 3420.0 4200.0
2 KAD 17700.0 21200.0
3 KDS 11460.0 12680.0
4 KUD 11460.0 12700.0
5 CU 5925.0 6665.0
6 KAU 27500.0 31000.0
7 KUS 13760.0 14500.0
8 KUU 13800.0 14500.0
TABLE 3 Beam sequence to meet user requirements
id beam_name freq_min freq_max
1 KAD 17700.0 21200.0
2 KDS 11460.0 12680.0
3 KUS 13760.0 14500.0
4 KUD 11460.0 12700.0
5 KUU 13800.0 14500.0
6 KAU 27500.0 31000.0
And 3, gridding the region of the earth needing to be subjected to support degree calculation into K grids.
1. According to a target area region calibrated on a map by a user, according to the longitude and latitude of the target area region on the map, a segmentation grid of an external matrix is calculated by using a graphical tool, the region is divided into 100 grid points, and each grid point is identified by a longitude and latitude sequence, as shown in fig. 3.
2. And judging a selection identifier method _ flag of the calculation method, wherein the method _ flag is 1, and then using a demodulation threshold to continue the judgment of the grid points.
And 4, calculating whether each grid point meets the minimum communication requirement of the satellite system downlink by using a demodulation threshold threesold.
And (3) testing whether each grid point meets the lowest communication requirement of a downlink according to the 100 grid points acquired in the step (3) in sequence, setting a mark ok _ flag field in each grid point to indicate whether the grid point meets the requirement, initializing the field to 0, and performing the following operation on each grid point.
For the grid points [ -178.6787,0.9793], for the beam cycle satisfying the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed, taking the beam KAD as an example:
1. it is determined whether the mesh point has reached a valid coverage requirement. And if the qualified flag ok _ flag is detected to be 0, continuing the following operation.
2. A path d of the satellite system to the grid point earth station is calculated. The calculation formula is as follows:
Figure BDA0003550554850000171
wherein
Figure BDA0003550554850000172
Is the difference between the longitude of the satellite and the grid point earth station and α is the grid point earth station latitude. The bore _ long of the s _ beam table is read by using the spatial spectrum resource identifier ntc _ id to obtain that the longitude of the satellite is 130 and the longitude of the earth station is-178.6787, so that the theta is 308.6787. The earth station latitude is 0.9793, from which a distance d of 40148.003(Km) is obtained.
3. The beam center frequency f is calculated. Freq _ min corresponding to beam _ name in freq _ list of beam: 17700 and freq _ max: 21200, i.e., the minimum frequency and the maximum frequency of the beam, the beam center frequency is calculated. The calculation formula is as follows:
f=(freq_min+freq_max)/(2×1000)
this results in a beam center frequency f of 19.45.
4. Calculating the path loss LF. From the path d and the beam center frequency f obtained above, the path loss L from the grid point earth station to the satellite system is obtainedF. The calculation formula is as follows:
LF=92.45+20lg(d×f)
thereby obtaining a path loss LFIs 210.301.
5. Calculating satellite transmission power Pi. All beam group identifications grp _ id under the beam are obtained by substituting space spectrum resource identifications ntc _ id into a grp table in a srs database, then the beam group identifications grp _ id are substituted into an emiss table in a srs database to obtain transmitting power pwr _ ds _ nbw, and the maximum transmitting power value under the beam is selected as satellite transmitting power PiIt was-34.1 dbm.
6. A grid point earth station receive gain is calculated. The contour of the beam, coArea, is calculated as "contour" [ [18.8,38.8], [19.00,40.00], [20.1,40.1], [20,39], [18.8,38.8] ], "gain":25}, { "contour" [ [17.8,37.8], [18.00,41.00], [21.1,41.1], [21,38], [17.8,37.8], [ gain ":5} ] by substituting ntc _ id of the spatial spectrum resource and beam _ name into the GIMs interface. The contour line coArea and the latitude and longitude coordinate [ -178.6787,0.9793] fields of the grid point earth stations are then substituted into the GIMs to obtain a gain value gain corresponding to the grid point earth stations of 8.93206882477.
7. Calculating satellite gain Gi. All beam group identities grp _ id are obtained by bringing the spatial spectrum resource identities ntc _ id into the grp table in the srs database. The maximum receiving antenna gain G corresponding to the satellite is obtained by being brought into an e _ as _ stn table in an srs database through grp _ idiWas 41.7.
8. Calculating grid point earth station received power Pr. From the path loss L obtained previouslyFTransmitting power PiAnd grid point geostationary gain GiThe received power P of the grid point earth station can be obtainedr. The calculation formula is as follows:
Pr=Gi+Pi-LF+gain
the received power at grid points-178.6787, 0.9793 is found to be-158.94.
9. It is detected whether the earth station at the mesh point meets the minimum communication requirements under the beam. Judging according to the demodulation threshold 18 input by the user, PrIf the grid point earth station is not reached yet, the grid point earth station is judged to be not reached yetAnd requesting to enter the next round of detection.
And step six, counting the number of the up-to-standard grid points in the grid point earth station, and calculating the support degree of the downlink. The ok flag of 100 mesh point earth stations is cyclically detected, and the num variable is set to record the number of mesh points which reach the downlink communication standard. If the ok _ flag of the grid point earth station is 1, adding one to the num variable; when ok _ flag is 0, num is unchanged, resulting in num of 36. Finally, the ratio of the number of the earth stations of the qualified grid point to the total number of the earth stations of the grid point is used as the support degree support for the downlinkdown. The calculation formula is as follows:
Figure BDA0003550554850000191
this gives a downlink support of 0.36.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims (8)

1. A method for calculating the support of spatial spectrum resources to a downlink of a satellite communication system, comprising the steps of:
s1, acquiring related calculation parameter information;
s2, searching beams which meet the target frequency band condition in the space spectrum resources to be calculated;
s3, dividing the region grid of the earth needing to be subjected to support degree calculation into K grid points; judging a selection identifier method _ flag of the calculation method, if the method _ flag is 0, judging the grid point earth station by using the transmitting power spectral density psd required by a user, and turning to step S4; if the method _ flag is 1, judging the grid point earth station by using a demodulation threshold required by a user, and turning to the step S5;
s4, judging whether each grid point earth station meets the lowest communication requirement of the satellite system downlink by using the transmitting power spectral density required by the user, and turning to the step S6;
s5, calculating whether each grid point earth station meets the lowest communication requirement of the satellite communication system downlink by using the demodulation threshold required by the user, and turning to the step S6;
and S6, counting the number of the mesh point earth stations meeting the minimum communication requirement of the satellite communication system downlink, and calculating the downlink support degree.
2. The method for calculating the downlink support of spatial spectrum resources in the satellite communication system according to claim 1, wherein the calculation parameter information in step S1 specifically includes a spatial spectrum resource identifier ntc _ id that needs to be calculated, an earth region that needs to be calculated for the support, a target frequency band freq _ need, and a selection identifier method _ flag for the calculation method, a user-required transmit power spectral density psd, and a user-required demodulation threshold.
3. The method for calculating the downlink support degree of a satellite communication system from spatial spectrum resources according to claim 1, wherein in step S2, according to the spatial spectrum resource identifier ntc _ id, all beams under the resource are obtained, the frequency range of each beam is sequentially read and compared with the target frequency band freq _ need, if there are overlapping frequency bands between the two, the beam information is stored in an array freq _ list, and the beam is used as the beam in the spatial spectrum resource that meets the target frequency band condition.
4. The method for calculating the downlink support of the satellite communication system from the spatial frequency spectrum resources as claimed in claim 1, wherein the step S3 is implemented by dividing the region into K grid points according to the longitude and latitude of the region on the map, each grid point being provided with an earth station of the satellite communication system, called a grid point earth station, and each grid point earth station being identified by a longitude and latitude sequence.
5. The method for calculating the support of spatial spectrum resources to the downlink of a satellite communication system according to claim 1, wherein in step S4, it is tested whether each of the grid point stas meets the minimum communication requirement of the downlink in sequence according to the K grid point stas obtained in step S3, and a flag ok _ flag field is set in each grid point sta to indicate whether the grid point sta meets the requirement, and the flag ok _ flag field is initialized to 0; for each mesh point earth station, performing the following operations on each beam meeting the target frequency band condition in each spatial frequency spectrum resource in freq _ list:
s41, judging whether the mesh point earth station meets the lowest communication requirement of the downlink, detecting the mark ok _ flag field, if the ok _ flag is 1, indicating that the mesh point earth station meets the lowest communication requirement of the downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating the beam; if ok _ flag is 0, continuing to step S42;
s42, calculating the path distance d between the satellite and the grid point earth station in the satellite communication system, wherein the calculation formula is as follows:
Figure FDA0003550554840000021
wherein, theta is the difference between the longitude of the satellite and the grid point earth station, and alpha is the latitude of the grid point earth station; acquiring longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id;
s43, calculating the central frequency f of the wave beam; according to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2;
s44, calculating path loss LF(ii) a Obtaining the path loss L from the earth station of the grid point to the satellite according to the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
s45, calculating satellite transmitting power P; according to a beam group identifier grp _ id corresponding to a beam, acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite; the beam group identification grp _ id is used for identifying a group of beams with certain same performance under the same spatial frequency spectrum resource identification;
s46, calculating satellite gain G; acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as the transmission gain G of the satellite;
s47, calculating the Effective Isotropic Radiated Power (EIRP) of the grid point earth station; according to the path loss LFThe effective omnidirectional radiation power EIRP is obtained by the satellite transmitting power P and the satellite gain G, and the calculation formula is as follows:
EIRP=P+G-LF
s48, detecting whether the mesh point earth station reaches the lowest communication requirement of the satellite communication system downlink under the beam; judging according to the transmitting power spectral density psd required by a user, if the EIRP is more than or equal to the psd, indicating that the earth station of the mesh point meets the lowest communication requirement of a satellite communication system downlink under the wave beam, and setting ok _ flag to be 1; if EIRP < psd, the mesh point earth station can not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
6. The method of calculating support of spatial spectrum resources for a downlink of a satellite communication system according to claim 1,
the step S5 specifically includes:
sequentially calculating whether each grid point earth station meets the lowest communication requirement of a downlink of a satellite communication system according to the K grid point earth stations obtained in the step S3, setting an identification ok _ flag field in each grid point earth station to indicate whether the grid point earth station meets the requirement, and initializing the ok _ flag field to be 0; for each grid point earth station, for each beam meeting the target frequency band condition in each spatial spectrum resource in freq _ list, the following operations are performed:
s51, judging whether the mesh point earth station meets the lowest communication requirement of the satellite communication system downlink; detecting an identification ok _ flag field of the mesh point earth station, if the ok _ flag field is 1, indicating that the mesh point earth station meets the lowest communication requirement of a satellite communication system downlink, directly ending the operation, and jumping to the next mesh point earth station to continue operating a beam of the mesh point earth station; if ok _ flag is 0, proceeding to step S52;
s52, calculating a path d between a satellite and a grid point earth station in the satellite system; the calculation formula is as follows:
Figure FDA0003550554840000041
wherein, theta is the difference between the longitude of the satellite and the grid point earth station, and alpha is the latitude of the grid point earth station; acquiring longitude and latitude coordinates of the satellite by using the space spectrum resource identifier ntc _ id;
s53, calculating the central frequency f of the wave beam; according to the frequency range { freq _ min, freq _ max } of the beam, the beam center frequency is calculated by the formula:
f=(freq_min+freq_max)/2,
s54, calculating path loss LF(ii) a Obtaining path loss L from the mesh point earth station to the satellite communication system from the path d and the beam center frequency fFThe calculation formula is as follows:
LF=92.45+20lg(d×f);
s55, calculating satellite transmitting power P; according to a beam group identifier grp _ id corresponding to a beam, acquiring the maximum transmitting power pwr _ ds _ nbw of the beam group through the grp _ id, and taking the maximum transmitting power pwr _ ds _ nbw as the transmitting power P of the satellite; the beam group identification grp _ id is used for identifying a group of beams with certain same performance under the same spatial frequency spectrum resource identification;
s56, calculating a grid point earth station receiving gain; according to the space spectrum resource identifier ntc _ id and the projection range of the beam on the earth, calculating the contour line coArea of the projection range of the beam, and according to the distance between the central point of the projection range and the grid point earth station, obtaining the receiving gain of the grid point earth station;
s57, calculating satellite receiving gain G; acquiring corresponding gain according to the beam group identifier grp _ id corresponding to the beam, and taking the gain as satellite receiving gain G;
s58, calculating the receiving power P of the grid point earth stationr(ii) a According to the path loss LFThe satellite transmitting power P and the grid point earth station gain are obtained to obtain the grid point earth station receiving power PrThe calculation formula is as follows:
Pr=G+P-LF+gain;
s59, detecting whether the mesh point earth station reaches the lowest communication requirement of the satellite communication system downlink under the beam; judging according to the demodulation threshold required by the user, if P isrIf the wave beam is greater than or equal to threshold, the earth station of the grid point meets the lowest communication requirement of a downlink of a satellite communication system under the wave beam, and ok _ flag is set to be 1; if PrIf < threshold, it means that the earth station of the grid point does not reach the minimum communication requirement of the downlink of the satellite communication system under the beam.
7. The method for calculating support of spatial spectrum resources for a downlink of a satellite communication system according to any one of claims 5 or 6,
the beam group identity grp _ id is used to identify a group of beams with certain same performance under the same spatial spectrum resource identity.
8. The method for calculating the support of spatial spectrum resources in the downlink of the satellite communication system according to claim 1, wherein in step S6, for K mesh point earth stations, the ok _ flag field is detected cyclically, and a num variable is set to record the number of mesh point earth stations that have reached the minimum communication requirement in the downlink of the satellite communication system; if the ok _ flag field of the grid point earth station is 1, adding 1 to the num variable; num is unchanged when the ok _ flag field of the grid point earth station is 0; taking the ratio of the number of the grid point earth stations meeting the minimum communication requirement to the total number of the grid point earth stations as a space frequency spectrumResource support for downlink of satellite communication systemupThe calculation formula is as follows:
Figure FDA0003550554840000051
wherein K is the total number of grid points earth stations.
CN202210262381.5A 2022-03-16 2022-03-16 Method for calculating support degree of space spectrum resources to satellite communication system downlink Active CN114630362B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210262381.5A CN114630362B (en) 2022-03-16 2022-03-16 Method for calculating support degree of space spectrum resources to satellite communication system downlink

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210262381.5A CN114630362B (en) 2022-03-16 2022-03-16 Method for calculating support degree of space spectrum resources to satellite communication system downlink

Publications (2)

Publication Number Publication Date
CN114630362A true CN114630362A (en) 2022-06-14
CN114630362B CN114630362B (en) 2022-10-21

Family

ID=81902981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210262381.5A Active CN114630362B (en) 2022-03-16 2022-03-16 Method for calculating support degree of space spectrum resources to satellite communication system downlink

Country Status (1)

Country Link
CN (1) CN114630362B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6069583A (en) * 1996-05-09 2000-05-30 Agence Spatiale Europeene Receiver for a navigation system, in particular a satellite navigation system
CN109672469A (en) * 2018-11-30 2019-04-23 航天科工空间工程发展有限公司 Avoiding angle calculation method coexists in a kind of GSO and NGSO satellite frequency spectrum
CN112835868A (en) * 2021-01-11 2021-05-25 北京新星宇航科技有限公司 Method for constructing space frequency track resource comprehensive database
CN113282645A (en) * 2021-07-23 2021-08-20 广东粤港澳大湾区硬科技创新研究院 Satellite time sequence parameter analysis method, system, terminal and storage medium
CN113627799A (en) * 2021-08-12 2021-11-09 天津(滨海)人工智能军民融合创新中心 Military benefit-based GSO frequency-rail resource efficiency evaluation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6069583A (en) * 1996-05-09 2000-05-30 Agence Spatiale Europeene Receiver for a navigation system, in particular a satellite navigation system
CN109672469A (en) * 2018-11-30 2019-04-23 航天科工空间工程发展有限公司 Avoiding angle calculation method coexists in a kind of GSO and NGSO satellite frequency spectrum
CN112835868A (en) * 2021-01-11 2021-05-25 北京新星宇航科技有限公司 Method for constructing space frequency track resource comprehensive database
CN113282645A (en) * 2021-07-23 2021-08-20 广东粤港澳大湾区硬科技创新研究院 Satellite time sequence parameter analysis method, system, terminal and storage medium
CN113627799A (en) * 2021-08-12 2021-11-09 天津(滨海)人工智能军民融合创新中心 Military benefit-based GSO frequency-rail resource efficiency evaluation method

Also Published As

Publication number Publication date
CN114630362B (en) 2022-10-21

Similar Documents

Publication Publication Date Title
US7035632B2 (en) Path loss data normalization for growth management of a cellular system
EP3890361B1 (en) Cell longitude and latitude prediction method and device, server, base station, and storage medium
CN101060689B (en) A method and equipment for planning the communication system network
US6282427B1 (en) Selection of location measurement units for determining the position of a mobile communication station
CA2358028C (en) Path loss data normalization for growth management of a cellular system
JP3972755B2 (en) Position measuring method, and terminal device and server used therefor
CN103945331B (en) A kind of localization method for the angle estimation that using WIFI field strength leave away
Ajose et al. Propagation measurements and modelling at 1800 MHz in Lagos Nigeria
CN101277539A (en) Method and system for estimating position of mobile equipment
CN107807346A (en) Adaptive WKNN outdoor positionings method based on OTT Yu MR data
CN108990074B (en) Map information-based power system wireless private network base station construction planning method
CN102480743A (en) Positioning method of cross-cell coverage, apparatus and system thereof
CN114520997B (en) Method, device, equipment and storage medium for positioning 5G network interference source
CN110831017A (en) Site selection method for power system wireless private network base station construction
CN113627799B (en) GSO frequency rail resource efficiency evaluation method based on military benefit
US10454597B1 (en) Systems and methods for locating telecommunication cell sites
CN106888504A (en) Indoor location fingerprint positioning method based on FM Yu DTMB signals
CN114630362B (en) Method for calculating support degree of space spectrum resources to satellite communication system downlink
CN114665950B (en) Method for calculating support degree of space spectrum resource to satellite communication system uplink
CN101848045A (en) Method for measuring reference direction of dual-polarization intelligent antenna business beam
CN111356168B (en) Network coverage simulation method and device
CN113395763B (en) NB terminal high-precision positioning method and system based on 5G/NB-IoT cluster node information fusion
CN115173974A (en) Downlink interference calculation method based on space resource database and contour line interpolation
CN114002717A (en) GSO rail position selection method based on operational risk
CN115378521A (en) Method, device and related equipment for determining wireless channel multipath information

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