Micro-deformation monitoring method, electronic device, and storage medium

Abstract

A micro-deformation monitoring method, comprising: on the basis of received data of a wireless communication network, determining a channel impulse response (110); on the basis of the channel impulse response, generating a range-angle map (120); on the basis of the range-angle map, determining a phase difference between a monitoring target point and a reference point (130); and on the basis of changes in the phase difference, determining the deformation of the monitoring target point (140). Further provided are an electronic device and a computer-readable storage medium.

Description

Micro-deformation monitoring method, electronic device and storage medium Technical Field

[0001] The present application relates to the field of wireless communication technology, for example, to a method for monitoring micro-deformation, an electronic device, and a storage medium. Background Art

[0002] Micro-deformation monitoring technology can monitor real-time changes in topography and provide early warning signals before geological disasters occur, helping people take preventative and response measures. This can reduce losses from natural disasters such as landslides, debris flows, and collapses. Micro-deformation monitoring technology can also monitor bridges, high-rise buildings, subways, and large venues, enabling timely identification of potential hazards and safeguarding public life and property.

[0003] Currently, micro-deformation monitoring technology is mainly implemented through sensors, such as electronic levels, dynamic weight measuring instruments, inclinometers, displacement meters, strain gauges, and dynamometers. These sensors monitor the micro-deformation of the monitored object through various sensors such as humidity, pressure, and stress. The greatest advantage of physical sensors is that they can accurately monitor some details, such as internal stress and pressure, and can monitor local deformation information with high precision. However, they are not accurate enough for monitoring overall deformation. Sensorless micro-deformation monitoring methods include optical and electronic measuring instruments, such as total stations, electromagnetic rangefinders, theodolites, and levels. These periodically measure monitoring points at specific distances and angles, obtain the three-dimensional coordinates of the measurement points, and determine whether the observed object has undergone horizontal or vertical displacement based on subsequent measurements. These methods offer high flexibility and measurement accuracy, but require specialized equipment for real-time monitoring of the monitored object, which is costly. Currently, monitoring can also be performed using Global Positioning System (GPS) signals. GPS measurement technology has the advantage of high precision in deformation monitoring, but GPS has poor penetration and is easily obstructed, which can affect monitoring accuracy. Micro-deformation monitoring using satellite remote sensing technology offers high timeliness and accuracy. However, common satellite remote sensing sensors, including aperture radar, lidar, and optical satellite sensors, are significantly affected by weather conditions. When obscured by clouds, fog, or water vapor, the quality of monitored images degrades significantly. Therefore, current micro-deformation monitoring methods urgently require a low-cost, high-quality, high-precision monitoring technology with extensive coverage.

[0004] Summary of the Invention

[0005] The embodiments of the present application provide a method, electronic device and storage medium for monitoring micro-deformations to solve the problem that the monitoring method is greatly affected by the environment, achieve high-precision detection of micro-deformations of the monitored target, increase the stability of micro-deformation monitoring, and improve the monitoring range to achieve overall detection of the deformation of the monitored object.

[0006] The present invention provides a method for monitoring micro deformation, wherein the method includes:

[0007] determining a channel impulse response based on received data from a wireless communication network;

[0008] generating a range angle graph according to the channel impulse response;

[0009] Determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram;

[0010] The deformation amount of the monitoring target point is determined according to the change of the phase difference.

[0011] An embodiment of the present application further provides an electronic device, wherein the electronic device includes:

[0012] one or more processors;

[0013] a memory for storing one or more programs;

[0014] When the one or more programs are executed by the one or more processors, the one or more processors implement the micro-deformation monitoring method as described in any one of the embodiments of the present application.

[0015] An embodiment of the present application further provides a computer-readable storage medium storing one or more programs, which are executed by one or more processors to implement the micro-deformation monitoring method as described in any one of the embodiments of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] FIG1 is a schematic diagram of a micro-deformation monitoring system provided in an embodiment of the present application;

[0017] FIG2 is a flow chart of a method for monitoring micro-deformation provided in an embodiment of the present application;

[0018] FIG3 is a flow chart of another method for monitoring micro-deformation provided in an embodiment of the present application;

[0019] FIG4 is a flow chart of another method for monitoring micro-deformation provided in an embodiment of the present application;

[0020] FIG5 is an example diagram of a method for monitoring micro-deformation provided in an embodiment of the present application;

[0021] FIG6 is a diagram of a channel impulse response provided in an embodiment of the present application;

[0022] FIG7 is an example diagram of a reference point for determining a LOS path on an RA diagram provided by an embodiment of the present application;

[0023] FIG8 is an RA diagram of an environmental air mining provided in an embodiment of the present application;

[0024] FIG9 is an RA diagram of a corner reflector placement provided in an embodiment of the present application;

[0025] FIG10 is an example diagram of a phase difference time curve provided in an embodiment of the present application;

[0026] FIG11 is an RA diagram of another environmental air mining provided in an embodiment of the present application;

[0027] FIG12 is an RA diagram of another embodiment of the present application for placing a corner reflector;

[0028] FIG13 is an example diagram of another phase difference time curve provided in an embodiment of the present application;

[0029] FIG14 is a schematic structural diagram of a micro-deformation monitoring device provided in an embodiment of the present application;

[0030] FIG15 is a schematic structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0031] It should be understood that the specific implementations described herein are only used to explain the present application and are not used to limit the present application.

[0032] In the subsequent description, suffixes such as "module", "component" or "unit" used to represent elements are only used to facilitate the description of this application and have no specific meaning. Therefore, "module", "component" or "unit" can be used interchangeably.

[0033] Figure 1 is a schematic diagram of a micro-deformation monitoring provided by an embodiment of the present application. Referring to Figure 1, the embodiment of the present application transmits a monitoring signal through a fifth-generation mobile communication technology (5th Generation, 5G) base station deployed on the ground, calculates the echo signal phase of the monitored target in real time, and determines the deformation amount of the detection target based on the echo signal phase change, thereby achieving millimeter-level deformation monitoring accuracy and having the advantages of low cost and wide coverage.

[0034] Figure 2 is a flow chart of a method for monitoring micro deformation provided by an embodiment of the present application. The embodiment of the present application is applicable to the case where micro deformation of an object is monitored based on a base station. The method can be executed by a micro deformation monitoring device, which can be implemented by software and / or hardware methods and can generally be integrated into a base station or terminal device. As shown in Figure 2, the method provided by the embodiment of the present application specifically includes the following steps:

[0035] Step 110: Determine a channel impulse response based on received data from the wireless communication network.

[0036] Among them, the wireless communication network can be a network that communicates with the help of wireless signals, and the wireless communication network includes but is not limited to 5G communication network, near-field communication network, etc. The received data can be a signal received through the wireless communication network. The received data can be used to determine the phase of the echo signal. The received data can include a data matrix determined by the number of subcarriers and the number of antennas. The channel impulse response can be the response of a unit pulse signal after being transmitted through the wireless communication network.

