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
environment
sensing
millimeter wave
channel
occlusion
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 sensing method under shielding effect
Technical Field
The invention relates to the field of wireless communication, in particular to a perception and communication integrated system design in the field of new generation wireless communication.
Background
Millimeter wave (mmWave) has become a research focus in the current wireless communication field due to its characteristics of high bandwidth, high reliability and high integration, combined with large-scale multiple-input multiple-output (MIMO) technology. With the rapid development of the wireless communication industry, in consideration of the rapid increase of connection devices and services, the application scenarios of increasingly complex wireless mobile communication are considered, and the propagation environment of wireless communication signals is increasingly complex. Meanwhile, as wireless communication base stations are more densely deployed, and a receiver and a sender have higher computing power and physical performance than before, the information processing capability is greatly improved, which puts higher requirements on a communication system. Not only do they want to continue to provide the original communication services, but they also want the communication devices to use their own powerful processing power to accomplish the perception of the environment. Specifically, in future wireless communication scenarios, not only wireless broadband connection but also accurate environmental information including, but not limited to, the position, shape, state, electromagnetic characteristics, and the like of background scatterers of stationary or moving objects in the environment are required for new technologies such as smart cities, autonomous driving, and drone positioning.
How to use wireless communication equipment to sense environment based on a wireless communication architecture and realize sensing and communication integration is an important research direction of a next generation wireless communication system. When the uplink is used for environment sensing, the base station receives an uplink communication signal sent by a user, performs data communication and realizes passive environment sensing. One of the challenges in the design of a perceptual-communicative integrated system is the potentially large number of unknown variables in the environment, and thus the sparsity of the target environment itself should be exploited. For example, in cellular communication networks, buildings are sparsely distributed over the wireless network coverage. In addition, the scatterers generally have a shielding effect, and the scatterers close to the user can shield the electromagnetic signals, so that the electromagnetic signals cannot reach the scatterers far away from each other in the same propagation direction. Therefore, not all scatterers in the sensing range will affect the multipath channel of the same user, and the mutual shielding causes different users to face different scatterer environments, and multiple users need to perform multi-view joint sensing. At present, the existing environment perception algorithm does not utilize the relation between the distribution of environment scatterers and the shielding effect, the problem is not considered in a targeted manner, and the performance is poor when an imaging model with the shielding effect is solved.
In conclusion, the problem of inconsistency of the user observation target under the millimeter wave environment perception problem and the occlusion effect is comprehensively considered, and how to jointly realize the separation of the environment information in the uplink data and the resolution of the occlusion effect has higher research difficulty and practical significance.
Disclosure of Invention
The invention aims to solve the problem of how a base station utilizes uplink data transmitted by multiple users to sense the environment in an uplink wireless communication scene. The invention uses the pilot signal or other known data sequence of the existing communication system to sense, separate the receiving and transmitting processing, can be compatible with the existing communication system, and realizes the integration of sensing and communication. Considering that different users face different scatterer environments due to the shielding effect existing between environment scatterers, a millimeter wave integrated communication and sensing method under the shielding effect is provided.
The invention adopts the following specific technical scheme:
a millimeter wave integrated communication and sensing method under a shielding effect comprises the following steps:
1) in the T-th time slot, the base station receives pilot sequence s signals with the length of L, which are sent by all active users in the space;
2) after receiving the signal, the base station converts the environment sensing problem into a compressed sensing reconstruction problem based on a millimeter wave channel model considering the shielding effect;
3) establishing a model of the distribution and shielding relation of the environment scatterers based on the relative position relation between the environment scatterers;
4) calculating the relationship among the received data, the distribution of the environmental scatterers and the shielding effect, and establishing a probabilistic reasoning model based on a factor graph;
5) and (3) solving the compressed sensing reconstruction problem in the step (3) based on a bilinear approximate message transmission method by using an occlusion relation by combining the models obtained in the step (3) and the step (4), thereby realizing the sensing of the environment.
In one or more embodiments, in order to convert the environmental perception problem in step 2) into a compressed perception reconstruction problem, the multipath channel H caused by the environmental scatterers needs to be subjected to a millimeter wave channel model based on consideration of the occlusion effectSEstimation is performed and the environmental information in the channel is separated.
In one or more embodiments, the millimeter wave channel model and compressed sensing reconstruction problem in step 2) is:
a) discretizing the environmental information, regarding the environmental information in the whole space as a point cloud, each point in the point cloud representing its surrounding size lr,wrAnd hrEnvironmental information of the microcubes, which are called pixels; the length, width and height of the perceived scatterer environment are L respectivelyr、WrAnd HrIf the number of the point clouds in the space is N ═ Lr/lr×Wr/wr×Hr/hr(ii) a Each pixel may be empty or have scatterers present; we use a scattering coefficient xnTo represent the scattering coefficient of the microcube in which the cloud point of the nth point is located, if the interior of the microcube is empty, x isn0; the environmental information of the whole space can be represented by the following variables:
x=[x1,x2,…,xN]T
2, b) representing the system model as follows, a plurality of users in the space share time-frequency resources, and the nth frequency resource blockRA BS receivesThe antenna receives signals as
Figure BDA0003284405150000031
Wherein
Figure BDA0003284405150000032
Represents NRThe received signal of the antenna of the individual BS,
Figure BDA0003284405150000033
represents NuEach user transmits a pilot of length L symbols, w is noise,
Figure BDA0003284405150000034
representing the channel coefficients of the direct-view path of the user to the receive antenna,
Figure BDA0003284405150000035
channel coefficients representing the user to spatial point cloud locations,
Figure BDA0003284405150000036
representing the channel coefficients of the user to the spatial point cloud locations, an
Figure BDA0003284405150000037
Wherein
Figure BDA0003284405150000038
Occlusion matrix representing the user to the spatial point cloud location, when Vs1When all zero columns appear, it indicates that all users cannot perceive the pixel point, i.e. the pixel point is out of the perception range, Vs2(nR)∈{0,1}N×1A shielding matrix for representing the direct projection of the spatial point cloud position to the receiving antenna when Vs2(nR) When zero element appears, it represents that the receiving antenna can not sense the pixel point, the pixel point is out of the sensing range, and V (n)R)=Vs1(nR)diag(Vs2(nR));
2, c) based on the system model, the optimization problem for solving the environmental information is expressed as,
Figure BDA0003284405150000041
during communication, channel uncertainty comes mainly from unknown environmental scatterer distributions, for the direct-view channel HLOSUsing a known statistical model description whereby the user transmits a known pilot s for channel HSThe estimation is performed so as to transform the solution problem into a compressed sensing reconstruction problem:
Figure BDA0003284405150000042
which will be NRThe channel estimates from the individual BS receive antennas are concatenated into a matrix, where HSIn the known manner, it is known that,
Figure BDA0003284405150000043
h is known, the environment information x needs to be solved under the condition that the occlusion matrix V is unknown.
In one or more embodiments, the model of the environment scatterer distribution and occlusion relationship in step 3) is:
a pixel point C with a scattering coefficient exists between pixel points A and B, a model of environment scatterer distribution and shielding relation is provided for shielding judgment, the position of the pixel point A is made to be a coordinate origin, coordinates of the pixel points B and C are respectively expressed as B and C, and the condition that the pixel point C shields a direct-view path between the pixel points A and B is met;
the distance d between the pixel point C and the direct-view path between the pixel points a and B is less than the threshold value β:
Figure BDA0003284405150000044
the angle between vector b and vector c is acute:
b·c>0;
the pixel point C is positioned between the pixel points A and B:
||c·b||<||b||2
in one or more embodiments, one factor graph-based probabilistic inference model in step 4) is:
will be provided with
Figure BDA0003284405150000051
The posterior probability with x is decomposed into:
Figure BDA0003284405150000052
constructing a factor graph model according to the compressed sensing reconstruction problem in the step 2) and the shielding relation in the step 3); the factor graph comprises variable nodes
Figure BDA0003284405150000053
And x, function node
Figure BDA0003284405150000054
px(x) And
Figure BDA0003284405150000055
and make HSLength M ═ NuNRWherein:
Figure BDA0003284405150000056
wherein HSIs a noisy observation of z, the noise is a variance of σwWhite gaussian noise of (1); p is a radical ofx(x) The prior distribution representing x is a gaussian-bernoulli distribution:
px(x)=(1-λ)δ(x)+λN(x;θx,σx);
where x denotes the sparsity of the environmental information x,
Figure BDA0003284405150000057
representing environment information x-pairThe blocking relationship of the road is as follows:
Figure BDA0003284405150000058
in one or more embodiments, the bilinear approximate message passing method using the occlusion relation in the step 5) is as follows:
a) initializing algorithm parameters, let t be the number of iterations, and for M-1, 2
Figure BDA0003284405150000059
For M1, 2,.. ang.m, N1, 2,. ang.n, p set in step 4) is usedx(x) And
Figure BDA0003284405150000061
generating an estimated mean of x
Figure BDA0003284405150000062
Sum variance
Figure BDA0003284405150000063
Is estimated as a mean value of
Figure BDA0003284405150000064
Sum variance
Figure BDA0003284405150000065
B) calculating for M ═ 1, 2
Figure BDA0003284405150000066
Is observed as a mean value
Figure BDA0003284405150000067
Sum variance
Figure BDA0003284405150000068
The method comprises the following specific steps:
Figure BDA0003284405150000069
Figure BDA00032844051500000610
wherein,
Figure BDA00032844051500000611
Figure BDA00032844051500000612
c) calculating z for M ═ 1, 2mIs estimated as a mean value of
Figure BDA00032844051500000613
Sum variance
Figure BDA00032844051500000614
The method comprises the following specific steps:
Figure BDA00032844051500000615
Figure BDA00032844051500000616
wherein,
Figure BDA00032844051500000617
wherein C is a normalization variable;
d) calculating the mean value of the residuals for M1, 2
Figure BDA00032844051500000618
Sum variance
Figure BDA00032844051500000619
The method comprises the following specific steps:
Figure BDA00032844051500000620
Figure BDA00032844051500000621
e) for M1, 2,.., M, N1, 2.., N, calculating
Figure BDA0003284405150000071
Is observed as a mean value
Figure BDA0003284405150000072
Sum variance
Figure BDA0003284405150000073
For N1, 2
Figure BDA0003284405150000074
Is observed as a mean value
Figure BDA0003284405150000075
Sum variance
Figure BDA0003284405150000076
The method comprises the following specific steps:
Figure BDA0003284405150000077
Figure BDA0003284405150000078
Figure BDA0003284405150000079
Figure BDA00032844051500000710
f) calculating for M1, 2
Figure BDA00032844051500000711
Is estimated mean value of
Figure BDA00032844051500000712
Sum variance
Figure BDA00032844051500000723
For N1, 2nIs observed as a mean value
Figure BDA00032844051500000713
Sum variance
Figure BDA00032844051500000714
The method comprises the following specific steps:
Figure BDA00032844051500000715
Figure BDA00032844051500000716
Figure BDA00032844051500000717
Figure BDA00032844051500000718
wherein,
Figure BDA00032844051500000719
Figure BDA00032844051500000720
5.g) repeating steps b) to f) until a convergence condition is reached
Figure BDA00032844051500000721
Obtaining an estimate of the environmental information x
Figure BDA00032844051500000722
The invention has the beneficial effects that: in a wireless communication uplink, in a scene where a user sends data for environment sensing (for example, one multi-antenna base station uses uplink data sent by a plurality of single-antenna users for environment sensing), the millimeter wave environment sensing method under the shielding effect, namely the millimeter wave integrated communication and sensing method under the shielding effect, uses pilot signals or other known data sequences of the existing communication system for sensing, can be compatible with the existing communication system, and realizes the integration of sensing and communication. Firstly, a millimeter wave channel model considering the occlusion effect is constructed, and the environment sensing problem is converted into a compressed sensing reconstruction problem. And then establishing a model of the distribution and shielding relation of the environment scatterers based on the relative position relation between the environment scatterers. 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 an occlusion 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.
Drawings
FIG. 1 is a diagram of a millimeter wave environment perception scene under an occlusion effect;
FIG. 2 is a schematic diagram of the relationship between the distribution of environmental scatterers and occlusion;
FIG. 3 is a factor graph based on occlusion relationships;
FIG. 4 is a graph of the environmental perception visual results comparing the present invention with other compressed perception reconstruction algorithms;
FIG. 5 is a graph of the number of users versus the accuracy of the environmental perception MSE when comparing the present invention with other compressed sensing reconstruction algorithms;
FIG. 6 is a graph of signal-to-noise ratio (SNR) versus environmental perceptual accuracy (MSE) when comparing the present invention to other compressed perceptual reconstruction algorithms.
Detailed Description
As shown in fig. 1, first we consider a scenario where one Base Station (BS) is deployed in an outdoor area and there are multiple active Users (UEs). In an uplink communication scenario, multiple single-antenna users simultaneously transmit uplink communication signals to the BS. The transmitted signal is scattered by scatterers and multipath propagates to the BS for reception. As shown in the above diagram, the transmission signal of the user 1 is received by the AP after being scattered only by the target scatterers 1 and 3, and the target scatterer 2 does not affect this process.
An embodiment of the invention provides a millimeter wave integrated communication and sensing method under a shielding effect, which comprises the following steps:
1) in the T-th time slot, the base station receives pilot sequence s signals with the length of L, which are sent by all active users in the space;
2) after receiving the signal, the base station converts the environment sensing problem into a compressed sensing reconstruction problem based on a millimeter wave channel model considering the shielding effect;
in an embodiment of the present invention, in order to convert the environmental sensing problem in step 2) into a compressive sensing reconstruction problem, it is necessary to apply a millimeter wave channel model considering an occlusion effect to a multipath channel H caused by an environmental scattererSEstimation is performed and the environmental information in the channel is separated.
Specifically, in an embodiment of the present invention, the millimeter wave channel model and the compressed sensing reconstruction problem in step 2) are as follows:
a) discretizing the environment information, regarding the environment information in the whole space as a point cloud, each point in the point cloud representing its surroundings with a size of lr,wrAnd hrThe microcubes are called pixels. The length, width and height of the sensed scatterer environment are respectively Lr,WrAnd HrIf the number of the point clouds in the space is N ═ Lr/lr×Wr/wr×Hr/hr. Each pixel may be empty or have scatterers present. We use a scattering coefficient xnTo represent the scattering coefficient of the microcube in which the cloud point of the nth point is located, if the interior of the microcube is empty, x isn0. The environmental information of the whole space can be represented by the following variables:
x=[x1,x2,…,xN]T
2.b) the system model is represented as follows. Multiple users in space share time frequency resource, on arbitrary frequency resource block, the n-thRThe BS receive antenna receive signal is represented as:
Figure BDA0003284405150000091
wherein
Figure BDA0003284405150000101
Represents NRThe received signal of the antenna of the individual BS,
Figure BDA0003284405150000102
represents NuEach user transmits a pilot of length L symbols and w is noise.
Figure BDA0003284405150000103
Representing the channel coefficients of the direct-view path of the user to the receive antenna,
Figure BDA0003284405150000104
channel coefficients representing the user to spatial point cloud locations,
Figure BDA0003284405150000105
representing the channel coefficients of the user to the spatial point cloud locations, an
Figure BDA0003284405150000106
Wherein
Figure BDA0003284405150000107
Occlusion matrix representing the user to the spatial point cloud location, when Vs1When all zero columns appear, all users cannot perceive the pixel point, namely the pixel point is out of the perception range. Vs2(nR)∈{0,1}N×1A shielding matrix for representing the direct projection of the spatial point cloud position to the receiving antenna when Vs2(nR) When zero element appears, it represents that the receiving antenna can not sense the pixel point, the pixel point is out of the sensing range, and V (n)R)=Vs1(nR)diag(Vs2(nR))。
2, c) based on the system model, the optimization problem for solving the environmental information is expressed as,
Figure BDA0003284405150000108
during communication, channel uncertainty comes mainly from unknown environmental scatterer distributions, for the direct-view channel HLOSDescribed using a known statistical model. Whereby the user transmits a known pilot s on channel HSThe estimation is performed so as to transform the solution problem into a compressed perceptual reconstruction problem, namely:
Figure BDA0003284405150000109
which will be NRObtained by a BS receiving antennaIs concatenated into a matrix, where HSIn the known manner, it is known that,
Figure BDA00032844051500001010
h is known. The environment information x needs to be solved under the condition that the occlusion matrix V is unknown.
3) Establishing a model of the distribution and shielding relation of the environment scatterers based on the relative position relation between the environment scatterers;
specifically, in an embodiment of the present invention, the model of the environment scatterer distribution and the occlusion relationship in step 3) is:
as shown in fig. 2, a pixel point C having a scattering coefficient exists between pixel points a and B, and we propose a model of the relationship between the distribution of the environmental scatterers and the occlusion to perform occlusion determination. Let the position of pixel point a be the origin of coordinates, then the coordinates of pixel points B and C are denoted as B and C, respectively. The pixel point C is judged to shield the direct-view path between the pixel points A and B by meeting the following three-point condition.
The distance d between pixel C and the direct-view path between pixels a and B is less than the threshold beta,
Figure BDA0003284405150000111
the angle between vector b and vector c is acute,
b·c>0;
pixel C is located between pixels a and B,
||c·b||<||b||2
4) and calculating the relation among the received data, the distribution of the environmental scatterers and the shielding effect, and establishing a probabilistic reasoning model based on a factor graph.
Specifically, in an embodiment of the present invention, a probabilistic inference model based on a factor graph in step 4) is specifically as follows.
Will be provided with
Figure BDA0003284405150000112
And the x posterior probability is decomposed into,
Figure BDA0003284405150000113
as shown in fig. 3, a factor graph model is constructed according to the compressed sensing reconstruction problem of step 2) and the occlusion relation in step 3). The factor graph comprises variable nodes
Figure BDA0003284405150000114
And x, function node
Figure BDA0003284405150000115
px(x) And
Figure BDA0003284405150000116
for the sake of brevity, let HSLength M ═ NuNR
Wherein
Figure BDA0003284405150000117
Represents HSIs a noisy observation of z, the noise is a variance of σwWhite gaussian noise. p is a radical ofx(x) The prior distribution representing x is a gaussian-bernoulli distribution,
px(x)=(1-λ)δ(x)+λN(x;θx,σx);
where λ represents the sparsity of the environmental information x.
Figure BDA0003284405150000121
Representing the occlusion relationship of the environment information x to the channel,
Figure BDA0003284405150000122
5) and (3) solving the compressed sensing reconstruction problem in the step (3) based on a bilinear approximate message transmission method by using an occlusion relation by combining the models obtained in the step (3) and the step (4), thereby realizing the sensing of the environment.
Specifically, in an embodiment of the present invention, the bilinear approximate message passing method using the occlusion relationship in this step is:
a) initializing algorithm parameters. Let t be the number of iterations, for M1, 2
Figure BDA0003284405150000123
For M1, 2,.. ang.m, N1, 2,. ang.n, p set in step 4) is usedx(x) And
Figure BDA0003284405150000124
generating an estimated mean of x
Figure BDA0003284405150000125
Sum variance
Figure BDA0003284405150000126
Is estimated as a mean value of
Figure BDA0003284405150000127
Sum variance
Figure BDA0003284405150000128
B) calculating for M ═ 1, 2
Figure BDA0003284405150000129
Is observed as a mean value
Figure BDA00032844051500001210
Sum variance
Figure BDA00032844051500001211
Figure BDA00032844051500001212
Figure BDA00032844051500001213
Wherein,
Figure BDA00032844051500001214
Figure BDA00032844051500001215
c) calculating z for M ═ 1, 2mIs estimated as a mean value of
Figure BDA00032844051500001216
Sum variance
Figure BDA00032844051500001217
Figure BDA0003284405150000131
Figure BDA0003284405150000132
Wherein,
Figure BDA0003284405150000133
c is a normalization variable.
D) calculating the mean value of the residuals for M1, 2
Figure BDA0003284405150000134
Sum variance
Figure BDA0003284405150000135
Figure BDA0003284405150000136
Figure BDA0003284405150000137
E) for M1, 2,.., M, N1, 2.., N, calculating
Figure BDA0003284405150000138
Is observed as a mean value
Figure BDA0003284405150000139
Sum variance
Figure BDA00032844051500001310
For N1, 2
Figure BDA00032844051500001311
Is observed as a mean value
Figure BDA00032844051500001312
Sum variance
Figure BDA00032844051500001313
Figure BDA00032844051500001314
Figure BDA00032844051500001315
Figure BDA00032844051500001316
Figure BDA00032844051500001317
F) forM1, 2, 1, N, and calculating
Figure BDA00032844051500001318
Is estimated mean value of
Figure BDA00032844051500001319
Sum variance
Figure BDA00032844051500001320
For N1, 2nIs observed as a mean value
Figure BDA00032844051500001321
Sum variance
Figure BDA00032844051500001322
Figure BDA00032844051500001323
Figure BDA00032844051500001324
Figure BDA00032844051500001325
Figure BDA0003284405150000141
Wherein
Figure BDA0003284405150000142
Figure BDA0003284405150000143
G) repeating the stepStep b) to step f) until a convergence condition is reached
Figure BDA0003284405150000144
Obtaining an estimate of the environmental information x
Figure BDA0003284405150000145
As can be seen by computer simulation: as shown in fig. 4, where the gamma algorithm completely ignores occlusion effects. The Biliner GAMP algorithm considers the occlusion matrix V and the environment information x as independent variables and does not use the occlusion relationship to respectively solve. Compared with the former two algorithms, the millimeter wave environment sensing algorithm based on the shielding effect obviously improves the accuracy of environment sensing. Fig. 5 shows that the environmental perception effect of the method of the present invention is gradually improved and superior to the existing algorithm as the number of users increases. Fig. 5 shows that the environmental perception effect of the method of the present invention is gradually improved and superior to the existing algorithm as the signal-to-noise ratio increases.
In summary, in a wireless communication uplink, in a scenario where a user transmits data for environment sensing (for example, one multi-antenna base station performs environment sensing by using uplink data transmitted by multiple single-antenna users), the millimeter wave environment sensing method under the occlusion effect proposed in the embodiment of the present invention is a millimeter wave integrated communication and sensing method under the occlusion effect, which senses by using a pilot signal or other known data sequence of an existing communication system, and is compatible with the existing communication system, thereby implementing integration of sensing and communication. Firstly, a millimeter wave channel model considering the occlusion effect is constructed in the embodiment of the invention, and the environment sensing problem is converted into a compressed sensing reconstruction problem. And then establishing a model of the distribution and shielding relation of the environment scatterers based on the relative position relation between the environment scatterers. 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 an occlusion 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.
The above examples are intended to illustrate the invention, but not to limit it. Any modification and variation of the present invention within the spirit of the present invention and the scope of the claims will fall within the scope of the present invention.

Claims (6)

1. A millimeter wave integrated communication and sensing method under a shielding effect is characterized by comprising the following steps:
1) in the T-th time slot, the base station receives pilot sequence s signals with the length of L, which are sent by all active users in the space;
2) after receiving the signal, the base station converts the environment sensing problem into a compressed sensing reconstruction problem based on a millimeter wave channel model considering the shielding effect;
3) establishing a model of the distribution and shielding relation of the environment scatterers based on the relative position relation between the environment scatterers;
4) calculating the relationship among the received data, the distribution of the environmental scatterers and the shielding effect, and establishing a probabilistic reasoning model based on a factor graph;
5) and (3) solving the compressed sensing reconstruction problem in the step (3) based on a bilinear approximate message transmission method by using an occlusion relation by combining the models obtained in the step (3) and the step (4), thereby realizing the sensing of the environment.
2. The millimeter wave integrated communication and sensing method under the occlusion effect according to claim 1, wherein in order to convert the environment sensing problem in step 2) into the compressive sensing reconstruction problem, the multipath channel H caused by the environment scatterer needs to be subjected to the multi-path channel H based on a millimeter wave channel model considering the occlusion effectSEstimation is performed and the environmental information in the channel is separated.
3. The millimeter wave integrated communication and sensing method under the occlusion effect according to claim 1 or 2, wherein the millimeter wave channel model and the compressive sensing reconstruction problem in step 2) are as follows:
a) discretizing the environmental information, regarding the environmental information in the whole space as a point cloud, each point in the point cloud representing its surrounding size lr,wrAnd hrEnvironmental information of the microcubes, which are called pixels; the length, width and height of the perceived scatterer environment are L respectivelyr、WrAnd HrIf the number of the point clouds in the space is N ═ Lr/lr×Wr/wr×Hr/hr(ii) a Each pixel may be empty or have scatterers present; we use a scattering coefficient xnTo represent the scattering coefficient of the microcube in which the cloud point of the nth point is located, if the interior of the microcube is empty, x isn0; the environmental information of the whole space can be represented by the following variables:
x=[x1,x2,…,xN]T
2, b) representing the system model as follows, a plurality of users in the space share time-frequency resources, and the nth frequency resource blockRThe BS receive antenna receive signal is represented as:
Figure FDA0003284405140000021
wherein
Figure FDA0003284405140000022
Represents NRThe received signal of the antenna of the individual BS,
Figure FDA0003284405140000023
represents NuEach user transmits a pilot of length L symbols, w is noise,
Figure FDA0003284405140000024
representing the channel coefficients of the direct-view path of the user to the receive antenna,
Figure FDA0003284405140000025
channel coefficients representing the user to spatial point cloud locations,
Figure FDA0003284405140000026
representing the channel coefficients of the user to the spatial point cloud locations, an
Figure FDA0003284405140000027
Wherein
Figure FDA0003284405140000028
An occlusion matrix representing a user to a spatial point cloud location when
Figure FDA0003284405140000029
When all zero columns appear, the pixel point can not be sensed by all users, namely the pixel point is out of the sensing range,
Figure FDA00032844051400000210
a shielding matrix for representing the direct projection of the spatial point cloud position to the receiving antenna when
Figure FDA00032844051400000211
When zero element appears, it represents that the receiving antenna can not sense the pixel point, and the pixel point is out of the sensing range, and
Figure FDA00032844051400000212
2, c) based on the system model, the optimization problem for solving the environmental information is expressed as,
Figure FDA00032844051400000213
during communication, channel uncertainty comes mainly from unknown environmental scatterer distributions, for the direct-view channel HLOSUsing a known statistical model description whereby the user transmits a known pilot s for channel HSThe estimation is performed so as to transform the solution problem into a compressed sensing reconstruction problem:
Figure FDA0003284405140000031
which will be NRThe channel estimates from the individual BS receive antennas are concatenated into a matrix, where HSIn the known manner, it is known that,
Figure FDA0003284405140000032
h is known, the environment information x needs to be solved under the condition that the occlusion matrix V is unknown.
4. The millimeter wave integrated communication and perception method under the occlusion effect according to claim 1 or 2, wherein the model of the environment scatterer distribution and occlusion relationship in step 3) is:
a pixel point C with a scattering coefficient exists between pixel points A and B, a model of environment scatterer distribution and shielding relation is provided for shielding judgment, the position of the pixel point A is made to be a coordinate origin, coordinates of the pixel points B and C are respectively expressed as B and C, and the condition that the pixel point C shields a direct-view path between the pixel points A and B is met;
the distance d between the pixel point C and the direct-view path between the pixel points a and B is less than the threshold value β:
Figure FDA0003284405140000033
the angle between vector b and vector c is acute:
b·c>0;
the pixel point C is positioned between the pixel points A and B:
||c·b||<||b||2
5. the millimeter wave environment sensing method according to claim 3, wherein the probabilistic inference model based on the factor graph in step 4) is:
will be provided with
Figure FDA0003284405140000034
The posterior probability with x is decomposed into:
Figure FDA0003284405140000041
constructing a factor graph model according to the compressed sensing reconstruction problem in the step 2) and the shielding relation in the step 3); the factor graph comprises variable nodes
Figure FDA0003284405140000042
And x, function node
Figure FDA0003284405140000043
And
Figure FDA0003284405140000044
and make HSLength M ═ NuNRWherein:
Figure FDA0003284405140000045
wherein HSIs a noisy observation of z, the noise is a variance of σwWhite gaussian noise of (1); p is a radical ofx(x) The prior distribution representing x is a gaussian-bernoulli distribution:
px(x)=(1-λ)δ(x)+λN(x;θx,σx);
where x denotes the sparsity of the environmental information x,
Figure FDA0003284405140000046
the occlusion relation of the environment information x to the channel is represented as follows:
Figure FDA0003284405140000047
6. the millimeter wave integrated communication and perception method under the occlusion effect according to claim 5, wherein the bilinear approximate message passing method using the occlusion relationship in step 5) is:
a) initializing algorithm parameters, let t be the number of iterations, and for M-1, 2
Figure FDA0003284405140000048
For M1, 2,.. ang.m, N1, 2,. ang.n, p set in step 4) is usedx(x) And
Figure FDA0003284405140000049
generating an estimated mean of x
Figure FDA00032844051400000410
Sum variance
Figure FDA00032844051400000411
Figure FDA00032844051400000412
Is estimated as a mean value of
Figure FDA00032844051400000413
Sum variance
Figure FDA00032844051400000414
B) calculating for M ═ 1, 2
Figure FDA00032844051400000415
Is observed as a mean value
Figure FDA00032844051400000416
Sum variance
Figure FDA00032844051400000417
The method comprises the following specific steps:
Figure FDA0003284405140000051
Figure FDA0003284405140000052
wherein,
Figure FDA0003284405140000053
Figure FDA0003284405140000054
c) calculating z for M ═ 1, 2mIs estimated as a mean value of
Figure FDA0003284405140000055
Sum variance
Figure FDA0003284405140000056
The method comprises the following specific steps:
Figure FDA0003284405140000057
Figure FDA0003284405140000058
wherein,
Figure FDA0003284405140000059
wherein C is a normalization variable;
d) calculating the mean value of the residuals for M1, 2
Figure FDA00032844051400000510
Sum variance
Figure FDA00032844051400000511
The method comprises the following specific steps:
Figure FDA00032844051400000512
Figure FDA00032844051400000513
e) for M1, 2,.., M, N1, 2.., N, calculating
Figure FDA00032844051400000514
Is observed as a mean value
Figure FDA00032844051400000515
Sum variance
Figure FDA00032844051400000516
For N1, 2
Figure FDA00032844051400000517
Is observed as a mean value
Figure FDA00032844051400000518
Sum variance
Figure FDA00032844051400000519
The method comprises the following specific steps:
Figure FDA00032844051400000520
Figure FDA0003284405140000061
Figure FDA0003284405140000062
Figure FDA0003284405140000063
f) calculating for M1, 2
Figure FDA0003284405140000064
Is estimated mean value of
Figure FDA0003284405140000065
Sum variance
Figure FDA0003284405140000066
For N1, 2nIs observed as a mean value
Figure FDA0003284405140000067
Sum variance
Figure FDA0003284405140000068
The method comprises the following specific steps:
Figure FDA0003284405140000069
Figure FDA00032844051400000610
Figure FDA00032844051400000611
Figure FDA00032844051400000612
wherein,
Figure FDA00032844051400000613
Figure FDA00032844051400000614
5.g) repeating steps b) to f) until a convergence condition is reached
Figure FDA00032844051400000615
Obtaining 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
WO2024021440A1 (en) * 2022-07-29 2024-02-01 浙江大学 Iterative focused 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 the communication sensory perceptual system of trackside environment sensing
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 the communication sensory perceptual system of trackside environment sensing
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 (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
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
Yuan et al. Reconfigurable-intelligent-surface empowered wireless communications: Challenges and opportunities
Yan et al. Passive beamforming and information transfer design for reconfigurable intelligent surfaces aided multiuser MIMO systems
Wang et al. Intelligent reflecting surface-assisted millimeter wave communications: Joint active and passive precoding design
EP3406059B1 (en) Method for decoding symbol and receiver for receiving and decoding symbol
Imoize et al. Standard propagation channel models for MIMO communication systems
CN113965881B (en) Millimeter wave integrated communication and sensing method under shielding effect
CN102457319B (en) Business data distribution method and device
CN112039564B (en) Satellite large-scale MIMO wide-coverage pre-coding transmission method
Wang et al. Sparse channel estimation in millimeter wave communications: Exploiting joint AoD-AoA angular spread
WO2024021440A1 (en) Iterative focused millimeter-wave integrated communication and sensing method
Iimori et al. A stochastic gradient descent approach for hybrid mmWave beamforming with blockage and CSI-error robustness
Chaccour et al. Joint sensing and communication for situational awareness in wireless THz systems
WO2024110362A1 (en) Method of joint communication and environment perception
Ebrahiem et al. A deep learning approach for channel estimation in 5G wireless communications
CN116743222A (en) Multiplexing beam forming method, communication perception calculation integrated system and related device
Chen et al. Performance evaluation of BER for an Massive-MIMO with M-ary PSK scheme over Three-Dimension correlated channel
Ashraf et al. Maximizing detection of target with multiple direction possibilities to support immersive communications in Metaverse
Teng Outage probability and ergodic capacity analysis of satellite–terrestrial NOMA system with mixed RF/mmWave relaying
Li et al. A GAN-GRU Based Space-Time Predictive Channel Model for 6G Wireless Communications
Sande et al. Fast convergence for efficient beamforming in massive MIMO heterogeneous networks
Yamada et al. Plane-wave and vector-rotation approximation technique for reducing computational complexity to simulate MIMO propagation channel using ray-tracing
Al-Tous et al. Static reflecting surface based on population-level optimization
Tong et al. Focused sensing in a wireless communication system
US11438044B1 (en) Interference aware eigen-beamforming based on second order statistics
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