CN105554808B - A kind of D2D in Cellular Networks is to resource allocation methods and device - Google Patents

Abstract

The embodiment of the invention discloses the D2D in a kind of Cellular Networks to resource allocation methods and device, this method are as follows: collect all phone users and D2D pairs of channel state information in base station range, it determines each channel gain, then calculates separately throughput gain of each D2D to the uplink and downlink channel for being multiplexed each phone user respectively again；For each D2D pairs, the comparison of throughput gain based on its uplink and downlink channel for being multiplexed each phone user, determine D2D pairs of the multiplexer mode, determine each D2D to corresponding multiple channels again, finally determine that each D2D and is each D2D to distribution channel and transmission power to the transmission power on multiple channels corresponding to it.The present invention compares selection multiplexer mode according to the analysis that each throughput gain of D2D user couple carries out, and increases D2D significantly resource reusable in system, throughput of system greatly increases.

Description

A kind of D2D in Cellular Networks is to resource allocation methods and device

Technical field

The present invention is more particularly directed to the D2D in a kind of Cellular Networks to resource allocation methods and device.

Background technique

In recent years, with the fast development of mobile multi-media service, requirement of the user to transmission rate is higher and higher, has The band resource of limit can not meet the growth requirement of wireless communication system, and D2D (Device-to-Device) is communicated as a result, Technology comes into being and causes domestic and international industry, the strong interest of academia.D2D communication is a kind of control in cellular system Under, allow adjacent terminal end user to carry out the new technology of direct communication by reuse pattern resource.D2D under cellular network covering is logical Letter can be improved the spectrum efficiency of network, increase throughput of system, expands cell coverage area, mitigates base station burden, reduces and pass Defeated time delay, equipment power dissipation meet user experience demand.

But on the other hand, since the D2D communication technology has been multiplexed the frequency spectrum resource of cellular cell, D2D user is to (referred to as D2D to) and phone user communicate on the same channel, co-channel interference can be generated between the two, in the case where serious interference The communication quality that the two can seriously be reduced, reduces the performance of system.Therefore, rationally the channel resource of effective D2D user couple divides It is the critical issue of D2D Communication Studies under cellular network with scheme.

In the prior art, there are mainly two types of the resource multiplex mode of D2D communication in cellular networks: up channel is multiplexed mould Formula, such as Chinese patent CN 103079262A；And down channel multiplexer mode.It is apparent that both aforementioned multiplexer modes are all It is that all D2D in base station range is allowed to be multiplexed up channel to all multiplexing down channels or all, D2D is to reusable Resource is few, causes D2D small to the handling capacity of system.

Summary of the invention

The embodiment of the invention discloses the D2D in a kind of Cellular Networks to resource allocation methods and device, existing for solving Resource Allocation Formula in few to reusable resource, D2D small to the handling capacity of the system problem of D2D.Technical solution is as follows:

In a first aspect, the embodiment of the invention provides the D2D in a kind of Cellular Networks to resource allocation methods, comprising steps of

All phone users and D2D pairs of channel state information in base station range are collected, determines each channel gain, institute Stating each channel gain includes that channel gain, each D2D pairs of channel between base station between each phone user and base station increase Beneficial, each D2D to respectively between each phone user channel gain and each D2D to internal channel gain；

According to each channel gain, calculate separately each D2D to the up channel for being multiplexed each phone user respectively and The throughput gain of down channel；

For each D2D pairs, throughput gain and multiplexing based on its up channel for being multiplexed each phone user are each The comparison of the throughput gain of the down channel of phone user determines D2D pairs of the multiplexer mode are as follows: be individually multiplexed up channel It carries out communication or individually multiplexing down channel is communicated；

According to each D2D to the throughput gain for the up channel and down channel for being multiplexed each phone user respectively and All D2D pairs of multiplexer mode determines each D2D to corresponding multiple channels；

According to D2D pairs of multiplexer mode and each D2D to corresponding multiple channels, determine each D2D to right in its institute The transmission power on multiple channels answered, and be each D2D to the multiple channels and more corresponding to it distributed corresponding to it Corresponding transmission power on a channel.

Preferably, described to calculate separately each D2D to the up channel and down channel for being multiplexed each phone user respectively Throughput gain used by formula include:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jMultiplexing Phone user C_iDown channel throughput gain,Indicate D2D pairs of maximum transmission power, P_BIndicate the transmitting function of base station Rate, P_CIndicate the transmission power of phone user, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D pairs D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jIt is interior The channel gain in portion, channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is path Loss index, d are distance of the transmitting terminal to receiving end, N₀Indicate noise figure.

Preferably, described to be directed to each D2D pairs, each throughput gain and multiplexing downlink of up channel are multiplexed based on it The comparison of each throughput gain of channel determines D2D pairs of the multiplexer mode are as follows: individually multiplexing up channel carry out communication or Individually multiplexing down channel is communicated, comprising:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed up channel Average throughput flow gain, in turn,

Each throughput gain of down channel is multiplexed respectively compared with the average throughput gain, and statistical value Greater than the quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain, if quantity is more than predetermined integers, Determine that the D2D communicates independent multiplexing down channel, otherwise the D2D communicates independent multiplexing up channel；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is phone user in the base station range Quantity, N is D2D pairs of quantity in the base station range.

Preferably, the gulping down to the up channel and down channel for being multiplexed each phone user respectively according to each D2D Flow gain and all D2D pairs of multiplexer modes are spat, determines each D2D to corresponding multiple channels, comprising:

A1: based on all D2D pairs of multiplexer mode, in uplink throughput gain matrix Ω^uIn, it will individually be multiplexed downlink letter The D2D that road is communicated is to corresponding elementIt is updated to any value less than or equal to zero；In downlink throughput capacity gain matrix Ω^dIn, the D2D that independent multiplexing up channel is communicated is to corresponding elementIt is updated to any less than or equal to zero Value；

A2: initialization matrix X=[x_ji]_N×M=0, define Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N }；

A3: it finds out and meets be expert at j ∈ W in throughput gain matrix Ω_j, column i ∈ Ф_iGreatest member, determine square Battle array X in the corresponding element x of row and column where the corresponding greatest member_ji=1；

A4: Ф is updated_i=Ф_iI, W_j=W_jJ,

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jMember Plain T_jiGreater than 0, return step A3；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiIt is greater than 0, update W_j={ 1,2 ..., N }, return step A3；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, It then determines output matrix X, and determines each D2D to corresponding multiple channels according to matrix X；

Wherein, the phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station D2D in coverage area is to D_jQuantity be N number of, j ∈ { 1,2 ..., N }, M > N, throughput gain matrix Ω=[T_ji]_N×M, square Battle array X characterizes up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when matrix X characterization up channel multiplexing Situation matrix X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to multiplexing phone user's up channel Throughput gain matrix Ω^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel throughput gain When matrix X characterizes down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D To the throughput gain matrix Ω of multiplexing phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUnder The throughput gain of row channelCharacterize phone user C_iUp channel distribute to D2D to D_j,Table Levy phone user C_iUp channel be not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D To D_j,Characterize phone user C_iDown channel be not yet assigned to D2D to D_j。

Preferably, it is described according to D2D pairs of multiplexer mode and each D2D to corresponding multiple channels, determine each D2D To the transmission power on multiple channels corresponding to it, comprising:

Based on D2D pairs of multiplexer mode, determine that the Signal to Interference plus Noise Ratio calculation formula of phone user and D2D pairs is respectively as follows:

Meanwhile determining the formula of the maximum of D2D system total rate and constraint condition are as follows:

s.t.

IfOtherwiseIfOtherwise

Based on each D2D to corresponding multiple channels, by the phase of the maximum of the D2D system total rate and constraint condition Answer simplified formula are as follows:

Lagrangian is converted by the total function of rate of the maximum of D2D system and the respective formula of constraint condition:

Its dual function are as follows:

The operational formula of power optimized solution is acquired according to KKT condition are as follows:

Based on the operational formula of power optimized solution, each power optimized solution P is found out_ji ^*, based on each power optimized acquired Solve P_ji ^*, it is each D2D to the transmission power on the multiple channels distributed corresponding to it；Wherein, P_ji ^*It is D2D to D_jIn honeycomb User C_iTransmission power on upstream or downstream channel；

Wherein,

Phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N },For phone user C_iThe letter of up channel dry make an uproar Than,For phone user C_iDown channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iUp channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iThe Signal to Interference plus Noise Ratio of down channel,Indicate D2D to D_jTo honeycomb User C_iUp channel multiplexed situation, if phone user C_iUp channel by D2D to D_jMultiplexing, thenOtherwiseIndicate D2D to D_jTo phone user C_iDown channel multiplexed situation, if phone user C_iDown channel quilt D2D is to D_jMultiplexing, thenOtherwiseChannel gain g=G β d^-α, G therein is path loss constant, and β represents letter Road decline constant, α are path loss index, and d is distance of the transmitting terminal to receiving end；ParameterIndicate D2D to D_jIt is used in honeycomb Family C_iTransmission power on up channel,Indicate D2D to D_jIn phone user C_iTransmission power on down channel, P_BIt indicates The transmission power of base station, P_CIndicate phone user C_iTransmission power, g_BiIndicate phone user C_iChannel between base station increases Benefit, g_BjIndicate D2D to D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel gain, N₀Indicate noise figure；max R_D2DIndicate all D2D to the maximum of the system of composition Rate, γ^uIndicate the uplink minimum signal-to-noise ratio of preset phone user, γ^dIndicate the downlink minimum noise of preset phone user Than,Indicate D2D pairs of maximum transmission power.

P_ji ^*D2D is characterized to D_jIn phone user C_iPower optimized solution on up channelOr the power on down channel Optimal solutionWork as P_ji ^*Characterize the power optimized solution on up channelWhen, P_jiD2D is characterized to D_jIn C_iOn up channel Transmission powerγ characterizes the uplink minimum signal-to-noise ratio γ of preset phone user^u；Work as P_ji ^*Characterize the power on down channel Optimal solutionWhen, P_jiD2D is characterized to D_jIn C_iTransmission power on down channelγ is characterized under preset phone user Row minimum signal-to-noise ratio γ^d；WithRespectively D2D is to D_jIn phone user C_iTransmitting function on up channel and down channel Rate；μ_jWithFor two multiplication factors in Lagrangian, μ_j ^*WithTo pass through pair acquired according to gradient search algorithm The optimal solution of even function；maxAnd maxRespectively all D2D are to composition in up channel and down channel The total rate of maximum of system.

Preferably, the gradient search algorithm comprising steps of

B1: initialization multiplication factorμ_j(0) and step parameter δ,

Wherein, i ∈ { 1,2 ..., M }, j ∈ { 1,2 ..., N }, and corresponding element x_ji≠0；

B2: it updates

B3: t=t+1, return step B2 are updated；Until(μ_i(0),μ_i(1),…, μ_i(t+1)) it restrains, enablesμ_j ^*=μ_i(t+1)。

Wherein, t indicates the number of iterations of the gradient search algorithm,And μ_j(t) when respectively indicating the t times iteration Multiplication factorμ_j, element x_jiCharacterize elementOr elementWork as P_ji ^*Characterize the power optimized solution on up channel When, element x_jiCharacterize elementWork as P_ji ^*Characterize the power optimized solution on down channelWhen, element x_jiCharacterize element

On the other hand, the embodiment of the present invention also provides the D2D in a kind of Cellular Networks to resource allocation device, comprising:

Channel gain determining module, for collecting all phone users and D2D pairs of channel status in base station range Information, determines each channel gain, and each channel gain includes channel gain between each phone user and base station, each Channel gain, each D2D between D2D pairs and base station are to channel gain respectively between each phone user and each D2D is to internal channel gain；

Throughput gain computing module, for it is every to multiplexing respectively to calculate separately each D2D according to each channel gain The up channel of a phone user and the throughput gain of down channel；

Multiplexer mode determining module, for being directed to each D2D pairs, based on its up channel for being multiplexed each phone user The comparison of the throughput gain of the down channel of throughput gain and each phone user of multiplexing, determines D2D pairs of the multiplexing mould Formula are as follows: individually multiplexing up channel carries out communication or individually multiplexing down channel is communicated；

Channel determination module, for being believed according to each D2D the up channel and downlink for being multiplexed each phone user respectively The throughput gain in road and all D2D pairs of multiplexer modes determine each D2D to corresponding multiple channels；

Resource distribution module is determined for the multiplexer mode and each D2D according to D2D pairs to corresponding multiple channels Each D2D to the transmission power on multiple channels corresponding to it, and for each D2D to distribute its corresponding to multiple channels And corresponding transmission power on multiple channels corresponding to it.

Preferably, each D2D is calculated separately in the throughput gain computing module to being multiplexed each phone user respectively Up channel and down channel throughput gain used by formula include:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jMultiplexing Phone user C_iDown channel throughput gain,Indicate D2D pairs of maximum transmission power, P_BIndicate the transmitting function of base station Rate, P_CIndicate the transmission power of phone user, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D pairs D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jIt is interior The channel gain in portion, channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is path Loss index, d are distance of the transmitting terminal to receiving end, N₀Indicate noise figure.

Preferably, the multiplexer mode determining module is specifically used for:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed up channel Average throughput flow gain, in turn,

Each throughput gain of down channel is multiplexed respectively compared with the average throughput gain, and statistical value Greater than the quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain, if quantity is more than predetermined integers, Determine that the D2D communicates independent multiplexing down channel, otherwise the D2D communicates independent multiplexing up channel；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is phone user in the base station range Quantity, N is D2D pairs of quantity in the base station range.

Preferably, the channel determination module includes:

Throughput gain updating unit, for the multiplexer mode based on all D2D pairs, in uplink throughput gain matrix Ω^uIn, the D2D that independent multiplexing down channel is communicated is to corresponding elementIt is updated to any less than or equal to zero Value；In downlink throughput capacity gain matrix Ω^dIn, the D2D that independent multiplexing up channel is communicated is to corresponding element It is updated to any value less than or equal to zero；

Initial cell, for initializing matrix X=[x_ji]_N×M=0, define Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N}；

Element x_jiDetermination unit meets be expert at j ∈ W for finding out in throughput gain matrix Ω_j, column i ∈ Ф_i Greatest member, determine in matrix X with the corresponding element x of row and column where the corresponding greatest member_ji=1；

Channels determining unit, for updating Ф_i=Ф_iI, W_j=W_jJ,

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jMember Plain T_jiGreater than 0, element x is returned_jiDetermination unit；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiIt is greater than 0, update W_j={ 1,2 ..., N } returns to element x_jiDetermination unit；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, It then determines output matrix X, and determines each D2D to corresponding multiple channels according to matrix X；

Wherein, the phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station D2D in coverage area is to D_jQuantity be N number of, j ∈ { 1,2 ..., N }, M > N, throughput gain matrix Ω=[T_ji]_N×M, square Battle array X characterizes up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when matrix X characterization up channel multiplexing Situation matrix X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to multiplexing phone user's up channel Throughput gain matrix Ω^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel throughput gainWhen Matrix X characterizes down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D pairs It is multiplexed the throughput gain matrix Ω of phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iDownlink The throughput gain of channelCharacterize phone user C_iUp channel distribute to D2D to D_j,Characterize bee Nest user C_iUp channel be not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D to D_j,Characterize phone user C_iDown channel be not yet assigned to D2D to D_j。

Preferably, the resource distribution module includes:

Transmission power determination unit determines phone user and D2D pairs of letter is dry makes an uproar for the multiplexer mode based on D2D pairs It is respectively as follows: than calculation formula

Meanwhile determining the formula of the maximum of D2D system total rate and constraint condition are as follows:

s.t.

IfOtherwiseIfOtherwise

Based on each D2D to corresponding multiple channels, by the phase of the maximum of the D2D system total rate and constraint condition Answer simplified formula are as follows:

Lagrangian is converted by the total function of rate of the maximum of D2D system and the respective formula of constraint condition:

Its dual function are as follows:

The operational formula of power optimized solution is acquired according to KKT condition are as follows:

Based on the operational formula of power optimized solution, each power optimized solution P is found out_ji ^*, based on each power optimized acquired Solve P_ji ^*, it is each D2D to the transmission power on the multiple channels distributed corresponding to it；Wherein, P_ji ^*It is D2D to D_jIn honeycomb User C_iTransmission power on upstream or downstream channel；

Wherein,

Phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N },For phone user C_iThe letter of up channel dry make an uproar Than,For phone user C_iDown channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iUplink letter The Signal to Interference plus Noise Ratio in road,It is D2D to D_jIt is multiplexed phone user C_iThe Signal to Interference plus Noise Ratio of down channel,Indicate D2D to D_jTo bee Nest user C_iUp channel multiplexed situation, if phone user C_iUp channel by D2D to D_jMultiplexing, thenOtherwiseIndicate D2D to D_jTo phone user C_iDown channel multiplexed situation, if phone user C_iDown channel quilt D2D is to D_jMultiplexing, thenOtherwiseChannel gain g=G β d^-α, G therein is path loss constant, and β represents letter Road decline constant, α are path loss index, and d is distance of the transmitting terminal to receiving end；ParameterIndicate D2D to D_jIt is used in honeycomb Family C_iTransmission power on up channel,Indicate D2D to D_jIn phone user C_iTransmission power on down channel, P_BIt indicates The transmission power of base station, P_CIndicate phone user C_iTransmission power, g_BiIndicate phone user C_iChannel between base station increases Benefit, g_BjIndicate D2D to D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel gain, N₀Indicate noise figure；max R_D2DIndicate all D2D to the maximum of the system of composition Rate, γ^uIndicate the uplink minimum signal-to-noise ratio of preset phone user, γ^dIndicate the downlink minimum noise of preset phone user Than,Indicate D2D pairs of maximum transmission power；

P_ji ^*D2D is characterized to D_jIn phone user C_iPower optimized solution on up channelOr the power on down channel Optimal solutionWork as P_ji ^*Characterize the power optimized solution on up channelWhen, P_jiD2D is characterized to D_jIn C_iOn up channel Transmission powerγ characterizes the uplink minimum signal-to-noise ratio γ of preset phone user^u；Work as P_ji ^*Characterize the power on down channel Optimal solutionWhen, P_jiD2D is characterized to D_jIn C_iTransmission power on down channelγ is characterized under preset phone user Row minimum signal-to-noise ratio γ^d；WithRespectively D2D is to D_jIn phone user C_iTransmitting function on up channel and down channel Rate；μ_jWithFor two multiplication factors in Lagrangian, μ_j ^*WithTo pass through pair acquired according to gradient search algorithm The optimal solution of even function；maxAnd maxSystem of all D2D to composition respectively in up channel and down channel The total rate of maximum；

Channel and power distributing unit, for for each D2D to distribute its corresponding to multiple channels and its corresponding to Multiple channels on corresponding transmission power.

Preferably, realize that device used by the gradient search algorithm includes:

Initial submodule: for initializing multiplication factorμ_j(0) and step parameter δ,

Wherein, i ∈ { 1,2 ..., M }, j ∈ { 1,2 ..., N }, and corresponding element x_ji≠0；

Update submodule: for updating

It determines submodule: for updating t=t+1, updated t being returned and updates submodule；Until(μ_i(0),μ_i(1),…,μ_i(t+1)) it restrains, enablesμ_j ^*=μ_i (t+1)；

Wherein, t indicates the number of iterations of the gradient search algorithm,And μ_j(t) when respectively indicating the t times iteration Multiplication factorμ_j, element x_jiCharacterize elementOr elementWork as P_ji ^*Characterize the power optimized solution on up channel When, element x_jiCharacterize elementWork as P_ji ^*Characterize the power optimized solution on down channelWhen, element x_jiCharacterize element

D2D in a kind of Cellular Networks provided in an embodiment of the present invention uses resource allocation methods and device using a D2D Family to the mode for being multiplexed multiple phone user's resources, meanwhile, according to handling capacity caused by the specific space diversity of D2D user couple Gain carries out analysis and compares the best multiplexer mode of selection, increases D2D significantly resource reusable in system, D2D is to system Handling capacity greatly increase, and then the problem of avoid serious interference between D2D couples and phone user, and use more D2D Family is docked in network, improves D2D user to the overall performance of the system of composition.

Detailed description of the invention

In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.

Fig. 1 is the applicable scene figure of the embodiment of the present invention；

Fig. 2 is flow diagram of the D2D in a kind of Cellular Networks provided in an embodiment of the present invention to resource allocation methods；

Fig. 3 is structural schematic diagram of the D2D in a kind of Cellular Networks provided in an embodiment of the present invention to resource allocation device；

Fig. 4 is the column diagram of different throughput of system provided in an embodiment of the present invention；

Fig. 5 is the line chart of the relationship between D2D handling capacity and phone user's quantity provided in an embodiment of the present invention；

Fig. 6 is the line chart of D2D handling capacity provided in an embodiment of the present invention and D2D to the relationship between quantity；

Fig. 7 be D2D handling capacity provided in an embodiment of the present invention at a distance from two users of each D2D centering between pass The line chart of system.

Specific embodiment

Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.

It should be noted that there are three types of the Channel distribution mode of D2D communication in cellular networks: one, a D2D user is multiple With a resource block, such as Chinese patent CN 103686743A, such channel distribution form has that the availability of frequency spectrum is low； Two, multiple D2D users are multiplexed a resource block, and such channel distribution form can generate serious interference；Three, a D2D user Multiple resource blocks are multiplexed, used in the embodiment of the present invention is exactly the third channel multiplexing mode as described above.Described herein Resource block refers to channel, up channel and down channel including phone user.

The method of the embodiment of the present invention can be applied in the server of base station, set it is of course also possible to apply with base station In standby associated device server.The application scenarios of the embodiment of the present invention as shown in Figure 1, in the coverage area of base station BS, There are three phone user C₁、C₁、C₃And two D2D are to D₁、D₂；In figure, straight solid line with the arrow indicates upstream data signals, Straight dotted line with the arrow indicates that downlink data signal, waveform solid line with the arrow indicate uplink interference signal, wave with the arrow Shape wave dotted line indicates downlink interference signal.

D2D in a kind of Cellular Networks provided by the embodiment of the present invention may include steps of resource allocation methods:

S1: collecting all phone users and D2D pairs of channel state information in base station range, determines that each channel increases Benefit, each channel gain include channel gain, each D2D pairs of letter between base station between each phone user and base station Road gain, each D2D to respectively between each phone user channel gain and each D2D to internal channel gain；

S2: according to each channel gain, each D2D is calculated separately to the up channel for being multiplexed each phone user respectively With the throughput gain of down channel；

In the embodiment of the present invention, each D2D is calculated to the up channel and down channel for being multiplexed each phone user respectively Throughput gain used by formula can be with are as follows:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jMultiplexing Phone user C_iDown channel throughput gain, it should explanation, D described herein_jAnd C_iBe respectively D2D pairs and The number of phone user；

Indicate D2D pairs of maximum transmission power, P_BIndicate the transmission power of base station, P_CIndicate the transmitting function of phone user Rate, it should which the transmission power of explanation, D2D pairs of maximum transmission power, the transmission power of base station and phone user is all It can be searched in communication standard determining, can be used directly, do not need to carry out additional calculating；

G in formula_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D to D_jBetween base station Channel gain, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel increases Benefit, the calculation formula of channel gain are g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is road Diameter loss index, d are distance of the transmitting terminal to receiving end, and the unit of d is km；N₀Indicate noise figure, above-mentioned N₀, G, β, α and d This four parameters just belong to what collection in step S1 obtained.

S3: being directed to each D2D pairs, and throughput gain and multiplexing based on its up channel for being multiplexed each phone user are every The comparison of the throughput gain of the down channel of a phone user determines D2D pairs of the multiplexer mode are as follows: individually multiplexing uplink letter Road carries out communication or individually multiplexing down channel is communicated；

It should be noted that some D2D make it be multiplexed any phone user's to due to self-position diversity Channel, throughput gain are both greater than other D2D pairs, will lead to base station in this way for more channel distribution to the same D2D It is right, and other D2D can not access network to dividing because throughput gain is smaller less than channel.

It is above-mentioned to be directed to each D2D pairs in the embodiment of the present invention in order to dock more D2D into network, based on its multiplexing The throughput gain of the down channel of the throughput gain and each phone user of multiplexing of the up channel of each phone user Comparison, determines D2D pairs of the multiplexer mode are as follows: individually multiplexing up channel carries out communication or individually multiplexing down channel is led to Believe, may include:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed up channel Average throughput flow gain, in turn, be multiplexed each throughput gain of down channel respectively with the average throughput flow gain Compare, and statistical value is greater than the quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain, if quantity More than predetermined integers, determine that the D2D communicates independent multiplexing down channel, otherwise the D2D is to independent multiplexing up channel It is communicated；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is phone user in the base station range Quantity, N is D2D pairs of quantity in the base station range.