[0037] In an embodiment of the present application, received data transmitted through a wireless communication network can be obtained, and a channel impulse response can be generated by processing the received data. Specifically, an inverse Fourier transform can be performed on the obtained received data, so that the result generated by processing the received data is used as the channel impulse response.

[0038] Step 120: Generate a range angle diagram according to the channel impulse response.

[0039] Among them, the distance angle diagram can be an image constructed with the signal arrival azimuth and distance as dimensions respectively. The distance angle diagram can include multiple dimensions, including but not limited to distance dimension, azimuth dimension, pitch angle dimension, etc.

[0040] In an embodiment of the present application, the channel impulse response may include signal paths at different distances, and the signal path at each distance is generated by superimposing signals at different angles. The channel impulse response can be separated in the angular dimension, and the signal paths at different distances and angles are determined to form a distance angle diagram. In an embodiment of the present application, the method of separating the channel impulse response in the angular dimension includes but is not limited to a spectrum peak search method, a super-resolution algorithm, and the like.

[0041] Step 130: Determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram.

[0042] The monitoring target point can be a target that requires micro-deformation monitoring. The monitoring target point can include large venues, mountains, bridges, subways, bridges, etc. The reference point can be a location that assists in determining the micro-deformation of the monitoring target point. The reference point can include a base station with a line of sight (LOS) path to the base station that receives the received data or a path in an environment with a high impulse response. The reference point can be pre-configured or automatically selected according to preset rules.

[0043] In an embodiment of the present application, corresponding position points can be determined in a distance angle diagram by respectively measuring the distance and phase of the monitoring target point and the reference point relative to the device that receives the received data, and the phase difference between the monitoring target point and the reference point can be determined by using the phase corresponding to the position point. It is understood that the distance and phase of the monitoring target point and the reference point relative to the base station that receives the received data can be known in advance, and the position points corresponding to the monitoring target point and the reference point can be obtained in the distance angle diagram by using the distance and phase.

[0044] Step 140: Determine the deformation amount of the monitoring target point according to the change of the phase difference.

[0045] The change condition may be the change of the phase difference over time, which may be determined by the difference in the phase difference at different moments, and the deformation variable may be the distance over which the monitoring target point is displaced, which may be determined by the change condition.

[0046] In an embodiment of the present application, the phase difference of the monitoring target point relative to the reference point obtained in the aforementioned step can be saved. Changes can be determined by comparing the phase differences saved at different moments. This change may include, but is not limited to, changes in phase differences between adjacent moments and changes in phase differences at each moment relative to the first recorded phase difference. The deformation of the monitoring target point can be calculated based on the changes in phase differences and actual physical quantities of the wireless communication network. For example, the displacement corresponding to the changes can be determined using the wavelength of the wireless communication network, and this displacement can be used as the deformation.

[0047] In an embodiment of the present application, received data is obtained in a wireless communication network, and a channel impulse response corresponding to the received data is determined, the channel impulse response is processed into a distance angle diagram, and the phase difference between the target point and the reference point is determined based on the distance angle diagram. The deformation of the monitored target point is determined by the change in the phase difference, thereby realizing high-precision deformation monitoring of the monitored target point. The change in the phase difference can be determined by the received data of the wireless communication network, which can reduce the impact of the environment on the monitoring of small deformations and enhance the robustness of the small deformation monitoring process. With the help of multi-directional signals in the received data, all-round monitoring of the monitored target point is realized, the monitoring range is expanded, and small deformations are monitored through existing commercial base stations without the need to add additional equipment, thereby reducing the monitoring cost.

[0048] FIG3 is a flow chart of another method for monitoring micro-deformation provided in an embodiment of the present application. The embodiment of the present application is a refinement of the above-mentioned embodiment. Referring to FIG3 , the method provided in the embodiment of the present application specifically includes the following steps:

[0049] Step 210: Perform delay calibration and phase calibration on the received data.

[0050] In an embodiment of the present application, after receiving received data from a wireless communication network, due to the existence of delay and phase drift in the wireless communication network, the received data can be subjected to delay calibration and phase calibration in sequence. The delay calibration can be to compensate for the delay difference of the data received by each receiving antenna in the received data using the data corresponding to each receiving antenna, and the phase calibration can be to compensate for the phase difference of the data corresponding to each receiving antenna in the received data.

[0051] In some application embodiments, phase calibration of received data includes: determining phase differences between other symbol data in the received data and the first symbol data; and compensating the phase differences to the corresponding other symbol data to complete the phase calibration of the received data. The other symbol data may be data transmitted by orthogonal frequency division multiplexing (OFDM) symbols other than the first OFDM symbol in the received data, the first symbol data may be data transmitted by the first OFDM symbol in the received data, and the phase difference may be the phase difference between data transmitted on the other OFDM symbols other than the first OFDM symbol and the first OFDM symbol.

[0052] In an embodiment of the present application, the data transmitted by each OFDM symbol in the received data can be determined, and the phase differences corresponding to the data transmitted on the other OFDM symbols except the first OFDM symbol and the first OFDM symbol can be determined respectively. Each other symbol data is compensated using the phase difference corresponding to its OFDM symbol, and the phase calibration of the received data is completed after all other symbol data are compensated for the phase difference.

[0053] In other application embodiments, performing delay calibration on received data includes:

[0054] Determine the time delay differences between other symbol data and the first symbol data in the received data; and compensate each time delay difference to the corresponding other symbol data to complete the time delay calibration of the received data.

[0055] In an embodiment of the present application, the data transmitted by each OFDM symbol in the received data can be determined, and the delay of other OFDM symbols except the first OFDM symbol can be compared with the data transmitted on the first OFDM symbol, so as to determine the delay difference between each other symbol data and the first symbol data, and the delay difference corresponding to each other symbol data is used to compensate for the delay difference of its OFDM symbol, thereby realizing delay calibration of the received data.

[0056] Step 220: Generate an inverse Fourier transform processing result of the received data, and use the processing result as a channel impulse response.

[0057] Among them, the Inverse Fast Fourier Transform (IFFT) can convert the received data from the signal frequency domain expression to the signal time domain expression.

[0058] In an embodiment of the present application, the received data after delay calibration and phase calibration can be subjected to an inverse Fourier transform so that the received data is converted into a signal time domain expression. The processing result generated after processing can reflect the strength of the path of different signals changing over time, and the processing result generated by the inverse Fourier transform can be used as a channel impulse response.

[0059] Step 230: Construct horizontal and vertical steering vectors for any direction angle and any pitch angle, and determine the Kronecker product of the horizontal steering vector and the vertical steering vector.

[0060] In an embodiment of the present application, since the channel impulse response is the superposition of signals at different angles, in order to detect the monitoring target point, the channel impulse response can be separated according to different angles, and the horizontal and vertical steering vectors can be constructed based on any direction angle and any pitch angle, and the horizontal and vertical steering vectors can be subjected to Kronecker product operations.

[0061] Step 240: Determine the modulus of the product of the Kronecker product and the channel impulse response, and construct a range-angle graph based on the modulus.

[0062] In an embodiment of the present application, the Kronecker product of the horizontal steering vector and the vertical steering vector can be multiplied by the channel impulse response to determine the product, thereby separating the radial strength superimposed at different angles in the channel impulse response into different directional angles and pitch angles. The modulus of the product may include the signal strength or radial strength corresponding to different directional angles, different pitch angles and different distances. Based on the correspondence between the above-mentioned directional angles, pitch angles, distances and signal strengths, a distance angle graph is constructed, and the distance angle graph may include at least three-dimensional information such as distance, directional angle and pitch angle.

[0063] Step 250: Determine first position information of the monitoring target point in the distance-angle graph according to the distance information and angle information of the monitoring target point.

[0064] Among them, the distance information and angle information can be the position information of the monitoring target point in space. The distance information and angle information of the monitoring target point can be pre-measured or pre-configured. It can be understood that the angle information can include direction angle and pitch angle, etc.

[0065] In an embodiment of the present application, the distance information and angle information of the pre-measured or pre-configured monitoring target point can be extracted, and the first position information corresponding to the monitoring target point can be determined in the distance-angle diagram according to the distance information and angle information. The first position information can include multiple dimensional coordinates, and these dimensional coordinates can correspond to the distance information and angle information of the monitoring target point.

[0066] Step 260: Determine the second position information of the reference point in the distance angle diagram, and determine the phase difference according to the first position information and the second position information.

[0067] The reference point may be a standard reference point selected in the distance angle diagram, and the reference point may include the LOS diameter in the distance angle diagram or a diameter whose diameter intensity is greater than the diameter intensity of the target detection point.

[0068] Specifically, the reference point can be determined according to the reference point in the distance angle diagram, and the second position information of the reference point can be extracted from the distance angle diagram. The difference between the angle in the first position information and the angle in the second position information can be used as the phase difference between the monitoring target point and the reference point.

[0069] Step 270: Determine the change in the phase difference between the two most recent moments.

[0070] In the embodiment of the present application, the phase difference between the monitoring target point and the reference point can be determined at different times, and the phase difference can be saved. When determining the change amount, the phase difference between the monitoring target point and the reference point determined at the two most recent times can be selected based on the current time, and the difference between the phase differences at the two most recent times can be used as the change amount. It is understood that the difference between the phase difference at the later time and the phase difference at the previous time can be used as the change amount.

[0071] Step 280: Determine a one-way distance change corresponding to the change amount according to the carrier wavelength and carrier phase change corresponding to the received data as a deformation amount.

[0072] Among them, the carrier wavelength and carrier phase change can be the physical quantities of the carrier signal when the wireless communication network transmits and receives data, the carrier wavelength can be the carrier wavelength of the beam transmitting and receiving data, and the carrier phase change can be the phase corresponding to the change of one wavelength of the beam transmitting and receiving data, and the one-way distance change can be the position offset corresponding to a change determined by the carrier wavelength and carrier phase change.

[0073] In an embodiment of the present application, the carrier wavelength and carrier phase change corresponding to the received data transmission can be obtained. The carrier wavelength and carrier phase change can be pre-configured. Since the change amount is actually the phase change amount, the one-way distance change amount corresponding to the phase change amount can be determined through the relationship between the carrier wavelength and the carrier phase change, and the one-way distance change amount can be used as a deformation variable.

[0074] In some application embodiments, determining a one-way distance change corresponding to a change in a carrier wavelength and a carrier phase corresponding to received data as a deformation variable includes:

[0075] Calling formula Determine the deformation amount, where Δd represents the deformation amount, φ represents the phase change, λ represents the carrier wavelength, and the carrier phase change is 2π.

[0076] In the embodiment of the present application, the change amount, carrier wavelength and carrier phase change can be substituted into the formula The one-way distance change is determined as the deformation variable. In the formula, Δd represents the deformation variable, φ represents the change, λ represents the carrier wavelength, and the carrier phase change is 2π.

[0077] In some application embodiments, determining the time delay differences between other symbol data and the first symbol data in the received data includes:

[0078] Calling formula

[0079] Determine the delay difference, where represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, Indicates the i-th symbol of the first sub +M subcarriers corresponding to the received data, Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, and conj() represents a function for determining the conjugate value.

[0080] In the embodiment of the present application, the formula can be called for each other symbol data in the received data Determine the time delay difference between the other symbol data and the first symbol data, where: represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, Indicates the i-th symbol of the first sub +M subcarriers corresponding to the received data, Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, conj() represents a function for determining the conjugate value, and i can be an integer greater than 1.

[0081] In some other application embodiments, compensating each delay difference to corresponding other symbol data to complete delay calibration of received data includes:

[0082] Calling formula Determine the received data after delay compensation, where R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i represents the received data corresponding to the i-th symbol before delay compensation, Indicates the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, and j represents a complex number.

[0083] In the embodiment of the present application, each delay difference is compensated to the corresponding other symbol data respectively, and the compensation method can be obtained by the formula Implementation, where R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i represents the received data corresponding to the i-th symbol before delay compensation, Indicates the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, j represents a complex number, and i can be an integer greater than 1.

[0084] In some other application embodiments, determining the phase difference between each of the other symbols in the received data and the first symbol data includes:

[0085] Calling formula Δφ i =angle(R i *conj(R i,τ )) Determine the phase difference, where Δφ represents the phase difference corresponding to the i-th symbol, R i Indicates the received data corresponding to the i-th symbol before delay compensation, R i,τ represents the received data corresponding to the i-th symbol after delay difference compensation, and conj() represents a function for determining the conjugate value.

[0086] In the embodiment of the present application, for each other symbol, the corresponding received data and the received data after delay calibration can be substituted into the formula angle(R i *conj(R i,τ )) Determine the phase difference, where Δφ represents the phase difference corresponding to the i-th symbol, R i Indicates the received data corresponding to the i-th symbol before delay compensation, R i,τ represents the received data corresponding to the i-th symbol after delay difference compensation, conj() represents a function for determining the conjugate value, and i can be an integer greater than 1.

[0087] In some application embodiments, compensating the phase difference to corresponding other symbol data to complete phase calibration of received data includes:

[0088] Call formula R i,τ,Δφ =R i,τ *exp(j*Δφ i ) determines the phase-calibrated received data, where R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol that has not been phase-calibrated, exp() represents an exponential function with e as the base, j represents a complex number, and i can be an integer greater than 1.

[0089] In the embodiment of the present application, performing phase calibration for each other symbol data may include substituting the received data corresponding to the i-th symbol and the phase difference into the formula R i,τ,Δφ =R i,τ *exp(j*Δφ i) determines the phase-calibrated received data corresponding to the i-th symbol, where R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol that has not been phase-calibrated, exp() represents an exponential function with e as the base, j represents a complex number, and i can be an integer greater than 1.

[0090] Furthermore, based on the above-mentioned application embodiment, the correspondence between the distance in the distance angle diagram and the sampling time of the channel impulse response satisfies the following formula:

[0091] Where d represents the distance, T s represents the sampling time of the channel impulse response, T LOS represents the line-of-sight propagation time, d LOS Indicates the line-of-sight propagation distance.

[0092] FIG4 is a flow chart of another method for monitoring micro-deformation provided in an embodiment of the present application. The embodiment of the present application is a refinement of the above-mentioned embodiment. Referring to FIG4 , the method provided in the embodiment of the present application specifically includes the following steps:

[0093] Step 310: Determine a channel impulse response based on received data from the wireless communication network.

[0094] Step 320: Generate a range angle diagram according to the channel impulse response.

[0095] Step 330: Determine the reference point corresponding to the received data.

[0096] In an embodiment of the present application, when the received data is obtained, a reference point can be determined through the received data to monitor the deformation of the monitoring target point. The reference point can determine the LOS diameter or the diameter energy greater than the diameter of the monitoring target point through the received data.

[0097] Step 340: Determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram.

[0098] Step 350: Determine the deformation amount of the monitoring target point according to the change of the phase difference.

[0099] Step 360: Generate a phase difference variation curve according to the phase differences at different moments.

[0100] In an embodiment of the present application, a phase difference change curve can be constructed with time and phase difference as two dimensions for the phase difference at different times. The phase difference change curve can reflect the changes of the reference point and the monitoring target point over time.

[0101] In some application embodiments, determining a reference point corresponding to the received data includes:

[0102] There is a line-of-sight propagation path corresponding to the received data, and the position point corresponding to the line-of-sight propagation path in the distance angle diagram is used as a reference point;

[0103] There is no line-of-sight propagation path corresponding to the received data, and the first path is selected as the reference point in the distance angle diagram, wherein the path energy value of the first path is greater than the monitoring target point.

[0104] The line-of-sight propagation path may be a channel model in which there is no obstruction between the two base stations, and the line of sight between the two base stations is not blocked.

[0105] In an embodiment of the present application, a distance angle diagram can be determined by receiving data, and the existence of a line-of-sight propagation path can be determined by the distance angle diagram. The position point of the line-of-sight propagation path corresponding to the distance angle diagram, that is, the position point where another base station is located in the line of sight, can be used as a reference point. If there is no line-of-sight propagation path, a first path can be selected in the distance angle diagram. The path energy value of the first path can be greater than the path energy value of the monitoring target point, and the first path can be used as a reference point.

[0106] FIG5 is an example diagram of a method for monitoring micro-deformation provided in an embodiment of the present application. The method for monitoring micro-deformation provided in an embodiment of the present application can be applied to the scenario shown in FIG1 . A transmitting base station sends a beam carrying data into space. The beam is reflected by the monitored target and reaches a receiving base station. The receiving base station can monitor the micro-deformation of the monitored target based on the received data and the physical parameters of the wireless communication network. Specifically, the monitoring method includes the following steps:

[0107] (1) Process the received data and calculate the channel impulse response.

[0108] In the embodiment of the present application, the received data R can be represented by a subcarrier N sub * Number of antennas N ant Matrix, where 1≤i sub ≤N sub , 1≤i ant ≤N ant The signal for transmitting and receiving data R is a 5G OFDM symbol, the subcarrier spacing is Δscs, the interval between two adjacent received symbols is T, and the i-th symbol is represented by i sym .

[0109] Since wireless communication networks suffer from delay and phase drift, this factor can be eliminated through delay-phase calibration. Starting from the second OFDM symbol, the data carried by all OFDM symbols is aligned with the data carried by the first OFDM symbol. The specific process of delay and phase calibration includes:

[0110] 1) Delay calibration:

[0111] i-th sym The delay difference between the first OFDM symbol carrying data and the first OFDM symbol carrying data is:

[0112] in, represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, Indicates the i-th symbol of the first sub +M subcarriers corresponding to the received data, Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, and conj() represents a function for determining the conjugate value.

[0113] In the embodiments of this application, Can be of length n ant Each element in the vector corresponds to the delay difference of each receiving antenna. The calculated delay difference can be used to compensate for the received data:

[0114] Among them, R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i represents the received data corresponding to the i-th symbol before delay compensation, Indicates the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, and j represents a complex number.

[0115] 2) Phase calibration:

[0116] i-th sym The phase difference between the first OFDM symbol carrying data and the first OFDM symbol carrying data is: Δφ i =angle(R i *conj(R i,τ ))

[0117] Among them, Δφ represents the phase difference corresponding to the i-th symbol, Ri Indicates the received data corresponding to the i-th symbol before delay compensation, R i,τ represents the received data corresponding to the i-th symbol after delay difference compensation, and conj() represents a function for determining the conjugate value.

[0118] Compensate the received data based on the phase difference: R i,τ,Δφ =R i,τ *exp(j*Δφ i )

[0119] Among them, R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol without phase calibration, exp() represents an exponential function with e as the base, and j represents a complex number.

[0120] In the embodiment of the present application, after the delay-phase calibration is performed on the received data, the received data can be subjected to IFFT to calculate the channel impulse response: h=ifft(R τ,Δφ ,N)

[0121] Among them, N is the number of points for IFFT, R τ,Δφ Indicates the received data after delay-phase calibration.

[0122] Specifically, one dimension of the channel impulse response is the sampling time Ts. The channel impulse response can represent the capabilities of signal paths at different distances. When a line-of-sight (LOS) path exists, the path with the highest capability in the channel impulse response is the LOS path, and information at a distance greater than the LOS path is the target's reflection path.

[0123] (2) Process the channel impulse response to generate a range-angle (RA) map.

[0124] In the embodiment of the present application, the channel impulse response includes signal paths at different distances, which is essentially the superposition of signals at different angles. In order to determine the path of the monitoring target, it is also necessary to separate the channel impulse response in the angle dimension. The method used is not limited to the peak search method, super-resolution algorithm and other technologies to extract the monitoring targets at different angles at a certain distance. Take the peak search method to separate the channel impulse response as an example:

[0125] 1) Constructing steering vectors of different angles, including but not limited to pitch angles and azimuth angles.

[0126] For any azimuth The pitch angle θ is used to construct the horizontal and vertical steering vectors respectively, and then the Kronecker product of the constructed steering vectors is calculated: a=kron(a H ,a V )

[0127] in, represents the azimuth angle, θ represents the elevation angle, j represents a complex number, Nz represents the number of antennas in the vertical direction of the receiving base station, Ny represents the number of antennas in the horizontal direction, T represents vector transpose, and λ represents the carrier wavelength.

[0128] 2) Calculate the signal strength of the received data at different angles at all distances:

[0129] For each steering vector, multiply it by the corresponding point of each antenna data of h at all distances, and obtain the modulus value. The range angle map (RA map) is obtained through the modulus value: RA = |h*a|

[0130] In an embodiment of the present application, the RA diagram is a multi-dimensional matrix of distance-elevation-azimuth. Furthermore, in actual testing, the horizontal or vertical angle can be fixed to a constant value to obtain a two-dimensional plane of distance azimuth or distance elevation. The signal energy (i.e., RA value) at different positions in the two-dimensional plane can be represented by different colors.

[0131] (3) Determine the coordinate positions of the target point and the reference point on the RA diagram.

[0132] In the embodiments of the present application, a reference point can be selected. The wireless communication system itself is affected by phase drift and factors such as temperature, humidity, and wind speed, and the system state changes in real time. A reference point can be determined first, and the target data result is subtracted from the reference point result to eliminate the inherent influence of the system. The method for selecting the reference point includes:

[0133] a) In the case where there is a LOS path between the transmitting base station and the receiving base station, the LOS path in the channel impulse response can be used as the reference point, and the coordinate value of the LOS path in the RA diagram can be used as the coordinate position P of the reference point. ref ;

[0134] b) If there is no LOS path between the transmitting base station and the receiving base station, a path in an environment with strong channel impulse response (such as a building) can be used as a reference point. The energy value of this path is at least not less than the target to be measured. The coordinate value of this reference point in the RA diagram can be used as the coordinate position P of the reference point. ref ;

[0135] c) For the case where there is neither a LOS path nor a path with strong energy between the transmitting base station and the receiving base station, a reference point can be constructed by adding a strong reflector in the environment as the reference point, and the coordinate value of the strong reflector in the RA diagram is used as the coordinate position P of the reference point. ref.

[0136] In the embodiment of the present application, the RA map may include multiple targets, which may include pedestrians, vehicles, buildings, etc. Taking the corner reflector as the target to be measured to test its real displacement as an example, the corner reflector can return the irradiated signal along the original path, and its radar cross section (RCS) is about 45m 2 .

[0137] a) Determine the position P of the corner reflector in the RA diagram by its distance and elevation and azimuth angles relative to the base station tar .

[0138] b) The corner reflector is not placed in the environment, and RA diagram A is generated by collecting data. The corner reflector is placed in the environment, and RA diagram B is generated by collecting data. Comparing the two RA diagrams, the target added to diagram B relative to diagram A can be determined to be the corner reflector. The coordinate position of the corner reflector in diagram B is P tar ′, through P tar ′=P tar , it can be seen that the position coordinates determined by the RA diagram are accurate.

[0139] (4) Calculate the phase difference between the target and the reference point before and after the target moves slightly.

[0140] In the embodiment of the present application, the phase before the target moves can be determined in the RA diagram as The phase of the reference point is After the target moves, the phase after the target moves is determined in the RA diagram as The phase of the reference point is Phase difference between the target and the reference point before the target undergoes slight deformation and movement The phase difference between the target and the reference point after the target moves slightly

[0141] (5) Calculate the change in the phase difference between the target and the reference point before and after the target moves slightly, and calculate the deformation amount based on the change.

[0142] In the embodiment of the present application, the change in the phase difference between the target and the reference point before and after the target moves slightly is: The carrier wavelength of the beam transmitting and receiving data is λ, and the phase change corresponding to one wavelength is 2π. Considering the round-trip signal, the one-way distance moved by the target can be expressed as:

[0143] The embodiment of the present application realizes the monitoring of tiny deformations through the 5G communication network, which can fully utilize the wide coverage advantage of the communication network to perform overall monitoring of the monitoring target, with millimeter-level monitoring accuracy, without the need for additional monitoring equipment, and reducing the monitoring cost of tiny deformations.

[0144] Furthermore, the method provided in the embodiment of the present application is verified by taking the slight deformation perception test with the reference point as the LOS path as an example. The test environment includes: the base station is a 4.9GHz 5G active antenna unit (Active Antenna Unit, AAU), 64 antennas, a bandwidth of 100MHz, a sampling rate of 122.88MHz, the base station is erected on the roof of a building, and a corner reflector is placed on the ground 139m away from the base station as a monitoring target. The pitch angle of the corner reflector relative to the base station is about -8 degrees, and the azimuth angle is about 1 degree. The corner reflector is placed on a guide rail, and the displacement of the corner reflector can be precisely controlled by the guide rail. There is a LOS path between the transceiver base station that does not pass through the target reflection.

[0145] Testing process:

[0146] (1) Without placing the corner reflector, collect data for a period of time, which is recorded as data 1.

[0147] (2) Place a corner reflector and collect data for a period of time. The guide rail controls the corner reflector to move 1 cm instantaneously, and collect data for another period of time, which is recorded as data 2.

[0148] Data processing process: process data one and data two separately.

[0149] (1) Perform delay phase calibration on the received data. Starting from the data of the second OFDM symbol, all OFDM symbols are aligned with the data of the first OFDM symbol. In the following formula, M can be 100, SCS = 30kHz = 30000;

[0150] in, represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, Indicates the i-th symbol of the first sub +M subcarriers corresponding to the received data, Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, and conj() represents a function for determining the conjugate value.

[0151] The calculated delay difference is used to compensate the original data to obtain the received data after delay compensation:

[0152] Among them, R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i represents the received data corresponding to the i-th symbol before delay compensation, Indicates the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, and j represents a complex number.

[0153] Calculate the phase difference:

[0154] The original data and the delay-compensated data are conjugate-multiplied and the angle is calculated to determine the phase difference: Δφ i =angloe(R i *conj(R i,τ ))

[0155] Use the calculated phase difference to compensate for the delay-compensated data to obtain the delay-phase aligned data: R i,τ,Δφ =R i,τ *exp(j*Δφ i )

[0156] Among them, R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol without phase calibration, exp() represents an exponential function with e as the base, and j represents a complex number.

[0157] A 4096-point IFFT is performed on the delay-phase aligned data to generate a channel impulse response. The horizontal axis of the channel impulse response, Ts = 1 / 122.88e, represents the sampling time. As shown in Figure 6, there are two strong signal paths in the environment: Ts = 1645, the LOS path between the transceiver base station, and Ts = 1701, the strong reflection path from the building.

[0158] (2) Set the direction angle The pitch angle θ = -8° constructs the horizontal and vertical steering vectors respectively, and then performs the Kronecker product: a=kron(a H ,a V )

[0159] in, represents the azimuth angle, θ represents the elevation angle, j represents a complex number, Nz represents the number of antennas in the vertical direction of the receiving base station, Ny represents the number of antennas in the horizontal direction, T represents vector transpose, and λ represents the carrier wavelength.

[0160] For each steering vector, perform a dot multiplication with the corresponding antenna data in h and take the modulus value. The resulting azimuth-filtered data can be represented by a range-angle graph, where the relationship between the range dimension and the Ts dimension in the channel impulse response is as follows: Where d represents the distance, T s represents the sampling time of the channel impulse response, T LOS represents the line-of-sight propagation time, d LOS Represents the line-of-sight distance, that is, the straight-line distance between the transmitting and receiving base stations. RA=|h*a|

[0161] (3) Find the reference point and the position of the target in the RA diagram. Use the LOS path as the reference point. The RA diagram of the data 1 obtained in step (1) is shown in Figure 7. In Figure 7, the position of LOS in the RA diagram is P ref (1.22,-4), which is used as a reference point.

[0162] Based on the distance of the corner reflector and its elevation and azimuth angles relative to the base station, its position can be determined in the RA diagram. Alternatively, it can be determined by the following method:

[0163] According to steps (1) and (2), the RA diagrams of the air-mined data and the data with the corner reflector are obtained in the range of 90m-165m, as shown in Figures 8 and 9, respectively. The RA diagram P of the corner reflector can be determined by comparing Figures 8 and 9. tar A target appears at the position (139, 1), which can be determined to be a corner reflector.

[0164] (4) Determine the phase difference between the corner reflector and the reference point for the RA diagram at different times of data 2, and record the phase of the corner reflector. Phase of the reference point The phase difference between the corner reflector and the reference point is:

[0165] The phase difference at different times can be A curve graph showing the phase difference changing with time is drawn, which is shown in FIG10 .

[0166] The 4.9GHz signal carrier wavelength is:

[0167] First deformation:

[0168] The average phase difference between the target and the reference point before deformation is After deformation The change in phase difference:

[0169] Deformation:

[0170] The actual measured deformation value is 1 cm, and the deformation monitoring error is 0.16 cm, or 1.6 mm.

[0171] Second deformation:

[0172] The average phase difference between the target and the reference point before deformation is After deformation The change in phase difference:

[0173] Deformation:

[0174] The actual measured deformation value is 1 cm, and the deformation monitoring error is 0.01 cm, or 0.1 mm.

[0175] In some other application embodiments, for example, two corner reflectors are set up in an environment where there is no LOS path, one corner reflector is used as a reference corner reflector, and the other corner reflector is used as a target corner reflector to perform micro-deformation monitoring.

[0176] Test environment: Suitable for environments with no LOS path or strong environmental path. Corner reflectors were placed on the ground at distances of 143m and 132m, respectively, as the target and reference corner reflectors. The base station was a 4.9GHz 5G AAU with 64 antennas, a bandwidth of 100MHz, and a sampling rate of 122.88MHz. The base station was erected on a building rooftop, and two corner reflectors were placed on the ground 139m away. The corner reflectors were placed on rails, which allowed for precise control of their displacement.

[0177] Specific testing process:

[0178] (1) Without placing the corner reflector, collect data for a period of time, which is recorded as data 1.

[0179] (2) Place corner reflectors at distances of 143m and 132m on the ground. Collect data for a period of time. Then, move the guide rail control corner reflectors 2mm instantaneously. Then, collect data for another period of time.

[0180] Process data 1 and data 2 separately:

[0181] (1) Perform delay phase calibration on the received data.

[0182] The specific method is the same as that in the above embodiment and will not be described again here.

[0183] A 4096-point IFFT is performed on the delay-phase-calibrated data to generate the channel impulse response. The horizontal axis, Ts = 1 / 122.88e6, represents the sampling time. There are no LOS paths in the environment, nor are there strong LOS paths in the environment.

[0184] (2) Set the direction angle At an elevation angle of θ = -8°, construct horizontal and vertical steering vectors. Then, perform a Kronecker product. For each steering vector, perform a dot-wise multiplication with the corresponding antenna data in h, and then take the modulus. This results in azimuth-filtered data, which can be represented as a range-angle graph.

[0185] (3) Determination of the location of the reference point:

[0186] The RA diagrams of data 1 and data 2 are shown in Figures 11 and 12, respectively. Comparing Figures 11 and 12, after placing the corner reflector, targets appear at positions (143, -2) and (132, -2), which are the positions P of the target corner reflector and the reference corner reflector, respectively. tar 、P ref .

[0187] (4) According to the RA diagram of data 2 at different times, the phase difference between the target corner reflector and the reference corner reflector can be determined. The phase of the target corner reflector can be recorded as Phase of the reference corner reflector The phase difference between the target corner reflector and the reference corner reflector is:

[0188] The phase difference at different times can be A curve graph showing the phase difference changing with time is drawn, which is shown in FIG13 .

[0189] First deformation:

[0190] The average phase difference between the target corner reflector and the reference corner reflector before deformation is After deformation The change in phase difference:

[0191] Deformation:

[0192] The actual measured deformation value is 2mm, and the deformation monitoring error is 0.7mm.

[0193] Second deformation:

[0194] The average phase difference between the target corner reflector and the reference corner reflector before deformation is After deformation The change in phase difference:

[0195] Deformation:

[0196] The actual measured deformation value is 2mm, and the deformation monitoring error is 0.1mm.

[0197] FIG14 is a schematic diagram of the structure of a micro-deformation monitoring device provided in an embodiment of the present application. The device can execute the micro-deformation monitoring method provided in any embodiment of the present application and has the corresponding functional modules and beneficial effects of the execution method. The device can be implemented by software and / or hardware. As shown in FIG14, the device provided in an embodiment of the present application specifically includes:

[0198] The data processing module 410 is configured to determine a channel impulse response based on received data from the wireless communication network.

[0199] The map generating module 420 is configured to generate a range angle map according to the channel impulse response.

[0200] The phase difference module 430 is configured to determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram.

[0201] The deformation variable module 440 is configured to determine the deformation variable of the monitoring target point according to the change of the phase difference.

[0202] Based on the above application embodiment, the data processing module 410 includes:

[0203] A data calibration unit is used to perform delay calibration and phase calibration on the received data.

[0204] A processing execution unit is used to generate a processing result of the inverse Fourier transform of the received data, and use the processing result as the channel impulse response.

[0205] Based on the above-mentioned application embodiment, the data calibration unit performs delay calibration on the received data, including:

[0206] Determining the time delay differences between other symbol data and the first symbol data in the received data;

[0207] Each of the time delay differences is compensated to the corresponding other symbol data to complete the time delay calibration of the received data.

[0208] Based on the above application embodiment, determining the time delay difference between other symbol data and the first symbol data in the received data includes:

[0209] Calling formula Determine the delay difference, wherein the represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, and Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, and conj() represents a function for determining the conjugate value.

[0210] Based on the above-mentioned embodiment of the application, each of the delay differences is compensated to the corresponding other symbol data to complete the delay calibration of the received data, including:

[0211] Calling formula Determine the received data after delay difference compensation, wherein the R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i represents the received data corresponding to the i-th symbol before delay difference compensation, represents the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, and j represents a complex number.

[0212] Based on the above application embodiment, the data calibration unit performs phase calibration on the received data, including:

[0213] Determining phase differences between other symbol data in the received data and the first symbol data;

[0214] The phase difference is compensated to the corresponding other symbol data to complete the phase calibration of the received data.

[0215] Based on the above application embodiment, determining the phase difference between each of the other symbols in the received data and the first symbol data includes:

[0216] Calling formula Δφ i =angle(R i *conj(R i,τ )) Determine the phase difference, wherein the Δφ i represents the phase difference corresponding to the i-th symbol, the R i represents the received data corresponding to the i-th symbol before delay difference compensation, the R i,τ represents the received data corresponding to the i-th symbol after delay difference compensation, and conj() represents a function for determining the conjugate value.

[0217] Based on the above-mentioned embodiment of the application, compensating the phase difference to the corresponding other symbol data to complete the phase calibration of the received data includes:

[0218] Call formula R i,τ,Δφ =R i,τ *exp(j*Δφ i ) determines the received data after phase calibration, wherein R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol that has not been phase-calibrated, exp() represents an exponential function with e as the base, and j represents a complex number.

[0219] Based on the above application embodiment, the graph generation module 420 includes:

[0220] The steering vector unit is configured to construct horizontal and vertical steering vectors for any direction angle and any pitch angle, and determine the Kronecker product of the horizontal steering vector and the vertical steering vector.

[0221] The map generating unit is configured to determine a modulus of a product of the Kronecker product and the channel impulse response, and construct the distance angle map based on the modulus.

[0222] Based on the above-mentioned application embodiment, the correspondence between the distance and the sampling time of the channel impulse response in the distance angle graph in the graph generation module 420 satisfies the following formula:

[0223] Wherein, d represents the distance, T s represents the sampling time of the channel impulse response, the T LOS represents the line-of-sight propagation time, and d LOS Indicates the line-of-sight propagation distance.

[0224] Based on the above application embodiment, the phase difference module 430 includes:

[0225] The position information unit is used to determine the first position information of the monitoring target point in the distance-angle graph according to the distance information and angle information of the monitoring target point.

[0226] A phase difference unit is used to determine the second position information of the reference point in the distance angle diagram, and determine the phase difference according to the first position information and the second position information.

[0227] Based on the above application embodiment, it also includes: a reference point module, used to determine the reference point corresponding to the received data.

[0228] Based on the above application embodiment, the reference point module is specifically configured to do at least one of the following: if there is a line-of-sight propagation path corresponding to the received data, the position point of the line-of-sight propagation path corresponding to the distance angle graph is used as the reference point;

[0229] There is no line-of-sight propagation path corresponding to the received data, and a first path is selected as the reference point in the distance angle diagram, wherein the path energy value of the first path is greater than the monitoring target point.

[0230] Based on the above-mentioned embodiment, the deformation variable module 440 includes:

[0231] The phase difference changing unit is used to determine the change amount of the phase difference between the two most recent moments.

[0232] The deformation amount determining unit is used to determine the one-way distance change amount corresponding to the change amount as the deformation amount according to the carrier wavelength and carrier phase change corresponding to the received data.

[0233] Based on the above application embodiment, the shape variable determination unit is specifically used to: call the formula Determine the deformation amount, wherein the Δd represents the deformation amount, the represents the variation, λ represents the carrier wavelength, and the carrier phase variation is 2π.

[0234] On the basis of the above application embodiment, it further includes: a change curve module, which is used to generate a phase difference change curve according to the phase difference at different moments.

[0235] Figure 15 is a structural diagram of an electronic device provided in an embodiment of the present application, which electronic device includes a processor 10, a memory 11, an input device 12 and an output device 13; the number of processors 10 in the electronic device can be one or more, and Figure 15 takes one processor 10 as an example; the processor 10, memory 11, input device 12 and output device 13 in the electronic device can be connected via a bus or other means, and Figure 15 takes connection via a bus as an example.

[0236] The memory 11, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, and modules, such as the modules corresponding to the micro-deformation monitoring device in the embodiment of the present application (data processing module 410, image generation module 420, phase difference module 430, and deformation amount module 440). The processor 10 executes the software programs, instructions, and modules stored in the memory 11 to execute various functional applications and data processing of the electronic device, thereby implementing the above-mentioned method.

[0237] The memory 11 may primarily include a program storage area and a data storage area. The program storage area may store an operating system and at least one application required for a function; the data storage area may store data generated based on the use of the electronic device. Furthermore, the memory 11 may include high-speed random access memory and may also include non-volatile memory, such as at least one disk storage device, flash memory device, or other non-volatile solid-state memory device. In some instances, the memory 11 may further include memory remotely located relative to the processor 10, and such remote memory may be connected to the electronic device via a network. Examples of such networks include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0238] The input device 12 may be used to receive input digital or character information and generate key signal input related to user settings and function control of the electronic device. The output device 13 may include a display device such as a display screen.

[0239] The present application also provides a storage medium containing computer-executable instructions. When the computer-executable instructions are executed by a computer processor, the computer-executable instructions are used to perform a method for monitoring micro-deformation. The method includes:

[0240] determining a channel impulse response based on received data from a wireless communication network;

[0241] generating a range angle graph according to the channel impulse response;

[0242] Determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram;

[0243] The deformation amount of the monitoring target point is determined according to the change of the phase difference.

[0244] Through the above description of the implementation methods, those skilled in the art can clearly understand that the present application can be implemented with the help of software and necessary general-purpose hardware, and of course can also be implemented by hardware. Based on this understanding, the technical solution of the present application can be embodied in the form of a software product, which can be stored in a computer-readable storage medium, such as a computer's floppy disk, read-only memory (ROM), random access memory (RAM), flash memory (FLASH), hard disk or optical disk, etc., including a number of instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute the methods described in each embodiment of the present application.

[0245] It is worth noting that in the embodiments of the above-mentioned device, the various units and modules included are only divided according to functional logic, but are not limited to the above-mentioned division, as long as the corresponding functions can be achieved; in addition, the specific names of the functional units are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application.

[0246] Those skilled in the art will appreciate that all or some of the steps, devices, and functional modules / units in the methods disclosed above may be implemented as software, firmware, hardware, or appropriate combinations thereof.

[0247] In a hardware implementation, the division between the functional modules / units mentioned in the above description does not necessarily correspond to the division of physical components; for example, a physical component may have multiple functions, or a function or step may be performed by several physical components in cooperation. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, a digital signal processor, or a microprocessor, or implemented as hardware, or implemented as an integrated circuit, such as an application-specific integrated circuit. The corresponding software may be distributed on a computer-readable medium, which may include computer storage media (or non-transitory media) and communication media (or temporary media). As is well known to those skilled in the art, the term computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program modules or other data). Computer storage media include, but are not limited to, RAM, ROM, Electrically Erasable Programmable Read-Only Memory (EEPROM), flash memory or other memory technology, Compact Disc Read Only Memory (CD-ROM), Digital Versatile Disk (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store the desired information and can be accessed by a computer. In addition, it is well known to those skilled in the art that communication media generally contain computer-readable instructions, data structures, program modules, or other data in a modulated data signal such as a carrier wave or other transport mechanism, and may include any information delivery media.

Claims

1. A method for monitoring micro-deformation, comprising: determining a channel impulse response based on received data from a wireless communication network; generating a range angle graph according to the channel impulse response; Determine the phase difference between the monitoring target point and the reference point according to the distance angle diagram; The deformation amount of the monitoring target point is determined according to the change of the phase difference.

2. The method according to claim 1, wherein: Determining the channel impulse response according to the received data of the wireless communication network includes: Performing delay calibration and phase calibration on the received data; A processing result of an inverse Fourier transform of the received data is generated, and the processing result is used as the channel impulse response.

3. The method according to claim 2, wherein: Performing delay calibration on the received data, comprising: Determining the time delay differences between other symbol data and the first symbol data in the received data; Each of the time delay differences is compensated to the corresponding other symbol data to complete the time delay calibration of the received data.

4. The method according to claim 2, wherein: Performing phase calibration on the received data, comprising: Determining phase differences between other symbol data in the received data and the first symbol data; The phase difference is compensated to the corresponding other symbol data to complete the phase calibration of the received data.

5. The method according to claim 3, wherein: The determining of the time delay differences between the other symbol data and the first symbol data in the received data includes: Calling formula Determine the delay difference, wherein the represents the delay difference of the i-th symbol data in the received data, M represents a constant value set based on experience, and Indicates the i-th symbol of the first sub +M subcarriers corresponding to the received data, the Represents all received data corresponding to the i-th symbol, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, angle[] represents a function for determining the phase, and conj() represents a function for determining the conjugate value.

6. The method according to claim 3, wherein: The compensating each of the delay differences to the corresponding other symbol data to complete the delay calibration of the received data includes: Calling formula Determine the received data after delay difference compensation, wherein the R i,τ represents the received data corresponding to the i-th symbol after delay compensation, R i Indicates the time before delay compensation The received data corresponding to the i-th symbol, represents the delay difference of the i-th symbol data in the received data, i sub represents the number of subcarriers on the i-th symbol, Δscs represents the subcarrier spacing, exp() represents an exponential function with e as the base, and j represents a complex number.

7. The method according to claim 4, wherein: The determining of the phase differences between the first symbol data and the other symbol data in the received data includes: Calling formula Δφ i =angle(R i *conj(R i,τ )) Determine the phase difference, wherein the Δφ i represents the phase difference corresponding to the i-th symbol, the R i represents the received data corresponding to the i-th symbol before delay compensation, the R i,τ represents the received data corresponding to the i-th symbol after delay difference compensation, and conj() represents a function for determining the conjugate value.

8. The method according to claim 4, wherein: The compensating the phase difference to the corresponding other symbol data to complete the phase calibration of the received data includes: Call formula R i,τ,Δφ =R i,τ *exp(j*Δφ i ) determines the received data after phase calibration, wherein R i,τ,Δφ represents the received data of the i-th symbol after phase alignment, R i,τ represents the received data corresponding to the i-th symbol that has not been phase-calibrated, exp() represents an exponential function with e as the base, and j represents a complex number.

9. The method according to claim 1, wherein: Generating a distance angle diagram according to the channel impulse response includes: constructing a horizontal steering vector and a vertical steering vector for any direction angle and any pitch angle, and determining a Kronecker product of the horizontal steering vector and the vertical steering vector; A modulus of a product of the Kronecker product and the channel impulse response is determined, and the range-angle map is constructed based on the modulus.

10. The method according to claim 9, wherein: The corresponding relationship between the distance in the distance angle diagram and the sampling time of the channel impulse response satisfies the following formula: Wherein, d represents the distance, T s represents the sampling time of the channel impulse response, the T LOS represents the line-of-sight propagation time, and d LOS Indicates the line-of-sight propagation distance.

11. The method according to claim 1, wherein: Determining the phase difference between the monitoring target point and the reference point according to the distance angle diagram includes: Determine first position information of the monitoring target point in the distance-angle graph according to the distance information and angle information of the monitoring target point; Determine second position information of the reference point in the distance angle graph, and determine the phase difference according to the first position information and the second position information.

12. The method according to claim 1, further comprising: A reference point corresponding to the received data is determined.

13. The method according to claim 12, wherein: Determining the reference point corresponding to the received data includes at least one of the following: There is a line-of-sight propagation path corresponding to the received data, and a position point corresponding to the distance angle graph of the line-of-sight propagation path is used as the reference point; There is no line-of-sight propagation path corresponding to the received data, and a first path is selected as the reference point in the distance angle diagram, wherein the path energy value of the first path is greater than the path energy value of the monitoring target point.

14. The method according to claim 1, wherein: The determining the deformation amount of the monitoring target point according to the change of the phase difference includes: Determining a change in the phase difference between the two most recent moments; According to the carrier wavelength and carrier phase changes corresponding to the received data, a one-way distance change corresponding to the change is determined as the deformation amount.

15. The method according to claim 14, wherein: The determining, according to the carrier wavelength and carrier phase changes corresponding to the received data, a one-way distance change corresponding to the change as the deformation amount, includes: Calling formula λ determines the deformation amount, wherein the Δd represents the deformation amount, the represents the variation, λ represents the carrier wavelength, and the carrier phase variation is 2π.

16. The method according to claim 1, further comprising: A phase difference variation curve is generated according to the phase difference at different moments.

17. An electronic device comprising: at least one processor; a memory configured to store at least one program; When the at least one program is executed by the at least one processor, the at least one processor implements the micro-deformation monitoring method according to any one of claims 1 to 16.

18. A computer-readable storage medium, wherein: The computer-readable storage medium stores at least one program, and the at least one program is executed by at least one processor to implement the micro-deformation monitoring method according to any one of claims 1 to 16.

Images