patents.google.com

CN110987150B - A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar - Google Patents

  • ️Fri Jul 23 2021
A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar Download PDF

Info

Publication number
CN110987150B
CN110987150B CN201911295312.9A CN201911295312A CN110987150B CN 110987150 B CN110987150 B CN 110987150B CN 201911295312 A CN201911295312 A CN 201911295312A CN 110987150 B CN110987150 B CN 110987150B Authority
CN
China
Prior art keywords
signal
radar
micro
channel
baseband signal
Prior art date
2019-12-16
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911295312.9A
Other languages
Chinese (zh)
Other versions
CN110987150A (en
Inventor
熊玉勇
彭志科
李松旭
任泽生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiao Tong University
Original Assignee
Shanghai Jiao Tong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
2019-12-16
Filing date
2019-12-16
Publication date
2021-07-23
2019-12-16 Application filed by Shanghai Jiao Tong University filed Critical Shanghai Jiao Tong University
2019-12-16 Priority to CN201911295312.9A priority Critical patent/CN110987150B/en
2020-04-10 Publication of CN110987150A publication Critical patent/CN110987150A/en
2021-07-23 Application granted granted Critical
2021-07-23 Publication of CN110987150B publication Critical patent/CN110987150B/en
Status Active legal-status Critical Current
2039-12-16 Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于单频连续波雷达的超微幅振动测量方法,涉及非接触式振动测量技术领域,包括如下步骤:雷达天线正对测量目标放置,然后控制前后移动雷达前端,同时同步采集雷达基带信号,包括I通道输出信号I0(t)和Q通道输出信号Q0(t);利用采集的基带信号I0(t)和Q0(t)进行圆拟合与圆心估计,得到初始化I/Q通道基带信号直流偏移量估计值;雷达的前端保持静止,同步采集微波雷达的I/Q通道基带信号,并估计I/Q通道基带信号随时间变化的漂移量;从I/Q通道基带信号中减去漂移量并进行低通滤波处理,采用反正切解调算法实现目标微振动位移时域信息提取。该方法能够解决超微幅振动测量中基带信号的直流偏移量估计问题,显著降低测量误差,减少计算量,提高测量精度。

Figure 201911295312

The invention discloses an ultra-micro-amplitude vibration measurement method based on a single-frequency continuous wave radar, and relates to the technical field of non-contact vibration measurement. Collect radar baseband signals, including I channel output signal I 0 (t) and Q channel output signal Q 0 (t); use the collected baseband signals I 0 (t) and Q 0 (t) to perform circle fitting and circle center estimation, The estimated value of DC offset of the baseband signal of the initialized I/Q channel is obtained; the front end of the radar is kept stationary, the baseband signal of the I/Q channel of the microwave radar is collected synchronously, and the drift of the baseband signal of the I/Q channel with time is estimated; The drift amount is subtracted from the baseband signal of the /Q channel and low-pass filtering is performed, and the arctangent demodulation algorithm is used to extract the time domain information of the target micro-vibration displacement. The method can solve the problem of estimating the DC offset of the baseband signal in the ultra-micro-amplitude vibration measurement, significantly reduce the measurement error, reduce the calculation amount, and improve the measurement accuracy.

Figure 201911295312

Description

Ultramicro amplitude vibration measurement method and system based on single-frequency continuous wave radar

Technical Field

The invention relates to the technical field of non-contact vibration measurement, in particular to a single-frequency continuous wave radar-based ultramicro amplitude vibration measurement method and system.

Background

Vibration measurement is an important way for realizing structural mechanical property test, health monitoring and fault diagnosis and other micro-motion characteristic extraction. There are two types of measurement, contact and noncontact vibration, which have significant advantages in terms of sensor installation simplicity, lightweight structural measurement, and the like. The vibration measurement method based on the single-frequency continuous wave microwave radar is that a continuous wave microwave signal with a certain frequency is transmitted, an electromagnetic echo signal reflected by a target is received, and vibration information extraction is realized through Doppler nonlinear phase demodulation. For single-frequency continuous wave microwave radar vibration measurement based on zero intermediate frequency orthogonal output, an arc tangent demodulation method is often adopted to extract vibration displacement time domain information. The method needs accurate direct current offset compensation, and in the prior art, the direct current offset is estimated by performing arc fitting and circle center estimation by utilizing orthogonal output baseband signals. However, for micro-amplitude, especially ultra-micro-amplitude (1/360 with a vibration amplitude less than the carrier wavelength) vibration measurement, because the phase change of the baseband signal is very small and the arc occupancy is very small, the direct current offset of the baseband signal cannot be accurately estimated by using the arc fitting and circle center estimation method, and the estimation error is very large, so that the ultra-micro-amplitude vibration measurement cannot be performed. In actual measurement, due to inevitable noise interference, the difficulty of accurately estimating the direct current offset of the radar baseband signal is further increased.

In addition, due to hardware performance imperfections, single frequency continuous wave radars in a zero if configuration can produce amplitude drift in the baseband signal due to frequency drift of the vco and output drift of the mixer. The drift is generally small and can be ignored for large amplitude vibration measurement, but for ultra-small amplitude vibration measurement, the amplitude drift of the baseband signal can have a large influence on the extraction of the small amplitude vibration information.

The invention aims to solve the problem of ultra-micro amplitude vibration measurement based on a single-frequency continuous wave microwave radar, solve the problems of radar baseband signal direct current offset estimation and baseband voltage signal amplitude drift in the ultra-micro amplitude vibration measurement, and accurately extract vibration displacement time domain information.

Therefore, those skilled in the art are dedicated to develop a method and a system for measuring ultra-micro amplitude vibration based on a single-frequency continuous wave radar, which can solve the problem of estimating the direct current offset of radar baseband I/Q signals in ultra-micro amplitude vibration measurement, significantly reduce the measurement error of an ultra-micro amplitude vibration displacement time sequence, and reduce the calculation amount. The problem of baseband signal drift caused by non-ideality of hardware performance is eliminated, and the precision of the ultra-micro amplitude vibration measurement is further improved. The method and the system solve the problem that the single-frequency continuous wave radar ultramicro amplitude vibration measurement is difficult to apply.

Disclosure of Invention

In view of the above defects in the prior art, the technical problem to be solved by the invention is to solve the problem of ultra-micro amplitude vibration measurement based on a single-frequency continuous wave microwave radar, solve the problems of radar baseband signal direct current offset estimation and baseband voltage signal amplitude drift in the ultra-micro amplitude vibration measurement, and accurately extract vibration displacement time domain information.

In order to achieve the aim, the invention provides an ultramicro amplitude vibration measurement method based on a single-frequency continuous wave radar, which comprises the following steps:

(1) the radar antenna is placed right opposite to the measurement target, then the front end of the radar is controlled to move back and forth, and simultaneously, radar baseband signals are synchronously acquired, including an I channel output signal I0(t) and Q channel output signal Q0(t);

(2) Using the acquired baseband signal0(t) and said Q0(t) carrying out circle fitting and circle center estimation to obtain an initialized I/Q channel baseband signal direct current offset estimation value;

(3) the front end of the radar keeps static, the I/Q channel baseband signal of the microwave radar is synchronously collected, and the drift amount of the I/Q channel baseband signal changing along with time is estimated;

(4) and subtracting the drift amount from the I/Q channel baseband signal, performing low-pass filtering, and extracting target micro-vibration displacement time domain information by adopting an arc tangent demodulation algorithm.

Further, the moving method of the radar front section in the step (1) is as follows:

firstly, controlling to move forwards at a constant speed by lambda/2, and then moving backwards at a constant speed by lambda/2 to return to the original position. Wherein λ is a transmission carrier wavelength of the radar.

Further onInitializing a DC offset DC of the I/Q channel baseband signal in the step (2)I0And DCQ0The estimation method comprises the following steps:

setting the baseband signal I of acquisition0(t) and Q0(t) discrete sequences are each I0[n]And Q0[n]The coordinate of the center of the fitting circle is (x)0,y0) And the radius is r, taking a matrix A, and the matrix m and the matrix B are respectively as follows:

Figure BDA0002320356370000021

obtaining x by parameter optimization estimation0And y0The value of (c).

Further, the method for optimizing the estimation comprises the following steps: min | | Am-B | non-conducting phosphorl2Or min | | Am-B | | non-woven phosphorl1Wherein min (·) represents a minimum operation, | |. luminance |)l2Represents 2-norm operation, | |l1Represents a 1-norm operation to find x0And y0After estimating the value of (2), take DCI0=x0,DCQ0=y0

Further, the drift amounts I _ drift (t) and Q _ drift (t) of the I/Q baseband signal varying with time in the step (3) are obtained by:

I_drift(t)=I_trend(t)-mean[I_trend(t)]

Q_drift(t)=Q_trend(t)-mean[Q_trend(t)]

wherein mean [. cndot. ] is an averaging operation, the calculation method of the trend components I _ trend (t) and Q _ trend (t) of the I/Q baseband signal is local weighted linear regression smoothing, and the number of smoothing points is preferably an integer equal to the sampling frequency.

Further, the cut-off frequency of the low-pass filtering in the step (4) is 5 fobjectWherein f isobjectIs the vibration frequency of the measured object.

Further, the calculation method of the time series of the micro-vibration displacement in the step (4) is that

Figure BDA0002320356370000031

In the formula

Figure BDA0002320356370000032

A phase time sequence of the I/Q baseband signal calculated for the arctan demodulation, wherein,

Figure BDA0002320356370000033

the calculation method is

Figure BDA0002320356370000034

Furthermore, the device comprises a radar front end, a translation sliding table, an acquisition and processing module, a controller and a vibration analysis and display module, wherein the radar front end is fixedly arranged on the translation sliding table.

Further, the radar front end comprises a microwave signal source, a power divider, a power amplifier, a quadrature phase shifter, a low noise amplifier, a frequency mixer, a low pass filter, a signal conditioning circuit, a transmitting antenna and a receiving antenna, wherein the microwave signal source is connected with the power divider, the power divider is connected with the power amplifier and the quadrature phase shifter, the power amplifier is connected with the transmitting antenna, the receiving antenna is connected with the low noise amplifier, two output ends of the low noise amplifier and the quadrature phase shifter are respectively connected with the frequency mixer, two output ends of the frequency mixer are respectively connected with the low pass filter, an output end of the low pass filter is connected with the signal conditioning circuit, a signal of the microwave signal source is divided into two paths by the power divider, and one path is connected with the transmitting antenna by the power amplifier, one path is connected with the frequency mixer through the quadrature phase shifter.

Furthermore, the radar front end is used for generating and transmitting a single-frequency continuous wave microwave signal, receiving an electromagnetic echo scattered by a target, and obtaining a baseband I/Q signal through amplification, frequency mixing, filtering and other processing;

the acquisition and processing module comprises a signal acquisition unit and a signal processing unit, wherein the signal acquisition unit is used for synchronously acquiring baseband I/Q signals output by the front end of the radar, and the signal processing unit is used for realizing the estimation of an initialized direct current offset value, the estimation of a baseband I/Q signal drift amount and the calculation of a vibration displacement time sequence in the single-frequency continuous wave radar-based ultramicro amplitude vibration measurement method;

the controller is used for controlling the sequential execution of all steps in the single-frequency continuous wave radar-based ultramicro amplitude vibration measurement method, controlling the movable sliding table, setting sampling parameters, executing signal acquisition, performing estimation operation of initialized direct current offset, estimating I/Q baseband signal drift amount and calculating and processing a vibration displacement time sequence;

and the vibration analysis and display module is used for extracting and analyzing vibration characteristics such as vibration amplitude, frequency and time-frequency distribution of the measured ultramicro amplitude vibration time domain information according to requirements and displaying signal waveforms and related measurement and characteristic extraction results.

The single-frequency continuous wave radar-based ultramicro amplitude vibration measurement method and system provided by the invention can solve the problem of direct current offset estimation of radar baseband I/Q signals in ultramicro amplitude vibration measurement, obviously reduce the measurement error of an ultramicro amplitude vibration displacement time sequence and reduce the calculated amount at the same time. The problem of baseband signal drift caused by non-ideality of hardware performance is eliminated, and the precision of the ultra-micro amplitude vibration measurement is further improved. The method and the system solve the problem that single-frequency continuous wave radar ultramicro amplitude vibration measurement is difficult to apply, and simultaneously provide a measuring method and a measuring system which are low in cost, high in precision and easy to operate.

The conception, the specific structure and the technical effects of the present invention will be further described with reference to the accompanying drawings to fully understand the objects, the features and the effects of the present invention.

Drawings

FIG. 1 is a flow chart of a method for measuring ultra-micro amplitude vibration based on a single-frequency continuous wave radar according to a preferred embodiment of the present invention;

FIG. 2 is a schematic diagram of the structure of the mobile radar front end of the present invention;

FIG. 3 is a diagram illustrating initial circle fitting and DC offset estimation results according to an embodiment of the present invention;

FIG. 4 is a diagram illustrating an example of the extraction of trend components of an I-channel baseband signal according to the present invention;

FIG. 5 is a diagram illustrating an example of the extraction of trend components of a Q-channel baseband signal according to the present invention;

FIG. 6 is a comparison of a prior art ultramicro amplitude vibration waveform measurement with a reference value in an embodiment of the present invention;

FIG. 7 shows the comparison of the measured ultramicro amplitude vibration waveform with a reference value according to an embodiment of the present invention;

FIG. 8 is a graph of prior art circle fitting and DC offset estimation results in an embodiment of the present invention;

FIG. 9 is a structural block diagram of an ultra-micro amplitude vibration measurement system based on a single-frequency continuous wave radar according to the present invention;

fig. 10 is a schematic structural diagram of a radar front end according to the present invention.

Detailed Description

The technical contents of the preferred embodiments of the present invention will be more clearly and easily understood by referring to the drawings attached to the specification. The present invention may be embodied in many different forms of embodiments and the scope of the invention is not limited to the embodiments set forth herein.

In the drawings, structurally identical elements are represented by like reference numerals, and structurally or functionally similar elements are represented by like reference numerals throughout the several views. The size and thickness of each component shown in the drawings are arbitrarily illustrated, and the present invention is not limited to the size and thickness of each component. The thickness of the components may be exaggerated where appropriate in the figures to improve clarity.

A method for measuring ultra-micro amplitude vibration based on a single-frequency continuous wave radar is shown in figure 1 and comprises the following steps:

step

1, placing a radar antenna over against a measurement target, then controlling to move the front end of a microwave radar back and forth, and synchronously acquiring radar baseband signals including an I channel output signal I0(t) and Q channel output signal Q0(t) of (d). In order to accurately estimate the direct current offset of the initialized I/Q channel baseband signal by using a circle fitting and circle center estimation algorithm, the phase of the baseband signal is changed approximately due to the movement of a radar antenna relative to a targetEqual to 2 pi, namely, the circular arc is distributed on the whole circumference in the circle fitting, and the front-back movement distance of the radar is preferably lambda/2. For example, as shown in fig. 2, a method of moving a microwave radar front end: firstly, controlling to move forwards at a constant speed by lambda/2, and then moving backwards at a constant speed by lambda/2 to return to the original position. Wherein, λ is the transmission carrier wavelength of the microwave radar.

Step

2, utilizing the collected baseband signal I0(t) and Q0And (t) carrying out circle fitting and circle center estimation on the signals to obtain an initialized I/Q channel baseband signal direct current offset estimation value.

DC offset DC of initialized I/Q channel baseband signalI0And DCQ0The estimation method comprises the following steps:

setting the baseband signal I of acquisition0(t) and Q0(t) discrete sequences are each I0[n]And Q0[n]The coordinate of the center of the fitting circle is (x)0,y0) And the radius is r, taking a matrix A, and the matrix m and the matrix B are respectively as follows:

Figure BDA0002320356370000051

obtaining x by parameter optimization estimation0And y0The optimized estimation method is min | | | Am-B | luminancel2Or min | | Am-B | | non-woven phosphorl1

In the formula, min (& lt|. & gt) represents the minimum value operation, | | | & lt| & gtLul2Represents 2-norm operation, | |l1Representing a 1-norm operation.

Finding x0And y0After estimating the value of (2), take DCI0=x0,DCQ0=y0

FIG. 3 shows an example of the estimation result of the measured initialized DC offset, in which the carrier frequency is 24.2GHz, the moving distance of the front end of the radar is 6.2mm, and the I/Q signal I is acquired for the baseband0(t) and Q0(t) performing circle fitting and circle center estimation to obtain initialized direct current offset values of the I/Q channels, wherein the initialized direct current offset values are respectively DCI0=4.2327,DCQ0=-0.8836。

And 3, keeping the front end of the radar still, carrying out ultramicro amplitude vibration measurement on the target to be measured, synchronously acquiring the I/Q signal of the radar baseband, and estimating the drift amount of the I/Q signal along with the change of time.

The method for obtaining the drift amounts I _ drift (t) and Q _ drift (t) of the IQ baseband signal along with the change of time comprises the following steps:

I_drift(t)=I_trend(t)-mean[I_trend(t)]

Q_drift(t)=Q_trend(t)-mean[Q_trend(t)]

in the formula, mean [. cndot. ] is the operation of taking the average value.

The trend components I _ trend (t) and Q _ trend (t) of the baseband IQ signal are preferably derived by locally weighted linear regression smoothing. Since the carrier and mixer output drift over time tends to exhibit a slow low frequency drift, the number of smoothing points for the trend component is preferably equal to an integer number of sampling frequencies. Fig. 4 and 5 show an example of the extraction result of the trend component of the measured I-channel and Q-channel baseband signals, respectively.

And 4, subtracting the drift amount from the baseband IQ signal, performing low-pass filtering, and extracting the target micro-vibration displacement time domain information by adopting an arc tangent demodulation algorithm.

The cut-off frequency of the low-pass filtering is determined according to the vibration frequency f of the measured targetobjectPreferably 5 x f, or a priori knowledge of the measured signal strengthobject. For example, the micro-vibration frequency of the object to be measured is 10Hz, and the cut-off frequency of the low-pass filtering is preferably set to 50 Hz.

Time series of vibratory displacements x [ n ]]The calculation method comprises the following steps:

Figure BDA0002320356370000061

in the formula

Figure BDA0002320356370000062

The phase time sequence of the calculated baseband signal is calculated for arctan demodulation.

Figure BDA0002320356370000063

The calculation method comprises the following steps:

Figure BDA0002320356370000064

in the measurement, when the phase time sequence of the baseband signal is calculated by adopting an arc tangent demodulation algorithm, the initialized direct current offset value of the baseband I/Q signal is directly used for demodulation. The calculation method solves the technical problem that in the ultramicro amplitude vibration measurement, the direct current offset cannot be estimated based on the collected baseband I/Q signals, so that the ultramicro amplitude vibration measurement cannot be carried out; on the other hand, the calculation of the time sequence of the ultra-micro amplitude vibration displacement directly adopts the initialized direct current offset for compensation calculation, and the direct current offset does not need to be estimated for multiple times based on a sliding window, so that the calculation amount of the algorithm is further reduced. Fig. 6 and fig. 7 show the measurement results of the time domain information of the ultra-micro amplitude vibration displacement measured based on the methods proposed by the prior art and the present invention, respectively. It can be seen that the prior art cannot realize ultra-micro amplitude vibration measurement, the measurement error is very large, the result of performing direct current offset estimation by using baseband I/Q signals is shown in fig. 8, and the error between circular arc fitting and circle center estimation is very large. Compared with the prior art, the method provided by the invention can be used for better extracting the ultramicro amplitude vibration information.

A single-frequency continuous wave radar-based ultramicro amplitude vibration measurement system is shown in a system structure block diagram in fig. 9, and comprises:

the radar front end is used for generating and transmitting single-frequency continuous wave microwaves, receiving electromagnetic echoes scattered by a target, adopting a direct down-conversion structure, and obtaining zero intermediate frequency baseband signals through amplification, frequency mixing, filtering and other processing. As shown in fig. 10, the radar front end includes a microwave signal source, a power divider, a power amplifier, a quadrature phase shifter, a low noise amplifier, a mixer, a low pass filter, a signal conditioning circuit, a transmitting antenna, and a receiving antenna, where the microwave signal source is connected to the power divider, the power divider is connected to the power amplifier and the quadrature phase shifter, the power amplifier is connected to the transmitting antenna, the receiving antenna is connected to the low noise amplifier, two output ends of the low noise amplifier and the quadrature phase shifter are respectively connected to the mixer, two output ends of the mixer are respectively connected to the low pass filter, and an output end of the low pass filter is connected to the signal conditioning circuit. The signal of the microwave signal source is divided into two paths by the power divider, one path is connected with the transmitting antenna by the power amplifier, and the other path is connected with the frequency mixer by the quadrature phase shifter.

And the translation sliding table is used for controlling the front end of the moving radar to move together with the sliding table according to a specific requirement, wherein the front end of the radar is fixedly arranged on the translation sliding table.

The acquisition and processing module comprises a signal acquisition unit and a signal processing unit, wherein the signal acquisition unit is used for synchronously acquiring baseband I/Q signals output by the front end of the radar. The signal processing unit is used for realizing the estimation of the initialized direct current offset value, the estimation of the baseband I/Q signal drift amount and the calculation of the vibration displacement time sequence in the ultramicro amplitude vibration measurement method based on the single-frequency continuous wave radar.

And the controller is used for controlling the movement, signal acquisition and processing operation of the translation sliding table. Preferably, the controller is configured to control to sequentially execute each step in the above method for measuring ultra-micro amplitude vibration based on the single-frequency continuous wave radar, and to control the moving of the sliding table, set sampling parameters, execute signal acquisition, perform estimation operation of an initialized dc offset, estimate a baseband I/Q signal drift amount, and calculate and process a vibration displacement time sequence.

And the vibration analysis and display module is used for extracting and analyzing vibration characteristics such as vibration amplitude, frequency and time-frequency distribution of the measured ultramicro amplitude vibration time domain information according to requirements and displaying signal waveforms and related measurement and characteristic extraction results.

The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (2)

1.一种基于单频连续波雷达的超微幅振动测量方法,其特征在于,包括如下步骤:1. an ultra-micro-amplitude vibration measuring method based on single frequency continuous wave radar, is characterized in that, comprises the steps: (1)雷达天线正对测量目标放置,然后控制前后移动雷达前端,同时同步采集雷达基带信号,包括I通道输出信号I0(t)和Q通道输出信号Q0(t);(1) The radar antenna is placed facing the measurement target, and then the radar front end is controlled to move forward and backward, and the radar baseband signal is collected synchronously, including the I channel output signal I 0 (t) and the Q channel output signal Q 0 (t); (2)利用采集的基带信号所述I0(t)和所述Q0(t)进行圆拟合与圆心估计,得到初始化的I/Q通道基带信号直流偏移量估计值;(2) using the I 0 (t) and the Q 0 (t) of the collected baseband signals to perform circle fitting and circle center estimation to obtain an initialized I/Q channel baseband signal DC offset estimate; (3)所述雷达前端保持静止,同步采集所述雷达的所述I/Q通道基带信号,并估计所述I/Q通道基带信号随时间变化的漂移量;(3) The radar front end remains stationary, synchronously collects the I/Q channel baseband signal of the radar, and estimates the drift amount of the I/Q channel baseband signal over time; (4)从所述I/Q通道基带信号中减去漂移量,并进行低通滤波处理,然后采用反正切解调算法实现目标微振动位移时域信息的提取;(4) subtract the drift amount from the baseband signal of the I/Q channel, and carry out low-pass filtering processing, and then adopt the arctangent demodulation algorithm to realize the extraction of target micro-vibration displacement time domain information; 所述步骤(1)中的所述雷达前段的移动方法为:The moving method of the front section of the radar in the step (1) is: 首先控制匀速向前移动λ/2,然后再匀速向后移动λ/2回到原始位置,其中,λ为所述雷达的发射载波波长;First, it is controlled to move forward by λ/2 at a constant speed, and then move backward by λ/2 at a constant speed to return to the original position, where λ is the wavelength of the transmitting carrier of the radar; 所述步骤(2)中的初始化所述I/Q通道基带信号的直流偏移量DCI0和DCQ0的估计方法为:The method for estimating the DC offset DC I0 and DC Q0 of the I/Q channel baseband signal in the step (2) is: 设采集的基带信号I0(t)和Q0(t)的离散序列分别为I0[n]和Q0[n],拟合圆的圆心坐标为(x0,y0),半径为r,取矩阵A,矩阵m和矩阵B分别为:Let the discrete sequences of the collected baseband signals I 0 (t) and Q 0 (t) be I 0 [n] and Q 0 [n] respectively, the center coordinates of the fitted circle are (x 0 , y 0 ), and the radius is r, take matrix A, matrix m and matrix B as:

Figure FDA0003075486650000011

Figure FDA0003075486650000011

通过参数优化估计得到x0和y0的值;The values of x 0 and y 0 are obtained by parameter optimization estimation; 所述参数优化估计的方法为:min||Am-B||l2或min||Am-B||l1,其中,min(.)表示最小值运算,||.||l2表示2-范数运算,||.||l1表示1-范数运算,求得x0和y0的估计值后,取DCI0=x0,DCQ0=y0The parameter optimization estimation method is: min||Am-B|| l2 or min||Am-B|| l1 , where min(.) represents the minimum value operation, and ||.|| l2 represents the 2-norm Number operation, ||.|| l1 represents 1-norm operation, after obtaining the estimated values of x 0 and y 0 , take DC I0 =x 0 , DC Q0 =y 0 ; 所述步骤(3)中的所述I/Q通道基带信号随时间变化的所述漂移量I_drift(t)和Q_drift(t)求取方法为:The method for obtaining the drifts I_drift(t) and Q_drift(t) of the I/Q channel baseband signal over time in the step (3) is: I_drift(t)=I_trend(t)-mean[I_trend(t)]I_drift(t)=I_trend(t)-mean[I_trend(t)] Q_drift(t)=Q_trend(t)-mean[Q_trend(t)]Q_drift(t)=Q_trend(t)-mean[Q_trend(t)] 式中,mean[·]为取平均值运算,所述I/Q通道基带信号的趋势分量I_trend(t)和Q_trend(t)的求取方法为局部加权线性回归平滑,平滑点数优选为等同于采样频率的整数;In the formula, mean[ ] is the average operation, and the method for obtaining the trend components I_trend(t) and Q_trend(t) of the I/Q channel baseband signal is local weighted linear regression smoothing, and the number of smoothing points is preferably equal to Integer of sampling frequency; 所述步骤(4)中的所述低通滤波的截止频率为5*fobject,其中,fobject是被测目标的振动频率;The cut-off frequency of the low-pass filtering in the step (4) is 5*f object , wherein f object is the vibration frequency of the measured target; 所述步骤(4)中的所述微振动位移的时间序列的计算方法为

Figure FDA0003075486650000021

式中

Figure FDA0003075486650000022

为所述反正切解调计算得到的所述I/Q通道基带信号的相位时间序列,其中,

Figure FDA0003075486650000023

计算方法为

Figure FDA0003075486650000024

The calculation method of the time series of the micro-vibration displacement in the step (4) is:

Figure FDA0003075486650000021

in the formula

Figure FDA0003075486650000022

The phase time series of the I/Q channel baseband signal calculated for the arctangent demodulation, wherein,

Figure FDA0003075486650000023

The calculation method is

Figure FDA0003075486650000024

2.一种基于单频连续波雷达的超微幅振动测量系统,其特征在于,包括雷达前端,平移滑台,采集与处理模块,控制器和振动分析与显示模块,所述雷达前端安装固定在平移滑台上;2. an ultra-micro-amplitude vibration measurement system based on single-frequency continuous wave radar, is characterized in that, comprises radar front-end, translational slide, acquisition and processing module, controller and vibration analysis and display module, described radar front-end is installed and fixed on the translation stage; 所述雷达前端包括微波信号源、功分器、功率放大器、正交移相器、低噪声放大器、混频器、低通滤波器、信号调理电路、发射天线和接收天线,其中,所述微波信号源与所述功分器相连,所述功分器与所述功率放大器和所述正交移相器相连,所述功率放大器与所述发射天线相连,所述接收天线与所述低噪声放大器相连,所述低噪声放大器与所述正交移相器的两个输出端分别与所述混频器相连,两个所述混频器输出端分别与所述低通滤波器相连,所述低通滤波器的输出端与所述信号调理电路相连,所述微波信号源的信号经过所述功分器分为两路,一路经过所述功率放大器连接所述发射天线,一路经过所述正交移相器与所述混频器相连;The radar front end includes a microwave signal source, a power divider, a power amplifier, a quadrature phase shifter, a low-noise amplifier, a frequency mixer, a low-pass filter, a signal conditioning circuit, a transmitting antenna and a receiving antenna, wherein the microwave The signal source is connected to the power divider, the power divider is connected to the power amplifier and the quadrature phase shifter, the power amplifier is connected to the transmitting antenna, and the receiving antenna is connected to the low noise The two output ends of the low-noise amplifier and the quadrature phase shifter are respectively connected to the mixer, and the two output ends of the mixer are respectively connected to the low-pass filter, so The output end of the low-pass filter is connected to the signal conditioning circuit, the signal of the microwave signal source is divided into two paths through the power divider, one path is connected to the transmitting antenna through the power amplifier, and the other a quadrature phase shifter is connected to the mixer; 所述雷达前端,用于产生并发射单频连续波微波信号,并接收目标散射的电磁回波,通过放大、混频及滤波等处理得到I/Q基带信号;The radar front-end is used to generate and transmit single-frequency continuous wave microwave signals, receive electromagnetic echoes scattered by the target, and obtain I/Q baseband signals through processing such as amplification, frequency mixing, and filtering; 所述采集与处理模块,包括信号采集单元和信号处理单元,信号采集单元用于对所述雷达前端输出的基带I/Q信号进行同步采集,信号处理单元用于实现对如权利要求1所述的基于单频连续波雷达的超微幅振动测量方法中初始化直流偏移值的估计,基带I/Q信号漂移量的估计及振动位移时间序列的计算;The acquisition and processing module includes a signal acquisition unit and a signal processing unit, the signal acquisition unit is used for synchronous acquisition of the baseband I/Q signal output by the radar front end, and the signal processing unit is used for The estimation of the initial DC offset value, the estimation of the baseband I/Q signal drift and the calculation of the vibration displacement time series in the ultra-micro-amplitude vibration measurement method based on the single-frequency continuous wave radar; 所述控制器,用于控制顺序执行如权利要求1所述的基于单频连续波雷达的超微幅振动测量方法中的各个步骤,控制所述平移滑台、设置采样参数及执行信号采集、进行初始化直流偏移量的估计运算、I/Q基带信号漂移量的估计及振动位移时间序列的计算与处理;Described controller, is used for controlling to carry out each step in the ultra-micro-amplitude vibration measuring method based on single frequency continuous wave radar as claimed in claim 1 in sequence, controls described translation slide table, sets sampling parameter and carries out signal acquisition, Carry out the estimation operation of the initial DC offset, the estimation of the I/Q baseband signal drift, and the calculation and processing of the vibration displacement time series; 振动分析与显示模块,用于根据需求对测量的超微幅振动时域信息进行振动幅值、频率及时频分布等振动特征的提取与分析,显示信号波形、相关测量及特征提取结果。The vibration analysis and display module is used to extract and analyze vibration characteristics such as vibration amplitude, frequency and frequency distribution of the measured ultra-micro-amplitude vibration time domain information according to the requirements, and display the signal waveform, related measurement and feature extraction results.
CN201911295312.9A 2019-12-16 2019-12-16 A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar Active CN110987150B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911295312.9A CN110987150B (en) 2019-12-16 2019-12-16 A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911295312.9A CN110987150B (en) 2019-12-16 2019-12-16 A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar

Publications (2)

Publication Number Publication Date
CN110987150A CN110987150A (en) 2020-04-10
CN110987150B true CN110987150B (en) 2021-07-23

Family

ID=70094206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911295312.9A Active CN110987150B (en) 2019-12-16 2019-12-16 A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar

Country Status (1)

Country Link
CN (1) CN110987150B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111609920B (en) * 2020-05-13 2021-06-11 上海交通大学 Hand-held microwave vibration measuring system
CN111624588B (en) * 2020-05-27 2023-04-18 上海交通大学 Active sound source positioning and sound signal reconstruction method and device
CN111624598B (en) * 2020-06-11 2023-07-18 石家庄铁道大学 A signal processing method for structural dynamic deformation measurement based on millimeter-wave radar
CN111708019B (en) * 2020-06-12 2023-04-18 上海交通大学 Engineering structure dynamic characteristic monitoring system based on microwave radar
CN112082638B (en) * 2020-10-14 2022-03-29 西安电子科技大学 Non-contact high-frequency micro vibration signal measuring method based on near-field microwave
CN113589249B (en) * 2021-07-16 2024-08-09 中山艾朗格科技有限公司 Signal processing method for calibrating direct current offset of single-frequency continuous wave Doppler radar
CN114608699B (en) * 2022-03-07 2023-02-24 中国矿业大学 Elevator head sheave vibration abnormity detection system and method based on millimeter wave radar

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9602571D0 (en) * 1996-02-08 1996-04-10 Mecon Limited Radar for vibration detection
CN102046076A (en) * 2008-04-03 2011-05-04 Kai医药公司 Non-contact physiological motion sensor and method of use thereof
WO2012152476A1 (en) * 2011-05-12 2012-11-15 Robert Bosch Gmbh Method for detecting gestures
CN203661075U (en) * 2013-12-18 2014-06-18 中国科学院微电子研究所 Non-contact human body surface micro-motion information detection device
CN103919527A (en) * 2013-01-14 2014-07-16 财团法人工业技术研究院 Motion/disturbance signal detection system and method
CN103948381A (en) * 2014-04-09 2014-07-30 浙江大学 Imaging method for large dynamic range Doppler biological signal
CN106644030A (en) * 2016-08-31 2017-05-10 上海交通大学 Doppler radar-based non-contact type vibration measuring method
CN106805940A (en) * 2015-12-02 2017-06-09 由国峰 A kind of continuous wave bioradar sign detection means

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9602571D0 (en) * 1996-02-08 1996-04-10 Mecon Limited Radar for vibration detection
CN102046076A (en) * 2008-04-03 2011-05-04 Kai医药公司 Non-contact physiological motion sensor and method of use thereof
WO2012152476A1 (en) * 2011-05-12 2012-11-15 Robert Bosch Gmbh Method for detecting gestures
CN103919527A (en) * 2013-01-14 2014-07-16 财团法人工业技术研究院 Motion/disturbance signal detection system and method
CN203661075U (en) * 2013-12-18 2014-06-18 中国科学院微电子研究所 Non-contact human body surface micro-motion information detection device
CN103948381A (en) * 2014-04-09 2014-07-30 浙江大学 Imaging method for large dynamic range Doppler biological signal
CN106805940A (en) * 2015-12-02 2017-06-09 由国峰 A kind of continuous wave bioradar sign detection means
CN106644030A (en) * 2016-08-31 2017-05-10 上海交通大学 Doppler radar-based non-contact type vibration measuring method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An Effective Accuracy Evaluation Method for LFMCW Radar Displacement Monitoring With Phasor Statistical Analysis;Xiong Yuyong, et al;《IEEE SENSORS JOURNAL》;20191215;第19卷(第24期);第12224-12234页 *
Evaluation of FMCW Radar for Vibration Sensing in Industrial Environments;Christian Zeintl, et al;《2019 29th International Conference Radioelektronika》;20190610;第1-5页 *
一种强杂波环境下雷达目标微振动信号提取方法;曾涛 等;《信号处理》;20190630;第35卷(第6期);第935-943页 *
基于参数化解调的旋转目标微多普勒频率提取方法;陈是扦 等;《上海航天》;20181031;第35卷(第5期);第39-44页 *

Also Published As

Publication number Publication date
CN110987150A (en) 2020-04-10

Similar Documents

Publication Publication Date Title
CN110987150B (en) 2021-07-23 A method and system for ultra-micro-amplitude vibration measurement based on single-frequency continuous wave radar
Scherr et al. 2014 An efficient frequency and phase estimation algorithm with CRB performance for FMCW radar applications
JP3222075U (en) 2019-07-11 Integrated multi-channel RF circuit with phase sensing capability
CN107607923B (en) 2021-07-30 Vibration monitoring system and signal processing method based on LFMCW radar
EP3051308B1 (en) 2020-07-29 Radar apparatus and object detection method
CN110095759B (en) 2023-07-07 Spectral Estimation of Noise in Radar Installations
JP5352959B2 (en) 2013-11-27 Direct sequence spread spectrum radar, method used in radar and computer program
US9998240B2 (en) 2018-06-12 Electric field intensity distribution measurement device and electric field intensity distribution measurement method
JP6573599B2 (en) 2019-09-11 Method and apparatus for measuring broadband measurement signals
CN103630888A (en) 2014-03-12 High-precision real-time microwave velocity and distance measuring device based on symmetrical triangle LFMCW (Linear Frequency Modulation Continuous Wave) radar
CN115508820A (en) 2022-12-23 Target detection method of linear frequency modulation pulse radar
US20160139264A1 (en) 2016-05-19 Method for Determining a Distance Between an FMCW Ranging Device and a Target
CN106199582A (en) 2016-12-07 The method being applied to the two-frequency CW radar human body locating and tracking of indoor
US20160135694A1 (en) 2016-05-19 Medical radar method and system
CN113433540B (en) 2023-12-12 Phase difference detection equipment for FMCW modulation
JP2016526150A5 (en) 2019-03-28
US11346921B2 (en) 2022-05-31 Hidden reception monitoring during a signal reception operation
JP6864112B2 (en) 2021-04-21 Exploration of concrete structure using electromagnetic waves
US11206090B2 (en) 2021-12-21 Error measurement method using a time-variant stopband test signal
CN103795675B (en) 2017-04-12 Method for carrier synchronization under low signal-to-noise ratio and high dynamic environment
JP5354126B2 (en) 2013-11-27 Direct sequence spread spectrum radar
CN114217266B (en) 2024-06-11 Processing method for reducing channel error based on antenna data time domain fitting
Strombeck et al. 2019 AMCW radar of micrometer accuracy distance measurement and monitoring
Kim et al. 2020 Fine Doppler frequency estimation of windowed complex sinusoidal signals for radar systems
Piccinni et al. 2019 Distance Estimation Algorithm for Wireless Localization Systems Based on Lyapunov Sensitivity Theory

Legal Events

Date Code Title Description
2020-04-10 PB01 Publication
2020-04-10 PB01 Publication
2020-05-05 SE01 Entry into force of request for substantive examination
2020-05-05 SE01 Entry into force of request for substantive examination
2021-07-23 GR01 Patent grant
2021-07-23 GR01 Patent grant