CN113965881A - Millimeter wave integrated communication and sensing method under shielding effect - Google Patents

Millimeter wave integrated communication and sensing method under shielding effect Download PDF

Info

Publication number
CN113965881A
CN113965881A CN202111156208.9A CN202111156208A CN113965881A CN 113965881 A CN113965881 A CN 113965881A CN 202111156208 A CN202111156208 A CN 202111156208A CN 113965881 A CN113965881 A CN 113965881A
Authority
CN
China
Prior art keywords
occlusion
perception
pixel
environmental
channel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111156208.9A
Other languages
Chinese (zh)
Other versions
CN113965881B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202111156208.9A priority Critical patent/CN113965881B/en
Publication of CN113965881A publication Critical patent/CN113965881A/en
Application granted granted Critical
Publication of CN113965881B publication Critical patent/CN113965881B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/18Network planning tools

Landscapes

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

Abstract

The invention provides a millimeter wave integrated communication and perception method under the shielding effect, which utilizes pilot signals or other known data sequences of the existing communication system to carry out perception. And establishing a probabilistic reasoning model based on a factor graph according to the relation among the received data, the distribution of the environmental scatterers and the shielding effect. And finally, solving a compressed sensing reconstruction problem based on a bilinear approximate message transmission method utilizing the shielding relation, thereby realizing the sensing of the environment. Compared with the existing environment perception reconstruction algorithm, the millimeter wave environment perception algorithm based on the occlusion effect obviously improves the accuracy of environment perception, and provides an efficient environment perception method for the design of a future perception communication integrated system.

Description

遮挡效应下的毫米波一体化通信与感知方法Millimeter-wave integrated communication and perception method under occlusion effect

技术领域technical field

本发明涉及无线通信领域,尤其新一代无线通信领域中的感知与通信一体化系统设计。The invention relates to the field of wireless communication, in particular to the design of an integrated system of perception and communication in the field of new generation wireless communication.

背景技术Background technique

毫米波(mmWave)由于其具有高带宽、高可靠性和高集成性特点,结合大规模多输入多输出(MIMO)技术,已经成为了当前无线通信领域中的研究重点。随着无线通信行业的快速发展,考虑到快速增多的连接设备和服务,越发复杂的无线移动通信的应用场景,无线通信信号的传播环境日益复杂。同时,随着无线通信基站越发密集的部署,并且接收方和发送方具有相比以前更加强大的计算能力和物理性能,信息处理能力大幅提高,这对于通信体制提出了更高的要求。我们不仅希望它们能够继续提供原本的通信服务,同时希望通信设备利用自己强大的处理能力完成对环境的感知。具体而言,在未来的无线通信场景中,智慧城市、自动驾驶和无人机定位等新型技术不仅需要无线宽带连接,还需要感知准确的环境信息,包括但不限于该环境中静止或移动物体的背景散射体的位置、形状、状态和电磁特性等。Millimeter wave (mmWave) has become a research focus in the current wireless communication field due to its high bandwidth, high reliability, and high integration, combined with massive multiple-input multiple-output (MIMO) technology. With the rapid development of the wireless communication industry, considering the rapidly increasing number of connected devices and services, and the increasingly complex application scenarios of wireless mobile communication, the propagation environment of wireless communication signals is increasingly complex. At the same time, with the increasingly dense deployment of wireless communication base stations, and the receiver and sender have more powerful computing power and physical performance than before, the information processing capacity has been greatly improved, which puts forward higher requirements for the communication system. We not only hope that they can continue to provide original communication services, but also hope that communication devices can use their powerful processing capabilities to complete the perception of the environment. Specifically, in future wireless communication scenarios, new technologies such as smart cities, autonomous driving, and drone positioning require not only wireless broadband connections, but also perception of accurate environmental information, including but not limited to stationary or moving objects in the environment The location, shape, state, and electromagnetic properties of the background scatterers.

如何基于无线通信架构,利用无线通信设备进行环境感知,实现感知通信一体化是下一代无线通信系统的重要研究方向。其中,在利用上行链路进行环境感知时,基站通过接收用户发送的上行链路通信信号,进行数据通信的同时,实现被动式环境感知。感知通信一体化系统设计的一大挑战在于环境中潜在的大量未知变量,因此应该利用目标环境本身的稀疏性。例如,在蜂窝通信网络中,建筑物稀疏地分布在无线网络覆盖范围内。除此之外,散射体之间普遍具有遮挡效应,距离用户较近位置的散射体会对电磁信号产生遮挡,使其不能到达至同一传播方向下距离较远的散射体。因此,并不是感知范围内的所有散射体都会对同一用户的多径信道产生影响,互相遮挡造成了不同用户面临着不同的散射体环境,多个用户需要进行多视角的联合感知。目前,现有的环境感知算法没有利用环境散射体分布和遮挡效应之间关系,没有对此问题进行针对性的考虑,在求解具有遮挡效应的成像模型时性能较差。How to use wireless communication equipment for environmental perception based on the wireless communication architecture and realize the integration of perception and communication is an important research direction of the next generation wireless communication system. Wherein, when using the uplink for environment perception, the base station implements passive environment perception while performing data communication by receiving the uplink communication signal sent by the user. One of the major challenges in the design of sensing-communication integrated systems lies in the potentially large number of unknown variables in the environment, so the sparsity of the target environment itself should be exploited. For example, in a cellular communication network, buildings are sparsely distributed within the wireless network coverage. In addition, there is generally a blocking effect between scatterers, and the scatterers located closer to the user will block the electromagnetic signal, so that they cannot reach the scatterers farther away in the same propagation direction. Therefore, not all scatterers within the sensing range will affect the multipath channel of the same user, and mutual occlusion causes different users to face different scatterer environments, and multiple users need to perform joint sensing from multiple perspectives. At present, the existing environment perception algorithms do not take advantage of the relationship between the distribution of environmental scatterers and the occlusion effect, and do not take this problem into account, and the performance is poor when solving the imaging model with the occlusion effect.

综上所述,综合考虑毫米波环境感知问题和遮挡效应下的用户观测目标不一致的问题,如何联合实现上行数据中环境信息的分离和遮挡效应的分辨具有较高的研究难度和现实意义。To sum up, considering the millimeter wave environment perception problem and the inconsistency of user observation targets under the occlusion effect, how to jointly realize the separation of environmental information in the uplink data and the resolution of the occlusion effect has high research difficulty and practical significance.

发明内容SUMMARY OF THE INVENTION

本发明的目的是为了解决在上行链路无线通信场景中,基站如何利用多用户发送的上行数据进行环境感知的问题。本发明利用现有通信系统的导频信号或其它已知数据序列来进行感知,收发处理分离,可与现有通信系统兼容,实现感知通信一体化。考虑到由于环境散射体之间存在遮挡效应,不同用户面临着不同的散射体环境,提出了一种遮挡效应下的毫米波一体化通信与感知方法。The purpose of the present invention is to solve the problem of how a base station uses uplink data sent by multiple users to perform environment perception in an uplink wireless communication scenario. The present invention utilizes the pilot signal or other known data sequence of the existing communication system for sensing, separates the sending and receiving processing, is compatible with the existing communication system, and realizes the integration of sensing and communication. Considering that due to the occlusion effect between ambient scatterers, different users face different scatterer environments, an integrated communication and perception method of mmWave under occlusion effect is proposed.

本发明所采用的具体技术方案如下:The concrete technical scheme adopted in the present invention is as follows:

一种遮挡效应下的毫米波一体化通信与感知方法,包括如下步骤:A millimeter wave integrated communication and perception method under occlusion effect, comprising the following steps:

1)在第T个时隙中,基站接收空间内的所有活跃用户发送的长度为L的导频序列s信号;1) In the T-th time slot, the base station receives a pilot sequence s signal of length L sent by all active users in the space;

2)基站接收到信号后,基于一种考虑遮挡效应的毫米波信道模型,将环境感知问题转换为压缩感知重构问题;2) After the base station receives the signal, based on a millimeter wave channel model considering the occlusion effect, the environment perception problem is converted into a compressed sensing reconstruction problem;

3)基于环境散射体之间的相对位置关系,建立一种环境散射体分布与遮挡关系的模型;3) Based on the relative positional relationship between environmental scatterers, a model of the relationship between the distribution and occlusion of environmental scatterers is established;

4)计算接收数据、环境散射体分布和遮挡效应之间的关系,建立一种基于因子图的概率推理模型;4) Calculate the relationship between the received data, the distribution of environmental scatterers and the occlusion effect, and establish a probabilistic inference model based on factor graphs;

5)结合步骤3)和步骤4)中得到的模型,基于一种利用遮挡关系的双线性近似消息传递方法,求解步骤3)中的压缩感知重构问题,实现对环境的感知。5) Combining the models obtained in steps 3) and 4), based on a bilinear approximation message passing method using occlusion relations, the compressive sensing reconstruction problem in step 3) is solved to realize the perception of the environment.

