CN105491510A - Service unloading method for resource sharing in dense heterogeneous cellular network - Google Patents

Abstract

The invention discloses a service unloading method for resource sharing in a dense heterogeneous cellular network. The method comprises the following steps: initializing a network resource state, forming a stable alliance partition for multiple channels, performing frequency spectrum leasing and service unloading in an alliance, and periodically repeating the previous steps by taking network environment change in to consideration at the same time. According to the invention, the service unloading method for resource sharing, provided by the invention, solves the problems of lack of frequency resources and unbalanced server loads existing in the dense heterogeneous cellular network; by taking such performance indexes as network performance, network time delays and the like into consideration, an effectiveness function is designed, a macro-micro base station cooperation framework based on an alliance gaming theory is brought forward, and high-efficiency server unloading can be realized for resource sharing; and the high-efficiency utilization of frequency spectrum resources in the alliance is realized based on the stable alliance partition, and the network overall performance is improved accordingly.

Description

The business discharging method that in a kind of intensive isomery cellular network, resourceoriented is shared

Technical field

The invention belongs to wireless communication technology field, particularly relate to the business discharging method that in a kind of intensive isomery cellular network, resourceoriented is shared.

Background technology

Along with mobile communication develop rapidly, the rare opportunities and challenges becoming future mobile communications and develop between the surge of business datum demand and frequency spectrum resource.But, under the environment of current heterogeneous network height fusion development, in macrocellular, dispose the various low power nodes (i.e. micro-base station) comprising microcellulor (Microcell), Home eNodeB (Femtocell) and via node (RelayNode) becomes one of key technology promoting network capacity.Due to constraints such as micro-base station range and antenna amount restrictions, cause existing in network the uneven and low inferior problem of the availability of frequency spectrum of serious load, hinder the further deployment of isomery cellular network and the lifting of the availability of frequency spectrum.

Business unloading, as a kind of scheme of the very promising lifting availability of frequency spectrum, is the technology realizing being unloaded by the partial data flow in network under given conditions other network.But in intensive isomery cellular network, isomerism and the complexity of wireless environment increase greatly, existing business Unloading Technology can not more efficiently Resources allocation and mitigate interference, therefore result in the problem of isomery cellular network degraded performance.Therefore, how in intensive isomery cellular network situation, to carry out the unloading of rational business in conjunction with other advanced technology and realize the efficient utilization of resources, to improve overall performance of network, become problem demanding prompt solution.Current in lifting, frequency spectrum resource utilizes, substantially comprise following several technology: the frequency spectrum leasing technology based on state information, the business Unloading Technology based on different Q os, based on game theoretic non-cooperation resource allocation methods.

Patent application document " a kind of hierarchical cellular network intermediate frequency spectrum rent method " (the publication number CN103561409A of Xian Electronics Science and Technology University, application number 201310385814.7, applying date 2013.08.29) in a kind of hierarchical cellular network intermediate frequency spectrum rent method is disclosed, the broadcast message that the method receives according to service provider adjusts frequency spectrum price and transmitting power in real time.But the weak point of the method is: time under intensive isomery honeycomb scene, to need between service provider interaction times and mutual information data comparatively large, and the requirement of real-time under the scene such as the wireless data service become when it can not ensure and dynamic wireless channel decline.

Patent application document " the business discharging method based on different QoS " (the publication number CN103327541A of Beijing University of Post & Telecommunication, application number 201310187980.6, applying date 2013.05.20) in a kind of business discharging method based on different QoS is disclosed, the method devises the business unloading relative influence factor and selects business unloading target BS accordingly, improves the utilance of Radio Resource.The method Shortcomings part is: the problem causing base station performance to decline after not considering target BS unloading business, and does not take suitable incentive measure to promote to cooperate between base station.

Patent application document " based on game theoretic radio resource optimizing method in LTE-A relay system " (the publication number CN103369568A of Xi'an Communications University, application number 201310291899.2, applying date 2013.07.11) in disclose based on game theoretic radio resource optimizing method in a kind of LTE-A relay system, the method mainly solves the compromise problem when introducing relaying in cellular cell between throughput of system and user fairness.The deficiency that the method exists is: need repeatedly to solve first-order partial derivative in its game resource optimization process, calculation of complex and convergence rate is slow.

Summary of the invention

The present invention is directed to above-mentioned the deficiencies in the prior art, provide the business discharging method that in a kind of intensive isomery cellular network, resourceoriented is shared, the method is by grand-micro-base station cooperation framework in research intensive isomery cellular network, and in conjunction with frequency spectrum leasing and business Unloading Technology, achieve the efficiency utilization of frequency spectrum resource, and improve overall performance of network.

For solving problems of the prior art, the concrete technical scheme of employing is:

The business discharging method that in intensive isomery cellular network, resourceoriented is shared, comprises the following steps:

S1, initialization network resource status;

S2, form Stable coalitions subregion towards multichannel;

Frequency spectrum leasing and business unloading in S3, alliance;

When S4, consideration change of network environment, periodically repeat above step, realize frequency spectrum leasing and business unloading in the alliance under multichannel Stable coalitions subregion.

Preferred scheme, in initialization network resource status described in step S1, the partial frequency spectrum resource of the multiplexing macro base station in micro-base station, and macro base station and micro-base station keep fixing descending power signal transmission.Described fixing descending power is maximum descending power.

Preferred scheme further, described in step S2 towards the method step of multichannel formation Stable coalitions subregion is:

S21, towards multichannel build interference list;

S22, formation Stable coalitions.

The described method towards multichannel structure interference list is:

S21a, based on receive channel strength instruction (RSSI), macro base station detects the micro-base station be active in all downlink sub-channels k ∈ Κ (Κ is multiplexing number of sub-channels) be re-used by affiliated macro base station user, micro-base station also perception macro base station user in multiplexing subchannel;

Macro base station user m ∈ k in S21b, all multiplexing downlink sub-channels _m(k _mmacro base station number of users in multiplexing subchannel k) the micro-base station detected is disturbed by by force to weak sequence, form the interference list of respective sub-channel respectively, and feed back to macro base station by wireless channel;

The macro base station user perceived in multiplexing subchannel sorts by distance by S21c, micro-base station n ∈ N (N is micro-base station number) from the near to the remote.

The method of described formation Stable coalitions is:

S22a, macro base station based on interference list, successively with the micro-base station n ∈ N in respective sub-channel k _k, N _k∈ N (N _kmicro-base station number for multiplexing subchannel k) consult in pairs, form new alliance, then as potential cooperation object when micro-base station n has a mind to cooperation;

The micro-base station n of S22b, macro base station and potential cooperation calculates the user utility x before alliance _i(S _n, Π _n), (i ∈ S _n, S _npotential alliance for this reason under channel, i is macro base station user m ∪ micro-base station user l ∈ L in potential alliance _n, L _nfor micro-base station user quantity that micro-base station is served, Π _nfor alliance's subregion) and alliance total utility ν (S _n, Π _n).Alliance total utility ν (S before alliance _n, Π _Ν) calculating formula is as follows:

Above formula represents the effectiveness sum of all users in potential cooperative alliances, | S _n| be number of users in alliance, represent the macro base station user utility in current potential cooperative alliances and the summation of micro-base station user effectiveness.Wherein, the user utility x before alliance _i(S _n, Π _n) macro base station user utility x can be divided into _m(S _n, Π _n) and micro-base station user effectiveness x _l(S _n, Π _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x}_m(S_n,{\Π}_N)=\frac{{\mathrm{\μ}}_m^{NC}{(S_n,{\Π}_N)}^{\mathrm{\δ}}}{{\left(D_m^{NC}\right)}^{(1-\mathrm{\δ})}}\\ x_l(S_n,{\Π}_N)=\frac{{\mathrm{\μ}}_l^{NC}{(S_n,{\Π}_N)}^{\mathrm{\δ}}}{{\left(D_l^{NC}\right)}^{(1-\mathrm{\δ})}}\end{array}\right.$

In formula, (S _n, Π _n) ^δwith (S _n, Π _n) ^δbe respectively reached at the information rate of macro base station user and reached at the information rate of micro-base station user of non-cooperation under current alliance subregion, S _npotential alliance for this reason under channel, Π _nfor current alliance subregion, with be respectively the average delay of macro base station user and the average delay of micro-base station user under non-cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, the macro base station user Ke Da information rate under non-cooperation (S _n, Π _n) ^δbe calculated as follows:

${\mathrm{\μ}}_m^{NC}(S_n,{\Π}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{{\mathrm{\Σ}}_{n\∈N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\σ}}^2})$

Reached at information rate when this formula represents that macro base station user m is subject to micro-base station interference in same sub-channel, wherein, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the micro-base station interference signal sum in the same sub-channel that macro base station user receives, represent macro base station user m Signal to Interference plus Noise Ratio now.

Micro-base station user of non-cooperation can reach information rate (l ∈ l _n, l _nmicro-base station user quantity for micro-base station is served) be calculated as follows:

${\mathrm{\μ}}_l^C(S_n,{\Π}_N)=B_{l}log(1+\frac{|H_{n,l}{|}^2{P}_n}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\≠n}{\overset{N_k}{\Σ}}|H_{i,l}{|}^2{P}_i+{\mathrm{\σ}}^2})$

