CN108988912B - Transceiver joint optimization method and device - Google Patents

Transceiver joint optimization method and device Download PDF

Info

Publication number
CN108988912B
CN108988912B CN201810993919.3A CN201810993919A CN108988912B CN 108988912 B CN108988912 B CN 108988912B CN 201810993919 A CN201810993919 A CN 201810993919A CN 108988912 B CN108988912 B CN 108988912B
Authority
CN
China
Prior art keywords
optimal
matrix
component
constraint condition
precoding matrix
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.)
Active
Application number
CN201810993919.3A
Other languages
Chinese (zh)
Other versions
CN108988912A (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.)
Sichuan University
Original Assignee
Sichuan 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 Sichuan University filed Critical Sichuan University
Priority to CN201810993919.3A priority Critical patent/CN108988912B/en
Publication of CN108988912A publication Critical patent/CN108988912A/en
Application granted granted Critical
Publication of CN108988912B publication Critical patent/CN108988912B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/022—Site diversity; Macro-diversity
    • H04B7/026—Co-operative diversity, e.g. using fixed or mobile stations as relays
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0426—Power distribution
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03891—Spatial equalizers
    • H04L25/03898—Spatial equalizers codebook-based design
    • H04L25/0391—Spatial equalizers codebook-based design construction details of matrices

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention provides a transceiver joint optimization method and a transceiver joint optimization device, which are applied to a wireless energy supply MIMO relay system based on a TS protocol. The method comprises the following steps: acquiring a first channel matrix, a second channel matrix, an energy conversion rate, a nominal power, a first peak power and a second peak power; performing singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix and a second diagonal matrix; and optimizing the time switching factor, the first precoding matrix, the second precoding matrix and the third precoding matrix of the system according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix and the second diagonal matrix and the energy constraint condition and the power constraint condition. The method can optimize the communication process of each node of the wireless energy supply MIMO relay system so as to enhance the signal transmission efficiency of the system.

Description

