CN110278024A - The power system capacity optimization method and device of Communication Satellite Constellation - Google Patents
The power system capacity optimization method and device of Communication Satellite Constellation Download PDFInfo
- Publication number
- CN110278024A CN110278024A CN201910727794.4A CN201910727794A CN110278024A CN 110278024 A CN110278024 A CN 110278024A CN 201910727794 A CN201910727794 A CN 201910727794A CN 110278024 A CN110278024 A CN 110278024A
- Authority
- CN
- China
- Prior art keywords
- transmission power
- mapping relations
- optimal solution
- constraint condition
- interference
- 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
Links
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/18513—Transmission in a satellite or space-based system
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Relay Systems (AREA)
Abstract
The present invention provides the power system capacity optimization methods and device of a kind of Communication Satellite Constellation, it is related to the technical field of satellite communication, Communication Satellite Constellation includes multiple native to this world geosynchronous satellites, the described method includes: obtaining the constraint condition of system capacity performance in Communication Satellite Constellation, constraint condition includes the constraint condition of the constraint condition of wave beam and earth station mapping relations, the constraint condition of transmission power and interference value, and interference value is characterized by mapping relations and transmission power;The optimal solution of transmission power and the optimal solution of mapping relations are calculated according to constraint condition;It is communicated according to the optimal solution of the optimal solution of transmission power and mapping relations, in the case that guarantee is not interfered same frequency geo-synchronous orbit satellite communication system, reaches maximum power system capacity.
Description
Technical field
The present invention relates to technical field of satellite communication, more particularly, to a kind of power system capacity optimization side of Communication Satellite Constellation
Method and device.
Background technique
Classify according to satellite orbit, satellite communication system includes geostationary orbit (Geostationary Earth
Orbit Satellite, GSO) satellite communication system and Non GEO (NGSO) satellite communication system.Compared to GSO satellite,
In being located at, the NGSO satellite system of the low orbit advantages that have that propagation delay is short, link load is small etc..Since single NGSO satellite covers
Lid range is extremely limited, and constellation can be made up of Track desigh and realizes the whole world, full-time covering.
With the increase of NGSO satellite, identical frequency, the satellite of adjacent tracks are more and more, this make NGSO constellation with
The problem of co-channel interference of GSO satellite system is got worse.Further, since needing consideration variable more in practical application and modulation system
It is complex, often there is the situation of resource allocation unevenness, so that NGSO constellation cannot reach optimal power system capacity.
Summary of the invention
The purpose of the present invention is to provide the power system capacity optimization methods and device of Communication Satellite Constellation, guarantee not to same frequency
In the case that geostationary orbit communication system is interfered, reach maximum power system capacity.
In a first aspect, the embodiment of the present invention provides a kind of power system capacity optimization method of Communication Satellite Constellation, the satellite
Communication satellite constellation includes multiple native to this world geosynchronous satellites, which comprises
The constraint condition of system capacity performance in the Communication Satellite Constellation is obtained, the constraint condition includes wave beam and ground
Constraint condition, the constraint condition of the constraint condition of transmission power and interference value of ball station mapping relations, the interference value pass through institute
Mapping relations and the transmission power are stated to characterize;
The optimal solution of the transmission power and the optimal solution of the mapping relations are calculated according to the constraint condition;
It is communicated according to the optimal solution of the optimal solution of the transmission power and the mapping relations.
In alternative embodiments, the constraint condition includes wave beam and earth station need to meet default mapping relations,
Interference value, which need to be less than interference threshold value, transmission power, within the scope of power threshold and each in the Communication Satellite Constellation to be defended
The space-time characterisation variation range of star is corresponding with predetermined trajectory configuration.
In alternative embodiments, the optimal solution of the transmission power and described is calculated according to the constraint condition
The optimal solution of mapping relations, comprising:
Solution calculating is carried out according to interference characteristic and binding feature, the transmission power and mapping for obtaining each chain road are closed
System;
The transmission power of each chain road and mapping relations are iterated calculating, respectively obtain the transmission power
Optimal solution and the mapping relations optimal solution.
In alternative embodiments, solution calculating is carried out according to interference characteristic and binding feature, obtains each chain road
Transmission power and mapping relations, comprising:
Pass through the cumulative and characterization of the product of the transmission power and mapping relations of the multiple native to this world geosynchronous satellite
Interference value calculates the interference value and interference threshold value, obtains the hair for each chain road for meeting communication service requirement
Penetrate power and mapping relations.
In alternative embodiments, the method also includes:
Maximum throughput is calculated according to the optimal solution of the optimal solution of the transmission power and the mapping relations.
In alternative embodiments, it is calculated according to the optimal solution of the optimal solution of the transmission power and the mapping relations
Obtain maximum throughput, further includes:
Maximum throughput is calculated according to the following formula:
Wherein, ai,jIndicate the mapping relations of i-th of wave beam and j-th of earth station, pi,jIndicate i-th of wave beam to j-th
The transmission power size of earth station, ci,jRepresent the maximum communication rate of each of the links, qi,jFor product coefficient, IthFor interference threshold
Value.
In alternative embodiments, the method also includes:
Adjust the modulation system of the Communication Satellite Constellation.
In alternative embodiments, the modulation system includes Adaptive Modulation and Coding mode, pulse-code modulation mode, phase
One of position keying modulation mode, quadrature amplitude modulation mode and minimum shift keying modulation system are a variety of.
In alternative embodiments, in optimal solution and the institute that the transmission power is calculated according to the constraint condition
Before the optimal solution for stating mapping relations, further includes:
Convex optimization problem conversion operation is carried out, so as to meet the canonical form of convex optimization problem after conversion.
Second aspect, the embodiment of the present invention provide a kind of power system capacity optimization device of Communication Satellite Constellation, the satellite
Communication satellite constellation includes multiple native to this world geosynchronous satellites, and described device includes:
Module is obtained, for obtaining the constraint condition of system capacity performance in the Communication Satellite Constellation, the constraint item
Part includes the constraint condition of the constraint condition of wave beam and earth station mapping relations, the constraint condition of transmission power and interference value, institute
Interference value is stated by the mapping relations and the transmission power to characterize;
Computing module, optimal solution and the mapping for the transmission power to be calculated according to the constraint condition are closed
The optimal solution of system;
Communication module, for being communicated according to the optimal solution of the transmission power and the optimal solution of the mapping relations.
The embodiment of the invention provides the power system capacity optimization method and device of a kind of Communication Satellite Constellation, satellite communication stars
Seat includes multiple native to this world geosynchronous satellites, calculates transmitting according to the constraint condition of system capacity performance in Communication Satellite Constellation
The optimal solution of power and optimal solution with the mapping relations of earth station realize the combined optimization of transmission power and mapping relations, can
In the case where guaranteeing not interfere same frequency geostationary orbit communication system, reach maximum power system capacity.
Other features and advantages of the present invention will illustrate in the following description, also, partly become from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention are in specification and attached drawing
Specifically noted structure is achieved and obtained.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below
Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor
It puts, is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of power system capacity optimization method flow chart of Communication Satellite Constellation provided in an embodiment of the present invention;
Fig. 2 is the power system capacity optimization method flow chart of another Communication Satellite Constellation provided in an embodiment of the present invention;
Fig. 3 is the throughput of system curve comparison of native to this world gso satellite communication constellation provided in an embodiment of the present invention
Change one of schematic diagram;
Fig. 4 is the receiving end lump interference curve of native to this world gso satellite communication constellation provided in an embodiment of the present invention
One of changes in contrast schematic diagram;
Fig. 5 is the throughput of system curve comparison of native to this world gso satellite communication constellation provided in an embodiment of the present invention
Change the two of schematic diagram;
Fig. 6 is the receiving end lump interference curve of native to this world gso satellite communication constellation provided in an embodiment of the present invention
The two of changes in contrast schematic diagram;
Fig. 7 is that a kind of power system capacity of Communication Satellite Constellation provided in an embodiment of the present invention optimizes the functional module of device
Figure.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with attached drawing to the present invention
Technical solution be clearly and completely described, it is clear that described embodiments are some of the embodiments of the present invention, rather than
Whole embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creative work premise
Under every other embodiment obtained, shall fall within the protection scope of the present invention.
With the increase of NGSO satellite, identical frequency, the satellite of adjacent tracks are more and more, this make NGSO constellation with
The problem of co-channel interference of GSO satellite system is got worse.According to the regulation of ITU, in, low-orbit satellite must take in communication
Measure avoids generating harmful interference to GSO satellite, and transmission power is controlled under the premise of meeting own communication needs, can be effective
Avoid potential co-channel interference.Research at present in terms of power control is less, and main control strategy has: optimal energy allocation
Technology, descending self-adapting FEC coding techniques, loop-back power control technology, power distributing technique etc..These schemes are substantially all
It is the link-quality such as signal-to-noise ratio, equivalent power flux-density etc. for monitoring NGSO satellite system receiver end, when there may be same frequencies
When interference, transmission power is adjusted by certain rule or method, to avoid potential interference.
Since these schemes acquire a certain degree of difficulty in Project Realization, in practice, relatively conventional way is to close its hair
Machine is penetrated to avoid interference, and other are not generated with the link of interference, transmitter is corresponded to and is generally in maximum power transmission shape
State.Since transmitter is closed, current certain links will be unable to continue as earth station offer communication service.Therefore, logical to guarantee
The case where that believes is continuous, and earth station is generally in the covering of more stars, when a certain link can not continue communication, earth station will be accessed newly
NGSO satellite link with guarantee communication it is continuous.At this point, will be directly influenced according to which kind of rule access NGSO satellite system
The system performances such as call blocking rate, handover overhead.The characteristics of in view of the communication of NGSO satellite system, more stars relatively conventional at present
Covering access strategy has distance priority strategy, cover time preference strategy, maximum elevation preference strategy etc..Further, since the earth
The problems such as fluctuating of topographical surface, the variation of atmospheric environment, atural object block also needs to adjust according to the actual situation in communication
Mode processed, to guarantee communication quality in different environments.
Therefore, in order to guarantee NGSO communication system continuous-stable, and interference is not generated to other GSO communication systems, set
Timing generally requires to consider many variables, meanwhile, a variety of different control strategies such as power control, selecting methods, modulation methods
Formula etc. is complex, it is likely that causes resource allocation on star uneven, the performance of entire NGSO constellation systems is unable to reach most preferably.
Based on this, a kind of the power system capacity optimization method and device of Communication Satellite Constellation provided in an embodiment of the present invention can
In the case where guaranteeing not interfere same frequency geostationary orbit communication system, reach maximum power system capacity.
For convenient for understanding the present embodiment, first to a kind of Communication Satellite Constellation disclosed in the embodiment of the present invention
Power system capacity optimization method describes in detail.
Fig. 1 provides a kind of power system capacity optimization method flow chart of Communication Satellite Constellation for the embodiment of the present invention.
Referring to Fig.1, the power system capacity optimization method of Communication Satellite Constellation is applied to mobile communication Satellite communication system skill
Art field, Communication Satellite Constellation include multiple native to this world geosynchronous satellites, comprising the following steps:
Step S102, obtains the constraint condition of system capacity performance in Communication Satellite Constellation, constraint condition include wave beam with
The constraint condition of the constraint condition of earth station mapping relations, the constraint condition of transmission power and interference value, interference value pass through mapping
Relationship and transmission power characterize.
The optimal solution of transmission power and the optimal solution of mapping relations are calculated according to constraint condition by step S104.
Step S106 is communicated according to the optimal solution of the optimal solution of transmission power and mapping relations.
In the preferred embodiment of practical application, calculated according to the constraint condition of system capacity performance in Communication Satellite Constellation
The optimal solution of transmission power and optimal solution with the mapping relations of earth station realize that transmission power and combining for mapping relations are excellent
Change, it is ensured that in the case where not interfering same frequency geostationary orbit communication system, reach maximum power system capacity.
Wherein, constraint condition includes wave beam and earth station need to meet default mapping relations, interference value need to be less than interference door
Limit value, transmission power need within the scope of power threshold and in Communication Satellite Constellation the space-time characterisation variation range of each satellite with
Predetermined trajectory configuration is corresponding.
It should be noted that constraint is evaded in the interference of NGSO satellite system, advised according to the existing protocol frame of International Telecommunication Association
Fixed, NGSO satellite system must not generate harmful interference to GSO satellite system, and interference evaluation index (interference threshold value) is usually
GSO satellite system receiver end lump interference-to-noise ratio or equivalent power flux-density, and it is directed to different frequency ranges, different interference
Evaluation index sets different interference protection standards, is not limited thereto.
The transmission power of NGSO satellite system constrains, and in order to avoid the interference to other systems in space, while meeting and being
It unites itself communication quality demand, the transmission power P of NGSO system needs to be limited in certain section: where upper limit PmaxFor hair
The maximum value that the machine of penetrating can emit, lower limit PminFor the transmission power minimum value for meeting receiving end carrier-to-noise ratio thresholding, to guarantee itself
Link quality;When transmission power P is equal to zero, then it represents that transmitter-off stops transmitting.
The orbital configuration of NGSO satellite system itself constrains, i.e., NGSO satellite system must be run on the track of its setting,
The space-time characterisation variation range of each NGSO satellite need to be corresponding with the scheduled orbital configuration of NGSO satellite.
In order to clearly know the optimal solution calculation of transmission power and mapping relations, in above-described embodiment
Step S104, comprising:
1, solution calculating is carried out according to interference characteristic and binding feature, the transmission power and mapping for obtaining each chain road are closed
System.
Wherein, pass through the cumulative and characterization of the product of the transmission power and mapping relations of multiple native to this world geosynchronous satellites
Interference value calculates interference value and interference threshold value, obtains the transmitting function for each chain road for meeting communication service requirement
Rate and mapping relations.
2, the transmission power of each chain road and mapping relations are iterated calculating, respectively obtain the optimal of transmission power
The optimal solution of solution and mapping relations.
Further, the above method further include: calculated according to the optimal solution of the optimal solution of transmission power and mapping relations
Power system capacity is optimized to maximum throughput, and according to maximum throughput.
Maximum throughput is calculated according to the following formula:
Wherein, ai,jIndicate the mapping relations of i-th of wave beam and j-th of earth station, pi,jIndicate i-th of wave beam to j-th
The transmission power size of earth station, ci,jRepresent the maximum communication rate of each of the links, qi,jFor product coefficient, IthFor interference threshold
Value.
The embodiment of the present invention needs while meeting above-mentioned constraint condition, on this basis, is closed by transmission power and mapping
The combined optimization of system, distributing system resource is to improve power system capacity.
A NGSO constellation systems and GSO satellite system are considered with the scene coexisted frequently, and NGSO system is by S NGSO satellite
It is formed with N number of earth station, K electricity is installed on every NGSO Satellite and sweeps phased array beam.
Define the transmission power matrix P=[p of NGSO constellation systemsi,j](S·K)×N, wave beam and earth station on NGSO Satellite
Mapping relations matrix A=[ai,j](S·K)×N, wherein pi,jIndicate i-th of wave beam to j-th of earth station transmission power size,
Each pi,jMeet Pmin≤pi,j≤PmaxOr pi,j=0;ai,jIndicate the mapping relations of i-th of wave beam and j-th of earth station, if
I-th of wave beam of NGSO constellation and j NGSO earth station establish link, then correspond to ai,jValue is 1, and otherwise value is 0.Specifically
Implementation method is as shown in Figure 2, comprising the following steps:
Step S201 determines optimization aim, the i.e. a certain system performance of NGSO constellation systems, is not limited thereto, this
Corresponding constraint condition can be arranged according to system capacity performance by taking system capacity performance as an example in inventive embodiments, utilize constraint item
Part constructs the constraint equation for meeting constraint condition, wherein constraint equation need to include interference constraints equation, to power matrix P's
Constraint equation and constraint equation to wave beam mapping relations matrix A.
Here, in order to more accurately calculate the optimal solution of transmission power and mapping relations, according to constraint condition meter
Before calculation obtains the optimal solution of transmission power and the optimal solution of mapping relations, also typically include:
Convex optimization problem conversion operation is carried out, so as to meet the canonical form of convex optimization problem after conversion.
Step S202 suitably converts system performance and constraint equation, so as to meet the standard of convex optimization problem after conversion
Form;
Step S203 solves convex optimization problem, analyzes the dry of all potential links between satellite and earth station
Characteristic and binding feature are disturbed, every latent transmission power p on the link is respectively obtainedi,jWith wave beam and the earth on NGSO Satellite
Mapping relations of standing ai,j。
Step S204 iterates to calculate pi,jAnd ai,jOptimal solutionWithObtain optimal dematrix PoptAnd Aopt。
In solution procedure, every potential link transmission power pi,jThe matrix P at place0It is special by interference characteristic and constraint
Property optimization after obtain transmission power matrix P1, wherein matrix P1Determine mapping relations matrix A1, pass through mapping relations matrix A1Come
Optimize transmission power matrix P1Obtain new emission matrix P2, P2Determination has obtained mapping relations matrix A again2, and so on, no
Disconnected ground mutually iteration to restraining to the endWith
Step S205 adjusts modulation system according to spread state, the link margin in each of the links, is meeting needed for engineering
Under conditions of surplus, link rate is maximized, power system capacity C is madeoptIt is optimal.
As a kind of possible embodiment, consider that power system capacity index is handling capacity T in this examplem.NGSO and GSO constellation
The simulation parameter of satellite system is respectively as shown in the following table 1, table 2:
Table 1
Table 2
Define transmission power matrix P=[pi,j](S·K)×N, wherein S, K, N are provided by upper table, each pi,jMeet 0≤pi,j
≤Pmax。
Define beam link allocation matrix A=[ai,j](S·K)×NIf ai,jI-th of wave beam of ∈ { 0,1 } expression NGSO constellation
With the mapping relations of j-th of NGSO earth station, this example assumes that a wave beam same time can only at most build with One Earth One Family station
Vertical communication link, vice versa, therefore beam allocation element ai,jMeet:
By defined above, we can be by handling capacity TmIt is indicated by following formula:
Wherein, ci,jRepresent the maximum communication rate of each of the links.
Since NGSO constellation systems and GSO satellite system use same communication frequency point, when NGSO constellation i-th of wave beam with
When j-th of earth station establishes downstream communications link, i-th of wave beam of NGSO constellation can generate with frequency GSO earth station receiving end
Interference, corresponding interference value Ii,jIt may be expressed as:
Wherein, Gngt(θ1) it is transmitting gain of the wave beam on GSO earth station direction on NGSO star, Ggr(θ2) it is the GSO earth
It stands in the reception gain where i-th of wave beam in NGSO satellite direction, di,jFor NGSO satellite where i-th of wave beam to the GSO earth
The distance stood, p 'i,jFor NGSO beam transmission power pi,jConvert the equivalent transmission power of overlapped frequency bands, relationship are as follows:
Wherein, WngFor NGSO constellation systems downstream communications bandwidth, WgIt is therefore done for GSO constellation systems downstream communications bandwidth
Disturb value Ii,jWith transmission power pi,jRelationship can be expressed as form:
Wherein, qi,jFor coefficient, then it is as follows to receive the interference I calculating of the lump from NGSO constellation for GSO earth station:
To sum up, NGSO constellation systems should meet following constraint condition:
System emission power pi,jMeet 0≤pi,j≤Pmax。
One wave beam of system can only establish communication link with One Earth One Family station, and vice versa.
System should meet interference and evade condition, and interference value is lower than interference discrimination standard Ith(interference-to-noise ratio threshold value be-
12.2dB)。
Using combined optimization method, to acquire system performance handling capacity TmMaximum value, the specific steps are as follows:
The optimization problem and constraint equation of this problem are constructed first, it may be assumed that
Convex optimization problem, the i.e. corresponding conditions of transformation model are converted by this problem, is allowed to meet the mark of convex optimization problem
Quasi- form.Conversion process is as follows: objective function is to seek handling capacity TmMaximum value, this maximization problems is converted to and seeks minimum
Problem;Optimized variable is transmission power pi,jWith beam link ai,j, wherein pi,jFor continuous variable, ai,jIt, will for discontinuous variable
The discontinuous variable a of 0-1i,jLoosen into the continuous variable that value range is [0,1], final result are as follows:
It solves the problem: listing Lagrangian:
Wherein, λ, μ, η are Lagrange multiplier, Lagrange multiplier are initialized, with Lagrange duality decomposition method
Solve the problem, the dual function of Central Plains problem are as follows:
This dual function is resolved into a primal problem and (SK) × N number of subproblem solves, wherein each subproblem
Optimal solutionWithMeet Karush-Kuhn-Tucker (KKT) optimal condition, iterates to calculate pi,jAnd ai,jOptimal solutionWithIt is iterating to calculateDuring, due to there is the relationship for being mutually coupled, restricting in constraint condition C2 and C3, because
This is introduced into two sub-models in graph theory, and the problem of seeking is converted into bipartite graph least weight match algorithm problem, is calculated using Kuhn-Munkres
Method completes least weight match algorithm, is calculated to update
Iteration calculates the optimal solution of convex optimization problem until convergenceWithObtain optimal dematrix PoptWith
Aopt, and calculate handling capacity TmOptimal value;
On the basis of above-mentioned throughput-optimal value calculates, the modulation system of Communication Satellite Constellation can be adjusted again, makes to handle up
Measure TmReach the more figure of merit.
Wherein, modulation system includes but is not limited to Adaptive Modulation and Coding mode (ACM), pulse-code modulation mode (PCM), phase
Position keying modulation mode (PSK/DPSK/QPSK), quadrature amplitude modulation mode (QAM) and minimum shift keying modulation system
One of (MSK) or it is a variety of.
For in GSO, with the NGSO constellation systems under frequency coexistence scenario, by combined optimization mode, realization has been evaded
The distribution of transmission power, optimization beam link in the satellite of interference, and the system resources such as modulation system are adjusted, with lifting system
Whole volume.
As a kind of optional embodiment, if NGSO constellation systems since interference value exceeds interference threshold, to guarantee NGSO
The normal communication service of constellation systems carries out, then corresponding link is carried out the operation such as on-off, switching, so that mapping relations or distribution
Power emission variation, the combined optimization side by transmission power and mapping relations then provided through the embodiment of the present invention at this time
Method again optimizes system performance.
For the optimization performance for highlighting real-time optimization method of the present invention, now by the embodiment of the present invention and conventional allocation strategy 1,2
It compares and analyzes respectively, wherein allocation strategy 1,2 is described as follows:
Allocation strategy 1: the mode for the airspace isolation that current NGSO constellation systems use is dry with frequency to the nocuousness of GSO to evade
It disturbs, setting link disjoint angle is 6 °, on this basis from the angle of earth station, with wave beam money on " maximum elevation criterion " distribution star
Source.
Allocation strategy 2: interference evasion tactics identical with allocation strategy 1, the airspace that setting link disjoint angle is 6 ° are used
Isolation;Under this constraint, with link communication rate ci,jAs weight, most authority assignment problem is solved using KM algorithm, according to asking
It solves result and completes wave beam resource allocation on corresponding star.
By Fig. 4,6 it is found that three kinds of methods can effectively evade harmful interference of the NGSO constellation systems to GSO satellite system.
As shown in figure 3, corresponding throughput of system is integrally higher than using distribution when NGSO constellation systems are using present invention method
Throughput of system when strategy 1, is averagely higher by 6.1%;As shown in Fig. 5,6, present invention method and allocation strategy 2 are collecting
The difference of corresponding NGSO constellation systems handling capacity is smaller in the case that total interference constraints are weaker;But when rail conditions are poor, GSO
The received lump interference constraints of earth station become strong, and present invention method compares allocation strategy 2, and the worst-case value of handling capacity improves
11.6%.
As the above analysis, the embodiment of the present invention by adjusting the transmission power of transmitter in NGSO constellation systems,
The system resources such as the mapping relations of wave beam and earth station, adjustment modulation system on NGSO Satellite, it is ensured that not to same frequency
In the case that ball geo-stationary orbit communication system is interfered, optimized allocation of resources is optimal power system capacity.
As shown in fig. 7, the embodiment of the present invention also provides a kind of power system capacity optimization device of Communication Satellite Constellation, satellite is logical
Believe that constellation includes multiple native to this world geosynchronous satellites, device includes:
Module is obtained, for obtaining the constraint condition of system capacity performance in Communication Satellite Constellation, constraint condition includes wave
Constraint condition, the constraint condition of the constraint condition of transmission power and interference value of beam and earth station mapping relations, interference value pass through
Mapping relations and transmission power characterize;
Computing module, for the optimal solution of transmission power and the optimal solution of mapping relations to be calculated according to constraint condition;
Communication module, for being communicated according to the optimal solution of transmission power and the optimal solution of mapping relations.
The power system capacity of Communication Satellite Constellation provided in an embodiment of the present invention optimizes device, defends with provided by the above embodiment
The power system capacity optimization method technical characteristic having the same of star communication satellite constellation reaches so also can solve identical technical problem
To identical technical effect.
The power system capacity optimization method of Communication Satellite Constellation provided by the embodiment of the present invention and the computer program of device
Product, the computer readable storage medium including storing program code, the instruction that said program code includes can be used for executing
Previous methods method as described in the examples, specific implementation can be found in embodiment of the method, and details are not described herein.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description
It with the specific work process of device, can refer to corresponding processes in the foregoing method embodiment, details are not described herein.
In addition, in the description of the embodiment of the present invention unless specifically defined or limited otherwise, term " installation ", " phase
Even ", " connection " shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or be integrally connected;It can
To be mechanical connection, it is also possible to be electrically connected;It can be directly connected, can also can be indirectly connected through an intermediary
Connection inside two elements.For the ordinary skill in the art, above-mentioned term can be understood at this with concrete condition
Concrete meaning in invention.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
In the description of the present invention, it should be noted that term " center ", "upper", "lower", "left", "right", "vertical",
The orientation or positional relationship of the instructions such as "horizontal", "inner", "outside" be based on the orientation or positional relationship shown in the drawings, merely to
Convenient for description the present invention and simplify description, rather than the device or element of indication or suggestion meaning must have a particular orientation,
It is constructed and operated in a specific orientation, therefore is not considered as limiting the invention.In addition, term " first ", " second ",
" third " is used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance.
The embodiment of the present invention also provides a kind of electronic equipment, including memory, processor and storage are on a memory and can
The computer program run on a processor, processor realize satellite communication provided by the above embodiment when executing computer program
The step of power system capacity optimization method of constellation.
The embodiment of the present invention also provides a kind of computer readable storage medium, and meter is stored on computer readable storage medium
Calculation machine program executes the power system capacity optimization side of the Communication Satellite Constellation of above-described embodiment when computer program is run by processor
The step of method.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features;And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention
Within the scope of.
Claims (10)
1. a kind of power system capacity optimization method of Communication Satellite Constellation, which is characterized in that the Communication Satellite Constellation includes multiple
Native to this world geosynchronous satellite, which comprises
The constraint condition of system capacity performance in the Communication Satellite Constellation is obtained, the constraint condition includes wave beam and earth station
The constraint condition of the constraint condition of mapping relations, the constraint condition of transmission power and interference value, the interference value are reflected by described
Relationship and the transmission power are penetrated to characterize;
The optimal solution of the transmission power and the optimal solution of the mapping relations are calculated according to the constraint condition;
It is communicated according to the optimal solution of the optimal solution of the transmission power and the mapping relations.
2. system according to claim 1 capacity optimization method, which is characterized in that the constraint condition includes wave beam and ground
Default mapping relations that ball station need to meet, interference value need to be less than interference threshold value, transmission power need within the scope of power threshold with
And in the Communication Satellite Constellation each satellite space-time characterisation variation range it is corresponding with predetermined trajectory configuration.
3. power system capacity optimization method according to claim 2, which is characterized in that be calculated according to the constraint condition
The optimal solution of the optimal solution of the transmission power and the mapping relations, comprising:
Solution calculating is carried out according to interference characteristic and binding feature, obtains the transmission power and mapping relations of each chain road;
The transmission power of each chain road and mapping relations are iterated calculating, respectively obtain the transmission power most
The optimal solution of excellent solution and the mapping relations.
4. power system capacity optimization method according to claim 3, which is characterized in that according to interference characteristic and binding feature into
Row, which solves, to be calculated, and the transmission power and mapping relations of each chain road are obtained, comprising:
Pass through the cumulative and characterization interference of the product of the transmission power and mapping relations of the multiple native to this world geosynchronous satellite
Value, the interference value and interference threshold value are calculated, the transmitting function for each chain road for meeting communication service requirement is obtained
Rate and mapping relations.
5. system according to claim 1 capacity optimization method, which is characterized in that the method also includes:
Maximum throughput is calculated according to the optimal solution of the optimal solution of the transmission power and the mapping relations.
6. power system capacity optimization method according to claim 3, which is characterized in that according to the optimal solution of the transmission power
Maximum throughput is calculated with the optimal solution of the mapping relations, further includes:
Maximum throughput is calculated according to the following formula:
Wherein, ai,jIndicate the mapping relations of i-th of wave beam and j-th of earth station, pi,jIndicate i-th of wave beam to j-th of earth
The transmission power size stood, ci,jRepresent the maximum communication rate of each of the links, qi,jFor product coefficient, IthFor interference threshold value.
7. system according to claim 1 capacity optimization method, which is characterized in that the method also includes:
Adjust the modulation system of the Communication Satellite Constellation.
8. power system capacity optimization method according to claim 7, which is characterized in that the modulation system includes adaptive adjusts
Coding mode, pulse-code modulation mode, phase keying modulation system, quadrature amplitude modulation mode and minimum shift keying modulation methods processed
One of formula is a variety of.
9. system according to claim 1 capacity optimization method, which is characterized in that calculated according to the constraint condition
To the transmission power optimal solution and the mapping relations optimal solution before, further includes:
Convex optimization problem conversion operation is carried out, so as to meet the canonical form of convex optimization problem after conversion.
10. a kind of power system capacity of Communication Satellite Constellation optimizes device, which is characterized in that the Communication Satellite Constellation includes multiple
Native to this world geosynchronous satellite, described device include:
Module is obtained, for obtaining the constraint condition of system capacity performance in the Communication Satellite Constellation, the constraint condition packet
The constraint condition of the constraint condition of wave beam and earth station mapping relations, the constraint condition of transmission power and interference value is included, it is described dry
Value is disturbed by the mapping relations and the transmission power to characterize;
Computing module, for the optimal solution and the mapping relations of the transmission power to be calculated according to the constraint condition
Optimal solution;
Communication module, for being communicated according to the optimal solution of the transmission power and the optimal solution of the mapping relations.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910727794.4A CN110278024B (en) | 2019-08-07 | 2019-08-07 | System capacity optimization method and device for satellite communication constellation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910727794.4A CN110278024B (en) | 2019-08-07 | 2019-08-07 | System capacity optimization method and device for satellite communication constellation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110278024A true CN110278024A (en) | 2019-09-24 |
CN110278024B CN110278024B (en) | 2020-07-24 |
Family
ID=67965868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910727794.4A Active CN110278024B (en) | 2019-08-07 | 2019-08-07 | System capacity optimization method and device for satellite communication constellation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110278024B (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111865398A (en) * | 2020-07-01 | 2020-10-30 | 哈尔滨工业大学(深圳) | Satellite-ground transmission method under large-scale LEO satellite deployment |
CN112994778A (en) * | 2021-02-07 | 2021-06-18 | 哈尔滨工业大学 | High-throughput satellite beam on-demand scheduling method based on service priority and rate self-adaption |
CN113536550A (en) * | 2021-06-29 | 2021-10-22 | 中国人民解放军63920部队 | Situation map generation method and device, nonvolatile storage medium and processor |
CN113595616A (en) * | 2021-07-27 | 2021-11-02 | 北京邮电大学 | NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization |
CN113691332A (en) * | 2021-09-10 | 2021-11-23 | 北京科技大学 | Co-channel interference characterization method and device of low-earth-orbit satellite communication system |
CN113825238A (en) * | 2021-09-03 | 2021-12-21 | 天地信息网络研究院(安徽)有限公司 | Flexible signaling beam multi-scene application method in satellite communication |
CN114301512A (en) * | 2021-12-10 | 2022-04-08 | 西安电子科技大学 | Multilayer low-orbit satellite network-oriented bearing capacity assessment method and system |
CN114553300A (en) * | 2022-02-18 | 2022-05-27 | 清华大学 | Method and system for rapidly allocating wireless resources of multi-antenna earth station of satellite system |
CN116865843A (en) * | 2023-09-05 | 2023-10-10 | 武汉能钠智能装备技术股份有限公司 | System and method for distributing links between satellites of satellite system |
CN117375706A (en) * | 2023-12-04 | 2024-01-09 | 成都本原星通科技有限公司 | Low-orbit inter-satellite interference optimization method and system for receiving end |
CN117544255A (en) * | 2023-11-09 | 2024-02-09 | 中国人民解放军军事科学院系统工程研究院 | NGSO satellite interference condition matrix manifold visualization method and device |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998020634A2 (en) * | 1996-11-07 | 1998-05-14 | Teledesic Llc | Frequency sharing for satellite communication system |
US6259913B1 (en) * | 1998-08-06 | 2001-07-10 | Alcatel | Method of allocating links between a set of areas and a set of satellites |
US20040157646A1 (en) * | 1995-02-22 | 2004-08-12 | Raleigh Gregory Gene | Method and apparatus for adaptive transmission beam forming in a wireless communication system |
CN102355294A (en) * | 2011-11-01 | 2012-02-15 | 东南大学 | Multipoint coordinated beam forming and power allocation method for single base station power constraint |
CN105744531A (en) * | 2016-02-04 | 2016-07-06 | 中国空间技术研究院 | Inline interference suppression based geostationary orbit (GEO) and non-geostationary orbit (NGEO) communication satellite spectrum sharing method |
CN106027138A (en) * | 2016-05-05 | 2016-10-12 | 清华大学 | Ground station system and method for avoiding collinear interference with geostationary satellite |
CN107276660A (en) * | 2017-06-22 | 2017-10-20 | 清华大学 | Resource allocation methods and device in non-orthogonal multiple air-ground coordination communication system |
CN107342804A (en) * | 2016-04-28 | 2017-11-10 | 北京大学 | Suppress the method for adjacent star interference in satellite communication system |
CN108419286A (en) * | 2018-01-18 | 2018-08-17 | 北京邮电大学 | A kind of allocation algorithm in face of 5G UAV Communications united beam and power |
CN108631840A (en) * | 2018-07-05 | 2018-10-09 | 东南大学 | Multi-beam satellite system robust pre-coding method based on outage probability constraint |
US10104641B2 (en) * | 2013-03-01 | 2018-10-16 | Intel IP Corporation | Network selection in a heterogeneous network |
CN108882245A (en) * | 2018-07-03 | 2018-11-23 | 中国人民解放军陆军工程大学 | GEO and L EO cognitive satellite network and dynamic frequency allocation method thereof |
CN109412682A (en) * | 2018-12-19 | 2019-03-01 | 北京卫星信息工程研究所 | Low rail constellation satellite beams frequency resource allocation method in one kind |
CN109743098A (en) * | 2018-12-30 | 2019-05-10 | 清华大学 | Frequency spectrum sharing method, device and the electronic equipment of spectra sharing system |
-
2019
- 2019-08-07 CN CN201910727794.4A patent/CN110278024B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040157646A1 (en) * | 1995-02-22 | 2004-08-12 | Raleigh Gregory Gene | Method and apparatus for adaptive transmission beam forming in a wireless communication system |
WO1998020634A2 (en) * | 1996-11-07 | 1998-05-14 | Teledesic Llc | Frequency sharing for satellite communication system |
US6259913B1 (en) * | 1998-08-06 | 2001-07-10 | Alcatel | Method of allocating links between a set of areas and a set of satellites |
CN102355294A (en) * | 2011-11-01 | 2012-02-15 | 东南大学 | Multipoint coordinated beam forming and power allocation method for single base station power constraint |
US10104641B2 (en) * | 2013-03-01 | 2018-10-16 | Intel IP Corporation | Network selection in a heterogeneous network |
CN105744531A (en) * | 2016-02-04 | 2016-07-06 | 中国空间技术研究院 | Inline interference suppression based geostationary orbit (GEO) and non-geostationary orbit (NGEO) communication satellite spectrum sharing method |
CN107342804A (en) * | 2016-04-28 | 2017-11-10 | 北京大学 | Suppress the method for adjacent star interference in satellite communication system |
CN106027138A (en) * | 2016-05-05 | 2016-10-12 | 清华大学 | Ground station system and method for avoiding collinear interference with geostationary satellite |
CN107276660A (en) * | 2017-06-22 | 2017-10-20 | 清华大学 | Resource allocation methods and device in non-orthogonal multiple air-ground coordination communication system |
CN108419286A (en) * | 2018-01-18 | 2018-08-17 | 北京邮电大学 | A kind of allocation algorithm in face of 5G UAV Communications united beam and power |
CN108882245A (en) * | 2018-07-03 | 2018-11-23 | 中国人民解放军陆军工程大学 | GEO and L EO cognitive satellite network and dynamic frequency allocation method thereof |
CN108631840A (en) * | 2018-07-05 | 2018-10-09 | 东南大学 | Multi-beam satellite system robust pre-coding method based on outage probability constraint |
CN109412682A (en) * | 2018-12-19 | 2019-03-01 | 北京卫星信息工程研究所 | Low rail constellation satellite beams frequency resource allocation method in one kind |
CN109743098A (en) * | 2018-12-30 | 2019-05-10 | 清华大学 | Frequency spectrum sharing method, device and the electronic equipment of spectra sharing system |
Non-Patent Citations (2)
Title |
---|
WEI ZHANG等: "Capacity Model and Constraints Analysis for Integrated Remote Wireless Sensor and Satellite Network in Emergency Scenarios", 《SENSORS》 * |
靳瑾等: "全球动态场景下非静止轨道通信星座干扰发生概率和系统可用性", 《清华大学学报(自然科学版)》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111865398A (en) * | 2020-07-01 | 2020-10-30 | 哈尔滨工业大学(深圳) | Satellite-ground transmission method under large-scale LEO satellite deployment |
CN112994778A (en) * | 2021-02-07 | 2021-06-18 | 哈尔滨工业大学 | High-throughput satellite beam on-demand scheduling method based on service priority and rate self-adaption |
CN113536550B (en) * | 2021-06-29 | 2024-01-12 | 中国人民解放军63920部队 | Situation map generation method and device, nonvolatile storage medium and processor |
CN113536550A (en) * | 2021-06-29 | 2021-10-22 | 中国人民解放军63920部队 | Situation map generation method and device, nonvolatile storage medium and processor |
CN113595616A (en) * | 2021-07-27 | 2021-11-02 | 北京邮电大学 | NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization |
CN113595616B (en) * | 2021-07-27 | 2022-05-17 | 北京邮电大学 | NGSO intersystem interference mitigation method based on satellite antenna beam pointing optimization |
CN113825238A (en) * | 2021-09-03 | 2021-12-21 | 天地信息网络研究院(安徽)有限公司 | Flexible signaling beam multi-scene application method in satellite communication |
CN113825238B (en) * | 2021-09-03 | 2024-02-20 | 天地信息网络研究院(安徽)有限公司 | Flexible signaling beam multi-scene application method in satellite communication |
CN113691332A (en) * | 2021-09-10 | 2021-11-23 | 北京科技大学 | Co-channel interference characterization method and device of low-earth-orbit satellite communication system |
CN114301512A (en) * | 2021-12-10 | 2022-04-08 | 西安电子科技大学 | Multilayer low-orbit satellite network-oriented bearing capacity assessment method and system |
CN114301512B (en) * | 2021-12-10 | 2022-10-14 | 西安电子科技大学 | Method and system for evaluating bearing capacity of multilayer low-orbit satellite network |
CN114553300A (en) * | 2022-02-18 | 2022-05-27 | 清华大学 | Method and system for rapidly allocating wireless resources of multi-antenna earth station of satellite system |
CN114553300B (en) * | 2022-02-18 | 2023-02-17 | 清华大学 | Method and system for rapidly allocating wireless resources of multi-antenna earth station of satellite system |
CN116865843B (en) * | 2023-09-05 | 2023-12-26 | 武汉能钠智能装备技术股份有限公司 | System and method for distributing links between satellites of satellite system |
CN116865843A (en) * | 2023-09-05 | 2023-10-10 | 武汉能钠智能装备技术股份有限公司 | System and method for distributing links between satellites of satellite system |
CN117544255A (en) * | 2023-11-09 | 2024-02-09 | 中国人民解放军军事科学院系统工程研究院 | NGSO satellite interference condition matrix manifold visualization method and device |
CN117544255B (en) * | 2023-11-09 | 2024-07-02 | 中国人民解放军军事科学院系统工程研究院 | NGSO satellite interference condition matrix manifold visualization method |
CN117375706A (en) * | 2023-12-04 | 2024-01-09 | 成都本原星通科技有限公司 | Low-orbit inter-satellite interference optimization method and system for receiving end |
CN117375706B (en) * | 2023-12-04 | 2024-03-12 | 成都本原星通科技有限公司 | Low-orbit inter-satellite interference optimization method and system for receiving end |
Also Published As
Publication number | Publication date |
---|---|
CN110278024B (en) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110278024A (en) | The power system capacity optimization method and device of Communication Satellite Constellation | |
US8902808B2 (en) | Interference mitigation with scheduling and dynamic power spectrum allocation for wireless networks | |
JP2012510749A (en) | Method and apparatus for power allocation in a multi-carrier system | |
US9338753B2 (en) | Method and apparatus for performance management in wireless backhaul networks via power control | |
Blough et al. | Optimal one-shot scheduling for MIMO networks | |
Wang et al. | Admission control and power allocation for NOMA-based satellite multi-beam network | |
Wang et al. | Spectrum optimization for cognitive satellite communications with cournot game model | |
Dimitrov et al. | Radio resource management techniques for high throughput satellite communication systems | |
Lee et al. | Interference mitigation for reverse spectrum sharing in B5G/6G satellite-terrestrial networks | |
Wang et al. | Adaptive beam pattern selection and resource allocation for NOMA-based LEO satellite systems | |
Barceló-Lladó et al. | Distributed power and carrier allocation in multibeam satellite uplink with individual SINR constraints | |
Khorov et al. | Joint power control and time division to improve spectral efficiency in dense Wi-Fi networks | |
Garnaev et al. | Bargaining over fair performing dual radar and communication task | |
Kai et al. | Optimal scheduling and power control for in-band full-duplex communication in WLANs | |
Long et al. | Throughput optimization in energy harvesting based cognitive iot with cooperative sensing | |
Zhang et al. | Power control for NOMA HAP-UAV integrated aerial access networks | |
Shaat et al. | Joint flow control and link scheduling in hybrid terrestrial-satellite wireless backhauling network | |
CN105409268B (en) | network capacity and coverage optimization method and device | |
Li et al. | Fixed-point algorithms for energy-efficient power allocation in spectrum-sharing wireless networks | |
Nakahira et al. | Capacity and quality enhancement using an adaptive resource allocation for multi-beam mobile satellite communication systems | |
Khorov et al. | Cloud-based Management of Energy-Efficient Dense IEEE 802.11 ax Networks | |
Guo et al. | Distributed power and rate allocation with fairness for cognitive radios in wireless ad hoc networks | |
Wang | Multi-objective Resource Allocation for UAV-assisted Air-Ground Integrated MC-NOMA Networks | |
Chai et al. | DRL-Based Dynamic Resource Allocation for Multi-Beam Satellite Systems | |
CN114710195B (en) | Low-orbit satellite energy-efficient resource allocation method based on beam hopping technology |
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 |