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 PDFInfo
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating 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射频源下的低复杂度的解调反射信号的方案,无需解调射频源信息。即在对反射设备定位的同时,实现反射设备与接收机的通信。反射设备通过发送导频,接收端事先已知导频与反射设备的对应关系,通过接收端已知的导频和接收到的信号的导频进行相关运算,可以判决该反射信号和反射设备的对应关系,最终实现对某一确定的反射设备的信息解调和定位。
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.
Description
技术领域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:
其中,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=[θ0,θ1,θ2,...,θ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 =[θ 0 ,θ 1 ,θ 2 ,...,θ 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:
S5、获得所有到达角的方向信号对应的功率以后,根据反向散射信号相比较于直射链路信号会经历两次衰落的特性,区分出直射信号功率和反射信号功率通过直射信号功率计算射频源到接收机的距离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. and reflected signal power through direct signal power Calculate the distance d f from the RF source to the receiver:
其中,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:
其中,为环境射频源到第l个反射设备的距离,为第l个反射设备到接收机的距离,αl为第l个反射设备的反射系数;in, is the distance from the ambient RF source to the lth reflecting device, 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:
其中,Δθ为直射径和反射径之间的夹角,Δθ=|θd-θl|,通过方程解得和 Among them, Δθ is the angle between the direct radiation diameter and the reflected diameter, Δθ=|θ d -θ l |, which can be obtained by solving the equation and
S7、利用反向散射信号到达角信息θl和反射设备到接收机的距离信息计算反射设备的坐标(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 Calculate the coordinates of the reflective device (x l , y l ):
实现定位;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:
其中,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)
其中,是一个常数,pl=1,s(n)是射频源信号,解调时对B(m)判决法则如下:in, 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:
从而恢复出反射设备传输的码元;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;环境射频源发射射频源信号其中s(t)是功率归一化射频源基带信号,发射功率为ps。信道是强LOS径,图中代表环境射频源到反向散射体的信道参数,代表环境射频源到接收机的信道参数,代表反向散射体到接收机的信道参数。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 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 represents the channel parameters from the ambient RF source to the backscatterer, represents the channel parameters from the ambient RF source to the receiver, 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的反射设备的反射信号到达角为θl。Taking 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
其中l=0,1,...L-1。是带通噪声,是射频源到编号为l的反射设备的距离,是编号为l的反射设备到接收机的距离。cl(t)为编号为l的反射设备的要传输的基带信号,αl为反射系数。和为方向矢量,其中为空间相位,fc是载波频率,d是天线间距,c是光速,λ是载波波长。where l=0,1,...L-1. is the bandpass noise, is the distance from the RF source to the reflecting device numbered l, 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. and is the direction vector, where 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
其中ω(n)是复基带噪声,服从循环对称复高斯分布σ2是噪声功率。where ω(n) is the complex baseband noise, which follows a cyclic symmetric complex Gaussian distribution σ 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
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),其中 1. According to the set sampling number N, perform N sampling to obtain the sample sequence y(n), where
2.计算y(n)自相关矩阵R,R=E{y(n)yH(n)},这里可以用时间平均代替统计平均,即 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
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)]T。4. Construct the vector a(z), a(z)=[1, z -1 ,...,z- (M-1) ] T .
构造如下函数: Construct the following function:
乘以自身共轭转置,得到 Multiplying by the self-conjugate transpose, we get
i=3,4,...,M。 i=3,4,...,M.
5.定义多项式5. Define Polynomials
令其为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=ejφ,φ是空间角频率,φ=-πsinθ。因此7. Since z=e jφ , φ is the spatial angular frequency, φ=−πsinθ. therefore
计算得到信号到达方向估计。The calculation obtains an estimate of the direction of arrival of the signal.
得到直射信号到达角度估计值和反射信号到达角度估计值后,进行空域滤波。目的是只保留某一个方向的信号,将其他方向信号置零。Obtain the estimated value of the angle of arrival of the direct signal and the estimated angle of arrival of the reflected signal 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.
空域滤波的算法有很多种,在此举例一种空域滤波的算法,以保留方向的反射信号,将其他估计出来的方向的信号置零为例说明:There are many kinds of spatial filtering algorithms. Here is an example of a spatial filtering algorithm to preserve the The reflected signal of the direction, and the signal of other estimated directions is set to zero as an example to illustrate:
写出矩阵代表要迫零的方向的集合。提出优化问题:write out the matrix Represents a collection of directions to zero forcing. Ask the optimization problem:
s.t wHA=0st w H A = 0
具体算法步骤如下:The specific algorithm steps are as follows:
1.对想要迫零的方向,通过方向向量构造矩阵 1. For the direction you want to force zero, construct a matrix through the direction vector
2.求解BHA=0可通过AHB=0求解2. Solving B H A=0 can be solved by A H B=0
3.其中B+是B矩阵的M-P广义逆矩阵。3. where B + is the MP generalized inverse of the B matrix.
这里以只保留到达角为的反射信号为例,空域滤波完成后,信号形式为:Here, only the arrival angle is retained as 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:
其中,w是空域滤波过程中设计的权向量,p0=wHa(θ0),εl=wHa(θl)εd=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(θ l )ε d =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)信号的功率为当p0≈1,反射系数α0已知,n和β实现进行实地测量后已知的情况下,则可以估计出的值,即射频源到反射设备的距离和反射设备到接收机的距离的乘积D。The power of the measured z 0 (n) signal is When p 0 ≈ 1, the reflection coefficient α 0 is known, and n and β are known after field measurements, then it can be estimated that 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.
保留直射链路的信号,抑制所有反射信号以估计射频源到接收机的距离,按照上述类似的方法进行空域滤波,可以得到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
式中pd=wHa(θd),pd≈1,εl是很小的值,约为0,可以归为噪声。因此,可以求出信号zd(n)的功率计算出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 Calculate the value of d f to estimate the distance from the RF source to the receiver.
因此对反射信号的定位问题可以建模成如下的几何问题:Therefore, the positioning problem of the reflected signal can be modeled as the following geometric problem:
θd和θ0在执行DOA估计以后可以得到,而df通过只保留直射路径的信号可以估计出来,作为两条径距离的乘积也可以估计出来。那么由三角形的余弦定理。直射径和反射径之间的夹角为Δθ=|θd-θ0|,可以建立如下方程:θ d and θ 0 can be obtained after performing DOA estimation, and d f can be estimated by keeping only the direct path signal, 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 Δθ=|θ d -θ 0 |, and the following equation can be established:
该方程中,D和Δθ和df均已知,两个方程可以解得和这样,反射设备到达接收机的距离即可估计出来,同时已知信号到达角θ0,则可以估计出反射设备的位置。In this equation, D and Δθ and d f are known, and the two equations can be solved as and In this way, the distance from the reflecting device to the receiver 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 θ
最终得到该反射信号的定位结果。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:
可以将模型简化为如下形式: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)
其中是一个常数,和由于空域滤波已经置零,将其归入ω(n)项内。in is a constant, and 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:
在一个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
对该信号进行处理,令:To process the signal, let:
则ωcross(n)为共轭相乘运算中产生的交叉项,也归于噪声项。but ω 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:
如此即可恢复出反射设备传输的码元。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:
在该仿真条件下,可以知道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)
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)
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)
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 |
-
2021
- 2021-02-05 CN CN202110159173.8A patent/CN112954792B/en active Active
Patent Citations (8)
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)
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 |