CN107070515A - A kind of D2D cooperation transmission methods under the conditions of rician fading channel - Google Patents
A kind of D2D cooperation transmission methods under the conditions of rician fading channel Download PDFInfo
- Publication number
- CN107070515A CN107070515A CN201710238088.4A CN201710238088A CN107070515A CN 107070515 A CN107070515 A CN 107070515A CN 201710238088 A CN201710238088 A CN 201710238088A CN 107070515 A CN107070515 A CN 107070515A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- msubsup
- 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/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
-
- 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/0417—Feedback systems
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/70—Services for machine-to-machine communication [M2M] or machine type communication [MTC]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/12—Wireless traffic scheduling
- H04W72/1215—Wireless traffic scheduling for collaboration of different radio technologies
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/12—Wireless traffic scheduling
- H04W72/1263—Mapping of traffic onto schedule, e.g. scheduled allocation or multiplexing of flows
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/50—Allocation or scheduling criteria for wireless resources
- H04W72/54—Allocation or scheduling criteria for wireless resources based on quality criteria
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Mobile Radio Communication Systems (AREA)
- Nitrogen And Oxygen Or Sulfur-Condensed Heterocyclic Ring Systems (AREA)
Abstract
The invention discloses the D2D cooperation transmission methods under the conditions of a kind of rician fading channel, comprise the following steps community user and D2D user that (1) scheduling participates in cooperation transmission；(2) estimation community user US0End subscriber D is received with D2D links1Doppler shift, obtain channel time delay coefficient correlation；(3) performance gain that two kinds of transmission modes are obtained is assessed, base station is determined the cooperation transmission mode of use；(4) user estimates respective channel information, and will feed back to base station after these information quantizations；(5) carry out data transmission under the pre-coding scheme for choosing pattern.The achievable rate of cooperation transmission method of the present invention can preferably be fitted the statistical property of real system, utilize this closed expression, can be before system actual feedback channel information, just base station is helped quickly and effectively to make judgement, then according to collaborative demand, user is only to channel condition information necessary to base station feedback, reduction feedback link transmission load.
Description
Technical field
The present invention relates to the cooperation transmission method in mobile communication multiple-input and multiple-output (MIMO) system, especially a kind of Lay
D2D cooperation transmission methods under this fading channel conditions.
Background technology
For in cell mobile communication systems, the available signal power that certain customers face is low and interfering signal power is high,
The difficult point of high data rate can not be effectively obtained, constantly having scholar to propose in recent years can be by equipment and equipment
Direct communication link is set up between (Device-to-Device, D2D) to effectively improve the data rate of user.
This D2D technologies are in existing cellular mobile communication environment so that the user positioned at cell edge can be by altogether
The frequency resource of cell other users is enjoyed, directly setting up link with neighbouring terminal carries out data transmission.Due to the two of D2D communications
The distance between individual user is very near, it is possible to use less transmission power obtains higher message transmission rate, using multiple
With the abundant excavated space resource of technology, the availability of frequency spectrum and power efficiency are improved to greatest extent.
Closely the Radio Link between D2D user largely there may be los path, and wireless channel is more likely
Rice channel model is fitted, the analysis of rician fading channel transmission performance is always the difficult point of wireless communication field, when working as sighting distance
When main signal amplitude peak is intended to ratio K →-∞ of 0, D2D link sighting distance energy and non line of sight energy, L-S distribution transformation
For its special case Rayleigh rayleigh distributeds.
The introducing of D2D links will inevitably bring significant co-channel interference (ICI), constrain system up to transmission
The lifting of speed.There are two kinds of conventional precoding processing patterns base station when sending signal, maximizes speed (maximum
Rate, MR) precoding and interference eliminate (interference cancellation, IC) precoding mode：Under IC patterns,
Base station utilizes D2D user D1The channel information of feedback does AF panel to it；Under MR patterns, base station needs to obtain community user
US0The channel information of feedback, maximization speed precoding processing is done to cell communication link.The difference of two kinds of collaboration modes is
Under IC patterns, base station not only needs community user US0The channel information of feedback, while needing D2D user D1The channel letter of feedback
Breath is, it is necessary to which extraly feedback channel information eliminates the interference that D2D user is subject to；And under MR patterns, base station only needs cell
User US0The channel condition information of feedback, thus feedback link load is obviously reduced, system delay can also be decreased.
Existing MR and IC mode evaluations scheme, which is all built upon hypothesis base station, can know the basis of ideal communication channel information
On, but in practical application scene, either time division duplex (TDD) or FDD (FDD) system, this hypothesis is all very
Difficulty is met.If system uses time division duplex (TDD), transmitting terminal can estimate uplink by reciprocity principle in base station
Path channels information and as the channel information of downlink, but up-link and downlink are not full symmetric
, thus just unavoidably there is evaluated error in the channel information that is obtained by reciprocity of base station.If system uses frequency division
Duplexing (FDD), is generally estimated downlink channel information by user terminal, then user by feedback link by downlink
Feedback of channel information to base station, but user can be limited with feedack amount in real system.User is in order to meet
The capacity requirement of feedback link, channel information is often first quantified, the information after quantization then is passed through into limited ratio
Spy feeds back to base station, and last base station recovers channel information as much as possible according to this feedback information.This scheme may result in base
Stand between the system channel information obtained and real channel information and there is a certain amount error, the size of this quantization error by
The capacity of feedback link is determined.There is the information delay problem that feedback link is brought when in addition using feedback mechanism, i.e. base station is obtained
The feedback information obtained is not the channel information at current time, but the channel information before one section is delayed.
Under the influence of above-mentioned factors, the downlink information and real channel information of base station acquisition are not fully
Matching, causes the adaptive transmission scheme performance based on ideal communication channel information to be remarkably decreased.
The content of the invention
Goal of the invention：It can be effectively applied to solve the deficiencies in the prior art there is provided one kind in actual scene, work as neighbour
Under conditions of there is line of sight between near D2D user, delay Limited Feedback channel information can only be obtained in base station end, obtained
In the case of there is the non-ideal factors such as quantization error, delay of feedback and path loss in channel information, it would be preferable to support equipment is to setting
Standby direct communication, is effectively improved adaptive in the reachable transmission rate of system and the cellular mobile communication D2D systems of the availability of frequency spectrum
Answer cooperation transmission method.
Technical scheme：D2D cooperation transmission methods under the conditions of a kind of rician fading channel of the present invention, including following step
Suddenly：
(1) scheduling participates in community user and the D2D user of cooperation transmission；
(2) estimation community user US0End subscriber D is received with D2D links1Doppler shift, obtain channel time delay phase relation
Number；
(3) performance gain that two kinds of transmission modes are obtained is assessed, base station is determined the cooperation transmission mode of use；
(4) user estimates respective channel condition information, and will feed back to base station after these information quantizations；
(5) carry out data transmission under the pre-coding scheme for choosing pattern.
The step (1) includes：Base station scheduling chooses a community user US0, its band resource is shared in the presence of straight
Up to the proximal subscribers Ds and D in footpath1D2D communicates between carrying out device-to-device, user Ds and D1Between wireless channel approximately obey Lay
This fading channel.
The step (2) includes：
(21) estimation community user US0End subscriber D is received with D2D links1Doppler shift
The average rate travel of the user terminal of cooperation, estimation community user US are participated in base station investigation cell0It is how general
Strangle frequency deviation f0, formula is：
f0=fcv0/c；
Estimate that D2D links receive end subscriber D1Doppler shift f1, formula is：
f1=fcv1/c；
Wherein, c represents the light velocity, fcRepresent carrier frequency, v0Represent community user US0Average rate travel, v1Represent
D2D links receive end subscriber D1Average rate travel.
(22) channel time delay coefficient correlation is obtained
Calculate community user US0Channel time delay correlation coefficient ρ0, formula is：
ρ0=J0(2πf0Ts)；
Calculate D2D links and receive end subscriber D1Channel time delay correlation coefficient ρ1, formula is：
ρ1=J0(2πf1Ts)；
Wherein, π represents pi, TsRepresent the length of each symbol period, J0() represents first kind zero-order Bessel letter
Number.
The step (3) includes：
(31) calculate and maximize under speed precoding mode, the closed solutions R of the reachable transmission rate of systemMR
First, community user US is calculated according to following formula0Reachable transmission rate in the case where maximizing speed precoding mode
Closed solutions
Wherein, subscript b represents cell base station, and subscript 0 represents the community user US of scheduling0, subscript s represents D2D links
Send end subscriber Ds；Lb0Represent cell base station b to community user US0Path loss, formula is：Lb0=(db0/d0)-l；Ls0Table
Show D2D link transmitting terminals s to community user US0Path loss, formula is：Ls0=(ds0/d0)-l；d0Represent transmit power ginseng
Distance is examined, l represents that decay factor, P is lost in signal propagation pathbRepresent cell base station b transmission power, PsRepresent D2D links hair
Sending end s transmission power, ρ0Represent community user US0Channel time delay coefficient correlation, B0Represent community user US0The letter of feedback
Channel state information quantizing bit number, M represents the quantity of Base Transmitter antenna.
Secondly, end subscriber D is received according to following formula computing device to device link1In the case where maximizing speed precoding mode
Up to the closed solutions of transmission rate
Wherein, subscript 1 represents the reception end subscriber D of D2D links1；Lb1Represent that cell base station b to D2D links receiving terminal is used
Family D1Path loss, formula is：Lb0=(db1/d0)-l；Ls1Represent D2D link transmitting terminals s to reception end subscriber D1Path damage
Consume, formula is：Ls1=(ds1/d0)-l。
Then, calculated and maximized under speed precoding mode according to following formula, the closed solutions R of the reachable transmission rate of systemMR：
(32) calculate interference to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC
First, community user US is calculated according to following formula0Closing for the reachable transmission rate under precoding mode is eliminated in interference
Formula solution
Secondly, end subscriber D is received according to following formula computing device to device link1Interference eliminate precoding mode under can
Up to the closed solutions of transmission rate
Wherein, ρ1Represent that device-to-device link receives end subscriber D1Channel time delay coefficient correlation；B1Represent equipment to setting
Standby terminating subscriber D1To the channel condition information bit number of base station feedback；K represents the Rice channel factor, is D2D link sighting distance energy
Amount and the ratio of non line of sight energy.
Then, calculate interference according to following formula to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC：
(33) performance gain obtained under both of which is calculated
Calculated according to following formula and maximize the performance gain R that system is obtained under speed precoding modeMR,↑：
Interference is calculated according to following formula and eliminates the performance gain R that system is obtained under precoding modeIC,↑：
Wherein, R0For community user US0Achievable rate in the case where not introducing device-to-device communication condition, according to following formula meter
Obtain：
Then the performance gain that system is obtained under the latter two transmission mode for introducing D2D communications is compared：If meeting RMR,↑
≥RIC,↑, then the transmission mode that selection maximization speed precoding mode is used as nearly base station user；Conversely, selection interference disappears
The transmission mode used except precoding mode as nearly base station user.
Wherein, function F1(x0,x1,x2, M) calculation formula be：
Wherein, x0, x1, x2, M is argument of function, and M is base station transmission antenna number；I, ζ, j are that summation operation accords with operation
The index subscript of sequence；j！Represent numeral j factorial；log2(e) it is that it is naturally right with 2 with 2 logarithms for being bottom natural number e
Number is reciprocal each other, i.e. log2(e)=1/ln2,
Function I (x0,x1,x2, m, n) calculation formula it is as follows：
Wherein, x0, x1, x2, m, n is argument of function, and q, k is the index subscript that summation operation accords with the sequence of operation,
To take out the number of combinations of q element from m different elements,For incomplete gamma functions,For gamma function, wherein, u, x is the independent variable of gamma function；
Work as x2When=1, m=0, n=1, function I (x0,x1,x2, m, n) it is reduced to：
Wherein,It is the first rank exponential integral function, k is argument of function.
Calculate R0Call function F1(x0,x1,x2, M) when, if there is x2=0 situation, it is impossible to directly calculated with above-mentioned formula,
Derive special case formula as follows：
Function F2(x0,x1,x2, M) calculation formula be：
F2(x0,x1,x2)=log2(e)(Iξ(x0,x1,x2,K)+Iξ(x0,x2,x1,K))；
Wherein,
Wherein, x0, x1, x2, K is argument of function, and K is the Rice channel factor of respective channels, and m, k, ζ transports for summation
The index subscript of the operator sequence of operation.
The step (4) includes：
(41) community user US0End subscriber D is received with D2D links1Channel estimation is carried out respectively, to obtain channel status letter
Breath and path loss；
(42) according to user and the known code book in base station, channel condition information is quantified in real time, and by after quantization
Information feedback is to base station.
The step (5) includes：
(51) base station calculates precoding vector according to the transmission mode of selection using the channel condition information of feedback；
(52) community user US will be sent to0Data carry out precoding, then carry out data transmission, D2D user utilize
Carry out data transmission with frequency resource.
Beneficial effect：Compared with prior art, the present invention has advantages below：(1) present invention can fully estimate reality
The penalty that the non-ideal factors such as path loss present in system, quantization error and delay of feedback error are caused, selection
The transmission mode being more suitable for, so as to obtain higher system transfer rate.(2) invention is applied and cellular mobile communication D2D systems
When in system, it is possible to use the closed solutions of achievable rate, before feedback link feedback instantaneous channel conditions information, estimate in advance
It is more excellent using systematic function during which kind of pre-coding scheme between base station and community user, therefore, in different system environments from
Different pre-coding schemes are adaptively selected, feedback channel can be fully utilized, implementation complexity is reduced, is preferably
System performance.
Brief description of the drawings
Fig. 1 is the system schematic of the inventive method；
Fig. 2 be under MR patterns system achievable rate with the relation curve changed with reference to signal to noise ratio；
Fig. 3 be under IC patterns system achievable rate with the relation curve changed with reference to signal to noise ratio；
Fig. 4 is user US0And user D1Achievable rate changes with the relation curve changed with reference to signal to noise ratio.
Embodiment
Explanation is further explained to technical solution of the present invention below in conjunction with accompanying drawing.
System schematic as shown in Figure 1, BSbFor base station, US0For the community user of selection, D1Presence with Ds to choose
The D2D link proximal subscribers of line of sight, user Ds and D1Between wireless channel approximately obey rician fading channel.
Solid arrow represents active link in Fig. 1, and dotted arrow represents interfering link, hb0For base station BSbAnd community user
US0Between channel condition information, hb1For base station BSbWith D2D link users D1Between channel condition information, hs1For D2D chains
Road user Ds and D1Between channel condition information, hs0For D2D link users Ds and community user US0Between channel status letter
Breath.
A kind of D2D cooperation transmission methods under the conditions of rician fading channel, comprise the following steps：
(1) scheduling participates in community user and the D2D user of cooperation transmission
Base station scheduling chooses a community user US0, its band resource is shared to the proximal subscribers Ds that there is line of sight
And D1D2D communicates between carrying out device-to-device, user Ds and D1Between wireless channel approximately obey rician fading channel.
(2) estimation community user US0End subscriber D is received with D2D links1Doppler shift, obtain channel time delay phase relation
Number
(21) estimation community user US0End subscriber D is received with D2D links1Doppler shift
The average rate travel of the user terminal of cooperation, estimation community user US are participated in base station investigation cell0It is how general
Strangle frequency deviation f0, formula is：
f0=fcv0/c；
Estimate that D2D links receive end subscriber D1Doppler shift f1, formula is：
f1=fcv1/c；
Wherein, c represents the light velocity, fcRepresent carrier frequency, v0Represent community user US0Average rate travel, v1Represent
D2D links receive end subscriber D1Average rate travel.
(22) channel time delay coefficient correlation is obtained
Calculate community user US0Channel time delay correlation coefficient ρ0, formula is：
ρ0=J0(2πf0Ts)；
Calculate D2D links and receive end subscriber D1Channel time delay correlation coefficient ρ1, formula is：
ρ1=J0(2πf1Ts)；
Wherein, π represents pi, TsRepresent the length of each symbol period, J0() represents first kind zero-order Bessel letter
Number.
(3) performance gain that two kinds of transmission modes are obtained is assessed, base station is determined the cooperation transmission mode of use
(31) calculate and maximize under speed precoding mode, the closed solutions R of the reachable transmission rate of systemMR
First, community user US is calculated according to following formula0Reachable transmission rate in the case where maximizing speed precoding mode
Closed solutions
Wherein, subscript b represents cell base station, and subscript 0 represents the community user US of scheduling0, subscript s represents D2D links
Send end subscriber Ds；
Lb0Represent cell base station b to community user US0Path loss, formula is：Lb0=(db0/d0)-l；
Ls0Represent D2D link transmitting terminals s to community user US0Path loss, formula is：Ls0=(ds0/d0)-l；
d0Represent transmit power reference distance；
L represents that decay factor is lost in signal propagation path；
PbRepresent cell base station b transmission power；
PsRepresent D2D link transmitting terminals s transmission power；
ρ0Represent community user US0Channel time delay coefficient correlation；
B0Represent community user US0The channel condition information quantizing bit number of feedback；
M represents the quantity of Base Transmitter antenna.
Secondly, end subscriber D is received according to following formula computing device to device link1In the case where maximizing speed precoding mode
Up to the closed solutions of transmission rate
Wherein, subscript 1 represents the reception end subscriber D of D2D links1；Lb1Represent that cell base station b to D2D links receiving terminal is used
Family D1Path loss, formula is Lb0=(db1/d0)-l；Ls1Represent D2D link transmitting terminals s to reception end subscriber D1Path damage
Consumption, formula is Ls1=(ds1/d0)-l。
Then, calculated and maximized under speed precoding mode according to following formula, the closed solutions R of the reachable transmission rate of systemMR：
(32) calculate interference to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC
First, community user US is calculated according to following formula0Closing for the reachable transmission rate under precoding mode is eliminated in interference
Formula solution
Secondly, end subscriber D is received according to following formula computing device to device link1Interference eliminate precoding mode under can
Up to the closed solutions of transmission rate
Wherein, ρ1Represent that device-to-device link receives end subscriber D1Channel time delay coefficient correlation；B1Represent equipment to setting
Standby terminating subscriber D1To the channel condition information bit number of base station feedback；K represents the Rice channel factor, is D2D link sighting distance energy
Amount and the ratio of non line of sight energy.
Then, calculate interference according to following formula to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC：
(33) performance gain obtained under both of which is calculated
Calculated according to following formula and maximize the performance gain R that system is obtained under speed precoding modeMR,↑：
Interference is calculated according to following formula and eliminates the performance gain R that system is obtained under precoding modeIC,↑：
Wherein, R0For achievable rates of the community user US0 in the case where not introducing device-to-device communication condition, according to following formula meter
Obtain：
Then the performance gain that system is obtained under the latter two transmission mode for introducing D2D communications is compared：If meeting RMR,↑
≥RIC,↑, then the transmission mode that selection maximization speed precoding mode is used as nearly base station user；Conversely, selection interference disappears
The transmission mode used except precoding mode as nearly base station user.
In addition, function F1(x0,x1,x2, M) calculation formula be：
Wherein, x0, x1, x2, M is argument of function, and M is base station transmission antenna number；I, ζ, j are that summation operation accords with behaviour
Make the index subscript of sequence；j！Represent numeral j factorial；log2(e) it is its nature with 2 with 2 logarithms for being bottom natural number e
Logarithm is reciprocal each other, i.e. log2(e)=1/ln2.
Function I (x0,x1,x2, m, n) calculation formula it is as follows：
Wherein, x0, x1, x2, m, n is argument of function；Q, k are the index subscript that summation operation accords with the sequence of operation；
To take out the number of combinations of q element from m different elements；For incomplete gamma functions；For gamma function, wherein, u, x is the independent variable of gamma function.
Work as x2When=1, m=0, n=1, function I (x0,x1,x2, m, n) it is reduced to：
Wherein,It is the first rank exponential integral function, k is argument of function.
As calculating R0Call function F1(x0,x1,x2, M) when, there is x2=0 situation, now can not directly be counted with above-mentioned formula
Calculate, we are that this special case derivation formula is as follows：
Function F2(x0,x1,x2, M) calculation formula be：
F2(x0,x1,x2)=log2(e)(Iξ(x0,x1,x2,K)+Iξ(x0,x2,x1,K))；
Wherein,
Wherein, x0, x1, x2, K is argument of function, and K is the Rice channel factor of respective channels；M, k, ζ transport for summation
The index subscript of the operator sequence of operation.
(4) user estimates respective channel information, and by these information by feeding back to base station after codebook quantification
(41) community user US0End subscriber D is received with D2D links1Channel estimation is carried out respectively, to obtain channel status letter
Breath and path loss；
(42) according to user and the known code book in base station, channel condition information is quantified in real time, and by after quantization
Information feedback is to base station.
Wherein, channel condition information h systems can use random vector quantization method, channel condition information be carried out limited anti-
Feedback.Channel norm | | h | | channel quality information is may be defined as, is the sign of channel quality, accounts for amount of channel feedback smaller.Generally
In feeding back channel state information, channel direction informationIt account for primary feedback amount.Also referred to as normalization channel to
Amount, average energy is normalized to 1, and receiving-transmitting sides random vector of making an appointment quantifies code book, and receiving terminal reflects channel direction information
A codebook space code book minimum with its angular separation is mapped to, code book label is fed back to by transmitting terminal by feedback link, is sent out
Sending end recovers channel direction information to quantify vector by code bookThe angle quantified between vector and actual channel vector is to have
Limit feedback quantization errorIf number of bits of feedback is B, the size of codebook space is 2B；B is smaller, required feedback band
It is wide smaller, but quantified precision is poorer；Conversely, B is bigger, required feedback bandwidth is bigger, but can improve quantified precision.
Under MR precoding modes, base station need to only obtain cell base station b to community user US from feedback link0Channel shape
State information hb0.Under IC precoding modes, base station need to obtain cell base station b to community user US simultaneously from feedback link0Letter
Channel state information hb0, and cell base station b to D2D links reception end subscriber D1Channel condition information hb1.In systems in practice,
Two kinds of precoding modes are each advantageous, and the present invention provides a kind of scheme of adaptively selected precoding mode.
(5) carry out data transmission under the pre-coding scheme for choosing pattern
(51) base station calculates precoding vector according to the transmission mode of selection using the channel condition information of feedback；
(52) community user US will be sent to0Data carry out precoding, then carry out data transmission, D2D user utilize
Carry out data transmission with frequency resource.
Embodiment：
For unified benchmark, noise variance normalization is transmitted a signal to up to the signal to noise ratio at reference distance with base station
To characterize base-station transmitting-power Pb, it is illustrated below in, assume reference distance d0=100m, signal propagation path loss decay because
Sub- l=3.5.
In the case where maximizing speed precoding mode, introduce the total reachable transmission rate of system and cell after D2D communications and use
Family US0With D2D user D1Each up to transmission rate, it is as shown in Figure 2 with signal to noise ratio situation of change.As shown in Figure 2：
Under MR patterns, by D2D user D1The curve changed using the inventive method Rate Theory result with signal to noise ratio with
The curve changed using numerical simulation rate values with signal to noise ratio can be seen that using the inventive method the calculated results with adopting
Extremely close with Numerical Simulation Results, i.e., the inventive method the calculated results being capable of preferable Curve fitting simulation result.
Under MR patterns, by community user US0The curve changed using the inventive method Rate Theory result with signal to noise ratio
With the curve that is changed with signal to noise ratio using numerical simulation rate values can be seen that using the inventive method the calculated results with
Consistent using Numerical Simulation Results, i.e., the inventive method the calculated results can be good at Curve fitting simulation result.
In addition, under MR patterns, by community user US0With D2D user D1And Rate Theory result and and rate values
As a result two kinds of curves can be seen that using the inventive method and Rate Theory result of calculation with using numerical simulation and rate values
As a result basically identical, i.e., the inventive method and Rate Theory result of calculation can be preferably fitted and rate simulating result.
Eliminated in interference under precoding mode, introduce the total reachable transmission rate of system and community user after D2D communications
US0With D2D user D1It is as shown in Figure 3 with signal to noise ratio change curve each up to transmission rate.Understand：
Under IC patterns, by D2D user D1The curve changed using the inventive method Rate Theory result with signal to noise ratio with
The curve changed using numerical simulation rate values with signal to noise ratio can be seen that using the inventive method the calculated results with adopting
Extremely close with Numerical Simulation Results, i.e., the inventive method the calculated results being capable of preferable Curve fitting simulation result.
Under IC patterns, by community user US0The curve changed using the inventive method Rate Theory result with signal to noise ratio
With the curve that is changed with signal to noise ratio using numerical simulation rate values can be seen that using the inventive method the calculated results with
Consistent using Numerical Simulation Results, i.e., the inventive method the calculated results can be good at Curve fitting simulation result.
In addition, under IC patterns, by community user US0With D2D user D1And Rate Theory result and and rate values
As a result two kinds of curves can be seen that using the inventive method and Rate Theory result of calculation with using numerical simulation and rate values
As a result basically identical, i.e., the inventive method and Rate Theory result of calculation can be preferably fitted and rate simulating result.
By the comparison of the invention described above the calculated results and Numerical Simulation Results, put forward expression formula is demonstrated in two
Can preferably Curve fitting simulation result under precoding mode.
Design parameter is：M=4, Ps=Pb/ 20, db0=50m, db1=450m, ds0=500m, ds1=10m, f0Ts=
0.1, f1Ts=0.01, B0=B1=8, K=3dB.
Under two kinds of collaboration modes, D2D user D1The achievable rate increment of acquisition, and community user US0Rate loss,
It is as shown in Figure 4 with signal to noise ratio change curve.As shown in Figure 4：
When transmit power is relatively low, when two kinds of precoding mode performances are suitable, from the point of view of system implementation complexity,
This is appropriate to the occasion simple, it is necessary to the less MR precoding modes of channel condition information amount using realizing.
When transmit power is larger, when the application condition for the channel condition information that feedback channel feeds back is small, quantify
When error and all smaller delay time error, the performance of IC precoding modes is substantially better than MR precoding modes, is adapted to pre- using IC
Coding mode.
Design parameter：M=4, Ps=Pb/ 20, db0=50m, db1=450m, ds0=500m, ds1=10m, f0Ts=0.1,
f1Ts=0.01, B0=B1=8, K=3dB.
In simulation process, it has been found that due to influenceing the factor of systematic function very many, relevant parameter has each signal source
Transmit power, the distance between each communication terminal, feedback channel bandwidth and delay of feedback information length etc..Therefore, exist
In real system, which kind of coordination strategy the selection of D2D system base-stations uses, and is the decision problem of a multidimensional.In systems in practice,
Precoding mode is sent except selection is suitable, base station and D2D can not also be kept to send the proportional relation of power of equipment, drop
The transimission power of low D2D links, because D2D communications are typical short-range communications, transmitting terminal can be sent with relatively low power
Signal transmits data.Also, in systems in practice, the feedback of channel condition information is also inevitably heavy in system opens
Pin.From simulation result, the achievable rate of cooperation transmission method of the present invention can preferably be fitted the statistics spy of real system
Property, using this closed expression, it can just help base station quickly and effectively to make before system actual feedback channel information and sentence
Certainly, then according to collaborative demand, user is only to channel condition information necessary to base station feedback, reduction feedback link transmission load.
Claims (7)
1. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel, it is characterised in that comprise the following steps：
(1) scheduling participates in community user and the D2D user of cooperation transmission；
(2) estimation community user US0End subscriber D is received with D2D links1Doppler shift, obtain channel time delay coefficient correlation；
(3) performance gain that two kinds of transmission modes are obtained is assessed, base station is determined the cooperation transmission mode of use；
(4) user estimates respective channel condition information, and will feed back to base station after these information quantizations；
(5) carry out data transmission under the pre-coding scheme for choosing pattern.
2. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 1, it is characterised in that
The step (1) includes：Base station scheduling chooses a community user US0, its band resource is shared to the neighbour that there is line of sight
Nearly user Ds and D1D2D communicates between carrying out device-to-device, user Ds and D1Between wireless channel approximately obey Lay this decline letter
Road.
3. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 1, it is characterised in that
The step (2) includes：
(21) estimation community user US0End subscriber D is received with D2D links1Doppler shift
The average rate travel of the user terminal of cooperation, estimation community user US are participated in base station investigation cell0Doppler shift
f0, formula is：
f0=fcv0/c；
Estimate that D2D links receive end subscriber D1Doppler shift f1, formula is：
f1=fcv1/c；
Wherein, c represents the light velocity, fcRepresent carrier frequency, v0Represent community user US0Average rate travel, v1Represent D2D links
Receive end subscriber D1Average rate travel；
(22) channel time delay coefficient correlation is obtained
Calculate community user US0Channel time delay correlation coefficient ρ0, formula is：
ρ0=J0(2πf0Ts)；
Calculate D2D links and receive end subscriber D1Channel time delay correlation coefficient ρ1, formula is：
ρ1=J0(2πf1Ts)；
Wherein, π represents pi, TsRepresent the length of each symbol period, J0() represents first kind zero Bessel function.
4. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 1, it is characterised in that
The step (3) includes：
(31) calculate and maximize under speed precoding mode, the closed solutions R of the reachable transmission rate of systemMR
First, community user US is calculated according to following formula0The closed solutions of reachable transmission rate in the case where maximizing speed precoding mode
<mrow>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, subscript b represents cell base station, and subscript 0 represents the community user US of scheduling0, subscript s represents the transmitting terminal of D2D links
User Ds；Lb0Represent cell base station b to community user US0Path loss, formula is Lb0=(db0/d0)-l；Ls0Represent D2D chains
Road transmitting terminal s to community user US0Path loss, formula is：Ls0=(ds0/d0)-l；d0Represent transmit power reference distance, l
Represent signal propagation path loss decay factor, PbRepresent cell base station b transmission power, PsRepresent D2D link transmitting terminals s's
Transmission power, ρ0Represent community user US0Channel time delay coefficient correlation, B0Represent community user US0The channel status letter of feedback
Quantizing bit number is ceased, M represents the quantity of Base Transmitter antenna；
Secondly, end subscriber D is received according to following formula computing device to device link1It is reachable in the case where maximizing speed precoding mode
The closed solutions of transmission rate
<mrow>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, subscript 1 represents the reception end subscriber D of D2D links1；Lb1Represent that cell base station b to D2D links receive end subscriber D1's
Path loss, formula is Lb0=(db1/d0)-l；Ls1Represent D2D link transmitting terminals s to reception end subscriber D1Path loss, it is public
Formula is：Ls1=(ds1/d0)-l；
Then, calculated and maximized under speed precoding mode according to following formula, the closed solutions R of the reachable transmission rate of systemMR：
<mrow>
<msup>
<mi>R</mi>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
(32) calculate interference to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC
First, community user US is calculated according to following formula0The closed solutions of the reachable transmission rate under precoding mode are eliminated in interference
<mrow>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<msup>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Secondly, end subscriber D is received according to following formula computing device to device link1The reachable biography under precoding mode is eliminated in interference
The closed solutions of defeated speed
<mrow>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<msubsup>
<mi>&rho;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msup>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&rho;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, ρ1Represent that device-to-device link receives end subscriber D1Channel time delay coefficient correlation；B1Represent that device-to-device is received
End subscriber D1To the channel condition information bit number of base station feedback；K represents the Rice channel factor, be D2D link sighting distance energy with
The ratio of non line of sight energy；
Then, calculate interference according to following formula to eliminate under precoding mode, the closed solutions R of the reachable transmission rate of systemIC：
<mrow>
<msup>
<mi>R</mi>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
(33) performance gain obtained under both of which is calculated
Calculated according to following formula and maximize the performance gain R that system is obtained under speed precoding modeMR,↑：
<mrow>
<msup>
<mi>R</mi>
<mrow>
<mi>M</mi>
<mi>R</mi>
<mo>,</mo>
<mo>&UpArrow;</mo>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
Interference is calculated according to following formula and eliminates the performance gain R that system is obtained under precoding modeIC,↑：
<mrow>
<msup>
<mi>R</mi>
<mrow>
<mi>I</mi>
<mi>C</mi>
<mo>,</mo>
<mo>&UpArrow;</mo>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
Wherein, R0For community user US0Achievable rate in the case where not introducing device-to-device communication condition, is calculated according to following formula
Arrive：
<mrow>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>b</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>b</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&rho;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Then the performance gain that system is obtained under the latter two transmission mode for introducing D2D communications is compared：If meeting RMR,↑≥RIC ,↑, then the transmission mode that selection maximization speed precoding mode is used as nearly base station user；Conversely, selection interference eliminates pre-
The transmission mode that coding mode is used as nearly base station user.
5. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 4, it is characterised in that：
Function F1(x0,x1,x2, M) calculation formula be：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>M</mi>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mfrac>
<msub>
<mi>I</mi>
<mn>01</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mfrac>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&zeta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>i</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mi>&zeta;</mi>
</mrow>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mi>j</mi>
</msup>
<msup>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mi>i</mi>
</msup>
<mi>j</mi>
<mo>!</mo>
</mrow>
</mfrac>
<mi>I</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>-</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, x0, x1, x2It is argument of function with M, M is respective base station transmission antenna number；I, ζ and j are that summation operation accords with behaviour
Make the index subscript of sequence；j！Represent numeral j factorial；log2(e) it is its nature with 2 with 2 logarithms for being bottom natural number e
Logarithm is reciprocal each other, i.e. log2(e)=1/ln2；
Function I (x0,x1,x2, m, n) calculation formula it is as follows：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mi>n</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mi>m</mi>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mi>q</mi>
</msup>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>m</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>q</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>q</mi>
<mo>-</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>&NotEqual;</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mi>n</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mi>m</mi>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mi>q</mi>
</msup>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>m</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>q</mi>
<mo>+</mo>
<mi>k</mi>
<mo>-</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>q</mi>
<mo>-</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>&NotEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>q</mi>
</msup>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>m</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>q</mi>
<mo>-</mo>
<mi>m</mi>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>q</mi>
<mo>-</mo>
<mi>n</mi>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, x0, x1, x2, m and n are argument of function；Q, k are the index subscript that summation operation accords with the sequence of operation；For from
The number of combinations of q element is taken out in m different elements；For incomplete gamma functions；For gamma function, wherein, u, x is the independent variable of gamma function；
Work as x2When=1, m=0, n=1, function I (x0,x1,x2, m, n) it is expressed as：
<mrow>
<msub>
<mi>I</mi>
<mn>01</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
</msup>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>&NotEqual;</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Wherein, E1(k)=∫1 ∞t-1e-ktDt is the first rank exponential integral function, and k is argument of function；
Calculate R0Call function F1(x0,x1,x2, M) when, if there is x2=0 situation, it is impossible to directly calculated with above-mentioned formula, is derived
Special case formula is as follows：
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>M</mi>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
</msup>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&zeta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>i</mi>
</munderover>
<mfrac>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mi>&zeta;</mi>
</msup>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>)</mo>
<mo>!</mo>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mi>I</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mo>;</mo>
</mrow>
Function F2(x0,x1,x2, M) calculation formula be：
F2(x0,x1,x2)=log2(e)(Iξ(x0,x1,x2,K)+Iξ(x0,x2,x1,K))；
Wherein,
<mrow>
<msub>
<mi>I</mi>
<mi>&xi;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</msup>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&zeta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mrow>
<mi>m</mi>
<mo>!</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>)</mo>
</mrow>
<mo>!</mo>
</mrow>
</mfrac>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
<mo>,</mo>
<mi>&zeta;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
3
<mrow>
<msub>
<mi>I</mi>
<mi>&xi;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</msup>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&zeta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mrow>
<mi>m</mi>
<mo>!</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
<mo>-</mo>
<mi>&zeta;</mi>
<mo>)</mo>
</mrow>
<mo>!</mo>
</mrow>
</mfrac>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
<mo>,</mo>
<mi>&zeta;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, x0, x1, x2It is argument of function with K, K is the Rice channel factor of respective channels；M, k and ζ are summation operation
Accord with the index subscript of the sequence of operation.
6. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 1, it is characterised in that
The step (4) includes：
(41) community user US0End subscriber D is received with D2D links1Carry out channel estimation respectively, with obtain channel condition information and
Path loss；
(42) according to user and the known code book in base station, channel condition information is quantified in real time, and by the channel after quantization
Status information feedback is to base station.
7. the D2D cooperation transmission methods under the conditions of a kind of rician fading channel according to claim 1, it is characterised in that
The step (5) includes：
(51) base station calculates precoding vector according to the transmission mode of selection using the channel condition information of feedback；
(52) community user US will be sent to0Data carry out precoding, then carry out data transmission, D2D user utilize with frequency provide
Source carries out data transmission.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710238088.4A CN107070515B (en) | 2017-04-13 | 2017-04-13 | D2D cooperative transmission method under Rice fading channel condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710238088.4A CN107070515B (en) | 2017-04-13 | 2017-04-13 | D2D cooperative transmission method under Rice fading channel condition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107070515A true CN107070515A (en) | 2017-08-18 |
| CN107070515B CN107070515B (en) | 2020-04-24 |
Family
ID=59601958
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710238088.4A Active CN107070515B (en) | 2017-04-13 | 2017-04-13 | D2D cooperative transmission method under Rice fading channel condition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107070515B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111614397A (en) * | 2020-05-22 | 2020-09-01 | 南京邮电大学 | High-altitude unmanned platform communication system and multi-antenna wireless transmission method |
| CN112601286A (en) * | 2020-11-12 | 2021-04-02 | 南京邮电大学 | User scheduling method based on channel estimation error |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2171949B1 (en) * | 2007-07-10 | 2011-10-26 | QUALCOMM Incorporated | Methods and apparatus for successive interference cancellation based on rate capping in peer-to-peer networks |
| CN103095423A (en) * | 2013-01-31 | 2013-05-08 | 西安电子科技大学 | Multi-user cooperative transmission method based on dimension to dimension (D2D) inter-group communication |
| CN103986558A (en) * | 2014-05-26 | 2014-08-13 | 东南大学 | Self-adaptation cooperative transmission method in cellular mobile communication D2D system |
| CN105071843A (en) * | 2015-07-29 | 2015-11-18 | 东南大学 | Large-scale MIMO system low-complexity polynomial expansion matrix inversion method and application thereof |
| CN105684545A (en) * | 2013-10-31 | 2016-06-15 | Lg电子株式会社 | Method by which terminal transmits and receives signal to/from base station by using device-to-device communication in wireless communication system, and device for same |
| US20170013495A1 (en) * | 2015-07-10 | 2017-01-12 | Lg Electronics Inc. | Method and apparatus for an input data processing via a local computing or offloading based on power harvesting in a wireless communication system |
| CN106550409A (en) * | 2016-12-12 | 2017-03-29 | 西安电子科技大学 | Compound relaying motivational techniques in coordination downlink communication |
-
2017
- 2017-04-13 CN CN201710238088.4A patent/CN107070515B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2171949B1 (en) * | 2007-07-10 | 2011-10-26 | QUALCOMM Incorporated | Methods and apparatus for successive interference cancellation based on rate capping in peer-to-peer networks |
| CN103095423A (en) * | 2013-01-31 | 2013-05-08 | 西安电子科技大学 | Multi-user cooperative transmission method based on dimension to dimension (D2D) inter-group communication |
| CN105684545A (en) * | 2013-10-31 | 2016-06-15 | Lg电子株式会社 | Method by which terminal transmits and receives signal to/from base station by using device-to-device communication in wireless communication system, and device for same |
| CN103986558A (en) * | 2014-05-26 | 2014-08-13 | 东南大学 | Self-adaptation cooperative transmission method in cellular mobile communication D2D system |
| US20170013495A1 (en) * | 2015-07-10 | 2017-01-12 | Lg Electronics Inc. | Method and apparatus for an input data processing via a local computing or offloading based on power harvesting in a wireless communication system |
| CN105071843A (en) * | 2015-07-29 | 2015-11-18 | 东南大学 | Large-scale MIMO system low-complexity polynomial expansion matrix inversion method and application thereof |
| CN106550409A (en) * | 2016-12-12 | 2017-03-29 | 西安电子科技大学 | Compound relaying motivational techniques in coordination downlink communication |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111614397A (en) * | 2020-05-22 | 2020-09-01 | 南京邮电大学 | High-altitude unmanned platform communication system and multi-antenna wireless transmission method |
| CN111614397B (en) * | 2020-05-22 | 2022-07-26 | 南京邮电大学 | High-altitude unmanned platform communication system and multi-antenna wireless transmission method |
| CN112601286A (en) * | 2020-11-12 | 2021-04-02 | 南京邮电大学 | User scheduling method based on channel estimation error |
| CN112601286B (en) * | 2020-11-12 | 2023-05-26 | 南京邮电大学 | User scheduling method based on channel estimation error |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107070515B (en) | 2020-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103986558B (en) | Adaptive cooperation transmission method in a kind of cellular mobile communication D2D systems | |
| CN102484543B (en) | Pilot-based Sinr Estimation For Mimo Systems | |
| CN101378280B (en) | Multi-input multi-output system based on antenna selection and signal processing method thereof | |
| US8036257B2 (en) | Method of determining at least one transmit mode parameter for a multiple-input multiple-output system | |
| CN105338609B (en) | Multiaerial system high energy efficiency dynamic power allocation method | |
| CN102983934A (en) | Method and device for linear precoding in multi-user multiple input multiple output system | |
| CN103312390B (en) | Based on the method for precoding of interference alignment, emitter and equipment | |
| CN101867462A (en) | Multi-base station cooperation linear precoding method based on minimum total bit error rate | |
| WO2011053501A1 (en) | Differential closed-loop transmission feedback in wireless communication systems | |
| CN101459495B (en) | Multi-user MIMO communication system, method and apparatus for transmitting signal | |
| CN102055563A (en) | Adaptive joint linear precoding method applicable to multi-base station coordination | |
| CN102545987A (en) | Multicell self-adaption cooperative transmission method on basis of delayed feedback | |
| CN107707284B (en) | Mixed precoding method based on channel statistic codebook quantization feedback | |
| CN101783776A (en) | Precoding feedback method, system, user equipment and base station | |
| CN102811491A (en) | Joint distribution method of limited feedback number of bits in multi-point cooperation system | |
| KR20060119144A (en) | Multiple users multiple input multiple output antenna commucation apparatus and method thereof | |
| Ekpenyong et al. | Feedback constraints for adaptive transmission | |
| CN102752037A (en) | Multi-antenna relay pre-code robust construction method under delayed feedback | |
| CN107070515A (en) | A kind of D2D cooperation transmission methods under the conditions of rician fading channel | |
| CN105978612A (en) | Pilot pollution-eliminated double-layer pre-coding design method in massive MIMO system | |
| Tien Ngoc et al. | Harvested energy and spectral efficiency trade-offs in multicell MIMO wireless networks | |
| CN107579762B (en) | Multi-cell cooperative precoding method based on quantification and channel information statistics | |
| CN1910830B (en) | Method for transmitting/receiving signal in mimo system | |
| CN101989867B (en) | Cooperative communication method and system, base station and mobile terminal device | |
| CN102056311B (en) | Base station scheduling method for improving MU-MIMO performance and device |
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 |