For example, in embodiments of the present invention, the quantity for setting phone user has M, and D2D is N number of to having, it should explanation It is that base station can distribute a up channel and down channel for each phone user, so being directed to each D2D pairs, can first counts Its sum of M throughput gain for being multiplexed up channel is calculated, then divided by M, obtains its average throughput for being multiplexed up channel and increase Benefit, then M throughput gain of down channel is multiplexed respectively compared with the average throughput gain, statistical value is flat greater than this The quantity of the throughput gain of the multiplexing down channel of equal throughput gain；

If quantity is more than predetermined integers, determine that the D2D communicates independent multiplexing down channel, otherwise this D2D pairs Individually multiplexing up channel is communicated, for example, M=20, N=15, in embodiments of the present invention, 2M/N=2.667, if on take Whole then value is 3, then value is 2 for lower rounding, it is assumed that rounding in use, then should statistical value greater than the average throughput flow gain Corresponding D2D determines the D2D under independent multiplexing the quantity of the throughput gain of multiplexing down channel if the quantity is more than 3 Row channel is communicated, and otherwise the D2D communicates independent multiplexing up channel.

The choosing value of above-mentioned 2M/N is mainly the upstream or downstream letter for the fairness and phone user for considering channel distribution It road at most can only be by a D2D to multiplexing.There are M up channel and M in the embodiment of the present invention, in base station range Down channel, then D2D is 2M to reusable total number of channels, it is contemplated that the fairness between D2D pairs makes each D2D pairs as far as possible The number of channel being assigned to is identical, then each D2D is 2M/N to the number of resources being assigned to.

S4: the throughput gain according to each D2D to the up channel and down channel that are multiplexed each phone user respectively And all D2D pairs of multiplexer mode, determine each D2D to corresponding multiple channels；

The throughput gain according to each D2D to the up channel and down channel that are multiplexed each phone user respectively And all D2D pairs of multiplexer mode, determine that each D2D to corresponding multiple channels, may include:

A1: based on all D2D pairs of multiplexer mode, in uplink throughput gain matrix Ω^uIn, it will individually be multiplexed downlink letter The D2D that road is communicated is to corresponding elementIt is updated to any value less than or equal to zero；In downlink throughput capacity gain matrix Ω^dIn, the D2D that independent multiplexing up channel is communicated is to corresponding elementIt is updated to any less than or equal to zero Value；

Any value less than or equal to zero mentioned in above-mentioned steps indicates, by D2D to corresponding elementOrMore Newly be less than or equal to zero value for any one, it should explanation, for different D2D pairs, this corresponding updated be less than or Null any value can be equal, can not also wait；Certainly, it in actual computerized algorithm, directly will can individually answer The D2D communicated with down channel is to corresponding elementIt is updated to zero, independent multiplexing up channel is communicated D2D to corresponding elementAlso it is updated to zero.

A2: initialization matrix X=[x_ji]_N×M=0, i.e., by the element whole zero setting in the matrix X of N row M column, make square Battle array X is null matrix, in addition defines two array Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N }.

A3: it finds out and meets be expert at j ∈ W in throughput gain matrix Ω_j, column i ∈ Ф_iGreatest member, determine square Battle array X in the corresponding element x of row and column where the corresponding greatest member_ji=1；

For example, it is assumed that in D2D to the throughput gain matrix Ω of multiplexing phone user's up channel^uIn, in the 5th row the The element of 7 columnFor current matrix Ω^uIn greatest member, be expert at j=5, column i=7 based on this greatest member, will The element x of be expert at j=5, column i=7 in matrix X₅₇It is assigned a value of 1；It is worth noting that, j=5 and i=7 described here 5 ∈ W should all be met_jAnd 7 ∈ Ф_i。

A4: Ф is updated_i=Ф_iI, W_j=W_jJ, i.e., updated array Ф_iTo remove element i from former array；It updates Array W afterwards_jTo remove element j from former array, for example, former array W_j={ 1,2,3,4,5,6 }, Ф_i=1,2,3,4,5, 6,7,8 }, in step A3, the element in the 5th row the 7th columnFor current matrix Ω^uIn greatest member, then from array W_j Middle removal element 5, from array Ф_iElement 7 is removed, therefore updated two numbers group is respectively as follows: W_j={ 1,2,3,4,6 }, Ф_i= {1,2,3,4,5,6,8}。

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jMember Plain T_jiGreater than 0, return step A3；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiIt is greater than 0, update W_j={ 1,2 ..., N }, return step A3；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, It then determines output matrix X, and determines each D2D to corresponding multiple channels according to matrix X.

Additionally it should illustrate, in the embodiment of the present invention, a channel is only capable of by a D2D to multiplexing, so, whenI.e. the M up channel or M down channel of phone user is all by D2D to multiplexing, the model of base station covering The channel that can be multiplexed in enclosing, should stop algorithm at this time, determine output matrix X；

And in throughput gain matrix Ω, all corresponding i ∈ Ф_iWith j ∈ W_jElement T_jiWhen being respectively less than equal to 0, this When no matter in the range of the covering of base station whether there is or not the channel that can be multiplexed, since D2D handles up to being multiplexed brought by remaining channel Flow gain be zero even negative value, the increase to the handling capacity of entire D2D system be without benefit or it is unfavorable, institute in this case It needs to stop algorithm, determines output matrix X.

In above-mentioned steps A1~A4, the quantity of phone user is M in base station range, i ∈ { 1,2 ..., M }, D2D Pair quantity be N number of, j ∈ { 1,2 ..., N }, and M > N, throughput gain matrix Ω=[T_ji]_N×M；

Matrix X characterizes up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when in matrix X characterization Row channel multiplexing situation matrix X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D and uses multiplexing honeycomb The throughput gain matrix Ω of family up channel^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel handle up Flow gainWhen matrix X characterizes down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to the throughput gain matrix Ω of multiplexing phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed honeycomb User C_iDown channel throughput gain

In the matrix X finally exported,Characterize phone user C_iUp channel distribute to D2D to D_j, Characterize phone user C_iUp channel be not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D is to D_j,Characterize phone user C_iDown channel be not yet assigned to D2D to D_j。

It is emphasized that step A1 to A4, even more can be effectively avoided base station preferable D2D pairs to throughput gain Distribute excessive channel, and the poor D2D of throughput gain is to can not be assigned to channel；Step A1 to A4 makes all D2D It is communicated to channel can be assigned to, therefore, the embodiment of the present invention can make more D2D dock into network.

S5: according to D2D pairs of multiplexer mode and each D2D to corresponding multiple channels, determine each D2D in its institute Transmission power on corresponding multiple channels, and for each D2D to distribute its corresponding to multiple channels and its corresponding to Corresponding transmission power on multiple channels.

It is described according to D2D pairs of multiplexer mode and each D2D to corresponding multiple channels, determine each D2D at it Transmission power on corresponding multiple channels, may is that

Based on D2D pairs of multiplexer mode, determine that the Signal to Interference plus Noise Ratio calculation formula of phone user and D2D pairs is respectively as follows:

Meanwhile determining the formula of the maximum of D2D system total rate and constraint condition are as follows:

The total rate calculations formula of the maximum of D2D system:

It s.t. is constraint condition are as follows:

This constraint condition indicates each phone user C_iMinimum channel quality requirement；

This constraint condition indicates D2D to D_jPower constraint；

In above-mentioned formulaExpression is directed to phone user C_iIt is allSum,Expression is directed to honeycomb use Family C_iIt is allSum, due toAndValue be 0 or 1, so constraint condition indicate each phone user uplink and Down channel at most can only be by a D2D to multiplexing；

In view of in the embodiment of the present invention, D2D pairs of multiplexer mode are as follows: individually multiplexing up channel carries out communication or individually multiple It is communicated with down channel, for D2D to D_jCannot concurrent multiplexing uplink and downlink channel, therefore, when D2D is to D_jIndividually in multiplexing When row channel is communicated, about the D2D to D_jIt is allIt is 0, i.e.,Likewise, when D2D is to D_jIt is individually multiple When being communicated with down channel, about the D2D to D_jIt is allIt is 0, i.e.,

IfOtherwiseIfOtherwise

This constraint condition indicates: D2D is to D_jIt is multiplexed phone user C_iUp channel when,ThenThere is value；D2D To D_jIt is multiplexed phone user C_iDown channel when,ThenThere is value.Obvious, D2D is to D_jIt is multiplexed phone user C_i's When up channel,D2D is to D_jIt is multiplexed phone user C_iDown channel when,

After each D2D determines the channel being multiplexed, D2D can to the calculation method of the transmission power in respective channel It is achieved by the prior art.

In order to optimize the calculation method of power, in embodiments of the present invention, firstly, based on each D2D to corresponding more A channel simplifies the respective formula of the total rate of the maximum of the D2D system and constraint condition are as follows:

Lagrangian is converted by the total function of rate of the maximum of D2D system and the respective formula of constraint condition:

Its dual function are as follows:

The operational formula of power optimized solution is acquired according to KKT condition are as follows:

Based on the operational formula of power optimized solution, each power optimized solution P is found out_ji ^*, based on each power optimized acquired Solve P_ji ^*, it is each D2D to the transmission power on the multiple channels distributed corresponding to it；Wherein, P_ji ^*It is D2D to D_jIn honeycomb User C_iTransmission power on upstream or downstream channel；

Wherein,

Phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N },For phone user C_iThe letter of up channel dry make an uproar Than,For phone user C_iDown channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iUplink letter The Signal to Interference plus Noise Ratio in road,It is D2D to D_jIt is multiplexed phone user C_iThe Signal to Interference plus Noise Ratio of down channel；

Indicate D2D to D_jTo phone user C_iUp channel multiplexed situation, if phone user C_iUp channel quilt D2D is to D_jOtherwise multiplexing then indicates D2D to D_jTo phone user C_iDown channel multiplexed situation, If phone user C_iDown channel by D2D to D_jMultiplexing, then otherwise

Channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α refers to for path loss Number, d are distance of the transmitting terminal to receiving end；ParameterIndicate D2D to D_jIn phone user C_iTransmission power on up channel,Indicate D2D to D_jIn phone user C_iTransmission power on down channel, P_BIndicate the transmission power of base station, P_CIndicate honeycomb User C_iTransmission power, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D to D_jBetween base station Channel gain, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel increases Benefit, N₀Indicate noise figure；

max R_D2DIndicate all D2D to the maximum rate of the system of composition, γ^uIndicate the uplink of preset phone user Minimum signal-to-noise ratio, γ^dIndicate the downlink minimum signal-to-noise ratio of preset phone user,Indicate D2D pairs of maximum transmission power；

P_ji ^*D2D is characterized to D_jIn phone user C_iPower optimized solution on up channelOr the power on down channel Optimal solutionWork as P_ji ^*Characterize the power optimized solution on up channelWhen, P_jiD2D is characterized to D_jIn C_iOn up channel Transmission powerγ characterizes the uplink minimum signal-to-noise ratio γ of preset phone user^u；Work as P_ji ^*Characterize the power on down channel Optimal solutionWhen, P_jiD2D is characterized to D_jIn C_iTransmission power on down channelγ is characterized under preset phone user Row minimum signal-to-noise ratio γ^d；WithRespectively D2D is to D_jIn phone user C_iTransmitting function on up channel and down channel Rate；

μ_jWithFor two multiplication factors in Lagrangian, μ_j ^*WithFor by being asked according to gradient search algorithm The optimal solution of the dual function obtained；maxAnd maxAll D2D are to composition respectively in up channel and down channel System the total rate of maximum.

In addition, above-mentioned gradient search algorithm comprising steps of

B1: initialization multiplication factorμ_j(0) and step parameter δ, in the embodiment of the present invention, assignment can be initializedμ_j(0)=0, δ=0.01；

Wherein, i ∈ { 1,2 ..., M }, j ∈ { 1,2 ..., N }, and corresponding element x_ji≠0；

B2: it updates

B3: t=t+1, return step B2 are updated；Until(μ_i(0),μ_i(1),…, μ_i(t+1)) it restrains, enablesμ_j ^*=μ_i(t+1)。

Wherein, t indicates the number of iterations of the gradient search algorithm,And μ_j(t) when respectively indicating the t times iteration Multiplication factorμ_j, element x_jiCharacterize elementOr elementWork as P_ji ^*Characterize the power optimized solution on up channel When, element x_jiCharacterize elementWork as P_ji ^*Characterize the power optimized solution on down channelWhen, element x_jiCharacterize element

It should be noted that working as x_jiWhen=0, corresponding P_ji ^*Value be necessarily zero, so being answered when using gradient search algorithm First determine whether x corresponding to corresponding i and j_jiValue whether be zero, only x_jiIt just needs to carry out gradient search algorithm when ≠ 0.

The embodiment of the present invention carries out emulation experiment in single base station cell that radius is 500M, wherein phone user and D2D are used Family is uniformly distributed, base station transmitting power 46dB, and the transmission power of phone user is 23dB, and D2D maximum transmission power is 23dB, The interference threshold value of the uplink and downlink channel of phone user is respectively 7dB and 13dB.

In Fig. 4~7, PCUOD, WU and WD characterize three kinds of channel allocation modes, and PCUOD characterizes institute of the embodiment of the present invention The Channel distribution mode of use, WU characterize D2D to the up channel of all multiplexing phone users, and WD characterizes D2D to whole multiplexings The down channel of phone user.CU in Fig. 3 and 4 indicates phone user.Fig. 5~7 and D2D described below handle up scale Show all D2D to the handling capacity of the system of composition.

When the quantity that the quantity of phone user is 60, D2D pairs is 20, the distance of two users of each D2D centering is At 20 meters, the case where above-mentioned three kinds of Channel distribution modes correspond to different throughput of system, is as shown in Figure 4, it is clear that, using this Channel distribution mode used by inventive embodiments, all D2D are maximum to the handling capacity of the system of composition, all D2D to and institute The handling capacity for the system for having phone user to constitute is maximum.

In Fig. 4, " CU up " characterizes the system that the up channel of all phone users is constituted, and " CU down " characterizes all bees The system that the down channel of nest user is constituted, " CU all " characterize the system that all phone users are constituted, and " D2D " characterization is all For D2D to the system of composition, " Total " characterizes all D2D to the system constituted with all phone users.

When D2D pairs of quantity is 20, when the distance of two users of each D2D centering is 20 meters, with phone user The variation of quantity, all D2D are as shown in Figure 5 to the variation of the handling capacity of the system of composition.As seen from the figure, either any letter The road method of salary distribution, all D2D can increase the handling capacity of the system of composition with the increase of phone user's quantity；Obviously , from figure 5 it can be seen that Channel distribution mode used by the embodiment of the present invention, is better than being obviously improved for D2D handling capacity Other Channel distribution modes.

It is corresponding using the Channel distribution mode of the embodiment of the present invention when the quantity of phone user is 100 such as Fig. 5 D2D throughput ratio uses D2D handling capacity when WU Channel distribution mode to be higher by 52%, when than using WD Channel distribution mode D2D handling capacity is higher by 225%.

When the quantity of phone user is 60, when the distance of two users of each D2D centering is 20 meters, with D2D pairs The variation of quantity, all D2D are as shown in Figure 6 to the variation of the handling capacity of the system of composition.As seen from the figure, either any letter The road method of salary distribution, all D2D can increase with D2D the increase of quantity the handling capacity of the system of composition, due to D2D pairs Reusable channel total amount is constant, therefore the increase of D2D handling capacity is all relatively slower.

From fig. 6 it can be seen that when the distance of two users of the quantity and each D2D centering of phone user is definite value When, the variation with D2D to quantity, using the embodiment of the present invention Channel distribution mode when corresponding D2D handling capacity it is always bright It is aobvious to be higher than the corresponding D2D handling capacity of other Channel distribution modes.When D2D to quantity be 30 when, using the letter of the embodiment of the present invention The road method of salary distribution, corresponding D2D throughput ratio use D2D handling capacity when WU Channel distribution mode to be higher by 45%, than using D2D handling capacity when WD Channel distribution mode is higher by 190%.

When the quantity that the quantity of phone user is 60, D2D pairs is 20, with each D2D centering two users away from From variation, all D2D are as shown in Figure 7 to the variation of the handling capacity of the system of composition.As seen from the figure, either any channel The method of salary distribution, all D2D handling capacities can all be reduced with the increase of the distance of two users of D2D centering.Abscissa in Fig. 7 " the distance between D2D " is the distance of two users of above-mentioned each D2D centering.

It can be seen from figure 7 that being implemented during the distance change of two users in D2D couples using the present invention Corresponding D2D handling capacity is apparently higher than the corresponding D2D of other Channel distribution modes always and handles up when the Channel distribution mode of example Amount.When the distance of two users of each D2D centering is 35 meters, using the Channel distribution mode of the embodiment of the present invention, Corresponding D2D throughput ratio uses D2D handling capacity when WU Channel distribution mode to be higher by 65%, than using WD Channel distribution mode When D2D handling capacity be higher by 186%.

