CN112954792B - Multi-reflection device joint positioning and communication method based on environment backscattering - Google Patents

Multi-reflection device joint positioning and communication method based on environment backscattering Download PDF

Info

Publication number
CN112954792B
CN112954792B CN202110159173.8A CN202110159173A CN112954792B CN 112954792 B CN112954792 B CN 112954792B CN 202110159173 A CN202110159173 A CN 202110159173A CN 112954792 B CN112954792 B CN 112954792B
Authority
CN
China
Prior art keywords
signal
reflection device
reflection
receiver
radio frequency
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
CN202110159173.8A
Other languages
Chinese (zh)
Other versions
CN112954792A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110159173.8A priority Critical patent/CN112954792B/en
Publication of CN112954792A publication Critical patent/CN112954792A/en
Application granted granted Critical
Publication of CN112954792B publication Critical patent/CN112954792B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明属于无线通信技术领域,具体的说是涉及一种基于环境反向散射的多反射设备联合定位与通信方法。本发明使用接收机接收信号,接收到的信号通常包含来自环境射频源的直射链路信号和L个反射设备反向散射信号。通过空域信号处理技术,实现同时对多个反射设备进行定位。定位的同时,提供了一种在OFDM射频源下的低复杂度的解调反射信号的方案,无需解调射频源信息。即在对反射设备定位的同时,实现反射设备与接收机的通信。反射设备通过发送导频,接收端事先已知导频与反射设备的对应关系,通过接收端已知的导频和接收到的信号的导频进行相关运算,可以判决该反射信号和反射设备的对应关系,最终实现对某一确定的反射设备的信息解调和定位。

Figure 202110159173

The invention belongs to the technical field of wireless communication, and in particular relates to a joint positioning and communication method for multi-reflection equipment based on environmental backscattering. The present invention uses a receiver to receive a signal, and the received signal usually includes a direct link signal from an ambient radio frequency source and a backscattered signal from L reflective devices. Through the spatial signal processing technology, the positioning of multiple reflection devices can be realized at the same time. At the same time of positioning, a low-complexity solution for demodulating reflected signals under an OFDM radio frequency source is provided without demodulating radio frequency source information. That is, while positioning the reflection device, the communication between the reflection device and the receiver is realized. The reflection device transmits the pilot frequency, and the receiving end knows the corresponding relationship between the pilot frequency and the reflection device in advance. By performing a correlation operation between the pilot frequency known at the receiving end and the pilot frequency of the received signal, the reflected signal and the reflection device can be determined. Corresponding relationship, and finally realize the information demodulation and positioning of a certain reflective device.

Figure 202110159173

Description

一种基于环境反向散射的多反射设备联合定位与通信方法A method for joint positioning and communication of multi-reflection devices based on environmental backscattering

技术领域technical field

本发明属于无线通信技术领域,具体的说是涉及一种基于环境反向散射的多反射设备联合定位与通信方法。The invention belongs to the technical field of wireless communication, and in particular relates to a joint positioning and communication method for multi-reflection equipment based on environmental backscattering.

背景技术Background technique

环境反向散射技术是近几年来提出来的物联网的解决方案之一,它使用了环境中已经存在的射频源作为信号源,不需要像传统反向散射技术单独配置射频信号源。反向散射技术本身具备功耗低、成本低的优点,基于环境反向散射技术的应用在此基础上,由于免去了部署专用射频信号源,进一步地降低了成本,提升了频谱效率。Environmental backscattering technology is one of the IoT solutions proposed in recent years. It uses the RF source already existing in the environment as the signal source, and does not need to configure the RF signal source separately like the traditional backscattering technology. Backscattering technology itself has the advantages of low power consumption and low cost. Based on the application of environmental backscattering technology, the deployment of dedicated RF signal sources is eliminated, which further reduces costs and improves spectral efficiency.

随着科技发展,对于物联网设备同时进行通信和定位的需求越来越广泛。设计一种可以同时实现对环境反向散射通信系统中多个反射设备的通信和定位的系统,具有广阔的应用前景。比如,在物流管理中,使用机器人代替人工操作的需求越来越大,当机器人知道了物流件的信息和位置以后,才能进行相应的行动,因此要求对于使用了物联网设备的物流件既需要信息被正确解码,同时也需要知道其位置,才能更加方便使用机器人等设备代替繁琐的人工操作。With the development of science and technology, the demand for simultaneous communication and positioning of IoT devices is becoming more and more extensive. Designing a system that can simultaneously realize the communication and positioning of multiple reflection devices in the environmental backscatter communication system has broad application prospects. For example, in logistics management, there is an increasing demand for using robots instead of manual operations. Only after robots know the information and location of the logistics parts can they take corresponding actions. Therefore, it is required that logistics parts using IoT devices need both The information is decoded correctly, and its location also needs to be known, so that it is more convenient to use equipment such as robots to replace tedious manual operations.

但是由于环境反向散射技术使用了环境中已经存在的射频信号源,相较于传统的反向散射技术,在射频源信号干扰消除、对反向散射信号进行解调等方面带来了新的挑战。OFDM技术广泛应用于当今的商用的通信系统,因此,OFDM信号是环境反向散射通信系统常用的环境信号,它具有循环前缀结构,因此,针对这种常见的射频源信号的特点设计基于环境反向散射的多反射设备联合定位与通信系统也具有重要的现实意义。However, since the environmental backscattering technology uses the radio frequency signal source that already exists in the environment, compared with the traditional backscattering technology, it brings new advantages in eliminating the interference of the radio frequency source signal and demodulating the backscattered signal. challenge. OFDM technology is widely used in today's commercial communication systems. Therefore, OFDM signals are commonly used environmental signals in environmental backscatter communication systems. It has a cyclic prefix structure. Therefore, according to the characteristics of this common RF source signal, the design based on environmental feedback The combined positioning and communication system of the multi-reflection equipment for the scattering also has important practical significance.

发明内容SUMMARY OF THE INVENTION

本发明的主要内容是提出一种基于环境反向散射的多反射设备联合定位与通信方法,实现了对多个反射设备同时进行定位和通。The main content of the present invention is to propose a joint positioning and communication method for multiple reflection devices based on environmental backscattering, which realizes the simultaneous positioning and communication of multiple reflection devices.

本发明使用接收机接收信号,接收到的信号通常包含来自环境射频源的直射链路信号和L个反射设备反向散射信号。通过空域信号处理技术,实现同时对多个反射设备进行定位。本发明在实现对多个反射信号的定位的同时,提供了一种在OFDM射频源下的低复杂度的解调反射信号的方案,无需解调射频源信息。即在对反射设备定位的同时,实现反射设备与接收机的通信。反射设备通过发送导频,接收端事先已知导频与反射设备的对应关系,通过接收端已知的导频和接收到的信号的导频进行相关运算,可以判决该反射信号和反射设备的对应关系,最终实现对某一确定的反射设备的信息解调和定位。The present invention uses a receiver to receive a signal, and the received signal usually includes a direct link signal from an ambient radio frequency source and a backscattered signal from L reflective devices. Through the spatial signal processing technology, the positioning of multiple reflection devices can be realized at the same time. The present invention provides a low-complexity solution for demodulating reflected signals under an OFDM radio frequency source while realizing the positioning of multiple reflected signals, without demodulating radio frequency source information. That is, while positioning the reflection device, the communication between the reflection device and the receiver is realized. The reflection device transmits the pilot frequency, and the receiving end knows the corresponding relationship between the pilot frequency and the reflection device in advance. By performing a correlation operation between the pilot frequency known at the receiving end and the pilot frequency of the received signal, the reflected signal and the reflection device can be determined. Corresponding relationship, and finally realize the information demodulation and positioning of a certain reflective device.