The reached information rate of this formula when micro-base station user l is subject to the interference of micro-base station in macro base station and same sub-channel before representing alliance, wherein, B _lrepresent the bandwidth of subchannel belonging to this micro-base station user l, | H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _{0, l}| ²represent the channel gain between macro base station and micro-base station user l, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, | H _n,l| ²p _nrepresent the available signal power that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that micro-base station user l receives the interference signal of the same sub-channel of macro base station, represent that micro-base station user l receives the interference signal of the micro-base station i under same sub-channel, represent micro-base station user l Signal to Interference plus Noise Ratio now;

with be respectively the average delay of macro base station user and the average delay of micro-base station user under non-cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^{NC}=\frac{{{\widetilde{\mathrm{\λ}}}_m}^{NC}}{2{\mathrm{\μ}}_m^{NC}({\mathrm{\μ}}_m^{NC}-{{\widetilde{\mathrm{\λ}}}_m}^{NC})}\\ D_l^{NC}=\frac{{{\widetilde{\mathrm{\λ}}}_l}^{NC}}{2{\mathrm{\μ}}_l^{NC}({\mathrm{\μ}}_l^{NC}-{{\widetilde{\mathrm{\λ}}}_l}^{NC})}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering that maximum retransmit number is D under non-cooperation with the practical communication load of micro-base station user be respectively

$\left\{\begin{array}{c}{{\widetilde{\mathrm{\λ}}}_m}^{NC}={\mathrm{\λ}}_m\underset{d=1}{\overset{D}{\Σ}}{\mathrm{Pt}}_m^{NC}{(1-{\mathrm{Pt}}_m^{NC})}^{d-1}\\ {{\widetilde{\mathrm{\λ}}}_m}^{NC}={\mathrm{\λ}}_l\underset{d=1}{\overset{D}{\Σ}}{\mathrm{Pt}}_l^{NC}{(1-{\mathrm{Pt}}_l^{NC})}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of micro-base station-micro-base station user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, with represent the probability of macro base station user transmission success and the probability of micro-base station user transmission success under non-cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, (i=l or m) represent that data the d time retransmit just successfully probability, calculating formula is distinguished as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^{NC}=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{{\Σ}_{n\∈N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\σ}}^2}\≥{\mathrm{\γ}}_m\}\\ {\mathrm{Pt}}_l^{NC}=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_l}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\≠n}{\overset{N_k}{\Σ}}|H_{i,l}{|}^2{P}_i+{\mathrm{\σ}}^2}\≥{\mathrm{\γ}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent that macro base station user receives the micro-base station interference signal sum in same sub-channel, represent macro base station user m Signal to Interference plus Noise Ratio now;

| H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, | H _{0, l}| ²represent the channel gain between macro base station and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, γ _lfor micro-base station user Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that micro-base station user receives the interference signal of the same sub-channel of macro base station, represent that micro-base station user receives the interference signal of the micro-base station under same sub-channel, represent micro-base station user l Signal to Interference plus Noise Ratio now.

The user utility x after alliance is estimated in S22c, the micro-base station of macro base station and potential alliance _i ^*(S _n, Π _Ν) and alliance total utility ν ^*(S _n, Π _Ν), the alliance total utility ν after alliance ^*(S _n, Π _Ν) calculating formula is as follows:

This formula represents alliance S _nin the effectiveness sum of all users, | S _n| be the number of users of current alliance, represent the macro base station user utility after to cooperation in current alliance subregion and the summation of micro-base station user effectiveness.Wherein, the user utility x after alliance _i ^*(S _n, Π _Ν) macro base station user utility x can be divided into _m ^*(S _n, Π _n) and micro-base station user effectiveness x _l ^*(S _n, Π _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x_m}^{*}(S_n,{\Π}_N)=\frac{{\mathrm{\μ}}_m^C{(S_n,{\Π}_N)}^{\mathrm{\δ}}}{{\left(D_m^C\right)}^{(1-\mathrm{\δ})}}\\ {x_l}^{*}(S_n,{\Π}_N)=\frac{{\mathrm{\μ}}_l^C{(S_n,{\Π}_N)}^{\mathrm{\δ}}}{{\left(D_l^C\right)}^{(1-\mathrm{\δ})}}\end{array}\right.$

In formula, (S _n, Π _n) and (S _n, Π _n) be respectively reached at the information rate of macro base station user and reached at the information rate of micro-base station user of cooperation under current alliance subregion, with be respectively the average delay of macro base station user and the average delay of micro-base station user under cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, reached at the information rate of macro base station user under cooperation (S _n, Π _n) be calculated as follows:

${\mathrm{\μ}}_m^C(S_n,{\Π}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\∈N_k\backslash S_n}{\Σ}|H_{n,m}{|}^2{P}_n+{\mathrm{\σ}}^2})$

This formula tabular form macro base station user m is subject to micro-base station n, the n ∈ N outside alliance after alliance _ks _ninterference time reached at information rate, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives, namely macro base station user receives interference signal and the white Gaussian noise sum of the micro-base station outside alliance, represent the Signal to Interference plus Noise Ratio under macro base station user cooperation;

Reached at the information rate of micro-base station user l under cooperation (S _n, Π _n) be calculated as follows:

${\mathrm{\μ}}_l^C(S_n,{\Π}_N)=B_{l}log(1+\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\≠n}{\overset{N_k\backslash S_n}{\Σ}}|H_{i,l}{|}^2{P}_i+{\mathrm{\σ}}^2})$

This formula represents that micro-base station user is subject to micro-base station n ∈ N outside alliance after alliance _ks _ninterference time reached at information rate, wherein, B _lrepresent the bandwidth of subchannel belonging to this micro-base station user l, | H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, | H _n,l| ²p _nrepresent the available signal power that micro-base station user receives, namely micro-base station user receives interference signal and the white Gaussian noise sum of the micro-base station outside alliance, for the Signal to Interference plus Noise Ratio of micro-base station user under cooperation;

with be respectively the average delay of macro base station user and the average delay of micro-base station user under cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^C=\frac{{\widetilde{\mathrm{\λ}}}_m^C}{2{\mathrm{\μ}}_m^C({\mathrm{\μ}}_m^C-{\widetilde{\mathrm{\λ}}}_m^C)}\\ D_l^C=\frac{{\widetilde{\mathrm{\λ}}}_l^C}{2{\mathrm{\μ}}_l^C({\mathrm{\μ}}_l^C-{\widetilde{\mathrm{\λ}}}_l^C)}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering under cooperation that maximum retransmit number is D with the practical communication load of micro-base station user be respectively