D2D in a kind of Cellular Networks provided in an embodiment of the present invention is to resource allocation methods, using a D2D user to multiple With the mode of multiple phone user's resources, meanwhile, according to throughput gain caused by the specific space diversity of D2D user couple into The best multiplexer mode of selection is compared in row analysis, increases D2D significantly resource reusable in system, D2D handles up to system Amount greatly increases, and then the problem of avoid serious interference between D2D couples and phone user, and docks more D2D users Enter in network, improves D2D user to the overall performance of the system of composition.

Corresponding to above method embodiment, as shown in figure 3, the embodiment of the invention also provides the D2D in a kind of Cellular Networks To resource allocation device, which includes:

Channel gain determining module 210, for collecting all phone users and D2D pairs of channel shape in base station range State information, determines each channel gain, and each channel gain includes channel gain between each phone user and base station, each Channel gain, each D2D between D2D pairs and base station are to channel gain respectively between each phone user and each D2D is to internal channel gain；

Throughput gain computing module 220, for calculating separately each D2D to multiple respectively according to each channel gain With the up channel of each phone user and the throughput gain of down channel；

Multiplexer mode determining module 230 is multiplexed the uplink letter of each phone user based on it for being directed to each D2D pairs The throughput gain in road and be multiplexed each phone user down channel throughput gain comparison, determine the D2D to answering With mode are as follows: individually multiplexing up channel carries out communication or individually multiplexing down channel is communicated；

Channel determination module 240, for according to each D2D to the up channel and downlink for being multiplexed each phone user respectively The throughput gain of channel and all D2D pairs of multiplexer modes determine each D2D to corresponding multiple channels；

Resource distribution module 250, for the multiplexer mode and each D2D according to D2D pairs to corresponding multiple channels, really Fixed each D2D is each D2D to the multiple letters distributed corresponding to it to the transmission power on multiple channels corresponding to it Corresponding transmission power on road and multiple channels corresponding to it.

It uses specifically, calculating separately each D2D in the throughput gain computing module 220 each honeycomb is multiplexed respectively Formula used by the up channel at family and the throughput gain of down channel includes:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jMultiplexing Phone user C_iDown channel throughput gain,Indicate D2D pairs of maximum transmission power, P_BIndicate the transmitting function of base station Rate, P_CIndicate the transmission power of phone user, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D pairs D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jIt is interior The channel gain in portion, channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is path Loss index, d are distance of the transmitting terminal to receiving end, N₀Indicate noise figure.

Specifically, the multiplexer mode determining module 230 is specifically used for:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed up channel Average throughput flow gain, in turn,

Each throughput gain of down channel is multiplexed respectively compared with the average throughput gain, and statistical value Greater than the quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain, if quantity is more than predetermined integers, Determine that the D2D communicates independent multiplexing down channel, otherwise the D2D communicates independent multiplexing up channel；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is phone user in the base station range Quantity, N is D2D pairs of quantity in the base station range.

Specifically, the channel determination module 240 includes:

Throughput gain updating unit, for the multiplexer mode based on all D2D pairs, in uplink throughput gain matrix Ω^uIn, the D2D that independent multiplexing down channel is communicated is to corresponding elementIt is updated to any less than or equal to zero Value；In downlink throughput capacity gain matrix Ω^dIn, the D2D that independent multiplexing up channel is communicated is to corresponding element It is updated to any value less than or equal to zero；

Initial cell, for initializing matrix X=[x_ji]_N×M=0, define Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N}；

Element x_jiDetermination unit meets be expert at j ∈ W for finding out in throughput gain matrix Ω_j, column i ∈ Ф_i Greatest member, determine in matrix X with the corresponding element x of row and column where the corresponding greatest member_ji=1；

Channels determining unit, for updating Ф_i=Ф_iI, W_j=W_jJ,

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jMember Plain T_jiGreater than 0, element x is returned_jiDetermination unit；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiIt is greater than 0, update W_j={ 1,2 ..., N } returns to element x_jiDetermination unit；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, It then determines output matrix X, and determines each D2D to corresponding multiple channels according to matrix X；

Wherein, the phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station D2D in coverage area is to D_jQuantity be N number of, j ∈ { 1,2 ..., N }, M > N, throughput gain matrix Ω=[T_ji]_N×M；

Matrix X characterizes up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when in matrix X characterization Row channel multiplexing situation matrix X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D and uses multiplexing honeycomb The throughput gain matrix Ω of family up channel^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel handle up Flow gainWhen matrix X characterizes down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to the throughput gain matrix Ω of multiplexing phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed honeycomb User C_iDown channel throughput gain

Characterize phone user C_iUp channel distribute to D2D to D_j,Characterize phone user C_iUplink Channel is not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D to D_j,Characterize bee Nest user C_iDown channel be not yet assigned to D2D to D_j。

Specifically, the resource distribution module 250 includes transmission power determination unit and channel and power determining unit.

Transmission power determination unit determines phone user and D2D pairs of letter is dry makes an uproar for the multiplexer mode based on D2D pairs It is respectively as follows: than calculation formula

Meanwhile determining the formula of the maximum of D2D system total rate and constraint condition are as follows:

s.t.

IfOtherwiseIfOtherwise

Based on each D2D to corresponding multiple channels, by the phase of the maximum of the D2D system total rate and constraint condition Answer simplified formula are as follows:

Lagrangian is converted by the total function of rate of the maximum of D2D system and the respective formula of constraint condition:

Its dual function are as follows:

The operational formula of power optimized solution is acquired according to KKT condition are as follows:

Based on the operational formula of power optimized solution, each power optimized solution P is found out_ji ^*, based on each power optimized acquired Solve P_ji ^*, it is each D2D to the transmission power on the multiple channels distributed corresponding to it；Wherein, P_ji ^*It is D2D to D_jIn honeycomb User C_iTransmission power on upstream or downstream channel；

Wherein,

Phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N },For phone user C_iThe letter of up channel dry make an uproar Than,For phone user C_iDown channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iUplink letter The Signal to Interference plus Noise Ratio in road,It is D2D to D_jIt is multiplexed phone user C_iThe Signal to Interference plus Noise Ratio of down channel；

Indicate D2D to D_jTo phone user C_iUp channel multiplexed situation, if phone user C_iUp channel quilt D2D is to D_jOtherwise multiplexing then indicates D2D to D_jTo phone user C_iDown channel multiplexed situation, If phone user C_iDown channel by D2D to D_jMultiplexing, then otherwise

Channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α refers to for path loss Number, d are distance of the transmitting terminal to receiving end；ParameterIndicate D2D to D_jIn phone user C_iTransmission power on up channel,Indicate D2D to D_jIn phone user C_iTransmission power on down channel, P_BIndicate the transmission power of base station, P_CIndicate honeycomb User C_iTransmission power, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D to D_jBetween base station Channel gain, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel increases Benefit, N₀Indicate noise figure；

max R_D2DIndicate all D2D to the maximum rate of the system of composition, γ^uIndicate the uplink of preset phone user Minimum signal-to-noise ratio, γ^dIndicate the downlink minimum signal-to-noise ratio of preset phone user,Indicate D2D pairs of maximum transmission power；

P_ji ^*D2D is characterized to D_jIn phone user C_iPower optimized solution on up channelOr the power on down channel Optimal solutionWork as P_ji ^*Characterize the power optimized solution on up channelWhen, P_jiD2D is characterized to D_jIn C_iOn up channel Transmission powerγ characterizes the uplink minimum signal-to-noise ratio γ of preset phone user^u；Work as P_ji ^*Characterize the power on down channel Optimal solutionWhen, P_jiD2D is characterized to D_jIn C_iTransmission power on down channelγ is characterized under preset phone user Row minimum signal-to-noise ratio γ^d；WithRespectively D2D is to D_jIn phone user C_iTransmitting function on up channel and down channel Rate；

μ_jWithFor two multiplication factors in Lagrangian, μ_j ^*WithFor by being asked according to gradient search algorithm The optimal solution of the dual function obtained；maxAnd maxAll D2D are to composition respectively in up channel and down channel System the total rate of maximum.

Channel and power distributing unit, for for each D2D to distribute its corresponding to multiple channels and its corresponding to Multiple channels on corresponding transmission power.

