CN107255513A - Skyscraper amplitude of oscillation measuring method and device - Google Patents
Skyscraper amplitude of oscillation measuring method and device Download PDFInfo
- Publication number
- CN107255513A CN107255513A CN201710340369.0A CN201710340369A CN107255513A CN 107255513 A CN107255513 A CN 107255513A CN 201710340369 A CN201710340369 A CN 201710340369A CN 107255513 A CN107255513 A CN 107255513A
- Authority
- CN
- China
- Prior art keywords
- mrow
- sequence
- amplitude
- imf components
- displacement signal
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Abstract
The embodiment of the present invention provides a kind of skyscraper amplitude of oscillation measuring method and device.Method includes：Gather amplitude of oscillation monitoring source displacement signal x (t)；Primary filtering is carried out to displacement signal x (t), to obtain displacement signal y (t) after filtering；Mode decomposition is carried out to displacement signal y (t), intrinsic mode function IMF components I is obtainedm(t), wherein, m=1,2 ..., M, M is IMF components Im(t) sum；IMF components I is calculated using presetting methodm(t) with displacement signal y (t) coefficient correlation, coefficient correlation sequence is obtained；Signal reconstruction is carried out according to coefficient correlation sequence, building amplitude of oscillation information is obtained.Skyscraper amplitude of oscillation measuring method provided in an embodiment of the present invention and device, by carrying out primary filtering and carrying out mode decomposition to filtered signal, and carry out signal reconstruction using the coefficient correlation of the IMF components after mode decomposition and filtered displacement signal, more accurate building amplitude of oscillation information is obtained, the precision of skyscraper amplitude of oscillation measurement is improved.
Description
Technical field
The present embodiments relate to structure detection technical field, more particularly to a kind of skyscraper amplitude of oscillation measuring method and dress
Put.
Background technology
With China's rapid development of economy, the skyscrapers of domestic each big and medium-sized cities setting up like the mushrooms after rain
Come, especially megalopolis such as Beijing, Shanghai, Guangzhou, Shenzhen etc..Skyscraper can be difficult to resist in construction period by various
The influence of factor and deform upon, add after the completion of construction throughout the year using the influences such as caused building materials aging, skyscraper deformation
Equally it is difficult to avoid that.Skyscraper can be produced under the bad weather such as typhoon and waved with the wind with beam wind, such as Guangzhou is built in high wind
Into skyscraper-Guangzhou international financial center, building is high 432 meters, and the maximum amplitude of oscillation in strong wind can reach 60 centimetres.If on
State deformation and the amplitude of oscillation and exceed safe range, gently then influence is used, the heavy then serious accident such as topple that collapses.Skyscraper
Different key positions the amplitude of oscillation monitoring have become skyscraper monitoring health status and safety monitoring important content.
At present, the common technology of skyscraper amplitude of oscillation monitoring includes displacement meter, and accelerometer, amesdial, linear differential becomes
The touch sensors such as depressor.But, touch sensor needs to be laid in the interested position of skyscraper, for building
High-rise is difficult to dispose touch sensor.In order to overcome the shortcoming of touch sensor, total powerstation, GPS, video measuring
Gradually it is applied in the monitoring of the skyscraper amplitude of oscillation Deng tether-free technologies.Total powerstation can monitor local engineering project deformation, Neng Goushi
Existing unattended round-the-clock automatic monitoring, has the disadvantage that measurement distance is limited.GPS absolute precisions can reach centimetre even grade,
But at the intensive place of skyscraper, gps satellite is blocked causes monitoring failure or influence monitoring accuracy with multipath effect.Video
Measurement is the set of the technologies such as Photogrammetry technology and computer vision, by way of under contactless state with image sequence
Moment records the 3 d space coordinate of moving object, and motion deformation of object etc. is analyzed.But it is due to image sensing
The limitation of device resolution ratio and frame frequency, it is difficult to the high-precision monitoring skyscraper amplitude of oscillation.The measurement of ground microwave interference has remote
Remote measurement, high-precision advantage, available for high level monitoring amplitude of oscillation monitoring, but are due to building oneself factor, environmental factor and set
, inevitably there is influence of noise, the precision of reduction skyscraper amplitude of oscillation monitoring in the influence of standby oneself factor.There is presently no
A kind of method can effectively improve the precision of skyscraper amplitude of oscillation measurement.
The content of the invention
In order to solve the above-mentioned technical problem, the embodiment of the present invention provides a kind of skyscraper amplitude of oscillation measuring method and device.
On the one hand, the embodiment of the present invention provides a kind of skyscraper amplitude of oscillation measuring method, including：
Gather amplitude of oscillation monitoring source displacement signal x (t)；
Primary filtering is carried out to institute displacement signal x (t), to obtain displacement signal y (t) after filtering；
Mode decomposition is carried out to institute displacement signal y (t), intrinsic mode function IMF components I is obtainedm(t), wherein, m=
1,2 ..., M, M are the IMF components Im(t) sum；
The IMF components I is calculated using presetting methodm(t) with institute displacement signal y (t) coefficient correlation, obtain related
Coefficient sequence；
Signal reconstruction is carried out according to the coefficient correlation sequence, building amplitude of oscillation information is obtained.
On the other hand, the embodiment of the present invention provides a kind of skyscraper amplitude of oscillation measurement apparatus, including：
Acquisition module, for gathering amplitude of oscillation monitoring source displacement signal x (t)；
Filtration module, for carrying out primary filtering to institute displacement signal x (t), to obtain displacement signal y (t) after filtering；
Mode decomposition module, for carrying out mode decomposition to institute displacement signal y (t), obtains IMF points of intrinsic mode function
Measure Im(t), wherein, m=1,2 ..., M, M is the IMF components Im(t) sum；
Computing module, for calculating the IMF components I using presetting methodm(t) it is related to institute displacement signal y (t)
Coefficient, obtains coefficient correlation sequence；
Reconstructed module, signal reconstruction is carried out according to the coefficient correlation sequence, obtains building amplitude of oscillation information.
Skyscraper amplitude of oscillation measuring method provided in an embodiment of the present invention and device, by carrying out primary filtering and to filter
Signal after ripple carries out mode decomposition, and utilizes the coefficient correlation of the IMF components after mode decomposition and filtered displacement signal
Signal reconstruction is carried out, more accurate building amplitude of oscillation information is obtained, improves the precision of skyscraper amplitude of oscillation measurement.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are this hairs
Some bright embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can be with root
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the schematic flow sheet of skyscraper amplitude of oscillation measuring method provided in an embodiment of the present invention；
Fig. 2 is skyscraper amplitude of oscillation measurement apparatus structural representation provided in an embodiment of the present invention；
Fig. 3 is another skyscraper amplitude of oscillation measurement apparatus structural representation provided in an embodiment of the present invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Fig. 1 is the schematic flow sheet of skyscraper amplitude of oscillation measuring method provided in an embodiment of the present invention, as shown in figure 1, institute
The method of stating includes：
Step 10, collection amplitude of oscillation monitoring source displacement signal x (t)；
Step 20, primary filtering is carried out to institute displacement signal x (t), to obtain displacement signal y (t) after filtering；
Step 30, to institute displacement signal y (t) carry out mode decomposition, obtain intrinsic mode function IMF components Im(t), its
In, m=1,2 ..., M, M is the IMF components Im(t) sum；
Step 40, presetting method is utilized to calculate the IMF components Im(t) with institute displacement signal y (t) coefficient correlation,
Obtain coefficient correlation sequence；
Step 50, according to the coefficient correlation sequence carry out signal reconstruction, obtain building amplitude of oscillation information.
When carrying out the skyscraper amplitude of oscillation using method provided in an embodiment of the present invention and measuring, the amplitude of oscillation to skyscraper first
Monitoring source displacement signal is acquired, and will each be gathered moment corresponding displacement and be recorded and come out with function representation, and be designated as x
(t)；Then primary filtering is carried out to displacement signal x (t), some obvious interference signals are filtered out, and filtered displacement is believed
Number it is designated as y (t)；After primary filtering is completed, using empirical mode decomposition (Empirical Mode Decomposition,
Abbreviation EMD) it is regular to filtered displacement signal y (t) progress mode decomposition, obtain M intrinsic mode function (intrinsic
Mode functions, abbreviation IMF) component Im(t), wherein, m=1,2 ..., M；Then according to obtained IMF components Im(t) and
Filtered displacement signal y (t) calculates IMF components I using default computational methodsm(t) with displacement signal y (t) phase
Relation number；Due to the frequency of M IMF component obtained when carrying out mode decomposition to y (t) be it is descending arrange,
That is I1(t) frequency highest, IM(t) frequency is maximum, and the actual amplitude of oscillation frequency of skyscraper is in certain frequency
Within the scope of, therefore according to IMF components Im(t) coefficient correlation with displacement signal y (t) can be divided to IMF components,
Swing signal is reconstructed, so as to obtain building amplitude of oscillation information.
Skyscraper amplitude of oscillation measuring method provided in an embodiment of the present invention, by carrying out primary filtering and to filtered
Signal carries out mode decomposition, and is believed using the coefficient correlation of IMF components and filtered displacement signal after mode decomposition
Number reconstruct, obtain more it is accurate building amplitude of oscillation information, improve the skyscraper amplitude of oscillation measure precision.
It is further, described that primary filtering is carried out to institute displacement signal x (t) on the basis of above-described embodiment, to obtain
Take filtering after displacement signal y (t) the step of be specially：
Institute displacement signal x (t) is decomposed using empirical mode decomposition EMD rules, K IMF components C is obtainedk(t)
With residual error R, wherein, k=1,2 ..., K；
According to default medium filtering window size respectively to the IMF components Ck(t) it is filtered, obtains after filtering
IMF components Dk(t)；
IMF components D after cumulative all filteringsk(t) filtered displacement signal y (t), is obtained.
Due to containing noise in the skyscraper amplitude of oscillation displacement signal that collects, had a certain impact for post-processing,
Therefore first filtering is first carried out before signal decomposition is carried out and removes larger interference signal.It is first sharp when being filtered for the first time
Displacement signal x (t) is subjected to empirical mode decomposition with EMD rules, K IMF components C is obtainedk(t) with residual error R, wherein, k=1,
2 ..., K, i.e.,：
Then according to the size of default medium filtering window to above-mentioned K IMF components Ck(t) it is filtered, is filtered
IMF components D after ripplek(t)
A kind of definition mode of optional medium filtering window size is：The medium filtering of first three IMF component decomposed
Window size is defined as 2+ (k+1) * 3, and the medium filtering window size of remaining IMF component decomposed is defined as 2+k2, wherein k is
The order of IMF components.
IMF components D after all median filter process that add upk(t) filtered displacement signal y (t), is obtained.
Then recycle filtered displacement signal y (t) to carry out mode decomposition, rebuild signal sequence, it is possible to
To more accurate amplitude of oscillation information.
Method provided in an embodiment of the present invention, is entered using the EMD methods with medium filtering to the displacement signal collected
The primary filtering of row, has filtered out larger noise, and it is more accurate reliable to make to rebuild obtained signal sequence, improves measurement essence
Degree.
On the basis of above-described embodiment, further, the utilization presetting method calculates the IMF components Im(t) with
The method of institute displacement signal y (t) coefficient correlation is specially：
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y
(t) average value,For Im(t) average value, n is number of samples.
Calculate the IMF components I that first noise reduction process is obtained by step 3m(t) with filtered displacement signal y (t)
Coefficient correlation, erased noise component, and carry out signal reconstruction, the skyscraper swing data for the influence of noise that is eliminated；
Specifically, mentioned before the frequency of M IMF component obtained when carrying out mode decomposition to y (t) be by greatly to
It is small to be arranged, that is to say, that I1(t) frequency highest, IM(t) frequency is maximum, and the actual amplitude of oscillation frequency of skyscraper
It is within certain frequency range, therefore according to IMF components Im(t) coefficient correlation with displacement signal y (t) can be right
IMF components are divided, and reconstruct swing signal, so as to obtain building amplitude of oscillation information.Understand, coefficient correlation is to reconstruct swing signal
Very big effect is played, the coefficient correlation can be calculated by following several formula：
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y
(t) average value,For Im(t) average value, n is number of samples.
M IMF components I can be obtained by calculatingm(t) with institute displacement signal y (t) coefficient R (y (t), Im
(t) division of signal sequence), is can be carried out according to obtained coefficient correlation, so that reconstruction signal sequence, obtains amplitude of oscillation letter
Breath.
The method provided using the present embodiment, can calculate IMF components Im(t) with institute displacement signal y (t) phase
Relation number, and using the coefficient of relationship calculated to IMF components Im(t) divided, reconstruction signal sequence is obtained after noise reduction
Amplitude of oscillation information, make measurement result more accurate.
It is further, described that signal reconstruction is carried out according to the coefficient correlation sequence on the basis of above-described embodiment, obtain
Take building the amplitude of oscillation information the step of be specially：
The IMF components I according to the coefficient correlation retrievalm(t) noise sequence and the boundary i of signal sequence with
And signal sequence and the boundary j of trend sequence,
Wherein, i is the IMF components corresponding to first local minimum coefficient correlation in the coefficient correlation sequence, and j is
The IMF components corresponding to coefficient correlation highest maximum in the coefficient correlation sequence after the corresponding coefficient correlations of i；
IMF components before j after i and j are reconstructed into signal sequence, the IMF components after j are reconstructed into trend sequence
Row；
Add up the signal sequence, obtains building amplitude of oscillation information；Add up the trend sequence, obtains building amplitude of oscillation trend letter
Breath.
The frequency of M IMF component obtained when carrying out mode decomposition to y (t) be it is descending arrange, it is and high
The actual amplitude of oscillation frequency of layer building is within certain frequency range, therefore according to IMF components Im(t) with displacement signal y
(t) coefficient correlation can be divided to IMF components, reconstruct swing signal, so as to obtain building amplitude of oscillation information.Specifically,
IMF components can be divided into noise sequence, signal sequence and the class of trend sequence three according to the characteristics of signal frequency, as its name suggests,
Noise sequence is the component of the interference signal caused by extraneous factors such as wind, earth shocks to signal；The building of signal sequence
The component of swing signal；Trend sequence is to build the component of amplitude of oscillation trend signal.By to IMF components Im(t) with displacement signal y
(t) analysis of coefficient correlation can obtain first local minimum R1, R1 correspondence IMF vector sequence of coefficient correlation sequence
In one-component i, we using i as noise sequence and signal sequence boundary；Same way, it is also possible to obtain coefficient correlation sequence
One-component j in highest maximum R2 in coefficient after middle R1, R2 correspondence IMF vector sequences, we regard j as signal
The boundary of sequence and trend sequence.Therefore, the IMF components before i and i are noise sequence, the IMF components after i before j and j
For signal sequence, the IMF components after j are trend sequence.
After noise sequence, signal sequence and trend sequence is divided, each IMF component in signal sequence is tired out
Plus can obtain removing the building swing signal after noise, each IMF component in trend sequence, which add up, to be obtained
The amplitude of oscillation general morphologictrend of building.
The method that the present embodiment is provided, utilizes IMF components Im(t) coefficient correlation with displacement signal y (t) is to IMF components
Divided, obtain removing the building swing signal after noise and the amplitude of oscillation general morphologictrend of building, not only increase high level
The precision of amplitude of oscillation measurement is built, the amplitude of oscillation general trend of building can also be obtained, subsequent analysis is convenient for.
On the basis of the various embodiments described above, further, the method for the collection amplitude of oscillation monitoring source displacement signal x (t)
Specially：Terrestrial microwave interferometer collection amplitude of oscillation monitoring source displacement signal x (t) is utilized in skyscraper lee side.
The method that the present embodiment is provided carries out amplitude of oscillation measurement using terrestrial microwave interferometer, the ground when carrying out data acquisition
Microwave interferometer needs to be placed on tripod, and wind-force can cause tripod to tremble, and terrestrial microwave interferometer measurement precision is very high,
Trembling for tripod can cause larger interference signal, so as to influence the precision of data acquisition.Gather high in skyscraper lee side
Layer building amplitude of oscillation monitoring source displacement signal can effectively reduce influence of the wind-force to ground microwave interferometer.
Fig. 2 is skyscraper amplitude of oscillation measurement apparatus structural representation provided in an embodiment of the present invention, as shown in Fig. 2 device
Including：Acquisition module 1, filtration module 2, mode decomposition module 3, computing module 4 and reconstructed module 5, wherein, acquisition module 1 is used
In collection amplitude of oscillation monitoring source displacement signal x (t)；Filtration module 2 is used to carry out institute displacement signal x (t) primary filtering, to obtain
Take displacement signal y (t) after filtering；Mode decomposition module 3 is used to carry out mode decomposition to institute displacement signal y (t), obtains intrinsic
Mode function IMF components Im(t), wherein, m=1,2 ..., M, M is the IMF components Im(t) sum；Computing module 4 is used for
The IMF components I is calculated using presetting methodm(t) with institute displacement signal y (t) coefficient correlation, phase relation number sequence is obtained
Row；Reconstructed module 5 carries out signal reconstruction according to the coefficient correlation sequence, obtains building amplitude of oscillation information.
It is 1 pair high by acquisition module first when carrying out the measurement of the skyscraper amplitude of oscillation using device provided in an embodiment of the present invention
The amplitude of oscillation monitoring source displacement signal of layer building is acquired, and the corresponding displacement of each collection moment is recorded and function table is used
Show to come, be designated as x (t)；Then filtration module 2 carries out primary filtering to displacement signal x (t), filters out some obvious interference letters
Number, and filtered displacement signal is designated as y (t)；After primary filtering is completed, EMD rules are utilized by mode decomposition module 3
Mode decomposition is carried out to filtered displacement signal y (t), M IMF components I is obtainedm(t), wherein, m=1,2 ..., M；Then
Computing module 4 is according to obtained IMF components Im(t) calculated with filtered displacement signal y (t) using default computational methods
Go out IMF components Im(t) with displacement signal y (t) coefficient correlation；Due to the M IMF obtained when carrying out mode decomposition to y (t)
The frequency of component is descending arranged, that is to say, that I1(t) frequency highest, IM(t) frequency is maximum, and high-rise
The actual amplitude of oscillation frequency of building is that within certain frequency range, therefore reconstructed module 5 is according to IMF components ImAnd position (t)
Shifting signal y (t) coefficient correlation can be divided to IMF components, reconstruct swing signal, so as to obtain building amplitude of oscillation information.
Skyscraper amplitude of oscillation measurement apparatus provided in an embodiment of the present invention, by carrying out primary filtering and to filtered
Signal carries out mode decomposition, and is believed using the coefficient correlation of IMF components and filtered displacement signal after mode decomposition
Number reconstruct, obtain more it is accurate building amplitude of oscillation information, improve the skyscraper amplitude of oscillation measure precision.
Fig. 3 is another skyscraper amplitude of oscillation measurement apparatus structural representation provided in an embodiment of the present invention, as shown in figure 3,
Device includes：Acquisition module 1, filtration module 2, mode decomposition module 3, computing module 4 and reconstructed module 5, wherein filtration module 2
Including：Resolving cell 21, filter unit 22 and summing elements 23；Wherein, the effect of acquisition module 1 and phase described in above-described embodiment
Together, here is omitted.Resolving cell 21 in filtration module 2 is used for using empirical mode decomposition EMD rules to the displacement
Signal x (t) is decomposed, and obtains K IMF components Ck(t) with residual error R, wherein, k=1,2 ..., K；Filter unit 22 is used for root
According to default medium filtering window size respectively to the IMF components Ck(t) it is filtered, obtains the IMF components D after filteringk
(t)；Summing elements 23 are for the IMF components D after all filterings that add upk(t) filtered displacement signal y (t), is obtained.Mode
The function of decomposing module 3, computing module 4 and reconstructed module 5 is identical with the embodiment, and here is omitted.
Contain noise in the skyscraper amplitude of oscillation displacement signal collected due to acquisition module 1, have one for post-processing
Fixed influence, therefore the first filtering of first progress removes larger interference signal before mode decomposition module 3 carries out signal decomposition.
Filtration module 2 when being filtered for the first time, and displacement signal x (t) is carried out empirical modal point by resolving cell 21 first with EMD rules
Solution, obtains K IMF components Ck(t) with residual error R, wherein, k=1,2 ..., K, i.e.,：
Then filter unit 22 according to the size of default medium filtering window to above-mentioned K IMF components Ck(t) filtered
Ripple, obtains filtered IMF components Dk(t)
A kind of definition mode of optional medium filtering window size is：The medium filtering of first three IMF component decomposed
Window size is defined as 2+ (k+1) * 3, and the medium filtering window size of remaining IMF component decomposed is defined as 2+k2, wherein k is
The order of IMF components.
Summing elements 23 add up the IMF components D after all median filter processk(t) filtered displacement signal y, is obtained
(t)。
Then recycle filtered displacement signal y (t) to carry out mode decomposition, rebuild signal sequence, it is possible to
To more accurate amplitude of oscillation information.
Device provided in an embodiment of the present invention, utilizes resolving cell 21,23 pairs of collection moulds of filter unit 22 and summing elements
Displacement signal that block 1 is collected carries out primary filtering, has filtered out larger noise, makes to rebuild obtained signal sequence more smart
It is really reliable, improve measurement accuracy.
On the basis of above-described embodiment, further, the method that the computing module obtains the coefficient correlation sequence
Specially：
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y
(t) average value,For Im(t) average value, n is number of samples.
Specifically, mentioning the M IMF component that mode decomposition module is obtained when carrying out mode decomposition to y (t) before
Frequency is descending arranged, that is to say, that I1(t) frequency highest, IM(t) frequency is maximum, and skyscraper
Actual amplitude of oscillation frequency is within certain frequency range, therefore according to IMF components Im(t) it is related to displacement signal y (t)
Coefficient can be divided to IMF components, reconstruct swing signal, so as to obtain building amplitude of oscillation information.Understand, coefficient correlation counterweight
Structure swing signal plays very big effect, and the coefficient correlation can be calculated by following several formula：
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y
(t) average value,For Im(t) average value, n is number of samples.
Computing module can obtain M IMF components I by calculatingm(t) with institute displacement signal y (t) coefficient R
(y(t),Im(t) division of signal sequence), is can be carried out according to obtained coefficient correlation, so that reconstruction signal sequence, is obtained
Amplitude of oscillation information.
The device provided using the present embodiment, computing module can calculate IMF components Im(t) with institute's displacement signal
Y (t) coefficient correlation, and the coefficient of relationship calculated is utilized to IMF components I by reconstructed modulem(t) divided, reconstructed
Signal sequence, obtains the amplitude of oscillation information after noise reduction, makes measurement result more accurate.
On the basis of above-described embodiment, further, the method that the reconstructed module obtains building amplitude of oscillation information is specific
For：
The IMF components I according to the coefficient correlation retrievalm(t) noise sequence and the boundary i of signal sequence with
And signal sequence and the boundary j of trend sequence,
Wherein, i is the IMF components corresponding to first local minimum coefficient correlation in the coefficient correlation sequence, and j is
The IMF components corresponding to coefficient correlation highest maximum in the coefficient correlation sequence after the corresponding coefficient correlations of i；
IMF components before j after i and j are reconstructed into signal sequence, the IMF components after j are reconstructed into trend sequence
Row；
Add up the signal sequence, obtains building amplitude of oscillation information；Add up the trend sequence, obtains building amplitude of oscillation trend letter
Breath.
When mode decomposition module carries out mode decomposition to y (t), the frequency of M obtained IMF component be it is descending enter
Row arrangement, and the actual amplitude of oscillation frequency of skyscraper is that within certain frequency range, therefore reconstructed module is according to meter
Calculate the IMF components I that module is calculatedm(t) coefficient correlation with displacement signal y (t) can be divided to IMF components, weight
Structure swing signal, so as to obtain building amplitude of oscillation information.
Specifically, IMF components can be divided into noise sequence, signal sequence and trend sequence according to the characteristics of signal frequency
Three classes, as its name suggests, noise sequence are the components of the interference signal caused by extraneous factors such as wind, earth shocks to signal；Letter
The component of the swing signal of the building of number sequence；Trend sequence is to build the component of amplitude of oscillation trend signal.By to IMF components Im
(t) analysis with displacement signal y (t) coefficient correlation can obtain first local minimum R1, R1 of coefficient correlation sequence
One-component i in correspondence IMF vector sequences, we using i as noise sequence and signal sequence boundary；Same way, it is also possible to
Obtain the one-component in highest maximum R2 in the coefficient in coefficient correlation sequence after R1, R2 correspondence IMF vector sequences
J, we regard j as signal sequence and the boundary of trend sequence.Therefore, the IMF components before i and i are j after noise sequence, i
And the IMF components before j are signal sequence, the IMF components after j are trend sequence.
After noise sequence, signal sequence and trend sequence is divided, each IMF component in signal sequence is tired out
Plus can obtain removing the building swing signal after noise, each IMF component in trend sequence, which add up, to be obtained
The amplitude of oscillation general morphologictrend of building.
The device that the present embodiment is provided, the IMF components I that reconstructed module is calculated using computing modulem(t) believe with displacement
Number y (t) coefficient correlation is divided to IMF components, and the amplitude of oscillation for obtaining removing building swing signal after noise and building is total
Body variation tendency, not only increases the precision of skyscraper amplitude of oscillation measurement, can also obtain the amplitude of oscillation general trend of building, be easy to
Carry out subsequent analysis.
On the basis of the various embodiments described above, the method tool of the acquisition module collection amplitude of oscillation monitoring source displacement signal x (t)
Body is：Terrestrial microwave interferometer collection amplitude of oscillation monitoring source displacement signal x (t) is utilized in skyscraper lee side.
In the device that the present embodiment is provided, acquisition module carries out amplitude of oscillation measurement using terrestrial microwave interferometer, is entering line number
Need to be placed on tripod according to terrestrial microwave interferometer during collection, wind-force can cause tripod to tremble, terrestrial microwave interferometer
Measurement accuracy is very high, and trembling for tripod can cause larger interference signal, so as to influence the precision of data acquisition.In skyscraper
Lee side's collection skyscraper amplitude of oscillation monitoring source displacement signal can effectively reduce influence of the wind-force to ground microwave interferometer.
Device embodiment described above is only schematical, wherein the unit illustrated as separating component can
To be or may not be physically separate, the part shown as unit can be or may not be physics list
Member, you can with positioned at a place, or can also be distributed on multiple NEs.It can be selected according to the actual needs
In some or all of module realize the purpose of this embodiment scheme.Those of ordinary skill in the art are not paying creativeness
Work in the case of, you can to understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
Realized by the mode of software plus required general hardware platform, naturally it is also possible to related work(is realized by hardware processor
Can module.Understood based on such, the part that above-mentioned technical proposal substantially contributes to prior art in other words can be with
The form of software product is embodied, and the computer software product can be stored in a computer-readable storage medium, such as ROM/
RAM, magnetic disc, CD etc., including some instructions to cause a computer equipment (can be personal computer, server, or
Person's network equipment etc.) perform method described in some parts of each embodiment or embodiment.
Finally it should be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
The present invention is described in detail with reference to the foregoing embodiments, it will be understood by those within the art that：It still may be used
To be modified to the technical scheme described in foregoing embodiments, or equivalent substitution is carried out to which part technical characteristic；
And these modification or replace, do not make appropriate technical solution essence depart from various embodiments of the present invention technical scheme spirit and
Scope.
Claims (10)
1. a kind of skyscraper amplitude of oscillation measuring method, it is characterised in that including：
Gather amplitude of oscillation monitoring source displacement signal x (t)；
Primary filtering is carried out to institute displacement signal x (t), to obtain displacement signal y (t) after filtering；
Mode decomposition is carried out to institute displacement signal y (t), intrinsic mode function IMF components I is obtainedm(t), wherein, m=1,
2 ..., M, M are the IMF components Im(t) sum；
The IMF components I is calculated using presetting methodm(t) with institute displacement signal y (t) coefficient correlation, coefficient correlation is obtained
Sequence；
Signal reconstruction is carried out according to the coefficient correlation sequence, building amplitude of oscillation information is obtained.
2. according to the method described in claim 1, it is characterised in that described that primary filtering is carried out to institute displacement signal x (t),
It is specially to obtain the step of displacement signal y (t) after filtering：
Institute displacement signal x (t) is decomposed using empirical mode decomposition EMD rules, K IMF components C is obtainedk(t) it is and residual
Poor R, wherein, k=1,2 ..., K；
According to default medium filtering window size respectively to the IMF components Ck(t) it is filtered, obtains IMF points after filtering
Measure Dk(t)；
IMF components D after cumulative all filteringsk(t) filtered displacement signal y (t), is obtained.
3. according to the method described in claim 1, it is characterised in that the utilization presetting method calculates the IMF components Im(t)
Method with institute displacement signal y (t) coefficient correlation is specially：
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
<mo>-</mo>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<msup>
<mo>&rsqb;</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y's (t)
Average value,For Im(t) average value, n is number of samples.
4. according to the method described in claim 1, it is characterised in that described that signal weight is carried out according to the coefficient correlation sequence
Structure, obtaining the step of building amplitude of oscillation information is specially：
The IMF components I according to the coefficient correlation retrievalm(t) noise sequence and the boundary i and letter of signal sequence
The boundary j of number sequence and trend sequence,
Wherein, i is the IMF components corresponding to first local minimum coefficient correlation in the coefficient correlation sequence, and j is described
The IMF components corresponding to coefficient correlation highest maximum in coefficient correlation sequence after the corresponding coefficient correlations of i；
IMF components before j after i and j are reconstructed into signal sequence, the IMF components after j are reconstructed into trend sequence；
Add up the signal sequence, obtains building amplitude of oscillation information；Add up the trend sequence, obtains building amplitude of oscillation tendency information.
5. according to any described method of Claims 1-4, it is characterised in that the collection amplitude of oscillation monitors source displacement signal x
(t) method is specially：Terrestrial microwave interferometer collection amplitude of oscillation monitoring source displacement signal x (t) is utilized in skyscraper lee side.
6. a kind of skyscraper amplitude of oscillation measurement apparatus, it is characterised in that including：
Acquisition module, for gathering amplitude of oscillation monitoring source displacement signal x (t)；
Filtration module, for carrying out primary filtering to institute displacement signal x (t), to obtain displacement signal y (t) after filtering；
Mode decomposition module, for carrying out mode decomposition to institute displacement signal y (t), obtains intrinsic mode function IMF components Im
(t), wherein, m=1,2 ..., M, M is the IMF components Im(t) sum；
Computing module, for calculating the IMF components I using presetting methodm(t) with institute displacement signal y (t) coefficient correlation,
Obtain coefficient correlation sequence；
Reconstructed module, signal reconstruction is carried out according to the coefficient correlation sequence, obtains building amplitude of oscillation information.
7. device according to claim 6, it is characterised in that the filtration module includes：
Resolving cell, for being decomposed using empirical mode decomposition EMD rules to institute displacement signal x (t), obtains K IMF
Component Ck(t) with residual error R, wherein, k=1,2 ..., K；
Filter unit, for according to default medium filtering window size respectively to the IMF components Ck(t) it is filtered, obtains
IMF components D after filteringk(t)；
Summing elements, for the IMF components D after all filterings that add upk(t) filtered displacement signal y (t), is obtained.
8. device according to claim 6, it is characterised in that the computing module obtains the side of the coefficient correlation sequence
Method is specially：
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
<mo>-</mo>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<msup>
<mo>&rsqb;</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, R (y (t), Im(t)) it is coefficient correlation, y (t) is displacement signal, Im(t) it is m rank IMF components,For y's (t)
Average value,For Im(t) average value, n is number of samples.
9. device according to claim 6, it is characterised in that the reconstructed module obtains the method tool of building amplitude of oscillation information
Body is：
The IMF components I according to the coefficient correlation retrievalm(t) noise sequence and the boundary i and letter of signal sequence
The boundary j of number sequence and trend sequence,
Wherein, i is the IMF components corresponding to first local minimum coefficient correlation in the coefficient correlation sequence, and j is described
The IMF components corresponding to coefficient correlation highest maximum in coefficient correlation sequence after the corresponding coefficient correlations of i；
IMF components before j after i and j are reconstructed into signal sequence, the IMF components after j are reconstructed into trend sequence；
Add up the signal sequence, obtains building amplitude of oscillation information；Add up the trend sequence, obtains building amplitude of oscillation tendency information.
10. according to any described device of claim 6 to 9, it is characterised in that acquisition module collection amplitude of oscillation monitoring source position
Shifting signal x (t) method is specially：Terrestrial microwave interferometer collection amplitude of oscillation monitoring source displacement letter is utilized in skyscraper lee side
Number x (t).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710340369.0A CN107255513B (en) | 2017-05-15 | 2017-05-15 | Skyscraper amplitude of oscillation measurement method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710340369.0A CN107255513B (en) | 2017-05-15 | 2017-05-15 | Skyscraper amplitude of oscillation measurement method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107255513A true CN107255513A (en) | 2017-10-17 |
| CN107255513B CN107255513B (en) | 2019-10-18 |
Family
ID=60027589
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710340369.0A Active CN107255513B (en) | 2017-05-15 | 2017-05-15 | Skyscraper amplitude of oscillation measurement method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107255513B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108647635A (en) * | 2018-05-09 | 2018-10-12 | 黑龙江大学 | The single channel image Blind extracting method for receiving model is supplied based on dividing equally to recycle |
| CN109031422A (en) * | 2018-08-09 | 2018-12-18 | 吉林大学 | A kind of seismic signal noise suppressing method based on CEEMDAN and Savitzky-Golay filtering |
| CN110456392A (en) * | 2019-08-23 | 2019-11-15 | 北京建筑大学 | A kind of tower crane beam position precise positioning reliability verification method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105547465A (en) * | 2015-12-08 | 2016-05-04 | 华北电力大学（保定） | Transformer vibration signal winding state feature extraction method |
-
2017
- 2017-05-15 CN CN201710340369.0A patent/CN107255513B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105547465A (en) * | 2015-12-08 | 2016-05-04 | 华北电力大学（保定） | Transformer vibration signal winding state feature extraction method |
Non-Patent Citations (4)
| Title |
|---|
| 刁建鹏: ""微波干涉测量在高层建筑动态监测中的应用"", 《振动与冲击》 * |
| 刘代志: "《地球物理与海洋安全》", 31 October 2009, 西安地图出版社 * |
| 文振华: "《航空发动机静电监测技术》", 30 September 2014, 知识产权出版社 * |
| 王欣 等: ""经验模态分解在高层建筑GPS变形监测分析中的应用"", 《城市勘测》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108647635A (en) * | 2018-05-09 | 2018-10-12 | 黑龙江大学 | The single channel image Blind extracting method for receiving model is supplied based on dividing equally to recycle |
| CN109031422A (en) * | 2018-08-09 | 2018-12-18 | 吉林大学 | A kind of seismic signal noise suppressing method based on CEEMDAN and Savitzky-Golay filtering |
| CN110456392A (en) * | 2019-08-23 | 2019-11-15 | 北京建筑大学 | A kind of tower crane beam position precise positioning reliability verification method |
| CN110456392B (en) * | 2019-08-23 | 2021-05-11 | 北京建筑大学 | Method for verifying accurate positioning reliability of position of cross arm of building tower crane |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107255513B (en) | 2019-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shen et al. | Remote sensing and evaluation of the wetland ecological degradation process of the Zoige Plateau Wetland in China | |
| CN103983273B (en) | A kind of real-time step-size estimation method based on acceleration transducer | |
| CN102937721B (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
| CN107255513A (en) | Skyscraper amplitude of oscillation measuring method and device | |
| CN110427857A (en) | A kind of transmission line of electricity geological disasters analysis method based on Remote Sensing Data Fusion Algorithm | |
| CN104459784B (en) | Two-dimensional Lg wave Q value tomographic imaging method based on single station data, double station data and double event data | |
| CN102023010B (en) | MEMS (micro-electromechanical system)-based wavelet field multisensor information fusion system and fusion method | |
| CN105426814A (en) | Old people stumbling detection method based on handset | |
| CN106990128A (en) | Plaque image and pollution characteristic coupling method in regional tailing pond environmental risk assessment | |
| CN105758404A (en) | Real-time positioning method and system of intelligent equipment | |
| CN106482917A (en) | A kind of detection method of cable-stayed bridge main-beam dynamic deflection | |
| CN110184885A (en) | A method of surface evenness is tested based on smart phone | |
| CN107862687A (en) | A kind of early warning system for being used to monitor agricultural pest | |
| CN102540252A (en) | High-precision median stacking method on basis of cross-correlation | |
| CN103473559A (en) | SAR image change detection method based on NSCT domain synthetic kernels | |
| CN109916485B (en) | Dynamic vehicle weighing method and device | |
| CN101487699B (en) | High-precision APS solar sensor and its implementing method | |
| CN106157616A (en) | A kind of magnitude of traffic flow short-term prediction device | |
| CN104036283A (en) | System and method for determining the activity of a mobile element | |
| CN112882032B (en) | Method and device for dynamically monitoring geological disaster SAR in key area of gas pipeline | |
| CN103064128B (en) | Based on the gravity field recover method of interstellar distance error model | |
| CN114812491B (en) | Transmission line earth surface deformation early warning method and device based on long-time sequence analysis | |
| Li et al. | Dynamic deformation monitoring of an offshore platform structure with accelerometers | |
| CN103376465A (en) | System and method for earthquake real-time analysis of building floor and storage medium | |
| CN105572077B (en) | A kind of gas Rayleigh Brillouin spectrum signal restoring method |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |