CN113286362A - Method for acquiring arrival time and arrival angle of multipath signal and related device - Google Patents
Method for acquiring arrival time and arrival angle of multipath signal and related device Download PDFInfo
- Publication number
- CN113286362A CN113286362A CN202110553482.3A CN202110553482A CN113286362A CN 113286362 A CN113286362 A CN 113286362A CN 202110553482 A CN202110553482 A CN 202110553482A CN 113286362 A CN113286362 A CN 113286362A
- Authority
- CN
- China
- Prior art keywords
- arrival
- matrix
- signals
- signal
- frequency response
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 239000011159 matrix material Substances 0.000 claims abstract description 159
- 230000004044 response Effects 0.000 claims abstract description 67
- 239000013598 vector Substances 0.000 claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 24
- 230000006870 function Effects 0.000 claims description 27
- 238000005314 correlation function Methods 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 13
- 238000005516 engineering process Methods 0.000 claims description 13
- 238000009499 grossing Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000000969 carrier Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 1
- GVVPGTZRZFNKDS-JXMROGBWSA-N geranyl diphosphate Chemical compound CC(C)=CCC\C(C)=C\CO[P@](O)(=O)OP(O)(O)=O GVVPGTZRZFNKDS-JXMROGBWSA-N 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
The utility model provides a method and a related device for acquiring the arrival time and the arrival angle of a multipath signal, which receive the signals which are sent by a base station and arrive through a plurality of paths through a plurality of antennas, and obtain the arrival time of the multipath signal according to the correlation between the sent signals and the received signals; the method comprises the steps of constructing a channel frequency response matrix according to signals of each subcarrier received by each antenna and signals of each subcarrier sent by a base station, constructing a noise matrix according to a smooth channel frequency response matrix, constructing a steering vector matrix according to phase differences of different subcarriers for the same antenna and phase differences of different antennas for the same subcarrier, constructing a 2D-MUSIC spectrum function by using the steering vector matrix and the noise matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function, thereby obtaining the arrival time and the arrival angle of matched multipath signals and providing accurate parameters for positioning wireless signals.
Description
Technical Field
The present disclosure relates to the field of wireless signal positioning technologies, and in particular, to a method and a related apparatus for obtaining an arrival time and an arrival angle of a multipath signal.
Background
With the development of the internet of things technology, the service based on the location awareness is greatly concerned by people, and how to acquire high-precision location information becomes a hotspot of research of people. The existing positioning methods are mostly based on sensor positioning, such as global navigation satellite system GNSS, lidar, inertial sensor IMU, and the like. Positioning technology based on satellite signals is relatively mature, but because satellite signals are easily shielded, the positioning mode can be more applied to open outdoor scenes, meanwhile, a laser radar is easily affected by light and bad weather, and an IMU (inertial measurement Unit) easily generates accumulated errors.
With the development of 5G mobile communication systems, the deployment of network infrastructures is more and more extensive, the low delay, high resolution, high reliability and the like of 5G networks make wireless signal positioning more attractive, wireless signal-based positioning can utilize the existing 5G network infrastructures to position networked vehicles, is not limited by city intensive environments and weather, and does not need to install additional sensor equipment in the vehicles. The parameter estimation is the basis of wireless signal positioning, and common parameters include signal arrival time, signal arrival angle and the like. However, in the related art, the time of arrival and the angle of arrival of the acquired signal are relatively single, and accurate parameters cannot be provided for wireless signal positioning.
Disclosure of Invention
In view of the above, the present disclosure is directed to a method and a related apparatus for acquiring time of arrival and angle of arrival of a multipath signal.
In view of the above, the present disclosure provides a method for acquiring time of arrival and angle of arrival of a multipath signal, including:
acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
constructing a correlation function of the transmitting signal and the receiving signal, and taking the time corresponding to a plurality of peak values which are higher than a peak value threshold value in the correlation function as a plurality of arrival times;
constructing a channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier sent by the base station, and constructing a noise matrix according to the channel frequency response matrix;
constructing a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and constructing a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
Based on the same inventive concept, the present disclosure provides a device for acquiring the arrival time and the arrival angle of a multipath signal, comprising:
the signal acquisition module is used for acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
an arrival time determining module, configured to construct a correlation function between the transmitted signal and the received signal, and use, as multiple arrival times, times corresponding to multiple peaks higher than a peak threshold in the correlation function;
a noise matrix constructing module, configured to construct a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station, and construct a noise matrix according to the channel frequency response matrix;
the steering vector matrix constructing module is used for constructing and obtaining a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and the arrival angle determining module is used for constructing and obtaining a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
Based on the same inventive concept, the present disclosure provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, which when executing the program implements the method as described above.
Based on the same inventive concept, the present disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to perform the above-described method.
As can be seen from the foregoing, the method and related apparatus for obtaining time of arrival and angle of arrival of multipath signals according to the present disclosure receive signals transmitted by a base station and arriving through multiple paths through multiple antennas, and obtain the time of arrival of multipath signals according to the correlation between the transmitted signals and the received signals; the method comprises the steps of constructing a channel frequency response matrix according to signals of each subcarrier received by each antenna and signals of each subcarrier sent by a base station, constructing a noise matrix according to a smooth channel frequency response matrix, constructing a steering vector matrix according to phase differences of different subcarriers for the same antenna and phase differences of different antennas for the same subcarrier, constructing a 2D-MUSIC spectrum function by using the steering vector matrix and the noise matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function, thereby obtaining the arrival time and the arrival angle of matched multipath signals and providing accurate parameters for positioning wireless signals.
Drawings
In order to more clearly illustrate the technical solutions in the present disclosure or related technologies, the drawings needed to be used in the description of the embodiments or related technologies are briefly introduced below, and it is obvious that the drawings in the following description are only embodiments of the present disclosure, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a schematic flowchart of a method for acquiring time of arrival and angle of arrival of a multipath signal according to an embodiment of the present disclosure;
fig. 2 is a schematic diagram of a scenario in which a transmission signal reaches multiple antennas through multiple different paths according to an embodiment of the present disclosure;
fig. 3 is a schematic structural diagram of a multipath signal time-of-arrival and angle-of-arrival obtaining apparatus according to an embodiment of the present disclosure;
fig. 4 is a more specific hardware structure diagram of an electronic device according to an embodiment of the present disclosure.
Detailed Description
For the purpose of promoting a better understanding of the objects, aspects and advantages of the present disclosure, reference is made to the following detailed description taken in conjunction with the accompanying drawings.
It is to be noted that technical terms or scientific terms used in the embodiments of the present disclosure should have a general meaning as understood by those having ordinary skill in the art to which the present disclosure belongs, unless otherwise defined. The use of "first," "second," and similar terms in the embodiments of the disclosure is not intended to indicate any order, quantity, or importance, but rather to distinguish one element from another. The word "comprising" or "comprises", and the like, means that the element or item listed before the word covers the element or item listed after the word and its equivalents, but does not exclude other elements or items. The terms "connected" or "coupled" and the like are not restricted to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "upper", "lower", "left", "right", and the like are used merely to indicate relative positional relationships, and when the absolute position of the object being described is changed, the relative positional relationships may also be changed accordingly.
With the development of the internet of things technology, the service based on the location awareness is greatly concerned by people, and how to acquire high-precision location information becomes a hotspot of research of people. The existing positioning methods are mostly based on sensor positioning, such as global navigation satellite system GNSS, lidar, inertial sensor IMU, and the like. Positioning technology based on satellite signals is relatively mature, but because satellite signals are easily shielded, the positioning mode can be more applied to open outdoor scenes, meanwhile, a laser radar is easily affected by light and bad weather, and an IMU (inertial measurement Unit) easily generates accumulated errors.
With the development of 5G mobile communication systems, the deployment of network infrastructures is more and more extensive, the low delay, high resolution, high reliability and the like of 5G networks make wireless signal positioning more attractive, wireless signal-based positioning can utilize the existing 5G network infrastructures to position networked vehicles, is not limited by city intensive environments and weather, and does not need to install additional sensor equipment in the vehicles. The parameter estimation is the basis for the positioning of the radio signal. Common parameters include signal time of arrival, signal angle of arrival, and the like. However, in the related art, the time of arrival and the angle of arrival of the acquired signal are relatively single, and accurate parameters cannot be provided for wireless positioning.
The inventor proposes that multipath propagation can be generated due to reflection of urban buildings, indoor walls and the like of wireless signals, multipath components comprise geometric information of the position of a target to be detected and the surrounding environment, positioning accuracy can be improved by using multipath signals, and therefore accurate estimation needs to be carried out on arrival time and arrival angle of the multipath signals, and accurate parameters can be provided for wireless signal positioning.
Referring to fig. 1, it is a schematic flowchart of a method for acquiring time of arrival and angle of arrival of a multipath signal according to an embodiment of the present disclosure. Multipath refers to a propagation phenomenon that a signal reaches a receiving antenna from a transmitting antenna through a plurality of paths, and multipath propagation is caused by scattering of the signal by the atmosphere, reflection and refraction of the signal by an ionosphere, and reflection of the signal by surface objects such as mountains and buildings. For a signal propagating through any of the multiple paths, the time of arrival refers to the time taken for the signal to reach the receiving antenna from the transmitting antenna, i.e., the propagation time; the angle of arrival refers to the angle of incidence of the signal when it reaches the receiving antenna, and the angle of incidence is the angle between the incident ray and the normal of the incident surface.
The method for acquiring the arrival time and the arrival angle of the multipath signal comprises the following steps:
s110, acquiring signals sent by a base station as sending signals, and acquiring signals received by a plurality of antennas as receiving signals.
In some embodiments, a transmitting antenna is arranged on the base station side, a receiving antenna is arranged on the target side to be measured, and the transmitting signal passes through a plurality of paths from the transmitting antenna to the receiving antenna. The present disclosure uses the side of the target to be measured as its own viewing angle, so the antenna is called a receiving antenna in the present disclosure.
The base station sends a 5G downlink signal to a target to be detected, the target to be detected utilizes a CSI-RS/PRS signal in the received downlink signal to perform multipath signal parameter estimation, and specifically, the arrival time and the arrival angle of the multipath signal are obtained. The base station uplink and downlink signals are transmitted by the base station by taking the base station as a reference, the signal received by the target to be detected is a downlink signal and is transmitted by the target to be detected, and the signal received by the base station is an uplink signal. The format and characteristics of the CSI-RS/PRS signal conform to the R16 standard of 3 GPP. The downlink signal transmitted by the base station is represented in the form of a time domain signal, and the transmitted signal is represented as s (t), wherein s (t) is a function which changes along with time, t refers to a time variable, and for the transmitted signal, the time variable is the transmission time.
In some embodiments, acquiring a signal sent by a base station as a transmission signal refers to directly acquiring data of the transmission signal from a base station side, rather than acquiring a signal received through an antenna from a target side to be measured, where the signal received through the antenna by the target side to be measured is a received signal.
In some embodiments, the plurality of antennas are arranged in a row, and the distance between every two adjacent antennas is the same and is greater than half of the wavelength of the transmission signal.
In some embodiments, the transmit signal travels multiple paths to multiple antennas; the received signal is a superposition of the transmitted signals that arrive at the multiple antennas via multiple paths.
Referring to fig. 2, where d is a distance between two adjacent antennas, and θ is an arrival angle of a signal of a certain path.
The transmission signals reach the multiple antennas through multiple paths, and the reception signals are the superposition of the transmission signals reaching the multiple antennas through the multiple paths, the reception signals can be expressed as:
where r (t) is the received signal, L is the number of paths, L is the first path, alFor the path gain of the l path, τlFor the path delay of the l-th path, n (t) is noise, t is time variable, and for the received signal, the time variable is the receiving time.
And S120, constructing a correlation function of the transmitted signal and the received signal, and taking the time corresponding to a plurality of peaks higher than a peak threshold in the correlation function as a plurality of arrival times.
In some embodiments, constructing a correlation function of a transmitted signal and a received signal comprises:
wherein R (τ) is a correlation function of the transmitted signal and the received signal; s (t) is a transmission signal, and r (t + τ) is a reception signal.
T in r (t + τ) is a time variable, and τ is a time delay.
The image of the correlation function R (τ) is a continuous curve, and τ where the peak occurs is the time of arrival of the multipath signal.
Due to multipath propagation, a plurality of peaks appear in the correlation function, and the time corresponding to the plurality of peaks of which the peaks exceed the peak threshold in the correlation function is extracted as the arrival time of the signal passing through each path.
S130, constructing according to the signals of each subcarrier received by each antenna and the signals of each subcarrier sent by the base station to obtain a channel frequency response matrix, and constructing according to the channel frequency response matrix to obtain a noise matrix.
In some embodiments, the signal of each subcarrier received by each antenna is divided by the signal of each subcarrier transmitted by the base station to construct a channel frequency response matrix, and the CFR represents the channel frequency response matrix:
wherein, the number of rows M is the number of antennas, and the number of columns K is the number of subcarriers occupied by CSI-RS/PRS signals.
The frequency response refers to a characteristic that the amplitude and phase of a signal are changed by frequency change, and the frequency response is composed of an amplitude-frequency characteristic and a phase-frequency characteristic. The amplitude-frequency characteristic represents the relationship between the increase and decrease of the gain and the signal frequency; the phase-frequency characteristic represents the phase distortion relationship at different signal frequencies.
The channel frequency response reflects the state of each channel, that is, how the delay of the channel affects the transmission signal, and is a real signal channel, so that the channel frequency response includes noise.
In some embodiments, an eigenvector and an eigenvalue of a product of a channel frequency response matrix and a conjugate transpose matrix of the channel frequency response matrix are obtained, and a noise matrix is constructed according to an eigenvector corresponding to an eigenvalue lower than an eigenvalue threshold.
In some embodiments, after constructing the channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier transmitted by the base station, the method further includes: processing the channel frequency response matrix by using a smoothing technology to obtain a smoothed channel frequency response matrix;
the method specifically comprises the following steps:
dividing a channel frequency response matrix into 2(M-1) submatrices; wherein M is the number of antennas;
and taking the elements in each sub-matrix as one column of the smoothed channel frequency response matrix to construct the smoothed channel frequency response matrix.
The dividing method comprises the following steps:
wherein, the first sub-matrix comprises all elements except the last column in the first row and the second row in the CFR matrix; the second sub-matrix comprises all elements except the first column in the first row and the second row in the CFR matrix; the third sub-matrix comprises all elements except the last column in the second row and the third row in the CFR matrix; the fourth sub-matrix comprises all elements except the first column in the second row and the third row in the CFR matrix; by analogy, … …, the 2(M-1) -1 sub-matrix includes all the elements except the last column in the (M-1) th row and the M-th row in the CFR matrix; the 2 nd (M-1) sub-matrix includes all elements except the first column in the (M-1) th row and the M th row in the CFR matrix.
The resulting smoothed channel frequency response matrix is:
the smoothing of the CFR is to increase the measurement, which is equivalent to extending the antenna array, and the increase of the measurement can improve the accuracy.
In some embodiments, an eigenvector and an eigenvalue of a product of a smoothed channel frequency response matrix and a conjugate transpose matrix of the smoothed channel frequency response matrix are obtained, and a noise matrix is constructed according to an eigenvector corresponding to an eigenvalue lower than an eigenvalue threshold.
Is CFRsmothConjugate transpose matrix of (CFR)smothWith its conjugate transpose matrixBy performing matrix multiplication, i.e.
ObtainingAccording to the eigenvector corresponding to the eigenvalue lower than the eigenvalue threshold, a noise matrix is constructed, and Q is usednRepresenting a noise matrix.
The characteristic value lower than the characteristic value threshold value corresponds to a noise space, and the characteristic value not lower than the characteristic value threshold value corresponds to a signal space.
And S140, constructing and obtaining a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier.
The subcarriers with different frequencies can generate phase difference when reaching the same antenna, and the subcarriers with the same frequency can also generate phase difference when reaching different antennas, and based on the phase difference, a steering vector matrix can be constructed to be:
wherein the content of the first and second substances,is a steering vector matrix, theta is an arrival angle, and tau is an arrival time; t is the phase difference of different sub-carriers to the same antenna, phi is the phase difference of different antennas to the same sub-carrier, K is the number of sub-carriers, and M is the number of antennas.
Where f is the frequency of the subcarrier, Δ f is the interval of the subcarrier, c is the speed of light, d is the distance between two adjacent antennas, and j is the imaginary unit.
The steering vector matrix is constructed from the antenna structure and the structure of the transmitted signal in combination with the time of arrival and the angle of arrival, and is understood to be a function related to θ, τ (θ, τ are independent variables), which is ideally noise-free.
S150, constructing according to the noise matrix and the steering vector matrix to obtain a 2D-MUSIC spectrum function, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
The 2D-MUSIC spectral function is:
wherein P (theta, tau) is a 2D-MUSIC spectral function;in order to be a matrix of steering vectors,as a conjugate transpose of the steering vector matrix, QnIn the form of a noise matrix, the noise matrix,is the conjugate transpose of the noise matrix.
And fixing the arrival time tau, wherein the fixed value is the signal arrival time, searching the arrival angle theta, obtaining the arrival angle matched with the path by fixing each arrival time, and finally obtaining the arrival time and the arrival angle matched with all paths.
The feature value lower than the feature value threshold value corresponds to a noise space, the feature value not lower than the feature value threshold value corresponds to a signal space, the noise space and the signal space are orthogonal, and the product of the orthogonal terms is zero, so when the denominator of the spectrum function tends to be zero due to a certain pair of arrival time theta and the arrival angle tau, namely a spectrum peak appears, that is, a certain pair of arrival time and arrival angle are solved.
As can be seen from the foregoing, the method and related apparatus for obtaining time of arrival and angle of arrival of multipath signals according to the present disclosure receive signals transmitted by a base station and arriving through multiple paths through multiple antennas, and obtain the time of arrival of multipath signals according to the correlation between the transmitted signals and the received signals; the method comprises the steps of constructing a channel frequency response matrix according to signals of each subcarrier received by each antenna and signals of each subcarrier sent by a base station, constructing a noise matrix according to a smooth channel frequency response matrix, constructing a steering vector matrix according to phase differences of different subcarriers for the same antenna and phase differences of different antennas for the same subcarrier, constructing a 2D-MUSIC spectrum function by using the steering vector matrix and the noise matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function, thereby obtaining the arrival time and the arrival angle of matched multipath signals and providing accurate parameters for positioning wireless signals.
The method is based on a 5G communication system, in the scene of signal multipath propagation such as indoor or urban canyons and the like, the antenna array is expanded by utilizing the subcarrier modulated by the 5G signal OFDM and the spatial smoothing technology, and the arrival time and the arrival angle of a downlink multipath signal are jointly estimated by utilizing the combination of a correlation method and a traditional 2D-MUSIC algorithm. The method and the device are different from the conventional direct path signal parameter estimation, can estimate the multipath signal parameters, and can match the estimated arrival time and angle with each other, thereby providing a signal parameter basis for a system for positioning by using the multipath wireless signals.
It should be noted that the method of the embodiments of the present disclosure may be executed by a single device, such as a computer or a server. The method of the embodiment can also be applied to a distributed scene and completed by the mutual cooperation of a plurality of devices. In such a distributed scenario, one of the devices may only perform one or more steps of the method of the embodiments of the present disclosure, and the devices may interact with each other to complete the method.
It should be noted that the above describes some embodiments of the disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments described above and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
Based on the same inventive concept, corresponding to the method of any embodiment, the disclosure also provides a device for acquiring the arrival time and the arrival angle of the multipath signal.
Referring to fig. 3, the apparatus for acquiring time of arrival and angle of arrival of a multipath signal includes:
a signal obtaining module 310, configured to obtain a signal sent by a base station as a sending signal, and obtain signals received by multiple antennas as receiving signals;
an arrival time determining module 320, configured to construct a correlation function between the transmitted signal and the received signal, and take times corresponding to a plurality of peaks higher than a peak threshold in the correlation function as a plurality of arrival times;
a noise matrix constructing module 330, configured to obtain a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station, and obtain a noise matrix according to the channel frequency response matrix;
a steering vector matrix constructing module 340, configured to construct a steering vector matrix according to a phase difference of different subcarriers for the same antenna and a phase difference of different antennas for the same subcarrier;
and an arrival angle determining module 350, configured to construct a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtain an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, the functionality of the various modules may be implemented in the same one or more software and/or hardware implementations of the present disclosure.
The apparatus in the foregoing embodiment is used to implement the method for acquiring time of arrival and angle of arrival of a multipath signal in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which are not described herein again.
Based on the same inventive concept, corresponding to the method of any embodiment described above, the present disclosure further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor executes the program to implement the method for acquiring the time of arrival and the angle of arrival of the multipath signal according to any embodiment described above.
Fig. 4 is a schematic diagram illustrating a more specific hardware structure of an electronic device according to this embodiment, where the electronic device may include: a processor 1010, a memory 1020, an input/output interface 1030, a communication interface 1040, and a bus 1050. Wherein the processor 1010, memory 1020, input/output interface 1030, and communication interface 1040 are communicatively coupled to each other within the device via bus 1050.
The processor 1010 may be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an Application Specific Integrated Circuit (ASIC), or one or more Integrated circuits, and is configured to execute related programs to implement the technical solutions provided in the embodiments of the present disclosure.
The Memory 1020 may be implemented in the form of a ROM (Read Only Memory), a RAM (Random Access Memory), a static storage device, a dynamic storage device, or the like. The memory 1020 may store an operating system and other application programs, and when the technical solution provided by the embodiments of the present specification is implemented by software or firmware, the relevant program codes are stored in the memory 1020 and called to be executed by the processor 1010.
The input/output interface 1030 is used for connecting an input/output module to input and output information. The i/o module may be configured as a component in a device (not shown) or may be external to the device to provide a corresponding function. The input devices may include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output devices may include a display, a speaker, a vibrator, an indicator light, etc.
The communication interface 1040 is used for connecting a communication module (not shown in the drawings) to implement communication interaction between the present apparatus and other apparatuses. The communication module can realize communication in a wired mode (such as USB, network cable and the like) and also can realize communication in a wireless mode (such as mobile network, WIFI, Bluetooth and the like).
It should be noted that although the above-mentioned device only shows the processor 1010, the memory 1020, the input/output interface 1030, the communication interface 1040 and the bus 1050, in a specific implementation, the device may also include other components necessary for normal operation. In addition, those skilled in the art will appreciate that the above-described apparatus may also include only those components necessary to implement the embodiments of the present description, and not necessarily all of the components shown in the figures.
The electronic device of the foregoing embodiment is used to implement the method for acquiring time of arrival and angle of arrival of a multipath signal in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiment, which are not described herein again.
Based on the same inventive concept, the present disclosure also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the multipath signal time-of-arrival and angle-of-arrival acquisition method according to any of the above embodiments, corresponding to any of the above embodiment methods.
Computer-readable media of the present embodiments, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device.
The computer instructions stored in the storage medium of the above embodiment are used to enable the computer to execute the method for acquiring time of arrival and angle of arrival of a multipath signal according to any of the above embodiments, and have the beneficial effects of corresponding method embodiments, and are not described herein again.
It should be noted that the embodiments of the present disclosure can be further described in the following ways:
a method for acquiring the time of arrival and the angle of arrival of a multipath signal comprises the following steps:
acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
constructing a correlation function of the transmitting signal and the receiving signal, and taking the time corresponding to a plurality of peak values which are higher than a peak value threshold value in the correlation function as a plurality of arrival times;
constructing a channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier sent by the base station, and constructing a noise matrix according to the channel frequency response matrix;
constructing a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and constructing a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
Optionally, the plurality of antennas are arranged in a row, a distance between every two adjacent antennas is the same, and the distance is greater than half of a wavelength of the transmission signal.
Optionally, wherein the transmission signal reaches the plurality of antennas through a plurality of paths;
the received signal is a superposition of the transmitted signals that arrive at the plurality of antennas via the plurality of paths.
Optionally, the constructing a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station includes:
and dividing the signal of each subcarrier received by each antenna with the signal of each subcarrier sent by the base station to construct a channel frequency response matrix.
Optionally, after constructing a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station, the method further includes: processing the channel frequency response matrix by using a smoothing technology to obtain a smoothed channel frequency response matrix;
the method specifically comprises the following steps:
dividing the channel frequency response matrix into 2(M-1) sub-matrices; wherein M is the number of the antennas;
and taking the elements in each sub-matrix as a column of the smooth channel frequency response matrix to construct the smooth channel frequency response matrix.
Optionally, the constructing a noise matrix according to the channel frequency response matrix includes:
and acquiring an eigenvector and an eigenvalue of a product of the channel frequency response matrix and a conjugate transpose matrix of the channel frequency response matrix, and constructing the noise matrix according to the eigenvector corresponding to the eigenvalue lower than an eigenvalue threshold.
Optionally, the constructing a noise matrix according to the channel frequency response matrix includes:
and acquiring an eigenvector and an eigenvalue of a product of the smoothed channel frequency response matrix and a conjugate transpose matrix of the smoothed channel frequency response matrix, and constructing and acquiring the noise matrix according to the eigenvector corresponding to the eigenvalue lower than an eigenvalue threshold.
A multipath signal time-of-arrival and angle-of-arrival acquisition apparatus, comprising:
the signal acquisition module is used for acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
an arrival time determining module, configured to construct a correlation function between the transmitted signal and the received signal, and use, as multiple arrival times, times corresponding to multiple peaks higher than a peak threshold in the correlation function;
a noise matrix constructing module, configured to construct a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station, and construct a noise matrix according to the channel frequency response matrix;
the steering vector matrix constructing module is used for constructing and obtaining a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and the arrival angle determining module is used for constructing and obtaining a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method as described above when executing the program.
A non-transitory computer readable storage medium storing computer instructions for causing a computer to perform the above method.
Those of ordinary skill in the art will understand that: the discussion of any embodiment above is meant to be exemplary only, and is not intended to intimate that the scope of the disclosure, including the claims, is limited to these examples; within the idea of the present disclosure, also technical features in the above embodiments or in different embodiments may be combined, steps may be implemented in any order, and there are many other variations of the different aspects of the embodiments of the present disclosure as described above, which are not provided in detail for the sake of brevity.
In addition, well-known power/ground connections to Integrated Circuit (IC) chips and other components may or may not be shown in the provided figures for simplicity of illustration and discussion, and so as not to obscure the embodiments of the disclosure. Furthermore, devices may be shown in block diagram form in order to avoid obscuring embodiments of the present disclosure, and this also takes into account the fact that specifics with respect to implementation of such block diagram devices are highly dependent upon the platform within which the embodiments of the present disclosure are to be implemented (i.e., specifics should be well within purview of one skilled in the art). Where specific details (e.g., circuits) are set forth in order to describe example embodiments of the disclosure, it should be apparent to one skilled in the art that the embodiments of the disclosure can be practiced without, or with variation of, these specific details. Accordingly, the description is to be regarded as illustrative instead of restrictive.
While the present disclosure has been described in conjunction with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art in light of the foregoing description. For example, other memory architectures (e.g., dynamic ram (dram)) may use the discussed embodiments.
The disclosed embodiments are intended to embrace all such alternatives, modifications and variances which fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalents, improvements, and the like that may be made within the spirit and principles of the embodiments of the disclosure are intended to be included within the scope of the disclosure.
Claims (10)
1. A method for acquiring the time of arrival and the angle of arrival of a multipath signal comprises the following steps:
acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
constructing a correlation function of the transmitting signal and the receiving signal, and taking the time corresponding to a plurality of peak values which are higher than a peak value threshold value in the correlation function as a plurality of arrival times;
constructing a channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier sent by the base station, and constructing a noise matrix according to the channel frequency response matrix;
constructing a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and constructing a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
2. The method of claim 1, wherein the plurality of antennas are arranged in a row, a spacing between every two adjacent antennas is the same, and the spacing is greater than half a wavelength of the transmission signal.
3. The method of claim 1, wherein the transmit signal travels a plurality of paths to the plurality of antennas;
the received signal is a superposition of the transmitted signals that arrive at the plurality of antennas via the plurality of paths.
4. The method of claim 1, wherein the constructing a channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier transmitted by the base station comprises:
and dividing the signal of each subcarrier received by each antenna with the signal of each subcarrier sent by the base station to construct a channel frequency response matrix.
5. The method of claim 1, wherein after constructing a channel frequency response matrix according to the signals of each subcarrier received by each antenna and the signals of each subcarrier transmitted by the base station, the method further comprises: processing the channel frequency response matrix by using a smoothing technology to obtain a smoothed channel frequency response matrix;
the method specifically comprises the following steps:
dividing the channel frequency response matrix into 2(M-1) sub-matrices; wherein M is the number of the antennas;
and taking the elements in each sub-matrix as a column of the smooth channel frequency response matrix to construct the smooth channel frequency response matrix.
6. The method of claim 1, wherein said constructing a noise matrix from said channel frequency response matrix comprises:
and acquiring an eigenvector and an eigenvalue of a product of the channel frequency response matrix and a conjugate transpose matrix of the channel frequency response matrix, and constructing the noise matrix according to the eigenvector corresponding to the eigenvalue lower than an eigenvalue threshold.
7. The method of claim 5, wherein said constructing a noise matrix from said channel frequency response matrix comprises:
and acquiring an eigenvector and an eigenvalue of a product of the smoothed channel frequency response matrix and a conjugate transpose matrix of the smoothed channel frequency response matrix, and constructing and acquiring the noise matrix according to the eigenvector corresponding to the eigenvalue lower than an eigenvalue threshold.
8. A multipath signal time-of-arrival and angle-of-arrival acquisition apparatus, comprising:
the signal acquisition module is used for acquiring signals sent by a base station as sending signals and acquiring signals received by a plurality of antennas as receiving signals;
an arrival time determining module, configured to construct a correlation function between the transmitted signal and the received signal, and use, as multiple arrival times, times corresponding to multiple peaks higher than a peak threshold in the correlation function;
a noise matrix constructing module, configured to construct a channel frequency response matrix according to the signal of each subcarrier received by each antenna and the signal of each subcarrier sent by the base station, and construct a noise matrix according to the channel frequency response matrix;
the steering vector matrix constructing module is used for constructing and obtaining a steering vector matrix according to the phase difference of different subcarriers to the same antenna and the phase difference of different antennas to the same subcarrier;
and the arrival angle determining module is used for constructing and obtaining a 2D-MUSIC spectrum function according to the noise matrix and the steering vector matrix, and obtaining an arrival angle corresponding to each arrival time according to a plurality of arrival times based on the 2D-MUSIC spectrum function.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method of any one of claims 1 to 7 when executing the program.
10. A non-transitory computer readable storage medium storing computer instructions for causing a computer to perform the method of any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110553482.3A CN113286362A (en) | 2021-05-20 | 2021-05-20 | Method for acquiring arrival time and arrival angle of multipath signal and related device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110553482.3A CN113286362A (en) | 2021-05-20 | 2021-05-20 | Method for acquiring arrival time and arrival angle of multipath signal and related device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113286362A true CN113286362A (en) | 2021-08-20 |
Family
ID=77280471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110553482.3A Pending CN113286362A (en) | 2021-05-20 | 2021-05-20 | Method for acquiring arrival time and arrival angle of multipath signal and related device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113286362A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113939014A (en) * | 2021-09-24 | 2022-01-14 | 北京邮电大学 | Wireless indoor positioning method and positioning device based on channel state information and related equipment |
CN114375038A (en) * | 2022-01-12 | 2022-04-19 | 网络通信与安全紫金山实验室 | Positioning method, positioning device, base station, storage medium and computer program product |
CN114679356A (en) * | 2022-03-17 | 2022-06-28 | 西安电子科技大学 | Channel full-dimensional parameter extraction method independent of likelihood function |
WO2024016708A1 (en) * | 2022-07-19 | 2024-01-25 | 中兴通讯股份有限公司 | Sensing method, receiver, and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387864A (en) * | 2018-03-02 | 2018-08-10 | 北京邮电大学 | A kind of angle of arrival computational methods and device |
CN108646213A (en) * | 2018-05-09 | 2018-10-12 | 华南理工大学 | Direct wave AOA determination methods under a kind of indoor multipath environment |
WO2021032307A1 (en) * | 2019-08-22 | 2021-02-25 | Huawei Technologies Co., Ltd. | Radio station for client localization in multipath indoor environment |
-
2021
- 2021-05-20 CN CN202110553482.3A patent/CN113286362A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387864A (en) * | 2018-03-02 | 2018-08-10 | 北京邮电大学 | A kind of angle of arrival computational methods and device |
CN108646213A (en) * | 2018-05-09 | 2018-10-12 | 华南理工大学 | Direct wave AOA determination methods under a kind of indoor multipath environment |
WO2021032307A1 (en) * | 2019-08-22 | 2021-02-25 | Huawei Technologies Co., Ltd. | Radio station for client localization in multipath indoor environment |
Non-Patent Citations (1)
Title |
---|
吴志国: "基于WiFi信道状态信息的单基站室内定位研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113939014A (en) * | 2021-09-24 | 2022-01-14 | 北京邮电大学 | Wireless indoor positioning method and positioning device based on channel state information and related equipment |
CN113939014B (en) * | 2021-09-24 | 2022-09-02 | 北京邮电大学 | Wireless indoor positioning method based on channel state information and related equipment |
CN114375038A (en) * | 2022-01-12 | 2022-04-19 | 网络通信与安全紫金山实验室 | Positioning method, positioning device, base station, storage medium and computer program product |
CN114375038B (en) * | 2022-01-12 | 2024-04-26 | 网络通信与安全紫金山实验室 | Positioning method, apparatus, base station, storage medium and computer program product |
CN114679356A (en) * | 2022-03-17 | 2022-06-28 | 西安电子科技大学 | Channel full-dimensional parameter extraction method independent of likelihood function |
WO2024016708A1 (en) * | 2022-07-19 | 2024-01-25 | 中兴通讯股份有限公司 | Sensing method, receiver, and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113286362A (en) | Method for acquiring arrival time and arrival angle of multipath signal and related device | |
KR102356954B1 (en) | Method and apparatus for estimating position of terminal | |
EP2798366B1 (en) | Method and system for localization | |
CN108387864B (en) | Method and device for calculating angle of arrival | |
US20150087331A1 (en) | Method and apparatus of wi-fi-based positioning | |
US9763216B2 (en) | Radiator localization | |
US10823835B2 (en) | Electromagnetic reflection profiles | |
WO2018119949A1 (en) | Channel state information phase correction method and apparatus | |
KR20210060844A (en) | Method and apparatus for measuting position of three dimensions using radar sensor | |
Chen et al. | Signal acquisition of Luojia-1A low earth orbit navigation augmentation system with software defined receiver | |
US7414582B1 (en) | Method and apparatus for all-polarization direction finding | |
CN114779297A (en) | Terminal positioning method and related device | |
JP6331072B2 (en) | White space detection device, white space detection method, and program | |
Yang et al. | Central DOA estimation of incoherently distributed noncircular sources with cross-correlation matrix | |
CN116033339A (en) | Information reporting method, device, equipment and readable storage medium | |
Tian et al. | SAPIL: single access point based indoor localisation using Wi‐Fi L‐shaped antenna array | |
CN113939014B (en) | Wireless indoor positioning method based on channel state information and related equipment | |
JP2004080604A (en) | Method for estimating space angle spread and receiver | |
CN113766424B (en) | Unmanned aerial vehicle indoor positioning method, unmanned aerial vehicle indoor positioning device, computer equipment and storage medium | |
JP6311198B2 (en) | White space detection device, white space detection method, and program | |
CN117538854B (en) | Ranging method, ranging apparatus, computer device, and computer-readable storage medium | |
CN116669181B (en) | Indoor personnel positioning method and system based on WiFi multi-reflection path image | |
WO2024131752A1 (en) | Sensing method, sensing apparatus, communication device and storage medium | |
Yildirim et al. | Impact of MU-MIMO on Passive Wi-Fi Sensing: Threat or Opportunity? | |
US20230082198A1 (en) | Positioning Measurements |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210820 |