CN110220586A - Vibration monitoring method and system based on millimeter wave - Google Patents
Vibration monitoring method and system based on millimeter wave Download PDFInfo
- Publication number
- CN110220586A CN110220586A CN201910530699.5A CN201910530699A CN110220586A CN 110220586 A CN110220586 A CN 110220586A CN 201910530699 A CN201910530699 A CN 201910530699A CN 110220586 A CN110220586 A CN 110220586A
- Authority
- CN
- China
- Prior art keywords
- signal
- measured
- vibration source
- vibration
- target composite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明实施例提供一种基于毫米波的振动监测方法及系统，包括：获取待测振动源的目标复合无线信号；对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。本发明实施例通过将待测振动源及周围物体发送的毫米波雷达信号中各种干扰信号去除，得到包含待测振动源振动信号相关信息的视距反射信号，从而得到待测振动源的振幅和频率，以用于对待测振动源进行振动监测，提高了振动监测的精准度，大幅度降低了振动监测的误差。
Embodiments of the present invention provide a vibration monitoring method and system based on millimeter waves, including: acquiring a target composite wireless signal of the vibration source to be tested; extracting and processing the target composite wireless signal, and obtaining a visual signal of the vibration source to be tested The distance reflection signal; according to the line-of-sight reflection signal, the amplitude and frequency of the vibration source to be measured are obtained, so as to wirelessly monitor the vibration source to be measured. In the embodiment of the present invention, by removing various interference signals from the millimeter-wave radar signal sent by the vibration source to be measured and surrounding objects, a line-of-sight reflection signal containing information related to the vibration signal of the vibration source to be measured is obtained, thereby obtaining the amplitude of the vibration source to be measured and frequency for vibration monitoring of the vibration source to be measured, which improves the accuracy of vibration monitoring and greatly reduces the error of vibration monitoring.
Description
技术领域technical field
本发明涉及机械振动无线监测技术领域，尤其涉及一种基于毫米波的振动监测方法及系统。The invention relates to the technical field of mechanical vibration wireless monitoring, in particular to a vibration monitoring method and system based on millimeter waves.
背景技术Background technique
在现代工业中，诸如振幅、振动频率等振动参数指标，反映了机器的运行工况，并能够辅助机械故障诊断，甚至提前进行故障预测与报警。因此，准确的振动监测一直是工业自动化场景中不可缺少的一个环节。工业环境中的振动监测有如下三个特点：首先，工业设备在正常工作时振幅极其微小，约毫米级甚至是亚毫米级，微小振幅给振动监测精度设立了很高的技术门槛；其次，现代工业设备通常24小时不停机工作，倘若出现任何机械故障可能导致整个产线系统停滞，因此，振动监测应该具备高实时性、全天候监测的特点；最后，现代工业中广泛部署的泵机、马达等均有振动监测需求，导致高昂监测成本。In modern industry, vibration parameter indicators such as amplitude and vibration frequency reflect the operating conditions of the machine, and can assist in the diagnosis of mechanical faults, and even perform fault prediction and alarm in advance. Therefore, accurate vibration monitoring has always been an indispensable link in industrial automation scenarios. Vibration monitoring in industrial environments has the following three characteristics: First, the amplitude of industrial equipment is extremely small during normal operation, about millimeters or even sub-millimeters, and the small amplitude sets a high technical threshold for vibration monitoring accuracy; secondly, modern Industrial equipment usually works 24 hours without stopping. Any mechanical failure may cause the entire production line system to stagnate. Therefore, vibration monitoring should have the characteristics of high real-time and all-weather monitoring; finally, pumps, motors, etc. widely deployed in modern industries Both have vibration monitoring requirements, resulting in high monitoring costs.
在目前普遍的工业场景中，诸多非核心的机械设备在原厂设计装配的环节并未安装振动监测传感器，额外加装自动化振动监测装置首先需要考虑的问题便是部署成本和维护成本；传统的基于压电陶瓷的振动监测传感器需要繁琐部署地部署电源线和通讯线等，并且需要固定在已经运转的机器表面，使得其难以适用于普遍的传感器加装场景；高精度的振动监测传感器原理一般基于激光测距或激光干涉，然而激光器件成本过高且对部署环境要求苛刻，难以适用于工业场景大规模部署；通过用新型无线感知技术，如基于无线射频识别RFID技术的振动测量方法，这类方法只能测量精度有限的振动频率，难以恢复出亚毫米级的振幅。因此，上述现有的自动化振动监测技术均未被广泛地应用在工业场景中。而目前绝大多数企业和工厂选择人工手持测振仪的方式进行振动监测，尽管手持振动仪价格低廉且监测准确度适中，人工巡检的方式却具有难以全天候监控以及实时报警，并且人为主观测量误差较大。In the current common industrial scenarios, many non-core mechanical equipment are not installed with vibration monitoring sensors in the original factory design and assembly process. The first issue that needs to be considered when installing additional automatic vibration monitoring devices is the deployment cost and maintenance cost; Piezoelectric ceramic vibration monitoring sensors require cumbersome deployment of power lines and communication lines, etc., and need to be fixed on the surface of already running machines, making it difficult to apply to common sensor installation scenarios; the principle of high-precision vibration monitoring sensors is generally based on Laser ranging or laser interference, however, the cost of laser devices is too high and the deployment environment is harsh, so it is difficult to apply to large-scale deployment in industrial scenarios; by using new wireless sensing technologies, such as vibration measurement methods based on radio frequency identification (RFID) technology, such The method can only measure the vibration frequency with limited accuracy, and it is difficult to recover the submillimeter amplitude. Therefore, none of the above-mentioned existing automated vibration monitoring technologies has been widely applied in industrial scenarios. At present, most enterprises and factories choose to carry out vibration monitoring by hand-held vibrometer. Although the hand-held vibrometer is cheap and the monitoring accuracy is moderate, the manual inspection method has difficulties in all-weather monitoring and real-time alarm, and human subjective measurement. The error is large.
因此，现在亟需一种基于毫米波的振动监测方法及系统来解决上述问题。Therefore, there is an urgent need for a vibration monitoring method and system based on millimeter waves to solve the above problems.
发明内容Contents of the invention
针对现有技术存在的问题，本发明实施例提供一种基于毫米波的振动监测方法及系统。In view of the problems existing in the prior art, embodiments of the present invention provide a vibration monitoring method and system based on millimeter waves.
第一方面，本发明实施例提供了一种基于毫米波的振动监测方法，包括：In the first aspect, an embodiment of the present invention provides a vibration monitoring method based on millimeter waves, including:
获取待测振动源的目标复合无线信号；Obtain the target composite wireless signal of the vibration source to be tested;
对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；Extracting and processing the target composite wireless signal to obtain a line-of-sight reflection signal of the vibration source to be measured;
根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。Acquiring the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be measured.
进一步地，在所述获取待测振动源的目标复合无线信号之前，所述方法还包括：Further, before the acquisition of the target composite wireless signal of the vibration source to be tested, the method also includes:
获取待测振动源和周围物体反射的毫米波雷达信号；Obtain millimeter-wave radar signals reflected by the vibration source to be tested and surrounding objects;
对所述毫米波雷达信号进行提取处理，获取相同传播距离且传播路径不同的复合信号，以得到所述待测振动源的目标复合无线信号。The millimeter-wave radar signal is extracted and processed to obtain composite signals with the same propagation distance and different propagation paths, so as to obtain the target composite wireless signal of the vibration source to be measured.
具体地，所述目标复合无线信号包括视距反射信号和多径反射干扰信号，其中，所述多径反射干扰信号包括目标相关多径反射信号和环境多径反射信号。Specifically, the target composite wireless signal includes a line-of-sight reflection signal and a multipath reflection interference signal, wherein the multipath reflection interference signal includes a target-related multipath reflection signal and an environment multipath reflection signal.
进一步地，所述对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号，包括：Further, the extracting and processing the target composite wireless signal to obtain the line-of-sight reflection signal of the vibration source to be measured includes:
若发送-接收天线为一组时，获取所述目标复合无线信号；If the transmitting-receiving antenna is a group, acquiring the target composite wireless signal;
将所述目标复合无线信号投射到二维复平面，并得到所述目标复合无线信号对应的圆弧；Projecting the target composite wireless signal to a two-dimensional complex plane, and obtaining an arc corresponding to the target composite wireless signal;
根据圆弧几何拟合算法，获取二维复平面上所述目标复合无线信号对应圆弧所在圆的圆心坐标和半径；According to the arc geometric fitting algorithm, the center coordinates and the radius of the circle where the target composite wireless signal corresponds to the arc are obtained on the two-dimensional complex plane;
根据所述圆心坐标，通过坐标轴平移，将所述目标复合无线信号中多径反射干扰信号去除，得到所述待测振动源的视距反射信号。According to the coordinates of the center of the circle, the multipath reflection interference signal in the target composite wireless signal is removed by translation of the coordinate axis, so as to obtain the line-of-sight reflection signal of the vibration source to be measured.
进一步地，所述根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测，包括：Further, the acquisition of the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal to wirelessly monitor the vibration source to be measured includes:
获取所述视距反射信号的相位角序列；Obtaining the phase angle sequence of the line-of-sight reflection signal;
对所述相位角序列进行低通滤波直流分量处理，得到所述视距反射信号的时域表示；performing low-pass filtering DC component processing on the phase angle sequence to obtain a time-domain representation of the line-of-sight reflection signal;
根据最大似然参数估计方法，对所述视距反射信号的时域表示进行估算，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。According to the maximum likelihood parameter estimation method, the time-domain representation of the line-of-sight reflection signal is estimated, and the amplitude and frequency of the vibration source to be measured are obtained, so as to wirelessly monitor the vibration source to be measured.
进一步地，所述对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号，所述方法还包括：Further, the extracting and processing the target composite wireless signal to obtain the line-of-sight reflection signal of the vibration source to be measured, the method also includes:
S1，若发送-接收天线为多组时，获取多个目标复合无线信号；S1, if there are multiple sets of sending-receiving antennas, acquiring multiple target composite wireless signals;
S2，将多个目标复合无线信号投射到二维复平面，并得到每个目标复合无线信号对应的圆弧；根据圆弧几何拟合算法，获取二维复平面上每个目标复合无线信号对应圆弧所在圆的圆心坐标和半径，并根据圆心坐标拟合得到多个目标复合无线信号的共圆参数，并将所述共圆参数作为半圆约束条件，所述共圆参数包括共圆圆心坐标和半径；S2. Project multiple target composite wireless signals onto the two-dimensional complex plane, and obtain the arc corresponding to each target composite wireless signal; according to the circular arc geometric fitting algorithm, obtain the corresponding arc of each target composite wireless signal on the two-dimensional complex plane. The center coordinates and radius of the circle where the arc is located, and according to the center coordinates, the co-circle parameters of multiple target composite wireless signals are obtained, and the co-circle parameters are used as semicircle constraints, and the co-circle parameters include the co-circle center coordinates and radius;
S3，根据半径约束条件和所述共圆约束条件，通过圆弧几何拟合算法，获取在约束条件向下每个目标复合无线信号对应圆弧所在圆的圆心参数和半径，以获取新的共圆参数，若新的共圆参数满足预设阈值，则通过坐标轴平移，将每个目标复合无线信号中多径反射干扰信号去除，以得到每个接收天线对应的所述待测振动源的视距反射信号；若新的共圆参数不满足预设阈值，则将新的共圆参数作为下一次迭代过程的共圆约束条件，再次执行步骤S3，直至迭代得到的共圆参数满足预设阈值。S3. According to the radius constraint condition and the co-circle constraint condition, through the arc geometric fitting algorithm, obtain the center parameter and radius of the circle where each target composite wireless signal corresponds to the arc under the constraint condition, so as to obtain a new co-circle constraint condition circle parameter, if the new co-circle parameter satisfies the preset threshold, the multipath reflection interference signal in each target composite wireless signal is removed by coordinate axis translation, so as to obtain the vibration source to be measured corresponding to each receiving antenna Line-of-sight reflection signal; if the new co-circle parameter does not meet the preset threshold, the new co-circle parameter is used as the co-circle constraint condition for the next iteration process, and step S3 is executed again until the co-circle parameter obtained by iteration meets the preset threshold.
进一步地，所述根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测，还包括：Further, the acquiring the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be measured, further includes:
获取每个接收天线对应的视距反射信号的相位角序列；Obtain the phase angle sequence of the line-of-sight reflection signal corresponding to each receiving antenna;
对每个相位角序列进行低通滤波直流分量处理，得到每个视距反射信号的时域表示；Perform low-pass filtering DC component processing on each phase angle sequence to obtain the time domain representation of each line-of-sight reflection signal;
根据四分位均值数据融合算法，对每个视距反射信号的时域表示进行融合，得到融合后的时域表示；According to the quartile mean data fusion algorithm, the time domain representation of each line-of-sight reflection signal is fused to obtain the fused time domain representation;
根据最大似然参数估计方法，对融合后的时域表示进行估算，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。According to the maximum likelihood parameter estimation method, the fused time domain representation is estimated, and the amplitude and frequency of the vibration source to be measured are obtained, so as to wirelessly monitor the vibration source to be measured.
第二方面，本发明实施例提供了一种基于毫米波的振动监测系统，包括：In the second aspect, an embodiment of the present invention provides a vibration monitoring system based on millimeter waves, including:
复合信号获取模块，用于获取待测振动源的目标复合无线信号；Composite signal acquisition module, used to obtain the target composite wireless signal of the vibration source to be measured;
视距反射信号提取模块，用于对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；The line-of-sight reflection signal extraction module is used to extract and process the target composite wireless signal, and obtain the line-of-sight reflection signal of the vibration source to be measured;
监测模块，用于根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。The monitoring module is configured to obtain the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be measured.
第三方面，本发明实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如第一方面所提供的方法的步骤。In a third aspect, an embodiment of the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and operable on the processor. When the processor executes the program, the computer program described in the first aspect is implemented. The steps of the provided method.
第四方面，本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如第一方面所提供的方法的步骤。In a fourth aspect, an embodiment of the present invention provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method provided in the first aspect are implemented.
本发明实施例提供的一种基于毫米波的振动监测方法及系统，通过将待测振动源及周围物体发送的毫米波雷达信号中各种干扰信号去除，得到包含待测振动源振动信号相关信息的视距反射信号，从而得到待测振动源的振幅和频率，以用于对待测振动源进行振动监测，提高了振动监测的精准度，大幅度降低了振动监测的误差。In the millimeter-wave-based vibration monitoring method and system provided by the embodiments of the present invention, various interference signals from the millimeter-wave radar signals sent by the vibration source to be measured and surrounding objects are removed to obtain relevant information including the vibration signal of the vibration source to be measured The line-of-sight reflection signal is used to obtain the amplitude and frequency of the vibration source to be measured, which is used for vibration monitoring of the vibration source to be measured, which improves the accuracy of vibration monitoring and greatly reduces the error of vibration monitoring.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description These are some embodiments of the present invention. Those skilled in the art can also obtain other drawings based on these drawings without creative work.
图1为本发明实施例提供的基于毫米波的振动监测方法的流程示意图；FIG. 1 is a schematic flow diagram of a vibration monitoring method based on millimeter waves provided by an embodiment of the present invention;
图2为本发明实施例提供的多径效应环境中毫米波雷达信号的反射示意图；FIG. 2 is a schematic diagram of the reflection of a millimeter-wave radar signal in a multipath environment provided by an embodiment of the present invention;
图3为本发明实施例提供的亚毫米级振动的信号提取示意图；Fig. 3 is a schematic diagram of signal extraction of submillimeter-level vibration provided by an embodiment of the present invention;
图4为本发明实施例提供的振幅与频率的测量误差示意图；FIG. 4 is a schematic diagram of measurement errors of amplitude and frequency provided by an embodiment of the present invention;
图5为本发明另一实施例提供的振幅与频率的测量误差示意图；Fig. 5 is a schematic diagram of measurement errors of amplitude and frequency provided by another embodiment of the present invention;
图6为本发明又一实施例提供的振幅与频率的测量误差示意图；Fig. 6 is a schematic diagram of measurement errors of amplitude and frequency provided by another embodiment of the present invention;
图7为本发明实施例提供的基于毫米波的振动监测系统的结构示意图；FIG. 7 is a schematic structural diagram of a vibration monitoring system based on millimeter waves provided by an embodiment of the present invention;
图8为本发明实施例提供的电子设备结构示意图。FIG. 8 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments It is a part of embodiments of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.
针对现有技术存在的问题，本发明实施例提出了一种基于毫米波的振动监测方法及系统，仅需对准待测机器(待测振动源)，通过远距离地发送毫米波频段的无线信号，接收由振动源及周围环境反射的无线信号，并提取其中包含的微小振动信号即可，需要说明的是，本发明实施例提供的方法及系统，对于具有毫米级以及亚毫米级振幅的微小振动信号均适用。本发明实施例具有如下多点应用优势：首先，本发明实施例为非侵入式部署且部署安装成本低廉，基于无线远程监测的特性使其无需额外固定在待测机器表面，也无需在待测机器附近布线；其次，本发明实施例能够准确地提取亚毫米级振动的振幅与频率，在两米的测量距离时，本发明实施例的振动频率和振幅测量的误差分别约1Hz和10um，能够满足工业振动全天候实时监测的实际需求；最后，本发明实施例基于无线信号空间传播特性，使得本发明实施例能够对多振动目标同时进行振动监测。Aiming at the problems existing in the prior art, the embodiment of the present invention proposes a vibration monitoring method and system based on millimeter wave, which only needs to align the machine under test (vibration source to be tested), and transmit wireless Signal, just receive the wireless signal reflected by the vibration source and the surrounding environment, and extract the tiny vibration signal contained therein. It should be noted that the method and system provided by the embodiment of the present invention are suitable for Small vibration signals are applicable. The embodiment of the present invention has the following multi-point application advantages: First, the embodiment of the present invention is a non-invasive deployment with low deployment and installation costs. Based on the characteristics of wireless remote monitoring, it does not need to be additionally fixed on the surface of the machine to be tested, nor does it need to be installed on the surface of the machine to be tested. Wiring near the machine; secondly, the embodiments of the present invention can accurately extract the amplitude and frequency of submillimeter-level vibrations. At a measurement distance of two meters, the errors of the vibration frequency and amplitude measurements of the embodiments of the present invention are about 1 Hz and 10 um respectively, which can It satisfies the actual demand for all-weather real-time monitoring of industrial vibration; finally, the embodiment of the present invention is based on the spatial propagation characteristics of wireless signals, so that the embodiment of the present invention can simultaneously monitor the vibration of multiple vibration targets.
图1为本发明实施例提供的基于毫米波的振动监测方法的流程示意图，如图1所示，本发明实施例提供了一种基于毫米波的振动监测方法，包括：Fig. 1 is a schematic flowchart of a vibration monitoring method based on millimeter waves provided by an embodiment of the present invention. As shown in Fig. 1, an embodiment of the present invention provides a vibration monitoring method based on millimeter waves, including:
步骤101，获取待测振动源的目标复合无线信号。Step 101, acquiring the target composite wireless signal of the vibration source to be tested.
在本发明实施例中，首先对待测振动源发送毫米波雷达信号，并接收到通过待测振动源以及周围物体反射的毫米波雷达信号。由于微小振动是一个往复运动，会周期性地改变无线信号的传播路径长度，从而改变无线信号的相位；并且，由于毫米波频段的无线信号波长约1-10mm，因此可以通过分析其相位变化提取出对应的振动信号。然而，由无线信号在空间传播的特性可知，发射出去的毫米波无线信号，会被周围所有的物体反射，导致接收到的无线信号是一个复合信号，该复合信号中包含了多路不同传播距离的信号，同时，在该复合信号中也包含了由多条传播路径产生的相同传播距离的信号。在本发明实施例中，对待测振动源和周围物体反射的毫米波雷达信号进行信号提取，首先将不同传播距离的干扰信号去除，从而提取得到相同传播距离不同传播路径的复合信号，即待测振动源的目标复合无线信号。In the embodiment of the present invention, the millimeter-wave radar signal is first sent to the vibration source to be tested, and the millimeter-wave radar signal reflected by the vibration source to be tested and surrounding objects is received. Since the tiny vibration is a reciprocating motion, it will periodically change the propagation path length of the wireless signal, thereby changing the phase of the wireless signal; and, since the wavelength of the wireless signal in the millimeter wave band is about 1-10mm, it can be extracted by analyzing its phase change. output the corresponding vibration signal. However, according to the characteristics of wireless signal propagation in space, the transmitted millimeter wave wireless signal will be reflected by all surrounding objects, resulting in the received wireless signal being a composite signal, which contains multiple channels with different propagation distances. At the same time, the composite signal also includes signals with the same propagation distance generated by multiple propagation paths. In the embodiment of the present invention, the signal extraction is performed on the millimeter-wave radar signal reflected by the vibration source to be tested and the surrounding objects. First, the interference signals with different propagation distances are removed, thereby extracting a composite signal with the same propagation distance and different propagation paths, that is, the signal to be tested The target composite wireless signal of the vibration source.
步骤102，对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号。Step 102, extracting and processing the target composite wireless signal to obtain a line-of-sight reflection signal of the vibration source to be tested.
在本发明实施例中，在得到目标复合无线信号之后，目标复合无线信号中携带有视距反射信号和多径反射干扰信号，因此，需要将多径反射干扰信号去除，从而得到待测振动源在视距路径上的反射信号，即待测振动源的视距反射信号。In the embodiment of the present invention, after the target composite wireless signal is obtained, the target composite wireless signal carries line-of-sight reflection signals and multipath reflection interference signals. Therefore, it is necessary to remove the multipath reflection interference signals to obtain the vibration source to be tested. The reflection signal on the line-of-sight path, that is, the line-of-sight reflection signal of the vibration source to be measured.
步骤103，根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。Step 103: Acquire the amplitude and frequency of the vibration source to be tested according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be tested.
在本发明实施例中，在去除了其他干扰信号之后，获取到待测振动源的视距反射信号，然后计算得到该视距反射信号的相位角序列信息，并对相位角序列信息进行低通滤波过滤直流分量处理，得到待测振动源对应振动信号的时域表示，最后通过对振动信号的时域表示进行最大似然参数估计，从而得到该振动信号的振幅和频率，以用于对待测振动源进行分析，得到对应的振动监测数据。In the embodiment of the present invention, after removing other interference signals, the line-of-sight reflection signal of the vibration source to be measured is obtained, and then the phase angle sequence information of the line-of-sight reflection signal is calculated, and the phase angle sequence information is low-passed Filter and filter the DC component to obtain the time-domain representation of the vibration signal corresponding to the vibration source to be tested. Finally, the maximum likelihood parameter estimation is performed on the time-domain representation of the vibration signal to obtain the amplitude and frequency of the vibration signal for use in the test. The vibration source is analyzed to obtain the corresponding vibration monitoring data.
本发明实施例提供的基于毫米波的振动监测方法，通过将待测振动源及周围物体发送的毫米波雷达信号中各种干扰信号去除，得到包含待测振动源振动信号相关信息的视距反射信号，从而得到待测振动源的振幅和频率，以用于对待测振动源进行振动监测，提高了振动监测的精准度，大幅度降低了振动监测的误差。The millimeter-wave-based vibration monitoring method provided by the embodiment of the present invention removes various interference signals from the millimeter-wave radar signals sent by the vibration source to be measured and surrounding objects, and obtains the line-of-sight reflection containing information related to the vibration signal of the vibration source to be measured. Signal, so as to obtain the amplitude and frequency of the vibration source to be measured, which is used for vibration monitoring of the vibration source to be measured, which improves the accuracy of vibration monitoring and greatly reduces the error of vibration monitoring.
在上述实施例的基础上，其特征在于，在所述获取待测振动源的目标复合无线信号之前，所述方法还包括：On the basis of the above-mentioned embodiments, it is characterized in that, before the acquisition of the target composite wireless signal of the vibration source to be tested, the method further includes:
获取待测振动源和周围物体反射的毫米波雷达信号；Obtain millimeter-wave radar signals reflected by the vibration source to be tested and surrounding objects;
对所述毫米波雷达信号进行提取处理，获取相同传播距离且传播路径不同的复合信号，以得到所述待测振动源对应的目标复合无线信号。The millimeter-wave radar signal is extracted and processed to obtain composite signals with the same propagation distance and different propagation paths, so as to obtain a target composite wireless signal corresponding to the vibration source to be measured.
图2为本发明实施例提供的多径效应环境中毫米波雷达信号的反射示意图，可参考图2所示，在本发明实施例中，向待测振动源发送毫米波雷达信号，并对接收到的振动源以及周围物体反射的毫米波雷达信号进行分析。由于接受到的反射信号是一个复合信号，在本发明实施例中，对待测振动源和周围物体反射的毫米波雷达信号中的各个信号进行了相应解释，参照图2所示，在待测振动源和周围物体反射的毫米波雷达信号中，包含了4种路径的信号，其中，路径1和路径2表示不同传播距离的信号复合；路径1、路径3和路径4表示相同传播距离且传播路径不同的信号复合。需要说明的是，在本发明实施例中，忽略因微小振动引起的信号传播距离的改变，路径1、路径3和路径4的传播距离相差非常微小，可以理解为这三种路径的传播距离相等，而路径2的传播距离，与路径1、路径3以及路径4的传播距离相差很大，因此，首先需要将路径2类型的信号去除。Figure 2 is a schematic diagram of the reflection of the millimeter-wave radar signal in the multipath effect environment provided by the embodiment of the present invention, as shown in Figure 2. In the embodiment of the present invention, the millimeter-wave radar signal is sent to the vibration source to be tested, and received The vibration source and the millimeter-wave radar signal reflected by the surrounding objects are analyzed. Since the received reflected signal is a composite signal, in the embodiment of the present invention, each signal in the millimeter-wave radar signal reflected by the vibration source to be measured and the surrounding objects is explained accordingly. Referring to Figure 2, in the vibration to be measured The millimeter-wave radar signal reflected by the source and surrounding objects contains signals of four paths, among which, path 1 and path 2 represent signal composites with different propagation distances; path 1, path 3 and path 4 represent the same propagation distance and the propagation path Composite of different signals. It should be noted that, in the embodiment of the present invention, the change of the signal propagation distance caused by tiny vibrations is ignored, and the propagation distances of path 1, path 3, and path 4 are very slightly different, which can be understood as the propagation distances of these three paths are equal , and the propagation distance of path 2 is very different from the propagation distances of path 1, path 3, and path 4. Therefore, it is first necessary to remove the signal of path 2 type.
进一步地，在本发明实施例中，首先需要对待测振动源和周围物体反射的毫米波雷达信号进行提取处理，将不同传播距离的信号分离去除，从而得到只包含相同传播距离且传播路径不同的复合信号，即目标复合无线信号。具体地，在本发明实施例中，将接收到的待测振动源和周围物体反射的毫米波雷达信号y(t)表示为：Furthermore, in the embodiment of the present invention, it is first necessary to extract and process the millimeter-wave radar signals reflected by the vibration source to be measured and surrounding objects, and separate and remove signals with different propagation distances, so as to obtain Composite signal, that is, the target composite wireless signal. Specifically, in the embodiment of the present invention, the millimeter-wave radar signal y(t) received by the vibration source to be measured and reflected by surrounding objects is expressed as:
其中，F和K均表示雷达硬件参数，F表示毫米波雷达线性调频连续波的起始频率，K表示频率线性变换的斜率，p表示不同传播的记号，q表示相同传播距离且不同传播路径的记号，Ap,q表示第p个传播距离中第q个传播路径上的信号强度衰减因子，Rp,q表示第p个传播距离中第q个传播路径的长度，e表示复数，t表示时间，c表示光速，j表示复数的虚数单位。在本发明实施例中，忽略微小振动引起的信号传播距离的改变，则Rp,1＝Rp,2＝…＝Rp,q＝…＝Rp，通过对接收到的待测振动源和周围物体反射的毫米波雷达信号y(t)进行快速傅里叶变换得到Y(ω),根据雷达信号原理，Y(ω)不同频率分量对应不同传播距离的信号复合，假设雷达和待测振动源之间的距离为Rp*，首先确定距离Rp*对应的频率分量然后，将该频率分量代入Y(ω)中，从而得到用于表征相同传播距离且不同传播路径的信号复合Y(ω*)，公式表示为：Among them, F and K both represent the radar hardware parameters, F represents the starting frequency of the millimeter-wave radar chirp continuous wave, K represents the slope of the frequency linear transformation, p represents the sign of different propagation, and q represents the same propagation distance and different propagation paths Notation, A p,q represents the signal strength attenuation factor on the qth propagation path in the pth propagation distance, R p,q represents the length of the qth propagation path in the pth propagation distance, e represents a complex number, and t represents Time, c represents the speed of light, and j represents the imaginary unit of complex numbers. In the embodiment of the present invention, ignoring the change of signal propagation distance caused by tiny vibrations, then R p,1 =R p,2 =...=R p,q =...=R p , by analyzing the received vibration source to be measured The millimeter-wave radar signal y(t) reflected by the surrounding objects is subjected to fast Fourier transformation to obtain Y(ω). According to the principle of radar signals, different frequency components of Y(ω) correspond to signals with different propagation distances. The distance between the vibration sources is R p* , first determine the frequency component corresponding to the distance R p* Then, the frequency component is substituted into Y(ω), so as to obtain the signal composite Y(ω*) used to characterize the same propagation distance and different propagation paths, the formula is expressed as:
具体地，在本发明实施例中，对每个雷达线性调频周期T内的所有采样点进行快速傅里叶变换，再根据对应的频率分量获取Y(ω*)的一个采样点，最后通过组合多个连续周期内分别提取得到采样点，从而得到关于Y(ω*)的连续采样，以用于后续的步骤处理。Specifically, in the embodiment of the present invention, fast Fourier transform is performed on all sampling points in each radar chirp period T, and then a sampling point of Y(ω*) is obtained according to the corresponding frequency component, and finally by combining Sampling points are respectively extracted in a plurality of consecutive periods, so as to obtain continuous sampling about Y(ω*), which is used for subsequent step processing.
在上述实施例的基础上，所述目标复合无线信号包括视距反射信号和多径反射干扰信号，其中，所述多径反射干扰信号包括目标相关多径反射信号和环境多径反射信号。On the basis of the above embodiments, the target composite wireless signal includes a line-of-sight reflection signal and a multipath reflection interference signal, wherein the multipath reflection interference signal includes a target-related multipath reflection signal and an environment multipath reflection signal.
在上述实施例的基础上，所述对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号，包括：On the basis of the above-mentioned embodiments, the extraction processing of the target composite wireless signal to obtain the line-of-sight reflection signal of the vibration source to be measured includes:
若发送-接收天线为一组时，获取所述目标复合无线信号；If the transmitting-receiving antenna is a group, acquiring the target composite wireless signal;
将所述目标复合无线信号投射到二维复平面，并得到所述目标复合无线信号对应的圆弧；Projecting the target composite wireless signal to a two-dimensional complex plane, and obtaining an arc corresponding to the target composite wireless signal;
根据圆弧几何拟合算法，获取二维复平面上所述目标复合无线信号对应圆弧所在圆的圆心坐标和半径；According to the arc geometric fitting algorithm, the center coordinates and the radius of the circle where the target composite wireless signal corresponds to the arc are obtained on the two-dimensional complex plane;
根据所述圆心坐标，通过坐标轴平移，将所述目标复合无线信号中多径反射干扰信号去除，得到所述待测振动源的视距反射信号。According to the coordinates of the center of the circle, the multipath reflection interference signal in the target composite wireless signal is removed by translation of the coordinate axis, so as to obtain the line-of-sight reflection signal of the vibration source to be measured.
在上述实施例的基础上，所述根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测，包括：On the basis of the above embodiments, the acquisition of the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal to wirelessly monitor the vibration source to be measured includes:
获取所述视距反射信号的相位角序列；Obtaining the phase angle sequence of the line-of-sight reflection signal;
对所述相位角序列进行低通滤波直流分量处理，得到所述视距反射信号的时域表示；performing low-pass filtering DC component processing on the phase angle sequence to obtain a time-domain representation of the line-of-sight reflection signal;
根据最大似然参数估计方法，对所述视距反射信号的时域表示进行估算，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。According to the maximum likelihood parameter estimation method, the time-domain representation of the line-of-sight reflection signal is estimated, and the amplitude and frequency of the vibration source to be measured are obtained, so as to wirelessly monitor the vibration source to be measured.
在本发明实施例中，需要将多径反射干扰信号从目标复合无线信号中分离去除，从而得到单一的视距反射信号，其中，多径反射干扰信号包括目标相关多径反射信号和环境多径反射信号。图3为本发明实施例提供的亚毫米级振动的信号提取示意图，可参考图3所示，在本发明实施例中，视距反射信号、目标相关多径反射信号和环境多径反射信号叠加在一起，三者叠加产生的复合信号在二维复平面上表示为向量的矢量相加，使得该目标复合无线信号投影到二维复平面上。需要说明的是，在二维复平面上，视距反射信号作为本申请实施例的目标提取信号，随待测振动源的振动绕其起始点做简谐旋转运动，其中，在简谐旋转运动过程中，旋转幅角包含了待测振动源的振幅信息，旋转速率包含了待测振动源的频率信息；目标相关多径反射信号为干扰信号，和待测振动源不相关，其中，对于相同发送天线发送的毫米波雷达信号，由不同接收天线接收到的目标相关多径反射信号是绕其起始点的旋转；环境多径反射信号也是干扰信号，并且与待测振动源以及接收天线均不相关，在二维复平面上属于非运动向量。In the embodiment of the present invention, it is necessary to separate and remove the multipath reflection interference signal from the target composite wireless signal, so as to obtain a single line-of-sight reflection signal, wherein the multipath reflection interference signal includes target-related multipath reflection signals and environmental multipath reflection signals. reflected signal. Figure 3 is a schematic diagram of submillimeter vibration signal extraction provided by the embodiment of the present invention, as shown in Figure 3, in the embodiment of the present invention, the line-of-sight reflection signal, the target-related multipath reflection signal and the environment multipath reflection signal are superimposed Together, the composite signal generated by the superposition of the three is expressed as a vector addition of vectors on the two-dimensional complex plane, so that the target composite wireless signal is projected onto the two-dimensional complex plane. It should be noted that on the two-dimensional complex plane, the line-of-sight reflection signal, as the target extraction signal of the embodiment of the present application, performs a simple harmonic rotational motion around its starting point with the vibration of the vibration source to be measured, wherein, in the simple harmonic rotational motion During the process, the rotation angle contains the amplitude information of the vibration source to be measured, and the rotation rate contains the frequency information of the vibration source to be measured; the target-related multipath reflection signal is an interference signal, which is not related to the vibration source to be measured. Among them, for the same For the millimeter-wave radar signal sent by the transmitting antenna, the target-related multipath reflection signals received by different receiving antennas rotate around its starting point; the environmental multipath reflection signal is also an interference signal, and it is different from the vibration source to be measured and the receiving antenna. Correlation, in the two-dimensional complex plane belongs to the non-motion vector.
进一步地，在本发明实施例中，对提取视距反射信号的方法进行说明，在去除不同传播距离的信号之后，假设待测振动源的振动信号为a sin(2πft)，则对应的传播距离可表示为R(t)＝R0+a sin(2πft)，其中，R0表示待测振动源和毫米波雷达信号发送天线之间的测量距离。在理想情况下，接收天线仅接收到包含待测振动源信息的信号Y，可以理解的是，在本发明实施例中，信号Y可作为视觉反射信号进行说明，即视距反射信号然后，计算视距反射信号Y的相位角序列，并通过低通滤波直流分量对相位角序列进行处理，得到视距反射信号的时域表示a sin(2πft)，最后，通过最大似然参数估计方法多时域表示进行估算，从而得到待测振动源的振幅a和频率f的值，以对待测振动源进行无线监测。进一步地，上述对视距反射信号的提取方法属于理想情况，在实际的亚毫米级振动的无线监测环境中，接收天线接收到的信号Y并不是单路信号，而是包含多种干扰信号的复合信号，其包含了上述实施例的目标相关多径反射信号和环境多径反射信号。因此，需要将这两类干扰信号去除，可参考图3所示，当发送-接收天线为一组时，对于同一接收天线，目标相关多径反射信号和环境多径反射信号是非运动向量，因此，由视觉反射信号、目标相关多径反射信号和环境叠加产生的复合信号，在二维复平面上的运动是根据视距反射信号的简谐旋转运动形成的。需要说明的是，圆弧拟合算法包括几何圆弧拟合算法和代数圆弧拟合算法，其中，几何圆弧拟合算法对短圆弧的拟合效果更好，代数圆弧拟合算法的运算速率更快，优选地，在本发明实施例中，采用几何圆弧拟合算法的文伯格-马夸特(Levenberg-Marquardt，简称LM)优化算法进行圆弧拟合。通过圆弧几何拟合算法，获取视距反射信号的旋转中心，并通过坐标轴平移的方法消除目标相关多径反射信号和环境多径反射信号的干扰。在本发明实施例，将上述实施例关于Y(ω*)的连续采样作为信号Y进行说明，即Y＝Y(ω*)，具体步骤如下：Further, in the embodiment of the present invention, the method of extracting the line-of-sight reflection signal is described. After removing the signals of different propagation distances, assuming that the vibration signal of the vibration source to be measured is a sin (2πft), the corresponding propagation distance It can be expressed as R(t)=R 0 +a sin(2πft), where R 0 represents the measurement distance between the vibration source to be measured and the millimeter-wave radar signal transmitting antenna. In an ideal situation, the receiving antenna only receives the signal Y containing the information of the vibration source to be measured. It can be understood that in the embodiment of the present invention, the signal Y can be described as a visual reflection signal, that is, a line-of-sight reflection signal Then, calculate the phase angle sequence of the line-of-sight reflection signal Y, and process the phase angle sequence by low-pass filtering the DC component to obtain the time-domain representation of the line-of-sight reflection signal a sin(2πft), and finally, estimate it through the maximum likelihood parameter The method uses multi-time domain representation to estimate, so as to obtain the values of amplitude a and frequency f of the vibration source to be measured, so as to carry out wireless monitoring of the vibration source to be measured. Furthermore, the above-mentioned method for extracting the line-of-sight reflection signal belongs to the ideal situation. In the actual wireless monitoring environment of submillimeter level vibration, the signal Y received by the receiving antenna is not a single-channel signal, but contains a variety of interference signals. The composite signal includes the target-related multipath reflection signal and the environment multipath reflection signal in the above embodiment. Therefore, these two types of interference signals need to be removed, as shown in Figure 3, when the transmitting-receiving antennas are a group, for the same receiving antenna, the target-related multipath reflection signal and the environment multipath reflection signal are non-motion vectors, so , the composite signal generated by superposition of visual reflection signal, target-related multipath reflection signal and environment, the motion on the two-dimensional complex plane is formed according to the simple harmonic rotational motion of the line-of-sight reflection signal. It should be noted that the arc fitting algorithm includes a geometric arc fitting algorithm and an algebraic arc fitting algorithm, wherein the geometric arc fitting algorithm has a better fitting effect on short arcs, and the algebraic arc fitting algorithm The calculation speed is faster. Preferably, in the embodiment of the present invention, the arc fitting is performed by using the Levenberg-Marquardt (LM) optimization algorithm of the geometric arc fitting algorithm. Through the arc geometric fitting algorithm, the rotation center of the line-of-sight reflection signal is obtained, and the interference of the target-related multipath reflection signal and the environment multipath reflection signal is eliminated by the method of coordinate axis translation. In the embodiment of the present invention, the continuous sampling of Y(ω*) in the above embodiment is described as the signal Y, that is, Y=Y(ω*), and the specific steps are as follows:
步骤01，获取目标复合无线信号Y(ω*)的实数部分Re(Y)和虚数部分Im(Y)，并以二维坐标形式将实数部分和虚数部分投影到二维复平面上，并在二维复平面上绘制出目标复合无线信号对应圆弧；Step 01, obtain the real part Re(Y) and the imaginary part Im(Y) of the target composite wireless signal Y(ω*), and project the real part and the imaginary part onto the two-dimensional complex plane in the form of two-dimensional coordinates, and Draw the arc corresponding to the target composite wireless signal on the two-dimensional complex plane;
步骤02，根据LM优化算法的圆弧几何拟合算法，计算目标复合无线信号对应圆弧所在圆的参数，即所在圆的圆心坐标(x0，y0)与半径r；Step 02, according to the arc geometry fitting algorithm of the LM optimization algorithm, calculate the parameters of the circle where the target composite wireless signal corresponds to the arc, that is, the center coordinates (x 0 , y 0 ) and radius r of the circle;
步骤03，基于上述步骤得到的圆心坐标(x0，y0)，通过坐标轴平移的方法，去除目标相关多径反射信号和环境多径反射信号，从而得到目标单源信号，即视距反射信号Y′＝[Re(Y)-x0]+j[Im(Y)-y0]；Step 03, based on the coordinates of the center of the circle (x 0 , y 0 ) obtained in the above steps, the target-related multipath reflection signal and the environmental multipath reflection signal are removed by the method of coordinate axis translation, so as to obtain the target single-source signal, that is, the line-of-sight reflection Signal Y'=[Re(Y)-x 0 ]+j[Im(Y)-y 0 ];
步骤04，获取视距反射信号Y＇的相位角序列；Step 04, obtaining the phase angle sequence of the line-of-sight reflection signal Y';
步骤05，通过对视距反射信号Y＇的相位角序列进行低通滤波过滤直流分量处理，得到视距反射信号Y＇对应的时域表示，即得到关于待测振动源振动信号的时域表示a sin(2πft)；Step 05, by low-pass filtering the DC component processing on the phase angle sequence of the line-of-sight reflection signal Y', the time-domain representation corresponding to the line-of-sight reflection signal Y' is obtained, that is, the time-domain representation of the vibration signal of the vibration source to be measured is obtained a sin(2πft);
步骤06，通过最大似然参数估计方法，对视距反射信号Y＇的时域表示进行估算，从而得到待测振动源中振动信号的振幅a和频率f，以对待测振动源进行无线监测。Step 06: Estimate the time-domain representation of the line-of-sight reflection signal Y' through the maximum likelihood parameter estimation method, so as to obtain the amplitude a and frequency f of the vibration signal in the vibration source to be measured, so as to conduct wireless monitoring of the vibration source to be measured.
在上述实施例的基础上，所述对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号，所述方法还包括：On the basis of the above-mentioned embodiments, the extraction processing of the target composite wireless signal is carried out to obtain the line-of-sight reflection signal of the vibration source to be measured, and the method further includes:
S1，若发送-接收天线为多组时，获取多个目标复合无线信号；S1, if there are multiple sets of sending-receiving antennas, acquiring multiple target composite wireless signals;
S2，将多个目标复合无线信号投射到二维复平面，并得到每个目标复合无线信号对应的圆弧；根据圆弧几何拟合算法，获取二维复平面上每个目标复合无线信号对应圆弧所在圆的圆心坐标和半径，并根据圆心坐标拟合得到多个目标复合无线信号的共圆参数，并将所述共圆参数作为半圆约束条件，所述共圆参数包括共圆圆心坐标和半径；S2. Project multiple target composite wireless signals onto the two-dimensional complex plane, and obtain the arc corresponding to each target composite wireless signal; according to the circular arc geometric fitting algorithm, obtain the corresponding arc of each target composite wireless signal on the two-dimensional complex plane. The center coordinates and radius of the circle where the arc is located, and according to the center coordinates, the co-circle parameters of multiple target composite wireless signals are obtained, and the co-circle parameters are used as semicircle constraints, and the co-circle parameters include the co-circle center coordinates and radius;
S3，根据半径约束条件和所述共圆约束条件，通过圆弧几何拟合算法，获取在约束条件向下每个目标复合无线信号对应圆弧所在圆的圆心参数和半径，以获取新的共圆参数，若新的共圆参数满足预设阈值，则通过坐标轴平移，将每个目标复合无线信号中多径反射干扰信号去除，以得到每个接收天线对应的所述待测振动源的视距反射信号；若新的共圆参数不满足预设阈值，则将新的共圆参数作为下一次迭代过程的共圆约束条件，再次执行步骤S3，直至迭代得到的共圆参数满足预设阈值。S3. According to the radius constraint condition and the co-circle constraint condition, through the arc geometric fitting algorithm, obtain the center parameter and radius of the circle where each target composite wireless signal corresponds to the arc under the constraint condition, so as to obtain a new co-circle constraint condition circle parameter, if the new co-circle parameter satisfies the preset threshold, the multipath reflection interference signal in each target composite wireless signal is removed by coordinate axis translation, so as to obtain the vibration source to be measured corresponding to each receiving antenna Line-of-sight reflection signal; if the new co-circle parameter does not meet the preset threshold, the new co-circle parameter is used as the co-circle constraint condition for the next iteration process, and step S3 is executed again until the co-circle parameter obtained by iteration meets the preset threshold.
在上述实施例的基础上，所述根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测，还包括：On the basis of the above embodiments, the acquisition of the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be measured, further includes:
获取每个接收天线对应的视距反射信号的相位角序列；Obtain the phase angle sequence of the line-of-sight reflection signal corresponding to each receiving antenna;
对每个相位角序列进行低通滤波直流分量处理，得到每个视距反射信号的时域表示；Perform low-pass filtering DC component processing on each phase angle sequence to obtain the time domain representation of each line-of-sight reflection signal;
根据四分位均值数据融合算法，对每个视距反射信号的时域表示进行融合，得到融合后的时域表示；According to the quartile mean data fusion algorithm, the time domain representation of each line-of-sight reflection signal is fused to obtain the fused time domain representation;
根据最大似然参数估计方法，对融合后的时域表示进行估算，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。According to the maximum likelihood parameter estimation method, the fused time domain representation is estimated, and the amplitude and frequency of the vibration source to be measured are obtained, so as to wirelessly monitor the vibration source to be measured.
在对单组发送-接收天线的信号进行提取时，得到的待测振动源对应的振动信号原始波形以及参数，还存在一定的误差，其误差的主要来源为圆弧拟合。优选地，在本发明实施例中，为了进一步降低信号提取误差，通过设置多组发送-接收天线，对多天线的设置来限制圆弧拟合的搜索空间，相比上述实施例提供的单组天线的步骤，本发明实施例引入两个约束条件：半径约束条件和共圆约束条件。其中，半径约束条件用于限定半径的范围，半径长度和天线-物体距离以及天线增益相关，在本发明实施例中，对于每一个发送-接收天线组成的对(排除天线增益影响)，通过多次测量统计出不同距离下半径的长度，半径约束条件为误差20％以内，例如，统计出天线至物体的距离为x米，则半径长度为y，半径约束条件为y×0.8～y×1.2，可参考图3所示，二维复平面中视距反射向量的长度是固定的，由物体到天线距离以及天线雷达散射截面积确定的，不受振动源的振动信号影响，因此，对应每一组发送-接收天线，可以统计出不同距离反射信号的视距反射向量的长度，从而在小范围内通过半径约束条件优化圆弧拟合过程。另外，视距反射向量的长度和目标相关多径反射向量的长度，均由物体到天线距离以及天线雷达散射截面积确定，由于在多组天线情况下，不同天线组的目标相关多径反射向量具有共圆性质，因此，通过将两个向量长度相除，从而排除该组天线被其他天线组的影响，得到一个与天线无关的约束比例参数，用于优化圆弧拟合过程，即共圆约束条件，需要说明的是，对于多组天线的情况(天线组数量大于两个时)，将约束比例参数组成的序列方差作为共圆约束条件，从而优化圆弧拟合过程。When extracting the signal of a single set of transmitting-receiving antennas, there are still certain errors in the original waveform and parameters of the vibration signal corresponding to the vibration source to be tested, and the main source of the error is arc fitting. Preferably, in the embodiment of the present invention, in order to further reduce the signal extraction error, by setting multiple sets of transmitting-receiving antennas, the setting of multiple antennas limits the search space of arc fitting, compared to the single set of antennas provided by the above-mentioned embodiment For the steps of the antenna, the embodiment of the present invention introduces two constraint conditions: a radius constraint condition and a co-circle constraint condition. Wherein, the radius constraint condition is used to limit the scope of the radius, and the radius length is related to the antenna-object distance and the antenna gain. The length of the radius at different distances can be counted for each measurement, and the radius constraint condition is within 20% of the error. For example, if the distance from the antenna to the object is calculated to be x meters, the radius length is y, and the radius constraint condition is y×0.8～y×1.2 , as shown in Figure 3, the length of the line-of-sight reflection vector in the two-dimensional complex plane is fixed, which is determined by the distance from the object to the antenna and the radar cross-sectional area of the antenna, and is not affected by the vibration signal of the vibration source. Therefore, corresponding to each A group of sending-receiving antennas can count the length of the line-of-sight reflection vector of different distance reflection signals, so as to optimize the arc fitting process in a small range through radius constraints. In addition, the length of the line-of-sight reflection vector and the length of the target-related multipath reflection vector are both determined by the distance from the object to the antenna and the radar cross-sectional area of the antenna. Because in the case of multiple antennas, the target-related multipath reflection vector of different antenna groups It has the property of co-circle, therefore, by dividing the length of the two vectors, the group of antennas is excluded from the influence of other antenna groups, and a constraint ratio parameter independent of the antenna is obtained, which is used to optimize the arc fitting process, that is, co-circle Constraint conditions, it should be noted that for the case of multiple antenna groups (when the number of antenna groups is greater than two), the sequence variance composed of constraint ratio parameters is used as the co-circle constraint condition to optimize the arc fitting process.
进一步地，在本发明实施例中，对于多组天线组情况下，信号提取的具体步骤如下：Further, in the embodiment of the present invention, in the case of multiple antenna groups, the specific steps of signal extraction are as follows:
步骤S11，获取每组接收天线接收到的对应目标复合无线信号；Step S11, obtaining the corresponding target composite wireless signal received by each group of receiving antennas;
步骤S12，获取每个目标复合无线信号的圆弧参数，具体步骤可参考根据上述实施例中步骤01至步骤02的方法，即将多个目标复合无线信号投射到二维复平面，并得到每个目标复合无线信号对应的圆弧；根据LM优化算法的圆弧几何拟合算法，获取二维复平面上每个目标复合无线信号对应圆弧所在圆的圆心坐标和半径；再根据圆心坐标拟合得到多个目标复合无线信号的共圆参数，并将该共圆参数作为半圆约束条件；Step S12, obtain the arc parameters of each target composite wireless signal, the specific steps can refer to the method according to step 01 to step 02 in the above embodiment, that is, project multiple target composite wireless signals onto the two-dimensional complex plane, and obtain each The arc corresponding to the target composite wireless signal; according to the circular arc geometric fitting algorithm of the LM optimization algorithm, obtain the center coordinates and radius of the circle where each target composite wireless signal corresponds to the arc on the two-dimensional complex plane; then fit according to the center coordinates Obtaining the co-circle parameters of multiple target composite wireless signals, and using the co-circle parameters as semicircle constraints;
步骤S13，在LM优化算法的圆弧几何拟合算法中添加半径约束条件和共圆约束条件，再次对每组接收天线接收到目标复合无线信号进行圆弧几何拟合，得到每个目标复合无线信号对应圆弧所在圆的圆心坐标和半径，再根据圆心坐标拟合得到新的共圆参数，并对新的共圆参数的收敛情况进行判断，在本发明实施例中，通过设置预设阈值判断是否收敛，若收敛，则进行下一步骤，若未收敛，则将本次得到的新的共圆参数作为下一次迭代过程的共圆参数，再次执行步骤S13，直至迭代得到的共圆参数满足预设阈值；Step S13, adding radius constraints and co-circle constraints to the circular arc geometric fitting algorithm of the LM optimization algorithm, and performing circular arc geometric fitting on the target composite wireless signals received by each group of receiving antennas again, to obtain each target composite wireless signal The signal corresponds to the center coordinates and radius of the circle where the arc is located, and then the new co-circle parameters are obtained by fitting according to the center coordinates, and the convergence of the new co-circle parameters is judged. In the embodiment of the present invention, by setting the preset threshold Judging whether it is converged, if converged, proceed to the next step, if not converged, use the new co-circle parameter obtained this time as the co-circle parameter of the next iteration process, and execute step S13 again until the co-circle parameter obtained by iteration Satisfy the preset threshold;
步骤S14，在确定共圆参数满足预设阈值之后，通过步骤S12中得到的每个目标复合无线信号对应圆弧所在圆的圆心坐标，获取每个视距反射信号对应的时域表示，即得到关于待测振动源振动信号的时域表示，具体步骤可参考根据上述实施例中步骤03至步骤05的方法，本发明实施例不再赘述；Step S14, after determining that the co-circle parameter satisfies the preset threshold, the time-domain representation corresponding to each line-of-sight reflection signal is obtained through the coordinates of the center of the circle where each target composite wireless signal corresponds to the arc obtained in step S12, that is, Regarding the time-domain representation of the vibration signal of the vibration source to be measured, the specific steps can refer to the method according to the steps 03 to 05 in the above-mentioned embodiments, and the embodiments of the present invention will not be repeated;
步骤S15，根据四分位均值数据融合算法(四分之一位数和四分之三位数的均值)，对每个视距反射信号的时域表示进行融合，得到融合后的时域表示；Step S15, according to the quartile mean data fusion algorithm (the mean value of the 1/4 digit and the 3/4 digit), the time domain representation of each line-of-sight reflection signal is fused to obtain the fused time domain representation ;
步骤S16，根据最大似然参数估计方法，对融合后的时域表示进行估算，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。Step S16 , estimating the fused time domain representation according to the maximum likelihood parameter estimation method, and acquiring the amplitude and frequency of the vibration source to be measured, so as to wirelessly monitor the vibration source to be measured.
为了验证上述各实施例的有效性，在本发明实施例中，通过使用一台振幅和频率连续可调的高精度振动传感器校准仪(频率控制范围1-1000Hz，精度1Hz，误差约±1Hz；振幅控制范围15-500um，精度1um，误差约±5um)进行评估。在本发明实施例中，分别选取测量距离为0.5m、1m、2m和3m，振动频率为20Hz、50Hz、100Hz、200Hz和300Hz，振幅为15um、30um、50um、100um、200um和400um的测试数据进行验证说明，其中，与本发明实施例进行对比的方法，是直接求取相位角提取的理论方法，与对比方法相比，本发明实施例提供的亚毫米级振动的无线监测方法，主要是对振动估计的精准度进行提升，因此，在频率估计上，两者之间的表现情况基本一致。In order to verify the effectiveness of the above-mentioned embodiments, in the embodiment of the present invention, by using a high-precision vibration sensor calibrator with continuously adjustable amplitude and frequency (frequency control range 1-1000Hz, accuracy 1Hz, error about ± 1Hz; The amplitude control range is 15-500um, the precision is 1um, and the error is about ±5um) for evaluation. In the embodiment of the present invention, the test data with measuring distances of 0.5m, 1m, 2m and 3m, vibration frequencies of 20Hz, 50Hz, 100Hz, 200Hz and 300Hz, and amplitudes of 15um, 30um, 50um, 100um, 200um and 400um were respectively selected. Verification instructions are given, wherein, the method for comparison with the embodiment of the present invention is a theoretical method for directly obtaining the phase angle extraction. Compared with the comparison method, the wireless monitoring method for submillimeter vibration provided by the embodiment of the present invention is mainly The accuracy of vibration estimation is improved, so in terms of frequency estimation, the performance between the two is basically the same.
具体地，图4为本发明实施例提供的振幅与频率的测量误差示意图，可参考图4所示，在本发明实施例中，在测量距离为100cm，振动频率为50Hz的条件下，当振幅较小时，本发明实施例的振幅测量误差小于10um；当振幅较大时，本发明实施例的振幅测量相对误差约10％，相比对比方法，本发明实施例的振幅测量误差减少了10倍左右。另外，本发明实施例和对比方法的频率测量的误差均较小，误差均小于1Hz。Specifically, Fig. 4 is a schematic diagram of the measurement error of amplitude and frequency provided by the embodiment of the present invention, as shown in Fig. 4 , in the embodiment of the present invention, under the condition that the measurement distance is 100 cm and the vibration frequency is 50 Hz, when the amplitude When it is small, the amplitude measurement error of the embodiment of the present invention is less than 10um; when the amplitude is large, the relative error of the amplitude measurement of the embodiment of the present invention is about 10%. Compared with the comparison method, the amplitude measurement error of the embodiment of the present invention is reduced by 10 times about. In addition, the frequency measurement errors of the embodiment of the present invention and the comparison method are both small, and the errors are both less than 1 Hz.
图5为本发明另一实施例提供的振幅与频率的测量误差示意图，可参考图5所示，在测量距离为100cm，振幅为15um的情况下，当振动频率较高时，本发明实施例的振幅测量误差更小，90％的振幅测量误差小于5um；当振动频率较低时，由于振动校准仪输出小频率振动会出现相对误差较大的情况，导致本发明实施例体现出振幅测量误差较大的现象。而对比方法的振幅测量误差虽然更加恒定且集中，但其误差是本发明实施例的2倍左右。另外，本发明实施例和对比方法的频率测量的误差同样较小，误差均小于5Hz。Fig. 5 is a schematic diagram of the measurement error of amplitude and frequency provided by another embodiment of the present invention, as shown in Fig. 5, in the case where the measurement distance is 100cm and the amplitude is 15um, when the vibration frequency is high, the embodiment of the present invention The amplitude measurement error is smaller, and 90% of the amplitude measurement error is less than 5um; when the vibration frequency is low, the relative error will be relatively large due to the small frequency vibration output by the vibration calibrator, resulting in the embodiment of the present invention reflecting the amplitude measurement error larger phenomenon. Although the amplitude measurement error of the comparative method is more constant and concentrated, its error is about twice that of the embodiment of the present invention. In addition, the frequency measurement errors of the embodiment of the present invention and the comparative method are also small, and the errors are both less than 5 Hz.
图6为本发明又一实施例提供的振幅与频率的测量误差示意图，可参考图6所示，在振幅为100um，振动频率为50Hz的情况下，当测量距离较短时，本发明实施例相比对比方法，振幅测量误差更小，特别在小于200cm的测量距离时，90％的振幅测量误差小于10um(10％相对误差)；当测量距离较长时，例如，测量距离为300cm时，由于信噪比更低，导致本发明实施例的振幅测量误差稍大，但相比对比方法，还是具有较低的误差，如图6所示，对比方法的整体振幅测量误差为本发明实施例的8倍左右，可以看出，对比方法的振幅测量误差非常恒定，测量距离50cm至300cm的振幅测量误差非常接近。另外，本发明实施例和对比方法的频率测量误差同样较小，两个方法在测量距离50cm至300cm的情况下，频率测量误差几乎相等，均小于1Hz，因此在图6中呈现出一条竖线的情况。Fig. 6 is a schematic diagram of the measurement error of the amplitude and frequency provided by another embodiment of the present invention, as shown in Fig. 6, when the amplitude is 100um and the vibration frequency is 50Hz, when the measurement distance is relatively short, the embodiment of the present invention Compared with the comparative method, the amplitude measurement error is smaller, especially when the measurement distance is less than 200cm, 90% of the amplitude measurement error is less than 10um (10% relative error); when the measurement distance is longer, for example, when the measurement distance is 300cm, Due to the lower signal-to-noise ratio, the amplitude measurement error of the embodiment of the present invention is slightly larger, but compared with the comparative method, it still has a lower error. As shown in Figure 6, the overall amplitude measurement error of the comparative method is the embodiment of the present invention. It can be seen that the amplitude measurement error of the comparison method is very constant, and the amplitude measurement error of the measurement distance 50cm to 300cm is very close. In addition, the frequency measurement errors of the embodiment of the present invention and the comparison method are also small. When the measurement distance is 50cm to 300cm, the frequency measurement errors of the two methods are almost equal, both less than 1Hz, so a vertical line appears in Figure 6 Case.
由上述的实施例的评估测试可知，本发明各实施例提供的方法，以无线且非侵入式(2m内的监测数据较准)监测机械振动频率低至15um振幅的微小振动，其振幅测量误差约10um，相对测量误差约为10％，相比对比方法，监测准确度效果提高了5倍左右，在毫米级以及亚毫米级的振动监测下具有显著的效果。From the evaluation tests of the above-mentioned embodiments, it can be seen that the methods provided by each embodiment of the present invention monitor the mechanical vibration frequency as low as 15um amplitude in a wireless and non-invasive way (the monitoring data within 2m is accurate), and the amplitude measurement error It is about 10um, and the relative measurement error is about 10%. Compared with the comparison method, the monitoring accuracy effect is increased by about 5 times, and it has a significant effect in the vibration monitoring of millimeter and submillimeter levels.
图7为本发明实施例提供的基于毫米波的振动监测系统的结构示意图，如图7所示，本发明实施例提供了一种基于毫米波的振动监测系统，包括复合信号获取模块701、视距反射信号提取模块702和监测模块703，其中，复合信号获取模块701用于获取待测振动源的目标复合无线信号；视距反射信号提取模块702用于对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；监测模块703用于根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。Fig. 7 is a schematic structural diagram of a vibration monitoring system based on millimeter waves provided by an embodiment of the present invention. As shown in Fig. 7, an embodiment of the present invention provides a vibration monitoring system based on millimeter waves, including a composite signal acquisition module 701, a Distance reflection signal extraction module 702 and monitoring module 703, wherein, composite signal acquisition module 701 is used for obtaining the target composite wireless signal of vibration source to be measured; Line-of-sight reflection signal extraction module 702 is used for extracting and processing described target composite wireless signal , to obtain the line-of-sight reflection signal of the vibration source to be measured; the monitoring module 703 is used to obtain the amplitude and frequency of the vibration source to be measured according to the line-of-sight reflection signal, so as to wirelessly monitor the vibration source to be measured .
在本发明实施例中，首先，基于毫米波的振动监测系统通过发送天线，向待测振动源发送毫米波雷达信号，并通过接收天线接收到待测振动源以及周围物体反射的毫米波雷达信号。由于微小振动是一个往复运动，会周期性地改变无线信号的传播路径长度，从而改变无线信号的相位；并且，由于毫米波频段的无线信号波长约1-10mm，因此可以通过分析其相位变化提取出对应的振动信号。然而，由无线信号在空间传播的特性可知，发射出去的毫米波无线信号，会被周围所有的物体反射，导致接收到的无线信号是一个复合信号，该复合信号中包含了多路不同传播距离的信号，同时，在该复合信号中也包含了由多条传播路径产生的相同传播距离的信号。在本发明实施例中，复合信号获取模块701对待测振动源和周围物体反射的毫米波雷达信号进行信号提取，将不同传播距离的干扰信号去除，从而提取得到相同传播距离不同传播路径的复合信号，即待测振动源的目标复合无线信号。由于目标复合无线信号中携带有视距反射信号和多径反射干扰信号，因此，视距反射信号提取模块702需要将多径反射干扰信号去除，从而得到待测振动源在视距路径上的反射信号，即待测振动源的视距反射信号。最后，监测模块703计算得到该视距反射信号的相位角序列信息，并对相位角序列信息进行低通滤波过滤直流分量处理，得到待测振动源对应振动信号的时域表示，并通过对振动信号的时域表示进行最大似然参数估计，从而得到该振动信号的振幅和频率，以用于对待测振动源进行分析，得到对应的振动监测数据。In the embodiment of the present invention, first, the millimeter-wave-based vibration monitoring system sends millimeter-wave radar signals to the vibration source to be measured through the transmitting antenna, and receives the millimeter-wave radar signals reflected by the vibration source to be measured and surrounding objects through the receiving antenna . Since the tiny vibration is a reciprocating motion, it will periodically change the propagation path length of the wireless signal, thereby changing the phase of the wireless signal; and, since the wavelength of the wireless signal in the millimeter wave band is about 1-10mm, it can be extracted by analyzing its phase change. output the corresponding vibration signal. However, according to the characteristics of wireless signal propagation in space, the transmitted millimeter wave wireless signal will be reflected by all surrounding objects, resulting in the received wireless signal being a composite signal, which contains multiple channels with different propagation distances. At the same time, the composite signal also includes signals with the same propagation distance generated by multiple propagation paths. In the embodiment of the present invention, the composite signal acquisition module 701 performs signal extraction on millimeter-wave radar signals reflected by the vibration source to be measured and surrounding objects, and removes interference signals with different propagation distances, thereby extracting composite signals with the same propagation distance and different propagation paths , that is, the target composite wireless signal of the vibration source to be tested. Since the target composite wireless signal carries line-of-sight reflection signals and multi-path reflection interference signals, the line-of-sight reflection signal extraction module 702 needs to remove the multi-path reflection interference signals, so as to obtain the reflection of the vibration source to be measured on the line-of-sight path The signal is the line-of-sight reflection signal of the vibration source to be measured. Finally, the monitoring module 703 calculates the phase angle sequence information of the line-of-sight reflection signal, and performs low-pass filtering on the phase angle sequence information to filter the DC component to obtain the time domain representation of the vibration signal corresponding to the vibration source to be measured, and through the vibration The time domain representation of the signal is used for maximum likelihood parameter estimation to obtain the amplitude and frequency of the vibration signal, which can be used to analyze the vibration source to be measured and obtain the corresponding vibration monitoring data.
本发明实施例提供的一种基于毫米波的振动监测系统，通过将待测振动源及周围物体发送的毫米波雷达信号中各种干扰信号去除，得到包含待测振动源振动信号相关信息的视距反射信号，从而得到待测振动源的振幅和频率，以用于对待测振动源进行振动监测，提高了振动监测的精准度，大幅度降低了振动监测的误差。A millimeter-wave-based vibration monitoring system provided by an embodiment of the present invention removes various interference signals from the millimeter-wave radar signals sent by the vibration source to be measured and surrounding objects, and obtains a visual image containing information related to the vibration signal of the vibration source to be measured. The amplitude and frequency of the vibration source to be measured can be obtained by reflecting the signal from the distance, which can be used for vibration monitoring of the vibration source to be measured, which improves the accuracy of vibration monitoring and greatly reduces the error of vibration monitoring.
在上述实施例的基础上，所述系统还包括第一处理模块和第二处理模块，其中，第一处理模块用于获取待测振动源和周围物体反射的毫米波雷达信号；第二处理模块用于对所述毫米波雷达信号进行提取处理，获取相同传播距离且传播路径不同的复合信号，以得到所述待测振动源的目标复合无线信号。On the basis of the above-mentioned embodiments, the system further includes a first processing module and a second processing module, wherein the first processing module is used to obtain millimeter-wave radar signals reflected by the vibration source to be measured and surrounding objects; the second processing module It is used for extracting and processing the millimeter-wave radar signal to obtain composite signals with the same propagation distance and different propagation paths, so as to obtain the target composite wireless signal of the vibration source to be measured.
本发明实施例提供的系统是用于执行上述各方法实施例的，具体流程和详细内容请参照上述实施例，此处不再赘述。The system provided by the embodiments of the present invention is used to execute the above-mentioned method embodiments. Please refer to the above-mentioned embodiments for specific procedures and details, and details will not be repeated here.
图8为本发明实施例提供的电子设备结构示意图，参照图8，该电子设备可以包括：处理器(processor)801、通信接口(Communications Interface)802、存储器(memory)803和通信总线804，其中，处理器801，通信接口802，存储器803通过通信总线804完成相互间的通信。处理器801可以调用存储器803中的逻辑指令，以执行如下方法：获取待测振动源的目标复合无线信号；对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。FIG. 8 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention. Referring to FIG. 8, the electronic device may include: a processor (processor) 801, a communication interface (Communications Interface) 802, a memory (memory) 803, and a communication bus 804, wherein , the processor 801 , the communication interface 802 , and the memory 803 communicate with each other through the communication bus 804 . The processor 801 can call the logic instructions in the memory 803 to perform the following method: obtain the target composite wireless signal of the vibration source to be tested; perform extraction processing on the target composite wireless signal, and obtain the line-of-sight reflection of the vibration source to be measured signal; according to the line-of-sight reflection signal, the amplitude and frequency of the vibration source to be tested are obtained, so as to wirelessly monitor the vibration source to be tested.
此外，上述的存储器803中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。In addition, the above logic instructions in the memory 803 may be implemented in the form of software function units and when sold or used as an independent product, may be stored in a computer-readable storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes. .
另一方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的基于毫米波的振动监测方法，例如包括：获取待测振动源的目标复合无线信号；对所述目标复合无线信号进行提取处理，获取所述待测振动源的视距反射信号；根据所述视距反射信号，获取所述待测振动源的振幅和频率，以对所述待测振动源进行无线监测。On the other hand, an embodiment of the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it is implemented to perform the millimeter-wave-based vibration provided by the above-mentioned embodiments. The monitoring method, for example, includes: acquiring the target composite wireless signal of the vibration source to be tested; extracting and processing the target composite wireless signal, and obtaining the line-of-sight reflection signal of the vibration source to be measured; according to the line-of-sight reflection signal, obtaining The amplitude and frequency of the vibration source to be tested are used for wireless monitoring of the vibration source to be tested.
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。The device embodiments described above are only illustrative, and the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in One place, or it can be distributed to multiple network elements. Part or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. It can be understood and implemented by those skilled in the art without any creative efforts.
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。Through the above description of the implementations, those skilled in the art can clearly understand that each implementation can be implemented by means of software plus a necessary general-purpose hardware platform, and of course also by hardware. Based on this understanding, the essence of the above technical solution or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic discs, optical discs, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute the methods described in various embodiments or some parts of the embodiments.
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still be Modifications are made to the technical solutions described in the foregoing embodiments, or equivalent replacements are made to some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present invention.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910530699.5A CN110220586A (en) | 2019-06-19 | 2019-06-19 | Vibration monitoring method and system based on millimeter wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910530699.5A CN110220586A (en) | 2019-06-19 | 2019-06-19 | Vibration monitoring method and system based on millimeter wave |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110220586A true CN110220586A (en) | 2019-09-10 |
Family
ID=67817805
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910530699.5A Pending CN110220586A (en) | 2019-06-19 | 2019-06-19 | Vibration monitoring method and system based on millimeter wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110220586A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110995626A (en) * | 2019-10-11 | 2020-04-10 | 清华大学 | Concurrent channel parameter estimation method, system, electronic device and storage medium |
| CN111609920A (en) * | 2020-05-13 | 2020-09-01 | 上海交通大学 | Handheld Microwave Vibration Measurement System |
| CN111983570A (en) * | 2020-07-16 | 2020-11-24 | 北京宏锐星通科技有限公司 | Interference method for SAR-GMTI mode |
| CN114608699A (en) * | 2022-03-07 | 2022-06-10 | 中国矿业大学 | Elevator head sheave vibration abnormity detection system and method based on millimeter wave radar |
| CN114689161A (en) * | 2022-03-30 | 2022-07-01 | 珠海格力电器股份有限公司 | Detection method, nonvolatile storage medium, and detection system |
| CN114966321A (en) * | 2022-05-27 | 2022-08-30 | 桂林电子科技大学 | An inductive path reversal measuring device |
| CN115436938A (en) * | 2022-09-22 | 2022-12-06 | 山东大学 | Passive signal measuring device and measuring method based on millimeter wave radar |
| CN115598644A (en) * | 2022-11-28 | 2023-01-13 | 北京木牛领航科技有限公司（Cn） | Collision detection method, device, medium and equipment |
| CN115956192A (en) * | 2020-09-11 | 2023-04-11 | 日本电信电话株式会社 | Vibration detection device and vibration detection method |
| CN119197737A (en) * | 2024-09-19 | 2024-12-27 | 淮北矿业股份有限公司 | Non-contact vibration detection device and method for hoist motor based on millimeter wave radar |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101050977A (en) * | 2007-04-28 | 2007-10-10 | 浙江大学 | Method and system for measuring low frequency vibration using millimeter wave |
| US20150030256A1 (en) * | 2012-05-09 | 2015-01-29 | Duke University | Multi-sensor compressive imaging |
| CN106644030A (en) * | 2016-08-31 | 2017-05-10 | 上海交通大学 | Doppler radar-based non-contact type vibration measuring method |
| CN107607924A (en) * | 2017-08-21 | 2018-01-19 | 上海交通大学 | The signal processing method that a kind of fmcw radar static state noise jamming eliminates |
| CN109521404A (en) * | 2018-10-12 | 2019-03-26 | 上海交通大学 | The evaluation of accuracy and system of vibration measurement based on fmcw radar |
| CN109883535A (en) * | 2019-01-25 | 2019-06-14 | 陕煤集团神木张家峁矿业有限公司 | A kind of vibration measurement device and method based on millimeter wave |
-
2019
- 2019-06-19 CN CN201910530699.5A patent/CN110220586A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101050977A (en) * | 2007-04-28 | 2007-10-10 | 浙江大学 | Method and system for measuring low frequency vibration using millimeter wave |
| US20150030256A1 (en) * | 2012-05-09 | 2015-01-29 | Duke University | Multi-sensor compressive imaging |
| CN106644030A (en) * | 2016-08-31 | 2017-05-10 | 上海交通大学 | Doppler radar-based non-contact type vibration measuring method |
| CN107607924A (en) * | 2017-08-21 | 2018-01-19 | 上海交通大学 | The signal processing method that a kind of fmcw radar static state noise jamming eliminates |
| CN109521404A (en) * | 2018-10-12 | 2019-03-26 | 上海交通大学 | The evaluation of accuracy and system of vibration measurement based on fmcw radar |
| CN109883535A (en) * | 2019-01-25 | 2019-06-14 | 陕煤集团神木张家峁矿业有限公司 | A kind of vibration measurement device and method based on millimeter wave |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110995626A (en) * | 2019-10-11 | 2020-04-10 | 清华大学 | Concurrent channel parameter estimation method, system, electronic device and storage medium |
| CN111609920A (en) * | 2020-05-13 | 2020-09-01 | 上海交通大学 | Handheld Microwave Vibration Measurement System |
| CN111983570A (en) * | 2020-07-16 | 2020-11-24 | 北京宏锐星通科技有限公司 | Interference method for SAR-GMTI mode |
| CN115956192A (en) * | 2020-09-11 | 2023-04-11 | 日本电信电话株式会社 | Vibration detection device and vibration detection method |
| CN114608699A (en) * | 2022-03-07 | 2022-06-10 | 中国矿业大学 | Elevator head sheave vibration abnormity detection system and method based on millimeter wave radar |
| CN114608699B (en) * | 2022-03-07 | 2023-02-24 | 中国矿业大学 | Elevator head sheave vibration abnormity detection system and method based on millimeter wave radar |
| CN114689161A (en) * | 2022-03-30 | 2022-07-01 | 珠海格力电器股份有限公司 | Detection method, nonvolatile storage medium, and detection system |
| CN114966321A (en) * | 2022-05-27 | 2022-08-30 | 桂林电子科技大学 | An inductive path reversal measuring device |
| CN115436938A (en) * | 2022-09-22 | 2022-12-06 | 山东大学 | Passive signal measuring device and measuring method based on millimeter wave radar |
| CN115436938B (en) * | 2022-09-22 | 2024-05-07 | 山东大学 | A passive signal measurement device and measurement method based on millimeter wave radar |
| CN115598644A (en) * | 2022-11-28 | 2023-01-13 | 北京木牛领航科技有限公司（Cn） | Collision detection method, device, medium and equipment |
| CN119197737A (en) * | 2024-09-19 | 2024-12-27 | 淮北矿业股份有限公司 | Non-contact vibration detection device and method for hoist motor based on millimeter wave radar |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110220586A (en) | Vibration monitoring method and system based on millimeter wave | |
| Jiang et al. | mmVib: micrometer-level vibration measurement with mmwave radar | |
| Guo et al. | Measuring micrometer-level vibrations with mmWave radar | |
| KR101767924B1 (en) | Method and system for detecting location of multi-target | |
| CN102607936A (en) | Method and system for online creep monitoring | |
| CN111337919B (en) | Method for measuring antenna pattern of secondary radar and secondary radar | |
| US9491587B2 (en) | Spatially mapping radio frequency data in 3-dimensional environments | |
| CN112630771A (en) | Multi-target tracking method and system based on millimeter wave equipment | |
| KR101331833B1 (en) | Method for positioning using time difference of arrival | |
| CN114119465A (en) | Point cloud data processing method and device | |
| CN113009432B (en) | Method, device and equipment for improving measurement accuracy and target detection accuracy | |
| JPWO2013136648A1 (en) | Wide area multilateration system, central office, and two-dimensional position calculation method used therefor | |
| CN112731293A (en) | Non-contact sound and vibration combined detection system and detection method | |
| CN105068053A (en) | Method for extracting background signal from radar cross section (RCS) measurement data | |
| CN111123246B (en) | Method for testing maximum range finding capability and angle measurement precision of airborne laser radar based on circular scanning | |
| CN113009486A (en) | Human body sensing method and system based on millimeter wave radar | |
| US20110153269A1 (en) | Methods and apparatus for filtering noise in a three-dimensional position measurement system | |
| CN108169634B (en) | Method for accurately obtaining partial discharge ultrahigh frequency signal time difference | |
| CN115278524A (en) | Positioning method, device, equipment and storage medium of sensor network node | |
| CN115061102A (en) | Calibration method, calibration device, electronic equipment and storage medium | |
| Sinitsyn et al. | Passive acoustic radar for aircraft trajectory tracking | |
| CN118759452B (en) | A UWB-based intelligent mold opening positioning control system and method | |
| US12270902B2 (en) | Object detection apparatus, object detection method, and non-transitory computer readable-medium | |
| CN107071730B (en) | Indoor radio frequency positioning environment parameter real-time automatic learning method and system applied by same | |
| CN113514809B (en) | Vector miss distance parameter measurement method and device, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190910 |
|
| RJ01 | Rejection of invention patent application after publication |