Transceiver joint optimization method and device
Technical Field
The invention relates to the technical field of data communication, in particular to a method and a device for joint optimization of transceivers.
Background
With the continuous development of data communication technology, MIMO (Multiple-Input Multiple-Output) communication technology is applied more and more widely. However, currently, the existing MIMO relay communication system implements signal transmission communication by performing power constraint on signal transmission based on rated power. When the MIMO relay communication system performs signal communication by such a communication method, the signal transmission efficiency of the communication system is lowered due to the limitation of the rated power.
Disclosure of Invention
In order to overcome the above-mentioned deficiencies in the prior art, the present invention provides a method and an apparatus for joint optimization of transceivers, which can optimize a communication process of each node of a wireless energy-supplying MIMO relay system, thereby enhancing a signal transmission efficiency of the relay system.
Regarding to the method, an embodiment of the present invention provides a transceiver joint optimization method, which is applied to a wireless-powered MIMO relay system based on a time switching ts (time switching) protocol, where the system includes a source node, a relay node, and a destination node, where the relay node is powered wirelessly by an energy signal sent by the source node to transmit an information signal from the source node to the destination node, and the method includes:
acquiring a first channel matrix between the source node and the relay node, a second channel matrix between the relay node and the destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power on the source node, a first peak power on the source node and a second peak power on the relay node;
performing singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix;
and optimizing a time switching factor of the system, a first precoding matrix used for transmitting energy signals on the source node, a second precoding matrix used for transmitting information signals on the source node and a third precoding matrix used for transmitting information signals from the source node on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix and the second diagonal matrix according to an energy constraint condition and a power constraint condition.
As for an apparatus, an embodiment of the present invention provides a transceiver joint optimization apparatus, applied to a wireless-powered MIMO relay system based on a time-switched TS protocol, the system including a source node, a relay node, and a destination node, wherein the relay node is wirelessly powered by an energy signal sent by the source node to transmit an information signal from the source node to the destination node, the apparatus includes:
an information obtaining module, configured to obtain a first channel matrix between the source node and the relay node, a second channel matrix between the relay node and the destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power at the source node, a first peak power at the source node, and a second peak power at the relay node;
the matrix decomposition module is used for respectively carrying out singular value decomposition on the first channel matrix and the second channel matrix to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix;
a joint optimization module, configured to optimize a time switching factor of the system, a first precoding matrix used on the source node for transmitting an energy signal, a second precoding matrix used on the source node for transmitting an information signal, and a third precoding matrix used on the relay node for transmitting an information signal from the source node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix, and the second diagonal matrix according to an energy constraint condition and a power constraint condition.
Compared with the prior art, the transceiver joint optimization method and device provided by the embodiment of the invention have the following beneficial effects: the method can optimize the communication process of each node of the relay system in a mode of carrying out joint optimization on the transceiver parameters corresponding to the wireless energy supply MIMO relay system, so as to enhance the signal transmission efficiency of the relay system, wherein the transceiver parameters comprise a time switching factor of the system, a first precoding matrix used for transmitting energy signals on the source node, a second precoding matrix used for transmitting information signals on the source node, and a third precoding matrix used for transmitting information signals from the source node on the relay node. The method is applied to a wireless energy supply MIMO relay system, the system comprises a source node, a relay node and a destination node, wherein the relay node sends an energy signal by the source node for wireless energy supply so as to transmit an information signal from the source node to the destination node, and the source node, the relay node and the destination node are in signal communication based on the DF (Decode-and-Forward) technology of TS protocol. Firstly, the method obtains a first channel matrix between the source node and the relay node, a second channel matrix between the relay node and the destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power at the source node, a first peak power at the source node, and a second peak power at the relay node. Then, the method performs singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix. And finally, the method optimizes a time switching factor of the system, a first precoding matrix used for transmitting energy signals on the source node, a second precoding matrix used for transmitting information signals on the source node and a third precoding matrix used for transmitting information signals from the source node on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix and the second diagonal matrix according to an energy constraint condition and a power constraint condition, so that the time switching factor and each precoding matrix of the system are adjusted and optimized in an energy constraint and power constraint mode to improve the signal transmission efficiency and the data transmission quantity of the relay system.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments are briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention, and therefore should not be considered as limiting the scope of the claims of the present invention, and it is obvious for those skilled in the art that other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 is a flowchart illustrating a joint optimization method for transceivers according to an embodiment of the present invention.
Fig. 2 is a schematic flow chart of the sub-steps included in step S230 shown in fig. 1.
Fig. 3 is another flowchart of a joint optimization method for transceivers according to an embodiment of the present invention.
Fig. 4 is a block diagram of a joint optimization apparatus for transceivers according to an embodiment of the present invention.
Fig. 5 is another block diagram of a joint optimization apparatus for transceivers according to an embodiment of the present invention.
Icon: 100-joint optimization means of the transceivers; 110-an information acquisition module; 120-a matrix decomposition module; 130-a joint optimization module; 140-configuration module.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, 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 some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. 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.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
In the description of the present invention, it should be noted that the terms "first", "second", "third", and the like are used only for distinguishing the description, and are not intended to indicate or imply relative importance.
Some embodiments of the invention are described in detail below with reference to the accompanying drawings. The embodiments described below and the features of the embodiments can be combined with each other without conflict.
Fig. 1 is a flowchart illustrating a method for joint optimization of transceivers according to an embodiment of the present invention. In the embodiment of the invention, the transceiver joint optimization method is applied to a wireless energy supply MIMO relay system based on a TS protocol, and is used for carrying out joint optimization on transceiver parameters corresponding to the wireless energy supply MIMO relay system so as to optimize the communication process of each node of the system. The wireless energy supply MIMO relay system carries out signal communication based on DF technology, the wireless energy supply MIMO relay system comprises a source node, a relay node and a destination node which are communicated with each other, the source node sends an energy signal to the relay node to carry out wireless function on the relay node, the source node sends an information signal to the relay node, and the relay node transmits the information signal to the destination node based on the energy corresponding to the energy signal, thereby realizing the communication process of the wireless energy supply MIMO relay system. The transceiver parameters corresponding to the system include a time switching factor of the system, a first precoding matrix used for transmitting an energy signal on the source node, a second precoding matrix used for transmitting an information signal on the source node, and a third precoding matrix used for transmitting an information signal from the source node on the relay node.
In this embodiment, if a signal transmission cycle of a wireless energy supply MIMO relay system for performing transmission communication on an information signal is T, the signal transmission cycle may be divided into three time periods, where a corresponding time duration of a first time period is α T, and at this time, the source node sends an energy signal to the relay node; the corresponding duration of the second time period is (1-alpha) T/2, and at the moment, the source node sends the information signal to the relay node; and the corresponding duration of the third time period is also (1-alpha) T/2, and at the moment, the relay node sends the information signal to the destination node by using the energy corresponding to the energy signal. Wherein α is a time switching factor of the system, and its value is greater than 0 and less than 1. The specific procedures and steps of the joint optimization method for transceivers shown in fig. 1 are described in detail below.
Step S210, a first channel matrix between a source node and a relay node, a second channel matrix between the relay node and a destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power at the source node, a first peak power at the source node, and a second peak power at the relay node are obtained.
In this embodiment, there is an additive white gaussian noise v with an average value of 0 between the source node and the relay noderAn additive white Gaussian noise v with an average value of 0 exists between the relay node and the destination nodedThus, for the relay node, additive white gaussian noise v between the source node and the relay node is availablerVariance of (2)
Figure BDA0001781432410000061
Representing the strength of the anti-interference capability when the relay node receives signals; for a destination node, an additive white Gaussian noise v between the relay node and the destination node can be useddVariance of (2)
Figure BDA0001781432410000062
And the strength of the anti-interference capability when the relay node receives the signal is represented.
In this embodiment, the source node, the relay node, and the destination node may all be provided with a plurality of antennas, where the number of antennas on the source node may be NsMeans that the number of antennas on the relay node is available NrMeans that the number of antennas at the destination node is available as NdIf the first channel matrix is represented by N, the first channel matrix is Nr×NsIs denoted by H(ii) a The second channel matrix is Nd×NrMay be denoted by G.
In this embodiment, the nominal power represents an average power that can be obtained by the source node, and P is availableSRepresents; the energy conversion rate represents the conversion efficiency of extracting corresponding energy from the energy signal by the relay node after receiving the energy signal from the source node, and can be represented by eta; the first peak power represents a maximum transmission power of the source node, available Pm,sRepresents; the second peak power represents a maximum transmission power, available P, of the relay nodem,rAnd (4) showing.
Step S220, performing singular value decomposition on the first channel matrix and the second channel matrix, respectively, to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix.
In this embodiment, the method may use the first channel matrix H and the second channel matrix G to be represented by the following formula by performing singular value decomposition on the first channel matrix H and the second channel matrix G, respectively:
Figure BDA0001781432410000071
wherein, ΛhRepresenting a first diagonal matrix, U, corresponding to the first channel matrix HhRepresenting the left sub-matrix, V, corresponding to the first channel matrix H after singular value decompositionhRepresenting the right sub-matrix, Λ, corresponding to the first channel matrix H after singular value decompositiongRepresenting a second diagonal matrix, U, corresponding to a second channel matrix GgRepresenting the left sub-matrix, V, corresponding to the second channel matrix G after singular value decompositiongRepresenting the right-hand submatrix of the second channel matrix G after singular value decomposition (·)HRepresenting the hermitian conjugate transpose. The diagonal elements in the first diagonal matrix corresponding to the first channel matrix H are sequentially arranged according to a decreasing order, and the diagonal elements in the second diagonal matrix corresponding to the second channel matrix G are sequentially arranged according to a decreasing orderThe sequence is arranged in sequence.
Step S230, optimizing a time switching factor of the system, a first precoding matrix on the source node, a second precoding matrix on the source node, and a third precoding matrix on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix, and the second diagonal matrix according to an energy constraint condition and a power constraint condition.
In this embodiment, the first precoding matrix is used to transmit an energy signal to the relay node, the second precoding matrix is used to transmit an information signal to the relay node, and the third precoding matrix is used to transmit an information signal from the source node to the destination node. When the system is used for carrying out transmission communication on an information signal, the energy signal is required to be precoded by the source node through the first precoding matrix and then sent to the relay node, and then the information signal is precoded by the source node through the second precoding matrix and then sent to the relay node. After receiving the encoded energy signal and the encoded information signal, the relay node decodes the encoded energy signal and the encoded information signal to obtain energy corresponding to the energy signal and information corresponding to the information signal, pre-encodes the information corresponding to the information signal through a third pre-encoding matrix, and then sends the information signal processed by the third pre-encoding matrix to the destination node based on the energy corresponding to the energy signal.
Wherein the energy signal transmitted by the source node is N1Signal vector s of x 11The corresponding first precoding matrix is B1N of the expressions×N1The coding matrix of (2). The information signal transmitted by the source node is N2Signal vector s of x 12The corresponding second precoding matrix is B2N of the expressions×N2The coding matrix of (2). The third precoding matrix is N represented by Fr×N2Of the coding matrix, correspondingThe information signal which is coded by the third precoding matrix and is directly transmitted to the destination node by the relay node is N2X 1 signal vector
Figure BDA0001781432410000084
In this embodiment, there are energy constraints and power constraints in the system. The power constraint includes a first peak power Pm,sCorresponding first peak constraint condition, and second peak power Pm,rA corresponding second peak constraint requiring that the transmission power at the source node is not greater than the first peak power Pm,sThe second peak constraint requires that the transmission power at the relay node is not greater than the second peak power Pm,r. The energy constraint may be represented by the following equation:
Figure BDA0001781432410000081
Figure BDA0001781432410000082
Figure BDA0001781432410000083
where tr (-) denotes the trace of the matrix, PSRepresents the current nominal power of the source node, and MI (-) represents the interaction information between the source node and the destination nodeHRepresenting the hermitian conjugate transpose. In this embodiment, the method performs a time switching factor α on the system and a first precoding matrix B on the source node1A second precoding matrix B at the source node2And when the optimal value of the third precoding matrix F on the relay node is optimized, the energy constraint condition and the power constraint condition are both required to be met, so that the communication process of the system can be connectedAnd (5) synthesizing and optimizing.
In the present embodiment, the mutual information MI (α, B) in the above energy constraint condition2F) can be represented by the following formula:
Figure BDA0001781432410000091
Figure BDA0001781432410000092
wherein the content of the first and second substances,
Figure BDA0001781432410000093
is a number N2×N2The identity matrix, |, represents the matrix determinant,
Figure BDA0001781432410000094
representing additive white Gaussian noise v between the source node and the relay noderThe variance of (a) is determined,
Figure BDA0001781432410000095
representing additive white Gaussian noise v between the relay node and the destination nodedVariance of, N2Is not greater than min { rank (h), rank (g) }, rank () represents the rank of the matrix, and the remaining characters represent meanings as described above.
In this embodiment, after obtaining the singular value decomposition expressions corresponding to the first channel matrix H and the second channel matrix G, the method may decompose the expressions according to the first precoding matrix B1The second precoding matrix B2And the incidence relation between the third precoding matrix F and the first channel matrix H and the second channel matrix G, and the first precoding matrix B1The second precoding matrix B2And the respective corresponding optimal structures of the third precoding matrices F are represented by the following equation:
Figure BDA0001781432410000096
wherein, (.)*Denotes the optimum value, λbRepresenting a first precoding matrix B1Corresponding positive definite scalar quantity, vh,1Indicating that the first channel matrix H corresponds to the right sub-matrix VhFirst row of (V)h,1Indicating that the first channel matrix H corresponds to the right sub-matrix VhN on the leftmost side2Column, Λ2Representing a second precoding matrix B2Corresponding to N2×N2Diagonal matrix of, Vg,1Indicating that the second channel matrix G corresponds to the right sub-matrix VgN on the leftmost side2Column, ΛfN corresponding to the third precoding matrix F2×N2The diagonal matrix of (a).
In this case, the above equation for the transceiver parameter optimization criterion with the matrix variable representing the energy constraint condition, in combination with the power constraint condition including the first peak constraint condition and the second peak constraint condition, may obtain the following power allocation optimization equation with the scalar variable representing the power allocation optimization criterion:
Figure BDA0001781432410000101
Figure BDA0001781432410000102
Figure BDA0001781432410000103
Figure BDA0001781432410000104
Figure BDA0001781432410000105
Figure BDA0001781432410000106
wherein the content of the first and second substances,
Figure BDA00017814324100001011
λ2,iis represented by2The ith diagonal element of (2), λf,iIs represented byfThe ith diagonal element of (1) (. C)TExpressing the transposition of the matrix, wherein the expression of the first peak value constraint condition corresponding to the source node in the power distribution optimization expression is 0 < lambdab≤Pm，s
Figure BDA0001781432410000107
The expression of the second peak value constraint condition corresponding to the relay node in the power distribution optimization expression is
Figure BDA0001781432410000108
In this example, the method introduces
Figure BDA0001781432410000109
Then the pair i is 1, N2The following parameters exist:
Figure BDA00017814324100001010
wherein x isiRepresenting a corresponding first basis component, y, of said systemiRepresenting a corresponding second basis component of said system, aiRepresenting a corresponding third basis component of said system, biRepresenting a corresponding fourth basis component, λ, of said systemh,iIs represented byhThe ith diagonal element of (2), λg,iIs represented bygThe ith diagonal element of (1).
Whereas for a given temporal switching factor a the relaxation component t can be expressed
Figure BDA0001781432410000111
At this time, the power distribution optimization equation can be converted into the following system component optimization equation for representation:
Figure BDA00017814324100001113
Figure BDA0001781432410000112
Figure BDA0001781432410000113
Figure BDA0001781432410000114
Figure BDA0001781432410000115
Figure BDA0001781432410000116
at this time, the first fundamental component xiThere is a first basis component x for solvingiThe secondary expression equation of the optimal structure of (1) is as follows:
Figure BDA0001781432410000117
Figure BDA0001781432410000118
the second basis component yiThere is a second basis component y for solvingiIs the most important ofThe secondary expression equation for the preferred structure is as follows:
Figure BDA0001781432410000119
Figure BDA00017814324100001110
in this embodiment, the method may be implemented by solving for an optimal time-switching factor α*Optimal relaxation component t*The optimal first basis component
Figure BDA00017814324100001111
And an optimal second basis component
Figure BDA00017814324100001112
To obtain the optimal first precoding matrix B corresponding to the system1 *The second precoding matrix B is optimal2 *The optimal third precoding matrix F*And realizing the joint optimization of the transceiver parameters of the system.
Fig. 2 is a schematic flow chart illustrating the sub-steps included in step S230 shown in fig. 1. In the embodiment of the present invention, the step S230 includes a sub-step S231, a sub-step S232, a sub-step S233, and a sub-step S234.
And a substep S231, calculating a current optimal time switching factor of the system meeting the power constraint condition based on a golden section search method according to the first diagonal matrix, the second diagonal matrix, the nominal power, the first peak power, the second peak power and the energy conversion rate, and correspondingly obtaining an optimal positive definite scalar quantity corresponding to the first precoding matrix meeting the power constraint condition, an optimal diagonal matrix corresponding to the second precoding matrix and an optimal diagonal matrix corresponding to the third precoding matrix.
In this embodiment, the method is based on the first diagonal matrix ΛhAnd the second diagonal matrix ΛgCalculating a corresponding third basic component a of the systemiAnd a fourth basis component b corresponding to said systemiWherein i is 1, N2. The method will then be based on the first peak power Pm,sSecond peak power Pm,rThe third base component aiAnd the fourth basic component biSolving the optimal time switching factor alpha according to the golden section searching method*And is in alpha*In the solving process, the first precoding matrix B is obtained based on the system parameter optimization expression1Corresponding optimal positive definite scalar lambdab *The second precoding matrix B2Corresponding optimal diagonal matrix Λ2 *And the optimal diagonal matrix Lambda corresponding to the third pre-coding matrix Ff *。
Wherein, the method can obtain the optimal time switching factor alpha through the following similar code programs*The first precoding matrix B1Corresponding optimal positive definite scalar lambdab *The second precoding matrix B2Corresponding optimal diagonal matrix Λ2 *And the optimal diagonal matrix Lambda corresponding to the third pre-coding matrix Ff *：
Initialization:αl＝0 andαu＝1；
While|αu-αl|>εdo
Defineν1＝(δ-1)αl+(2-δ)αu andν2＝(2-δ)αl+(δ-1)αu；
Based on alpha ═ v1Solving the system component optimization formula to obtain alpha ═ v1Optimal relaxation component t of time*The optimal first basis component
Figure BDA0001781432410000121
And an optimal second basis component
Figure BDA0001781432410000122
And are basedAt α ═ v1Calculating F (v)1)；
Based on alpha ═ v2Solving the system component optimization formula to obtain alpha ═ v2Optimal relaxation component t of time*The optimal first basis component
Figure BDA0001781432410000131
And an optimal second basis component
Figure BDA0001781432410000132
And based on alpha ═ v2Calculating F (v)2)；
Figure BDA0001781432410000139
α*＝(αu+αl)/2；
Figure BDA0001781432410000133
Figure BDA0001781432410000134
Wherein ε is a normal number close to 0, and δ is a reduction factor equal to 1.618 (.)*Representing the optimum value, αlFor time switching a lower value of the factor alpha, alphauIs an upper limit value, v, of the time-switching factor alpha1And v2Solving an optimal time-switching factor alpha for a golden section search method*A switching factor, lambda, occurring iteratively in the process of (2)b *For the first precoding matrix B1The corresponding optimal positive definite scalar quantity,
Figure BDA0001781432410000135
for the second precoding matrix B2Corresponding optimal diagonal matrix Λ2 *The (i) th diagonal element of (a),
Figure BDA0001781432410000136
is the optimal diagonal matrix Lambda corresponding to the third pre-coding matrix Ff *The function F (alpha) is available for the ith diagonal element of (1)
Figure BDA0001781432410000137
That M (α) is α corresponds to the system component optimization formula
Figure BDA0001781432410000138
The value of (c).
It should be understood that the above-mentioned code program is only one implementation manner of the embodiment of the present invention, and should not be construed as limiting the scope of the present invention. In this embodiment, the step of correspondingly obtaining the optimal positive definite scalar corresponding to the first precoding matrix, the optimal diagonal matrix corresponding to the second precoding matrix, and the optimal diagonal matrix corresponding to the third precoding matrix, which satisfy the power constraint condition in the sub-step S231, includes:
in the process of calculating the optimal time switching factor, calculating a first basic component which corresponds to the switching factor occurring in iteration and satisfies a first peak constraint condition and an optimal second basic component which satisfies a second peak constraint condition according to a binary search method based on the nominal power, the diagonal element set of the first diagonal matrix, the diagonal element set of the second diagonal matrix, the energy conversion rate and the switching factor occurring in each iteration;
according to the optimal first basic component and the optimal second basic component corresponding to the optimal time switching factor obtained through final calculation and the incidence relation among the first basic component, the second basic component, the positive definite scalar corresponding to the first pre-coding matrix, the diagonal matrix corresponding to the second pre-coding matrix and the diagonal matrix corresponding to the third pre-coding matrix, the optimal positive definite scalar corresponding to the first pre-coding matrix, the optimal diagonal matrix corresponding to the second pre-coding matrix and the optimal diagonal matrix corresponding to the third pre-coding matrix are obtained through calculation.
The execution process of calculating the optimal first basis component meeting the first peak constraint condition and the optimal second basis component meeting the second peak constraint condition corresponding to the switching factor which appears in the iteration according to the binary search method is that the execution process is based on alpha-v in the code program1Solving the system component optimization formula to obtain alpha ═ v1Optimal relaxation component t of time*The optimal first basis component
Figure BDA0001781432410000143
And an optimal second basis component
Figure BDA0001781432410000144
"and" is based on α ═ v2Solving the system component optimization formula to obtain alpha ═ v2Optimal relaxation component t of time*The optimal first basis component
Figure BDA0001781432410000145
And an optimal second basis component
Figure BDA0001781432410000146
"wherein v is1And v2Solving an optimal time-switching factor alpha for a golden section search method*Iteratively occurring switching factors. The implementation can be represented by the following similar code programs:
Initialization:tl and tu；
While|tu-tl|>εdo
t*＝(αu+αl)/2；
Figure BDA0001781432410000141
wherein the value of beta is solved for the first basis component x by binary searchiComplementary equation of relaxation
Figure BDA0001781432410000142
Obtaining;
Figure BDA0001781432410000151
wherein the value of gamma is solved for the second basis component y by binary searchiComplementary equation of relaxation
Figure BDA0001781432410000152
Obtaining;
Figure BDA0001781432410000153
judgment of x0Whether the first peak value constraint condition is met or not is judged0Whether a second peak constraint condition is satisfied;
according to the judgment result and the first basic component xiIs expressed in terms of the secondary expression equation and the second basis component yiIs calculated to obtain the current optimal relaxation component t*The optimal first basis component
Figure BDA0001781432410000154
And an optimal second basis component
Figure BDA0001781432410000155
Figure BDA0001781432410000156
Wherein, tlIs the lower limit value of the relaxation component t, tuIs the upper limit value of the relaxation component t, t0For calculating a corresponding calculated relaxation component, x, within an execution loop of an optimal relaxation component based on said iteratively occurring switching factors0Is given as0Corresponding first basis component x basediThe first basis component, y, calculated from the relaxation complementary equation (as shown above)0Is given as0Corresponding based on the second basis component yiComplementary equations of relaxation (see above)Shown) the calculated second basis component.
In this embodiment, the calculation of the switching factor v occurring in the iteration according to the binary search method1Or v2Corresponding optimal first basis component satisfying first peak constraint condition
Figure BDA0001781432410000159
And an optimal second basis component satisfying a second peak constraint
Figure BDA0001781432410000158
Comprises the following steps:
at a switching factor v occurring based on said iteration1Or v2Calculating an optimal relaxation component t*In the loop (2), the calculated relaxation components t are respectively compared with each other according to a binary search method0The matched first basis component xiAnd the second basis component yiThe relaxation complementary equation of (a) is solved to obtain a corresponding first basis component x0And a second basis component y0；
Judging the calculated first basis component x0Judging whether the first peak value constraint condition is met or not, and judging the calculated second basic component y0Whether a second peak constraint condition is satisfied;
if the first basis component x0Satisfies the first peak constraint condition, and the second basis component y0If the second peak value constraint condition is satisfied, the calculated relaxation component t is used0The first basis component x0And the second basis component y0Respectively as the optimal relaxation component t of the current cycle*The optimal first basis component
Figure BDA0001781432410000161
And an optimal second basis component
Figure BDA0001781432410000162
Wherein, the method is as followsWill x0Substitution into
Figure BDA0001781432410000163
The mode of checking is carried out, and x is judged0Whether the first peak constraint is satisfied. The method comprises mixing a compound of formula (I) and (II)0Substitution into
Figure BDA0001781432410000164
The mode of checking is carried out, and y is judged0Whether the second peak constraint is satisfied.
In this embodiment, the calculation of the switching factor v occurring in the iteration according to the binary search method1Or v2Corresponding optimal first basis component satisfying first peak constraint condition
Figure BDA0001781432410000165
And an optimal second basis component satisfying a second peak constraint
Figure BDA0001781432410000166
Further comprising the steps of:
if the first basis component x0Satisfies the first peak constraint condition, and the second basis component y0If the second peak constraint condition is not satisfied, the first base component x is calculated0First basis component as optimum for the current cycle
Figure BDA0001781432410000167
Solving the second basis component y according to water injectioniThe corresponding secondary expression equation obtains the optimal second basic component of the current cycle
Figure BDA0001781432410000168
And based on the optimal second basis component
Figure BDA0001781432410000169
Calculating and obtaining the optimal relaxation component t of the current cycle*。
Wherein, the methodMethod for solving second basis component y according to water injection methodiThe corresponding secondary expression equation obtains the optimal second basic component of the current cycle
Figure BDA00017814324100001610
Then, will be according to the formula
Figure BDA00017814324100001611
Calculating and obtaining the optimal relaxation component t of the current cycle*。
In this embodiment, the calculation of the switching factor v occurring in the iteration according to the binary search method1Or v2Corresponding optimal first basis component satisfying first peak constraint condition
Figure BDA00017814324100001612
And an optimal second basis component satisfying a second peak constraint
Figure BDA00017814324100001613
Further comprising the steps of:
if the first basis component x0The first peak constraint is not satisfied, and the second basis component y0If the second peak value constraint condition is satisfied, the calculated second basic component y is used0Second basis component as optimum for the current cycle
Figure BDA0001781432410000171
Solving the first basis component x according to a water injection methodiThe corresponding secondary expression equation obtains the optimal first basic component of the current cycle
Figure BDA0001781432410000172
And based on the optimal first basis component
Figure BDA0001781432410000173
Calculating and obtaining the optimal relaxation component t of the current cycle*。
Wherein the method comprises the following steps of obtaining the target product according to a water injection methodSolving the first basis component xiThe corresponding secondary expression equation obtains the optimal first basic component of the current cycle
Figure BDA0001781432410000174
Then, will be according to the formula
Figure BDA00017814324100001718
Calculating and obtaining the optimal relaxation component t of the current cycle*。
In this embodiment, the calculation of the switching factor v occurring in the iteration according to the binary search method1Or v2Corresponding optimal first basis component satisfying first peak constraint condition
Figure BDA0001781432410000176
And an optimal second basis component satisfying a second peak constraint
Figure BDA0001781432410000177
Further comprising the steps of:
if the first basis component x0The first peak constraint is not satisfied, and the second basis component y0If the second peak value constraint condition is not met, the first basic component x is respectively solved according to a water injection methodiCorresponding secondary representation equation and second basis component yiCorresponding secondary expression equation to obtain the optimal first basic component of the current cycle
Figure BDA0001781432410000178
And an optimal second basis component
Figure BDA0001781432410000179
Based on the optimal first basis components, respectively
Figure BDA00017814324100001710
And the optimal second basis component
Figure BDA00017814324100001711
Calculating to obtain two optimal relaxation components corresponding to the two basic components, and selecting the relaxation component with smaller value in the two calculated relaxation components as the optimal relaxation component t of the current cycle*。
Wherein the optimal first basis component is based on
Figure BDA00017814324100001712
The calculated value of the relaxation component is
Figure BDA00017814324100001713
Based on the second basis component of the optimum
Figure BDA00017814324100001714
The calculated value of the relaxation component is
Figure BDA00017814324100001715
Then the optimal relaxation component t of the current sub-loop*I.e. the slack component with the smaller value of the two slack components.
In this embodiment, the method determines the optimal time-shift factor α*And obtaining the optimal first basic component matched with the switching factor appearing in the last iteration in the implementation process of the golden section search method
Figure BDA00017814324100001716
And an optimal second basis component
Figure BDA00017814324100001717
Then, the optimal first basic component is used
Figure BDA0001781432410000181
And the optimal second basis component
Figure BDA0001781432410000182
Respectively as a time switching factor alpha with respect to said optimum*Corresponding optimumFirst basis component of
Figure BDA0001781432410000183
And the optimal second basis component
Figure BDA0001781432410000184
And calculating to obtain the first precoding matrix B1Corresponding optimal positive definite scalar lambdab *The second precoding matrix B2Corresponding optimal diagonal matrix Λ2 *And the optimal diagonal matrix Lambda corresponding to the third pre-coding matrix Ff *。
And a substep S232, calculating to obtain the currently optimal first precoding matrix according to the currently corresponding optimal positive definite scalar of the first precoding matrix and the incidence relation between the first precoding matrix and the first channel matrix.
And a substep S233, calculating to obtain the currently optimal second precoding matrix according to the currently corresponding optimal diagonal matrix of the second precoding matrix and the incidence relation between the second precoding matrix and the first channel matrix.
And a substep S234, calculating to obtain the currently optimal third precoding matrix according to the currently corresponding optimal diagonal matrix of the third precoding matrix and the incidence relation between the third precoding matrix and the second channel matrix.
In this embodiment, the method obtains the optimal time switching factor α*The first precoding matrix B1Corresponding optimal positive definite scalar lambdab *The second precoding matrix B2Corresponding optimal diagonal matrix Λ2 *And the optimal diagonal matrix Lambda corresponding to the third pre-coding matrix Ff *Then, according to the first precoding matrix B1The second precoding matrix B2And an optimal structural expression of the third precoding matrix F
Figure BDA0001781432410000185
And calculating to obtain three optimal precoding matrixes corresponding to the system.
Step S240, configuring a signal transmission cycle of the system by using the optimized time-switching factor, configuring the source node by using the optimized first precoding matrix and the optimized second precoding matrix, and configuring the relay node by using the optimized third precoding matrix, so that the system performs communication based on the optimized precoding matrix and the optimized time-switching factor.
Fig. 4 is a block diagram of a joint optimization apparatus 100 for transceivers according to an embodiment of the present invention. In the embodiment of the present invention, the transceiver joint optimization apparatus 100 is applied to the above-mentioned wireless energy-supplying MIMO relay system based on the TS protocol, and the transceiver joint optimization apparatus 100 includes an information obtaining module 110, a matrix decomposition module 120, and a joint optimization module 130.
The information obtaining module 110 is configured to obtain a first channel matrix between a source node and a relay node, a second channel matrix between the relay node and a destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power at the source node, a first peak power at the source node, and a second peak power at the relay node.
In this embodiment, the information obtaining module 110 may perform step S210 shown in fig. 1, and the detailed description may refer to the above detailed description of step S210.
The matrix decomposition module 120 is configured to perform singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix.
In this embodiment, the matrix decomposition module 120 may execute step S220 shown in fig. 1, and the detailed description may refer to the above detailed description of step S220.
The joint optimization module 130 is configured to optimize a time switching factor of the system, a first precoding matrix used for transmitting an energy signal on the source node, a second precoding matrix used for transmitting an information signal on the source node, and a third precoding matrix used for transmitting an information signal from the source node on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix, and the second diagonal matrix according to an energy constraint condition and a power constraint condition.
In this embodiment, the joint optimization module 130 may perform step S230 shown in fig. 1 and sub-steps S231, S232, S233, and S234 shown in fig. 2, and the detailed description may refer to the detailed description of step S220, S231, S232, S233, and S234.
Fig. 5 is a block diagram of another transceiver joint optimization apparatus 100 according to an embodiment of the present invention. In the embodiment of the present invention, the transceiver joint optimization device 100 may further include a configuration module 140.
The configuration module 140 is configured to configure a signal transmission cycle of the system by using the optimized time switching factor, configure the source node by using the optimized first precoding matrix and the optimized second precoding matrix, and configure the relay node by using the optimized third precoding matrix, so that the system performs communication based on the optimized precoding matrix and the optimized time switching factor.
In summary, in the transceiver joint optimization method and apparatus provided by the present invention, the method can optimize the communication process of each node of the relay system in a manner of performing joint optimization on transceiver parameters corresponding to the wireless energy supply MIMO relay system, so as to enhance the signal transmission efficiency of the relay system, where the transceiver parameters include a time switching factor of the system, a first precoding matrix used for transmitting an energy signal at the source node, a second precoding matrix used for transmitting an information signal at the source node, and a third precoding matrix used for transmitting an information signal from the source node at the relay node. The method is applied to the wireless energy supply MIMO relay system based on the TS protocol. Firstly, the method obtains a first channel matrix between a source node and a relay node, a second channel matrix between the relay node and a destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power on the source node, a first peak power on the source node and a second peak power on the relay node. Then, the method performs singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix. And finally, the method optimizes a time switching factor of the system, a first precoding matrix used for transmitting energy signals on the source node, a second precoding matrix used for transmitting information signals on the source node and a third precoding matrix used for transmitting information signals from the source node on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix and the second diagonal matrix according to an energy constraint condition and a power constraint condition, so that the time switching factor and each precoding matrix of the system are adjusted and optimized in an energy constraint and power constraint mode to improve the signal transmission efficiency and the data transmission quantity of the relay system.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (6)