Specifically, realizing that device used by the gradient search algorithm includes:

Initial submodule: for initializing multiplication factorμ_j(0) and step parameter δ,

Wherein, i ∈ { 1,2 ..., M }, j ∈ { 1,2 ..., N }, and corresponding element x_ji≠0；

Update submodule: for updating

It determines submodule: for updating t=t+1, updated t being returned and updates submodule；UntilWithIt restrains, enablesμ_j ^*= μ_i(t+1)；

Wherein, t indicates the number of iterations of the gradient search algorithm,And μ_j(t) when respectively indicating the t times iteration Multiplication factorμ_j, element x_jiCharacterize elementOr elementWork as P_ji ^*Characterize the power optimized solution on up channel When, element x_jiCharacterize elementWork as P_ji ^*Characterize the power optimized solution on down channelWhen, element x_jiCharacterize element

D2D in a kind of Cellular Networks provided in an embodiment of the present invention is to resource allocation device, using a D2D user to multiple With the mode of multiple phone user's resources, meanwhile, according to throughput gain caused by the specific space diversity of D2D user couple into The best multiplexer mode of selection is compared in row analysis, increases D2D significantly resource reusable in system, D2D handles up to system Amount greatly increases, and then the problem of avoid serious interference between D2D couples and phone user, and docks more D2D users Enter in network, improves D2D user to the overall performance of the system of composition.

It should be noted that, in this document, relational terms such as first and second and the like are used merely to a reality Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those Element, but also including other elements that are not explicitly listed, or further include for this process, method, article or equipment Intrinsic element.In the absence of more restrictions, the element limited by sentence " including one ... ", it is not excluded that There is also other identical elements in the process, method, article or apparatus that includes the element.

Each embodiment in this specification is all made of relevant mode and describes, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for device reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.

Those of ordinary skill in the art will appreciate that all or part of the steps in realization above method embodiment is can It is completed with instructing relevant hardware by program, the program can store in computer-readable storage medium, The storage medium designated herein obtained, such as: ROM/RAM, magnetic disk, CD.

The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the scope of the present invention.It is all Any modification, equivalent replacement, improvement and so within the spirit and principles in the present invention, are all contained in protection scope of the present invention It is interior.

Claims (10)

1. the D2D in a kind of Cellular Networks is to resource allocation methods, which is characterized in that the method includes the steps:

All phone users and D2D pairs of channel state information in base station range are collected, determines each channel gain, it is described each Channel gain includes channel gain between each phone user and base station, each D2D pairs of channel gain between base station, every A D2D to respectively between each phone user channel gain and each D2D to internal channel gain；

According to each channel gain, each D2D is calculated separately to the up channel and downlink for being multiplexed each phone user respectively The throughput gain of channel；

For each D2D pairs, throughput gain and each honeycomb of multiplexing based on its up channel for being multiplexed each phone user The comparison of the throughput gain of the down channel of user determines D2D pairs of the multiplexer mode are as follows: individually multiplexing up channel carries out Communication or independent multiplexing down channel are communicated；

To the throughput gain for the up channel and down channel for being multiplexed each phone user respectively and owned according to each D2D D2D pairs of multiplexer mode determines each D2D to corresponding multiple channels；

According to D2D pairs of multiplexer mode and each D2D to corresponding multiple channels, determine each D2D to corresponding to it Transmission power on multiple channels, and for each D2D to distribute its corresponding to multiple channels and multiple letters corresponding to it Corresponding transmission power on road.

2. the method according to claim 1, wherein described calculate separately each D2D to being multiplexed each bee respectively Formula used by the up channel of nest user and the throughput gain of down channel includes:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jHoneycomb is multiplexed to use Family C_iDown channel throughput gain,Indicate D2D pairs of maximum transmission power, P_BIndicate the transmission power of base station, P_C Indicate the transmission power of phone user, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D to D_jWith base Channel gain between standing, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal Channel gain, channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is path loss Index, d are distance of the transmitting terminal to receiving end, N₀Indicate noise figure.

3. being multiplexed up channel based on it the method according to claim 1, wherein described be directed to each D2D pairs Each throughput gain and multiplexing down channel each throughput gain comparison, determine D2D pairs of the multiplexer mode are as follows: Individually multiplexing up channel carries out communication or individually multiplexing down channel is communicated, comprising:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed being averaged for up channel Throughput gain, in turn,

Each throughput gain of down channel is multiplexed respectively compared with the average throughput gain, and statistical value is greater than The quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain determines if quantity is more than predetermined integers The D2D communicates independent multiplexing down channel, and otherwise the D2D communicates independent multiplexing up channel；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is the number of phone user in the base station range Amount, the quantity that N is D2D pairs in the base station range.

4. the method according to claim 1, wherein described use according to each D2D each honeycomb is multiplexed respectively The up channel at family and the throughput gain of down channel and all D2D pairs of multiplexer modes, determine each D2D to corresponding Multiple channels, comprising:

A1: based on all D2D pairs of multiplexer mode, in uplink throughput gain matrix Ω^uIn, will individually be multiplexed down channel into The D2D of row communication is to corresponding elementIt is updated to any value less than or equal to zero；In downlink throughput capacity gain matrix Ω^d In, the D2D that independent multiplexing up channel is communicated is to corresponding elementIt is updated to any value less than or equal to zero；

A2: initialization matrix X=[x_ji]_N×M=0, define Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N }；

A3: it finds out and meets be expert at j ∈ W in throughput gain matrix Ω_j, column i ∈ Ф_iGreatest member, determine matrix X In with the corresponding element x of row and column where the corresponding greatest member_ji=1；

A4: Ф is updated_i=Ф_iI, W_j=W_jJ,

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jElement T_ji Greater than 0, return step A3；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiGreater than 0, more New W_j={ 1,2 ..., N }, return step A3；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, it is determined that Output matrix X, and determine each D2D to corresponding multiple channels according to matrix X；

Wherein, the phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N }, M > N, throughput gain matrix Ω=[T_ji]_N×M, matrix X table Levy up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when matrix X characterizes up channel multiplexed situation square Battle array X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D and handles up multiplexing phone user's up channel Flow gain matrix Ω^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel throughput gainAs matrix X Characterize down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to multiplexing The throughput gain matrix Ω of phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iDown channel Throughput gain Characterize phone user C_iUp channel distribute to D2D to D_j,Honeycomb is characterized to use Family C_iUp channel be not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D to D_j,Characterize phone user C_iDown channel be not yet assigned to D2D to D_j。

5. the method according to claim 1, wherein it is described according to D2D pairs of multiplexer mode and each D2D to institute Corresponding multiple channels determine each D2D to the transmission power on multiple channels corresponding to it, comprising:

Based on D2D pairs of multiplexer mode, determine that the Signal to Interference plus Noise Ratio calculation formula of phone user and D2D pairs is respectively as follows:

Meanwhile determining the formula of the maximum of D2D system total rate and constraint condition are as follows:

s.t.

IfOtherwiseIfOtherwise

Based on each D2D to corresponding multiple channels, by the corresponding public affairs of the maximum of the D2D system total rate and constraint condition Formula simplifies are as follows:

Lagrangian is converted by the total function of rate of the maximum of D2D system and the respective formula of constraint condition:

Its dual function are as follows:

The operational formula of power optimized solution is acquired according to KKT condition are as follows:

Based on the operational formula of power optimized solution, each power optimized solution P is found out_ji ^*, based on each power optimized solution acquired P_ji ^*, it is each D2D to the transmission power on the multiple channels distributed corresponding to it；Wherein, P_ji ^*It is D2D to D_jIt is used in honeycomb Family C_iTransmission power on upstream or downstream channel；

Wherein,

Phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, in the base station range D2D to D_jQuantity be N number of, j ∈ { 1,2 ..., N },For phone user C_iUp channel Signal to Interference plus Noise Ratio,For phone user C_iDown channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iUp channel Signal to Interference plus Noise Ratio,It is D2D to D_jIt is multiplexed phone user C_iThe Signal to Interference plus Noise Ratio of down channel,Indicate D2D to D_jHoneycomb is used Family C_iUp channel multiplexed situation, if phone user C_iUp channel by D2D to D_jMultiplexing, thenOtherwise Indicate D2D to D_jTo phone user C_iDown channel multiplexed situation, if phone user C_iDown channel quilt D2D is to D_jMultiplexing, thenOtherwiseChannel gain g=G β d^-α, G therein is path loss constant, and β represents letter Road decline constant, α are path loss index, and d is distance of the transmitting terminal to receiving end；ParameterIndicate D2D to D_jIt is used in honeycomb Family C_iTransmission power on up channel,Indicate D2D to D_jIn phone user C_iTransmission power on down channel, P_BIt indicates The transmission power of base station, P_CIndicate phone user C_iTransmission power, g_BiIndicate phone user C_iChannel between base station increases Benefit, g_BjIndicate D2D to D_jChannel gain between base station, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal channel gain, N₀Indicate noise figure；max R_D2DIndicate all D2D to the maximum of the system of composition Rate, γ^uIndicate the uplink minimum signal-to-noise ratio of preset phone user, γ^dIndicate the downlink minimum noise of preset phone user Than,Indicate D2D pairs of maximum transmission power；

P_ji ^*D2D is characterized to D_jIn phone user C_iPower optimized solution on up channelOr the power optimized on down channel SolutionWork as P_ji ^*Characterize the power optimized solution on up channelWhen, P_jiD2D is characterized to D_jIn C_iTransmitting on up channel Powerγ characterizes the uplink minimum signal-to-noise ratio γ of preset phone user^u；Work as P_ji ^*Characterize the power optimized on down channel SolutionWhen, P_jiD2D is characterized to D_jIn C_iTransmission power on down channelγ characterizes the downlink of preset phone user most Small signal-to-noise ratio γ^d；WithRespectively D2D is to D_jIn phone user C_iTransmission power on up channel and down channel；μ_j WithFor two multiplication factors in Lagrangian, μ_j ^*WithTo pass through the antithesis letter acquired according to gradient search algorithm Several optimal solutions；maxAnd maxAll D2D are to the system of composition respectively in up channel and down channel Maximum total rate.

6. according to the method described in claim 5, it is characterized in that, the gradient search algorithm comprising steps of

B1: initialization multiplication factorμ_j(0) and step parameter δ,

Wherein, i ∈ { 1,2 ..., M }, j ∈ { 1,2 ..., N }, and corresponding element x_ji≠0；

B2: it updates

B3: t=t+1, return step B2 are updated；Until(μ_i(0),μ_i(1),…,μ_i(t+ 1) it) restrains, enablesμ_j ^*=μ_i(t+1)；

Wherein, t indicates the number of iterations of the gradient search algorithm,And μ_j(t) multiplication when the t times iteration is respectively indicated The factorμ_j, element x_jiCharacterize elementOr elementWork as P_ji ^*Characterize the power optimized solution on up channelWhen, member Plain x_jiCharacterize elementWork as P_ji ^*Characterize the power optimized solution on down channelWhen, element x_jiCharacterize element It indicates to functionSeek multiplication factorGradient,It indicates to functionSeek multiplication factor μ_jGradient.

7. the D2D in a kind of Cellular Networks is to resource allocation device, which is characterized in that described device includes:

Channel gain determining module, for collecting all phone users and D2D pairs of channel state information in base station range, Determine each channel gain, each channel gain include channel gain between each phone user and base station, it is each D2D pairs with Channel gain, each D2D between base station to respectively between each phone user channel gain and each D2D to inside Channel gain；

Throughput gain computing module, for according to each channel gain, calculating separately each D2D to being multiplexed each bee respectively The up channel of nest user and the throughput gain of down channel；

Multiplexer mode determining module, for being directed to each D2D pairs, based on handling up for its up channel for being multiplexed each phone user The comparison of the throughput gain of flow gain and the down channel of each phone user of multiplexing, determines D2D pairs of the multiplexer mode are as follows: Individually multiplexing up channel carries out communication or individually multiplexing down channel is communicated；

Channel determination module, for according to each D2D to the up channel and down channel for being multiplexed each phone user respectively Throughput gain and all D2D pairs of multiplexer modes determine each D2D to corresponding multiple channels；

Resource distribution module determines each for the multiplexer mode and each D2D according to D2D pairs to corresponding multiple channels D2D to the transmission power on multiple channels corresponding to it, and for each D2D to distribute its corresponding to multiple channels and Corresponding transmission power on multiple channels corresponding to it.

8. device according to claim 7, which is characterized in that calculated separately in the throughput gain computing module each Formula used by throughput gain of the D2D to the up channel and down channel that are multiplexed each phone user respectively includes:

Wherein,It is D2D to D_jIt is multiplexed phone user C_iUp channel throughput gain,It is D2D to D_jHoneycomb is multiplexed to use Family C_iDown channel throughput gain,Indicate D2D pairs of maximum transmission power, P_BIndicate the transmission power of base station, P_C Indicate the transmission power of phone user, g_BiIndicate phone user C_iChannel gain between base station, g_BjIndicate D2D to D_jWith base Channel gain between standing, g_ijIndicate D2D to D_jWith phone user C_iBetween channel gain, g_jjIndicate D2D to D_jInternal Channel gain, channel gain g=G β d^-α, G therein is path loss constant, and β represents channel fading constant, and α is path loss Index, d are distance of the transmitting terminal to receiving end, N₀Indicate noise figure.

9. device according to claim 7, which is characterized in that the multiplexer mode determining module is specifically used for:

For each D2D pairs, it is multiplexed each throughput gain of up channel according to it, determines that it is multiplexed being averaged for up channel Throughput gain, in turn,

Each throughput gain of down channel is multiplexed respectively compared with the average throughput gain, and statistical value is greater than The quantity of the throughput gain of the multiplexing down channel of the average throughput flow gain determines if quantity is more than predetermined integers The D2D communicates independent multiplexing down channel, and otherwise the D2D communicates independent multiplexing up channel；

Wherein, the predetermined integers are the upper rounding or lower rounding of 2M/N, and M is the number of phone user in the base station range Amount, the quantity that N is D2D pairs in the base station range.

10. device according to claim 7, which is characterized in that the channel determination module includes:

Throughput gain updating unit, for the multiplexer mode based on all D2D pairs, in uplink throughput gain matrix Ω^uIn, The D2D that independent multiplexing down channel is communicated is to corresponding elementIt is updated to any value less than or equal to zero；Under Row throughput gain matrix Ω^dIn, the D2D that independent multiplexing up channel is communicated is to corresponding elementIt is updated to small In null any value；

Initial cell, for initializing matrix X=[x_ji]_N×M=0, define Ф_i={ 1,2 ..., M }, W_j=1,2 ..., N }；

Element x_jiDetermination unit meets be expert at j ∈ W for finding out in throughput gain matrix Ω_j, column i ∈ Ф_iMost Big element, determine in matrix X with the corresponding element x of row and column where the corresponding greatest member_ji=1；

Channels determining unit, for updating Ф_i=Ф_iI, W_j=W_jJ,

IfAndAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iWith j ∈ W_jElement T_ji Greater than 0, element x is returned_jiDetermination unit；

IfAnd in throughput gain matrix Ω, there are corresponding i ∈ Ф_iElement T_jiGreater than 0, more New W_j={ 1,2 ..., N } returns to element x_jiDetermination unit；

IfOr in throughput gain matrix Ω, all corresponding i ∈ Ф_iElement T_jiRespectively less than it is equal to 0, it is determined that Output matrix X, and determine each D2D to corresponding multiple channels according to matrix X；

Wherein, the phone user C in the base station range_iQuantity be M, i ∈ { 1,2 ..., M }, the base station covering D2D in range is to D_jQuantity be N number of, j ∈ { 1,2 ..., N }, M > N, throughput gain matrix Ω=[T_ji]_N×M, matrix X table Levy up channel multiplexed situation matrix X^uOr down channel multiplexed situation matrix X^d, when matrix X characterizes up channel multiplexed situation square Battle array X^uWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D and handles up multiplexing phone user's up channel Flow gain matrix Ω^u, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iUp channel throughput gainAs matrix X Characterize down channel multiplexed situation matrix X^dWhen, element x_jiCharacterize elementThroughput gain matrix Ω characterizes D2D to multiplexing The throughput gain matrix Ω of phone user's down channel^d, element T_jiD2D is characterized to D_jIt is multiplexed phone user C_iDown channel Throughput gain Characterize phone user C_iUp channel distribute to D2D to D_j,Honeycomb is characterized to use Family C_iUp channel be not yet assigned to D2D to D_j,Characterize phone user C_iDown channel distribute to D2D to D_j,Characterize phone user C_iDown channel be not yet assigned to D2D to D_j。