在一个或多个实施例中,为了将步骤2)中环境感知问题转换为压缩感知重构问题,需要在基于一种考虑遮挡效应的毫米波信道模型的基础上,对由环境散射体所造成多径信道HS进行估计,并分离出信道中的环境信息。In one or more embodiments, in order to convert the environmental perception problem in step 2) into a compressive sensing reconstruction problem, it is necessary to analyze the problems caused by environmental scatterers on the basis of a millimeter-wave channel model considering the occlusion effect. The multipath channel H S is estimated and the environmental information in the channel is separated.

在一个或多个实施例中,步骤2)中的毫米波信道模型和压缩感知重构问题为:In one or more embodiments, the millimeter-wave channel model and compressive sensing reconstruction problem in step 2) are:

2.a)环境信息离散化,将整个空间内的环境信息视作为点云,点云中的每个点代表其周围大小为lr,wr和hr的小立方体的环境信息,这些小立方体称为像素;感知的散射体环境的长、宽、高分别为Lr、Wr和Hr,则空间内点云的数量为N=Lr/lr×Wr/wr×Hr/hr;每一个像素内部可能是空的,也可能是存在散射体;我们使用一个散射系数xn来表示第n个点云点的所在小立方体的散射系数,若小立方体内部是空的,则xn=0;所以整个空间的环境信息可以用如下变量表示:2.a) The environmental information is discretized, and the environmental information in the entire space is regarded as a point cloud. Each point in the point cloud represents the environmental information of small cubes with sizes l r , wr and hr around them. The cube is called a pixel; the length, width and height of the perceived scatterer environment are L r , W r and H r respectively, then the number of point clouds in the space is N=L r /l r ×W r /w r ×H r /h r ; each pixel may be empty or there may be scatterers; we use a scattering coefficient x n to represent the scattering coefficient of the small cube where the nth point cloud point is located, if the interior of the small cube is empty , then x n =0; so the environmental information of the entire space can be represented by the following variables:

x=[x1,x2,…,xN]Tx=[x 1 , x 2 , . . . , x N ] T ;

2.b)将系统模型表示如下,在空间内的多个用户共享时频资源,在任意频率资源块上,第nR个BS接收天线接收信号表示为2.b) The system model is expressed as follows, multiple users in space share time-frequency resources, and on any frequency resource block, the signal received by the n- th BS receiving antenna is expressed as

Figure BDA0003284405150000031
Figure BDA0003284405150000031

其中

Figure BDA0003284405150000032
表示NR个BS天线的接收信号,
Figure BDA0003284405150000033
表示Nu个用户发送长度为L码元的导频,w为噪声,
Figure BDA0003284405150000034
表示用户到接收天线直视路径的信道系数,
Figure BDA0003284405150000035
表示用户到空间点云位置的信道系数,
Figure BDA0003284405150000036
表示用户到空间点云位置的信道系数,且
Figure BDA0003284405150000037
其中
Figure BDA0003284405150000038
表示用户至空间点云位置的遮挡矩阵,当Vs1出现全零列时,表示所有用户都无法感知到该像素点,即该像素点在感知范围外,Vs2(nR)∈{0,1}N×1表示空间点云位置直射至接收天线的遮挡矩阵,当Vs2(nR)出现零元素时,表示接收天线无法感知到该像素点,该像素点也在感知范围外,且V(nR)=Vs1(nR)diag(Vs2(nR));in
Figure BDA0003284405150000032
represents the received signals of NR BS antennas,
Figure BDA0003284405150000033
Indicates that Nu users send pilots of length L symbols, w is noise,
Figure BDA0003284405150000034
represents the channel coefficient of the direct line of sight from the user to the receiving antenna,
Figure BDA0003284405150000035
represents the channel coefficient from the user to the position of the spatial point cloud,
Figure BDA0003284405150000036
represents the channel coefficient from the user to the spatial point cloud location, and
Figure BDA0003284405150000037
in
Figure BDA0003284405150000038
Represents the occlusion matrix from the user to the position of the spatial point cloud. When all zero columns appear in V s1 , it means that all users cannot perceive the pixel, that is, the pixel is outside the perception range, V s2 (n R )∈{0, 1} N×1 represents the occlusion matrix that the spatial point cloud position directly hits the receiving antenna. When a zero element appears in V s2 (n R ), it means that the receiving antenna cannot perceive the pixel, and the pixel is also outside the perception range, and V(n R )=V s1 (n R )diag(V s2 (n R ));

2.c)基于上述系统模型,求解环境信息的优化问题表达为,2.c) Based on the above system model, the optimization problem of solving environmental information is expressed as,

Figure BDA0003284405150000041
Figure BDA0003284405150000041

在通信过程中,信道不确定性主要来自于未知的环境散射体分布,对于直视信道HLOS采用已知的统计模型描述,由此用户发送已知的导频s对信道HS进行估计,从而将求解问题转换为压缩感知重构问题:In the communication process, the channel uncertainty mainly comes from the unknown distribution of environmental scatterers. For the direct-view channel H LOS , a known statistical model is used to describe it, and the user sends the known pilot frequency s to estimate the channel H S. This transforms the solving problem into a compressive sensing reconstruction problem:

Figure BDA0003284405150000042
Figure BDA0003284405150000042

其将NR个BS接收天线所得的信道估计值拼接为矩阵,其中HS已知,

Figure BDA0003284405150000043
H已知,需要在遮挡矩阵V未知的条件下求解环境信息x。It splices the channel estimates obtained by the N R BS receiving antennas into a matrix, where H S is known,
Figure BDA0003284405150000043
H is known, and the environment information x needs to be solved under the condition that the occlusion matrix V is unknown.

在一个或多个实施例中,步骤3)中的环境散射体分布与遮挡关系的模型为:In one or more embodiments, the model of the relationship between ambient scatterer distribution and occlusion in step 3) is:

像素点A和B之间存在一个存在散射系数的像素点C,我们提出了一种环境散射体分布与遮挡关系的模型进行遮挡判断,令像素点A位置为坐标原点,则像素点B和C的坐标分别表示为b和c,通过满足以下三点条件判断像素点C对像素点A和B之间的直视路径产生了遮挡;There is a pixel point C with a scattering coefficient between pixel points A and B. We propose a model of the relationship between the distribution of environmental scatterers and occlusion for occlusion judgment. Let the position of pixel point A be the coordinate origin, then pixel points B and C The coordinates of , are expressed as b and c, respectively, and it is judged that pixel C occludes the direct line of sight between pixels A and B by satisfying the following three conditions;

像素点C与像素点A和B之间的直视路径的距离d小于阈值β:The distance d of the direct path between pixel C and pixels A and B is less than the threshold β:

Figure BDA0003284405150000044
Figure BDA0003284405150000044

向量b与向量c之间的夹角为锐角:The angle between the vector b and the vector c is an acute angle:

b·c>0;b·c>0;

像素点C位于像素点A和B之间:Pixel C is located between pixels A and B:

||c·b||<||b||2||c·b||<||b|| 2 .

在一个或多个实施例中,步骤4)中的一种基于因子图的概率推理模型为:In one or more embodiments, a factor graph-based probabilistic inference model in step 4) is:

Figure BDA0003284405150000051
与x后验概率分解为:Will
Figure BDA0003284405150000051
The posterior probability with x is decomposed as:

Figure BDA0003284405150000052
Figure BDA0003284405150000052

根据步骤2)的压缩感知重构问题和步骤3)中的遮挡关系构建因子图模型;因子图包含变量节点

Figure BDA0003284405150000053
和x,函数节点
Figure BDA0003284405150000054
px(x)和
Figure BDA0003284405150000055
并令HS长度M=NuNR,其中:A factor graph model is constructed according to the compressed sensing reconstruction problem in step 2) and the occlusion relationship in step 3); the factor graph contains variable nodes
Figure BDA0003284405150000053
and x, the function node
Figure BDA0003284405150000054
p x (x) and
Figure BDA0003284405150000055
And let H S length M=N u N R , where:

Figure BDA0003284405150000056
Figure BDA0003284405150000056

其中,HS是z的加噪观测,噪声是方差为σw的高斯白噪声;px(x)表示x的先验分布为高斯-伯努利分布:where H S is the noised observation of z, and the noise is Gaussian white noise with variance σ w ; p x (x) indicates that the prior distribution of x is a Gauss-Bernoulli distribution:

px(x)=(1-λ)δ(x)+λN(x;θx,σx);p x (x)=(1-λ)δ(x)+λN(x; θ x , σ x );

其中λ表示环境信息x的稀疏度,

Figure BDA0003284405150000057
表示环境信息x对信道的遮挡关系,具体如下:where λ represents the sparsity of the environmental information x,
Figure BDA0003284405150000057
Represents the occlusion relationship of the environmental information x to the channel, as follows:

Figure BDA0003284405150000058
Figure BDA0003284405150000058

在一个或多个实施例中,步骤5)中的利用遮挡关系的双线性近似消息传递方法为:In one or more embodiments, the bilinear approximation message passing method utilizing the occlusion relationship in step 5) is:

5.a)初始化算法参数,令t为迭代次数,对于m=1,2,...,M,残差均值为

Figure BDA0003284405150000059
对于m=1,2,...,M,n=1,2,...,N,根据步骤4)中设定的px(x)和
Figure BDA0003284405150000061
生成x的估计的均值
Figure BDA0003284405150000062
和方差
Figure BDA0003284405150000063
的估计的均值
Figure BDA0003284405150000064
和方差
Figure BDA0003284405150000065
5.a) Initialize the algorithm parameters, let t be the number of iterations, and for m = 1, 2, ..., M, the residual mean is
Figure BDA0003284405150000059
For m = 1, 2, ..., M, n = 1, 2, ..., N, according to the p x (x) and
Figure BDA0003284405150000061
generate the estimated mean of x
Figure BDA0003284405150000062
and variance
Figure BDA0003284405150000063
the estimated mean of
Figure BDA0003284405150000064
and variance
Figure BDA0003284405150000065

5.b)对于m=1,2,...,M,计算

Figure BDA0003284405150000066
的观测均值
Figure BDA0003284405150000067
和方差
Figure BDA0003284405150000068
具体如下:5.b) For m=1,2,...,M, calculate
Figure BDA0003284405150000066
The observed mean of
Figure BDA0003284405150000067
and variance
Figure BDA0003284405150000068
details as follows:

Figure BDA0003284405150000069
Figure BDA0003284405150000069

Figure BDA00032844051500000610
Figure BDA00032844051500000610

其中,in,

Figure BDA00032844051500000611
Figure BDA00032844051500000611

Figure BDA00032844051500000612
Figure BDA00032844051500000612

5.c)对于m=1,2,...,M,计算zm的估计的均值

Figure BDA00032844051500000613
和方差
Figure BDA00032844051500000614
具体如下:5.c) For m = 1, 2, ..., M, compute the mean of the estimates of z m
Figure BDA00032844051500000613
and variance
Figure BDA00032844051500000614
details as follows:

Figure BDA00032844051500000615
Figure BDA00032844051500000615

Figure BDA00032844051500000616
Figure BDA00032844051500000616

其中,in,

Figure BDA00032844051500000617
Figure BDA00032844051500000617

其中,C是归一化变量;where C is the normalization variable;

5.d)对于m=1,2,...,M,计算残差的均值

Figure BDA00032844051500000618
和方差
Figure BDA00032844051500000619
具体如下:5.d) For m = 1, 2, ..., M, calculate the mean of the residuals
Figure BDA00032844051500000618
and variance
Figure BDA00032844051500000619
details as follows:

Figure BDA00032844051500000620
Figure BDA00032844051500000620

Figure BDA00032844051500000621
Figure BDA00032844051500000621

5.e)对于m=1,2,...,M,n=1,2,...,N,计算

Figure BDA0003284405150000071
的观测均值
Figure BDA0003284405150000072
和方差
Figure BDA0003284405150000073
对于n=1,2,...,N,计算
Figure BDA0003284405150000074
的观测均值
Figure BDA0003284405150000075
和方差
Figure BDA0003284405150000076
具体如下:5.e) For m=1,2,...,M,n=1,2,...,N, compute
Figure BDA0003284405150000071
The observed mean of
Figure BDA0003284405150000072
and variance
Figure BDA0003284405150000073
For n=1,2,...,N, compute
Figure BDA0003284405150000074
The observed mean of
Figure BDA0003284405150000075
and variance
Figure BDA0003284405150000076
details as follows:

Figure BDA0003284405150000077
Figure BDA0003284405150000077

Figure BDA0003284405150000078
Figure BDA0003284405150000078

Figure BDA0003284405150000079
Figure BDA0003284405150000079

Figure BDA00032844051500000710
Figure BDA00032844051500000710

5.f)对于m=1,2,...,M,n=1,2,...,N,计算

Figure BDA00032844051500000711
的估计均值
Figure BDA00032844051500000712
和方差
Figure BDA00032844051500000723
对于n=1,2,...,N,计算xn的观测均值
Figure BDA00032844051500000713
和方差
Figure BDA00032844051500000714
具体如下:5.f) For m=1,2,...,M,n=1,2,...,N, calculate
Figure BDA00032844051500000711
the estimated mean of
Figure BDA00032844051500000712
and variance
Figure BDA00032844051500000723
For n = 1, 2, ..., N, compute the observed mean of x n
Figure BDA00032844051500000713
and variance
Figure BDA00032844051500000714
details as follows:

Figure BDA00032844051500000715
Figure BDA00032844051500000715

Figure BDA00032844051500000716
Figure BDA00032844051500000716

Figure BDA00032844051500000717
Figure BDA00032844051500000717

Figure BDA00032844051500000718
Figure BDA00032844051500000718

其中,in,

Figure BDA00032844051500000719
Figure BDA00032844051500000719

Figure BDA00032844051500000720
Figure BDA00032844051500000720

5.g)重复执行步骤b)至步骤f)直到达到收敛条件

Figure BDA00032844051500000721
得到环境信息x的估计值
Figure BDA00032844051500000722
5.g) Repeat step b) to step f) until the convergence condition is reached
Figure BDA00032844051500000721
Get an estimate of the environmental information x
Figure BDA00032844051500000722

本发明具有的有益效果是:在无线通信上行链路中,利用用户发送数据做环境感知的场景下(例如,一个多天线基站利用多个单天线用户发送的上行数据进行环境感知),本发明提出的遮挡效应下的毫米波环境感知方法,即一种遮挡效应下的毫米波一体化通信与感知方法,利用现有通信系统的导频信号或其它已知数据序列来进行感知,可与现有通信系统兼容,实现感知通信一体化。首先,本发明构建了一种考虑遮挡效应的毫米波信道模型,将环境感知问题转换为压缩感知重构问题。随后基于环境散射体之间的相对位置关系,建立一种环境散射体分布与遮挡关系的模型。并根据接收数据、环境散射体分布和遮挡效应之间的关系,建立一种基于因子图的概率推理模型。最后基于一种利用遮挡关系的双线性近似消息传递方法,求解压缩感知重构问题,实现对环境的感知。其利用多个用户观察视角,对环境进行了多视角的重构,相比现有的环境感知重构算法,本发明的基于遮挡效应的毫米波环境感知算法显著提升了环境感知的准确度,为未来感知通信一体化系统设计提供了一种高效的环境感知方法。The beneficial effects of the present invention are: in the uplink of wireless communication, in a scenario where data sent by users is used for environmental perception (for example, a multi-antenna base station uses uplink data sent by multiple single-antenna users for environmental perception), the present invention has the following advantages: The proposed millimeter-wave environment perception method under the occlusion effect, that is, a millimeter-wave integrated communication and perception method under the occlusion effect, uses the pilot signal of the existing communication system or other known data sequences for perception, which can be compared with the existing communication system. Compatible with communication systems to realize the integration of perception and communication. First, the present invention constructs a millimeter wave channel model considering the occlusion effect, and converts the environmental perception problem into a compressed sensing reconstruction problem. Then, based on the relative positional relationship between ambient scatterers, a model of the relationship between ambient scatterer distribution and occlusion is established. And according to the relationship between the received data, the distribution of environmental scatterers and the occlusion effect, a probabilistic inference model based on factor graph is established. Finally, based on a bilinear approximate message passing method using the occlusion relationship, the reconstruction problem of compressed sensing is solved to realize the perception of the environment. It uses multiple viewing angles of users to reconstruct the environment from multiple viewing angles. Compared with the existing environment perception reconstruction algorithm, the millimeter wave environment perception algorithm based on the occlusion effect of the present invention significantly improves the accuracy of environment perception, It provides an efficient environment perception method for the design of the future integrated system of perception and communication.

附图说明Description of drawings

图1是遮挡效应下的毫米波环境感知场景示意图;Figure 1 is a schematic diagram of a millimeter-wave environment perception scene under the occlusion effect;

图2是环境散射体分布与遮挡关系示意图;Figure 2 is a schematic diagram of the relationship between environmental scatterer distribution and occlusion;

图3基于遮挡关系的因子图;Fig. 3 factor graph based on occlusion relationship;

图4是将本发明与其他压缩感知重构算法相比较的环境感知直观结果图;Fig. 4 is the environment perception intuitive result diagram comparing the present invention with other compressed sensing reconstruction algorithms;

图5是将本发明与其他压缩感知重构算法相比较时,用户数量与环境感知准确度MSE的关系图;Fig. 5 is when comparing the present invention with other compressed sensing reconstruction algorithms, the relationship diagram of the number of users and the accuracy of environmental perception MSE;

