CN103517438A - Power distribution method and device in cooperative communication system - Google Patents

Abstract

The invention provides a power distribution method and device in a cooperative communication system. The method comprises the steps of acquiring a first relational expression between the total sending power of a source node and each relay node, the ratio between the sending power of each relay node and the sending power of the source node, and the channel state parameter of the system; calculating the ratio between the sending power of each relay node and the sending power of the source node when the total sending power is a minimum value according to the first relational expression as the distribution information of the sending power of each relay node and the sending power of the source node; and performing power distribution by using the distribution information. According to the invention, the flexible configuration of the sending power at each sending end of a multi-relay cooperative communication system.

Description

The distribution method of power and device in Turbo Detection for Cooperative Communication

Technical field

The present invention relates to wireless communication technology field, refer to especially distribution method and the device of power in a kind of Turbo Detection for Cooperative Communication.

Background technology

Multiple-input and multiple-output (MIMO) system is because the transmitting terminal in system and receiving terminal have configured many antennas, thereby obtained obvious space diversity and spatial multiplexing gain, improved the capacity of system.Yet along with the development of integrated circuit technique and manufacturing process, it is more and more microminiaturized that mobile terminal becomes, and due to the restriction of size and manufacturing cost, will overcome the impact of coherence distance, becomes very difficult at mobile terminal configuration multiple antenna simultaneously.

In recent years, in multi-user communication environment, adopt communication for coordination technology, make to configure respectively closing on mobile subscriber and can sharing antenna each other, collaborative transmitted signal according to certain way of single slave antenna, produce a kind of virtual environment that many antennas send that is similar to, obtain space diversity gain, improve the transmission performance of system.

Amplification forwarding (amplify-and-forward, AF) is the simplest a kind of cooperative transmission mode in Turbo Detection for Cooperative Communication.In amplification forwarding mode, via node amplifies to the received signal, is then transmitted to destination node, and destination node merges the signal from source node and each via node receiving, and with this, obtains diversity gain.Amplification forwarding mode complexity is lower, be easy to realize, more practical in actual LTE-A/4G communication system.

Yet the power distribution methods such as existing Turbo Detection for Cooperative Communication employing distribute the transmitted power equating between source node and each via node.Although the method is relatively simple,, owing to not considering each node each opposite sex in network, can not make full use of system resource, cannot realize transmitted power in the flexible configuration of each transmitting terminal of many relay cooperative communication systems.

Summary of the invention

The technical problem to be solved in the present invention is to provide distribution method and the device of power in a kind of Turbo Detection for Cooperative Communication, can realize transmitted power in the flexible configuration of each transmitting terminal of many relay cooperative communication systems.

For solving the problems of the technologies described above, embodiments of the invention provide technical scheme as follows:

On the one hand, provide the distribution method of power in a kind of Turbo Detection for Cooperative Communication, comprising:

Obtain the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

According to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node;

Use described assignment information to carry out power division.

The channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{\mathrm{\α}}_k){[\frac{3}{4m^2}\·\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\Π}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\α}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

A _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

Described according to described the first relational expression, when calculating described total transmitted power and being minimum value described in the step of ratio between the transmitted power of each via node and the transmitted power of described source node comprise:

Set total transmitted power P _{total, PAp}minimum border I He Gao border u, the primary vector α of a K dimension and the secondary vector α of K dimension ⁰;

The random initial value α that produces primary vector α ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

According to described the first expression formula, calculate corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Get current iteration frequency n=1;

Carry out iteration, produce at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

According to the 3rd vectorial β, calculate the first variable Δ α, be specially: $\mathrm{\Δ\α}=(u-1)*g_{\mathrm{\μ}}^{-1}\left(\mathrm{\β}\right);$ Wherein, $g_{\mathrm{\μ}}^{-1}\left(\mathrm{\β}\right)=\frac{{(1+\mathrm{\μ})}^{\left|\mathrm{\β}\right|}-1}{\mathrm{\μ}}\mathrm{sign}\left(\mathrm{\β}\right);$ Sign (β) is sign function; The second variable μ=10 ¹⁰⁰(n/N) ^q;

According to the first value of variable Δ α and the currency of primary vector α, generate the four-way amount α ' in interval (I, u), be specially: α '=α+Δ α;

According to described the first relational expression, calculate the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

If Δ P=P (α ')-P (α) < 0 or z < p=exp are ((n/N) ^qΔ P/|p (α) |/ε _f), use the currency of the 3rd vectorial α ' to upgrade primary vector α; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0;

If P (α) < P is (α ⁰), use the currency of primary vector α to upgrade secondary vector α ⁰; ,

Whether the currency that judges current iteration frequency n equals maximum iteration time N;

If YES, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node; If NO, current iteration frequency n adds 1, proceeds iteration.

Described maximum iteration time N equates with the total quantity K of described via node.

Described according to described the first expression formula, calculate corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰) step comprise:

By measuring, obtain power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

According to modulation system, obtain the value of modulation system constant m;

According to target error rate, error rate of system P is set _evalue;

According to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

On the other hand, provide the distributor of power in a kind of Turbo Detection for Cooperative Communication, comprising:

Acquiring unit, obtains the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

Computing unit, according to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node;

Allocation units, are used described assignment information to carry out power division.

The channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

Described computing unit comprises:

First sets subelement, sets total transmitted power P _{tatal, PAp}minimum border I He Gao border u;

Second sets subelement, sets the primary vector α of K dimension and the secondary vector α of K dimension ⁰,

The first numerical value produces subelement, produces at random the initial value α of primary vector α ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

The first computation subunit, according to described the first expression formula, calculates corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Value subelement, gets current iteration frequency n=1;

Second value produces subelement, while carrying out iteration, produces at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

The second computation subunit, calculates the first variable Δ α according to the 3rd vectorial β, is specially: $\mathrm{\&Delta;\&alpha;}=(u-1)*g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right);$ Wherein, $g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)=\frac{{(1+\mathrm{\&mu;})}^{\left|\mathrm{\&beta;}\right|}-1}{\mathrm{\&mu;}}\mathrm{sign}\left(\mathrm{\&beta;}\right);$ Sign (β) is sign function; The second variable μ=10 ¹⁰⁰(n/N) ^q;

The 3rd computation subunit, generates the four-way amount α ' in interval (I, u) according to the first value of variable Δ α and the currency of primary vector α, is specially: α '=α+Δ α;

The 4th computation subunit, according to described the first relational expression, calculates the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

First upgrades subelement, when Δ P=P (α ')-P (α) < 0 or z < p=exp ((n/N) ^qΔ P/|P (α) |/ε _f) time, use the currency of the 3rd vectorial α ' to upgrade primary vector α; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0;

Second upgrades subelement, as P (α) < P (α ⁰) time, use the currency of primary vector α to upgrade secondary vector α ⁰;

Judgment sub-unit, judges whether the currency of current iteration frequency n equals maximum iteration time N;

The 5th operator unit, when the currency of current iteration frequency n equals maximum iteration time N, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node; If NO, current iteration frequency n adds 1, proceeds iteration.

Described maximum iteration time N equates with the total quantity K of described via node.

Described the first computation subunit comprises:

Measurement module, by measuring, obtains power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

Acquisition module, according to modulation system, obtains the value of modulation system constant m;

Module is set, according to target error rate, error rate of system P is set _evalue;

Computing module, according to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

Embodiments of the invention have following beneficial effect:

In such scheme, according to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, assignment information as the transmitted power of described each via node and the transmitted power of described source node, distribute the constant power distribution method of identical transmitted power to compare with source node in present technology and each via node, consider each node each opposite sex in network, realized the flexible configuration of transmitted power at each transmitting terminal of many relay cooperative communication systems.

Accompanying drawing explanation

Fig. 1 is the schematic flow sheet of the embodiment of the distribution method of power in Turbo Detection for Cooperative Communication of the present invention;

Fig. 2 is the schematic flow sheet of another embodiment of the distribution method of power in Turbo Detection for Cooperative Communication of the present invention;

Fig. 3 is the schematic flow sheet of another embodiment of the distribution method of power in Turbo Detection for Cooperative Communication of the present invention;

Fig. 4 is the structural representation of the distributor of power in Turbo Detection for Cooperative Communication of the present invention.

Embodiment

For technical problem, technical scheme and advantage that embodiments of the invention will be solved are clearer, be described in detail below in conjunction with the accompanying drawings and the specific embodiments.

As shown in Figure 1, the schematic flow sheet for the embodiment of the distribution method of power in a kind of Turbo Detection for Cooperative Communication of the present invention, comprising:

Step 11, obtains the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

Step 12, according to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node;

Step 13, is used described assignment information to carry out power division.

Wherein, the channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

A _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

Step 12 comprises:

Step 121, sets total transmitted power P _{total, PAp}minimum border I He Gao border u, the primary vector α of a K dimension and the secondary vector α of K dimension ⁰;

Step 122, produces the initial value α of primary vector α at random ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

Step 123, according to described the first expression formula, calculates corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Step 124, gets current iteration frequency n=1;

Step 125, carries out iteration, produces at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

Step 126, calculates the first variable Δ α according to the 3rd vectorial β, is specially: $\mathrm{\&Delta;\&alpha;}=(u-1)*g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right);$ Wherein, $g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)=\frac{{(1+\mathrm{\&mu;})}^{\left|\mathrm{\&beta;}\right|}-1}{\mathrm{\&mu;}}\mathrm{sign}\left(\mathrm{\&beta;}\right);$ Sign (β) is sign function; The second variable μ=10 ¹⁰⁰(n/N) ^q;

Step 127, generates the four-way amount α ' in interval (I, u) according to the first value of variable Δ α and the currency of primary vector α, is specially: α '=α+Δ α;

Step 128, according to described the first relational expression, calculates the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

Step 129, judges for the first time;

Step 130, if Δ P=P (α ')-P (α) < 0 or z < p=exp are ((n/N) ^qΔ P/|P (α) |/ε _f), use the currency of the 3rd vectorial α ' to upgrade primary vector α; Otherwise, primary vector α is not upgraded; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0; Described maximum iteration time N equates with the total quantity K of described via node.

Step 131, judges for the second time;

Step 132, if P (α) < P is (α ⁰), use the currency of primary vector α to upgrade secondary vector α ⁰; Otherwise, not to secondary vector α ⁰upgrade;

Step 133, judges whether the currency of current iteration frequency n equals maximum iteration time N;

Step 134, if YES, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node;

Step 135, if NO, current iteration frequency n adds 1, proceeds iteration, performs step 125.

In above-mentioned treatment step, treat estimation model P _{total, PAp}minimum border I He Gao border u, maximum iterations N, annealing factor q, error margins ε _fvalue can set in advance.

Wherein, step 123 comprises:

Step 1231, by measuring, obtains power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

Step 1232, according to modulation system, obtains the value of modulation system constant m;

Step 1233, arranges error rate of system P according to target error rate _evalue;

Step 1234, according to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

The present invention proposes in a kind of Turbo Detection for Cooperative Communication based on simulated annealing (Simulated Annealing, SA) power distribution method, the method is in the situation that guaranteeing that systematic function reaches target error rate requirement, total transmitted power (the transmitted power sum of the transmitted power of source node and each via node) of many relayings amplification forwarding Turbo Detection for Cooperative Communication is minimized, there is relatively low computation complexity and network signal load simultaneously.That is to say, when communication for coordination, the method distributes the constant power distribution method of identical transmitted power to compare with source node in present technology and each via node, in the situation that reaching same error rate of system requirement, can obviously reduce total transmitted power of many relayings amplification forwarding Turbo Detection for Cooperative Communication, improve the power system capacity of Turbo Detection for Cooperative Communication.The universality that the present invention has and robustness make it at LTE-A/4G(Long Term Evolution, Long Term Evolution) there is very high using value in system.

The distribution method of transmitted power in communication for coordination of the present invention, the requirement of the Offered target error rate in Turbo Detection for Cooperative Communication, obtains source node and the transmitted power sum of each via node and the relation between channel state parameter and error rate of system in system; Then adopt simulated annealing, when error rate of system performance reaches target error rate and requires, the transmit power assignment situation of source node and each via node when iteration is obtained the total transmitted power of system and reached minimum value.Then, each via node and the source node transmitted power to be assigned to, sends information by the mode of direct mode or amplification forwarding to destination node.

Fig. 3 shows another implementing procedure figure of the power distribution method in communication for coordination of the present invention, comprising:

Step 1 according to the requirement of Turbo Detection for Cooperative Communication, first, arranges source node (S), destination node (D), a K via node (R in this system _k); K via node all participates in cooperative transmission; Then, according to the positional information of source node, via node and destination node, measure the channel condition information between source node and each via node, each via node and destination node, source node and destination node; Obtain the calculating formula of the error rate of system:

$P_e\&ap;\frac{3}{4m^2}\&CenterDot;\frac{1}{{\stackrel{\&OverBar;}{\mathrm{\&gamma;}}}_f}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{{\stackrel{\&OverBar;}{\mathrm{\&gamma;}}}_{g_k}}+\frac{1}{{\stackrel{\&OverBar;}{\mathrm{\&gamma;}}}_{h_k}}),$

Wherein, P _eerror rate of system for many relayings of Turbo Detection for Cooperative Communication AF cooperative transmission;

M is constant, by the modulation system adopting, is determined, for example, adopt BPSK(Binary Phase Shift Keying, biphase phase shift keying) during modulation, m=2;

for the instantaneous signal-to-noise ratio γ of source node to channel between destination node _fmean value;

instantaneous signal-to-noise ratio for channel between source node to a k via node

mean value;

be that k via node is to the instantaneous signal-to-noise ratio of channel between destination node mean value;

K is the total quantity of via node; K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K.

Step 2, according to the said system error rate, obtains the instantaneous signal-to-noise ratio of each branch road of direct transmission and cooperative transmission; Obtain the relation between the channel state parameter between the Turbo Detection for Cooperative Communication error rate and source node and destination node, source node and k via node, a k via node and destination node; Thereby obtain source node and the transmitted power sum of each via node and the calculation relational expression between channel state parameter and error rate of system in system;

Step 3, transmitted power sum based on source node and each via node in above-mentioned Turbo Detection for Cooperative Communication and the calculation relational expression between channel state parameter and error rate of system, according to the relational expression between the error rate of system calculating in step 2 and each channel state parameter, adopt SA algorithm, source node when iterative computation goes out the total transmitted power of system and reaches minimum value and the transmit power assignment situation of each via node;

Step 4, according to the power division situation calculating in step 3, source node and each via node according to the transmitted power being assigned to directly to transmit or the mode of amplification forwarding sends information to destination node.

The instantaneous signal-to-noise ratio of each branch road that obtains direct transmission and cooperative transmission described in above-mentioned steps 2, the concrete steps of obtaining the relation between the channel state parameter between the Turbo Detection for Cooperative Communication error rate and source node and destination node, source node and k via node, a k via node and destination node are as follows:

Step 21, obtains the calculation expression of the instantaneous signal-to-noise ratio of direct transmission and each transmission branch of cooperative transmission, and it is specially: suppose that directly transmission equates with the noise power spectral density of each transmission branch of cooperative transmission, its calculation expression is respectively:

γ _f＝|f| ² P ₀/N ₀； （2）

${\mathrm{\&gamma;}}_{g_k}={\left|g_k\right|}^2{P}_0/N_0;---\left(3\right)$

${\mathrm{\&gamma;}}_{h_k}={\left|h_k\right|}^2{P}_k/N_0;---\left(4\right)$

Wherein, γ _f,

be respectively the instantaneous signal-to-noise ratio of source node on to destination node, source node to a k via node and k via node to the communication leg of destination node;

| f| ², | g _k| ²with | h _k| ²be respectively the channel gain between channel gain, a k via node and the destination node between channel gain, source node and k the via node between source node and destination node;

P ₀and P _kbe respectively the transmitted power of source node and k via node;

N ₀power spectral density for noise in channel;

K is the sequence number of via node, k=1 ..., K;

K is the total quantity of via node; K is greater than 0 natural number.

Step 22, obtains the relational expression between the channel state parameter between the Turbo Detection for Cooperative Communication error rate and source node and destination node, source node and k via node, a k via node and destination node. and it is specific as follows:

By γ _f,

expression formula (2) (3) (4) be updated to above-mentioned AF collaborative transmission system error rate P _eexpression formula (1) is inner, and the expression formula obtaining is as follows:

$P_e\&ap;\frac{3}{4m^2}\&CenterDot;\frac{N_0^{K+1}}{P_0}\&CenterDot;\frac{1}{F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{P_0{G}_k}+\frac{1}{P_k{H}_k}),---\left(5\right)$

Wherein, P _efor error rate of system;

M is constant, by the modulation system adopting, is determined;

N ₀power spectral density for noise in channel;

P ₀and P _kbe respectively the transmitted power of source node and k via node;

F, G _kwith H _kbe respectively the channel gain between source node and destination node, source node and k via node, a k via node and destination node | f| ², | g _k| ²with | h _k| ²mean value;

K is the sequence number of via node, k=1 ..., K;

K is the total quantity of via node; K is greater than 0 natural number.

Step 23, according to AF collaborative transmission system error rate P obtained above _eand the relation between system each several part parameter, obtains source node and the transmitted power sum of each via node and the calculation relational expression between channel state parameter and error rate of system in system.It is specific as follows:

Suppose error rate of system P _emeet given target error rate requirement, and total transmitted power of system is variable, total transmitted power of the transmitted power of k via node and system just can be expressed as so:

P _k＝α _kP ₀； （6）

$P_{\mathrm{total},\mathrm{RAP}}=P_0+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{P}_K;---\left(7\right)$

Wherein, P ₀and P _kbe respectively the transmitted power of source node and k via node;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node;

P _{toatl, RAP}total transmitted power for system;

K is the sequence number of via node, k=1 ..., K;

K is the total quantity of via node; K is greater than 0 natural number.

Then, by P _kcalculating formula (6) be updated to the error rate of system P in step 22 _ecalculating formula (5), obtain error rate P _eand the relation between the transmitted power of system source node, its expression formula is as follows:

$P_0{\&ap;[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_e}\&CenterDot;\frac{1}{F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}}---\left(8\right)$

Wherein, P ₀transmitted power for source node;

M is constant, by the modulation system adopting, is determined;

N ₀power spectral density for noise in channel;

P _efor error rate of system;

F, G _kwith H _kbe respectively the channel gain between source node and destination node, source node and k via node, a k via node and destination node | f| ², | g _k| ²with | h _k| ²mean value;

α _kbe the ratio between the transmitted power of k via node and the transmitted power of source node;

K is the sequence number of via node, k=1 ..., K;

K is the total quantity of via node; K is greater than 0 natural number.

Finally, by P ₀expression formula (8) be updated in the calculating formula (7) of the total transmitted power of system, obtain source node and the transmitted power sum of each via node and the calculating formula between channel state parameter and error rate of system in system, as follows:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}}---\left(9\right)$

Wherein, P _{toatl, RAP}total transmitted power for system;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node;

M is constant, by the modulation system adopting, is determined;

N ₀power spectral density for noise in channel;

P _eit is error rate of system;

F, G _kwith H _kbe respectively the channel gain between source node and destination node, source node and k via node, a k via node and destination node | f| ², | g _k| ²with | h _k| ²mean value;

K is the sequence number of via node, k=1 ..., K;

K is the total quantity of via node; K is greater than 0 natural number.

Described in above-mentioned steps 3 according to the calculating formula between the transmitted power sum of source node and each via node in the system obtaining and channel state parameter and error rate of system, adopt SA algorithm, thereby when the total transmitted power of the system that calculates reaches minimum value, the implementation step of the power allocation scheme of source node and each via node is as follows:

Step 31, sets primary vector α=[α ₁, α ₂..., α _k], secondary vector α ⁰, secondary vector and secondary vector are all K dimensional vector; α _kbe the ratio between the transmitted power of k via node and the transmitted power of source node; K is the sequence number of via node, k=1 ..., K; K is the total quantity of via node; K is greater than 0 natural number.

Initialization.The random initial value α that produces α ₀, determine to treat estimation model P _{total, PAp}minimum (I) and the highest (u) border, iterations N, annealing factor q are set, (q > 0), error margins ε _f;

Step 32, to primary vector and secondary vector assignment, is specially α=α ₀, α ⁰=α, obtaining optimum solution can be expressed as follows:

P ⁰＝f(α ⁰)

Wherein, α ⁰optimum solution, P ⁰to take under the prerequisite of optimal distributing scheme, the minimum value obtaining, f is SA function, f function is above-mentioned P _{total, PAp}calculation relational expression.

Step 33, iteration.From n=0, start until n=N.

First, producing at random one has the vectorial β of equal dimension, β ∈ U[-1,1 with α], then calculate Δ α, its expression formula is as follows:

$\mathrm{\&Delta;\&alpha;}=(u-1)*g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)$

Wherein, $g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)=\frac{{(1+\mathrm{\&mu;})}^{\left|\mathrm{\&beta;}\right|}-1}{\mathrm{\&mu;}}\mathrm{sign}\left(\mathrm{\&beta;}\right),$ μ=10 ¹⁰⁰(n/N) ^q, β is the vector that random of producing and α have equal dimension, u is model P _{total, PAp}gao border, n is current iteration number of times, N is maximum iteration time, q is annealing factor.

Then a point is set in ， interval (I, u), specifically can represents as follows: α '=α+Δ α

If Δ P=P (α ')-P (α) < 0 or z < p=exp are ((n/N) ^qΔ P/|P (α) |/ε _f) (z is in interval [0,1] the random number producing), α=α ' and P (α)=P (α '); Wherein, P (α ') and P (α) are according to P _{total, PAp}the value calculated of calculation relational expression.

If P (α) < is P ⁰, α ⁰=α and P ⁰=P (α ⁰).

Finally, when iteration finishes, α ⁰be exactly the optimal solution of P (α), P ⁰it is exactly resulting minimum value.

As shown in Figure 4, the distributor for power in a kind of Turbo Detection for Cooperative Communication shown in the present, comprising:

Acquiring unit 41, obtains the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

Computing unit 42, according to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node.

Allocation units 43, are used described assignment information to carry out power division.

The channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

Described computing unit 42 comprises:

First sets subelement, sets total transmitted power P _{total, PAp}minimum border I He Gao border u;

Second sets subelement, sets the primary vector α of K dimension and the secondary vector α of K dimension ⁰,

The first numerical value produces subelement, produces at random the initial value α of primary vector α ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

The first computation subunit, according to described the first expression formula, calculates corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Value subelement, gets current iteration frequency n=1;

Second value produces subelement, while carrying out iteration, produces at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

The second computation subunit, calculates the first variable Δ α according to the 3rd vectorial β, is specially: $\mathrm{\&Delta;\&alpha;}=(u-1)*g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right);$ Wherein, $g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)=\frac{{(1+\mathrm{\&mu;})}^{\left|\mathrm{\&beta;}\right|}-1}{\mathrm{\&mu;}}\mathrm{sign}\left(\mathrm{\&beta;}\right);$ Sign (β) is sign function; The second variable μ=10 ¹⁰⁰(n/N) ^q;

The 3rd computation subunit, generates the four-way amount α ' in interval (I, u) according to the first value of variable Δ α and the currency of primary vector α, is specially: α '=α+Δ α;

The 4th computation subunit, according to described the first relational expression, calculates the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

First upgrades subelement, when Δ P=P (α ')-P (α) < 0 or z < p=exp ((n/N) ^qΔ P/|P (α) |/ε _f) time, use the currency of the 3rd vectorial α ' to upgrade primary vector α; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0; Described maximum iteration time N can equate with the total quantity K of described via node.

Second upgrades subelement, as P (α) < P (α ⁰) time, use the currency of primary vector α to upgrade secondary vector α ⁰;

Judgment sub-unit, judges whether the currency of current iteration frequency n equals maximum iteration time N;

The 5th operator unit, when the currency of current iteration frequency n equals maximum iteration time N, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node; If NO, current iteration frequency n adds 1, proceeds iteration.

Described the first computation subunit comprises:

Measurement module, by measuring, obtains power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

Acquisition module, according to modulation system, obtains the value of modulation system constant m;

Module is set, according to target error rate, error rate of system P is set _evalue;

Computing module, according to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

The above is the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, do not departing under the prerequisite of principle of the present invention; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (10)

1. a distribution method for power in Turbo Detection for Cooperative Communication, is characterized in that, comprising:

Obtain the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

According to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node;

Use described assignment information to carry out power division.

2. the distribution method of power in Turbo Detection for Cooperative Communication according to claim 1, is characterized in that,

The channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

3. the distribution method of transmitted power according to claim 2, it is characterized in that, described according to described the first relational expression, when calculating described total transmitted power and being minimum value described in the step of ratio between the transmitted power of each via node and the transmitted power of described source node comprise:

Set total transmitted power P _{total, PAp}minimum border I He Gao border u, the primary vector α of a K dimension and the secondary vector α of K dimension ⁰;

The random initial value α that produces primary vector α ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

According to described the first expression formula, calculate corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Get current iteration frequency n=1;

Carry out iteration, produce at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

According to the 3rd vectorial β, calculate the first variable Δ α;

According to the first value of variable Δ α and the currency of primary vector α, generate the four-way amount α ' in interval (I, u);

According to described the first relational expression, calculate the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

If Δ P=P (α ')-P (α) < 0 or z < p=exp are ((n/N) ^qΔ P/|P (α) |/ε _f), use the currency of the 3rd vectorial α ' to upgrade primary vector α; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0;

If P (α) < P is (α ⁰), use the currency of primary vector α to upgrade secondary vector α ⁰;

Whether the currency that judges current iteration frequency n equals maximum iteration time N;

If YES, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node; If NO, current iteration frequency n adds 1, proceeds iteration.

4. the distribution method of power in Turbo Detection for Cooperative Communication according to claim 3, is characterized in that,

The described step of calculating the first variable Δ α according to the 3rd vectorial β is specially:

$\mathrm{\&Delta;\&alpha;}=(u-1)*g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right);$ Wherein, $g_{\mathrm{\&mu;}}^{-1}\left(\mathrm{\&beta;}\right)=\frac{{(1+\mathrm{\&mu;})}^{\left|\mathrm{\&beta;}\right|}-1}{\mathrm{\&mu;}}\sin \mathrm{gn}\left(\mathrm{\&beta;}\right);$ Sign (β) is sign function; The second variable μ=10 ¹⁰⁰(n/N) ^q; Or

The described step according to the four-way amount α ' in the first value of variable Δ α and the currency of primary vector α generation interval (I, u) is specially: α '=α+Δ α.

5. the distribution method of power in Turbo Detection for Cooperative Communication according to claim 3, is characterized in that,

Described maximum iteration time N equates with the total quantity K of described via node.

6. the distribution method of power in Turbo Detection for Cooperative Communication according to claim 3, is characterized in that, described according to described the first expression formula, calculates corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰) step comprise:

By measuring, obtain power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

According to modulation system, obtain the value of modulation system constant m;

According to target error rate, error rate of system P is set _evalue;

According to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

7. a distributor for power in Turbo Detection for Cooperative Communication, is characterized in that, comprising:

Acquiring unit, obtains the first relational expression between the channel state parameter three of ratio between the transmitted power of total transmitted power of source node and each via node, described each via node and the transmitted power of described source node and system;

Computing unit, according to described the first relational expression, ratio described in when calculating described total transmitted power and being minimum value between the transmitted power of each via node and the transmitted power of described source node, as the assignment information of the transmitted power of described each via node and the transmitted power of described source node;

Allocation units, are used described assignment information to carry out power division.

8. the distributor of power in Turbo Detection for Cooperative Communication according to claim 7, is characterized in that,

The channel state parameter of described system comprises: the power spectral density N of noise in modulation system constant m, channel ₀, error rate of system P _e, the channel gain between source node and destination node the first mean value F, source node and k via node between the second mean value G of channel gain _k, the channel gain between a k via node and destination node the 3rd mean value H _k;

Described the first relational expression is:

$P_{\mathrm{total},\mathrm{PAp}}=(1+\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_k){[\frac{3}{4m^2}\&CenterDot;\frac{N_0}{P_{e}F}\underset{k=1}{\overset{K}{\mathrm{\&Pi;}}}(\frac{1}{G_k}+\frac{1}{{\mathrm{\&alpha;}}_k{H}_k})]}^{\frac{1}{K+1}};$

Wherein, P _{tOTAL, RAP}for described total transmitted power;

K is the total quantity of via node, and K is greater than 0 natural number;

K is the sequence number of via node, k=1 ..., K;

α _kthe ratio between the transmitted power of k via node and the transmitted power of source node.

9. the distributor of transmitted power according to claim 8, is characterized in that, described computing unit comprises:

First sets subelement, sets total transmitted power P _{total, PAp}minimum border I He Gao border u;

Second sets subelement, sets the primary vector α of K dimension and the secondary vector α of K dimension ⁰,

The first numerical value produces subelement, produces at random the initial value α of primary vector α ₀, and by initial value α ₀assignment is given described secondary vector α ⁰;

The first computation subunit, according to described the first expression formula, calculates corresponding first total transmission power value P (α), the described secondary vector α of currency of described primary vector α ⁰second total transmission power value P (α corresponding to currency ⁰);

Value subelement, gets current iteration frequency n=1;

Second value produces subelement, while carrying out iteration, produces at random the 3rd vectorial β of K dimension, β ∈ U[-1,1];

The second computation subunit, calculates the first variable Δ α according to the 3rd vectorial β;

The 3rd computation subunit, generates the four-way amount α ' in interval (I, u) according to the first value of variable Δ α and the currency of primary vector α, is specially:;

The 4th computation subunit, according to described the first relational expression, calculates the 4th total transmission power value P (α ') corresponding to currency of four-way amount α ';

First upgrades subelement, when Δ P=P (α ')-P (α) < 0 or z < p=exp ((n/N) ^qΔ P/|P (α) |/ε _f) time, use the currency of the 3rd vectorial α ' to upgrade primary vector α; Wherein, ternary z is the random number producing in interval [0,1]; Exp function is exponential function; N is maximum iteration time, ε _ffor error margins, q is annealing factor, q>0;

Second upgrades subelement, as P (α) < P (α ⁰) time, use the currency of primary vector α to upgrade secondary vector α ⁰;

Judgment sub-unit, judges whether the currency of current iteration frequency n equals maximum iteration time N;

The 5th operator unit, when the currency of current iteration frequency n equals maximum iteration time N, iteration finishes, by secondary vector α ⁰currency while being minimum value as described total transmitted power described in ratio between the transmitted power of each via node and the transmitted power of described source node; If NO, current iteration frequency n adds 1, proceeds iteration.

10. the distributor of power in Turbo Detection for Cooperative Communication according to claim 9, is characterized in that, described the first computation subunit comprises:

Measurement module, by measuring, obtains power spectral density N ₀, the first mean value F, the second mean value G _kand the 3rd mean value H _kvalue;

Acquisition module, according to modulation system, obtains the value of modulation system constant m;

Module is set, according to target error rate, error rate of system P is set _evalue;

Computing module, according to described power spectral density N ₀, the first mean value F, the second mean value G _k, the 3rd mean value H _k, modulation system constant m and error rate of system P _evalue, by described the first expression formula, calculate respectively corresponding the first total power value P (α), the described secondary vector α of currency of described primary vector α ⁰the second total power value P (α corresponding to currency ⁰).

Images