$\left\{\begin{array}{c}{\widetilde{\mathrm{\λ}}}_m^C={\mathrm{\λ}}_m^C\underset{d=1}{\overset{D}{\Σ}}{\mathrm{Pt}}_m^C{(1-{\mathrm{Pt}}_m^C)}^{d-1}\\ {\widetilde{\mathrm{\λ}}}_l^C={\mathrm{\λ}}_l^C\underset{d=1}{\overset{D}{\Σ}}{\mathrm{Pt}}_l^C{(1-{\mathrm{Pt}}_l^C)}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of micro-base station-micro-base station user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, with the probability of macro base station user transmission success and the probability of micro-base station user transmission success under expression cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, calculating formula is respectively as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^C=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\∈N_k\backslash S_n}{\Σ}|H_{n,m}{|}^2{P}_n+{\mathrm{\σ}}^2}\≥{\mathrm{\γ}}_m\}\\ {\mathrm{Pt}}_l^C=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\≠n}{\overset{N_k\backslash S_n}{\Σ}}|H_{i,l}{|}^2{P}_i+{\mathrm{\σ}}^2}\≥{\mathrm{\γ}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives, represent that macro base station user receives whole interference signals of micro-base station in same sub-channel outside alliance, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under macro base station user m cooperation;

| H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, γ _lfor micro-base station user Successful transmissions desired value, | H _n,l| ²p represents the available signal power that micro-base station user receives, represent the interference signal of micro-base station in same sub-channel outside the alliance that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under micro-base station user l cooperation.

S22d, as the user utility x after macro base station and micro-base station cooperation _i ^*(S _n, Π _n) and alliance total utility ν ^*(S _n, Π _n) be all not less than cooperation before effectiveness time, i.e. x _i ^*(S _n, Π _n) >x _i(S _n, Π _n) and ν ^*(S _n, Π _n) > ν (S _n, Π _n), then fed back by wire message way between macro base station and micro-base station and formally form alliance; Otherwise, judge whether macro base station and next micro-base station of potential cooperation form new alliance;

S22e, macro base station sequentially perform S22a, S22b, S22c, S22d with the micro-base station in all interference lists according to the order of sequence, until Stable coalitions is formed in all multiplexing subchannels.Alliance's total utility that in arbitrary subchannel, micro-base station departs from the new alliance subregion of current coalition formation is less than alliance's total utility of current steady alliance subregion.

Further preferred scheme, in alliance described in step S3, the method for frequency spectrum leasing and business unloading is:

S31, based on the Stable coalitions subregion formed, according to frequency spectrum leasing model, solved the optimal solution of optimum frequency spectrum leasing coefficient by convex optimization tool;

S32, macro base station lease to corresponding micro-BTS channel frequency spectrum resource according to the optimum frequency spectrum leasing coefficient calculating gained, the macro base station user simultaneously notifying to take this channel disconnects former link and sets up new downlink transfer with corresponding micro-base station and links, and the frequency spectrum resource of the relaying window coefficient ratio in gained frequency spectrum resource is distributed to macro base station user and is used for business transmission by corresponding micro-base station.

Wherein, the method solving the optimal solution of optimum frequency spectrum leasing coefficient described in step S31 is:

Sub-channel spectra multiplexing for micro-base station is normalized to a unit length by frequency spectrum leasing model, each unit is divided into three parts simultaneously:

Part I unit length is 1-α, and macro base station signal transmission gives affiliated macro base station user;

Part II unit length is α β, and unloaded macro base station user is given as relay transmission signal in micro-base station;

Part III unit length is α (1-β), and micro-base-station transmission signal gives affiliated micro-base station user;

Wherein, α is that frequency spectrum leasing coefficient represents that macro base station leases to the channel spectrum resource of micro-base station, and β is that relay transmission window coefficient represents that micro-base station is the relaying frequency spectrum resource of the macro base station user distribution of unloading;

According to frequency spectrum leasing model, the utility function under this model is proposed, and by the convex optimization tool optimal solution that to solve with frequency spectrum leasing factor alpha and relay transmission window coefficient β be variable:

$\left\{\begin{array}{cc}\underset{\mathrm{\&alpha;},\mathrm{\&beta;}}{\max }& x_i^{\&prime;}(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})\\ \underset{\mathrm{\&alpha;},\mathrm{\&beta;}}{\max }& v^{\&prime;}(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})\\ s.t.& 0<\mathrm{\&alpha;},\mathrm{\&beta;}<1\end{array}\right.$

Above formula represents and is meeting 0< α, during the condition of β <1, solves the maximum user utility x based on frequency spectrum leasing model in Stable coalitions subregion _i' (S' _n, Π ' _n) and maximum alliance total utility ν ' (S _n', Π ' _n), S' _nfor the alliance under current steady alliance subregion, Π ' _nfor current steady alliance subregion.Wherein, user utility x _i＇ (S _n', Π ' _n) macro base station user utility x can be divided into _m＇ (S' _n, Π ' _n) and micro-base station user effectiveness x ＇ _l(S' _n, Π ' _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x}_m^{\&prime;}(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})=\frac{{\mathrm{\&alpha;\&beta;\&mu;}}_m^{\&prime;}{(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\&Pi;}_N^{\&prime;})}^{\mathrm{\&delta;}}}{D_m^{\&prime;(1-\mathrm{\&delta;})}})\\ x_l^{\&prime;}(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})=\frac{\mathrm{\&alpha;}(1-\mathrm{\&beta;}){\mathrm{\&mu;}}_l^{\&prime;}{(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\&Pi;}_N^{\&prime;})}^{\mathrm{\&delta;}}}{D_l^{\&prime;(1-\mathrm{\&delta;})}})\end{array}\right.$

With

δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, α and β is respectively frequency spectrum leasing coefficient and relay transmission window coefficient.Wherein, reached at the information rate μ ' of the macro base station user in optimization problem _m(α, β, S' _n, Π ' _n) and reached at the information rate μ of micro-base station user _l' (α, β, S' _n, Π ' _n) calculating formula is as follows:

$\left\{\begin{array}{c}{\mathrm{\&mu;}}_m^{\&prime;}(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\&Pi;}_N^{\&prime;})={\mathrm{\&alpha;\&beta;\&mu;}}_m^C(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})\\ {\mathrm{\&mu;}}_l^{\&prime;}(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\&Pi;}_N^{\&prime;})=\mathrm{\&alpha;}(1-\mathrm{\&beta;}){\mathrm{\&mu;}}_l^C(S_n^{\&prime;},{\&Pi;}_N^{\&prime;})\end{array}\right.$

Represent macro base station user and the reached information rate of micro-base station user under frequency spectrum leasing model under current alliance subregion, (S' _n, Π ' _n) and (S' _n, Π ' _n) be that macro base station user Ke Da information rate under current alliance subregion under cooperation and micro-base station user can reach information rate respectively, calculating formula is as shown in S22c.

Further preferred scheme again, processing method when considering change of network environment described in step S4 is: when change of network environment, periodically repeat step S2 towards frequency spectrum leasing and business unloading in the alliance that multichannel forms Stable coalitions subregion and step S3; When network environment does not change, periodically repeat the initialization network resource status of step S1, S2 towards frequency spectrum leasing and business unloading in the alliance that multichannel forms Stable coalitions subregion and step S3.

By adopting above technical scheme, the business discharging method that in a kind of intensive isomery cellular network of the present invention, resourceoriented is shared is compared with the prior art, and its technique effect is:

The first, the present invention is directed in prior art and be not suitable for processing the resource allocation problem in dense network, provide the business discharging method that in a kind of intensive isomery cellular network, resourceoriented is shared, it efficiently solves the unbalanced problem of business, and then improves the availability of frequency spectrum.

The second, the present invention is directed in prior art and lack suitable incentive measure and to impel between base station the problems such as enforcement business unloading, by frequency spectrum leasing compensating technique partial frequency spectrum resource, win-win situation between base station is achieved to the base station of unloading business.

Three, the present invention is directed to calculation of complex in prior art and the slow problem of convergence rate, theoretical based on Game with Coalitions, define the Stable coalitions under multichannel, in alliance, solve optimum utility problem, while decreasing amount of calculation, accelerate convergence rate.

Accompanying drawing explanation

Fig. 1 is the schematic diagram of collaborative program and non-collaborative program difference under isomery honeycomb scene provided by the invention;

Fig. 2 is the business discharging method that in a kind of intensive isomery cellular network provided by the invention, resourceoriented is shared;

Fig. 3 is formation multichannel Stable coalitions subregion flow chart provided by the invention;

Fig. 4 is frequency spectrum leasing model schematic provided by the invention;

Fig. 5 the invention provides embodiment with Home eNodeB gain effect during Home eNodeB number change;

Fig. 6 the invention provides embodiment macro base station user's gain effect and deployed position relation.

Embodiment

Clearly understand to make object of the present invention, technical scheme and advantage, in conjunction with the accompanying drawings and embodiments the present invention is described in further detail, should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.

Below in conjunction with Fig. 1 to Fig. 6 and specific embodiment, application principle of the present invention is further described.

The embodiment of the present invention is mainly described based on home base station network (FemtocellNetworks), aims to provide the business discharging method that in a kind of intensive isomery cellular network, resourceoriented is shared.

Fig. 1 is the schematic diagram of collaborative program and non-collaborative program difference under isomery honeycomb scene provided by the invention, it should be noted that the present invention is not limited to home base station network in current schematic diagram and Home eNodeB (FBS) quantity and number of users.

The business discharging method that in intensive isomery cellular network, resourceoriented is shared, comprises the steps:

Step 1, initialization network resource status

Initialization network state is Home eNodeB multiplexing macro base station partial frequency spectrum resource (available subchannels), and macro base station and Home eNodeB keep stablizing descending power signal transmission.

Step 2, form Stable coalitions subregion towards multichannel

In units of each independently orthogonal sub-channels, macro base station and Home eNodeB are respectively based on channel strength instruction (RSSI) received, formed and disturb list accordingly, macro base station consults (whether having cooperative desire) in pairs to find potential alliance object according to the order of sequence from interference list top with Home eNodeB, and periodically to judge in current investigation subchannel whether cooperation forms new alliance, until towards multi channel Stable coalitions subregion formation for macro base station and potential Home eNodeB.

Frequency spectrum leasing and business unloading in step 3, alliance

Based on the Stable coalitions subregion formed, propose the utility function under this model according to frequency spectrum leasing model, and solve utility function optimal solution by convex optimization tool.Macro base station leases to a certain proportion of channel spectrum resource of corresponding Home eNodeB according to the optimum frequency spectrum leasing coefficient calculating gained, the macro base station user simultaneously notifying to take this channel disconnects former link and sets up new downlink transfer with corresponding Home eNodeB and links, and the frequency spectrum resource of the relaying window coefficient ratio in gained frequency spectrum resource is distributed to macro base station user and is used for business transmission by corresponding Home eNodeB.

When step 4, consideration change of network environment, as user mobility and new base station deployment etc., periodically repeat above step, realize frequency spectrum leasing and business unloading in the alliance under multichannel Stable coalitions subregion.

On the basis of technique scheme, as shown in Figure 2, its concrete steps comprise described step 2:

(1), interference list is built towards multichannel

1a), based on the channel strength instruction (ReceivedSignalStrengthIndicator received; RSSI), macro base station detects the Home eNodeB be active in all downlink sub-channels k ∈ Κ (Κ is multiplexing number of sub-channels) be re-used by affiliated macro base station user, Home eNodeB also perception macro base station user in multiplexing subchannel.

1b), the macro base station user m ∈ k in all multiplexing downlink sub-channels _m(k _mmacro base station number of users in multiplexing subchannel k) Home eNodeB detected is disturbed by the interference list forming by force respective sub-channel to weak sequence respectively, and feed back to macro base station by wireless channel;

1c), the macro base station user perceived in multiplexing subchannel sorts by distance by Home eNodeB n ∈ N (N is Home eNodeB quantity) from the near to the remote;

(2), Stable coalitions is formed

2a), macro base station based on interference list, successively with the Home eNodeB n ∈ N in respective sub-channel k _k, N _k∈ N (N _khome eNodeB quantity for multiplexing subchannel k) consult in pairs, form new alliance, then as potential cooperation object when Home eNodeB n has a mind to cooperation;

2b), macro base station and potential cooperation Home eNodeB n calculate the user utility x before alliance _i(S _n, Π _n), (i ∈ S _n, S _npotential alliance for this reason under channel, i is macro base station user m ∪ femtocell user l ∈ L in potential alliance _n, L _nfor the femtocell user quantity of Home eNodeB service, Π _nfor alliance's subregion) and alliance total utility ν (S _n, Π _n).Alliance total utility ν (S before alliance _n, Π _Ν) calculating formula is as follows:

Above formula represents the effectiveness sum of all users in potential cooperative alliances, | S _n| be number of users in alliance, represent the macro base station user utility in current potential cooperative alliances and the summation of femtocell user effectiveness.Wherein, the user utility x before alliance _i(S _n, Π _n) macro base station user utility x can be divided into _m(S _n, Π _n) and femtocell user effectiveness x _l(S _n, Π _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x}_m(S_n,{\&Pi;}_N)=\frac{{\mathrm{\&mu;}}_m^{NC}{(S_n,{\&Pi;}_N)}^{\mathrm{\&delta;}}}{{\left(D_m^{NC}\right)}^{(1-\mathrm{\&delta;})}}\\ x_l(S_n,{\&Pi;}_N)=\frac{{\mathrm{\&mu;}}_l^{NC}{(S_n,{\&Pi;}_N)}^{\mathrm{\&delta;}}}{{\left(D_l^{NC}\right)}^{(1-\mathrm{\&delta;})}}\end{array}\right.$

In formula, (S _n, Π _n) ^δwith (S _n, Π _n) ^δbe respectively reached at the information rate of macro base station user and reached at the information rate of femtocell user of non-cooperation under current alliance subregion, S _npotential alliance for this reason under channel, Π _nfor current alliance subregion, with be respectively the average delay of macro base station user and the average delay of femtocell user under non-cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, the macro base station user Ke Da information rate under non-cooperation (S _n, Π _n) ^δbe calculated as follows:

${\mathrm{\&mu;}}_m^{NC}(S_n,{\&Pi;}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{{\mathrm{\&Sigma;}}_{n\&Element;N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2})$

Reached at information rate when this formula represents that macro base station user m is subject to the interference of Home eNodeB in same sub-channel, wherein, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between Home eNodeB and macro base station user m, P _nrepresent the transmitting power of Home eNodeB n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the Home eNodeB interference signal sum in the same sub-channel that macro base station user receives, represent macro base station user m Signal to Interference plus Noise Ratio now.

The femtocell user of non-cooperation can reach information rate (l ∈ l _n, l _nfemtocell user quantity for Home eNodeB service) be calculated as follows:

${\mathrm{\&mu;}}_l^{NC}(S_n,{\&Pi;}_N)=B_l\log (1+\frac{|H_{n,l}{|}^2{P}_n}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\&NotEqual;n}{\overset{N_k}{\&Sigma;}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2})$

Reached at information rate when femtocell user l is subject to the interference of Home eNodeB in macro base station and same sub-channel before this formula represents alliance, wherein, B _lrepresent the bandwidth of subchannel belonging to this femtocell user l, | H _n,l| ²represent the channel gain between Home eNodeB n and femtocell user l, P _nrepresent the transmitting power of Home eNodeB n, | H _{0, l}| ²represent the channel gain between macro base station and femtocell user l, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between Home eNodeB i and femtocell user l, P _ifor the transmitting power of Home eNodeB, | H _n,l| ²p _nrepresent the available signal power that femtocell user receives, represent the interference signal that femtocell user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that femtocell user l receives the interference signal of the same sub-channel of macro base station, p _iexpression femtocell user l receives the interference signal of the Home eNodeB i under same sub-channel, represent femtocell user l Signal to Interference plus Noise Ratio now;

with be respectively the average delay of macro base station user and the average delay of femtocell user under non-cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^{NC}=\frac{{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC}}{2{\mathrm{\&mu;}}_m^{NC}({\mathrm{\&mu;}}_m^{NC}-{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC})}\\ D_l^{NC}=\frac{{{\widetilde{\mathrm{\&lambda;}}}_l}^{NC}}{2{\mathrm{\&mu;}}_l^{NC}({\mathrm{\&mu;}}_l^{NC}-{{\widetilde{\mathrm{\&lambda;}}}_l}^{NC})}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering that maximum retransmit number is D under non-cooperation with the practical communication load of femtocell user be respectively

$\left\{\begin{array}{c}{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC}={\mathrm{\&lambda;}}_m\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_m^{NC}{(1-{\mathrm{Pt}}_m^{NC})}^{d-1}\\ {{\widetilde{\mathrm{\&lambda;}}}_l}^{NC}={\mathrm{\&lambda;}}_l\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_l^{NC}{(1-{\mathrm{Pt}}_l^{NC})}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of Home eNodeB-femtocell user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, with represent the probability of macro base station user transmission success and the probability of femtocell user transmission success under non-cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, (i=l or m) represent that data the d time retransmit just successfully probability, calculating formula is distinguished as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^{NC}=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{{\&Sigma;}_{n\&Element;N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_m\}\\ {\mathrm{Pt}}_l^{NC}=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_l}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\&NotEqual;n}{\overset{N_k}{\&Sigma;}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between Home eNodeB and macro base station user m, P _nrepresent the transmitting power of Home eNodeB n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent that macro base station user receives the Home eNodeB interference signal sum in same sub-channel, represent macro base station user m Signal to Interference plus Noise Ratio now;

| H _n,l| ²represent the channel gain between Home eNodeB n and femtocell user l, | H _{0, l}| ²represent the channel gain between macro base station and femtocell user l, P _nrepresent the transmitting power of Home eNodeB n, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between Home eNodeB i and femtocell user l, P _ifor the transmitting power of Home eNodeB, γ _lfor femtocell user Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that femtocell user receives, represent the interference signal that femtocell user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that femtocell user receives the interference signal of the same sub-channel of macro base station, expression femtocell user receives the interference signal of the Home eNodeB under same sub-channel, represent femtocell user l Signal to Interference plus Noise Ratio now.

2c), macro base station and potential alliance Home eNodeB estimate the user utility x after alliance _i ^*(S _n, Π _Ν) and alliance total utility ν ^*(S _n, Π _Ν), the alliance total utility ν after alliance ^*(S _n, Π _Ν) calculating formula is as follows:

This formula represents alliance S _nin the effectiveness sum of all users, | S _n| be the number of users of current alliance, represent the macro base station user utility after to cooperation in current alliance subregion and the summation of femtocell user effectiveness.Wherein, the user utility x after alliance _i ^*(S _n, Π _Ν) macro base station user utility x can be divided into _m* (S _n, Π _n) and femtocell user effectiveness x _l* (S _n, Π _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x_m}^{*}(S_n,{\&Pi;}_N)=\frac{{\mathrm{\&mu;}}_m^C{(S_n,{\&Pi;}_N)}^{\mathrm{\&delta;}}}{{\left(D_m^C\right)}^{(1-\mathrm{\&delta;})}}\\ {x_l}^{*}(S_n,{\&Pi;}_N)=\frac{{\mathrm{\&mu;}}_l^C{(S_n,{\&Pi;}_N)}^{\mathrm{\&delta;}}}{{\left(D_l^C\right)}^{(1-\mathrm{\&delta;})}}\end{array}\right.$

In formula, (S _n, Π _n) and (S _n, Π _n) be respectively reached at the information rate of macro base station user and reached at the information rate of femtocell user of cooperation under current alliance subregion, with be respectively the average delay of macro base station user and the average delay of femtocell user under cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, reached at the information rate of macro base station user under cooperation (S _n, Π _n) be calculated as follows:

${\mathrm{\&mu;}}_m^C(S_n,{\&Pi;}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\&Element;N_k\backslash S_n}{\&Sigma;}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2})$

This formula tabular form macro base station user m is subject to the Home eNodeB n outside alliance after alliance, n ∈ N _ks _ninterference time reached at information rate, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between Home eNodeB and macro base station user m, P _nrepresent the transmitting power of Home eNodeB n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives, namely macro base station user receives interference signal and the white Gaussian noise sum of the Home eNodeB outside alliance, represent the Signal to Interference plus Noise Ratio under macro base station user cooperation;

Reached at the information rate of femtocell user l under cooperation (S _n, Π _n) be calculated as follows:

${\mathrm{\&mu;}}_l^C(S_n,{\&Pi;}_N)=B_{l}log(1+\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\&NotEqual;n}{\overset{N_k\backslash S_n}{\&Sigma;}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2})$

This formula represents that femtocell user is subject to Home eNodeB n ∈ N outside alliance after alliance _ks _ninterference time reached at information rate, wherein, B _lrepresent the bandwidth of subchannel belonging to this femtocell user l, | H _n,l| ²represent the channel gain between Home eNodeB n and femtocell user l, P _nrepresent the transmitting power of Home eNodeB n, | H _i,l| ²represent the channel gain between Home eNodeB i and femtocell user l, P _ifor the transmitting power of Home eNodeB, | H _n,l| ²p _nrepresent the available signal power that femtocell user receives, namely femtocell user receives interference signal and the white Gaussian noise sum of the Home eNodeB outside alliance, for the Signal to Interference plus Noise Ratio of femtocell user under cooperation;

with be respectively the average delay of macro base station user and the average delay of femtocell user under cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^C=\frac{{\widetilde{\mathrm{\&lambda;}}}_m^C}{2{\mathrm{\&mu;}}_m^C({\mathrm{\&mu;}}_m^C-{\widetilde{\mathrm{\&lambda;}}}_m^C)}\\ D_l^C=\frac{{\widetilde{\mathrm{\&lambda;}}}_l^C}{2{\mathrm{\&mu;}}_l^C({\mathrm{\&mu;}}_l^C-{\widetilde{\mathrm{\&lambda;}}}_l^C)}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering under cooperation that maximum retransmit number is D with the practical communication load of femtocell user be respectively

$\left\{\begin{array}{c}{\widetilde{\mathrm{\&lambda;}}}_m^C={\mathrm{\&lambda;}}_m^C\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_m^C{(1-{\mathrm{Pt}}_m^C)}^{d-1}\\ {\widetilde{\mathrm{\&lambda;}}}_l^C={\mathrm{\&lambda;}}_l^C\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_l^C{(1-{\mathrm{Pt}}_l^C)}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of Home eNodeB-femtocell user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, with the probability of macro base station user transmission success and the probability of femtocell user transmission success under expression cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, calculating formula is respectively as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^C=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\&Element;N_k\backslash S_n}{\&Sigma;}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_m\}\\ {\mathrm{Pt}}_l^C=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\&NotEqual;n}{\overset{N_k\backslash S_n}{\&Sigma;}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between Home eNodeB and macro base station user m, P _nrepresent the transmitting power of Home eNodeB n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives,

Represent that macro base station user receives whole interference signals of Home eNodeB in same sub-channel outside alliance, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under macro base station user m cooperation;

| H _n,l| ²represent the channel gain between Home eNodeB n and femtocell user l, P _nrepresent the transmitting power of Home eNodeB n, | H _i,l| ²represent the channel gain between Home eNodeB i and femtocell user l, P _ifor the transmitting power of Home eNodeB, γ _lfor femtocell user Successful transmissions desired value, | H _n,l| ²p represents the available signal power that femtocell user receives, represent the interference signal of Home eNodeB in same sub-channel outside the alliance that femtocell user receives, represent the interference signal that femtocell user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under femtocell user l cooperation.

2d), as the user utility x after macro base station and Home eNodeB cooperation _i ^*(S _n, Π _Ν) and alliance total utility ν ^*(S _n, Π _Ν) be all not less than cooperation before effectiveness time, i.e. x _i ^*(S _n, Π _Ν) >x _i(S _n, Π _Ν) and ν ^*(S _n, Π _Ν) > ν (S _n, Π _Ν), then fed back by wire message way between macro base station and Home eNodeB and formally form alliance; Otherwise, judge whether macro base station and next potential cooperation Home eNodeB form new alliance;

2e), macro base station sequentially performs 2a, 2b, 2c, 2d with the Home eNodeB in all interference lists, according to the order of sequence until Stable coalitions is formed in all multiplexing subchannels.

Fig. 4 is frequency spectrum leasing model schematic provided by the invention, and sub-channel spectra multiplexing for Home eNodeB is normalized to a unit length by frequency spectrum leasing model, each unit is divided into three parts simultaneously:

Part I unit length is 1-α, and macro base station signal transmission gives affiliated macro base station user,

Part II unit length is α β, and unloaded macro base station user is given as relay transmission signal in micro-base station,

Part III unit length is α (1-β), and micro-base-station transmission signal gives affiliated micro-base station user,

Wherein, for frequency spectrum leasing coefficient represents that macro base station leases to the channel spectrum resource of micro-base station, for relay transmission window coefficient represents that micro-base station is the relaying frequency spectrum resource of the macro base station user distribution of unloading.

Effect of the present invention can be further illustrated by emulation:

1. simulated conditions:

Simulating scenes of the present invention covers by a macro base station and multiple Home eNodeB overlap the intensive isomery cellular network formed, the hexagon of macro base station coverage to be radius be 1Km, wherein dispose N number of Home eNodeB and M macro base station user, the border circular areas of Home eNodeB coverage to be radius be 20m, comprise one family base station user, we remain unchanged by default transmit power in the middle of whole coalition formation process.In down direction transmitting procedure, signal is mainly subject to by the impact about the factor of distance such as path loss, shadow fading.In addition, when Home eNodeB and macro base station user set up down direction be connected time, need the wall penetration loss considering 12dB more.In network, macro base station disposes 500 available downlink sub-carrier, each sub carries allocation 180KHz bandwidth, and belongs to different OFDMA subchannels.The simulation parameter according to 3GPP standard formulation is listed in table 1.Because Home eNodeB and all use are random placement per family, for eliminating the impact producing channel randomness, all emulated datas are all averaged through 10 circulations.

Table 1 isomery cellular network simulated environment optimum configurations

2. emulate content and result

As seen from Figure 5, increase with Home eNodeB density, within the specific limits, Home eNodeB is increased by cooperation related gain, but when exceeding certain value, relevant benefit reduces on the contrary, tracing it to its cause is because the Home eNodeB that nearly all macro base station downlink sub-channels have selected current optimum all is formed alliances, macro base station frequency range becomes crowded, no longer has unnecessary user and unloads to Home eNodeB.Generally speaking, adopt proposed alliance cooperative model in the region that family's base station deployment rate is high, higher income will be obtained.

Fig. 6 provides macro base station user gain effect and deployed position relation, draw the border area covered at macro base station, network condition is poor and interruption rate is high, if macro base station can find the Home eNodeB that can form alliance in this case, will obtain better effectiveness, we also can see that macro base station fringe region alliance quantity is maximum (namely having unloaded maximum macro base station users) simultaneously.

The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1. the business discharging method that in intensive isomery cellular network, resourceoriented is shared, it is characterized in that, it comprises the following steps:

S1, initialization network resource status;

S2, form Stable coalitions subregion towards multichannel;

Frequency spectrum leasing and business unloading in S3, alliance;

When S4, consideration change of network environment, periodically repeat above step, realize frequency spectrum leasing and business unloading in the alliance under multichannel Stable coalitions subregion.

2. the business discharging method that in a kind of intensive isomery cellular network according to claim 1, resourceoriented is shared, it is characterized in that, in initialization network resource status described in step S1, the partial frequency spectrum resource of the multiplexing macro base station in micro-base station, and macro base station and micro-base station keep fixing descending power signal transmission.

3. the business discharging method that in a kind of intensive isomery cellular network according to claim 2, resourceoriented is shared, it is characterized in that, described fixing descending power is maximum descending power.

4. the business discharging method that in a kind of intensive isomery cellular network according to claim 1, resourceoriented is shared, is characterized in that, described in step S2 towards the method step of multichannel formation Stable coalitions subregion is:

S21, towards multichannel build interference list;

S22, formation Stable coalitions.

5. the business discharging method that in a kind of intensive isomery cellular network according to claim 4, resourceoriented is shared, is characterized in that, described in step S21 towards the method for multichannel structure interference list is:

S21a, based on receive channel strength instruction (RSSI), macro base station detects the micro-base station be active in all downlink sub-channels k ∈ Κ (Κ is multiplexing number of sub-channels) be re-used by affiliated macro base station user, micro-base station also perception macro base station user in multiplexing subchannel;

Macro base station user m ∈ k in S21b, all multiplexing downlink sub-channels _m(k _mmacro base station number of users in multiplexing subchannel k) the micro-base station detected is disturbed by by force to weak sequence, form the interference list of respective sub-channel respectively, and feed back to macro base station by wireless channel;

The macro base station user perceived in multiplexing subchannel sorts by distance by S21c, micro-base station n ∈ N (N is micro-base station number) from the near to the remote.

6. the business discharging method that in a kind of intensive isomery cellular network according to claim 4, resourceoriented is shared, it is characterized in that, the method forming Stable coalitions described in step S22 is:

S22a, macro base station based on interference list, successively with the micro-base station n ∈ N in respective sub-channel k _k, N _k∈ N (N _kmicro-base station number for multiplexing subchannel k) consult in pairs, form new alliance, then as potential cooperation object when micro-base station n has a mind to cooperation;

The micro-base station n of S22b, macro base station and potential cooperation calculates the user utility x before alliance _i(S _n, Π _n), (i ∈ S _n, S _npotential alliance for this reason under channel, i is macro base station user m ∪ micro-base station user l ∈ L in potential alliance _n, L _nfor micro-base station user quantity that micro-base station is served, Π _nfor alliance's subregion) and alliance total utility ν (S _n, Π _n); Alliance total utility ν (S before alliance _n, Π _Ν) calculating formula is as follows:

Above formula represents the effectiveness sum of all users in potential cooperative alliances, | S _n| be number of users in alliance, represent the macro base station user utility in current potential cooperative alliances and the summation of micro-base station user effectiveness; Wherein, the user utility x before alliance _i(S _n, Π _n) macro base station user utility x can be divided into _m(S _n, Π _n) and micro-base station user effectiveness x _l(S _n, Π _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x}_m(S_n,{\mathrm{\&Pi;}}_N)=\frac{{\mathrm{\&mu;}}_m^{NC}{(S_n,{\mathrm{\&Pi;}}_N)}^{\mathrm{\&delta;}}}{{\left(D_m^{NC}\right)}^{(1-\mathrm{\&delta;})}}\\ x_l(S_n,{\mathrm{\&Pi;}}_N)=\frac{{\mathrm{\&mu;}}_l^{NC}{(S_n,{\mathrm{\&Pi;}}_N)}^{\mathrm{\&delta;}}}{{\left(D_l^{NC}\right)}^{(1-\mathrm{\&delta;})}}\end{array}\right.$

In formula, with be respectively reached at the information rate of macro base station user and reached at the information rate of micro-base station user of non-cooperation under current alliance subregion, S _npotential alliance for this reason under channel, Π _nfor current alliance subregion, with be respectively the average delay of macro base station user and the average delay of micro-base station user under non-cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, the macro base station user Ke Da information rate under non-cooperation be calculated as follows:

${\mathrm{\&mu;}}_m^{NC}(S_n,{\mathrm{\&Pi;}}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{{\mathrm{\&Sigma;}}_{n\&Element;N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2})$

Reached at information rate when this formula represents that macro base station user m is subject to micro-base station interference in same sub-channel, wherein, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the micro-base station interference signal sum in the same sub-channel that macro base station user receives, represent macro base station user m Signal to Interference plus Noise Ratio now;

Micro-base station user of non-cooperation can reach information rate micro-base station user quantity for micro-base station is served) be calculated as follows:

${\mathrm{\&mu;}}_l^{NC}(S_n,{\mathrm{\&Pi;}}_N)=B_{l}log(1+\frac{|H_{n,l}{|}^2{P}_n}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\&NotEqual;n}{\overset{N_k}{\mathrm{\&Sigma;}}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2})$

The reached information rate of this formula when micro-base station user l is subject to the interference of micro-base station in macro base station and same sub-channel before representing alliance, wherein, B _lrepresent the bandwidth of subchannel belonging to this micro-base station user l, | H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _{0, l}| ²represent the channel gain between macro base station and micro-base station user l, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, | H _n,l| ²p _nrepresent the available signal power that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that micro-base station user l receives the interference signal of the same sub-channel of macro base station, represent that micro-base station user l receives the interference signal of the micro-base station i under same sub-channel, represent micro-base station user l Signal to Interference plus Noise Ratio now;

with be respectively the average delay of macro base station user and the average delay of micro-base station user under non-cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^{NC}=\frac{{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC}}{2{\mathrm{\&mu;}}_m^{NC}({\mathrm{\&mu;}}_m^{NC}-{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC})}\\ D_l^{NC}=\frac{{{\widetilde{\mathrm{\&lambda;}}}_l}^{NC}}{2{\mathrm{\&mu;}}_l^{NC}({\mathrm{\&mu;}}_l^{NC}-{{\widetilde{\mathrm{\&lambda;}}}_l}^{NC})}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering that maximum retransmit number is D under non-cooperation with the practical communication load of micro-base station user be respectively

$\left\{\begin{array}{c}{{\widetilde{\mathrm{\&lambda;}}}_m}^{NC}={\mathrm{\&lambda;}}_m\underset{d=1}{\overset{D}{\&Sigma;}}P{t}_m^{NC}{\left(1-{\mathrm{Pt}}_m^{NC}\right)}^{d-1}\\ {{\widetilde{\mathrm{\&lambda;}}}_l}^{NC}={\mathrm{\&lambda;}}_l\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_l^{NC}{(1-{\mathrm{Pt}}_l^{NC})}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of micro-base station-micro-base station user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, Pt _m ^nCand Pt _l ^nCrepresent the probability of macro base station user transmission success and the probability of micro-base station user transmission success under non-cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, (i=l or m) represent that data the d time retransmit just successfully probability, calculating formula is distinguished as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^{NC}=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{{\mathrm{\&Sigma;}}_{n\&Element;N_k}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_m\}\\ {\mathrm{Pt}}_l^{NC}=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_l}{|H_{0,l}{|}^2{P}_0+\underset{i=1,i\&NotEqual;n}{\overset{N_k}{\mathrm{\&Sigma;}}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that macro base station user receives, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent that macro base station user receives the micro-base station interference signal sum in same sub-channel, represent macro base station user m Signal to Interference plus Noise Ratio now;

| H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, | H _{0, l}| ²represent the channel gain between macro base station and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, P ₀represent the transmitting power of macro base station, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, γ _lfor micro-base station user Successful transmissions desired value, | H _{0, m}| ²p ₀represent the available signal power that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, | H _{0, l}| ²p ₀represent that micro-base station user receives the interference signal of the same sub-channel of macro base station, represent that micro-base station user receives the interference signal of the micro-base station under same sub-channel, represent micro-base station user l Signal to Interference plus Noise Ratio now;

The user utility x after alliance is estimated in S22c, the micro-base station of macro base station and potential alliance _i ^*(S _n, Π _Ν) and alliance total utility ν ^*(S _n, Π _Ν); Alliance total utility ν after alliance ^*(S _n, Π _Ν) calculating formula is as follows:

This formula represents alliance S _nin the effectiveness sum of all users, | S _n| be the number of users of current alliance, represent the macro base station user utility after to cooperation in current alliance subregion and the summation of micro-base station user effectiveness; Wherein, the user utility after alliance macro base station user utility can be divided into with micro-base station user effectiveness be calculated as follows respectively:

$\left\{\begin{array}{c}{x_m}^{*}(S_n,{\mathrm{\&Pi;}}_N)=\frac{{\mathrm{\&mu;}}_m^C{(S_n,{\mathrm{\&Pi;}}_N)}^{\mathrm{\&delta;}}}{{\left(D_m^C\right)}^{(1-\mathrm{\&delta;})}}\\ {x_l}^{*}(S_n,{\mathrm{\&Pi;}}_N)=\frac{{\mathrm{\&mu;}}_l^C{(S_n,{\mathrm{\&Pi;}}_N)}^{\mathrm{\&delta;}}}{{\left(D_l^C\right)}^{(1-\mathrm{\&delta;})}}\end{array}\right.$

In formula, with be respectively reached at the information rate of macro base station user and reached at the information rate of micro-base station user of cooperation under current alliance subregion, with be respectively the average delay of macro base station user and the average delay of micro-base station user under cooperation, δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, wherein, reached at the information rate of macro base station user under cooperation be calculated as follows:

${\mathrm{\&mu;}}_m^C(S_n,{\mathrm{\&Pi;}}_N)=Blog(1+\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\&Element;N_k\backslash S_n}{\&Sigma;}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2})$

This formula tabular form macro base station user m is subject to micro-base station n, the n ∈ N outside alliance after alliance _ks _ninterference time reached at information rate, B represents the bandwidth of subchannel belonging to this macro base station user m, the logarithm operation at log to be 10 the be end, | H _{0, m}| ²represent the channel gain between macro base station (subscript 0 represents macro base station) and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives, namely macro base station user receives interference signal and the white Gaussian noise sum of the micro-base station outside alliance, represent the Signal to Interference plus Noise Ratio under macro base station user cooperation;

Reached at the information rate of micro-base station user l under cooperation be calculated as follows:

${\mathrm{\&mu;}}_l^C(S_n,{\mathrm{\&Pi;}}_N)=B_{l}log(1+\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\&NotEqual;n}{\overset{N_k\backslash S_n}{\&Sigma;}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2})$

This formula represents that micro-base station user is subject to micro-base station n ∈ N outside alliance after alliance _ks _ninterference time reached at information rate, wherein, B _lrepresent the bandwidth of subchannel belonging to this micro-base station user l, | H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, | H _n,l| ²p _nrepresent the available signal power that micro-base station user receives, namely micro-base station user receives interference signal and the white Gaussian noise sum of the micro-base station outside alliance, for the Signal to Interference plus Noise Ratio of micro-base station user under cooperation;

with be respectively the average delay of macro base station user and the average delay of micro-base station user under cooperation, calculating formula is as follows:

$\left\{\begin{array}{c}{D}_m^C=\frac{{\widetilde{\mathrm{\&lambda;}}}_m^C}{2{\mathrm{\&mu;}}_m^C({\mathrm{\&mu;}}_m^C-{\widetilde{\mathrm{\&lambda;}}}_m^C)}\\ D_l^C=\frac{{\widetilde{\mathrm{\&lambda;}}}_l^C}{2{\mathrm{\&mu;}}_l^C({\mathrm{\&mu;}}_l^C-{\widetilde{\mathrm{\&lambda;}}}_l^C)}\end{array}\right.$

Wherein, the practical communication load of macro base station user when considering under cooperation that maximum retransmit number is D with the practical communication load of micro-base station user be respectively

$\left\{\begin{array}{c}{\widetilde{\mathrm{\&lambda;}}}_m^C={\mathrm{\&lambda;}}_m^C\underset{d=1}{\overset{D}{\&Sigma;}}P{t}_m^C{\left(1-{\mathrm{Pt}}_m^C\right)}^{d-1}\\ {\widetilde{\mathrm{\&lambda;}}}_l^C={\mathrm{\&lambda;}}_l^C\underset{d=1}{\overset{D}{\&Sigma;}}{\mathrm{Pt}}_l^C{(1-{\mathrm{Pt}}_l^C)}^{d-1}\end{array}\right.$

λ _mand λ _l(bits/s) represent the average arrival rate of macro base station-macro base station user and the average arrival rate of micro-base station-micro-base station user by being determined grouping during M/D/1 queuing model by poisson arrival process respectively, d is current re-transmission number, with the probability of macro base station user transmission success and the probability of micro-base station user transmission success under expression cooperation respectively, namely when signal to noise ratio (SINR) is higher than set respective objects value γ _m, γ _lprobability, calculating formula is respectively as follows:

$\left\{\begin{array}{c}{\mathrm{Pt}}_m^{NC}=\mathrm{Pr}\{\frac{|H_{0,m}{|}^2{P}_0}{\underset{n\&Element;N_k\backslash S_n}{\mathrm{\&Sigma;}}|H_{n,m}{|}^2{P}_n+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_m\}\\ {\mathrm{Pt}}_l^C=\mathrm{Pr}\{\frac{|H_{n,l}{|}^2{P}_n}{\underset{i=1,i\&NotEqual;n}{\overset{N_k\backslash S_n}{\mathrm{\&Sigma;}}}|H_{i,l}{|}^2{P}_i+{\mathrm{\&sigma;}}^2}\&GreaterEqual;{\mathrm{\&gamma;}}_l\}\end{array}\right.$

Wherein, Pr{SINR>=γ } when representing that Signal to Interference plus Noise Ratio is greater than certain target value gamma, the probability distribution of data Successful transmissions, | H _{0, m}| ²represent the channel gain between macro base station and macro base station user m, P ₀represent the transmitting power of macro base station for its user m, | H _n,m| ²represent the channel gain between micro-base station and macro base station user m, P _nrepresent the transmitting power of micro-base station n, σ ²for white Gaussian noise (AWGN) mean-square value, γ _mfor grand user's Successful transmissions desired value, | H _{0, m}| ²p ₀for the available signal power that macro base station user receives, represent that macro base station user receives whole interference signals of micro-base station in same sub-channel outside alliance, represent the interference signal that macro base station user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under macro base station user m cooperation;

| H _n,l| ²represent the channel gain between micro-base station n and micro-base station user l, P _nrepresent the transmitting power of micro-base station n, | H _i,l| ²represent the channel gain between micro-base station i and micro-base station user l, P _ifor the transmitting power of micro-base station, γ _lfor micro-base station user Successful transmissions desired value, | H _n,l| ²p represents the available signal power that micro-base station user receives, represent the interference signal of micro-base station in same sub-channel outside the alliance that micro-base station user receives, represent the interference signal that micro-base station user receives and white Gaussian noise sum, represent the Signal to Interference plus Noise Ratio under micro-base station user l cooperation;

S22d, as the user utility x after macro base station and micro-base station cooperation _i ^*(S _n, Π _n) and alliance total utility ν ^*(S _n, Π _n) be all not less than cooperation before effectiveness time, i.e. x _i ^*(S _n, Π _n) >x _i(S _n, Π _n) and ν ^*(S _n, Π _n) > ν (S _n, Π _n), then fed back by wire message way between macro base station and micro-base station and formally form alliance; Otherwise, judge whether macro base station and next micro-base station of potential cooperation form new alliance;

S22e, macro base station sequentially perform S22a, S22b, S22c, S22d with the micro-base station in all interference lists according to the order of sequence, until Stable coalitions is formed in all multiplexing subchannels.

7. the business discharging method that in a kind of intensive isomery cellular network according to claim 6, resourceoriented is shared, it is characterized in that, alliance's total utility that in arbitrary subchannel, micro-base station departs from the new alliance subregion of current coalition formation is less than alliance's total utility of current steady alliance subregion.

8. the business discharging method that in a kind of intensive isomery cellular network according to claim 1, resourceoriented is shared, is characterized in that, in alliance described in step S3, the method for frequency spectrum leasing and business unloading is:

S31, based on the Stable coalitions subregion formed, according to frequency spectrum leasing model, solved the optimal solution of optimum frequency spectrum leasing coefficient by convex optimization tool;

S32, macro base station lease to corresponding micro-BTS channel frequency spectrum resource according to the optimum frequency spectrum leasing coefficient calculating gained, the macro base station user simultaneously notifying to take this channel disconnects former link and sets up new downlink transfer with corresponding micro-base station and links, and the frequency spectrum resource of the relaying window coefficient ratio in gained frequency spectrum resource is distributed to macro base station user and is used for business transmission by corresponding micro-base station.

9. the business discharging method that in a kind of intensive isomery cellular network according to claim 8, resourceoriented is shared, it is characterized in that, the method solving the optimal solution of optimum frequency spectrum leasing coefficient described in step S31 is:

Sub-channel spectra multiplexing for micro-base station is normalized to a unit length by frequency spectrum leasing model, each unit is divided into three parts simultaneously:

Part I unit length is 1-α, and macro base station signal transmission gives affiliated macro base station user;

Part II unit length is α β, and unloaded macro base station user is given as relay transmission signal in micro-base station;

Part III unit length is α (1-β), and micro-base-station transmission signal gives affiliated micro-base station user;

Wherein, α is that frequency spectrum leasing coefficient represents that macro base station leases to the channel spectrum resource of micro-base station, and β is that relay transmission window coefficient represents that micro-base station is the relaying frequency spectrum resource of the macro base station user distribution of unloading;

According to frequency spectrum leasing model, the utility function under this model is proposed, and by the convex optimization tool optimal solution that to solve with frequency spectrum leasing factor alpha and relay transmission window coefficient β be variable:

$\left\{\begin{array}{cc}\underset{\mathrm{\&alpha;},\mathrm{\&beta;}}{\max }& x_i^{\&prime;}(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})\\ \underset{\mathrm{\&alpha;},\mathrm{\&beta;}}{\max }& v^{\&prime;}(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})\\ s.t.& 0<\mathrm{\&alpha;},\mathrm{\&beta;}<1\end{array}\right.$

Above formula represents and is meeting 0< α, during the condition of β <1, solves the maximum user utility x based on frequency spectrum leasing model in Stable coalitions subregion _i' (S ' _n, Π ' _n) and maximum alliance total utility ν ' (S ' _n, Π ' _n), S ' _nfor the alliance under current steady alliance subregion, Π ' _nfor current steady alliance subregion; Wherein, user utility x ' _i(S ' _n, Π ' _n) macro base station user utility x ' can be divided into _m(S ' _n, Π ' _n) and micro-base station user effectiveness x ' _l(' ' _n, Π ' _n), be calculated as follows respectively:

$\left\{\begin{array}{c}{x}_m^{\&prime;}(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})=\frac{{\mathrm{\&alpha;\&beta;\&mu;}}_m^{\&prime;}{(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})}^{\mathrm{\&delta;}}}{{D_m^{\&prime;}}^{(1-\mathrm{\&delta;})}})\\ x_l^{\&prime;}(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})=\frac{\mathrm{\&alpha;}(1-\mathrm{\&beta;}){\mathrm{\&mu;}}_l^{\&prime;}{(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})}^{\mathrm{\&delta;}}}{{D_l^{\&prime;}}^{(1-\mathrm{\&delta;})}}\end{array}\right.$

With

δ ∈ (0,1) is transmission capacity-time delay balance coefficient, and namely system is to the tolerance of propagation delay time, α and β is respectively frequency spectrum leasing coefficient and relay transmission window coefficient; Wherein, reached at the information rate μ ' of the macro base station user in optimization problem _m(α, β, S ' _n, Π ' _n) and reached at the information rate μ ' of micro-base station user _l(α, β, S ' _n, Π ' _n) calculating formula is as follows:

$\left\{\begin{array}{c}{\mathrm{\&mu;}}_m^{\&prime;}(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})=\mathrm{\&alpha;}\mathrm{\&beta;}{\mathrm{\&mu;}}_m^C(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})\\ {\mathrm{\&mu;}}_l^{\&prime;}(\mathrm{\&alpha;},\mathrm{\&beta;},S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})=\mathrm{\&alpha;}(1-\mathrm{\&beta;}){\mathrm{\&mu;}}_l^C(S_n^{\&prime;},{\mathrm{\&Pi;}}_N^{\&prime;})\end{array}\right.$

Represent macro base station user and the reached information rate of micro-base station user under frequency spectrum leasing model under current alliance subregion, with be that macro base station user Ke Da information rate under current alliance subregion under cooperation and micro-base station user can reach information rate respectively, calculating formula is as shown in S22c.

10. the business discharging method that in a kind of intensive isomery cellular network according to claim 1, resourceoriented is shared, it is characterized in that, processing method when considering change of network environment described in step S4 is: when change of network environment, periodically repeat step S2 towards frequency spectrum leasing and business unloading in the alliance that multichannel forms Stable coalitions subregion and step S3; When network environment does not change, periodically repeat the initialization network resource status of step S1, S2 towards frequency spectrum leasing and business unloading in the alliance that multichannel forms Stable coalitions subregion and step S3.