图6是将本发明与其他压缩感知重构算法相比较时,信噪比SNR与环境感知准确度MSE的关系图。FIG. 6 is a graph showing the relationship between the signal-to-noise ratio (SNR) and the environment perception accuracy (MSE) when the present invention is compared with other compressed sensing reconstruction algorithms.

具体实施方式Detailed ways

如图1所示,首先我们考虑如下场景,在室外区域中部署有一个基站(BS),存在多个活跃用户(UE)。上行通信场景中,多个单天线用户同时发送上行通信信号至BS。发送信号由散射体散射,多径传播至BS接收。环境中散射体之间存在遮挡效应,如上图所示,用户1的发送信号仅由目标散射体1、3散射后被AP接收,目标散射体2则不影响这一过程。As shown in Figure 1, first we consider the following scenario, where a base station (BS) is deployed in an outdoor area and there are multiple active users (UE). In an uplink communication scenario, multiple single-antenna users send uplink communication signals to the BS simultaneously. The transmitted signal is scattered by the scatterer and multipath propagated to the BS for reception. There is an occlusion effect between scatterers in the environment. As shown in the figure above, the signal sent by user 1 is only scattered by target scatterers 1 and 3 and then received by the AP, while target scatterer 2 does not affect this process.

本发明一实施例提供了一种遮挡效应下的毫米波一体化通信与感知方法,其包括如下步骤:An embodiment of the present invention provides a millimeter-wave integrated communication and sensing method under occlusion effect, which includes the following steps:

1)在第T个时隙中,基站接收空间内的所有活跃用户发送的长度为L的导频序列s信号;1) In the T-th time slot, the base station receives a pilot sequence s signal of length L sent by all active users in the space;

2)基站接收到信号后,基于一种考虑遮挡效应的毫米波信道模型,将环境感知问题转换为压缩感知重构问题;2) After the base station receives the signal, based on a millimeter wave channel model considering the occlusion effect, the environment perception problem is converted into a compressed sensing reconstruction problem;

其中,在本发明一实施例中,为了将步骤2)中环境感知问题转换为压缩感知重构问题,需要在基于一种考虑遮挡效应的毫米波信道模型的基础上,对由环境散射体所造成多径信道HS进行估计,并分离出信道中的环境信息。Among them, in an embodiment of the present invention, in order to convert the environmental perception problem in step 2) into a compressive sensing reconstruction problem, it is necessary to use a millimeter wave channel model considering the occlusion effect. The multipath channel H S is estimated, and the environmental information in the channel is separated.

具体地,在本发明一实施例中,本步骤2)中所述的毫米波信道模型和压缩感知重构问题为:Specifically, in an embodiment of the present invention, the millimeter wave channel model and the compressive sensing reconstruction problem described in step 2) are:

2.a)将环境信息离散化,将整个空间内的环境信息视作为点云,点云中的每个点代表其周围大小为lr,wr和hr的小立方体的环境信息,这些小立方体称为像素。感知的散射体环境的长宽高分别为Lr,Wr和Hr,则空间内点云的数量为N=Lr/lr×Wr/wr×Hr/hr。每一个像素内部可能是空的,也可能是存在散射体。我们使用一个散射系数xn来表示第n个点云点的所在小立方体的散射系数,若小立方体内部是空的,则xn=0。所以整个空间的环境信息可以用如下变量表示:2.a) Discretize the environmental information, regard the environmental information in the whole space as a point cloud, each point in the point cloud represents the environmental information of the small cubes around it with sizes l r , wr and hr , these Small cubes are called pixels. The length, width and height of the perceived scatterer environment are L r , W r and H r respectively, then the number of point clouds in the space is N= L r /l r ×W r / wr ×H r /hr . The interior of each pixel may be empty, or there may be scatterers. We use a scattering coefficient x n to represent the scattering coefficient of the small cube where the nth point cloud point is located. If the interior of the small cube is empty, then x n =0. Therefore, the environmental information of the entire space can be represented by the following variables:

x=[x1,x2,…,xN]Tx=[x 1 , x 2 , . . . , x N ] T ;

2.b)将系统模型表示如下。在空间内的多个用户共享时频资源,在任意频率资源块上,第nR个BS接收天线接收信号表示为:2.b) Express the system model as follows. Multiple users in the space share time-frequency resources, and on any frequency resource block, the signal received by the n Rth BS receiving antenna is expressed as:

Figure BDA0003284405150000091
Figure BDA0003284405150000091

其中

Figure BDA0003284405150000101
表示NR个BS天线的接收信号,
Figure BDA0003284405150000102
表示Nu个用户发送长度为L码元的导频,w为噪声。
Figure BDA0003284405150000103
表示用户到接收天线直视路径的信道系数,
Figure BDA0003284405150000104
表示用户到空间点云位置的信道系数,
Figure BDA0003284405150000105
表示用户到空间点云位置的信道系数,且
Figure BDA0003284405150000106
其中
Figure BDA0003284405150000107
表示用户至空间点云位置的遮挡矩阵,当Vs1出现全零列时,表示所有用户都无法感知到该像素点,即该像素点在感知范围外。Vs2(nR)∈{0,1}N×1表示空间点云位置直射至接收天线的遮挡矩阵,当Vs2(nR)出现零元素时,表示接收天线无法感知到该像素点,该像素点也在感知范围外,且V(nR)=Vs1(nR)diag(Vs2(nR))。in
Figure BDA0003284405150000101
represents the received signals of NR BS antennas,
Figure BDA0003284405150000102
Indicates that Nu users send pilots with a length of L symbols, and w is noise.
Figure BDA0003284405150000103
represents the channel coefficient of the direct line of sight from the user to the receiving antenna,
Figure BDA0003284405150000104
represents the channel coefficient from the user to the position of the spatial point cloud,
Figure BDA0003284405150000105
represents the channel coefficient from the user to the spatial point cloud location, and
Figure BDA0003284405150000106
in
Figure BDA0003284405150000107
It represents the occlusion matrix from the user to the position of the spatial point cloud. When all zero columns appear in V s1 , it means that all users cannot perceive the pixel, that is, the pixel is outside the perception range. V s2 (n R )∈{0,1} N×1 represents the occlusion matrix that the spatial point cloud position directly hits the receiving antenna. When a zero element appears in V s2 (n R ), it means that the receiving antenna cannot perceive the pixel, This pixel is also outside the perception range, and V(n R )=V s1 (n R )diag(V s2 (n R )).

2.c)基于上述系统模型,求解环境信息的优化问题表达为,2.c) Based on the above system model, the optimization problem of solving environmental information is expressed as,

Figure BDA0003284405150000108
Figure BDA0003284405150000108

在通信过程中,信道不确定性主要来自于未知的环境散射体分布,对于直视信道HLOS采用已知的统计模型描述。由此用户发送已知的导频s对信道HS进行估计,从而将求解问题转换为压缩感知重构问题,即:In the communication process, the channel uncertainty mainly comes from the unknown distribution of environmental scatterers, and a known statistical model is used to describe the direct-view channel H LOS . Therefore, the user sends the known pilot frequency s to estimate the channel H S , thus transforming the solving problem into a compressed sensing reconstruction problem, namely:

Figure BDA0003284405150000109
Figure BDA0003284405150000109

其将NR个BS接收天线所得的信道估计值拼接为矩阵,其中HS已知,

Figure BDA00032844051500001010
H已知。需要在遮挡矩阵V未知的条件下求解环境信息x。It splices the channel estimates obtained by the N R BS receiving antennas into a matrix, where H S is known,
Figure BDA00032844051500001010
H is known. The environment information x needs to be solved under the condition that the occlusion matrix V is unknown.

3)基于环境散射体之间的相对位置关系,建立一种环境散射体分布与遮挡关系的模型;3) Based on the relative positional relationship between environmental scatterers, a model of the relationship between the distribution and occlusion of environmental scatterers is established;

具体地,在本发明一实施例中,本步骤3)中的环境散射体分布与遮挡关系的模型为:Specifically, in an embodiment of the present invention, the model of the relationship between ambient scatterer distribution and occlusion in step 3) is:

如图2所示,像素点A和B之间存在一个存在散射系数的像素点C,我们提出了一种环境散射体分布与遮挡关系的模型进行遮挡判断。令像素点A位置为坐标原点,则像素点B和C的坐标分别表示为b和c。通过满足以下三点条件判断像素点C对像素点A和B之间的直视路径产生了遮挡。As shown in Figure 2, there is a pixel C with a scattering coefficient between pixels A and B. We propose a model of the relationship between the distribution of ambient scatterers and occlusion for occlusion judgment. Let the position of pixel A be the origin of coordinates, then the coordinates of pixels B and C are represented as b and c, respectively. By satisfying the following three conditions, it is determined that the pixel point C occludes the direct view path between the pixel points A and B.

像素点C与像素点A和B之间的直视路径的距离d小于阈值β,The distance d of the direct path between the pixel point C and the pixel points A and B is less than the threshold β,

Figure BDA0003284405150000111
Figure BDA0003284405150000111

向量b与向量c之间的夹角为锐角,The angle between the vector b and the vector c is an acute angle,

b·c>0;b·c>0;

像素点C位于像素点A和B之间,Pixel C is located between pixels A and B,

||c·b||<||b||2||c·b||<||b|| 2 ;

4)计算接收数据、环境散射体分布和遮挡效应之间的关系,建立一种基于因子图的概率推理模型。4) Calculate the relationship between the received data, the distribution of environmental scatterers and the occlusion effect, and establish a probabilistic inference model based on factor graphs.

具体地,在本发明一实施例中,本步骤4)中的一种基于因子图的概率推理模型具体如下。Specifically, in an embodiment of the present invention, a factor graph-based probabilistic inference model in step 4) is as follows.

Figure BDA0003284405150000112
与x后验概率分解为,Will
Figure BDA0003284405150000112
The posterior probability with x is decomposed into,

Figure BDA0003284405150000113
Figure BDA0003284405150000113

如图3所示,根据步骤2)的压缩感知重构问题和步骤3)中的遮挡关系构建因子图模型。因子图包含变量节点

Figure BDA0003284405150000114
和x,函数节点
Figure BDA0003284405150000115
px(x)和
Figure BDA0003284405150000116
为简洁表达,令HS长度M=NuNR。As shown in Figure 3, a factor graph model is constructed according to the compressed sensing reconstruction problem in step 2) and the occlusion relationship in step 3). Factor graph contains variable nodes
Figure BDA0003284405150000114
and x, the function node
Figure BDA0003284405150000115
p x (x) and
Figure BDA0003284405150000116
For brevity, let HS length M=N u N R .

其中in

Figure BDA0003284405150000117
Figure BDA0003284405150000117

表示HS是z的加噪观测,噪声是方差为σw的高斯白噪声。px(x)表示x的先验分布为高斯-伯努利分布,Denote that H S is the noised observation of z, and the noise is white Gaussian noise with variance σ w . p x (x) indicates that the prior distribution of x is a Gauss-Bernoulli distribution,

px(x)=(1-λ)δ(x)+λN(x;θx,σx);p x (x)=(1-λ)δ(x)+λN(x; θ x , σ x );

其中λ表示环境信息x的稀疏度。

Figure BDA0003284405150000121
表示环境信息x对信道的遮挡关系,where λ represents the sparsity of the environmental information x.
Figure BDA0003284405150000121
represents the occlusion relationship of the environmental information x to the channel,

Figure BDA0003284405150000122
Figure BDA0003284405150000122

5)结合步骤3)和步骤4)中得到的模型,基于一种利用遮挡关系的双线性近似消息传递方法,求解步骤3)中的压缩感知重构问题,实现对环境的感知。5) Combining the models obtained in steps 3) and 4), based on a bilinear approximation message passing method using occlusion relations, the compressive sensing reconstruction problem in step 3) is solved to realize the perception of the environment.

具体地,在本发明一实施例中,本步骤中的利用遮挡关系的双线性近似消息传递方法为:Specifically, in an embodiment of the present invention, the bilinear approximation message passing method using the occlusion relationship in this step is:

5.a)初始化算法参数。令t为迭代次数,对于m=1,2,...,M,残差均值为

Figure BDA0003284405150000123
对于m=1,2,...,M,n=1,2,...,N,根据步骤4)中设定的px(x)和
Figure BDA0003284405150000124
生成x的估计的均值
Figure BDA0003284405150000125
和方差
Figure BDA0003284405150000126
的估计的均值
Figure BDA0003284405150000127
和方差
Figure BDA0003284405150000128
5.a) Initialize the algorithm parameters. Let t be the number of iterations, and for m = 1, 2, ..., M, the residual mean is
Figure BDA0003284405150000123
For m = 1, 2, ..., M, n = 1, 2, ..., N, according to the p x (x) and
Figure BDA0003284405150000124
generate the estimated mean of x
Figure BDA0003284405150000125
and variance
Figure BDA0003284405150000126
the estimated mean of
Figure BDA0003284405150000127
and variance
Figure BDA0003284405150000128

5.b)对于m=1,2,...,M,计算

Figure BDA0003284405150000129
的观测均值
Figure BDA00032844051500001210
和方差
Figure BDA00032844051500001211
5.b) For m=1,2,...,M, calculate
Figure BDA0003284405150000129
The observed mean of
Figure BDA00032844051500001210
and variance
Figure BDA00032844051500001211

Figure BDA00032844051500001212
Figure BDA00032844051500001212

Figure BDA00032844051500001213
Figure BDA00032844051500001213

其中,in,

Figure BDA00032844051500001214
Figure BDA00032844051500001214

Figure BDA00032844051500001215
Figure BDA00032844051500001215

5.c)对于m=1,2,...,M,计算zm的估计的均值

Figure BDA00032844051500001216
和方差
Figure BDA00032844051500001217
5.c) For m = 1, 2, ..., M, compute the mean of the estimates of z m
Figure BDA00032844051500001216
and variance
Figure BDA00032844051500001217

Figure BDA0003284405150000131
Figure BDA0003284405150000131

Figure BDA0003284405150000132
Figure BDA0003284405150000132

其中,in,

Figure BDA0003284405150000133
Figure BDA0003284405150000133

C是归一化变量。C is the normalization variable.

5.d)对于m=1,2,...,M,计算残差的均值

Figure BDA0003284405150000134
和方差
Figure BDA0003284405150000135
5.d) For m = 1, 2, ..., M, calculate the mean of the residuals
Figure BDA0003284405150000134
and variance
Figure BDA0003284405150000135

Figure BDA0003284405150000136
Figure BDA0003284405150000136

Figure BDA0003284405150000137
Figure BDA0003284405150000137

5.e)对于m=1,2,...,M,n=1,2,...,N,计算

Figure BDA0003284405150000138
的观测均值
Figure BDA0003284405150000139
和方差
Figure BDA00032844051500001310
对于n=1,2,...,N,计算
Figure BDA00032844051500001311
的观测均值
Figure BDA00032844051500001312
和方差
Figure BDA00032844051500001313
5.e) For m=1,2,...,M,n=1,2,...,N, compute
Figure BDA0003284405150000138
The observed mean of
Figure BDA0003284405150000139
and variance
Figure BDA00032844051500001310
For n=1,2,...,N, compute
Figure BDA00032844051500001311
The observed mean of
Figure BDA00032844051500001312
and variance
Figure BDA00032844051500001313

Figure BDA00032844051500001314
Figure BDA00032844051500001314

Figure BDA00032844051500001315
Figure BDA00032844051500001315

Figure BDA00032844051500001316
Figure BDA00032844051500001316

Figure BDA00032844051500001317
Figure BDA00032844051500001317

5.f)对于m=1,2,...,M,n=1,2,...,N,计算

Figure BDA00032844051500001318
的估计均值
Figure BDA00032844051500001319
和方差
Figure BDA00032844051500001320
对于n=1,2,...,N,计算xn的观测均值
Figure BDA00032844051500001321
和方差
Figure BDA00032844051500001322
5.f) For m=1,2,...,M,n=1,2,...,N, calculate
Figure BDA00032844051500001318
the estimated mean of
Figure BDA00032844051500001319
and variance
Figure BDA00032844051500001320
For n = 1, 2, ..., N, compute the observed mean of x n
Figure BDA00032844051500001321
and variance
Figure BDA00032844051500001322

Figure BDA00032844051500001323
Figure BDA00032844051500001323

Figure BDA00032844051500001324
Figure BDA00032844051500001324

Figure BDA00032844051500001325
Figure BDA00032844051500001325

Figure BDA0003284405150000141
Figure BDA0003284405150000141

其中in

Figure BDA0003284405150000142
Figure BDA0003284405150000142

Figure BDA0003284405150000143
Figure BDA0003284405150000143

5.g)重复执行步骤b)至步骤f)直到达到收敛条件

Figure BDA0003284405150000144
得到环境信息x的估计值
Figure BDA0003284405150000145
5.g) Repeat step b) to step f) until the convergence condition is reached
Figure BDA0003284405150000144
Get an estimate of the environmental information x
Figure BDA0003284405150000145

通过计算机仿真可以看出:如图4所示,其中GAMP算法完全忽视遮挡效应。Bilinear GAMP算法将遮挡矩阵V和环境信息x视作为独立变量,不利用遮挡关系分别进行求解。相比前两种算法,本发明的基于遮挡效应的毫米波环境感知算法显著提升了环境感知的准确度。图5表明随着用户数量的增加,本发明的方法环境感知效果逐渐提升并优于现有的算法。图5表明随着信噪比的增加,本发明的方法环境感知效果逐渐提升并优于现有的算法。It can be seen from the computer simulation that as shown in Figure 4, the GAMP algorithm completely ignores the occlusion effect. The Bilinear GAMP algorithm regards the occlusion matrix V and the environmental information x as independent variables, and does not use the occlusion relationship to solve them separately. Compared with the first two algorithms, the millimeter wave environment perception algorithm based on the occlusion effect of the present invention significantly improves the accuracy of environment perception. FIG. 5 shows that with the increase of the number of users, the environment perception effect of the method of the present invention is gradually improved and is superior to the existing algorithm. FIG. 5 shows that with the increase of the signal-to-noise ratio, the environmental perception effect of the method of the present invention is gradually improved and is superior to the existing algorithm.

综上所述,在无线通信上行链路中,利用用户发送数据做环境感知的场景下(例如,一个多天线基站利用多个单天线用户发送的上行数据进行环境感知),本发明实施例中提出的遮挡效应下的毫米波环境感知方法,即一种遮挡效应下的毫米波一体化通信与感知方法,利用现有通信系统的导频信号或其它已知数据序列来进行感知,可与现有通信系统兼容,实现感知通信一体化。首先,本发明实施例中构建了一种考虑遮挡效应的毫米波信道模型,将环境感知问题转换为压缩感知重构问题。随后基于环境散射体之间的相对位置关系,建立一种环境散射体分布与遮挡关系的模型。并根据接收数据、环境散射体分布和遮挡效应之间的关系,建立一种基于因子图的概率推理模型。最后基于一种利用遮挡关系的双线性近似消息传递方法,求解压缩感知重构问题,实现对环境的感知。其利用多个用户观察视角,对环境进行了多视角重构,相比现有的环境感知重构算法,本发明实施例中的基于遮挡效应的毫米波环境感知算法显著提升了环境感知的准确度,为未来感知通信一体化系统设计提供了一种高效的环境感知方法。To sum up, in the uplink of wireless communication, in a scenario where data sent by users is used for environmental perception (for example, a multi-antenna base station uses uplink data sent by multiple single-antenna users for environmental perception), in this embodiment of the present invention, The proposed millimeter-wave environment perception method under the occlusion effect, that is, a millimeter-wave integrated communication and perception method under the occlusion effect, uses the pilot signal of the existing communication system or other known data sequences for perception, which can be compared with the existing communication system. Compatible with communication systems to realize the integration of perception and communication. First, in the embodiment of the present invention, a millimeter-wave channel model considering the occlusion effect is constructed, and the environment perception problem is converted into a compressed sensing reconstruction problem. Then, based on the relative positional relationship between ambient scatterers, a model of the relationship between ambient scatterer distribution and occlusion is established. And according to the relationship between the received data, the distribution of environmental scatterers and the occlusion effect, a probabilistic inference model based on factor graph is established. Finally, based on a bilinear approximate message passing method using the occlusion relationship, the reconstruction problem of compressed sensing is solved to realize the perception of the environment. It uses multiple user viewing angles to reconstruct the environment from multiple perspectives. Compared with the existing environment perception reconstruction algorithm, the millimeter wave environment perception algorithm based on the occlusion effect in the embodiment of the present invention significantly improves the accuracy of environment perception. It provides an efficient environment perception method for the design of the integrated system of perception and communication in the future.

上述实施例用来解释说明本发明,而不是对本发明进行限制。在本发明的精神和权利要求的保护范围内,对本发明做出任何修改和改变,都落入本发明的保护范围。The above-mentioned embodiments are used to illustrate the present invention, but not to limit the present invention. Any modification or change made to the present invention within the spirit of the present invention and the protection scope of the claims shall fall within the protection scope of the present invention.

Claims (6)

1.一种遮挡效应下的毫米波一体化通信与感知方法,其特征在于,包括如下步骤:1. a millimeter-wave integrated communication and perception method under occlusion effect, is characterized in that, comprises the steps: 1)在第T个时隙中,基站接收空间内的所有活跃用户发送的长度为L的导频序列s信号;1) In the T-th time slot, the base station receives a pilot sequence s signal of length L sent by all active users in the space; 2)基站接收到信号后,基于一种考虑遮挡效应的毫米波信道模型,将环境感知问题转换为压缩感知重构问题;2) After the base station receives the signal, based on a millimeter wave channel model considering the occlusion effect, the environment perception problem is converted into a compressed sensing reconstruction problem; 3)基于环境散射体之间的相对位置关系,建立一种环境散射体分布与遮挡关系的模型;3) Based on the relative positional relationship between environmental scatterers, a model of the relationship between the distribution and occlusion of environmental scatterers is established; 4)计算接收数据、环境散射体分布和遮挡效应之间的关系,建立一种基于因子图的概率推理模型;4) Calculate the relationship between the received data, the distribution of environmental scatterers and the occlusion effect, and establish a probabilistic inference model based on factor graphs; 5)结合步骤3)和步骤4)中得到的模型,基于一种利用遮挡关系的双线性近似消息传递方法,求解步骤3)中的压缩感知重构问题,实现对环境的感知。5) Combining the models obtained in steps 3) and 4), based on a bilinear approximation message passing method using occlusion relations, the compressive sensing reconstruction problem in step 3) is solved to realize the perception of the environment. 2.根据权利要求1所述的遮挡效应下的毫米波一体化通信与感知方法,其特征在于,为了将步骤2)中环境感知问题转换为压缩感知重构问题,需要在基于一种考虑遮挡效应的毫米波信道模型的基础上,对由环境散射体所造成多径信道HS进行估计,并分离出信道中的环境信息。2. the millimeter-wave integrated communication and perception method under the occlusion effect according to claim 1, is characterized in that, in order to convert the environmental perception problem in step 2) into the compressive sensing reconstruction problem, need to consider occlusion based on a kind of On the basis of the millimeter-wave channel model of the effect, the multipath channel H S caused by the environmental scatterers is estimated, and the environmental information in the channel is separated. 3.根据权利要求1或2所述的遮挡效应下的毫米波一体化通信与感知方法,其特征在于,步骤2)中的毫米波信道模型和压缩感知重构问题为:3. the millimeter-wave integrated communication and sensing method under the occlusion effect according to claim 1 and 2, is characterized in that, the millimeter-wave channel model in step 2) and compressive sensing reconstruction problem are: 2.a)环境信息离散化,将整个空间内的环境信息视作为点云,点云中的每个点代表其周围大小为lr,wr和hr的小立方体的环境信息,这些小立方体称为像素;感知的散射体环境的长、宽、高分别为Lr、Wr和Hr,则空间内点云的数量为N=Lr/lr×Wr/wr×Hr/hr;每一个像素内部可能是空的,也可能是存在散射体;我们使用一个散射系数xn来表示第n个点云点的所在小立方体的散射系数,若小立方体内部是空的,则xn=0;所以整个空间的环境信息可以用如下变量表示:2.a) The environmental information is discretized, and the environmental information in the entire space is regarded as a point cloud. Each point in the point cloud represents the environmental information of small cubes with sizes l r , wr and hr around them. The cube is called a pixel; the length, width and height of the perceived scatterer environment are L r , W r and H r respectively, then the number of point clouds in the space is N=L r /l r ×W r /w r ×H r /h r ; each pixel may be empty or there may be scatterers; we use a scattering coefficient x n to represent the scattering coefficient of the small cube where the nth point cloud point is located, if the interior of the small cube is empty , then x n =0; so the environmental information of the entire space can be represented by the following variables: x=[x1,x2,…,xN]Tx=[x 1 , x 2 , . . . , x N ] T ; 2.b)将系统模型表示如下,在空间内的多个用户共享时频资源,在任意频率资源块上,第nR个BS接收天线接收信号表示为:2.b) The system model is expressed as follows, multiple users in space share time-frequency resources, and on any frequency resource block, the signal received by the n Rth BS receiving antenna is expressed as:
Figure FDA0003284405140000021
Figure FDA0003284405140000021
其中
Figure FDA0003284405140000022
表示NR个BS天线的接收信号,
Figure FDA0003284405140000023
表示Nu个用户发送长度为L码元的导频,w为噪声,
Figure FDA0003284405140000024
表示用户到接收天线直视路径的信道系数,
Figure FDA0003284405140000025
表示用户到空间点云位置的信道系数,
Figure FDA0003284405140000026
表示用户到空间点云位置的信道系数,且
Figure FDA0003284405140000027
其中
Figure FDA0003284405140000028
表示用户至空间点云位置的遮挡矩阵,当
Figure FDA0003284405140000029
出现全零列时,表示所有用户都无法感知到该像素点,即该像素点在感知范围外,
Figure FDA00032844051400000210
表示空间点云位置直射至接收天线的遮挡矩阵,当
Figure FDA00032844051400000211
出现零元素时,表示接收天线无法感知到该像素点,该像素点也在感知范围外,且
Figure FDA00032844051400000212
in
Figure FDA0003284405140000022
represents the received signals of NR BS antennas,
Figure FDA0003284405140000023
Indicates that Nu users send pilots of length L symbols, w is noise,
Figure FDA0003284405140000024
represents the channel coefficient of the direct line of sight from the user to the receiving antenna,
Figure FDA0003284405140000025
represents the channel coefficient from the user to the position of the spatial point cloud,
Figure FDA0003284405140000026
represents the channel coefficient from the user to the spatial point cloud location, and
Figure FDA0003284405140000027
in
Figure FDA0003284405140000028
represents the occlusion matrix from the user to the position of the spatial point cloud, when
Figure FDA0003284405140000029
When an all-zero column appears, it means that all users cannot perceive the pixel, that is, the pixel is outside the perception range.
Figure FDA00032844051400000210
represents the occlusion matrix that the position of the spatial point cloud directly hits the receiving antenna, when
Figure FDA00032844051400000211
When a zero element appears, it means that the receiving antenna cannot perceive the pixel, the pixel is also outside the perception range, and
Figure FDA00032844051400000212
2.c)基于上述系统模型,求解环境信息的优化问题表达为,2.c) Based on the above system model, the optimization problem of solving environmental information is expressed as,
Figure FDA00032844051400000213
Figure FDA00032844051400000213
在通信过程中,信道不确定性主要来自于未知的环境散射体分布,对于直视信道HLOS采用已知的统计模型描述,由此用户发送已知的导频s对信道HS进行估计,从而将求解问题转换为压缩感知重构问题:In the communication process, the channel uncertainty mainly comes from the unknown distribution of environmental scatterers. For the direct-view channel H LOS , a known statistical model is used to describe it, and the user sends the known pilot frequency s to estimate the channel H S. This transforms the solving problem into a compressive sensing reconstruction problem:
Figure FDA0003284405140000031
Figure FDA0003284405140000031
其将NR个BS接收天线所得的信道估计值拼接为矩阵,其中HS已知,
Figure FDA0003284405140000032
H已知,需要在遮挡矩阵V未知的条件下求解环境信息x。
It splices the channel estimates obtained by the N R BS receiving antennas into a matrix, where H S is known,
Figure FDA0003284405140000032
H is known, and the environment information x needs to be solved under the condition that the occlusion matrix V is unknown.
4.根据权利要求1或2所述的遮挡效应下的毫米波一体化通信与感知方法,其特征在于,步骤3)中的环境散射体分布与遮挡关系的模型为:4. The millimeter-wave integrated communication and perception method under occlusion effect according to claim 1 and 2, is characterized in that, the model of environmental scatterer distribution and occlusion relation in step 3) is: 像素点A和B之间存在一个存在散射系数的像素点C,我们提出了一种环境散射体分布与遮挡关系的模型进行遮挡判断,令像素点A位置为坐标原点,则像素点B和C的坐标分别表示为b和c,通过满足以下三点条件判断像素点C对像素点A和B之间的直视路径产生了遮挡;There is a pixel point C with a scattering coefficient between pixel points A and B. We propose a model of the relationship between the distribution of environmental scatterers and occlusion for occlusion judgment. Let the position of pixel point A be the coordinate origin, then pixel points B and C The coordinates of , are expressed as b and c, respectively, and it is judged that pixel C occludes the direct line of sight between pixels A and B by satisfying the following three conditions; 像素点C与像素点A和B之间的直视路径的距离d小于阈值β:The distance d of the direct path between pixel C and pixels A and B is less than the threshold β:
Figure FDA0003284405140000033
Figure FDA0003284405140000033
向量b与向量c之间的夹角为锐角:The angle between the vector b and the vector c is an acute angle: b·c>0;b·c>0; 像素点C位于像素点A和B之间:Pixel C is located between pixels A and B: ||c·b||<||b||2||c·b||<||b|| 2 .
5.根据权利要求3所述的一种毫米波环境感知方法,其特征在于,步骤4)中的一种基于因子图的概率推理模型为:5. a kind of millimeter wave environment perception method according to claim 3, is characterized in that, a kind of probability inference model based on factor graph in step 4) is:
Figure FDA0003284405140000034
与x后验概率分解为:
Will
Figure FDA0003284405140000034
The posterior probability with x is decomposed as:
Figure FDA0003284405140000041
Figure FDA0003284405140000041
根据步骤2)的压缩感知重构问题和步骤3)中的遮挡关系构建因子图模型;因子图包含变量节点
Figure FDA0003284405140000042
和x,函数节点
Figure FDA0003284405140000043
Figure FDA0003284405140000044
并令HS长度M=NuNR,其中:
A factor graph model is constructed according to the compressed sensing reconstruction problem in step 2) and the occlusion relationship in step 3); the factor graph contains variable nodes
Figure FDA0003284405140000042
and x, the function node
Figure FDA0003284405140000043
and
Figure FDA0003284405140000044
And let H S length M=N u N R , where:
Figure FDA0003284405140000045
Figure FDA0003284405140000045
其中,HS是z的加噪观测,噪声是方差为σw的高斯白噪声;px(x)表示x的先验分布为高斯-伯努利分布:where H S is the noised observation of z, and the noise is Gaussian white noise with variance σ w ; p x (x) indicates that the prior distribution of x is a Gauss-Bernoulli distribution: px(x)=(1-λ)δ(x)+λN(x;θx,σx);p x (x)=(1-λ)δ(x)+λN(x; θ x , σ x ); 其中λ表示环境信息x的稀疏度,
Figure FDA0003284405140000046
表示环境信息x对信道的遮挡关系,具体如下:
where λ represents the sparsity of the environmental information x,
Figure FDA0003284405140000046
Represents the occlusion relationship of the environmental information x to the channel, as follows:
Figure FDA0003284405140000047
Figure FDA0003284405140000047
6.根据权利要求5所述的遮挡效应下的毫米波一体化通信与感知方法,其特征在于,步骤5)中的利用遮挡关系的双线性近似消息传递方法为:6. millimeter wave integrated communication and perception method under occlusion effect according to claim 5, is characterized in that, the bilinear approximate message passing method utilizing occlusion relation in step 5) is: 5.a)初始化算法参数,令t为迭代次数,对于m=1,2,...,M,残差均值为
Figure FDA0003284405140000048
对于m=1,2,...,M,n=1,2,...,N,根据步骤4)中设定的px(x)和
Figure FDA0003284405140000049
生成x的估计的均值
Figure FDA00032844051400000410
和方差
Figure FDA00032844051400000411
Figure FDA00032844051400000412
的估计的均值
Figure FDA00032844051400000413
和方差
Figure FDA00032844051400000414
5.a) Initialize the algorithm parameters, let t be the number of iterations, and for m = 1, 2, ..., M, the residual mean is
Figure FDA0003284405140000048
For m = 1, 2, ..., M, n = 1, 2, ..., N, according to the p x (x) and
Figure FDA0003284405140000049
generate the estimated mean of x
Figure FDA00032844051400000410
and variance
Figure FDA00032844051400000411
Figure FDA00032844051400000412
the estimated mean of
Figure FDA00032844051400000413
and variance
Figure FDA00032844051400000414
5.b)对于m=1,2,...,M,计算
Figure FDA00032844051400000415
的观测均值
Figure FDA00032844051400000416
和方差
Figure FDA00032844051400000417
具体如下:
5.b) For m=1,2,...,M, calculate
Figure FDA00032844051400000415
The observed mean of
Figure FDA00032844051400000416
and variance
Figure FDA00032844051400000417
details as follows:
Figure FDA0003284405140000051
Figure FDA0003284405140000051
Figure FDA0003284405140000052
Figure FDA0003284405140000052
其中,in,
Figure FDA0003284405140000053
Figure FDA0003284405140000053
Figure FDA0003284405140000054
Figure FDA0003284405140000054
5.c)对于m=1,2,...,M,计算zm的估计的均值
Figure FDA0003284405140000055
和方差
Figure FDA0003284405140000056
具体如下:
5.c) For m = 1, 2, ..., M, compute the mean of the estimates of z m
Figure FDA0003284405140000055
and variance
Figure FDA0003284405140000056
details as follows:
Figure FDA0003284405140000057
Figure FDA0003284405140000057
Figure FDA0003284405140000058
Figure FDA0003284405140000058
其中,in,
Figure FDA0003284405140000059
Figure FDA0003284405140000059
其中,C是归一化变量;where C is the normalization variable; 5.d)对于m=1,2,...,M,计算残差的均值
Figure FDA00032844051400000510
和方差
Figure FDA00032844051400000511
具体如下:
5.d) For m = 1, 2, ..., M, calculate the mean of the residuals
Figure FDA00032844051400000510
and variance
Figure FDA00032844051400000511
details as follows:
Figure FDA00032844051400000512
Figure FDA00032844051400000512
Figure FDA00032844051400000513
Figure FDA00032844051400000513
5.e)对于m=1,2,...,M,n=1,2,...,N,计算
Figure FDA00032844051400000514
的观测均值
Figure FDA00032844051400000515
和方差
Figure FDA00032844051400000516
对于n=1,2,...,N,计算
Figure FDA00032844051400000517
的观测均值
Figure FDA00032844051400000518
和方差
Figure FDA00032844051400000519
具体如下:
5.e) For m=1,2,...,M,n=1,2,...,N, compute
Figure FDA00032844051400000514
The observed mean of
Figure FDA00032844051400000515
and variance
Figure FDA00032844051400000516
For n=1,2,...,N, compute
Figure FDA00032844051400000517
The observed mean of
Figure FDA00032844051400000518
and variance
Figure FDA00032844051400000519
details as follows:
Figure FDA00032844051400000520
Figure FDA00032844051400000520
Figure FDA0003284405140000061
Figure FDA0003284405140000061
Figure FDA0003284405140000062
Figure FDA0003284405140000062
Figure FDA0003284405140000063
Figure FDA0003284405140000063
5.f)对于m=1,2,...,M,n=1,2,...,N,计算
Figure FDA0003284405140000064
的估计均值
Figure FDA0003284405140000065
和方差
Figure FDA0003284405140000066
对于n=1,2,...,N,计算xn的观测均值
Figure FDA0003284405140000067
和方差
Figure FDA0003284405140000068
具体如下:
5.f) For m=1,2,...,M,n=1,2,...,N, calculate
Figure FDA0003284405140000064
The estimated mean of
Figure FDA0003284405140000065
and variance
Figure FDA0003284405140000066
For n = 1, 2, ..., N, compute the observed mean of x n
Figure FDA0003284405140000067
and variance
Figure FDA0003284405140000068
details as follows:
Figure FDA0003284405140000069
Figure FDA0003284405140000069
Figure FDA00032844051400000610
Figure FDA00032844051400000610
Figure FDA00032844051400000611
Figure FDA00032844051400000611
Figure FDA00032844051400000612
Figure FDA00032844051400000612
其中,in,
Figure FDA00032844051400000613
Figure FDA00032844051400000613
Figure FDA00032844051400000614
Figure FDA00032844051400000614
5.g)重复执行步骤b)至步骤f)直到达到收敛条件
Figure FDA00032844051400000615
得到环境信息x的估计值
Figure FDA00032844051400000616
5.g) Repeat step b) to step f) until the convergence condition is reached
Figure FDA00032844051400000615
Get an estimate of the environmental information x
Figure FDA00032844051400000616
CN202111156208.9A 2021-09-28 2021-09-28 Millimeter wave integrated communication and sensing method under shielding effect Active CN113965881B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111156208.9A CN113965881B (en) 2021-09-28 2021-09-28 Millimeter wave integrated communication and sensing method under shielding effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111156208.9A CN113965881B (en) 2021-09-28 2021-09-28 Millimeter wave integrated communication and sensing method under shielding effect

Publications (2)

Publication Number Publication Date
CN113965881A true CN113965881A (en) 2022-01-21
CN113965881B CN113965881B (en) 2022-11-01

Family

ID=79462827

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111156208.9A Active CN113965881B (en) 2021-09-28 2021-09-28 Millimeter wave integrated communication and sensing method under shielding effect

Country Status (1)

Country Link
CN (1) CN113965881B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115243311A (en) * 2022-07-29 2022-10-25 浙江大学 Iterative focusing type millimeter wave integrated communication and sensing method
WO2024032009A1 (en) * 2022-08-10 2024-02-15 浙江大学 Model evolution-based environment sensing method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2020647A1 (en) * 2007-08-03 2009-02-04 Insigna Security Srl Automatic multi-user system for localization, alarm and personal emergency, operating in multi-standard mode in aerial environment
US20140140375A1 (en) * 2012-11-19 2014-05-22 King Fahd University Of Petroleum And Minerals Method for compressive sensing , reconstruction, and estimation of ultra-wideband channels
CN106772365A (en) * 2016-11-25 2017-05-31 南京理工大学 A kind of multipath based on Bayes's compressed sensing utilizes through-wall radar imaging method
CN109738898A (en) * 2019-01-15 2019-05-10 西安电子科技大学 Transmitter, collector, receiver and communication perception system for trackside environment perception
WO2021027305A1 (en) * 2019-08-12 2021-02-18 华为技术有限公司 Method for determining perception information during communication transmission and related device
CN112769461A (en) * 2020-12-11 2021-05-07 华南理工大学 Large-scale antenna channel estimation method based on millimeter wave intelligent reflector communication

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2020647A1 (en) * 2007-08-03 2009-02-04 Insigna Security Srl Automatic multi-user system for localization, alarm and personal emergency, operating in multi-standard mode in aerial environment
US20140140375A1 (en) * 2012-11-19 2014-05-22 King Fahd University Of Petroleum And Minerals Method for compressive sensing , reconstruction, and estimation of ultra-wideband channels
CN106772365A (en) * 2016-11-25 2017-05-31 南京理工大学 A kind of multipath based on Bayes's compressed sensing utilizes through-wall radar imaging method
CN109738898A (en) * 2019-01-15 2019-05-10 西安电子科技大学 Transmitter, collector, receiver and communication perception system for trackside environment perception
WO2021027305A1 (en) * 2019-08-12 2021-02-18 华为技术有限公司 Method for determining perception information during communication transmission and related device
CN112769461A (en) * 2020-12-11 2021-05-07 华南理工大学 Large-scale antenna channel estimation method based on millimeter wave intelligent reflector communication

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIN TONG ET.: "Joint Multi-User Communication and Sensing", 《IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115243311A (en) * 2022-07-29 2022-10-25 浙江大学 Iterative focusing type millimeter wave integrated communication and sensing method
WO2024021440A1 (en) * 2022-07-29 2024-02-01 浙江大学 Iterative focused millimeter-wave integrated communication and sensing method
US12019181B2 (en) 2022-07-29 2024-06-25 Zhejiang University Iterative focused millimeter wave integrated communication and sensing method
CN115243311B (en) * 2022-07-29 2025-01-14 浙江大学 Iterative focusing millimeter wave integrated communication and sensing method
WO2024032009A1 (en) * 2022-08-10 2024-02-15 浙江大学 Model evolution-based environment sensing method

Also Published As

Publication number Publication date
CN113965881B (en) 2022-11-01

Similar Documents

Publication Publication Date Title
Yan et al. Passive beamforming and information transfer design for reconfigurable intelligent surfaces aided multiuser MIMO systems
CN110463147B (en) Method for decoding symbols and receiver for receiving and decoding symbols
US11467248B2 (en) Position estimation device and communication device
Yuan et al. Joint 3D beamforming and trajectory design for UAV-enabled mobile relaying system
CN106909779B (en) Calculation method of Cramero bound for MIMO radar based on distributed processing
CN115243311B (en) Iterative focusing millimeter wave integrated communication and sensing method
CN113965881B (en) Millimeter wave integrated communication and sensing method under shielding effect
Chaccour et al. Joint sensing and communication for situational awareness in wireless THz systems
EP4266603A1 (en) Space division-based data processing method and communication device
Jiang et al. Novel 3-D irregular-shaped geometry-based channel modeling for semi-ellipsoid vehicle-to-vehicle scattering environments
WO2024110362A1 (en) Method of joint communication and environment perception
Roy et al. Going beyond RF: A survey on how AI-enabled multimodal beamforming will shape the NextG standard
Zhou et al. Path loss model based on cluster at 28 GHz in the indoor and outdoor environments
Shaddad et al. Channel estimation for intelligent reflecting surface in 6G wireless network via deep learning technique
Colpaert et al. 3D non-stationary channel measurement and analysis for MaMIMO-UAV communications
CN114900904B (en) Unlicensed random access method in LEO satellite communication scene
Liu et al. Networked sensing in 6G cellular networks: Opportunities and challenges
CN114690175A (en) Target direct detection and tracking method based on passive external radiation source radar
Ashraf et al. Maximizing detection of target with multiple direction possibilities to support immersive communications in metaverse
Alexandropoulos et al. Joint localization and channel estimation for UAV-assisted millimeter wave communications
Kim et al. CNN-based Doppler shift estimation for low earth orbit satellites
Xu et al. Deep learning based channel covariance matrix estimation with scene images
Tong et al. Focused sensing in a wireless communication system
Oliveira et al. Machine learning-based mmwave mimo beam tracking in v2i scenarios: Algorithms and datasets
Li et al. Exploiting cell-free massive mimo for computational imaging

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