1. A method for joint optimization of transceivers, applied to a wireless-powered MIMO relay system based on a time-switched TS protocol, the system comprising a source node, a relay node and a destination node, wherein the relay node is wirelessly powered by an energy signal sent by the source node to transmit an information signal from the source node to the destination node, the method comprising:
acquiring a first channel matrix between the source node and the relay node, a second channel matrix between the relay node and the destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power on the source node, a first peak power on the source node and a second peak power on the relay node;
performing singular value decomposition on the first channel matrix and the second channel matrix respectively to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix;
optimizing a time switching factor of the system, a first precoding matrix used for transmitting energy signals on the source node, a second precoding matrix used for transmitting information signals on the source node and a third precoding matrix used for transmitting information signals from the source node on the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix and the second diagonal matrix according to an energy constraint condition and a power constraint condition;
wherein the step of optimizing a time switching factor of the system, a first precoding matrix used for transmitting an energy signal at the source node, a second precoding matrix used for transmitting an information signal at the source node, and a third precoding matrix used for transmitting an information signal from the source node at the relay node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix, and the second diagonal matrix according to an energy constraint condition and a power constraint condition comprises:
calculating a current optimal time switching factor of the system meeting the power constraint condition according to the first diagonal matrix, the second diagonal matrix, the nominal power, the first peak power, the second peak power and the energy conversion rate based on a golden section search method, and correspondingly obtaining an optimal positive definite scalar corresponding to the first precoding matrix meeting the power constraint condition, an optimal diagonal matrix corresponding to the second precoding matrix and an optimal diagonal matrix corresponding to the third precoding matrix;
calculating to obtain a currently optimal first precoding matrix according to a currently corresponding optimal positive definite scalar quantity of the first precoding matrix and an incidence relation between the first precoding matrix and the first channel matrix;
calculating to obtain a currently optimal second precoding matrix according to the currently corresponding optimal diagonal matrix of the second precoding matrix and the incidence relation between the second precoding matrix and the first channel matrix;
calculating to obtain a currently optimal third precoding matrix according to a currently corresponding optimal diagonal matrix of the third precoding matrix and an incidence relation between the third precoding matrix and the second channel matrix;
in this process, the power constraint condition includes a first peak constraint condition corresponding to a first peak power and a second peak constraint condition corresponding to a second peak power, and the step of correspondingly obtaining an optimal positive scalar corresponding to the first precoding matrix, an optimal diagonal matrix corresponding to the second precoding matrix, and an optimal diagonal matrix corresponding to the third precoding matrix, which satisfy the power constraint condition, includes:
in the process of calculating the optimal time switching factor, calculating a first basic component which corresponds to the switching factor occurring in iteration and satisfies a first peak constraint condition and an optimal second basic component which satisfies a second peak constraint condition according to a binary search method based on the nominal power, the diagonal element set of the first diagonal matrix, the diagonal element set of the second diagonal matrix, the energy conversion rate and the switching factor occurring in each iteration;
according to the optimal first basic component and the optimal second basic component corresponding to the optimal time switching factor obtained through final calculation and the incidence relation among the first basic component, the second basic component, the positive definite scalar corresponding to the first pre-coding matrix, the diagonal matrix corresponding to the second pre-coding matrix and the diagonal matrix corresponding to the third pre-coding matrix, calculating to obtain the optimal positive definite scalar corresponding to the first pre-coding matrix, the optimal diagonal matrix corresponding to the second pre-coding matrix and the optimal diagonal matrix corresponding to the third pre-coding matrix;
in addition, the step of calculating an optimal first basis component satisfying a first peak constraint condition and an optimal second basis component satisfying a second peak constraint condition corresponding to the switching factor occurring in the iteration according to a binary search method includes:
in a cycle of calculating an optimal relaxation component based on the switching factor appearing in iteration, solving a relaxation complementary equation of the first basic component and a relaxation complementary equation of the second basic component which are matched with the calculated relaxation component according to a binary search method to obtain a corresponding first basic component and a corresponding second basic component;
judging whether the calculated first basic component meets the first peak value constraint condition or not, and judging whether the calculated second basic component meets the second peak value constraint condition or not;
if the first basis component meets the first peak constraint condition and the second basis component meets a second peak constraint condition, taking the calculated slack component, the calculated first basis component and the calculated second basis component as an optimal slack component, an optimal first basis component and an optimal second basis component of the current cycle respectively;
if the first basic component meets the first peak value constraint condition and the second basic component does not meet the second peak value constraint condition, taking the calculated first basic component as the optimal first basic component of the current cycle, solving a secondary expression equation corresponding to the second basic component according to a water injection method to obtain the optimal second basic component of the current cycle, and calculating the optimal relaxation component of the current cycle based on the optimal second basic component;
meanwhile, in the process, the interaction information between the source node and the destination node in the energy constraint condition is represented by the following formula:
Figure FDA0003073879380000041
wherein MI (-) represents the mutual information,
Figure FDA00030738793800000411
is a number N2×N2The identity matrix, |, represents the matrix determinant,
Figure FDA0003073879380000049
representing additive white Gaussian noise v between the source node and the relay noderThe variance of (a) is determined,
Figure FDA00030738793800000410
representing additive white Gaussian noise v between the relay node and the destination nodedVariance of, N2Is not greater than min { rank (H), rank (G) }, rank (. cndot.) represents the rank of the matrix, (. cndot.)HRepresenting a hermitian conjugate transpose, H for representing the first channel matrix, a being a time-switching factor of the system, G for representing the second channel matrix; f is used to represent the third precoding matrix;
the system component optimization formula involved in calculating the optimal relaxation component is as follows:
Figure FDA0003073879380000042
Figure FDA0003073879380000043
Figure FDA0003073879380000044
Figure FDA0003073879380000045
Figure FDA0003073879380000046
Figure FDA0003073879380000047
xi≥0，yi≥0，i＝1，…，N2
wherein the content of the first and second substances,
Figure FDA00030738793800000412
xi＝λ2，i，
Figure FDA0003073879380000048
xifor representing a corresponding first basis component, y, of said systemiFor representing a corresponding second basis component, a, of said systemiFor representing a corresponding third basis component of the system, biFor representing a corresponding fourth basis component, λ, of said systemh,iFor representing the first diagonal matrix ΛhThe ith diagonal element of (2), λg,iFor representing the second diagonal matrix ΛgThe ith diagonal element of (2), λ2,iFor representing the second precoding matrix B2Of the diagonal matrix of (a) is the ith diagonal element of the diagonal matrix of (b), λf,iThe ith diagonal element, P, of the diagonal matrix used to represent the third precoding matrix FSFor indicating the current nominal power of the source node, and t for indicating
Figure FDA0003073879380000051
A relaxation component for effecting expression;
wherein for solving the first basis component xiThe secondary expression equation of the optimal structure of (1) is as follows:
Figure FDA0003073879380000052
Figure FDA0003073879380000053
and for solving the second basis component yiThe secondary expression equation of the optimal structure of (1) is as follows:
Figure FDA0003073879380000054
Figure FDA0003073879380000055
yi≥0，i＝1，…，N2.
wherein, Pm,sFor representing said first peak power, Pm,rFor representing said second peak power and η for representing said energy conversion rate.
2. The method of claim 1, wherein the step of computing an optimal first basis component satisfying a first peak constraint and an optimal second basis component satisfying a second peak constraint corresponding to the iteratively occurring switching factors according to a binary search further comprises:
if the first basic component does not meet the first peak value constraint condition and the second basic component meets the second peak value constraint condition, the calculated second basic component is used as the optimal second basic component of the current cycle, a secondary expression equation corresponding to the first basic component is solved according to a water injection method to obtain the optimal first basic component of the current cycle, and the optimal relaxation component of the current cycle is calculated and solved based on the optimal first basic component.
3. The method of claim 1, wherein the step of computing an optimal first basis component satisfying a first peak constraint and an optimal second basis component satisfying a second peak constraint corresponding to the iteratively occurring switching factors according to a binary search further comprises:
if the first basic component does not meet the first peak value constraint condition and the second basic component does not meet the second peak value constraint condition, respectively solving a secondary expression equation corresponding to the first basic component and a secondary expression equation corresponding to the second basic component according to a water injection method to obtain an optimal first basic component and an optimal second basic component of the current cycle;
and calculating relaxation components corresponding to the two optimal basic components respectively based on the optimal first basic component and the optimal second basic component, and selecting the relaxation component with the smaller value in the two calculated relaxation components as the optimal relaxation component of the current cycle.
4. The method according to any one of claims 1-3, further comprising:
and configuring a signal transmission period of the system by using the optimized time switching factor, configuring the source node by using the optimized first precoding matrix and the optimized second precoding matrix, and configuring the relay node by using the optimized third precoding matrix, so that the system performs communication based on the optimized precoding matrix and the optimized time switching factor.
5. A joint optimization device for transceiver applied to a wireless-powered MIMO relay system based on a time-switched TS protocol, the system comprising a source node, a relay node and a destination node, wherein the relay node is wirelessly powered by the source node to transmit an energy signal for transmitting an information signal from the source node to the destination node, the device comprising:
an information obtaining module, configured to obtain a first channel matrix between the source node and the relay node, a second channel matrix between the relay node and the destination node, an energy conversion rate of an energy signal received by the relay node from the source node, a nominal power at the source node, a first peak power at the source node, and a second peak power at the relay node;
the matrix decomposition module is used for respectively carrying out singular value decomposition on the first channel matrix and the second channel matrix to obtain a first diagonal matrix corresponding to the first channel matrix and a second diagonal matrix corresponding to the second channel matrix;
a joint optimization module, configured to optimize a time switching factor of the system, a first precoding matrix used on the source node for transmitting an energy signal, a second precoding matrix used on the source node for transmitting an information signal, and a third precoding matrix used on the relay node for transmitting an information signal from the source node according to the nominal power, the first peak power, the second peak power, the energy conversion rate, the first diagonal matrix, and the second diagonal matrix according to an energy constraint condition and a power constraint condition;
wherein the joint optimization module is specifically configured to:
calculating a current optimal time switching factor of the system meeting the power constraint condition according to the first diagonal matrix, the second diagonal matrix, the nominal power, the first peak power, the second peak power and the energy conversion rate based on a golden section search method, and correspondingly obtaining an optimal positive definite scalar corresponding to the first precoding matrix meeting the power constraint condition, an optimal diagonal matrix corresponding to the second precoding matrix and an optimal diagonal matrix corresponding to the third precoding matrix;
calculating to obtain a currently optimal first precoding matrix according to a currently corresponding optimal positive definite scalar quantity of the first precoding matrix and an incidence relation between the first precoding matrix and the first channel matrix;
calculating to obtain a currently optimal second precoding matrix according to the currently corresponding optimal diagonal matrix of the second precoding matrix and the incidence relation between the second precoding matrix and the first channel matrix;
calculating to obtain a currently optimal third precoding matrix according to a currently corresponding optimal diagonal matrix of the third precoding matrix and an incidence relation between the third precoding matrix and the second channel matrix;
in this process, the power constraint condition includes a first peak constraint condition corresponding to a first peak power and a second peak constraint condition corresponding to a second peak power, and the manner in which the joint optimization module correspondingly obtains an optimal positive definite scalar corresponding to the first precoding matrix, an optimal diagonal matrix corresponding to the second precoding matrix, and an optimal diagonal matrix corresponding to the third precoding matrix that satisfy the power constraint condition includes:
in the process of calculating the optimal time switching factor, calculating a first basic component which corresponds to the switching factor occurring in iteration and satisfies a first peak constraint condition and an optimal second basic component which satisfies a second peak constraint condition according to a binary search method based on the nominal power, the diagonal element set of the first diagonal matrix, the diagonal element set of the second diagonal matrix, the energy conversion rate and the switching factor occurring in each iteration;
according to the optimal first basic component and the optimal second basic component corresponding to the optimal time switching factor obtained through final calculation and the incidence relation among the first basic component, the second basic component, the positive definite scalar corresponding to the first pre-coding matrix, the diagonal matrix corresponding to the second pre-coding matrix and the diagonal matrix corresponding to the third pre-coding matrix, calculating to obtain the optimal positive definite scalar corresponding to the first pre-coding matrix, the optimal diagonal matrix corresponding to the second pre-coding matrix and the optimal diagonal matrix corresponding to the third pre-coding matrix;
at this time, the way that the joint optimization module calculates the optimal first basis component meeting the first peak constraint condition and the optimal second basis component meeting the second peak constraint condition corresponding to the switching factor occurring in the iteration according to the binary search method includes:
in a cycle of calculating an optimal relaxation component based on the switching factor appearing in iteration, solving a relaxation complementary equation of the first basic component and a relaxation complementary equation of the second basic component which are matched with the calculated relaxation component according to a binary search method to obtain a corresponding first basic component and a corresponding second basic component;
judging whether the calculated first basic component meets the first peak value constraint condition or not, and judging whether the calculated second basic component meets the second peak value constraint condition or not;
if the first basis component meets the first peak constraint condition and the second basis component meets a second peak constraint condition, taking the calculated slack component, the calculated first basis component and the calculated second basis component as an optimal slack component, an optimal first basis component and an optimal second basis component of the current cycle respectively;
if the first basic component meets the first peak value constraint condition and the second basic component does not meet the second peak value constraint condition, taking the calculated first basic component as the optimal first basic component of the current cycle, solving a secondary expression equation corresponding to the second basic component according to a water injection method to obtain the optimal second basic component of the current cycle, and calculating the optimal relaxation component of the current cycle based on the optimal second basic component;
meanwhile, in the process, the interaction information between the source node and the destination node in the energy constraint condition is represented by the following formula:
Figure FDA0003073879380000091
wherein MI (-) represents the mutual information,
Figure FDA0003073879380000094
is a number N2×N2The identity matrix, |, represents the matrix determinant,
Figure FDA0003073879380000092
representing additive white Gaussian noise v between the source node and the relay noderThe variance of (a) is determined,
Figure FDA0003073879380000093
representing additive white Gaussian noise v between the relay node and the destination nodedVariance of, N2Is not greater than min { rank (H), rank (G) }, rank (. cndot.) represents the rank of the matrix, (. cndot.)HRepresenting a hermitian conjugate transpose, H for representing the first channel matrix, a being a time-switching factor of the system, G for representing the second channel matrix; f is used to represent the third precoding matrix;
the system component optimization formula involved in calculating the optimal relaxation component is as follows:
Figure FDA0003073879380000101
Figure FDA0003073879380000102
Figure FDA0003073879380000103
Figure FDA0003073879380000104
Figure FDA0003073879380000105
Figure FDA0003073879380000106
xi≥0，yi≥0，i＝1，…，N2
wherein the content of the first and second substances,
Figure FDA0003073879380000107
xifor representing a corresponding first basis component, y, of said systemiFor representing a corresponding second basis component, a, of said systemiFor representing a corresponding third basis component of the system, biFor representing a corresponding fourth basis component, λ, of said systemh,iFor representing the first diagonal matrix ΛhThe ith diagonal element of (2), λg,iFor representing the second diagonal matrix ΛgThe ith diagonal element of (2), λ2,iFor representing the second precoding matrix B2Of the diagonal matrix of (a) is the ith diagonal element of the diagonal matrix of (b), λf,iThe ith diagonal element, P, of the diagonal matrix used to represent the third precoding matrix FSFor indicating the current nominal power of the source node, and t for indicating
Figure FDA0003073879380000108
A relaxation component for effecting expression;
wherein for solving the first basis component xiThe secondary expression equation of the optimal structure of (1) is as follows:
Figure FDA0003073879380000109
Figure FDA00030738793800001010
and for solving the second basis component yiThe secondary expression equation of the optimal structure of (1) is as follows:
Figure FDA00030738793800001011
Figure FDA00030738793800001012
yi≥0，i＝1，…，N2.
wherein, Pm,sFor representing said first peak power, Pm,rFor representing said second peak power and η for representing said energy conversion rate.
6. The apparatus of claim 5, further comprising:
and the configuration module is used for configuring the signal transmission period of the system by adopting the optimized time switching factor, configuring the source node by adopting the optimized first precoding matrix and the optimized second precoding matrix, and configuring the relay node by adopting the optimized third precoding matrix, so that the system performs communication based on the optimized precoding matrix and the optimized time switching factor.
CN201810993919.3A 2018-08-29 2018-08-29 Transceiver joint optimization method and device Active CN108988912B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810993919.3A CN108988912B (en) 2018-08-29 2018-08-29 Transceiver joint optimization method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810993919.3A CN108988912B (en) 2018-08-29 2018-08-29 Transceiver joint optimization method and device