本发明的技术方案是:The technical scheme of the present invention is:

一种基于环境反向散射的多反射设备联合定位与通信方法,包括环境射频源、L个反射设备、一个有M根天线的接收机,反射设备的编号为l=0,1,...L-1,接收机坐标为(xr,yr),M≥L+1;其特征在于,所述联合定位与通信方法包括以下步骤:A method for joint positioning and communication of multi-reflection devices based on environmental backscatter, comprising an ambient radio frequency source, L reflective devices, and a receiver with M antennas, and the number of the reflective devices is l=0,1,... L-1, the receiver coordinates are (x r , y r ), M≥L+1; it is characterized in that the joint positioning and communication method includes the following steps:

S1、环境射频源发射射频源信号,反射设备对环境射频源信号进行反向散射,其中,反射设备的波形设计方法为:S1. The ambient radio frequency source transmits the radio frequency source signal, and the reflection device backscatters the ambient radio frequency source signal. The waveform design method of the reflection device is as follows:

假设反射设备符号周期是OFDM射频源符号周期的K倍,每个反射设备发出的二元信息用Bl(m)表示,Bl(m)=±1,当传输码片中的某一个符号B(m)=-1时,在对应的第k个OFDM符号周期内,设计如下的反射设备波形:Assuming that the symbol period of the reflection device is K times the symbol period of the OFDM radio frequency source, the binary information sent by each reflection device is represented by B l (m), B l (m)=±1, when a certain symbol in the transmission chip is transmitted When B(m)=-1, in the corresponding kth OFDM symbol period, design the following reflection device waveform:

Figure GDA0003476261240000021
Figure GDA0003476261240000021

其中,Nc是OFDM射频源符号的循环前缀的长度,N是数据部分长度,N+Nc构成了一个完整的OFDM符号的长度,k=1,...K代表在同一个反射设备符号周期内的OFDM射频源符号的序号;Among them, N c is the length of the cyclic prefix of the OFDM radio frequency source symbol, N is the length of the data part, N+N c constitutes the length of a complete OFDM symbol, k=1,...K represents the symbol in the same reflection device The sequence number of the OFDM radio source symbol in the period;

当传输码片的某一个符号B(m)=1时,反射设备波形恒定为1;When a certain symbol B(m)=1 of the transmission chip, the waveform of the reflection device is always 1;

S2、接收机接收信号,进行DOA估计,估计得到直射链路信号和多个反向散射信号的到达角θd和θL=[θ012,...,θl,...θL-1];S2. The receiver receives the signal, performs DOA estimation, and obtains the estimated angles of arrival of the direct link signal and multiple backscattered signals θ d and θ L =[θ 012 ,...,θ l , ...θ L-1 ];

S3、对所有估计出来的到达角的信号各执行一次空域滤波,每一次空域滤波时,设计权向量进行波束成形,处理接收的信号,只保留一个角度的信号,抑制其他方向信号的功率,获得所有到达角对应的方向信号;S3. Perform spatial filtering once for all the estimated angle of arrival signals. In each spatial filtering, design a weight vector to perform beamforming, process the received signals, keep only the signals of one angle, suppress the power of signals in other directions, and obtain Direction signals corresponding to all angles of arrival;

S4、计算空域滤波后的信号z(n)的功率,设信号有N个采样点,则信号的功率Pz计算方法为:S4. Calculate the power of the signal z (n) after spatial filtering. Suppose the signal has N sampling points, and the calculation method of the power Pz of the signal is:

Figure GDA0003476261240000031
Figure GDA0003476261240000031

S5、获得所有到达角的方向信号对应的功率以后,根据反向散射信号相比较于直射链路信号会经历两次衰落的特性,区分出直射信号功率

Figure GDA0003476261240000032
和反射信号功率
Figure GDA0003476261240000033
通过直射信号功率
Figure GDA0003476261240000034
计算射频源到接收机的距离df:S5. After obtaining the power corresponding to the direction signals of all the angles of arrival, according to the characteristic that the backscattered signal will experience two fading compared with the direct link signal, the power of the direct signal is distinguished.
Figure GDA0003476261240000032
and reflected signal power
Figure GDA0003476261240000033
through direct signal power
Figure GDA0003476261240000034
Calculate the distance d f from the RF source to the receiver:

Figure GDA0003476261240000035
Figure GDA0003476261240000035

其中,n代表该环境的衰减系数,β是一个常数,由载波频率和环境决定,n和β均可在实际测量中获得,ps为发射功率;Among them, n represents the attenuation coefficient of the environment, β is a constant, which is determined by the carrier frequency and the environment, both n and β can be obtained in the actual measurement, and ps is the transmit power;

其余的信号为反射链路信号,通过功率估计射频源到反射设备距离和反射设备到接收机距离的乘积D:The rest of the signals are reflected link signals, and the power is used to estimate the product D of the distance from the RF source to the reflecting device and the distance from the reflecting device to the receiver:

Figure GDA0003476261240000036
Figure GDA0003476261240000036

其中,

Figure GDA0003476261240000037
为环境射频源到第l个反射设备的距离,
Figure GDA0003476261240000038
为第l个反射设备到接收机的距离,αl为第l个反射设备的反射系数;in,
Figure GDA0003476261240000037
is the distance from the ambient RF source to the lth reflecting device,
Figure GDA0003476261240000038
is the distance from the lth reflection device to the receiver, α l is the reflection coefficient of the lth reflection device;

S6、通过几何原理估计反射设备到接收机的距离,具体为根据三角形的余弦定理建立如下方程:S6. Estimate the distance from the reflection device to the receiver through geometric principles, specifically establishing the following equation according to the cosine theorem of the triangle:

Figure GDA0003476261240000039
Figure GDA0003476261240000039

其中,Δθ为直射径和反射径之间的夹角,Δθ=|θdl|,通过方程解得

Figure GDA00034762612400000310
Figure GDA00034762612400000311
Among them, Δθ is the angle between the direct radiation diameter and the reflected diameter, Δθ=|θ dl |, which can be obtained by solving the equation
Figure GDA00034762612400000310
and
Figure GDA00034762612400000311

S7、利用反向散射信号到达角信息θl和反射设备到接收机的距离信息

Figure GDA00034762612400000312
计算反射设备的坐标(xl,yl):S7. Use the angle of arrival information θ l of the backscattered signal and the distance information from the reflection device to the receiver
Figure GDA00034762612400000312
Calculate the coordinates of the reflective device (x l , y l ):

Figure GDA0003476261240000041
Figure GDA0003476261240000041

实现定位;achieve positioning;

S8、解调出反射信号的信息:S8, demodulate the information of the reflected signal:

假设K=1,即反射设备波形c(n)的一个周期的持续时长和一个OFDM符号持续时长相同,同时忽略噪声,在一个OFDM符号周期内,建立:Assuming K=1, that is, the duration of one cycle of the reflection device waveform c(n) is the same as the duration of one OFDM symbol, while ignoring noise, within one OFDM symbol cycle, establish:

Figure GDA0003476261240000042
Figure GDA0003476261240000042

