CN111030795A - Transmission method for time-reversed wireless energy-carrying communication system under non-ideal channel - Google Patents

Transmission method for time-reversed wireless energy-carrying communication system under non-ideal channel Download PDF

Info

Publication number
CN111030795A
CN111030795A CN201911188194.1A CN201911188194A CN111030795A CN 111030795 A CN111030795 A CN 111030795A CN 201911188194 A CN201911188194 A CN 201911188194A CN 111030795 A CN111030795 A CN 111030795A
Authority
CN
China
Prior art keywords
user terminal
base station
channel
time
kth
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.)
Granted
Application number
CN201911188194.1A
Other languages
Chinese (zh)
Other versions
CN111030795B (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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN201911188194.1A priority Critical patent/CN111030795B/en
Publication of CN111030795A publication Critical patent/CN111030795A/en
Application granted granted Critical
Publication of CN111030795B publication Critical patent/CN111030795B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00Arrangements affording multiple use of the transmission path
    • H04L5/003Arrangements for allocating sub-channels of the transmission path
    • H04L5/0048Allocation of pilot signals, i.e. of signals known to the receiver
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00Arrangements affording multiple use of the transmission path
    • H04L5/003Arrangements for allocating sub-channels of the transmission path
    • H04L5/0053Allocation of signalling, i.e. of overhead other than pilot signals
    • H04L5/0055Physical resource allocation for ACK/NACK
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00Arrangements affording multiple use of the transmission path
    • H04L5/003Arrangements for allocating sub-channels of the transmission path
    • H04L5/0053Allocation of signalling, i.e. of overhead other than pilot signals
    • H04L5/0057Physical resource allocation for CQI

Landscapes

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

Abstract

本发明公开了一种非理想信道下时间反演无线携能通信系统的传输方法,该方法中充分考虑了时间反演特性对非理想信道状态信息下多用户多发单收无线携能通信系统可靠性和有效性的影响,利用时间反演时空聚焦特性增强了系统的效性,通过含有信道估计误差的能量‑速率闭合表达式,在非理想CSI的条件下规划优化问题,通过本发明提供的方法可以在保证一定信息可达速率条件下提升能量收集效率,相比于传统MISO‑SWIPT,提高了系统的可靠性和有效性。

Figure 201911188194

The invention discloses a transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel. In the method, the time-reversal characteristics are fully considered for the reliability of a multi-user multi-transmit single-receive wireless energy-carrying communication system under non-ideal channel state information. The efficiency of the system is enhanced by using the time-reversal space-time focusing characteristics, and the optimization problem is planned under the condition of non-ideal CSI through the energy-rate closed expression containing the channel estimation error. The method can improve the energy harvesting efficiency under the condition of ensuring a certain information reachability rate, and improve the reliability and effectiveness of the system compared with the traditional MISO‑SWIPT.

Figure 201911188194

Description

非理想信道下时间反演无线携能通信系统的传输方法Transmission method for time-reversed wireless energy-carrying communication system under non-ideal channel

技术领域technical field

本发明涉及通信技术领域,尤其涉及一种非理想信道下时间反演无线携能通信系统的传输方法。The present invention relates to the field of communication technologies, and in particular, to a transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel.

背景技术Background technique

随着通信技术的不断演进,绿色通信和万物互联成为未来通信网络主要发展趋势。手机、传感器、家用电器在内的海量设备被连接到通信网络中,实现远程控制、数据收集和数据交换等功能。低成本、低功耗的设备也是未来通信的一类典型设备,这类设备将在物联网(Internet of Things,IoT)、无线传感网(Wireless Sensor Network,WSN)应用领域愈加广泛。然而,这类设备运行所需要的能量往往来自于预先满电的电池,由于成本、尺寸和应用场景的限制,很难给设备电池充电或更换电池。因此从环境中的射频信号获取能量给各设备充能,同时实现安全接收信息数据,这是可行且高效的解决方案。这也为无线携能通信的概念奠定了基础。With the continuous evolution of communication technology, green communication and the Internet of Everything have become the main development trends of future communication networks. Massive devices, including mobile phones, sensors, and household appliances, are connected to communication networks for functions such as remote control, data collection, and data exchange. Low-cost and low-power devices are also typical devices for future communications, and such devices will be more and more widely used in the fields of Internet of Things (Internet of Things, IoT) and Wireless Sensor Network (WSN). However, the energy required for the operation of such devices often comes from pre-charged batteries, and it is difficult to charge or replace the device batteries due to cost, size, and application scenarios. Therefore, it is a feasible and efficient solution to obtain energy from the radio frequency signal in the environment to charge each device, and at the same time realize the safe reception of information and data. This also laid the foundation for the concept of wireless energy-carrying communication.

无线携能通信(Simultaneous Wireless Information and Power Transfer,SWIPT)系统中,通常采用EH路径进行能量采集,ID路径进行信息传输。目前,接收到的射频信号功率可用两种方式实现能量收集和信息传输,分别为:时间分配方案和功率分割方案,时间分配(Time Sharing,TS)方案以交替时间间隔方式分别向每条路径提供射频信号,功率分割(PowerSplitter,PS)方案则是根据功率分割比同时向各路径提供射频信号。In a Simultaneous Wireless Information and Power Transfer (SWIPT) system, an EH path is usually used for energy collection, and an ID path is used for information transmission. At present, the received RF signal power can be used for energy collection and information transmission in two ways: time allocation scheme and power division scheme. Time sharing (TS) scheme provides each path with alternate time intervals. For the radio frequency signal, the Power Splitter (PS) scheme simultaneously provides the radio frequency signal to each path according to the power split ratio.

但是现有的无线携能通信系统的传输方案中均是基于理想信道假设,即发射端可得到准确的信道状态信息,再通过优化功率分割实现能量、信息的有效收集。但在实际中,由于信道估计误差等因素的影响,同时无线信道的状态是时刻变化的,完全能够与信道匹配的理想信道状态信息(Channel State Information,CSI)的提取是非常困难的。所以在目前的方案中存在以下问题:一是普遍忽略了更符合实际的非理想信道对SWIPT系统的影响;二是现有的方案中普遍忽略了各种干扰对多用户的影响。这会导致SWIPT系统的服务质量QoS约束条件放宽,从而使得SWIPT系统的传输可能无法达到预期的效果。However, the transmission schemes of the existing wireless energy-carrying communication systems are all based on the ideal channel assumption, that is, the transmitter can obtain accurate channel state information, and then realize the effective collection of energy and information by optimizing power division. However, in practice, due to the influence of channel estimation errors and other factors, and the state of the wireless channel changes from time to time, it is very difficult to extract the ideal channel state information (Channel State Information, CSI) that can fully match the channel. Therefore, there are the following problems in the current scheme: First, the influence of more realistic non-ideal channels on the SWIPT system is generally ignored; second, the existing scheme generally ignores the influence of various interferences on multi-users. This will lead to relaxation of the QoS constraints of the quality of service of the SWIPT system, so that the transmission of the SWIPT system may not achieve the expected effect.

发明内容SUMMARY OF THE INVENTION

为解决上述技术问题,本发明提供一种非理想信道下时间反演无线携能通信系统的传输方法。In order to solve the above technical problems, the present invention provides a transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel.

本发明采用的技术方案是:The technical scheme adopted in the present invention is:

一种非理想信道下时间反演无线携能通信系统的传输方法,包括如下步骤:A transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel, comprising the following steps:

S1:用户终端向基站发送探测导频信号;S1: the user terminal sends a sounding pilot signal to the base station;

S2:所述基站根据接收到的信号估计对应用户终端的信道增益;S2: the base station estimates the channel gain of the corresponding user terminal according to the received signal;

S3:通过求解第一问题模型计算出基站在非理想信道下向第k个用户终端发送数据的最优发射功率pk,所述第一问题模型为:S3: Calculate the optimal transmit power pk of the base station sending data to the kth user terminal under the non-ideal channel by solving the first problem model, the first problem model is:

Figure BDA0002292913470000021
Figure BDA0002292913470000021

s.t.C1:Rk≥Rk * stC1:R k ≥R k *

C2:Ek≥Ek * C2:E k ≥E k *

Figure BDA0002292913470000022
Figure BDA0002292913470000022

C4:0≤pk≤ppeak C4: 0≤p k ≤p peak

Figure BDA0002292913470000023
Figure BDA0002292913470000024
且k∈{1,2…K}
Figure BDA0002292913470000023
Figure BDA0002292913470000024
and k∈{1,2…K}

其中,P1表示所述第一问题模型中待求解的目标函数,C1、C2、C3、C4以及C5表示所述第一问题模型中目标函数P1的约束条件,pk表示基站在第k时隙向对应的第k个用户终端发送数据的发射功率,αk表示第k个用户终端对应的非负模型变量,K表示用户终端数,

Figure BDA0002292913470000031
表示非理想信道下第k个用户终端的信道增益,
Figure BDA0002292913470000032
表示第k个用户终端的天线噪声,
Figure BDA0002292913470000033
表示系统噪声,Rk *表示预先设置的第k个时隙对应的用户终端的最小信息可达速率阈值,Ek *表示预先设置的第k个时隙对应的用户终端的最低接收能量阈值,Rk表示第k个时隙对应的用户终端在非理想信道下的实际信息可达速率,Ek表示第k个时隙对应的用户终端在非理想信道下的实际能量值,P表示预先设置的基站的最大发射功率,
Figure BDA0002292913470000034
表示第k个用户终端在非理性信道下接收到的符号间干扰功率,
Figure BDA0002292913470000035
表示k个用户终端在非理想信道下接收到的码间干扰功率,表示通信过程中偏向信息传输的加权因子,ppeak表示每个时隙的峰值功率,0<β<1,Rk、Ek
Figure BDA0002292913470000036
以及
Figure BDA0002292913470000037
均通过含有信道估计误差的表达式计算得到;Among them, P1 represents the objective function to be solved in the first problem model, C1, C2, C3, C4 and C5 represent the constraints of the objective function P1 in the first problem model, and p k represents the base station in the kth time slot is the transmit power of sending data to the corresponding kth user terminal, α k represents the non-negative model variable corresponding to the kth user terminal, K represents the number of user terminals,
Figure BDA0002292913470000031
represents the channel gain of the kth user terminal under the non-ideal channel,
Figure BDA0002292913470000032
represents the antenna noise of the kth user terminal,
Figure BDA0002292913470000033
represents the system noise, R k * represents the preset minimum information reachable rate threshold of the user terminal corresponding to the k th time slot, E k * represents the preset minimum receiving energy threshold of the user terminal corresponding to the k th time slot, R k represents the actual information reachable rate of the user terminal corresponding to the kth time slot in the non-ideal channel, E k represents the actual energy value of the user terminal corresponding to the kth time slot in the non-ideal channel, P represents the preset the maximum transmit power of the base station,
Figure BDA0002292913470000034
represents the intersymbol interference power received by the kth user terminal under the irrational channel,
Figure BDA0002292913470000035
Represents the inter-symbol interference power received by k user terminals in a non-ideal channel, represents the weighting factor biased towards information transmission in the communication process, p peak represents the peak power of each time slot, 0<β<1, R k , E k ,
Figure BDA0002292913470000036
as well as
Figure BDA0002292913470000037
are calculated by expressions containing channel estimation errors;

S4:基站将第k个时隙的发射功率设置为pk,并基于该发射功率发送经过时间反转处理的信息数据给对应的用户终端;S4: the base station sets the transmit power of the kth time slot to p k , and sends the information data subjected to time reversal processing to the corresponding user terminal based on the transmit power;

S5:用户终端在接收到信息数据后,将β倍的信号功率传递给信号解码器,将1-β倍的信号功率传递给能量收集器。S5: After receiving the information data, the user terminal transmits β times the signal power to the signal decoder, and transmits 1-β times the signal power to the energy collector.

进一步地,步骤S3中通过将所述第一问题模型转化成第二问题模型来计算基站第k个时隙的发射功率pk,所述第二问题模型为:Further, in step S3, the transmit power pk of the kth time slot of the base station is calculated by converting the first problem model into a second problem model, and the second problem model is:

Figure BDA0002292913470000038
Figure BDA0002292913470000038

s.t.C6:Rk≥Rk * stC6:R k ≥R k *

C7:Ek≥Ek * C7:E k ≥E k *

Figure BDA0002292913470000039
Figure BDA0002292913470000039

C9:0≤pk≤ppeak C9: 0≤p k ≤p peak

Figure BDA00022929134700000310
且k∈{1,2…K}
Figure BDA00022929134700000310
and k∈{1,2…K}

其中,P2表示所述第二问题模型中待求解的目标函数,C6、C7、C8、C9以及C10表示所述第二问题模型中目标函数P2的约束条件,qk表示非负的中间变量,qk=αkpkWherein, P2 represents the objective function to be solved in the second problem model, C6, C7, C8, C9 and C10 represent the constraints of the objective function P2 in the second problem model, q k represents a non-negative intermediate variable, q kk p k .

进一步地,步骤S3中通过CVX凸优化工具箱求解目标函数P2从而得到目标解pkFurther, in step S3, the objective function P2 is solved by using the CVX convex optimization toolbox to obtain the objective solution p k .

进一步地,所述用户终端为单天线用户终端。Further, the user terminal is a single-antenna user terminal.

进一步地,

Figure BDA0002292913470000041
Rk以及Ek分别通过以下公式计算得到:further,
Figure BDA0002292913470000041
R k and E k are calculated by the following formulas respectively:

Figure BDA0002292913470000042
Figure BDA0002292913470000042

Figure BDA0002292913470000043
Figure BDA0002292913470000043

Figure BDA0002292913470000044
Figure BDA0002292913470000044

Figure BDA0002292913470000045
Figure BDA0002292913470000045

Figure BDA0002292913470000046
Figure BDA0002292913470000046

其中,

Figure BDA0002292913470000047
in,
Figure BDA0002292913470000047

Figure BDA0002292913470000051
Figure BDA0002292913470000051

Figure BDA0002292913470000052
表示第k个用户终端在非理想信道下的接收功率,
Figure BDA0002292913470000053
表示能量转化效率,L表示多径总条数,M表示基站的发射天线总数,ψ表示预先设置的信道误差,p'k表示用户终端k的发射功率,
Figure BDA0002292913470000054
表示基站天线x与用户终端y之间的噪声功率,
Figure BDA0002292913470000055
表示基站天线i与用户终端j之间第r条多径的高斯白噪声,
Figure BDA0002292913470000052
represents the received power of the kth user terminal under the non-ideal channel,
Figure BDA0002292913470000053
represents the energy conversion efficiency, L represents the total number of multipaths, M represents the total number of transmit antennas of the base station, ψ represents the preset channel error, p' k represents the transmit power of the user terminal k,
Figure BDA0002292913470000054
represents the noise power between the base station antenna x and the user terminal y,
Figure BDA0002292913470000055
represents the white Gaussian noise of the r-th multipath between the base station antenna i and the user terminal j,

其中,(x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r)∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')},m,m'∈{1,2…M},k,k'∈{1,2…K},l,l'∈{1,2…L},

Figure BDA0002292913470000056
表示基站天线m与基站天线m'之间的相关性矩阵,
Figure BDA0002292913470000057
表示用户终端k的传输天线与用户终端k'的传输天线之间的相关性矩阵,(RU)kk'表示用户终端k与用户终端k'之间的相关性矩阵。where, (x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r) ∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')}, m,m'∈{1 ,2...M}, k,k'∈{1,2...K}, l,l'∈{1,2...L},
Figure BDA0002292913470000056
represents the correlation matrix between the base station antenna m and the base station antenna m',
Figure BDA0002292913470000057
represents the correlation matrix between the transmission antenna of the user terminal k and the transmission antenna of the user terminal k', and (R U ) kk' represents the correlation matrix between the user terminal k and the user terminal k'.

