CN104346528B - Effective vibration data interception method based on waveform characteristic statistics - Google Patents
Effective vibration data interception method based on waveform characteristic statistics Download PDFInfo
- Publication number
- CN104346528B CN104346528B CN201410568873.2A CN201410568873A CN104346528B CN 104346528 B CN104346528 B CN 104346528B CN 201410568873 A CN201410568873 A CN 201410568873A CN 104346528 B CN104346528 B CN 104346528B
- Authority
- CN
- China
- Prior art keywords
- dcs
- data
- value
- wave
- vibration data
- 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.)
- Active
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides an effective vibration data interception method based on waveform characteristic statistics. Filtering of a noise peak and interception of a vibration signal are realized by performing statistics on a waveform characteristic of acquired data based on a waveform characteristic of the vibration signal. With the adoption of the scheme, interference of a useless signal is effectively avoided, the operation times and a false alarm rate of a system are reduced, and the method is low in algorithm operand, high in timeliness, high in reliability and low in false alarm rate.
Description
Technical field
The invention belongs to effectively vibration data intercepts technical field, more particularly to a kind of counted based on wave character
Effective vibration data intercept method.
Background technology
MZ optical fiber perimeter intrusion detections instrument is based on Mach-Zehnder principle of interferences, is shaken using the optical fiber cable detection external world
Illegal invasion is moved and differentiates, electromagnetism interference is strong, be capable of achieving the long range real-time monitoring to the caused vibration of circumference invasion and fixed
Position.Used as the key technology of perimeter security system, the real-time and precise positioning of invasion is the main research of perimeter detection early warning system
One of direction.And in actual applications, signal Real-time Collection enters system, if all carried out to all input signals a series of
Computing, system cannot necessarily meet the demand of real-time；Meanwhile, garbage signal is incorporated in the middle of computing, by the accurate of system
Degree produces deleterious effect.Therefore it is necessary to detect to the vibration fragment in input signal, effectively vibration letter is only processed
Number, so as to mitigate pressure of the system on time complexity, the real-time of strengthening system improves the positioning precision of system.It is existing
Vibration signal intercept algorithm mainly using short-time energy method and short-time zero-crossing rate method, both algorithms all have algorithm and answer
Miscellaneous, operation efficiency is low, and real-time is not high and is vulnerable to the interference of noise, the defect such as reliability is not strong, it is difficult to meet MZ optical fiber week
Real-time monitoring and high reliability of boundary's intrusion detection instrument etc. are required.
Therefore, prior art existing defects, need to improve.
The content of the invention
The technical problem to be solved is for the deficiencies in the prior art, there is provided a kind of to be counted based on wave character
Effective vibration data intercept method.
Technical scheme is as follows：
A kind of effective vibration data intercept method counted based on wave character, wherein, comprise the following steps：
Step 101：Obtain replicate data DC [N] of original vibration data D [N]；
Step 102：Analysis DC [N], calculates the Wave crest and wave trough in DC [N]；
Step 103：Will be all in DC [N] neither the value of the crest again point of non-trough is set to 0；
Step 104：The value of point of all Wave crest and wave trough point absolute values in DC less than filtering threshold TW is set to into 0；
Step 105：It is M continuous fragment by DC point, the data length of each fragment is L=N/M；
Step 106：DC [1] is calculated to having an absolute value in DC [M] each fragment and DCS [1]～DCS [M]；
Step 107：Find maximum DCSMAX in DCS [1]～DCS [M]；
Step 108：Whether DCSMAX is judged more than peak threshold TM, if DCSMAX is not less than TM, into step
109, otherwise into step 113；
Step 109：Calculate the crest location PS of left side half and the crest location PE of right side half of DCS [1]～DCS [M]；
Step 110：Judge that PD=PE-PS, whether not less than data length threshold value TL, if PD is not less than TL, enters step
Rapid 111, otherwise into step 113；
Step 111：D [(PS-1) * L]～D [PE*L] is taken as effective vibration data；
Step 112：Submit effective vibration data to, go to step 114；
Step 113：Without effective vibration data, step 114 is gone to；
Step 114：Calculating terminates.
The described effective vibration data intercept method counted based on wave character, wherein in the step 101：It is original to shake
Dynamic data D [N] length N is 2 integral number power, if integral number power of the D data lengths less than 2, using the polishing of data 0.
The described effective vibration data intercept method counted based on wave character, wherein, in the step 102, calculate
During DC [N] medium wave peak trough, the continuous data more than 0 are just considered wave crest portions in one piece of data, and its maximum is this hop count
According to crest；The continuous data less than 0 are just considered wave trough portion in one piece of data, and its minimum of a value is the trough of this segment data.
The described effective vibration data intercept method counted based on wave character, wherein, in the step 105, by DC
It is divided into M continuous fragment, the data length of each fragment is L=N/M, and N, M are 2 integral number power.
The described effective vibration data intercept method counted based on wave character, wherein, the data of each fragment
Length is L=N/M, the value of M be less than N 2 integral number power, representative value is 32,64,128 or 256.
A kind of described effective vibration data intercept method counted based on wave character, wherein, in the step 106,
The value of i-th element DCS [i] (1≤i ＜ M) should be in DCS [1]～DCS [M]
A kind of described effective vibration data intercept method counted based on wave character, wherein, in the step 109,
The finding method of half crest location PS positions of left side is as follows：Centered on the PDCSMAX of DCSMAX positions, DCS is successively read
The value of [PDCSMAX-1]～DCS [1], first value less than DCSMAX/2 is then PSL, PS=PSL+1；In the step 109
The finding method of the crest location PE of right side half of DCS is as follows：Centered on the PDCSMAX of DCSMAX positions, DCS is successively read
The value of [PDCSMAX+1]～DCS [M], first value less than DCSMAX/2 is then PER, PE=PER+1.
Using such scheme, based on the wave character of vibration signal, carried out by the wave character to gathered data
Statistics, realizes filtering and intercepting to vibration signal to noise peak, effectively avoids the interference of garbage signal, reduces fortune
The rate of false alarm of number of times and system is calculated, algorithm operation quantity is little, real-time is high, reliability is high, false alarm rate is low.
Description of the drawings
Fig. 1 is the structural representation of the present invention.
Specific embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.
Embodiment 1
As shown in figure 1, the present invention provides a kind of effective vibration data intercept method counted based on wave character, wherein,
Comprise the following steps：
Step 101：Obtain replicate data DC [N] of original vibration data D [N]；
Step 102：Analysis DC [N], calculates the Wave crest and wave trough in DC [N]；
Step 103：Will be all in DC [N] neither the value of the crest again point of non-trough is set to 0；
Step 104：The value of point of all Wave crest and wave trough point absolute values in DC less than filtering threshold TW is set to into 0；
Step 105：It is M continuous fragment by DC point, the data length of each fragment is L=N/M；
Step 106：DC [1] is calculated to having an absolute value in DC [M] each fragment and DCS [1]～DCS [M]；
Step 107：Find maximum DCSMAX in DCS [1]～DCS [M]；
Step 108：Whether DCSMAX is judged more than peak threshold TM, if DCSMAX is not less than TM, into step
109, otherwise into step 113；
Step 109：Calculate the crest location PS of left side half and the crest location PE of right side half of DCS [1]～DCS [M]；
Step 110：Judge that PD=PE-PS, whether not less than data length threshold value TL, if PD is not less than TL, enters step
Rapid 111, otherwise into step 113；
Step 111：D [(PS-1) * L]～D [PE*L] is taken as effective vibration data；
Step 112：Submit effective vibration data to, go to step 114；
Step 113：Without effective vibration data, step 114 is gone to；
Step 114：Calculating terminates.
In the step 101：Original vibration data D [N] length N is 2 integral number power, if D data lengths are less than 2
Integral number power, using the polishing of data 0.
In the step 102, when calculating DC [N] medium wave peak trough, the continuous data more than 0 are just considered in one piece of data
Wave crest portions, its maximum is the crest of this segment data；The continuous data less than 0 are just considered wave trough portion in one piece of data,
Its minimum of a value is the trough of this segment data.
It is M continuous fragment by DC point in the step 105, the data length of each fragment is L=N/M, and N, M are 2
Integral number power.
The data length of each fragment is L=N/M, the value of M be less than N 2 integral number power, representative value is 32,
64th, 128 or 256.
In the step 106, the value of i-th element DCS [i] (1≤i ＜ M) should be in DCS [1]～DCS [M]
In the step 109, the finding method of the crest location PS positions of left side half is as follows：With DCSMAX positions
Centered on PDCSMAX, the value of DCS [PDCSMAX-1]～DCS [1] is successively read, first value less than DCSMAX/2 is then
PSL, PS=PSL+1；The finding method of the crest location PE of right side half of DCS is as follows in the step 109：It is in place with DCSMAX institutes
Centered on putting PDCSMAX, the value of DCS [PDCSMAX+1]～DCS [M] is successively read, first value less than DCSMAX/2 is then
PER, PE=PER+1.
Using technical scheme, based on the wave character of vibration signal, by the waveform to gathered data
Feature is counted, and realizes filtering and intercepting to vibration signal to noise peak, effectively avoids the interference of garbage signal,
The rate of false alarm of operation times and system is reduced, algorithm operation quantity is little, real-time is high, reliability is high, false alarm rate is low.
It should be appreciated that for those of ordinary skills, can according to the above description be improved or be converted,
And all these modifications and variations should all belong to the protection domain of claims of the present invention.
Claims (7)
1. it is a kind of based on wave character count effective vibration data intercept method, it is characterised in that comprise the following steps：
Step 101：Obtain replicate data DC [N] of original vibration data D [N]；
Step 102：Analysis DC [N], calculates the Wave crest and wave trough in DC [N]；
Step 103：Will be all in DC [N] neither the value of the crest again point of non-trough is set to 0；
Step 104：The value of point of all Wave crest and wave trough point absolute values in DC less than filtering threshold TW is set to into 0；
Step 105：It is M continuous fragment by DC point, the data length of each fragment is L=N/M；
Step 106：DC [1] is calculated to having an absolute value in DC [M] each fragment and DCS [1]～DCS [M]；
Step 107：Find maximum DCSMAX in DCS [1]～DCS [M]；
Step 108：Whether DCSMAX is judged more than peak threshold TM, it is no into step 109 if DCSMAX is not less than TM
Then enter step 113；
Step 109：Calculate the crest location PS of left side half and the crest location PE of right side half of DCS [1]～DCS [M]；
Step 110：Whether PD=PE-PS is judged not less than data length threshold value TL, if PD is not less than TL, into step
111, otherwise into step 113；
Step 111：D [(PS-1) * L]～D [PE*L] is taken as effective vibration data；
Step 112：Submit effective vibration data to, go to step 114；
Step 113：Without effective vibration data, step 114 is gone to；
Step 114：Calculating terminates.
2. it is as claimed in claim 1 to be based on effective vibration data intercept method that wave character is counted, it is characterised in that described
In step 101：Original vibration data D [N] length N is 2 integral number power, if integral number power of the D data lengths less than 2, makes
With the polishing of data 0.
3. it is as claimed in claim 1 to be based on effective vibration data intercept method that wave character is counted, it is characterised in that described
In step 102, when calculating DC [N] medium wave peak trough, the continuous data more than 0 are just considered wave crest portions in one piece of data, its
Maximum is the crest of this segment data；The continuous data less than 0 are just considered wave trough portion in one piece of data, and its minimum of a value is this
The trough of segment data.
4. it is as claimed in claim 1 to be based on effective vibration data intercept method that wave character is counted, it is characterised in that described
It is M continuous fragment by DC point in step 105, the data length of each fragment is L=N/M, and N, M are 2 integral number power.
5. it is as claimed in claim 1 to be based on effective vibration data intercept method that wave character is counted, it is characterised in that described
The data length of each fragment is L=N/M, the value of M be less than N 2 integral number power, representative value is 32,64,128 or 256.
6. a kind of effective vibration data intercept method counted based on wave character as claimed in claim 1, it is characterised in that
In the step 106, the value of i-th element DCS [i] (1≤i ＜ M) should be in DCS [1]～DCS [M]
7. a kind of effective vibration data intercept method counted based on wave character as claimed in claim 1, it is characterised in that
In the step 109, the finding method of the crest location PS positions of left side half is as follows：In being with DCSMAX positions PDCSMAX
The heart, is successively read the value of DCS [PDCSMAX-1]～DCS [1], and first value less than DCSMAX/2 is then PSL, PS=PSL+
1；The finding method of the crest location PE of right side half of DCS is as follows in the step 109：It is with DCSMAX positions PDCSMAX
Center, is successively read the value of DCS [PDCSMAX+1]～DCS [M], and first value less than DCSMAX/2 is then PER, PE=PER
+1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410568873.2A CN104346528B (en) | 2014-10-22 | 2014-10-22 | Effective vibration data interception method based on waveform characteristic statistics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410568873.2A CN104346528B (en) | 2014-10-22 | 2014-10-22 | Effective vibration data interception method based on waveform characteristic statistics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104346528A CN104346528A (en) | 2015-02-11 |
| CN104346528B true CN104346528B (en) | 2017-05-10 |
Family
ID=52502113
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410568873.2A Active CN104346528B (en) | 2014-10-22 | 2014-10-22 | Effective vibration data interception method based on waveform characteristic statistics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104346528B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105337259B (en) * | 2015-11-24 | 2018-09-07 | 珠海格力电器股份有限公司 | A kind of motor rotation blockage guard method and device |
| CN107024268B (en) * | 2017-05-19 | 2019-06-07 | 中国电子科技集团公司第四十一研究所 | A kind of effective vibration determination method based on threshold number excessively |
| CN107067608B (en) * | 2017-05-19 | 2019-03-05 | 中国电子科技集团公司第四十一研究所 | A kind of effective vibrational waveform intercept method based on three-level threshold determination |
| CN107886661B (en) * | 2017-11-10 | 2019-04-23 | 中国电子科技集团公司第四十一研究所 | A kind of intrusion event determination method based on vibrational energy accumulation |
| CN112562237B (en) * | 2020-12-11 | 2022-05-17 | 无锡科晟光子科技有限公司 | Optical fiber vibration system positioning accuracy optimization method applied to fence |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102435554A (en) * | 2011-09-08 | 2012-05-02 | 北京师范大学 | Method for acquiring farmland multiple-cropping index |
| CN103169456A (en) * | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
| JP5540959B2 (en) * | 2010-07-15 | 2014-07-02 | 横河電機株式会社 | Waveform measuring device |
-
2014
- 2014-10-22 CN CN201410568873.2A patent/CN104346528B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5540959B2 (en) * | 2010-07-15 | 2014-07-02 | 横河電機株式会社 | Waveform measuring device |
| CN102435554A (en) * | 2011-09-08 | 2012-05-02 | 北京师范大学 | Method for acquiring farmland multiple-cropping index |
| CN103169456A (en) * | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
Non-Patent Citations (2)
| Title |
|---|
| 基于AT89C51单片机的RS232串行数据截取器设计;郭付才，等.;《现代电子技术》;20120401;第35卷(第7期);第95-97页 * |
| 基于拉伸装置的光纤应变系数标定及误差分析;闫继送，等.;《光电工程》;20140131;第41卷(第1期);第23-28页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104346528A (en) | 2015-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104346528B (en) | Effective vibration data interception method based on waveform characteristic statistics | |
| KR102386155B1 (en) | How to protect your voice assistant from being controlled by machine learning-based silent commands | |
| CN104282103B (en) | A kind of vibrational waveform end-point detection new method based on short-time energy | |
| CN103234767A (en) | Nonlinear fault detection method based on semi-supervised manifold learning | |
| CN105067101A (en) | Fundamental tone frequency characteristic extraction method based on vibration signal for vibration source identification | |
| Cai et al. | A novel improved local binary pattern and its application to the fault diagnosis of diesel engine | |
| CN108090462A (en) | A kind of Emitter Fingerprint feature extracting method based on box counting dimension | |
| WO2019192252A1 (en) | Method and system for monitoring personal computer usage behaviour based on acoustic channel | |
| CN102624464A (en) | Method for decreasing influence of malicious users on cooperative spectrum sensing based on historical performance | |
| Wei et al. | A novel intrusion detection model for the CAN bus packet of in-vehicle network based on attention mechanism and autoencoder | |
| CN111598039A (en) | Abnormal sound monitoring method and monitoring device for petroleum pipeline monitoring | |
| Sun et al. | Non-contact diagnosis for gearbox based on the fusion of multi-sensor heterogeneous data | |
| CN103778916A (en) | Method and system for monitoring environmental sound | |
| CN103995969A (en) | Configurable optical fiber invasion event occurring end point detecting method and detector | |
| Petropulu et al. | Contactless monitoring of critical infrastructure [from the guest editors] | |
| CN103956013B (en) | A kind of wind and rain disturbing signal real-time judgment method | |
| CN110363033A (en) | A kind of chip security appraisal procedure and device | |
| CN107024268B (en) | A kind of effective vibration determination method based on threshold number excessively | |
| Zou et al. | OutletGuarder: detecting DarkSide ransomware by power factor correction signals in an electrical outlet | |
| CN107067608A (en) | A kind of effective vibrational waveform intercept method based on three-level threshold determination | |
| Li et al. | Research and simulation of network security situation prediction algorithm | |
| CN114186223A (en) | System and method for improving measurements of an intrusion detection system by transforming one-dimensional measurements into a multi-dimensional image | |
| CN108257364B (en) | Method for improving alarm reliability of distributed optical fiber monitoring system | |
| CN109945075A (en) | A kind of water supply line leakiness detection method and device | |
| CN113723200B (en) | Method for extracting time spectrum structural features of non-stationary signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190301 Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: China Electronics Technology Instrument and Meter Co., Ltd. Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: The 41st Institute of CETC |
|
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: CLP kesiyi Technology Co.,Ltd. Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd. |