其中,Lh是信道的时延扩展,代表被码间串扰污染的样点数,将这些样点丢弃;则有:Among them, L h is the delay spread of the channel, which represents the number of samples polluted by inter-symbol crosstalk, and these samples are discarded; then there are:

u(m)=|Kl|2|s(n)|2B(m)u(m)=|K l | 2 |s(n)| 2 B(m)

其中,

Figure GDA0003476261240000043
是一个常数,pl=1,s(n)是射频源信号,解调时对B(m)判决法则如下:in,
Figure GDA0003476261240000043
is a constant, p l =1, s(n) is the radio frequency source signal, the decision rule for B(m) during demodulation is as follows:

Figure GDA0003476261240000044
Figure GDA0003476261240000044

从而恢复出反射设备传输的码元;Thereby, the symbols transmitted by the reflection device are recovered;

S9、使用已知的导频序列与解码后的反射信号的导频部分进行相关运算,当解调后的信号与某个已知导频序列进行相关运算得到的结果最大时,即确定信号与已知导频的对应关系,因此在传输信息的同时,判断出该信号属于哪一个反射设备。S9. Use the known pilot sequence to perform a correlation operation with the pilot part of the decoded reflected signal. When the result of the correlation operation between the demodulated signal and a certain known pilot sequence is the largest, it is determined that the signal is related to The corresponding relationship of the pilot frequency is known, so it is determined which reflection device the signal belongs to while transmitting the information.

本发明的有益效果是:本发明提出一种基于环境反向散射的多反射设备联合定位与通信方法,在普遍使用的OFDM环境射频源场景下,同时对多个反射设备定位的基础上解码出反射设备的信息,具有很强的应用价值。The beneficial effects of the present invention are as follows: the present invention proposes a joint positioning and communication method for multi-reflection devices based on environmental backscattering. It reflects the information of the equipment and has strong application value.

附图说明Description of drawings

图1示出了本发明的系统构成Fig. 1 shows the system composition of the present invention

图2示出了本发明所考虑的一种反射设备的波形设计Fig. 2 shows the waveform design of a reflection device considered by the present invention

图3示出了本发明所考虑的一种信号结构说明Figure 3 shows an illustration of a signal structure considered by the present invention

图4示出了本发明所考虑的技术方案的信号被正确判决到正确反射设备仿真结果Fig. 4 shows the simulation result of the signal of the technical solution considered by the present invention being correctly judged to the correct reflection device

图5示出了本发明所考虑的技术方案的反射设备信息解调性能仿真结果Fig. 5 shows the simulation result of the information demodulation performance of the reflection device of the technical solution considered by the present invention

具体实施方式Detailed ways

下面结合附图对本发明进行详细的描述。The present invention will be described in detail below with reference to the accompanying drawings.

如图1所示,本发明基于的环境反向散射系统,包括环境射频源、L个反射设备、一个有M根天线的接收机,M≥L+1;环境射频源发射射频源信号

Figure GDA0003476261240000051
其中s(t)是功率归一化射频源基带信号,发射功率为ps。信道是强LOS径,图中
Figure GDA0003476261240000052
代表环境射频源到反向散射体的信道参数,
Figure GDA0003476261240000053
代表环境射频源到接收机的信道参数,
Figure GDA0003476261240000054
代表反向散射体到接收机的信道参数。n代表该环境的衰减系数,可以通过在某一具体环境实测而得,一般在2到4之间。df为环境射频源到接收机的距离,dh为环境射频源到反向散射体的距离,dg为反向散射体到接收机的距离。As shown in FIG. 1, the environmental backscattering system based on the present invention includes an environmental radio frequency source, L reflection devices, and a receiver with M antennas, M≥L+1; the environmental radio frequency source transmits the radio frequency source signal
Figure GDA0003476261240000051
where s(t) is the power-normalized RF source baseband signal, and the transmit power is p s . The channel is a strong LOS path, as shown in the figure
Figure GDA0003476261240000052
represents the channel parameters from the ambient RF source to the backscatterer,
Figure GDA0003476261240000053
represents the channel parameters from the ambient RF source to the receiver,
Figure GDA0003476261240000054
Represents the channel parameters from the backscatterer to the receiver. n represents the attenuation coefficient of the environment, which can be obtained by actual measurement in a specific environment, generally between 2 and 4. d f is the distance from the ambient RF source to the receiver, dh is the distance from the ambient RF source to the backscatterer, and dg is the distance from the backscatterer to the receiver.

以均匀线性阵列为例说明接收机接收的模型,反射设备的个数是L,对反射设备进行编号,为l=0,1,...L-1,设直射信号到达角为θd,编号为l的反射设备的反射信号到达角为θlTaking a uniform linear array as an example to illustrate the model of receiver reception, the number of reflection devices is L, and the number of reflection devices is l=0,1,...L-1, and the arrival angle of the direct signal is θ d , The angle of arrival of the reflected signal from the reflecting device numbered l is θ l .

接收到的信号可以表示为The received signal can be expressed as

Figure GDA0003476261240000055
Figure GDA0003476261240000055

其中l=0,1,...L-1。

Figure GDA0003476261240000056
是带通噪声,
Figure GDA0003476261240000057
是射频源到编号为l的反射设备的距离,
Figure GDA0003476261240000058
是编号为l的反射设备到接收机的距离。cl(t)为编号为l的反射设备的要传输的基带信号,αl为反射系数。
Figure GDA0003476261240000059
Figure GDA00034762612400000510
为方向矢量,其中
Figure GDA00034762612400000511
为空间相位,fc是载波频率,d是天线间距,c是光速,λ是载波波长。where l=0,1,...L-1.
Figure GDA0003476261240000056
is the bandpass noise,
Figure GDA0003476261240000057
is the distance from the RF source to the reflecting device numbered l,
Figure GDA0003476261240000058
is the distance from the reflective device numbered l to the receiver. c l (t) is the baseband signal to be transmitted by the reflection device numbered l, and α l is the reflection coefficient.
Figure GDA0003476261240000059
and
Figure GDA00034762612400000510
is the direction vector, where
Figure GDA00034762612400000511
is the spatial phase, f c is the carrier frequency, d is the antenna spacing, c is the speed of light, and λ is the carrier wavelength.

接收信号的数字基带形式可以表示为The digital baseband form of the received signal can be expressed as

Figure GDA0003476261240000061
Figure GDA0003476261240000061

其中ω(n)是复基带噪声,服从循环对称复高斯分布

Figure GDA0003476261240000062
σ2是噪声功率。where ω(n) is the complex baseband noise, which follows a cyclic symmetric complex Gaussian distribution
Figure GDA0003476261240000062
σ 2 is the noise power.

本发明设计了反射信号的波形,具体波形设计的方案如下:The present invention designs the waveform of the reflected signal, and the specific waveform design scheme is as follows:

假设反射设备符号周期是OFDM射频源符号周期的K倍。每个反射设备发出的二元信息用Bl(m)表示,Bl(m)=±1。当传输码片中的某一个符号B(m)=-1时,设计如下的反射设备波形It is assumed that the reflecting device symbol period is K times the OFDM radio source symbol period. The binary information emitted by each reflective device is denoted by B1(m), where B1 ( m)=± 1 . When a certain symbol B(m)=-1 in the transmission chip, the following reflection device waveform is designed

Figure GDA0003476261240000063
Figure GDA0003476261240000063