本发明提供的非理想信道下时间反演无线携能通信系统的传输方法充分考虑了时间反演特性对非理想信道状态信息下多用户多发单收无线携能通信系统可靠性和有效性的影响,利用时间反演时空聚焦特性增强了系统的效性,通过含有信道估计误差的能量-速率闭合表达式,在非理想CSI的条件下规划优化问题,通过本发明提供的方法可以在保证一定信息可达速率条件下提升能量收集效率,相比于传统MISO-SWIPT,提高了系统的可靠性和有效性。The transmission method of the time-reversed wireless energy-carrying communication system under the non-ideal channel provided by the present invention fully considers the influence of the time-reversal characteristics on the reliability and effectiveness of the multi-user multi-transmit single-receive wireless energy-carrying communication system under the non-ideal channel state information , the effectiveness of the system is enhanced by using the time-reversal space-time focusing characteristics, and the optimization problem is planned under the condition of non-ideal CSI through the energy-rate closed expression containing the channel estimation error. The method provided by the present invention can guarantee certain information. Compared with traditional MISO-SWIPT, the energy harvesting efficiency is improved under the condition of achievable rate, which improves the reliability and effectiveness of the system.

附图说明Description of drawings

下面将结合附图及实施例对本发明作进一步说明,附图中:The present invention will be further described below in conjunction with the accompanying drawings and embodiments, in which:

图1为本实施例提供的非理想信道下时间反演无线携能通信系统的传输方法的流程示意图;1 is a schematic flowchart of a transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel provided by the present embodiment;

图2为本实施例所提供的时间反演无线携能通信系统的系统框图;2 is a system block diagram of a time-reversed wireless energy-carrying communication system provided by the present embodiment;

图3为实验过程中SINR随SNR变化的关系图。;Figure 3 is a graph showing the relationship between SINR and SNR during the experiment. ;

图4为实验过程中用户终端的能量收集器接收到的能量随用基站发射功率变化的示意图;Fig. 4 is the schematic diagram that the energy received by the energy collector of the user terminal varies with the transmission power of the base station during the experiment;

图5为实验过程中用户终端的能量收集器接收到的能量随用户终端的可达信息速率变化的示意图。FIG. 5 is a schematic diagram illustrating the variation of the energy received by the energy harvester of the user terminal with the achievable information rate of the user terminal during the experiment.

具体实施方式Detailed ways

为了使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the technical problems, technical solutions and advantages to be solved by the present invention clearer, the following will be described in detail with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, not for The invention is limited.

本实施例提供了一种非理想信道下时间反演无线携能通信系统的传输方法,请参见图1所示,包括如下步骤:This embodiment provides a transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel, as shown in FIG. 1, including the following steps:

S1:用户终端向基站发送探测导频信号。S1: The user terminal sends a sounding pilot signal to the base station.

S2:基站根据接收到的信号估计对应用户终端的信道增益。S2: The base station estimates the channel gain of the corresponding user terminal according to the received signal.

为了保证用户终端在一定的信息可达速率条件下最大化能量获取,预先设置了第一问题模型。In order to ensure that the user terminal maximizes energy acquisition under the condition of a certain information reachability rate, a first problem model is preset.

S3:通过求解第一问题模型计算出基站在非理想信道下向第k个用户终端发送数据的最优发射功率pk,第一问题模型为:S3: Calculate the optimal transmit power pk of the base station sending data to the kth user terminal under the non-ideal channel by solving the first problem model. The first problem model is:

Figure BDA0002292913470000061
Figure BDA0002292913470000061

s.t.C1:Rk≥Rk * stC1:R k ≥R k *

C2:Ek≥Ek * C2:E k ≥E k *

Figure BDA0002292913470000062
Figure BDA0002292913470000062

C4:0≤pk≤ppeak C4: 0≤p k ≤p peak

Figure BDA0002292913470000063
且k∈{1,2…K}
Figure BDA0002292913470000063
and k∈{1,2…K}

其中,P1表示第一问题模型中待求解的目标函数,C1、C2、C3、C4以及C5表示第一问题模型中目标函数P1的约束条件,pk表示基站在第k时隙向对应的第k个用户终端发送数据的发射功率,αk表示第k个用户终端对应的非负模型变量,K表示用户终端数,

Figure BDA0002292913470000064
表示非理想信道下第k个用户终端的信道增益,
Figure BDA0002292913470000071
表示第k个用户终端的天线噪声,
Figure BDA0002292913470000072
表示系统噪声,Rk *表示预先设置的第k个时隙对应的用户终端的最小信息可达速率阈值,Ek *表示预先设置的第k个时隙对应的用户终端的最低接收能量阈值,Rk表示第k个时隙对应的用户终端在非理想信道下的实际信息可达速率,Ek表示第k个时隙对应的用户终端在非理想信道下的实际能量值,P表示预先设置的基站的最大发射功率,
Figure BDA0002292913470000073
表示第k个用户终端在非理性信道下接收到的符号间干扰功率,
Figure BDA0002292913470000074
表示k个用户终端在非理想信道下接收到的码间干扰功率,表示通信过程中偏向信息传输的加权因子,ppeak表示每个时隙的峰值功率,0<β<1,Rk、Ek
Figure BDA0002292913470000075
以及
Figure BDA0002292913470000076
均通过含有信道估计误差的表达式计算得到。Among them, P1 represents the objective function to be solved in the first problem model, C1, C2, C3, C4 and C5 represent the constraints of the objective function P1 in the first problem model, p k represents the base station in the kth time slot corresponding to the The transmit power of data sent by k user terminals, α k represents the non-negative model variable corresponding to the kth user terminal, K represents the number of user terminals,
Figure BDA0002292913470000064
represents the channel gain of the kth user terminal under the non-ideal channel,
Figure BDA0002292913470000071
represents the antenna noise of the kth user terminal,
Figure BDA0002292913470000072
represents the system noise, R k * represents the preset minimum information reachable rate threshold of the user terminal corresponding to the k th time slot, E k * represents the preset minimum receiving energy threshold of the user terminal corresponding to the k th time slot, R k represents the actual information reachable rate of the user terminal corresponding to the kth time slot in the non-ideal channel, E k represents the actual energy value of the user terminal corresponding to the kth time slot in the non-ideal channel, P represents the preset the maximum transmit power of the base station,
Figure BDA0002292913470000073
represents the intersymbol interference power received by the kth user terminal under the irrational channel,
Figure BDA0002292913470000074
Represents the inter-symbol interference power received by k user terminals in a non-ideal channel, represents the weighting factor biased towards information transmission in the communication process, p peak represents the peak power of each time slot, 0<β<1, R k , E k ,
Figure BDA0002292913470000075
as well as
Figure BDA0002292913470000076
Both are calculated by expressions containing channel estimation errors.