Publications (2)

Publication Number Publication Date
CN108988912A CN108988912A (en) 2018-12-11
CN108988912B true CN108988912B (en) 2021-06-29

Family

ID=64546848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810993919.3A Active CN108988912B (en) 2018-08-29 2018-08-29 Transceiver joint optimization method and device

Country Status (1)

Country Link
CN (1) CN108988912B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102882570A (en) * 2012-09-28 2013-01-16 东南大学 Optimum transceiving combined processing method for communication among equipment in mobile communication network
CN103181097A (en) * 2010-09-29 2013-06-26 Lg电子株式会社 Method and apparatus for efficient feedback in wireless communication system that supports multiple antennas
CN107332595A (en) * 2017-05-22 2017-11-07 华南理工大学 A kind of MIMO wireless energies communication network maximize handling capacity method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103181097A (en) * 2010-09-29 2013-06-26 Lg电子株式会社 Method and apparatus for efficient feedback in wireless communication system that supports multiple antennas
CN102882570A (en) * 2012-09-28 2013-01-16 东南大学 Optimum transceiving combined processing method for communication among equipment in mobile communication network
CN107332595A (en) * 2017-05-22 2017-11-07 华南理工大学 A kind of MIMO wireless energies communication network maximize handling capacity method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Joint Transceiver Optimization for Wireless Information and Energy Transfer in Nonregenerative MIMO Relay Systems";Bin Li著;《IEEE Transactions on Vehicular Technology》;20180712;第1-15页 *

Also Published As

Publication number Publication date
CN108988912A (en) 2018-12-11

Similar Documents

Publication Publication Date Title
CN102725967B (en) For the method and apparatus of information feed back and precoding
US10250311B2 (en) Method and system for MIMO communication
CN109644039A (en) Method and apparatus for downlink and uplink channel state acquisition of information
CN103975537B (en) For being modulated the adaptive method and apparatus of encoding scheme in mimo systems
WO2017124967A1 (en) Multi-antenna transmission method, base station, and user terminal
HUE032168T2 (en) Method for precoding using a block diagonal matrix
CN101425830B (en) Pre-encoding codeword determining method, system and apparatus
CN106571858B (en) Hybrid beam forming transmission system
WO2008011451A2 (en) Reduced feedback transmit beamforming
CN106464331B (en) Multi-antenna data transmission method, base station, user equipment and system
CN101110623A (en) Part user feedback method for limited feedback of multiple-input multiple-output multi-user system
CN103825678A (en) Three-dimensional multi-user multi-input and multi-output (3D MU-MIMO) precoding method based on Khatri-Rao product
CN101094022B (en) Transmitter, communication system, and communication method
CN106911376B (en) Wireless energy drive transmission method and device
US10659130B2 (en) Terminal, base station, and channel information obtaining method
EP2487950B1 (en) Wireless communication system, mobile terminal and collaborative communication method
CN114900398A (en) IRS (intelligent resilient framework) assisted cloud access network downlink beam forming method for non-ideal CSI (channel state information)
CN108988912B (en) Transceiver joint optimization method and device
CN102571172B (en) User scheduling method and equipment in MIMO wireless communication system
KR20120131224A (en) Transformation device and method
CN109120319B (en) Transceiver joint optimization method and device
CN102377526B (en) Method and device for transmitting feedback information in MIMO (multiple input multiple output) communication system
CN101645757B (en) Method and device for selecting precoding matrix
WO2022184010A1 (en) Information reporting method and apparatus, first device, and second device
CN102594489B (en) Pre-coded vector determination method, device, base station and user terminal

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