Nc是OFDM射频源符号的循环前缀的长度,N是数据部分长度,N+Nc构成了一个完整的OFDM符号的长度。N c is the length of the cyclic prefix of the OFDM radio frequency source symbol, N is the length of the data part, and N+N c constitutes the length of a complete OFDM symbol.

当传输码片的某一个符号B(m)=1时,反射设备波形恒定为1.When a certain symbol of the transmission chip B(m)=1, the reflection device waveform is constant at 1.

如图2所示,当传输-1的时候,波形设计为在一个OFDM射频源符号中间进行一次跳变,从1变为-1,而传输1的时候,波形设计为没有跳变。As shown in Figure 2, when -1 is transmitted, the waveform is designed to perform a transition in the middle of an OFDM radio frequency source symbol, from 1 to -1, while when 1 is transmitted, the waveform is designed to have no transition.

使用DOA估计算法可以估计出L+1个到达角,这里以Root-MUSIC算法(Barabell,"Improving the resolution performance of eigenstructure-based direction-finding algorithms,"ICASSP'83.IEEE International Conference on Acoustics,Speech,and Signal Processing,Boston,Massachusetts,USA,1983,pp.336-339,doi:10.1109/ICASSP.1983.1172124)为例说明估计算法:Using the DOA estimation algorithm, L+1 angles of arrival can be estimated. Here, the Root-MUSIC algorithm (Barabell, "Improving the resolution performance of eigenstructure-based direction-finding algorithms," ICASSP'83. IEEE International Conference on Acoustics, Speech, and Signal Processing, Boston, Massachusetts, USA, 1983, pp.336-339, doi: 10.1109/ICASSP.1983.1172124) as an example to illustrate the estimation algorithm:

1.根据设定的采样数目N,进行N次采样获得样本序列y(n),其中

Figure GDA0003476261240000065
1. According to the set sampling number N, perform N sampling to obtain the sample sequence y(n), where
Figure GDA0003476261240000065

2.计算y(n)自相关矩阵R,R=E{y(n)yH(n)},这里可以用时间平均代替统计平均,即

Figure GDA0003476261240000064
2. Calculate the y(n) autocorrelation matrix R, R=E{y(n)y H (n)}, where time average can be used instead of statistical average, that is
Figure GDA0003476261240000064

3.对R进行特征值分解,由于目标信号数目有L+1个,因此得到最小的M-L-1个特征值对应的归一化特征向量ui,其中i=L+2,L+3…,M。3. Perform eigenvalue decomposition on R. Since the number of target signals is L+1, the normalized eigenvector ui corresponding to the smallest ML-1 eigenvalue is obtained, where i =L+2, L+3... , M.

ui=[ui0,ui1,...ui(m-1)]T u i =[u i0 ,u i1 ,...u i(m-1) ] T

4.构造向量a(z),a(z)=[1,z-1,...,z-(M-1)]T4. Construct the vector a(z), a(z)=[1, z -1 ,...,z- (M-1) ] T .

构造如下函数:

Figure GDA0003476261240000071
Construct the following function:
Figure GDA0003476261240000071

Figure GDA0003476261240000072
乘以自身共轭转置,得到
Figure GDA0003476261240000072
Multiplying by the self-conjugate transpose, we get

Figure GDA0003476261240000073
i=3,4,...,M。
Figure GDA0003476261240000073
i=3,4,...,M.

5.定义多项式5. Define Polynomials

Figure GDA0003476261240000074
Figure GDA0003476261240000074

令其为0,解方程。求得z的值。Set it to 0 and solve the equation. Find the value of z.

6.目标信号数目为L+1个,对z值求模减去1后再取绝对值从小到大进行排序。这样会找到2L+2个最接近模为1的根(会出现重根)。由于反向散射体反射的信号经历了两次衰减,其信号强度远远小于直射信号。最接近单位圆的根属于直射链路信号。其余的根去重后可以得到反射链路信号到达角。6. The number of target signals is L+1, and the z value is modulo subtracted by 1, and then the absolute value is sorted from small to large. This will find the 2L+2 closest roots modulo 1 (multiple roots will occur). Since the signal reflected by the backscatterer undergoes two attenuations, its signal strength is much smaller than that of the direct signal. The roots closest to the unit circle belong to the direct link signal. The angle of arrival of the reflected link signal can be obtained after the rest of the roots are deduplicated.

7.由于z=e,φ是空间角频率,φ=-πsinθ。因此7. Since z=e , φ is the spatial angular frequency, φ=−πsinθ. therefore

Figure GDA0003476261240000075
Figure GDA0003476261240000075

计算得到信号到达方向估计。The calculation obtains an estimate of the direction of arrival of the signal.

得到直射信号到达角度估计值

Figure GDA0003476261240000076
和反射信号到达角度估计值
Figure GDA0003476261240000077
后,进行空域滤波。目的是只保留某一个方向的信号,将其他方向信号置零。Obtain the estimated value of the angle of arrival of the direct signal
Figure GDA0003476261240000076
and the estimated angle of arrival of the reflected signal
Figure GDA0003476261240000077
Afterwards, spatial filtering is performed. The purpose is to keep only the signals in a certain direction, and set the signals in other directions to zero.

空域滤波的算法有很多种,在此举例一种空域滤波的算法,以保留

Figure GDA0003476261240000078
方向的反射信号,将其他估计出来的方向的信号置零为例说明:There are many kinds of spatial filtering algorithms. Here is an example of a spatial filtering algorithm to preserve the
Figure GDA0003476261240000078
The reflected signal of the direction, and the signal of other estimated directions is set to zero as an example to illustrate:

写出矩阵

Figure GDA0003476261240000081
代表要迫零的方向的集合。提出优化问题:write out the matrix
Figure GDA0003476261240000081
Represents a collection of directions to zero forcing. Ask the optimization problem:

Figure GDA0003476261240000082
Figure GDA0003476261240000082

s.t wHA=0st w H A = 0

具体算法步骤如下:The specific algorithm steps are as follows:

1.对想要迫零的方向,通过方向向量构造矩阵

Figure GDA0003476261240000083
1. For the direction you want to force zero, construct a matrix through the direction vector
Figure GDA0003476261240000083

2.求解BHA=0可通过AHB=0求解2. Solving B H A=0 can be solved by A H B=0

3.

Figure GDA0003476261240000084
其中B+是B矩阵的M-P广义逆矩阵。3.
Figure GDA0003476261240000084
where B + is the MP generalized inverse of the B matrix.

这里以只保留到达角为

Figure GDA0003476261240000085
的反射信号为例,空域滤波完成后,信号形式为:Here, only the arrival angle is retained as
Figure GDA0003476261240000085
The reflected signal of , for example, after the spatial filtering is completed, the signal form is:

z(n)=wHy(n)z(n)=w H y(n)

即该信号为:That is, the signal is:

Figure GDA0003476261240000086
Figure GDA0003476261240000086

Figure GDA0003476261240000087
Figure GDA0003476261240000087

