CN104320174B - A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge - Google Patents
A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge Download PDFInfo
- Publication number
- CN104320174B CN104320174B CN201410471326.2A CN201410471326A CN104320174B CN 104320174 B CN104320174 B CN 104320174B CN 201410471326 A CN201410471326 A CN 201410471326A CN 104320174 B CN104320174 B CN 104320174B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- user
- msup
- 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.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0426—Power distribution
- H04B7/043—Power distribution using best eigenmode, e.g. beam forming or beam steering
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
- H04B7/046—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account
- H04B7/0465—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account taking power constraints at power amplifier or emission constraints, e.g. constant modulus, into account
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0619—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal using feedback from receiving side
- H04B7/0621—Feedback content
- H04B7/063—Parameters other than those covered in groups H04B7/0623 - H04B7/0634, e.g. channel matrix rank or transmit mode selection
Abstract
The present invention relates to a kind of satellite multi-beam cooperation transmission method.Effective utilization of multi-beam satellite mobile communication system medium wave interfascicular AF panel and frequency spectrum resource is solved, and is only used for the confinement problems of fixed satellite system.Method and step includes：A. user in the cooperative beam area of coverage is grouped at random, the user in each collaboration user group distributes a frequency range respectively from the different beams area of coverage, every group of user；B. satellite is calculated to the link effective transmission gain of user according to known satellite subscriber channel partial information；C. optimal group pre-coding matrix is calculated；D. precoding is carried out to the transmission signal of cooperative beam according to the optimal pre-coding matrix calculated.It is an advantage of the invention that efficiently solving the limitation of prior art, it can apply in satellite mobile communication system.Power gain is obtained, compared to traditional frequency multiplexing method, systematic function is greatly improved.
Description
Technical field
The present invention relates to a kind of technical field of satellite communication, more particularly, to it is a kind of be related to multibeam satellite system based on
The satellite multi-beam cooperation transmission method of partial channel knowledge.
Background technology
Multi-beam antenna technology provide not only flexible ground communication covering, and can realize higher data transfer rate,
Obtained a wide range of applications in modern satellite communications system.However, in multi-beam satellite system, due to same between wave beam
Channel disturbance and rare frequency spectrum resource limitation, hinder the further lifting of multibeam satellite system performance.
Suppress for inter-beam interference and resource allocation problem, the multibeam satellite system solution applied at present is,
Utilize traditional frequency multiplexing method；But, the system spectral efficiency based on this method is low.In order to solve this problem,
Newest investigative technique mainly employs the method for suppressing inter-cell interference in land honeycomb mobile communication, including cooperation transmission skill
Art and interference coordination technique.But, a basic assumption condition of these methods is that satellite transmitter or control juncture station need
There is accurate terrestrial user channel information；For mobile satellite communication system, the mobile communication system of synchronous satellite is based particularly on
System, due to propagation delay, this condition is very inappeasable.Therefore, these technologies and method proposed at present are only used for solid
Determine satellite system.
The content of the invention
The present invention mainly solves having for multi-beam satellite mobile communication system medium wave interfascicular AF panel and frequency spectrum resource
Utilizing question is imitated, and the confinement problems for being only used for fixed satellite system improve transmission speed there is provided one kind suppression interference
Rate, it is not limited only to the satellite multi-beam cooperation transmission method based on partial channel knowledge of fixed satellite system.
The above-mentioned technical problem of the present invention is mainly what is be addressed by following technical proposals：One kind is based on local channel
The satellite multi-beam cooperation transmission method of information, comprises the following steps：
A. user in the cooperative beam area of coverage is grouped at random, the user in each collaboration user group is not respectively from
Same footprint of a beam, every group of user distributes a frequency range；
B. satellite is calculated to the link effective transmission gain of user according to known satellite subscriber channel partial information：Using Lay
This shade channel model exports the outage probability expression formula of satellite user link, according to User Part channel information and user link
Outage probability requirement, tries to achieve the effective transmission gain of subscriber channel；
C. optimal group pre-coding matrix is calculated：According to partial channel knowledge, system power constraints and user rate
Request obtains the optimization problem construction of optimal precoding and power distribution object function, and Solve problems construction obtains optimal precoding
Matrix；
D. precoding is carried out to the transmission signal of cooperative beam according to the optimal pre-coding matrix calculated.The present invention is based on
The multi-beam cooperation transmission method of partial channel knowledge, it is excellent according to the impaction of partial channel state information and quality of service requirement of user
Change the precoding vectors of control signal and the distribution of resource, efficiently against the interference and the influence of channel fading between wave beam.By
In the slow time-varying characteristic of partial channel knowledge, proposition method efficiently solves the limitation of prior art, is multi-beam satellite
The performance improvement of GSM provides a kind of solution.
As a kind of preferred scheme, step a detailed processes are：Multi-beam cooperation transmission method realized by juncture station,
Critical point station administration in K wave beam, K footprint of a beam and is uniform-distribution with N number of user, using randomly selected method by user
Packet, every group of K user is respectively from the K different beams area of coverage.
As a kind of preferred scheme, step b detailed processes are：
The signal of every group of user is cooperated using multi-beam and sent, then terrestrial user k reception signal can be expressed as：
Wherein k is k-th of user, akIt is user k fading channel coefficients；wjFor user k transmitting precoding weight vector；
SjIt is the predetermined data-signal for being sent to user；It is the direction vector of user, K represents juncture station pipe
K wave beam of reason, or it is expressed as G=[g with radiation gain matrix1 g2…gk], it reflects satellite antenna to terrestrial user
Radiation gain characteristic；nkIt is to receive the Gaussian noise in signal；For the beam gain of k-th of beam center, it is contained
Transmitter antenna gain (dBi), path loss, the influence of receiving antenna gain and noise power, are defined as
L in formula (2)kFor the distance between user k and satellite, λ is carrier frequency, kBFor Boltzmann constant,It is user
Receiving antenna gain,It is maximum satellite transmitting antenna gain, depending on the mode parameter α, T of transmitting antennaRAnd BWPoint
Not Wei receiver noise temperature and transmission link bandwidth；
Then cooperative transmission system downlink user k reception signal Signal to Interference plus Noise Ratio can be expressed as
Therefore, user k channel capacity is Ck=log2(1+Γk)；
Satellite is calculated to user link outage probability：
Wherein COut, kIt is the target capacity for meeting outage probability, as outage capacity；
Using this shadow fading channel model of Lay, the probability density function of fading channel is：
Wherein channel parameterΩ is the flat of sight component
Equal power, m is Nakagami distributed constants, 2b0It is the mean power of scattering component,1F1(m, 1；c0γk) it is interflow hypergeometry letter
Number, thresholdingFor：
Solution formula (4) is obtained
Pout=Φ (uk)=1-exp (- uk)Q(uk) (7)
In formula (7)Q(uk) it is a convergent unlimited item value of series,
Using interior value of series in limited item number and approximant (8), obtain
Satellite user partial channel knowledge includes the directional information { g of subscriber channelj, j=1,2 ..., K } and statistic letter
Breath(the wherein subscriber channel under the User Part channel information known case
Directional information { gj, j=1,2 ..., K } it can be tried to achieve using the channel estimation of up-link, statistic information { Ω ', b '0, m } and can be
User terminal is estimated, and feeds back to satellite launch pole or juncture station), the user terminal probability P based on determinationoutIt is required that, using iteration
Dichotomy calculating formula (9), obtain the effective transmission gain of subscriber channel
According to ukDefinition, user's outage capacity is represented by
Assuming that the band bandwidth B of every group of user's distributionW, then user k signal transmission rate can be expressed as Rk=
BWCOut, k。
As a kind of preferred scheme, step c detailed processes are：
Based on the rate request of User Part channel information, system power constraint and user, construction precoding and power point
The optimization problem matched somebody with somebody, the target of the optimization problem is：Under fairness condition between ensureing user, maximize the speed of user,
Problem is：
P in formula (11)maxConstrained for system peak power, FkIt is user k rate request；
Using ZF pre-coding power optimized algorithm, orderAnd by ZF pre-coding matrix
Applied in (11), formula (11) is rewritten as
Solution formula (12), obtains optimal power allocation { pk, k=1,2 ..., K }, wherein user k distribution power is
[(GGT)-1]K, kpk；
Formula (11) is the feasibility Solve problems of a convex optimization problem, has many ripe fast algorithms to use.
Using precoding optimized algorithm, formula (11) is rewritten as
C in formula (13)2It is a convex function, C3It is a convex Second-order cone programming constraint, is solved using the dichotomy of iteration
Obtain optimal precoding vectors { wk, k=1,2 ..., K }.
As a kind of preferred scheme, step d detailed process is：
According to the optimal precoding vectors obtainedThe transmission signal of cooperative system is：
Therefore, it is an advantage of the invention that：
1. due to the slowly varying behavior that partial channel knowledge relative user position changes, the method for proposition is efficiently solved now
There is the limitation of technology, can apply in satellite mobile communication system.
2. although the perfectly correlated property of satellite multiple antennas, system can not obtain fading diversity gain；But, due to using association
Make precoding transmissions, the signal of each user is sent by multiple wave beams, and system can be obtained a kind of equivalent ' array gain ', or
Person, which says, obtains a kind of power gain, and compared to traditional frequency multiplexing method, systematic function is greatly improved.
Brief description of the drawings
Accompanying drawing 1 is a kind of structural representation of multibeam satellite system of the present invention；
Accompanying drawing 2 is a kind of method flow schematic diagram of the invention；
Accompanying drawing 3 is the first position distribution schematic diagram of collaboration user in the present invention；
Accompanying drawing 4 is second of position distribution schematic diagram of collaboration user in the present invention；
Accompanying drawing 5 is the third position distribution schematic diagram of collaboration user in the present invention；
Accompanying drawing 6 is cooperation transmission method and legacy frequencies multiplexing method performance comparision in collaboration user the first position distribution
Schematic diagram；
Accompanying drawing 7 is cooperation transmission method and legacy frequencies multiplexing method performance comparision in second of position distribution of collaboration user
Schematic diagram；
Accompanying drawing 8 is cooperation transmission method and legacy frequencies multiplexing method performance comparision in collaboration user the third position distribution
Schematic diagram.
Embodiment
Below by embodiment, and with reference to accompanying drawing, technical scheme is described in further detail.
Embodiment：
A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge of the present embodiment, runs in system, many
Beam satellite system is as shown in figure 1, include juncture station, multi-beam satellite and ground based mobile subscriber, and wherein black round dot is represented
Mobile subscriber.Setting juncture station has obtained the partial channel knowledge of user, and these channel informations are slow time-varyings, local channel letter
Breath includes the directional information { g of userj, j=1,2 ..., K } and statistic information
Wherein direction main information relies on the position of user and the radiation mode of antenna, the influence very little of frequency, therefore can utilize upper
The channel estimation of line link is obtained；The mean power Ω ' of sight component, the mean power b ' of scattering component0With Nakagami points
Cloth parameter m can estimate in user terminal, and feed back to juncture station；The flow of this method is as shown in Figure 2.
The present embodiment method includes following steps：
A. user in the cooperative beam area of coverage is grouped at random, the user in each collaboration user group is not respectively from
Same footprint of a beam, every group of user distributes a frequency range；Assuming that have N number of user in all footprint of a beams of critical point station administration,
User is divided into by M groups using random packet, every group includes K user, each user is from different footprint of a beams, so
Afterwards, every group of user distributes a frequency range, a width of B of band of each frequency rangeW.Each satellite critical point station administration K is assumed in the present embodiment
There are a user, therefore N=7 in=7 wave beams, each footprint of a beam.
B. the effective transmission gain of link is calculated according to known satellite subscriber channel partial information.Effective transmission gain coefficient
u*Satellite is described to a kind of Statistical Parameters of user link, it relies on fading characteristic and the interruption of link of association's channel
Probability level, in the present embodiment, it is contemplated that a kind of simple scene：The user of cooperation is in identical shade environment.
According to this shade channel model of Lay, export effective transmission gain is with systematic parameter and the relation of channel parameter, process：
The signal of every group of user is cooperated using multi-beam and sent, then terrestrial user k reception signal can be expressed as：
Wherein k is k-th of user, akIt is user k fading channel coefficients；wjFor user k transmitting precoding weight vector；
SjIt is the predetermined data-signal for being sent to user；It is the direction vector of user, K represents juncture station pipe
K wave beam of reason, or it is expressed as G=[g with radiation gain matrix1 g2…gK], it reflects satellite antenna to terrestrial user
Radiation gain characteristic；nkIt is to receive the Gaussian noise in signal；For the beam gain of k-th of beam center, it is contained
Transmitter antenna gain (dBi), path loss, the influence of receiving antenna gain and noise power, are defined as
L in formula (2)kFor the distance between user k and satellite, λ is carrier frequency, kBFor Boltzmann constant,It is user
Receiving antenna gain,It is maximum satellite transmitting antenna gain, depending on the mode parameter α, T of transmitting antennaRAnd BWRespectively
For receiver noise temperature and the bandwidth of transmission link；
Then cooperative transmission system downlink user k reception signal Signal to Interference plus Noise Ratio can be expressed as
Therefore, user k channel capacity is Ck=log2(1+Γk)；
Satellite is calculated to user link outage probability：
Wherein COut, kIt is the target capacity for meeting outage probability, as outage capacity；
Using this shadow fading channel model of Lay, the probability density function of fading channel is：
Wherein channel parameterΩ is sight component
Mean power, m be Nakagami distributed constants, 2b0It is the mean power of scattering component,1F1(m, 1；c0γk) super several to collaborate
What function, thresholdingFor：
Solution formula (4) is obtained
Pout=Φ (uk)=1-exp (- uk)Q(uk) (7)
In formula (7)Q(uk) it is a convergent unlimited item value of series,
Using interior value of series in limited item number and approximant (8), effective transmission gain and systematic parameter and letter are obtained
The relation of road parameter
In this example, it is assumed that outage probability requirement is Pout=0.01, it is average that shade channel circumstance employs Loo ' s
Shade ambient measurements：Mean power Ω=0.835, m=10.1 of sight component, the mean power 2b of scattering component0=
2.52.According to (9) formula, and dichotomy is used, the equivalent gain u of link can be tried to achieve*=0.01511.
For multi-beam cooperation transmission, the outage capacity of each user is
Its transmission rate can be expressed as Rk=BWCOut, k；Obviously, the speed of user depends on the position distribution of user, defended
The selection of the radiation mode of star antenna, channel statistics and precoder.In this example we assume that α=0.5, radiation is increased
Beneficial matrix G matrix element gK, nIt is modeled asWherein RbeamIt is terrestrial beam covering radius (- 3dB profiles), dM, nIt is
On ground with a distance from from n-th beam center of m-th of wave beam user.Also, we consider three kinds of different customer locations point
Cloth, analyzes the performance characteristics of proposition method；
The first user location distribution is as shown in figure 3, the radiation gain matrix of collaboration user can be modeled as：
Second of user location distribution is as shown in Figure 4：The radiation gain matrix of collaboration user can be modeled as：
The third user location distribution is as shown in figure 5, the radiation gain matrix of collaboration user can be modeled as：
C. optimal group pre-coding matrix is calculated.Partial channel knowledge based on user, system power constraint and user
Rate request, constructs the optimization problem of precoding and power distribution, and the target of the optimization problem is：The fairness between user is ensured
Under the conditions of, the speed of user is maximized, is specifically expressed as follows：
P in formula (11)maxConstrained for the peak power of system, FkIt is user k rate request；
Using ZF pre-coding power optimized algorithm, orderAnd by ZF pre-coding matrix
Applied in (11), formula (11) is rewritten as
Solution formula (12), obtains optimal power allocation { pk, k=1,2 ..., K }, wherein user k distribution power is
[(GGT)-1]K, kpk；
Using precoding optimized algorithm, formula (11) is rewritten as
C in formula (13)2It is a convex function, C3It is a convex Second-order cone programming constraint, is solved using the dichotomy of iteration
Obtain optimal precoding vectors { wk, k=1,2 ..., K }.
D. precoding is carried out to the transmission signal of cooperative beam according to the optimal pre-coding matrix calculated.According to what is obtained
Optimal precoding vectorsThe transmission signal of cooperative system is：
In this example, it is assumed that system is operated in L-band, frequency fc=2GHz, the distribution of each frequency range is with a width of BW=7
×104Hz, satellite antenna is in the emission maximum antenna gain of beam centerThe product of mobile subscriber's receiver
Matter parameter isThe covering diameter of terrestrial beam is
Fig. 3 is a kind of random user distribution environment, and Fig. 6 compares the algorithm proposed in the environment and traditional frequency is answered
With the rate capability of joint Power allocation algorithm；To channeling algorithm, the distribution of each user is with a width of BW/fR, wherein, fR
It is frequency reuse, the present embodiment selection fR=4.Fig. 6 shows that two kinds of cooperation transmission methods of proposition are multiplexed than legacy frequencies
Joint Power optimization method realizes bigger rate gain, also, the increase for the peak power that can be provided with system, this
Gain further increases；For example, when system emission power is 8dBW, precoding optimization method realizes gain 87.7%, when
When transmission power is 20dBW, this gain reaches 165.5%.When we can also notice that system is operated in low transmitting power,
Zero forcing joint power optimization method performance is significantly lower than precoding optimization method.
Fig. 7 compares the method proposed in Fig. 4 position distribution environment and is multiplexed joint Power optimization method with legacy frequencies
Performance；The user location distribution feature of this environment is, the direction vector { g of collaboration userj, j=1,2 ..., K } and have best
Orthogonality；Fig. 7 also show the performance of cooperation transmission method realization than embodiment 1 more preferably, for example, being 20dBW in transmission power
When, the gain that cooperation transmission method is realized than conventional method is up to 200%；Also, the two kinds of cooperation transmission method performances proposed are all
It is very good.
Fig. 8 compares the method proposed in Fig. 5 position distribution environment and is multiplexed joint Power optimization method with legacy frequencies
Performance；Because the center of the relative cooperative beam of the position distribution of user has suitable symmetry, the direction vector of collaboration user
{gj, j=1,2 ..., K } and there is a very big correlation, the distribution of this customer location is a kind of very unfavorable environment；Fig. 8 tables
Bright, the performance of collaboration method has very big decline, especially ZF precoding joint Power optimized algorithm；But, precoding is excellent
Change method still realizes very big performance gain, for example, when transmission power is 20dBW, precoding optimization method is than tradition
Method still realize 76% performance boost.Indicate robustness of the precoding optimization method to user location distribution.
Specific embodiment described herein is only to spirit explanation for example of the invention.Technology neck belonging to of the invention
The technical staff in domain can be made various modifications or supplement to described specific embodiment or be replaced using similar mode
Generation, but without departing from the spiritual of the present invention or surmount scope defined in appended claims.
Claims (4)
1. a kind of satellite multi-beam cooperation transmission method based on partial channel knowledge, it is characterised in that：Comprise the following steps：
A. user in the cooperative beam area of coverage is grouped at random, the user in each collaboration user group is respectively from different ripples
The beam area of coverage, every group of user distributes a frequency range；
B. satellite is calculated to the link effective transmission gain of user according to known satellite subscriber channel partial information：Using Lay, this is cloudy
Shadow channel model exports the outage probability expression formula of satellite user link, is interrupted according to User Part channel information and user link
Probability demands, try to achieve the effective transmission gain of subscriber channel；Detailed process is：
The signal of every group of user is cooperated using multi-beam and sent, then terrestrial user k reception signal can be expressed as：
<mrow>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein k is k-th of user, akIt is user k fading channel coefficients；wjFor user k transmitting precoding weight vector；sjIt is
The predetermined data-signal for being sent to user；It is the direction vector of user, K represents the K of critical point station administration
Individual wave beam, or it is expressed as G=[g with radiation gain matrix1g2…gK], it reflects radiation of the satellite antenna to terrestrial user and increased
Beneficial characteristic；nkIt is to receive the Gaussian noise in signal；For the beam gain of k-th of beam center, it contains transmitting day
Line gain, path loss, the influence of receiving antenna gain and noise power, are defined as
<mrow>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>G</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>G</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mfrac>
<msub>
<mi>L</mi>
<mi>k</mi>
</msub>
<mi>&lambda;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>k</mi>
<mi>B</mi>
</msub>
<msub>
<mi>T</mi>
<mi>R</mi>
</msub>
<msub>
<mi>B</mi>
<mi>W</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
L in formula (2)kFor the distance between user k and satellite, λ is carrier frequency, kBFor Boltzmann constant,It is that user connects
Receive antenna gain,It is maximum satellite transmitting antenna gain, depending on the mode parameter α, T of transmitting antennaRAnd BWRespectively
For receiver noise temperature and the bandwidth of transmission link；
Then cooperative transmission system downlink user k reception signal Signal to Interference plus Noise Ratio can be expressed as
<mrow>
<msub>
<mi>&Gamma;</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>k</mi>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Therefore, user k channel capacity is Ck=log2(1+Γk)；
Satellite is calculated to user link outage probability：
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo><</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mo><</mo>
<msubsup>
<mi>&gamma;</mi>
<mi>k</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein Cout,kIt is the target capacity for meeting outage probability, as outage capacity；
Using this shadow fading channel model of Lay, the probability density function of fading channel is：
<mrow>
<msub>
<mi>f</mi>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mmultiscripts>
<mi>F</mi>
<mn>1</mn>
</mmultiscripts>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
<mo>;</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
1
Wherein channel parameterΩ is averaged for sight component
Power, m is Nakagami distributed constants, 2b0It is the mean power of scattering component,1F1(m,1；c0γk) it is interflow hypergeometry letter
Number, thresholdingFor：
<mrow>
<msubsup>
<mi>&gamma;</mi>
<mi>k</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>/</mo>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>k</mi>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Solution formula (4) is obtained
Pout=Φ (uk)=1-exp (- uk)Q(uk) (7)
In formula (7)Q(uk) it is a convergent unlimited item value of series,
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Pi;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mi>n</mi>
<mo>!</mo>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo>/</mo>
<mi>k</mi>
<mo>!</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
Using interior value of series in limited item number and approximant (8), obtain
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&ap;</mo>
<mn>1</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>u</mi>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</mfrac>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Satellite user partial channel knowledge includes the directional information { g of subscriber channelj, j=1,2 ..., K } and statistic informationUnder User Part channel information known case, the user based on determination is whole
Hold probability PoutIt is required that, using the dichotomy calculating formula (9) of iteration, obtain the effective transmission gain of subscriber channel
According to ukDefinition, user's outage capacity is represented by
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&ap;</mo>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mo>*</mo>
</msubsup>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>k</mi>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
Assuming that the band bandwidth B of every group of user's distributionW, then user k signal transmission rate can be expressed as Rk=BWCout,k；
C. optimal group pre-coding matrix is calculated：Asked according to partial channel knowledge, system power constraints and user rate
The optimization problem construction of optimal precoding and power distribution object function is obtained, Solve problems construction obtains optimal precoding square
Battle array；
D. precoding is carried out to the transmission signal of cooperative beam according to the optimal pre-coding matrix calculated.
2. a kind of satellite multi-beam cooperation transmission method based on partial channel knowledge according to claim 1, its feature
It is that step a detailed processes are：Multi-beam cooperation transmission method realizes that critical point station administration K wave beam, K by juncture station
N number of user is uniform-distribution with footprint of a beam, using randomly selected method by user grouping, every group of K user comes respectively
From the K different beams area of coverage.
3. a kind of satellite multi-beam cooperation transmission method based on partial channel knowledge according to claim 1, its feature
It is that step c detailed processes are：
Based on the rate request of User Part channel information, system power constraint and user, construction precoding and power distribution
Optimization problem, problem is：
<mrow>
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<munder>
<mi>min</mi>
<mrow>
<mo>{</mo>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mo>}</mo>
</mrow>
</munder>
</mtd>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>:</mo>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
P in formula (11)maxConstrained for system peak power, FkIt is user k rate request；
Using ZF pre-coding power optimized algorithm, orderAnd by ZF pre-coding matrix
Applied in (11), formula (11) is rewritten as
<mrow>
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mo>{</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>}</mo>
</mrow>
</munder>
</mtd>
<mtd>
<mi>t</mi>
</mtd>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>GG</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>:</mo>
<mfrac>
<mrow>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>&GreaterEqual;</mo>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
</mrow>
<msub>
<mi>B</mi>
<mi>w</mi>
</msub>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<mo>:</mo>
<mi>t</mi>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
Solution formula (12), obtains optimal power allocation { pk, k=1,2 ..., K }, wherein user k distribution power is [(GGT)-1]k, kpk；
Using precoding optimized algorithm, formula (11) is rewritten as
<mrow>
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mo>{</mo>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>}</mo>
</mrow>
</munder>
</mtd>
<mtd>
<mi>t</mi>
</mtd>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<mi>t</mi>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>:</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<mo>:</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msubsup>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>k</mi>
</mrow>
<mi>K</mi>
</munderover>
<mo>|</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&GreaterEqual;</mo>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
</mrow>
<msub>
<mi>B</mi>
<mi>W</mi>
</msub>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
C in formula (13)2It is a convex function, C3It is a convex Second-order cone programming constraint, is solved and obtained using the dichotomy of iteration
Optimal precoding vectors { wk, k=1,2 ..., K }.
4. a kind of satellite multi-beam cooperation transmission method based on partial channel knowledge according to claim 1, its feature
The detailed process for being step d is：According to the optimal precoding vectors obtainedThe hair of cooperative system
Penetrating signal is：
<mrow>
<mi>x</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mo>*</mo>
</msubsup>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>K</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410471326.2A CN104320174B (en) | 2014-09-16 | 2014-09-16 | A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410471326.2A CN104320174B (en) | 2014-09-16 | 2014-09-16 | A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104320174A CN104320174A (en) | 2015-01-28 |
| CN104320174B true CN104320174B (en) | 2017-08-18 |
Family
ID=52375363
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410471326.2A Expired - Fee Related CN104320174B (en) | 2014-09-16 | 2014-09-16 | A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104320174B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518961A (en) * | 2019-08-29 | 2019-11-29 | 东南大学 | Extensive MIMO satellite mobile communication method and system |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190044591A1 (en) * | 2016-02-04 | 2019-02-07 | Ntt Docomo, Inc. | Base station, user equipment, precoding matrix application method, and precoding matrix acquisition method |
| CN107425892B (en) * | 2017-06-22 | 2019-02-26 | 清华大学 | Air-ground coordination communication means and device based on multi-beam multicast |
| CN107517088B (en) * | 2017-08-24 | 2020-09-11 | 北京邮电大学 | Interference analysis method of satellite-ground integrated system based on mixed channel fading |
| CN108809390B (en) * | 2018-05-18 | 2020-06-16 | 东南大学 | Robust transmission method suitable for multicast multi-beam satellite mobile communication system |
| CN109547075B (en) * | 2018-12-21 | 2021-09-07 | 东南大学 | Method and device for grouping users of multi-beam satellite communication system |
| CN109862584B (en) * | 2018-12-30 | 2020-10-23 | 中国科学院软件研究所 | Power resource allocation method and multi-beam satellite communication system applying same |
| CN110266362B (en) * | 2019-07-04 | 2020-08-18 | 北京邮电大学 | Millimeter wave based interference suppression method for constellation multi-beam reception |
| CN111162824B (en) * | 2019-12-27 | 2021-06-29 | 中国人民解放军国防科技大学 | Multi-beam high-flux satellite communication system based on MIMO and implementation method |
| CN111510870B (en) * | 2020-04-29 | 2020-12-11 | 清华大学 | Phased array satellite multicast method and device |
| CN111416648A (en) * | 2020-05-18 | 2020-07-14 | 北京邮电大学 | Multi-beam adaptive management method and device for low-earth-orbit satellite system |
| CN113037340B (en) * | 2021-03-26 | 2022-04-05 | 东方红卫星移动通信有限公司 | Safe precoding method of multi-beam low-orbit satellite communication system |
| CN113194484B (en) * | 2021-03-31 | 2022-05-20 | 浙江大学 | Large-scale access method based on inter-satellite cooperation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101039140A (en) * | 2007-03-16 | 2007-09-19 | 北京航空航天大学 | Wireless resource management system for space-based mobile communication |
| CN101656600A (en) * | 2009-09-14 | 2010-02-24 | 西安交通大学 | Implementation method of opportunistic cooperation diversity protocol based on reinforced selective amplification-forwarding |
| CN103796319A (en) * | 2014-01-16 | 2014-05-14 | 北京大学 | Frequency reuse method for multi-beam satellite mobile communication down link |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140135016A1 (en) * | 2012-11-09 | 2014-05-15 | Qualcomm Incorporated | Method and apparatus for optimizing the frequency of autonomous search functions for discovering csg cells |
-
2014
- 2014-09-16 CN CN201410471326.2A patent/CN104320174B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101039140A (en) * | 2007-03-16 | 2007-09-19 | 北京航空航天大学 | Wireless resource management system for space-based mobile communication |
| CN101656600A (en) * | 2009-09-14 | 2010-02-24 | 西安交通大学 | Implementation method of opportunistic cooperation diversity protocol based on reinforced selective amplification-forwarding |
| CN103796319A (en) * | 2014-01-16 | 2014-05-14 | 北京大学 | Frequency reuse method for multi-beam satellite mobile communication down link |
Non-Patent Citations (2)
| Title |
|---|
| 基于部分信道信息的MIMO系统预编码技术研究;冯佳;《中国优秀硕士学位论文全文数据库信息科技辑》;20120515(第5期);正文第21、30页 * |
| 莱斯信道下跨层系统的频谱效率分析;郭丽丽，岳殿武;《大连海事大学学报》;20071130;第33卷(第4期);第19-23页 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518961A (en) * | 2019-08-29 | 2019-11-29 | 东南大学 | Extensive MIMO satellite mobile communication method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104320174A (en) | 2015-01-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104320174B (en) | A kind of satellite multi-beam cooperation transmission method based on partial channel knowledge | |
| Guidotti et al. | Geographical scheduling for multicast precoding in multi-beam satellite systems | |
| US20170026093A1 (en) | Base station and terminal for distributed array massive multiple-input and multiple-output (mimo) communication antenna system | |
| Shibata et al. | A study on cell configuration for HAPS mobile communications | |
| KR102306100B1 (en) | Apparatus for processing transmission/reception signal for interference alignment in a mu-mimo interference broadcasting channel and method thereof | |
| Chatzinotas et al. | Energy-efficient MMSE beamforming and power allocation in multibeam satellite systems | |
| CN108449150A (en) | Star based on radio map information ground disturbance coordination method, device and equipment | |
| Chatzinotas et al. | Joint precoding with flexible power constraints in multibeam satellite systems | |
| CN101702700A (en) | Method for allocating minimum power of MIMO-OFDM multi-user system based on allelism | |
| CN109547075B (en) | Method and device for grouping users of multi-beam satellite communication system | |
| Dimitrov et al. | Radio resource management techniques for high throughput satellite communication systems | |
| CN105049166B (en) | Pilot distribution method based on user's geographical location information in extensive antenna cells | |
| CN105721029A (en) | Multi-user scheduling method based on double-codebook limited feedback in 3D MU-MIMO FDD system | |
| KR101646864B1 (en) | System and method of wireless fixed access using a multiple antenna array | |
| Zhang et al. | Trajectory and power optimization for multi-UAV enabled emergency wireless communications networks | |
| Park et al. | MIMO beamforming for QoS enhancement via analog, digital and hybrid relaying | |
| CN103560856B (en) | The part method for precoding and system of satellite ground synthesis mobile communication | |
| Okuyama et al. | Flexible antenna deployment for 5G distributed Massive MIMO in low SHF bands | |
| CN101783699B (en) | Signal emission method and device of down link in multi-user MIMO (Multiple Input Multiple Output) system | |
| Kaya et al. | Dense distributed massive MIMO: Precoding and power control | |
| Christopoulos et al. | Multibeam joint precoding: Frame-based design | |
| CN107947910A (en) | Ka multi-beam satellite resource allocation methods based on artificial fish-swarm algorithm | |
| CN113938179A (en) | Joint beam forming and power control method for interference of 5G base station to satellite user | |
| Dimitrov et al. | Radio resource management for forward and return links in high throughput satellite systems | |
| Wang et al. | Beamforming and power allocation in noma-based multibeam satellite systems with outage constraint |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170818 Termination date: 20180916 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |