CN109890039B - Method for allocating SWIPT relay resources under general interference - Google Patents

Method for allocating SWIPT relay resources under general interference Download PDF

Info

Publication number
CN109890039B
CN109890039B CN201910092634.7A CN201910092634A CN109890039B CN 109890039 B CN109890039 B CN 109890039B CN 201910092634 A CN201910092634 A CN 201910092634A CN 109890039 B CN109890039 B CN 109890039B
Authority
CN
China
Prior art keywords
relay
signal
probability
interruption
receiver
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201910092634.7A
Other languages
Chinese (zh)
Other versions
CN109890039A (en
Inventor
李建雄
赵可
丁学龙
史伟光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Polytechnic University
Original Assignee
Tianjin Polytechnic University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin Polytechnic University filed Critical Tianjin Polytechnic University
Priority to CN201910092634.7A priority Critical patent/CN109890039B/en
Publication of CN109890039A publication Critical patent/CN109890039A/en
Application granted granted Critical
Publication of CN109890039B publication Critical patent/CN109890039B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radio Relay Systems (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses a method for allocating SWIPT relay resources under general interference, which is based on a DF relay system, establishes an optimization problem by taking the maximum non-interruption probability as a target, and solves the problem by adopting a golden section method to obtain the optimal h2And maximum probability of non-interruption. Simulation experiments show that the method can effectively improve the non-interruption probability of the relay system under different parameters.

Description

Method for allocating SWIPT relay resources under general interference
Technical Field
The invention relates to the technical field of relay resource allocation, in particular to an SWIPT relay resource allocation method under general interference.
Background
With the advent of the 5G era and the expansion of the scale of wireless networks, the number of nodes in the network is increasing dramatically. Conventional network nodes use batteries to supply energy, which means that the network nodes have a limited life cycle. Faced with more and more nodes, replacing batteries for them obviously consumes a lot of manpower and material resources. However, the emergence and development of the Simultaneous transfer (SWIPT) technology provides a new idea for solving the problem. The SWIPT technology can enable nodes in a wireless network to collect energy through radio frequency signals, and receive information by using the collected energy. The technology avoids the disadvantage of frequent battery replacement, and can greatly prolong the service life of nodes in the network. Therefore, the application of the SWIPT technology in the wireless network is of great interest.
The SWIPT technology can effectively improve the frequency spectrum utilization rate of the network, reduce delay and reduce power consumption, so that a lot of students consider applying the SWIPT technology to a relay communication system. And when the SWIPT is applied to the relay system, a reasonable resource allocation strategy needs to be adopted to improve the performance of the relay system. The existing method for allocating SWIPT relay resources generally only considers the situation that interference is periodic information, however, in practical application, the interference often appears not periodically but randomly.
Therefore, how to provide a method for allocating SWIPT relay resources for general interference is a problem that needs to be solved urgently by those skilled in the art.
Disclosure of Invention
In view of this, the present invention provides a method for allocating SWIPT relay resources under general interference, which can effectively improve the non-interruption probability of the system.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for allocating SWIPT relay resources under general interference comprises the following steps:
(1) establishing a relay system model based on DF;
(2) defining the non-interruption probability p of the whole link based on the model established in the step (1)non＝p1p2(ii) a Wherein the content of the first and second substances,
Figure GDA0002043211370000021
is the probability of non-interruption of the transmitter to the receiver;
Figure GDA0002043211370000022
a non-outage probability for a relay to a receiver; wherein, PsFor the average power of the transmitted signal, PRTo forward the power of the signal, r0For the transmission rate threshold, ε is a straight line h2Point of intersection with the h-axis, h1The DF region and the EH region are the dividing lines,
Figure GDA0002043211370000023
h2line for resource distribution
Figure GDA0002043211370000024
I is the power value I, theta of the interference IIIs the expectation of I, thetahTheta is the expectation of the channel gain hgFor a channelExpectation of gain g, and channel gains h and g both obey an exponential;
(3) the optimization problem P5 is established with the goal of maximizing the probability of non-interruption:
(P5)max pnon(ε)＝p1p2
Figure GDA0002043211370000025
(4) solving the optimization problem P5 in the step (3) by adopting a golden section method to obtain an optimal resource allocation section line h2And maximum probability of non-interruption.
In order to further optimize the above technical solution, the step (4) specifically includes:
step 1: setting a value interval [ a, b ] and precision e of the initialization epsilon;
step 2: solving section golden section point a1＝a+(1-0.618)(b-a)，a2＝a+0.618(b-a)；
And 3, step 3: if the probability of non-interruption pnon(a1)＜pnon(a2) Jumping to the 4 th step, otherwise jumping to the 5 th step;
and 4, step 4: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a1Outputting the maximum non-interrupt probability
Figure GDA0002043211370000026
Otherwise, let a be a1，a1＝a2，a2Jumping to step 3 as a +0.618 (b-a);
and 5, step 5: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a2Outputting the maximum non-interrupt probability
Figure GDA0002043211370000031
Otherwise, let b be a2，a2＝a1，a1A + (1-0.618) (b-a), jump to step 3.
The technical proposal is similar to the prior artCompared with the prior art, the invention discloses and provides a method for allocating SWIPT relay resources under general interference, a DF-based relay system establishes an optimization problem by taking the maximum non-interruption probability as a target, and the optimal h is calculated by adopting a golden section method2And maximum probability of non-interruption. Simulation experiments show that the method can effectively improve the non-interruption probability of the relay system under different parameters.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
Fig. 1 is a schematic structural diagram of a method for allocating SWIPT relay resources under general interference provided by the present invention;
fig. 2 is a schematic diagram of a DF-based relay system model provided in the present invention;
FIG. 3 is a schematic diagram of relay timeslot division provided by the present invention;
FIG. 4 is a h-I resource allocation coordinate system provided by the present invention;
FIG. 5 shows a graph h according to the present invention2A schematic diagram of method resource allocation;
FIG. 6 shows the non-interruption probability of the relay system using the TS method according to the present invention;
FIG. 7 shows the present invention using h2The method relays the probability of system non-interruption;
FIG. 8 shows different interrupt rates r provided by the present invention0Maximum non-interruption probability of the lower relay system;
FIG. 9 shows different power consumptions P of the circuit according to the present inventioncMaximum non-interruption probability of the lower relay system;
FIG. 10 shows different interference parameters θ provided by the present inventionIMaximum non-interruption probability of the lower relay system.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, an embodiment of the present invention discloses a method for allocating SWIPT relay resources under general interference, including:
(1) establishing a relay system model based on DF;
(2) defining the non-interruption probability p of the whole link based on the model established in the step (1)non＝p1p2(ii) a Wherein the content of the first and second substances,
Figure GDA0002043211370000041
is the probability of non-interruption of the transmitter to the receiver;
Figure GDA0002043211370000042
a non-outage probability for a relay to a receiver; wherein, PsFor the average power of the transmitted signal, PRTo forward the power of the signal, r0For the transmission rate threshold, ε is a straight line h2Point of intersection with the h-axis, h1The DF region and the EH region are the dividing lines,
Figure GDA0002043211370000043
h2line for resource distribution
Figure GDA0002043211370000044
I is the power value I, theta of the interference IIIs the expectation of I, thetahTheta is the expectation of the channel gain hgIs the expectation of channel gain g, and both channel gains h and g obey an index;
(3) the optimization problem P5 is established with the goal of maximizing the probability of non-interruption:
(P5)max pnon＝p1p2
Figure GDA0002043211370000045
(4) solving the optimization problem P5 in the step (3) by adopting a golden section method to obtain an optimal resource allocation section line h2And maximum probability of non-interruption.
In order to further disclose the technical scheme, the step (4) specifically comprises the following steps:
the step (4) specifically comprises the following steps:
step 1: setting a value interval [ a, b ] and precision e of the initialization epsilon;
step 2: solving section golden section point a1＝a+(1-0.618)(b-a)，a2＝a+0.618(b-a)；
And 3, step 3: if the probability of non-interruption pnon(a1)＜pnon(a2) Jumping to the 4 th step, otherwise jumping to the 5 th step;
and 4, step 4: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a1Outputting the maximum non-interrupt probability
Figure GDA0002043211370000051
Otherwise, let a be a1，a1＝a2，a2Jumping to step 3 as a +0.618 (b-a);
and 5, step 5: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a2Outputting the maximum non-interrupt probability
Figure GDA0002043211370000052
Otherwise, let b be a2，a2＝a1，a1A + (1-0.618) (b-a), jump to step 3.
Under the condition that the existing interference is periodic information, the relay mostly adopts an AF forwarding mode, and the relay has the advantages of simple system structure and low circuit consumption energy. However, under general interference, the AF forwarding scheme does not perform well. Since the interference exists all the time and occurs non-periodically in a general interference environment, the relay amplifies the signal and also amplifies the interference signal and noise, which obviously reduces the signal to noise ratio of the signal and is not beneficial to the transmission of system information. The DF has the advantage that the relay decodes the received information before forwarding the signal, which allows the interference and noise in the original signal to be filtered out and avoids further transmission of the interference and noise. Therefore, in a general interference environment, the relay is more suitable to adopt the DF forwarding method.
A relay system model using the DF forwarding scheme is shown in fig. 2. It should be noted that, unlike the AF forwarding method, the DF forwarding method needs to decode information and consumes a certain amount of energy. The energy collected by the system is divided into two parts, one part is used for forwarding information, and the power of the forwarded signal is assumed to be PR(ii) a Another part is consumed by the decoding circuit, assuming that the power consumed is PC。
In FIG. 2, the transmitter is for transmitting signals
Figure GDA0002043211370000053
Denotes that its average power is PS. The signal reaches the relay through a wireless channel, the channel gain is h, and meanwhile, a general interference signal i also reaches the relay. y isRIndicating relayed received signal, nRIs the relay antenna noise. The relay received signal is divided into two paths by the time switcher, the upper path signal is used for decoding and forwarding, the lower path signal is used for collecting energy and providing the size P for the forwarded signalRThe signal after the forward amplification power is amplified is xR. The signal arrives at the receiver via a radio channel with a channel gain g, the received signal being denoted yD. At the same time, receiver antenna noise
Figure GDA0002043211370000054
Is introduced.
In a general interference environment, the signal-to-noise ratio of signals received by the relay and the receiver is affected by the channel state, and when the channel state is good, the signal-to-noise ratio of the received signals is high, and conversely, the signal-to-noise ratio of the signals is poor. And with interferenceThe signal-to-noise ratio of the signal fluctuates with time. Considering the channel gains h, g as mutually independent random variables and respectively subject to the expectation of theta in conjunction with fig. 2hAnd thetagIs used as the index distribution of (1). The relay adopts a TS operation strategy, and at the moment, the relay has two working modes which are respectively defined as a Decoding Forwarding (DF) mode and an energy collecting (EH) mode. The relay can judge the optimal working state of the relay according to h, I and g, wherein I is an interference power value. As shown in fig. 3, the communication process is divided into several slots and it is assumed that the relay can acquire channel state information, i.e., h, I and g. When a time slot arrives, the relay determines to allocate the time slot to EH or DF by judging the channel state information.
Considering the communication process from the transmitter to the relay, referring to fig. 2, in a certain channel state, the received signal at the relay is represented as:
Figure GDA0002043211370000061
the signal to noise ratio is:
Figure GDA0002043211370000062
at this time, the instantaneous information rate can be expressed as
Figure GDA0002043211370000063
The relay energy collects power of
PEH＝η(hPS+I)
Where η is energy conversion efficiency, and for convenience, η is 1.
Transmission being limited by taking into account delay, i.e. transmission rate below threshold r0The communication interruption is caused by time, and the probability of interruption of the communication from the transmitter to the relay is expressed as:
Figure GDA0002043211370000064
defining the non-outage probability p1＝1-pout，p1Then this can be expressed as:
Figure GDA0002043211370000071
considering the communication process from the relay to the receiver in fig. 2, in a certain channel state, the receiver receives signals as follows:
Figure GDA0002043211370000072
in the formula, xRThe forwarded signal representing the relay, assuming the relay can decode the information correctly, then
Figure GDA0002043211370000073
The signal-to-noise ratio of the received signal at the receiver is therefore:
Figure GDA0002043211370000074
at this time, the instantaneous information rate relayed to the receiver can be expressed as:
Figure GDA0002043211370000075
also consider a threshold value of r0Delay limited transmission, non-interruption probability p of relaying to receiver2Can be expressed as follows:
Figure GDA0002043211370000076
with h as the vertical axis and I as the horizontal axis, an h-I coordinate system, referred to as a resource allocation coordinate system, can be obtained, as shown in fig. 4. If it is required that the transmitter-to-relay transmission cannot be interrupted, then h and I must satisfy:
Figure GDA0002043211370000077
the method is simplified and can be obtained:
Figure GDA0002043211370000078
order to
Figure GDA0002043211370000079
The condition that no interruption occurs in the transmitter-to-relay transmission is mapped to the h-I coordinate system, so that the non-interruption region is h in FIG. 41The hatched portion above the straight line. In this area, no interruption occurs in the transmission of information, and the relay should operate in the DF mode, so this area is called the DF area. In contrast, in h1The lower zone for information transmission causes an interruption, and therefore, h1The lower zone relay should operate in the EH mode, and this zone is referred to as an EH zone. The energy collected by the relay is divided into two parts, and can be expressed by the formula:
PEH＝PC+PR
PEHfor collecting power, PCA fixed power consumed for the circuit, therefore, PRDepending on the amount of relay energy harvesting power. For relaying to the receiver, the uninterrupted communication needs to satisfy RRD≥r0I.e. by
Figure GDA0002043211370000081
For different channel states g, the relays require different forwarding powers to ensure uninterrupted communication. Therefore, in fig. 4, a certain area needs to be partitioned from the decode-and-forward area for energy collection. Defining resource allocation lines h2Is a straight line perpendicular to the h-axis, and the intersection with the h-axis is (0,. epsilon.) as shown in FIG. 5. h is2And h1The common upper region is divided for energy harvesting, h-axis,h1And h2The jointly formed triangular regions are used for decoding and forwarding. h is1The intersection point with the h axis is
Figure GDA0002043211370000082
Figure GDA0002043211370000083
The transmitter always transmits information to the relay, so the shaded triangular area cannot be 0, and the resource allocation line needs to be satisfied
Figure GDA0002043211370000084
I.e. the value range of epsilon is
Figure GDA0002043211370000085
For the DF and EH region allocation case of fig. 5, the transmitter-to-receiver non-outage probability can be expressed as:
Figure GDA0002043211370000086
the above equation is expressed as an integral:
Figure GDA0002043211370000087
the average power of the relay energy collection is:
Figure GDA0002043211370000088
after the average power of the relay energy collection is obtained, according to PEH＝PC+PRThe average forwarding power of the relay can be obtained as follows:
E[PR]＝E[PEH]-PC
the probability of non-interruption relayed to the receiver can be expressed as an integral:
Figure GDA0002043211370000091
defining the non-interruption probability of the whole link as pnonThe communication of the whole link is not interrupted, and it must be satisfied that the communication from the transmitter to the relay is not interrupted and the communication from the relay to the receiver is not interrupted at the same time, so the probability of non-interruption of the whole link can be expressed as:
pnon＝p1p2
at this time, the optimization problem P5 is established with the goal of maximizing the non-interruption probability:
(P5)max pnon＝p1p2
Figure GDA0002043211370000092
the optimization problem can be solved by adopting a golden section method.
The technical scheme provided by the invention is further explained by combining specific experimental simulation.
For convenience of subsequent description, the method for allocating SWIPT relay resources under general interference provided by the invention is simply referred to as h2A method.
Invention pair h2The method is subjected to simulation verification and is compared with the method for optimizing the TS coefficient. The simulation parameter settings used are as follows, PS＝5，PC＝0.5，r0＝0.2，θh＝1，θI＝5，θg＝1，
Figure GDA0002043211370000093
Figure GDA0002043211370000094
Referring to fig. 6, fig. 6 is a non-interruption probability of a relay system by using a TS coefficient method, which can be seen visually: when the TS coefficient is increased, the probability of system non-interruption is increased firstly and then reduced, and when alpha is 0.22, the rate of system non-interruption reachesMaximum, pnon＝0.59。
FIG. 7 shows relay adoption h2In the method, the non-interrupt probability of the system is distributed along with the resource2Simulation diagram of the change of (1). It can be seen from the figure that the probability of non-interruption p from the transmitter to the relay1Increases gradually with increasing epsilon because increasing epsilon increases the DF area in FIG. 5, while p increases1Is integrated over this area, and thus, p1It will increase. But when epsilon increases to a certain value, the probability of non-interruption p1Remains substantially unchanged because the integral of the probability density function approaches 0 in the DF integration region expanded by the increase in epsilon. Probability of non-interruption p relayed to receiver2Decreasing with increasing epsilon, because the EH area in fig. 5 decreases, the energy received by the repeater decreases, resulting in a decrease in the repeating power, and therefore a decrease in the signal-to-noise ratio of the repeating signal and a decrease in the probability of non-interruption. The best dividing line h can be obtained by golden section method26.43, the system non-outage probability is pnon0.68. It can be seen that employing h2The method can improve the non-interruption probability of the system by 15.3 percent.
Under different parameters, h2The process exhibits different properties. As shown in FIG. 8, the interrupt rate r is adjusted by fixing other parameters of the system0. With r0Is increased by h2The method and the relay non-interruption probability adopting the TS coefficient method are both gradually reduced, because the interruption rate is improved, the signal-to-noise ratio requirement of the system on the signal is improved, therefore, under the condition that other conditions are not changed, the signal with lower signal-to-noise ratio can generate communication interruption, and the non-interruption probability is reduced. However, it is apparent that h is used2The relay non-interruption probability of the method is larger than h which is not adopted2The relay non-interruption probability of the method is low, and the difference between the two is along with r0Is increased. Therefore, the relay system with high requirement on delay transmission is more suitable for adopting h2A method.
Fixing other parameters of the system, regulating the power consumption P of the circuitCTo obtain different non-interruption probabilities, as shown in FIG. 9. With PCIs increased by h2The method and the relay non-interruption probability adopting the TS coefficient method are both gradually reduced because the circuit consumption power is increased, on one hand, the relay needs to allocate more resources for energy collection, thereby reducing the area of the DF area in fig. 5. On the other hand, the power P for the retransmitted signalRThe signal-to-noise ratio of the retransmitted signal is reduced. This therefore leads to a reduction in the probability of relay non-interruption, otherwise unchanged. However, as can be seen from FIG. 9, h is used2The relay non-interruption probability of the method is obviously greater than that of the TS coefficient method.
As shown in fig. 10, the other parameters of the system are fixed, and the disturbance parameter theta is adjustedI. With thetaIIs increased by h2The method and the relay non-interruption probability adopting the TS coefficient method are both gradually reduced, because the signal-to-noise ratio of the relay receiving signal of the system is reduced due to the improvement of the interference power, and therefore, the non-interruption probability is reduced under the condition that other conditions are not changed. However, as is apparent from FIG. 10, h is used2The relay non-interruption probability of the method is greater than that of the TS coefficient method, and h is adopted2The non-outage probability of the method drops very slowly. This is because, when the interference is strong or the channel state is poor, the relay switches to the energy receiving mode; when the interference is weak or the channel state is good, the relay is switched to a decoding forwarding mode. However, for the relay adopting the TS coefficient method, regardless of the interference magnitude and the channel state, the relay can only perform the conversion of energy collection and decoding forwarding according to the fixed time switching coefficient, so for some states with large interference power or poor channel state, the relay system adopting the TS coefficient method may be in a decoding forwarding mode, which causes the signal-to-noise ratio of the relay received signal or the signal received by the receiver to be very low, and finally causes the non-interruption probability of the relay to be reduced.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (2)

1. A method for allocating SWIPT relay resources under general interference is characterized by comprising the following steps:
(1) establishing a relay system model based on DF;
(2) defining the non-interruption probability p of the whole link based on the model established in the step (1)non＝p1p2；
Wherein the content of the first and second substances,
Figure FDA0002976459090000011
is the probability of non-interruption of the transmitter to the receiver;
Figure FDA0002976459090000012
a non-outage probability for a relay to a receiver; wherein, PsFor the average power of the transmitted signal, PRTo forward the power of the signal, r0For the transmission rate threshold, ε is a straight line h2Point of intersection with the h-axis, h1The DF region and the EH region are the dividing lines,
Figure FDA0002976459090000013
h2line for resource distribution
Figure FDA0002976459090000014
I is the power value of interference II，θIIs the expectation of I, thetahTheta is the expectation of the channel gain hgIs the expectation of channel gain g, and both channel gains h and g obey an index;
(3) the optimization problem P5 is established with the goal of maximizing the probability of non-interruption:
(P5)max pnon(ε)＝p1p2
Figure FDA0002976459090000015
(4) solving the optimization problem P5 in the step (3) by adopting a golden section method to obtain an optimal resource allocation section line h2And a maximum probability of non-interruption;
the DF-based relay system model adopts a DF forwarding mode; the collected energy is divided into two parts, one part is used for forwarding information, and the power of the forwarded signal is assumed to be PR(ii) a Another part is consumed by the decoding circuit, assuming that the power consumed is PC；
For transmitting signals from transmitters
Figure FDA0002976459090000016
Denotes that its average power is PS(ii) a The signal reaches the relay through a wireless channel, the channel gain is h, and a general interference signal i also reaches the relay; y isRIndicating relayed received signal, nRIs the relay antenna noise; the relay received signal is divided into two paths by the time switcher, the upper path signal is used for decoding and forwarding, the lower path signal is used for collecting energy and providing the size P for the forwarded signalRThe signal after the forward amplification power is amplified is xR(ii) a The signal arrives at the receiver via a radio channel with a channel gain g, the received signal being denoted yD(ii) a At the same time, receiver antenna noise
Figure FDA0002976459090000021
Is introduced;
in an interference environment, the signal-to-noise ratio of signals received by a relay and a receiver is influenced by a channel state, when the channel state is good, the signal-to-noise ratio of the received signals is high, otherwise, the signal-to-noise ratio of the signals is poor; moreover, with the continuous change of the interference, the signal-to-noise ratio of the signal also fluctuates; considering the channel gains h, g as mutually independent random variables and subject to the expectation of θhAnd thetagThe distribution of indices; the relay adopts a TS operation strategy, and at the moment, the relay has two working modes which are respectively defined as a Decoding Forwarding (DF) mode and an energy collecting (EH) mode; the relay judges the optimal working state of the relay according to h, I and g, wherein I is an interference power value; dividing the communication process into a plurality of time slots, and assuming that a relay acquires channel state information, namely h, I and g; when a time slot arrives, the relay determines to distribute the time slot to EH or DF by judging the channel state information;
considering the communication process from the transmitter to the relay, in a certain channel state, the received signal of the relay is represented as:
Figure FDA0002976459090000022
the signal to noise ratio is:
Figure FDA0002976459090000023
at this time, the instantaneous information rate is expressed as
Figure FDA0002976459090000024
The relay energy collects power of
PEH＝η(hPS+I)
Wherein η is energy conversion efficiency, and for convenience, η is 1;
transmission with delay constraints taken into account, i.e. transmission rates lower thanThreshold value r0The communication interruption is caused by time, and the probability of interruption of the communication from the transmitter to the relay is expressed as:
Figure FDA0002976459090000031
defining the non-outage probability p1＝1-pout，p1Then the expression is:
Figure FDA0002976459090000032
considering the communication process from the relay to the receiver, in a certain channel state, the receiver receives the signal:
Figure FDA0002976459090000033
in the formula, xRThe forwarded signal representing the relay, assuming the relay correctly decoded the information, then
Figure FDA0002976459090000034
The signal-to-noise ratio of the received signal at the receiver is therefore:
Figure FDA0002976459090000035
at this point, the instantaneous information rate relayed to the receiver is expressed as:
Figure FDA0002976459090000036
also consider a threshold value of r0Delay limited transmission, non-interruption probability p of relaying to receiver2Is represented as follows:
Figure FDA0002976459090000037
taking h as a vertical axis and I as a horizontal axis, obtaining an h-I coordinate system, called a resource allocation coordinate system, and if the transmitter-to-relay transmission is required to be uninterrupted, h and I must satisfy:
Figure FDA0002976459090000041
simplifying to obtain:
Figure FDA0002976459090000042
order to
Figure FDA0002976459090000043
Mapping the condition that the transmission from a transmitter to a relay is not interrupted to an h-I coordinate system, wherein the information transmission in a non-interruption area is not interrupted, and the relay is supposed to work in a DF mode, so that the area is called as a DF area; in contrast, in h1The lower zone for information transmission causes an interruption, and therefore, h1The lower zone relay should operate in EH mode, and this zone is called EH zone; the energy collected by the relay is divided into two parts, and is expressed by the formula:
PEH＝PC+PR
PEHfor collecting power, PCA fixed power consumed for the circuit, therefore, PRThe magnitude of (d) depends on the magnitude of the relay energy harvesting power; for relaying to the receiver, the uninterrupted communication needs to satisfy RRD≥r0I.e. by
Figure FDA0002976459090000044
For different channel states g, the relay needs different forwarding powers to ensure uninterrupted communication; it is necessary to partition a certain area from a decode-and-forward area for energy collection(ii) a Defining resource allocation lines h2Is a straight line perpendicular to the h-axis, and the intersection point of the straight line and the h-axis is (0, epsilon); h is2And h1The common upper region is divided for energy harvesting, h-axis, h1And h2The triangular areas which are formed together are used for decoding and forwarding; h is1The intersection point with the h axis is
Figure FDA0002976459090000045
The transmitter always transmits information to the relay, so the shaded triangular area cannot be 0, and the resource allocation line needs to be satisfied
Figure FDA0002976459090000046
I.e. the value range of epsilon is
Figure FDA0002976459090000047
For the DF and EH region allocation cases, the transmitter-to-receiver non-outage probability is expressed as:
Figure FDA0002976459090000048
2. the method for allocating SWIPT relay resources under general interference according to claim 1, wherein the step (4) specifically comprises:
step 1: setting a value interval [ a, b ] and precision e of the initialization epsilon;
step 2: solving section golden section point a1＝a+(1-0.618)(b-a)，a2＝a+0.618(b-a)；
And 3, step 3: if the probability of non-interruption pnon(a1)＜pnon(a2) Jumping to the 4 th step, otherwise jumping to the 5 th step;
and 4, step 4: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a1Outputting the maximum non-interrupt probability
Figure FDA0002976459090000051
Otherwise, let a be a1，a1＝a2，a2Jumping to step 3 as a +0.618 (b-a);
and 5, step 5: if a is2-a1E, stopping iteration and outputting an optimal solution x*＝a2Outputting the maximum non-interrupt probability
Figure FDA0002976459090000052
Otherwise, let b be a2，a2＝a1，a1A + (1-0.618) (b-a), jump to step 3.
CN201910092634.7A 2019-01-30 2019-01-30 Method for allocating SWIPT relay resources under general interference Expired - Fee Related CN109890039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910092634.7A CN109890039B (en) 2019-01-30 2019-01-30 Method for allocating SWIPT relay resources under general interference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910092634.7A CN109890039B (en) 2019-01-30 2019-01-30 Method for allocating SWIPT relay resources under general interference

Publications (2)

Publication Number Publication Date
CN109890039A CN109890039A (en) 2019-06-14
CN109890039B true CN109890039B (en) 2021-06-25

Family

ID=66927454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910092634.7A Expired - Fee Related CN109890039B (en) 2019-01-30 2019-01-30 Method for allocating SWIPT relay resources under general interference

Country Status (1)

Country Link
CN (1) CN109890039B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111225399B (en) * 2020-02-26 2021-03-05 电子科技大学 Relay forwarding and resource allocation method in wireless data energy simultaneous transmission cooperative communication

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104301977A (en) * 2014-09-22 2015-01-21 西安交通大学 Transmission method in SWIPT system under interference channel
CN105610485A (en) * 2015-12-21 2016-05-25 东南大学 Wireless relay communication system SWIPT (Simultaneous Wireless Information and Power Transfer) method
KR20170120811A (en) * 2016-04-22 2017-11-01 성균관대학교산학협력단 Full-duplex relay communication method based on energy harvesting and full-duplex relay communication system based on energy harvesting
CN108039942A (en) * 2017-12-11 2018-05-15 天津工业大学 Passive relay collects interfering energy to improve optimal system Speed method in SWIPT systems
CN109088686A (en) * 2018-09-28 2018-12-25 西北工业大学 One kind being based on wireless messages and energy transmission method while 5G height frequency range

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104301977A (en) * 2014-09-22 2015-01-21 西安交通大学 Transmission method in SWIPT system under interference channel
CN105610485A (en) * 2015-12-21 2016-05-25 东南大学 Wireless relay communication system SWIPT (Simultaneous Wireless Information and Power Transfer) method
KR20170120811A (en) * 2016-04-22 2017-11-01 성균관대학교산학협력단 Full-duplex relay communication method based on energy harvesting and full-duplex relay communication system based on energy harvesting
CN108039942A (en) * 2017-12-11 2018-05-15 天津工业大学 Passive relay collects interfering energy to improve optimal system Speed method in SWIPT systems
CN109088686A (en) * 2018-09-28 2018-12-25 西北工业大学 One kind being based on wireless messages and energy transmission method while 5G height frequency range

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于携能通信的大规模无线协作网络终端性能分析;姚媛媛，李学华，黄赛，朱政宇，向维;《北京邮电大学学报》;20181231;全文 *
无线携能OFDM中继***的联合资源分配研究;何丽君，张广驰，黄高飞，万林青，崔苗，刘怡俊，林凡;《广东工业大学学报》;20180131;全文 *

Also Published As

Publication number Publication date
CN109890039A (en) 2019-06-14

Similar Documents

Publication Publication Date Title
CN107508628B (en) Cooperative transmission method in radio frequency energy collection relay network
CN101977097B (en) Method for combined adjustment of adaptive rate and power of remote WiFi link
CN102573033B (en) Multi-Femtocell downlink power interference control method based on game theory
CN108811025B (en) Wireless energy-carrying communication cooperative transmission scheme based on relay energy storage
CN110267294B (en) Random relay selection method based on energy cooperation
CN110784882B (en) Energy acquisition D2D communication resource allocation method based on reinforcement learning
CN105610485A (en) Wireless relay communication system SWIPT (Simultaneous Wireless Information and Power Transfer) method
Coutinho et al. Modeling and analysis of opportunistic routing in low duty-cycle underwater sensor networks
Liu et al. Energy-efficient transmission scheme for mobile data gathering in wireless sensor networks
CN110518993B (en) Resource allocation method of wireless energy collection type full-duplex cooperative relay system
CN108540247B (en) Spectrum sensing and information transmission method with energy collection in cognitive wireless network
CN109151775B (en) Relay-based D2D communication transmission mode switching method between cells
CN103096334A (en) Method and device for determining community extension bias boundary
CN110602758A (en) Cognitive energy-carrying relay communication method based on multi-slot wireless energy collection
CN109890039B (en) Method for allocating SWIPT relay resources under general interference
Iqbal et al. GWINs: Group-based medium access for large-scale wireless powered IoT networks
CN110366225B (en) Wireless energy supply multi-hop communication system node selection method
Zhang et al. Full-duplex cooperative cognitive radio networks with wireless energy harvesting
CN102196587B (en) Wireless-resource-dispatching method during multi-cell cooperation in relay-aided communication system
CN102395188A (en) Communication relay terminal and combined optimization method for relay selection and power distribution of communication relay terminal
KR20130142426A (en) Multihop transmission method for increasing node's lifetime in wireless ad hoc network
Nguyen et al. Wireless powered D2D communications underlying cellular networks: design and performance of the extended coverage
CN112073094A (en) TS (transport stream) strategy-based method for optimizing throughput in SWIPT (switched wire power over IP) relay system including direct link
CN107359926B (en) Full-duplex relay transmission method based on energy state
Hou et al. Bidirectional wireless information and power transfer with an energy accumulating relay

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210625

Termination date: 20220130

CF01 Termination of patent right due to non-payment of annual fee