其中,w是空域滤波过程中设计的权向量,p0=wHa(θ0),εl=wHa(θld=wHa(θd)代表对各个信号分量功率的抑制程度,取决于DOA估计的精确度。如果DOA估计很精确,则p0≈1,εl和εd都是很小的值,约为0,可以归为噪声,以下基于这样的情况进行讨论,即p0=1,ε=0。Among them, w is the weight vector designed in the spatial filtering process, p 0 =w H a(θ 0 ), ε l =w H a(θ ld =w H a(θ d ) represents the power of each signal component The degree of inhibition depends on the accuracy of the DOA estimation. If the DOA estimation is very accurate, then p 0 ≈1, ε l and ε d are both small values, about 0, which can be classified as noise. The following discussion is based on such a situation, that is, p 0 =1, ε = 0 .

测量z0(n)信号的功率为

Figure GDA0003476261240000088
当p0≈1,反射系数α0已知,n和β实现进行实地测量后已知的情况下,则可以估计出
Figure GDA0003476261240000089
的值,即射频源到反射设备的距离和反射设备到接收机的距离的乘积D。The power of the measured z 0 (n) signal is
Figure GDA0003476261240000088
When p 0 ≈ 1, the reflection coefficient α 0 is known, and n and β are known after field measurements, then it can be estimated that
Figure GDA0003476261240000089
The value of , that is, the product D of the distance from the RF source to the reflective device and the distance from the reflective device to the receiver.

Figure GDA00034762612400000810
Figure GDA00034762612400000810

保留直射链路的信号,抑制所有反射信号以估计射频源到接收机的距离,按照上述类似的方法进行空域滤波,可以得到Retain the signal of the direct link, suppress all reflected signals to estimate the distance from the RF source to the receiver, and perform spatial filtering in a similar way to the above, we can get

Figure GDA0003476261240000091
Figure GDA0003476261240000091

式中pd=wHa(θd),pd≈1,εl是很小的值,约为0,可以归为噪声。因此,可以求出信号zd(n)的功率

Figure GDA0003476261240000092
计算出df的值,以估计出射频源到接收机的距离。In the formula, p d =w H a(θ d ), p d ≈1, and ε l is a very small value, about 0, which can be classified as noise. Therefore, the power of the signal z d (n) can be found
Figure GDA0003476261240000092
Calculate the value of d f to estimate the distance from the RF source to the receiver.

Figure GDA0003476261240000093
Figure GDA0003476261240000093

因此对反射信号的定位问题可以建模成如下的几何问题:Therefore, the positioning problem of the reflected signal can be modeled as the following geometric problem:

θd和θ0在执行DOA估计以后可以得到,而df通过只保留直射路径的信号可以估计出来,

Figure GDA0003476261240000094
作为两条径距离的乘积也可以估计出来。那么由三角形的余弦定理。直射径和反射径之间的夹角为Δθ=|θd0|,可以建立如下方程:θ d and θ 0 can be obtained after performing DOA estimation, and d f can be estimated by keeping only the direct path signal,
Figure GDA0003476261240000094
It can also be estimated as the product of the two radial distances. Then by the cosine law of triangles. The angle between the direct diameter and the reflected diameter is Δθ=|θ d0 |, and the following equation can be established:

Figure GDA0003476261240000095
Figure GDA0003476261240000095

该方程中,D和Δθ和df均已知,两个方程可以解得

Figure GDA0003476261240000096
Figure GDA0003476261240000097
这样,反射设备到达接收机的距离
Figure GDA0003476261240000098
即可估计出来,同时已知信号到达角θ0,则可以估计出反射设备的位置。In this equation, D and Δθ and d f are known, and the two equations can be solved as
Figure GDA0003476261240000096
and
Figure GDA0003476261240000097
In this way, the distance from the reflecting device to the receiver
Figure GDA0003476261240000098
can be estimated, and at the same time the signal arrival angle θ 0 is known, the position of the reflecting device can be estimated.

反射设备坐标计算方法如下,设已知接收机的坐标为(xr,yr),反射设备的坐标为(x,y),到达角为θThe calculation method of the coordinates of the reflection device is as follows. Assume that the coordinates of the known receiver are (x r , y r ), the coordinates of the reflection device are (x, y), and the angle of arrival is θ

Figure GDA0003476261240000099
Figure GDA0003476261240000099

最终得到该反射信号的定位结果。Finally, the positioning result of the reflected signal is obtained.

定位同时解调反射信号的方法如下:The method to locate and demodulate the reflected signal at the same time is as follows:

为了使得分析方便,假设K=1,即反射设备波形c(n)的一个周期的持续时长和一个OFDM符号持续时长相同,该假设同样不会影响理论推导的正确性。即一个完整OFDM符号对应着一个反射设备的符号Bl(m),即一个Bl(m)符号周期内有N+Nc个c(n)样点。For the convenience of analysis, it is assumed that K=1, that is, the duration of one cycle of the reflection device waveform c(n) is the same as the duration of one OFDM symbol. This assumption also does not affect the correctness of the theoretical derivation. That is, a complete OFDM symbol corresponds to a symbol B l (m) of a reflection device, that is, there are N+N c c(n) samples in a B l (m) symbol period.

空域滤波后的信号如下:The spatially filtered signal is as follows:

Figure GDA0003476261240000101
Figure GDA0003476261240000101

Figure GDA0003476261240000102
Figure GDA0003476261240000102

可以将模型简化为如下形式:The model can be simplified to the following form:

z0(n)=K0s(n)c0(n)+ω(n)z 0 (n)=K 0 s(n)c 0 (n)+ω(n)

其中

Figure GDA0003476261240000103
是一个常数,
Figure GDA0003476261240000104
Figure GDA0003476261240000105
由于空域滤波已经置零,将其归入ω(n)项内。in
Figure GDA0003476261240000103
is a constant,
Figure GDA0003476261240000104
and
Figure GDA0003476261240000105
Since the spatial filtering has been zeroed, it is included in the ω(n) term.

如图3所示,假设信道h的时延扩展为Lh,那么将会对CP的符号造成一定的码间串扰,CP的长度大于时延扩展,因此还有一部分CP的符号不会造成码间串扰,就利用没有产生码间串扰的部分的样点进行处理。即ISI不会影响推导的正确性。As shown in Fig. 3 , assuming that the delay spread of channel h is L h , a certain amount of intersymbol interference will be caused to the symbols of the CP, and the length of the CP is greater than the delay spread, so there are still some symbols of the CP that will not cause code If the inter-symbol crosstalk is detected, the sample points of the part where the inter-symbol crosstalk is not generated are used for processing. That is, ISI does not affect the correctness of the derivation.

由于循环前缀的性质:Due to the nature of cyclic prefixes:

s(n)=s(n+N),n=0,...,Nc-1s(n)=s(n+N), n=0,...,N c -1

由于设计的波形的性质:Due to the nature of the designed waveform:

Figure GDA0003476261240000106
Figure GDA0003476261240000106

在一个OFDM符号周期,即一个B0(m)符号周期内,忽略噪声的影响,如图3所示,z0(n)信号有如下性质:当n=Lh-1,...,Nc-1时In one OFDM symbol period, that is, one B 0 (m) symbol period, ignoring the influence of noise, as shown in Figure 3, the z 0 (n) signal has the following properties: when n=L h -1,..., When N c -1

Figure GDA0003476261240000107
Figure GDA0003476261240000107

对该信号进行处理,令:To process the signal, let:

Figure GDA0003476261240000111
Figure GDA0003476261240000111

Figure GDA0003476261240000112
ωcross(n)为共轭相乘运算中产生的交叉项,也归于噪声项。but
Figure GDA0003476261240000112
ω cross (n) is the cross term produced in the conjugate multiplication operation, which is also attributed to the noise term.

因此,如果忽略噪声,则u(m)=|Ko|2|s(n)|2B(m),即处理后的u(m)是原发射符号B(m)乘上一个常数。Therefore, if noise is ignored, u(m)=|K o | 2 |s(n)| 2 B(m), that is, the processed u(m) is the original transmitted symbol B(m) multiplied by a constant.

解调时对B(m)判决法则如下:The decision rule for B(m) during demodulation is as follows:

Figure GDA0003476261240000113
Figure GDA0003476261240000113

如此即可恢复出反射设备传输的码元。In this way, the symbols transmitted by the reflection device can be recovered.

使用已知的导频序列与解码后的反射信号的导频部分进行相关运算,当解调后的信号与某个已知导频序列进行相关运算得到的结果最大时,即可确定信号与已知导频的对应关系,因此可以在传输信息的同时,知道该信号属于哪一个反射设备。Use the known pilot sequence to perform the correlation operation with the pilot part of the decoded reflected signal. When the result of the correlation operation between the demodulated signal and a certain known pilot sequence is the largest, it can be determined that the signal is related to the already decoded signal. The corresponding relationship of the pilot frequency is known, so it is possible to know which reflecting device the signal belongs to while transmitting the information.

仿真分析:Simulation analysis:

反射设备设为4个,定位区域为一个10米×10米的空间,接收机在该空间内,环境射频源在这个空间以外。The reflective devices are set to 4, the positioning area is a space of 10 meters × 10 meters, the receiver is in this space, and the ambient radio frequency source is outside this space.

环境射频源的OFDM符号的Nc设为16,N设为64,N+Nc为80,同时也是对一个OFDM符号的采样点数。反射设备符号周期和OFDM符号周期相同。一共采集80000个样点,即1000个符号序列做DOA估计。信源发射功率固定为1,信源坐标设置为(-5,15),接收机接收天线数目为8。信道参数中β设为1,距离衰减系数设为2.5。接收机固定坐标为(5,0).直射信号到达角为-33.7°,信源发射机到接收机的距离为18.03米,由之前信号模型,直射链路信号信噪比与发射信号信噪比之比为Δγd=-31.4dB。N c of the OFDM symbol of the ambient radio frequency source is set to 16, N is set to 64, and N+N c is 80, which is also the number of sampling points for one OFDM symbol. The reflector symbol period is the same as the OFDM symbol period. A total of 80,000 samples are collected, that is, 1,000 symbol sequences for DOA estimation. The transmit power of the source is fixed to 1, the coordinate of the source is set to (-5, 15), and the number of receiver antennas is 8. In the channel parameters, β is set to 1, and the distance attenuation coefficient is set to 2.5. The fixed coordinate of the receiver is (5,0). The arrival angle of the direct signal is -33.7°, and the distance from the source transmitter to the receiver is 18.03 meters. According to the previous signal model, the signal-to-noise ratio of the direct-link signal is related to the signal-to-noise ratio of the transmitted signal. The ratio of the ratios is Δγ d = -31.4 dB.

反射设备的反射系数设为α=0.2+0.3j。仿真中固定反射设备的坐标,编号为0的反射设备坐标为(2,3),编号为1的反射设备坐标为(4,7),编号为2的反射设备坐标为(6,1),编号为3的反射设备坐标为(7,9)。即,4个反射设备到达接收机的距离分别为4.24米,7.07米,1.41米和9.22米,4个反射设备到达接收机的真实方向分别为45.0°,8.1301°,-45.0°,-12.53°。反射设备的导频序列长度设置为64个反射设备符号。The reflection coefficient of the reflecting device is set to α=0.2+0.3j. In the simulation, the coordinates of the reflection device are fixed. The coordinates of the reflection device numbered 0 are (2,3), the coordinates of the reflection device numbered 1 are (4,7), and the coordinates of the reflection device numbered 2 are (6,1). The reflection device with number 3 has coordinates (7,9). That is, the distances from the four reflection devices to the receiver are 4.24 meters, 7.07 meters, 1.41 meters, and 9.22 meters, respectively, and the true directions of the four reflection devices to the receiver are 45.0°, 8.1301°, -45.0°, and -12.53°, respectively. . The length of the pilot sequence of the reflection device is set to 64 symbols of the reflection device.

由之前的信号模型,可以计算固定坐标的第i个反射设备的反射信号的功率和发射信号的信噪比之比为:From the previous signal model, the ratio of the power of the reflected signal of the ith reflecting device with fixed coordinates to the signal-to-noise ratio of the transmitted signal can be calculated as:

Figure GDA0003476261240000121
Figure GDA0003476261240000121

在该仿真条件下,可以知道4个反射设备的Δγi分别为-53.1dB,-57.1dB,-43.9dB,-61.2dB.可见不同位置的反射设备的反射信号功率相差可以达到约20dB.Under the simulation conditions, it can be known that the Δγ i of the four reflection devices are -53.1dB, -57.1dB, -43.9dB, and -61.2dB respectively. It can be seen that the reflected signal power difference of the reflection devices at different positions can reach about 20dB.

更改发射信号功率信噪比γs从30dB到60dB,直射链路接收信号信噪比为γs+Δγd,4个反射设备接收信号信噪比为γs+Δγi。每一轮进行若干次蒙特卡罗实验,每一次蒙特卡罗实验中,统计反射信号判决为相应反射设备的准确率,统计每个反射设备的误码率。Change the signal-to-noise ratio γ s of the transmit signal power from 30dB to 60dB, the signal-to-noise ratio of the received signal of the direct link is γ s +Δγ d , and the signal-to-noise ratio of the received signal of the four reflection devices is γ s +Δγ i . Several Monte Carlo experiments are carried out in each round. In each Monte Carlo experiment, the statistical reflection signal is determined as the accuracy of the corresponding reflection device, and the bit error rate of each reflection device is counted.

如图4所示,这是多反射设备的信号同时接入,对反射信号进行判决为反射设备的仿真结果。可以看出发射信号信噪比大于40dB时,判决正确的概率很高。说明该方案是有效的。同时,最低判决概率在20%左右,这是因为对于一个方向的信号,假设完全随机判决,则可能判决为4个反射设备反射信号或者直射信号,恰好判决正确的概率正是20%。As shown in Figure 4, this is the simultaneous access of signals of multiple reflection devices, and the reflection signal is judged as the simulation result of the reflection device. It can be seen that when the signal-to-noise ratio of the transmitted signal is greater than 40dB, the probability of a correct decision is very high. Show that the program is effective. At the same time, the lowest decision probability is about 20%. This is because for a signal in one direction, assuming a completely random decision, it may be decided as four reflected signals or direct signals, and the probability of the correct decision is exactly 20%.

图5是多反射设备信号解调的误码率仿真结果。可以看出,误码率性能受接收信号信噪比影响比较大。同时,当信噪比足够高时,反射设备的解调的误码率也是令人满意的,说明了该信号解调方案的正确性。Fig. 5 is the simulation result of the bit error rate of the signal demodulation of the multi-reflection device. It can be seen that the bit error rate performance is greatly affected by the signal-to-noise ratio of the received signal. At the same time, when the signal-to-noise ratio is high enough, the bit error rate of the demodulation of the reflection device is also satisfactory, which shows the correctness of the signal demodulation scheme.

Claims (1)

1.一种基于环境反向散射的多反射设备联合定位与通信方法,包括环境射频源、L个反射设备、一个有M根天线的接收机,反射设备的编号为l=0,1,...L-1,接收机坐标为(xr,yr),M≥L+1;其特征在于,所述联合定位与通信方法包括以下步骤:1. A multi-reflection device joint positioning and communication method based on environmental backscatter, comprising an ambient radio frequency source, L reflection devices, and a receiver with M antennas, and the number of the reflection devices is l=0,1,. ..L-1, the receiver coordinates are (x r , y r ), M≥L+1; it is characterized in that the joint positioning and communication method comprises the following steps: S1、环境射频源发射射频源信号,反射设备对环境射频源信号进行反向散射,其中,反射设备的波形设计方法为:S1. The ambient radio frequency source transmits the radio frequency source signal, and the reflection device backscatters the ambient radio frequency source signal. The waveform design method of the reflection device is as follows: 假设反射设备符号周期是OFDM射频源符号周期的K倍,每个反射设备发出的二元信息用Bl(m)表示,Bl(m)=±1,当传输码片中的某一个符号B(m)=-1时,在对应的第k个OFDM符号周期内,设计如下的反射设备波形:Assuming that the symbol period of the reflection device is K times the symbol period of the OFDM radio frequency source, the binary information sent by each reflection device is represented by B l (m), B l (m)=±1, when a certain symbol in the transmission chip is transmitted When B(m)=-1, in the corresponding kth OFDM symbol period, design the following reflection device waveform:
Figure FDA0003476261230000011
Figure FDA0003476261230000011
其中,Nc是OFDM射频源符号的循环前缀的长度,N是数据部分长度,N+Nc构成了一个完整的OFDM符号的长度,k=1,...K代表在同一个反射设备符号周期内的OFDM射频源符号的序号;Among them, N c is the length of the cyclic prefix of the OFDM radio frequency source symbol, N is the length of the data part, N+N c constitutes the length of a complete OFDM symbol, k=1,...K represents the symbol in the same reflection device The sequence number of the OFDM radio source symbol in the period; 当传输码片的某一个符号B(m)=1时,反射设备波形恒定为1;When a certain symbol B(m)=1 of the transmission chip, the waveform of the reflection device is always 1; S2、接收机接收信号,进行DOA估计,估计得到直射链路信号和多个反向散射信号的到达角θd和θL=[θ012,...,θl,...θL-1];S2. The receiver receives the signal, performs DOA estimation, and obtains the estimated angles of arrival of the direct link signal and multiple backscattered signals θ d and θ L =[θ 012 ,...,θ l , ...θ L-1 ]; S3、对所有估计出来的到达角的信号各执行一次空域滤波,每一次空域滤波时,设计权向量进行波束成形,处理接收的信号,只保留一个角度的信号,抑制其他方向信号的功率,获得所有到达角对应的方向信号;S3. Perform spatial filtering once for all the estimated angle of arrival signals. In each spatial filtering, design a weight vector to perform beamforming, process the received signals, keep only the signals of one angle, suppress the power of signals in other directions, and obtain Direction signals corresponding to all angles of arrival; S4、计算空域滤波后的信号z(n)的功率,设信号有N个采样点,则信号的功率Pz计算方法为:S4. Calculate the power of the signal z (n) after spatial filtering. Suppose the signal has N sampling points, and the calculation method of the power Pz of the signal is:
Figure FDA0003476261230000012
Figure FDA0003476261230000012
S5、获得所有到达角的方向信号对应的功率以后,根据反向散射信号相比较于直射链路信号会经历两次衰落的特性,区分出直射信号功率
Figure FDA0003476261230000021
和反射信号功率
Figure FDA0003476261230000022
通过直射信号功率
Figure FDA0003476261230000023
计算射频源到接收机的距离df
S5. After obtaining the power corresponding to the direction signals of all the angles of arrival, according to the characteristic that the backscattered signal will experience two fading compared with the direct link signal, distinguish the power of the direct signal
Figure FDA0003476261230000021
and reflected signal power
Figure FDA0003476261230000022
through direct signal power
Figure FDA0003476261230000023
Calculate the distance d f from the RF source to the receiver:
Figure FDA0003476261230000024
Figure FDA0003476261230000024
其中,n代表该环境的衰减系数,β是一个常数,由载波频率和环境决定,n和β均可在实际测量中获得,ps为发射功率;Among them, n represents the attenuation coefficient of the environment, β is a constant, determined by the carrier frequency and the environment, both n and β can be obtained in actual measurement, and ps is the transmit power; 其余的信号为反射链路信号,通过功率估计射频源到反射设备距离和反射设备到接收机距离的乘积D:The rest of the signals are reflected link signals, and the power is used to estimate the product D of the distance from the RF source to the reflecting device and the distance from the reflecting device to the receiver:
Figure FDA0003476261230000025
Figure FDA0003476261230000025
其中,
Figure FDA0003476261230000026
为环境射频源到第l个反射设备的距离,
Figure FDA0003476261230000027
为第l个反射设备到接收机的距离,αl为第l个反射设备的反射系数;
in,
Figure FDA0003476261230000026
is the distance from the ambient RF source to the lth reflecting device,
Figure FDA0003476261230000027
is the distance from the lth reflection device to the receiver, α l is the reflection coefficient of the lth reflection device;
S6、通过几何原理估计反射设备到接收机的距离,具体为根据三角形的余弦定理建立如下方程:S6. Estimate the distance from the reflection device to the receiver through geometric principles, specifically establishing the following equation according to the cosine theorem of the triangle:
Figure FDA0003476261230000028
Figure FDA0003476261230000028
其中,Δθ为直射径和反射径之间的夹角,Δθ=|θdl|,通过方程解得
Figure FDA0003476261230000029
Figure FDA00034762612300000210
Among them, Δθ is the angle between the direct radiation diameter and the reflected diameter, Δθ=|θ dl |, which can be obtained by solving the equation
Figure FDA0003476261230000029
and
Figure FDA00034762612300000210
S7、利用反向散射信号到达角信息θl和反射设备到接收机的距离信息
Figure FDA00034762612300000211
计算反射设备的坐标(xl,yl):
S7. Use the angle of arrival information θ l of the backscattered signal and the distance information from the reflection device to the receiver
Figure FDA00034762612300000211
Calculate the coordinates of the reflective device (x l , y l ):
Figure FDA00034762612300000212
Figure FDA00034762612300000212
实现定位;achieve positioning; S8、解调出反射信号的信息:S8, demodulate the information of the reflected signal: 假设K=1,即反射设备波形c(n)的一个周期的持续时长和一个OFDM符号持续时长相同,同时忽略噪声,在一个OFDM符号周期内,建立:Assuming K=1, that is, the duration of one cycle of the reflection device waveform c(n) is the same as the duration of one OFDM symbol, while ignoring noise, within one OFDM symbol cycle, establish:
Figure FDA0003476261230000031
Figure FDA0003476261230000031
其中,Lh是信道的时延扩展,代表被码间串扰污染的样点数,将这些样点丢弃;则有:Among them, L h is the delay spread of the channel, which represents the number of samples polluted by inter-symbol crosstalk, and these samples are discarded; then there are: u(m)=|Kl|2|s(n)|2B(m)u(m)=|K l | 2 |s(n)| 2 B(m) 其中,
Figure FDA0003476261230000032
是一个常数,pl=1,s(n)是射频源信号,解调时对B(m)判决法则如下:
in,
Figure FDA0003476261230000032
is a constant, p l =1, s(n) is the radio frequency source signal, the decision rule for B(m) during demodulation is as follows:
Figure FDA0003476261230000033
Figure FDA0003476261230000033
从而恢复出反射设备传输的码元;Thereby, the symbols transmitted by the reflection device are recovered; S9、使用已知的导频序列与解码后的反射信号的导频部分进行相关运算,当解调后的信号与某个已知导频序列进行相关运算得到的结果最大时,即确定信号与已知导频的对应关系,因此在传输信息的同时,判断出该信号属于哪一个反射设备。S9. Use the known pilot sequence to perform a correlation operation with the pilot part of the decoded reflected signal. When the result of the correlation operation between the demodulated signal and a certain known pilot sequence is the largest, it is determined that the signal is related to The corresponding relationship of the pilot frequency is known, so it is determined which reflection device the signal belongs to while transmitting the information.
CN202110159173.8A 2021-02-05 2021-02-05 Multi-reflection device joint positioning and communication method based on environment backscattering Active CN112954792B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110159173.8A CN112954792B (en) 2021-02-05 2021-02-05 Multi-reflection device joint positioning and communication method based on environment backscattering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110159173.8A CN112954792B (en) 2021-02-05 2021-02-05 Multi-reflection device joint positioning and communication method based on environment backscattering

Publications (2)

Publication Number Publication Date
CN112954792A CN112954792A (en) 2021-06-11
CN112954792B true CN112954792B (en) 2022-03-08

Family

ID=76242472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110159173.8A Active CN112954792B (en) 2021-02-05 2021-02-05 Multi-reflection device joint positioning and communication method based on environment backscattering

Country Status (1)

Country Link
CN (1) CN112954792B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116419149A (en) * 2021-12-31 2023-07-11 维沃移动通信有限公司 Positioning method, device and communication equipment
CN118542035A (en) * 2022-01-14 2024-08-23 上海诺基亚贝尔股份有限公司 Joint communication and awareness mechanism
CN115694686A (en) * 2022-09-27 2023-02-03 中国电信股份有限公司 Mixed signal processing method, device, apparatus, medium, and program
CN115616600A (en) * 2022-10-12 2023-01-17 中国电信股份有限公司 Positioning method, device, device, medium and program for backscattering equipment
CN116087875B (en) * 2023-04-07 2023-06-30 电子科技大学(深圳)高等研究院 Indoor non-line-of-sight scene passive tag three-dimensional positioning method based on environment information
CN117353804B (en) * 2023-12-05 2024-06-04 国网信息通信产业集团有限公司 A method and device for environmental backscatter communication

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6259991B1 (en) * 1999-02-10 2001-07-10 X-Cyte Inc. Environmental location system
CN106549692A (en) * 2016-12-12 2017-03-29 电子科技大学 A kind of signal acceptance method for backscatter communication system
CN106921415A (en) * 2017-03-09 2017-07-04 电子科技大学 A kind of signal acceptance method for environment backscatter communication system
CN107544054A (en) * 2017-08-15 2018-01-05 西京学院 A kind of indoor orientation method and device based on environment backscattering
CN107861100A (en) * 2017-11-10 2018-03-30 西京学院 A kind of indoor orientation method based on three-dimensional wave beam
CN109150253A (en) * 2018-08-07 2019-01-04 浙江大学 The information symbol detection method of environment backscatter system based on multiple antennas
CN110224965A (en) * 2019-06-17 2019-09-10 电子科技大学 One kind being based on half blind receiver design method of OFDM backscatter communication system
CN209486287U (en) * 2019-01-25 2019-10-11 西京学院 Indoor Positioning Device Based on Environmental Backscatter

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6259991B1 (en) * 1999-02-10 2001-07-10 X-Cyte Inc. Environmental location system
CN106549692A (en) * 2016-12-12 2017-03-29 电子科技大学 A kind of signal acceptance method for backscatter communication system
CN106921415A (en) * 2017-03-09 2017-07-04 电子科技大学 A kind of signal acceptance method for environment backscatter communication system
CN107544054A (en) * 2017-08-15 2018-01-05 西京学院 A kind of indoor orientation method and device based on environment backscattering
CN107861100A (en) * 2017-11-10 2018-03-30 西京学院 A kind of indoor orientation method based on three-dimensional wave beam
CN109150253A (en) * 2018-08-07 2019-01-04 浙江大学 The information symbol detection method of environment backscatter system based on multiple antennas
CN209486287U (en) * 2019-01-25 2019-10-11 西京学院 Indoor Positioning Device Based on Environmental Backscatter
CN110224965A (en) * 2019-06-17 2019-09-10 电子科技大学 One kind being based on half blind receiver design method of OFDM backscatter communication system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Gang Yang ; Ying-Chang Liang.Modulation in the Air: Backscatter Communication Over Ambient OFDM Carrier.《IEEE Transactions on Communications》.2017, *
Ying-Chang Liang ; Ruizhe Long.Cooperative Ambient Backscatter System: A Symbiotic Radio Paradigm for Passive IoT.《IEEE Wireless Communications Letters》.2019, *
反向散射通信技术与物联网;王公仆等;《物联网学报》;20170630(第01期);全文 *
智能无线通信技术研究概况;梁应敞;《通信学报》;20200720;全文 *

Also Published As

Publication number Publication date
CN112954792A (en) 2021-06-11

Similar Documents

Publication Publication Date Title
CN112954792B (en) Multi-reflection device joint positioning and communication method based on environment backscattering
CN112986905B (en) Multi-reflection equipment positioning method based on environmental backscattering
CN107132505B (en) Multi-target direct positioning method in direct and indirect mixed scene
CN106909779B (en) Calculation method of Cramero bound for MIMO radar based on distributed processing
CN104469941B (en) Indoor wireless positioning method based on WLAN WLAN ofdm signal cyclic prefix
CN105738865B (en) The direct localization method of multiple target under signal waveform known conditions
CN111123197A (en) TDOA-based target radiation source positioning method
CN109239683B (en) Clarmero boundary analysis method of broadband passive MIMO radar
CN108092926B (en) Parameter Estimation Method for Passive Backscatter Communication Channel
CN106950532B (en) Based on cross ambiguity function compared with wave-form similarity united delay time estimation method
CN102043143A (en) Simulation method for statistical MIMO (multiple input multiple output) radar target detection
CN104730517A (en) Bistatic MIMO radar multi-target tracking method
CN105939299A (en) Channel parameter estimation method based on improved SAGE algorithm
JP2018535585A (en) Method and computer-implemented method for decoding symbols transmitted via mm-wave channel, receiver for receiving and decoding symbols transmitted via mm-wave channel
CN110927660A (en) Mixed signal direction of arrival estimation method based on co-prime array
CN109768816B (en) non-Gaussian noise 3D-MIMO system data detection method
CN114415210A (en) A direct positioning method of satellite navigation deceptive interference source based on multi-dimensional information
CN105974362B (en) A kind of high-precision Passive Location of Combined estimator signal parameter and position
CN106125059B (en) Nonparametric Combined estimator signal and the Passive Location of position
CN108828505A (en) Angle-of- arrival estimation algorithm research and application based on machine learning
CN116405348A (en) Angle Estimation Method for Array Antennas and 5G New Radio Interface in Communication-Sensing Integrated System
CN114268902B (en) Pulse ultra-wideband direction finding method based on PDOA
CN112505624B (en) A Positioning System Based on Environmental Backscatter
CN117538821A (en) Channel estimation-based multi-path interference resistance method for underground long and narrow space of coal mine
CN113938172B (en) Method for aligning signal beam forming with passive sensor

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