第一问题模型为非凸优化问题,在计算过程中,可以通过引入非负变量qk来保证其连续性将非凸优化问题转化为凸优化问题,其中,qk=αkpk(k=1,···,k),也即是在步骤S3中可以将第一问题模型转化成第二问题模型来计算基站第k个时隙的发射功率pk,其中,第二问题模型为:The first problem model is a non-convex optimization problem. In the calculation process, the non-convex optimization problem can be transformed into a convex optimization problem by introducing a non-negative variable q k to ensure its continuity, where q kk p k (k =1,...,k), that is, in step S3, the first problem model can be transformed into a second problem model to calculate the transmit power p k of the kth time slot of the base station, where the second problem model is :

Figure BDA0002292913470000077
Figure BDA0002292913470000077

s.t.C6:Rk≥Rk * stC6:R k ≥R k *

C7:Ek≥Ek * C7:E k ≥E k *

Figure BDA0002292913470000078
Figure BDA0002292913470000078

C9:0≤pk≤ppeak C9: 0≤p k ≤p peak

Figure BDA0002292913470000079
且k∈{1,2…K}
Figure BDA0002292913470000079
and k∈{1,2…K}

其中,P2表示第二问题模型中待求解的目标函数,C6、C7、C8、C9以及C10表示第二问题模型中目标函数P2的约束条件,qk表示非负的中间变量,qk=αkpkAmong them, P2 represents the objective function to be solved in the second problem model, C6, C7, C8, C9 and C10 represent the constraints of the objective function P2 in the second problem model, q k represents a non-negative intermediate variable, q kk p k .

步骤S3中可以通过CVX凸优化工具箱求解目标函数P2从而得到目标解pk。应当说明的是,本实施例中的步骤S3可以由基站执行,也即可以由基站计算出pk,但是在其他的实施例中pk也可以由另外专门设置的计算终端计算得到,计算终端在计算出pk后将pk发送给基站以供基站设置发射功率。In step S3, the objective function P2 can be solved by the CVX convex optimization toolbox to obtain the objective solution p k . It should be noted that step S3 in this embodiment may be performed by the base station, that is, the base station may calculate p k , but in other embodiments, p k may also be calculated by a specially set computing terminal, and the computing terminal After the pk is calculated, the pk is sent to the base station for the base station to set the transmit power.

S4:基站将第k个时隙的发射功率设置为pk,并基于该发射功率发送经过时间反转处理的信息数据给对应的用户终端。S4: The base station sets the transmit power of the kth time slot to p k , and sends the information data subjected to time reversal processing to the corresponding user terminal based on the transmit power.

由于TR(Time Reversal,时间反演)的时空聚焦特性,基站的发送信号能在很短的时间间隔内汇集大部分有用信号功率,并在目标处形成聚。Due to the space-time focusing characteristic of TR (Time Reversal, time inversion), the transmitted signal of the base station can gather most of the useful signal power in a very short time interval, and form a cluster at the target.

S5:用户终端在接收到信息数据后,将β倍的信号功率传递给信号解码器,将1-β倍的信号功率传递给能量收集器。S5: After receiving the information data, the user terminal transmits β times the signal power to the signal decoder, and transmits 1-β times the signal power to the energy collector.

本实施例中的用户终端可以采用PS接收机。The user terminal in this embodiment may adopt a PS receiver.

本实施例所提供的时间反演无线携能通信系统的系统框图可以参见图2所示,图2中用户终端在接收到基站发送的经TR调制的信号功率后,将信号功率的β部分传递给信号接收机(也即信号解码器),将信号功率的1-β部分传递给能量接收机(也即能量收集器)。For the system block diagram of the time-reversed wireless energy-carrying communication system provided in this embodiment, see FIG. 2 . In FIG. 2 , after receiving the TR-modulated signal power sent by the base station, the user terminal transmits the β part of the signal power. To the signal receiver (ie, the signal decoder), the 1-beta portion of the signal power is passed to the energy receiver (ie, the energy harvester).

应当说明的是,优选的,本实施例中的用户终端为单天线用户终端。It should be noted that, preferably, the user terminal in this embodiment is a single-antenna user terminal.

假设本实施例中的发送端(基站)采用M根发送天线,合法接收端(用户终端侧)采用K根接收天线,也即系统中具有K个用户终端,每一用户终端均为单天线用户终端,发送端发射天线m∈(1,M)与合法接收端接收天线k∈(1,K)之间的信道冲击响应(channelimpulse response,CIR)可通过信道估计获得,本实施例中的CIR可记为:It is assumed that the transmitting end (base station) in this embodiment adopts M transmitting antennas, and the legitimate receiving end (user terminal side) adopts K receiving antennas, that is, there are K user terminals in the system, and each user terminal is a single-antenna user In the terminal, the channel impulse response (CIR) between the transmitting antenna m∈(1,M) of the transmitting end and the receiving antenna k∈(1,K) of the legitimate receiving end can be obtained by channel estimation. The CIR in this embodiment Can be written as:

Figure BDA0002292913470000081
Figure BDA0002292913470000081

其中,

Figure BDA0002292913470000082
Figure BDA0002292913470000083
分别是第l个抽头的振幅和延迟。CSI在时域内离散的矢量hmk∈CL ×1,特别的E[hmk[l]]=0,
Figure BDA0002292913470000084
E[·]表示信号期望。in,
Figure BDA0002292913470000082
and
Figure BDA0002292913470000083
are the amplitude and delay of the lth tap, respectively. CSI is a discrete vector h mkC L ×1 in the time domain, especially E[h mk [l]]=0,
Figure BDA0002292913470000084
E[·] represents the signal expectation.

单天线用户终端接收到的信号可以写作:The signal received by a single-antenna user terminal can be written as:

Figure BDA0002292913470000091
Figure BDA0002292913470000091

对应上式中各部分详解,其中sk和sk'分别代表第k个用户和第k'个用户的所要发送的符号(E(|s|2)=1)。将多径信道模型化为抽头延迟线模型,即假设hmk∈CL×1是第m根发射天线和第k个用户的目标天线之间的信道状态矩阵E(|hmk[l]|)=0,

Figure BDA0002292913470000092
hmk∈CL×1是引入系统的天线噪声
Figure BDA0002292913470000093
又由于hmk是属于复数的向量集合,则对复信道hm[l]进行复共轭和时序反转操作可以获得
Figure BDA0002292913470000094
因此,gm每个抽头的值可以写作:Corresponding to the detailed explanation of each part in the above formula, where sk and sk' represent the symbols to be sent by the kth user and the k'th user respectively (E(|s| 2 )=1). The multipath channel is modeled as a tapped delay line model, that is, it is assumed that h mkC L×1 is the channel state matrix E(|h mk [l]| )=0,
Figure BDA0002292913470000092
h mkC L×1 is the antenna noise introduced into the system
Figure BDA0002292913470000093
And since h mk is a vector set of complex numbers, the complex conjugation and timing inversion operations on the complex channel h m [l] can be obtained.
Figure BDA0002292913470000094
Therefore, the value of each tap of g m can be written as:

Figure BDA0002292913470000095
Figure BDA0002292913470000095

为使现在非理想CSI条件下达到目标能量收集效果,本实施例所提供的方法利用TR时空聚焦特性弥补SWIPT技术劣势。In order to achieve the target energy harvesting effect under the current non-ideal CSI conditions, the method provided in this embodiment makes use of the TR spatiotemporal focusing characteristic to make up for the technical disadvantage of SWIPT.

本实施例所提供的无线携能通信系统由具有M根发射天线的基站和N根单天线用户终端组成。在多径信道中,每个信道冲击响应的最大长度为L。在信息发送阶段,发送方首先经过时间反演镜,完成信号在非理想CSI条件下的时间反演过程,以下对上述过程展开分析:The wireless energy-carrying communication system provided in this embodiment is composed of a base station having M transmitting antennas and N single-antenna user terminals. In a multipath channel, the maximum length of the impulse response of each channel is L. In the information transmission stage, the sender first passes through the time inversion mirror to complete the time inversion process of the signal under non-ideal CSI conditions. The above process is analyzed as follows:

由于在实际应用中,真正的信道是一个未知的参数,因此将信道误差估计的影响进行如下建模:Since the real channel is an unknown parameter in practical applications, the effect of channel error estimation is modeled as follows:

Figure BDA0002292913470000096
Figure BDA0002292913470000096

Figure BDA0002292913470000097
和emk∈CL×1分别表示估计通道和误差向量,它们相互独立分布且带有一个非负的误差因子ψ,可以发现具有以下特征:
Figure BDA0002292913470000097
and emkC L×1 represent the estimated channel and the error vector, respectively, which are distributed independently of each other and have a non-negative error factor ψ, which can be found to have the following characteristics:

E[|emk[l]|2]=ψE[||hmk[l]|2]E[|e mk [l]| 2 ]=ψE[||h mk [l]| 2 ]

Figure BDA0002292913470000098
Figure BDA0002292913470000098

基于非理想信道的TR技术,结合上面两式得到

Figure BDA0002292913470000099
表示的是含有发送功率p'k的非理想信道TR预滤波向量,每个抽头的值定义为:Based on the TR technology of the non-ideal channel, combining the above two formulas, we can get
Figure BDA0002292913470000099
represents the non-ideal channel TR pre-filter vector with transmit power p' k , and the value of each tap is defined as:

Figure BDA00022929134700000910
Figure BDA00022929134700000910

Figure BDA0002292913470000101
是hmk[l]的在非理想信道条件下的共轭和时序反转操作。
Figure BDA0002292913470000101
is the conjugation and timing reversal operation of h mk [l] under nonideal channel conditions.

由此,第k个用户终端接收到的信号为Thus, the signal received by the kth user terminal is

Figure BDA0002292913470000102
Figure BDA0002292913470000102

本系统的SINR表达式为:The SINR expression of this system is:

Figure BDA0002292913470000103
Figure BDA0002292913470000103

各用户终端的接收功率可表示为:The received power of each user terminal can be expressed as:

Figure BDA0002292913470000104
Figure BDA0002292913470000104

符号间干扰的功率可表示为:The power of ISI can be expressed as:

Figure BDA0002292913470000105
Figure BDA0002292913470000105

用户间干扰的功率可表示为:The power of inter-user interference can be expressed as:

Figure BDA0002292913470000106
Figure BDA0002292913470000106

可得到SINR(Signal to Interference plus Noise Ratio,信号与干扰加噪声比)具体推导过程如下The specific derivation process of SINR (Signal to Interference plus Noise Ratio, signal to interference plus noise ratio) can be obtained as follows

Figure BDA0002292913470000107
Figure BDA0002292913470000107

本实施例所提供的方法充分考虑了信道误差对传输的影响,在推导过程中引入了信道误差,进一步推导各个表达式,各部分具体推导如下:The method provided by this embodiment fully considers the influence of channel error on transmission, introduces channel error in the derivation process, and further derives each expression. The specific derivation of each part is as follows:

Figure BDA0002292913470000108
Figure BDA0002292913470000108

Figure BDA0002292913470000109
Figure BDA0002292913470000109

Figure BDA0002292913470000111
Figure BDA0002292913470000111

所以,在多径误差信道和上面给出的相关性的假设下,本实施例中的

Figure BDA0002292913470000112
Figure BDA0002292913470000113
Rk以及Ek可以分别通过以下公式计算得到:Therefore, under the assumption of the multipath error channel and the correlation given above, the
Figure BDA0002292913470000112
Figure BDA0002292913470000113
R k and E k can be calculated by the following formulas respectively:

Figure BDA0002292913470000114
Figure BDA0002292913470000114

Figure BDA0002292913470000115
Figure BDA0002292913470000115

Figure BDA0002292913470000116
Figure BDA0002292913470000116

Figure BDA0002292913470000117
Figure BDA0002292913470000117

Figure BDA0002292913470000118
Figure BDA0002292913470000118

其中,

Figure BDA0002292913470000119
Figure BDA00022929134700001110
in,
Figure BDA0002292913470000119
Figure BDA00022929134700001110

Figure BDA00022929134700001111
表示第k个用户终端在非理想信道下的接收功率,
Figure BDA00022929134700001112
表示能量转化效率,L表示多径总条数,M表示基站的发射天线总数,ψ表示预先设置的信道误差,p'k表示用户终端k的发射功率,
Figure BDA0002292913470000121
表示基站天线x与用户终端y之间的噪声功率,
Figure BDA0002292913470000122
表示基站天线i与用户终端j之间第r条多径的高斯白噪声,
Figure BDA00022929134700001111
represents the received power of the kth user terminal under the non-ideal channel,
Figure BDA00022929134700001112
represents the energy conversion efficiency, L represents the total number of multipaths, M represents the total number of transmit antennas of the base station, ψ represents the preset channel error, p' k represents the transmit power of the user terminal k,
Figure BDA0002292913470000121
represents the noise power between the base station antenna x and the user terminal y,
Figure BDA0002292913470000122
represents the white Gaussian noise of the r-th multipath between the base station antenna i and the user terminal j,

其中,(x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r)∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')},m,m'∈{1,2…M},k,k'∈{1,2…K},l,l'∈{1,2…L},

Figure BDA0002292913470000123
表示基站天线m与基站天线m'之间的相关性矩阵,
Figure BDA0002292913470000124
表示用户终端k的传输天线与用户终端k'的传输天线之间的相关性矩阵,(RU)kk'表示用户终端k与用户终端k'之间的相关性矩阵。where, (x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r) ∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')}, m,m'∈{1 ,2...M}, k,k'∈{1,2...K}, l,l'∈{1,2...L},
Figure BDA0002292913470000123
represents the correlation matrix between the base station antenna m and the base station antenna m',
Figure BDA0002292913470000124
represents the correlation matrix between the transmission antenna of the user terminal k and the transmission antenna of the user terminal k', and (R U ) kk' represents the correlation matrix between the user terminal k and the user terminal k'.

以上计算

Figure BDA0002292913470000125
Rk以及Ek的计算公式都充分考虑了信道误差的影响,因此提高了系统的可靠性和有效性。above calculation
Figure BDA0002292913470000125
The calculation formulas of R k and E k fully consider the influence of the channel error, thus improving the reliability and effectiveness of the system.

为了验证本实施例所提供的方法的有效性,对不同条件下基于TR的无线携能通信系统的传输方法进行理论推导值的仿真验证,并分析系统性能在这些条件下的变化趋势,最后与传统无线携能通信系统传输方案进行性能对比,研究加入TR及凸优化后的鲁棒性方案对系统性能的影响。实验过程中的参数为:当m=m'时,(RT)mm'=1,当m≠m'时,(RT)mm'=ρT,ρT表示第m根天线和第m'根天线之间的相关程度,其它通用设置为抽头数L=110,基站的平均发射功率p=1W。In order to verify the effectiveness of the method provided in this embodiment, the theoretical derivation value of the transmission method of the TR-based wireless energy-carrying communication system under different conditions is simulated and verified, and the change trend of the system performance under these conditions is analyzed. The performance of the traditional wireless energy-carrying communication system transmission scheme is compared, and the impact of the robust scheme after adding TR and convex optimization on the system performance is studied. The parameters during the experiment are: when m=m', (R T ) mm' = 1, when m≠m', (R T ) mm' = ρ T , ρ T represents the mth antenna and the mth antenna 'The degree of correlation between the antennas, other common settings are the number of taps L=110, and the average transmit power of the base station p=1W.

进行了相关实验仿真。The related experimental simulations were carried out.

图3为在时间反演过程中理想信道和非理想信道这两种条件下SINR随SNR变化的关系图,图3中的信道误差因子为ψ=0.5,也就是说,在实际应用中的非理想信道条件下,会给系统带来一定程度上的削弱。Figure 3 shows the relationship between SINR and SNR under two conditions of ideal channel and non-ideal channel in the time inversion process. The channel error factor in Figure 3 is ψ=0.5, that is, the non-ideal channel in practical applications Under ideal channel conditions, the system will be weakened to a certain extent.

图4是在基站发射天线相关性为0(天线之间相互独立),传统的理想信道条件下的无线携能通信系统传输方法和本实施例所提供的非理想信道条件下基于TR的无线携能通信系统传输方法在实验过程中的仿真图,示出了用户终端的能量收集器接收到的能量随用户终端的发送功率变化的示意图,从图中可以看到仿真值和理论值相符,验证了本实施例中理论推导值在功率分割比例β、基站天线数量、基站发射功率变化时的正确性。其中具体的变化趋势有:基站的发送功率增大时,用户终端的能量接收器接收到的能量也增大,但通过本实施例提供的方法,在任一发送功率下,能量效率都优于传统方案的能量效率,证实了本实施例所提供方法的有效性。Fig. 4 shows the transmission method of the wireless energy-carrying communication system under the traditional ideal channel condition and the TR-based wireless portability under the non-ideal channel condition provided by the present embodiment when the correlation between the transmit antennas of the base station is 0 (the antennas are independent of each other). The simulation diagram of the transmission method of the energy communication system in the experimental process shows the schematic diagram of the energy received by the energy harvester of the user terminal changing with the transmission power of the user terminal. It can be seen from the figure that the simulation value is consistent with the theoretical value, and the verification The correctness of the theoretically derived value in this embodiment when the power division ratio β, the number of base station antennas, and the base station transmit power change. The specific change trend is as follows: when the transmission power of the base station increases, the energy received by the energy receiver of the user terminal also increases, but with the method provided in this embodiment, under any transmission power, the energy efficiency is better than the traditional The energy efficiency of the scheme confirms the effectiveness of the method provided in this example.

图5是在基站发射天线相关性为0(相互独立)、功率分割比例、天线数量发生变化时,通过仿真得到的。图中有两种传输方案,一种是本实施例所提供的基于TR的无线携能通信系统传输方法(非理想信道鲁棒性方法),另一种是传统的理想信道无线携能通信系统传输方法,从图中可以看到,在功率分割比例较小时(可达信息速率较小),采用鲁棒性方案的基于TR的无线携能通信系统收集到的能量已经接近增加一定数量天线的传统方案,而在天线数量不变时,鲁棒性方案不论是可达信息速率还是收集到的能量都比通过传统理想信道无线携能通信系统传输方案收集到的能量大。当需要从可达信息速率和收集到的能量两个方面整体的考虑系统性能时,鲁棒性方案中平均速率-能量区的大小接近于一个正方形框,且其面积大于传统方案,此时系统既能满足一定可达信息速率,也能兼顾收集一定大小的能量。本实施例所提供的传输方法相较于传统的传输方法明显提升了能量收集效率。Figure 5 is obtained through simulation when the base station transmit antenna correlation is 0 (independent of each other), the power division ratio and the number of antennas change. There are two transmission schemes in the figure, one is the TR-based wireless energy-carrying communication system transmission method (non-ideal channel robustness method) provided in this embodiment, and the other is the traditional ideal channel wireless energy-carrying communication system Transmission method, it can be seen from the figure that when the power division ratio is small (the achievable information rate is small), the energy collected by the TR-based wireless energy-carrying communication system using the robust scheme is close to that of adding a certain number of antennas. The traditional scheme, and when the number of antennas is constant, the robust scheme can achieve greater information rate and collected energy than the traditional ideal channel wireless energy-carrying communication system transmission scheme. When the overall system performance needs to be considered from the two aspects of the achievable information rate and the collected energy, the size of the average rate-energy region in the robust scheme is close to a square box, and its area is larger than that of the traditional scheme. At this time, the system It can not only meet a certain reachable information rate, but also take into account the collection of a certain amount of energy. Compared with the traditional transmission method, the transmission method provided in this embodiment obviously improves the energy collection efficiency.

针对SWIPT-MISO功率分割和假设中天线相互独立的局限性,本文将非理想信道条件下TR与SWIPT-MISO相结合并利用凸优化算法优化其能量收集,提出了一种基于TR的新型SWIPT-MISO传输方法。通过该方法,可以采用PS接收机,在考虑天线具有相关性的情况下,引入克罗内克模型推导多径信道下的平均速率-能量区理论表达式。仿真结果显示,本实施例提供的方法在发射端天线数量、功率分割比例等条件发生变化时,得到的理论值与仿真值始终处于同一条曲线上,验证了本实施例所提供方法的推导过程的正确性,而在随后与传统SWIPT-MISO传输方案的性能对比中,充分显示了引入TR技术对系统性能的改善作用。这也从另一方面说明,将TR和SWIPT-MISO结合,可以实现两者优势互补,达到理想效果。上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。Aiming at the limitations of SWIPT-MISO power splitting and the assumption that the antennas are independent of each other, this paper combines TR with SWIPT-MISO under non-ideal channel conditions and uses a convex optimization algorithm to optimize its energy harvesting. A new TR-based SWIPT-MISO is proposed. MISO transfer method. With this method, PS receivers can be used, and the Kronecker model is introduced to derive the theoretical expression of the average rate-energy region under multipath channels, considering the correlation of the antennas. The simulation results show that when the method provided in this embodiment changes the number of transmitter antennas, the power split ratio and other conditions, the obtained theoretical value and the simulated value are always on the same curve, which verifies the derivation process of the method provided by this embodiment. In the subsequent performance comparison with the traditional SWIPT-MISO transmission scheme, it fully shows the improvement of the system performance by introducing the TR technology. This also shows from another aspect that the combination of TR and SWIPT-MISO can complement each other's advantages and achieve ideal results. The above-mentioned serial numbers of the embodiments of the present invention are only for description, and do not represent the advantages or disadvantages of the embodiments.

通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。From the description of the above embodiments, those skilled in the art can clearly understand that the method of the above embodiment can be implemented by means of software plus a necessary general hardware platform, and of course can also be implemented by hardware, but in many cases the former is better implementation. Based on this understanding, the technical solutions of the present invention can be embodied in the form of software products in essence or the parts that make contributions to the prior art, and the computer software products are stored in a storage medium (such as ROM/RAM, magnetic disk, CD), including several instructions to make a terminal (which may be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) execute the methods described in the various embodiments of the present invention.

上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。The embodiments of the present invention have been described above in conjunction with the accompanying drawings, but the present invention is not limited to the above-mentioned specific embodiments, which are merely illustrative rather than restrictive. Under the inspiration of the present invention, without departing from the scope of protection of the present invention and the claims, many forms can be made, which all belong to the protection of the present invention.

Claims (5)

1.一种非理想信道下时间反演无线携能通信系统的传输方法,其特征在于,包括如下步骤:1. a transmission method for time inversion wireless energy-carrying communication system under a non-ideal channel, is characterized in that, comprises the steps: S1:用户终端向基站发送探测导频信号;S1: the user terminal sends a sounding pilot signal to the base station; S2:所述基站根据接收到的信号估计对应用户终端的信道增益;S2: the base station estimates the channel gain of the corresponding user terminal according to the received signal; S3:通过求解第一问题模型计算出基站在非理想信道下向第k个用户终端发送数据的最优发射功率pk,所述第一问题模型为:S3: Calculate the optimal transmit power pk of the base station sending data to the kth user terminal under the non-ideal channel by solving the first problem model, the first problem model is: P1:
Figure FDA0002292913460000011
P1:
Figure FDA0002292913460000011
s.t.C1:Rk≥Rk * stC1:R k ≥R k * C2:Ek≥Ek * C2:E k ≥E k * C3:
Figure FDA0002292913460000012
C3:
Figure FDA0002292913460000012
C4:0≤pk≤ppeak C4: 0≤p k ≤p peak C5:
Figure FDA0002292913460000013
且k∈{1,2…K}
C5:
Figure FDA0002292913460000013
and k∈{1,2…K}
其中,P1表示所述第一问题模型中待求解的目标函数,C1、C2、C3、C4以及C5表示所述第一问题模型中目标函数P1的约束条件,pk表示基站在第k时隙向对应的第k个用户终端发送数据的发射功率,αk表示第k个用户终端对应的非负模型变量,K表示用户终端数,
Figure FDA0002292913460000014
表示非理想信道下第k个用户终端的信道增益,
Figure FDA0002292913460000015
表示第k个用户终端的天线噪声,
Figure FDA0002292913460000016
表示系统噪声,Rk *表示预先设置的第k个时隙对应的用户终端的最小信息可达速率阈值,Ek *表示预先设置的第k个时隙对应的用户终端的最低接收能量阈值,Rk表示第k个时隙对应的用户终端在非理想信道下的实际信息可达速率,Ek表示第k个时隙对应的用户终端在非理想信道下的实际能量值,P表示预先设置的基站的最大发射功率,
Figure FDA0002292913460000017
表示第k个用户终端在非理性信道下接收到的符号间干扰功率,
Figure FDA0002292913460000018
表示k个用户终端在非理想信道下接收到的码间干扰功率,表示通信过程中偏向信息传输的加权因子,ppeak表示每个时隙的峰值功率,0<β<1,Rk、Ek
Figure FDA0002292913460000021
以及
Figure FDA0002292913460000022
均通过含有信道估计误差的表达式计算得到;
Among them, P1 represents the objective function to be solved in the first problem model, C1, C2, C3, C4 and C5 represent the constraints of the objective function P1 in the first problem model, and p k represents the base station in the kth time slot is the transmit power of sending data to the corresponding kth user terminal, α k represents the non-negative model variable corresponding to the kth user terminal, K represents the number of user terminals,
Figure FDA0002292913460000014
represents the channel gain of the kth user terminal under the non-ideal channel,
Figure FDA0002292913460000015
represents the antenna noise of the kth user terminal,
Figure FDA0002292913460000016
represents the system noise, R k * represents the preset minimum information reachable rate threshold of the user terminal corresponding to the k th time slot, E k * represents the preset minimum receiving energy threshold of the user terminal corresponding to the k th time slot, R k represents the actual information reachable rate of the user terminal corresponding to the kth time slot in the non-ideal channel, E k represents the actual energy value of the user terminal corresponding to the kth time slot in the non-ideal channel, P represents the preset the maximum transmit power of the base station,
Figure FDA0002292913460000017
represents the intersymbol interference power received by the kth user terminal under the irrational channel,
Figure FDA0002292913460000018
Represents the inter-symbol interference power received by k user terminals in a non-ideal channel, represents the weighting factor biased towards information transmission in the communication process, p peak represents the peak power of each time slot, 0<β<1, R k , E k ,
Figure FDA0002292913460000021
as well as
Figure FDA0002292913460000022
are calculated by expressions containing channel estimation errors;
S4:基站将第k个时隙的发射功率设置为pk,并基于该发射功率发送经过时间反转处理的信息数据给对应的用户终端;S4: the base station sets the transmit power of the kth time slot to p k , and sends the information data subjected to time reversal processing to the corresponding user terminal based on the transmit power; S5:用户终端在接收到信息数据后,将β倍的信号功率传递给信号解码器,将1-β倍的信号功率传递给能量收集器。S5: After receiving the information data, the user terminal transmits β times the signal power to the signal decoder, and transmits 1-β times the signal power to the energy collector.
2.如权利要求1所述的非理想信道下时间反演无线携能通信系统的传输方法,其特征在于,步骤S3中通过将所述第一问题模型转化成第二问题模型来计算基站第k个时隙的发射功率pk,所述第二问题模型为:2. The transmission method of the time-reversed wireless energy-carrying communication system under a non-ideal channel as claimed in claim 1, wherein in step S3, the first problem model is converted into the second problem model to calculate the first problem of the base station. The transmit power pk of k time slots, the second problem model is: P2:
Figure FDA0002292913460000023
P2:
Figure FDA0002292913460000023
s.t.C6:Rk≥Rk * stC6:R k ≥R k * C7:Ek≥Ek * C7:E k ≥E k * C8:
Figure FDA0002292913460000024
C8:
Figure FDA0002292913460000024
C9:0≤pk≤ppeak C9: 0≤p k ≤p peak C10:
Figure FDA0002292913460000025
且k∈{1,2…K}
C10:
Figure FDA0002292913460000025
and k∈{1,2…K}
其中,P2表示所述第二问题模型中待求解的目标函数,C6、C7、C8、C9以及C10表示所述第二问题模型中目标函数P2的约束条件,qk表示非负的中间变量,qk=αkpkWherein, P2 represents the objective function to be solved in the second problem model, C6, C7, C8, C9 and C10 represent the constraints of the objective function P2 in the second problem model, q k represents a non-negative intermediate variable, q kk p k .
3.如权利要求2所述的非理想信道下时间反演无线携能通信系统的传输方法,其特征在于,步骤S3中通过CVX凸优化工具箱求解目标函数P2从而得到目标解pk3. The transmission method of the time-reversed wireless energy-carrying communication system under a non-ideal channel as claimed in claim 2, wherein in step S3, the objective function P2 is obtained by solving the objective function P2 through the CVX convex optimization toolbox to obtain the objective solution p k . 4.如权利要求1所述的非理想信道下时间反演无线携能通信系统的传输方法,其特征在于,所述用户终端为单天线用户终端。4 . The transmission method for a time-reversed wireless energy-carrying communication system under a non-ideal channel according to claim 1 , wherein the user terminal is a single-antenna user terminal. 5 . 5.如权利要求1-4任一项所述的非理想信道下时间反演无线携能通信系统的传输方法,其特征在于,
Figure FDA0002292913460000031
Rk以及Ek分别通过以下公式计算得到:
5. The transmission method of the time-reversed wireless energy-carrying communication system under a non-ideal channel according to any one of claims 1-4, wherein,
Figure FDA0002292913460000031
R k and E k are calculated by the following formulas respectively:
Figure FDA0002292913460000032
Figure FDA0002292913460000032
Figure FDA0002292913460000033
Figure FDA0002292913460000033
Figure FDA0002292913460000034
Figure FDA0002292913460000034
Figure FDA0002292913460000035
Figure FDA0002292913460000035
Figure FDA0002292913460000036
Figure FDA0002292913460000036
其中,
Figure FDA0002292913460000037
Figure FDA0002292913460000038
in,
Figure FDA0002292913460000037
Figure FDA0002292913460000038
Figure FDA0002292913460000039
表示第k个用户终端在非理想信道下的接收功率,
Figure FDA00022929134600000310
表示能量转化效率,L表示多径总条数,M表示基站的发射天线总数,ψ表示预先设置的信道误差,p'k表示用户终端k的发射功率,
Figure FDA0002292913460000041
表示基站天线x与用户终端y之间的噪声功率,
Figure FDA0002292913460000042
表示基站天线i与用户终端j之间第r条多径的高斯白噪声,
Figure FDA0002292913460000039
represents the received power of the kth user terminal under the non-ideal channel,
Figure FDA00022929134600000310
represents the energy conversion efficiency, L represents the total number of multipaths, M represents the total number of transmit antennas of the base station, ψ represents the preset channel error, p' k represents the transmit power of the user terminal k,
Figure FDA0002292913460000041
represents the noise power between the base station antenna x and the user terminal y,
Figure FDA0002292913460000042
represents the white Gaussian noise of the r-th multipath between the base station antenna i and the user terminal j,
其中,(x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r)∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')},m,m'∈{1,2…M},k,k'∈{1,2…K},l,l'∈{1,2…L},
Figure FDA0002292913460000043
表示基站天线m与基站天线m'之间的相关性矩阵,
Figure FDA0002292913460000044
表示用户终端k的传输天线与用户终端k'的传输天线之间的相关性矩阵,(RU)kk'表示用户终端k与用户终端k'之间的相关性矩阵。
where, (x,y)∈{(m,l),(m',l),(m,k+1-l),(m,L+1-l)},(i,j,r) ∈{(m,k,L-1-l),(m,k,l),(m,k',l),(m',k',l')}, m,m'∈{1 ,2...M}, k,k'∈{1,2...K}, l,l'∈{1,2...L},
Figure FDA0002292913460000043
represents the correlation matrix between the base station antenna m and the base station antenna m',
Figure FDA0002292913460000044
represents the correlation matrix between the transmission antenna of the user terminal k and the transmission antenna of the user terminal k', and (R U ) kk' represents the correlation matrix between the user terminal k and the user terminal k'.
CN201911188194.1A 2019-11-28 2019-11-28 Transmission method of time reversal wireless energy-carrying communication system under non-ideal channel Active CN111030795B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911188194.1A CN111030795B (en) 2019-11-28 2019-11-28 Transmission method of time reversal wireless energy-carrying communication system under non-ideal channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911188194.1A CN111030795B (en) 2019-11-28 2019-11-28 Transmission method of time reversal wireless energy-carrying communication system under non-ideal channel

Publications (2)

Publication Number Publication Date
CN111030795A true CN111030795A (en) 2020-04-17
CN111030795B CN111030795B (en) 2022-03-22

Family

ID=70202896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911188194.1A Active CN111030795B (en) 2019-11-28 2019-11-28 Transmission method of time reversal wireless energy-carrying communication system under non-ideal channel

Country Status (1)

Country Link
CN (1) CN111030795B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180026481A1 (en) * 2016-05-03 2018-01-25 Origin Wireless, Inc. Method, system, and apparatus for wireless power transmission based on power waveforming
CN109005551A (en) * 2018-07-10 2018-12-14 南京邮电大学 A kind of multi-user's NOMA downlink power distributing method of imperfect channel state information
CN109379154A (en) * 2018-10-11 2019-02-22 重庆邮电大学 A secure transmission scheme based on time inversion technology
CN109617590A (en) * 2019-01-11 2019-04-12 华南理工大学 Physical layer security communication method for multiple-input single-output wireless energy-carrying communication system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180026481A1 (en) * 2016-05-03 2018-01-25 Origin Wireless, Inc. Method, system, and apparatus for wireless power transmission based on power waveforming
CN109005551A (en) * 2018-07-10 2018-12-14 南京邮电大学 A kind of multi-user's NOMA downlink power distributing method of imperfect channel state information
CN109379154A (en) * 2018-10-11 2019-02-22 重庆邮电大学 A secure transmission scheme based on time inversion technology
CN109617590A (en) * 2019-01-11 2019-04-12 华南理工大学 Physical layer security communication method for multiple-input single-output wireless energy-carrying communication system

Also Published As

Publication number Publication date
CN111030795B (en) 2022-03-22

Similar Documents

Publication Publication Date Title
CN107483088B (en) Massive MIMO Robust Precoding Transmission Method
CN101771509B (en) Orthogonal network space-time coding method and relay transmission system
CN103067134B (en) Channel feedback method in ofdm system
CN114285444B (en) Power optimization method for large-scale de-cellular MIMO system
Han et al. Time-reversal division multiple access in multi-path channels
CN102291842B (en) Virtual MIMO pairing method taking user QoS into account
CN106533590A (en) Uplink channel quality measurement method based on receiver EVM
CN103716263A (en) Parallel-cancellation-based high-efficient pilot frequency method in large-scale multi-antenna system
CN106534008B (en) The power compensation MMSE equalization methods of wireless multi-path channels
CN114665930B (en) Downlink Blind Channel Estimation Method for Decellularized Massive MIMO Systems
CN118590107B (en) Terminal direct connection and non-cellular heterogeneous network access mode selection method
CN107154818A (en) Co-channel full duplex bi-directional relaying transmission method while based on single carrier frequency domain equalization
CN117560043B (en) Non-cellular network power control method based on graph neural network
CN111030795B (en) Transmission method of time reversal wireless energy-carrying communication system under non-ideal channel
Ma et al. Time-reversal tunneling effects for cloud radio access network
CN101207595B (en) A synchronization sequence transmission method and sending end device
CN101521531B (en) A method of enhancing TD-LTE downlink performance by using antenna selection
Salehi et al. Channel estimation based on learning automata for OFDM systems
CN111756418B (en) Wireless data energy simultaneous transmission method based on receiving space modulation
Ganhão et al. Analytical evaluation of iterative packet combining and multipacket detection schemes for SC-FDE
Hu et al. Orthogonal-time-frequency-space signal design for integrated data and energy transfer: Benefits from Doppler offsets
US8848773B2 (en) Rate control for a virtual diversity receiver
CN107425887B (en) A beamforming method in a multi-antenna untrusted relay network
CN103560987A (en) Cross-layer design method in multi-user MIMO based on multiple time delay estimations
CN111404588A (en) Physical layer secure transmission method of full-duplex cognitive eavesdropping network

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