CN104990717B - A kind of magnetic-levitation train sensor signal processing method - Google Patents

A kind of magnetic-levitation train sensor signal processing method Download PDF

Info

Publication number
CN104990717B
CN104990717B CN201510444147.4A CN201510444147A CN104990717B CN 104990717 B CN104990717 B CN 104990717B CN 201510444147 A CN201510444147 A CN 201510444147A CN 104990717 B CN104990717 B CN 104990717B
Authority
CN
China
Prior art keywords
sampled value
flexible beam
sensor
virtual
value
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
Application number
CN201510444147.4A
Other languages
Chinese (zh)
Other versions
CN104990717A (en
Inventor
周丹峰
余佩倡
崔鹏
李�杰
王连春
李金辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510444147.4A priority Critical patent/CN104990717B/en
Publication of CN104990717A publication Critical patent/CN104990717A/en
Application granted granted Critical
Publication of CN104990717B publication Critical patent/CN104990717B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

A kind of magnetic-levitation train sensor signal processing method, the sampled value of this method continuous some sampling instants including acquisition sensor current time sampled value and before current time, is used as constraints to build virtual flexible beam sampled value, forms optimum fit curve;The sampled value of next sampling instant is predicted according to the extended position of virtual flexible beam;Next sampling instant sampled value is obtained, if the sampled value of next sampling instant and the difference of predicted value exceed the threshold value of setting, then it is assumed that spurious signal occur in sensor, it is necessary to be rejected to it, and sampled value is replaced with predicted value, it is on the contrary then retain sampled value.Effectively the burr that mixes in sensor signal, outlier can be rejected using this method, i.e., the spurious signal to sensor is corrected;The processing method provided using the present invention both can have been judged and corrected to spurious signal, the form fit sensing data of virtual flexible beam can also be utilized, so as to reach the effect of smooth noise reduction.

Description

A kind of magnetic-levitation train sensor signal processing method
Technical field
It is more particularly to a kind of suitable for medium-and low-speed maglev train the invention mainly relates to medium-and low-speed maglev train field The processing method of sensor signal.
Background technology
Medium-and low-speed maglev train is the Vehicular system that a kind of dependence electromagnetic force realizes contactless suspension, and its core system is outstanding Floating control system.Suspension control system is used for various sensors, and such as levitation gap sensor, electromagnet vertical acceleration are sensed The signals such as device, current sensor are acquired, and suitable voltage is exported after certain computing to control the electric current of electromagnet. However, due to working environment, the signal that these sensors are obtained there may come a time when there is deceptive information.For example, floating in magnetic When the traction electric machine of train and the work of other high-power electrical equipments, the interference of electromagnetic field stronger in short-term may be produced, may Gap sensor, current sensor etc. can be caused to export of short duration pulse signal, this wrong pulse can be produced to suspension system Raw impact even results in levitating electromagnet and collided with track;For another example, when the gap sensor of magnetic-levitation train is by track seam, survey Amount gap can be because track be discontinuous and exports larger value, and this spurious signal can equally disturb the normal value of clearance measurement.
Therefore, how to reduce influence of this spurious signal to suspension control system, to ensure that levitation gap is steady, turn into The problem of those skilled in the art's urgent need to resolve.
The content of the invention
The technical problem to be solved in the present invention, which is to provide one kind, to be believed in real time sensor in magnetic-levitation train running Number detected and judged, the spurious signal detected is rejected, so as to reduce sensor spurious signal to suspension control The influence of system, it is ensured that train levitation gap is steady.
The present invention solve its technical problem use technical scheme be, a kind of magnetic-levitation train sensor signal processing method, Comprise the following steps:
Step 1:Obtain sensor current time sampled value and before current time continuous some sampling instants sampling Value, virtual flexible beam is built using sampled value as constraints, forms optimum fit curve;
Wherein, the sampled value is as the virtual flexible beam of constraints structure:
Sampled value is connected on the longitudinal axis between virtual flexible beam by virtual linear spring, and virtual linear spring is according to adopting Sample value and the position difference of virtual flexible beam produce pulling force effect on virtual flexible beam so that virtual flexible beam be subjected to displacement and Bending, multiple sampled values are acted on simultaneously, the stress and moment of flexure of beam is reached balance, are determined that deflection of beam is deformed, are built in flat Virtual flexible equalizer bar under weighing apparatus state;
Step 2:The sampled value of next sampling instant is predicted according to the extended position of virtual flexible beam;
Step 3:Next sampling instant sampled value is obtained, if the sampled value of next sampling instant and the difference of predicted value exceed The threshold value of setting, then it is assumed that occur spurious signal in sensor, it is on the contrary then think that sensor signal is normal.
It is preferred that, when spurious signal occurs in the sensor, sampled value, on the contrary then reservation sampling are replaced with predicted value Value.
It is preferred that, the bending rigidity EI and virtual linear rigidity k s of virtual flexible beam can adjust, to be adapted to different frequencies Sensor signal with scope.
It is preferred that, the sampling instant number can be determined according to the type and purposes of sensor.
It is preferred that, time interval is identical in once virtual flexible beam building process between the neighbouring sample moment.
It is preferred that, the time interval between the sampling instant can be determined according to the type and purposes of sensor.
It is preferred that, in the sampling interval of the next sampling instant of prediction, it can be more than or equal between the sampling of sampling instant Every.
The magnetic-levitation train sensor signal processing method of the present invention can effectively to mix in sensor signal burr, Outlier is rejected, i.e., the spurious signal to sensor is corrected;Sampled value is determined and virtual soft using virtual linear spring Restraining force between property beam, makes these sampled values and virtual flexible beam reach the balance on mechanics, therefore, this balance can have The sampled data of effect fitting sensor, obtains sampled data variation tendency more smooth credible;The bending rigidity of virtual flexible beam Can arbitrarily it be adjusted with virtual linear spring rate, to be adapted to the sensor signal of different frequency bands scope.Using present invention processing Method both can be judged and be corrected to spurious signal, can also utilize the form fit sensing data of virtual flexible beam, So as to reach the effect of smooth noise reduction.
Brief description of the drawings
Fig. 1 is a kind of flow chart of embodiment of magnetic-levitation train sensor signal processing method of the present invention;
Fig. 2 is restriction relation coordinate schematic diagram between virtual flexible beam and sampled value;
Fig. 3 is the coordinate schematic diagram that spurious signal is judged and corrected with virtual flexible beam;
Fig. 4 is schematic diagram of the magnetic-levitation train gap sensor by track seam;
Fig. 5 is that gap sensor passes through spurious signal correction principle schematic diagram during track seam.
Embodiment
In order that those skilled in the art more fully understands technical scheme, below in conjunction with the accompanying drawings to the present invention It is described in further detail.
Reference picture 1, magnetic-levitation train sensor signal processing method of the present invention, comprises the following steps:
Step 1:Obtain sensor current time sampled value and before current time continuous some sampling instants sampling Value, virtual flexible beam is built using sampled value as constraints, forms optimum fit curve;
Obtain sensor current time sampled value and before current time some sampling instants sampled value, by sampled value Virtual flexible beam is built as constraints, the position of virtual flexible beam is fixed by multiple sampled values, virtual flexible beam Shape is the Smoothing fit of these sampled values, so as to form optimum fit curve.
It is preferred that, the sampling period at neighbouring sample moment is identical in once virtual flexible beam building process.
Step 2:The sampled value of next sampling instant is predicted according to the extended position of virtual flexible beam;
Virtual flexible beam change in shape trend represents the real change trend of signal, and overhanging section of virtual flexible beam can be with As the Approximate prediction of future signal, thus the sampling of next sampling instant can be predicted according to the extended position of virtual flexible beam Value.
Step 3:Next sampling instant sampled value is obtained, if the sampled value of next sampling instant and the difference of predicted value exceed The threshold value of setting, then it is assumed that occur spurious signal in sensor, it is on the contrary then think that sensor signal is normal.
Due to the controlled device as magnetic-levitation train, the electromagnet and bogie inertia of magnetic-levitation train are larger, therefore suspend The true measurement signal of system can not possibly have larger mutation.If next sampling instant sensor signal value and virtual flexible beam The difference of prediction signal exceedes the threshold value of setting, then it is considered that the sensor signal collected is false.
It is preferred that, in the sampling interval of the next sampling instant of prediction, it can be more than or equal between the sampling of sampling instant Every.
By above-mentioned embodiment, medium-and low-speed maglev train is obtained after sensor signal without using wave filter, can be to void Glitch is judged and corrected, and the form fit sensing data of virtual flexible beam can also be utilized, effectively to sensor The burr that mixes in signal, outlier are rejected, so as to reach the effect of smooth noise reduction.
Alternatively, the sampled value is built as constraints in the step 1 of above-mentioned embodiment empty The step of intending flexible beam can further comprise:The sampled value passes through virtual linear spring on the longitudinal axis between virtual flexible beam Connection, virtual linear spring produces pulling force effect on virtual flexible beam according to sampled value and the position difference of virtual flexible beam, So that virtual flexible beam is subjected to displacement and bent.Multiple sampled values are acted on simultaneously, may finally reach the stress and moment of flexure of beam Balance, deflection of beam deformation also can determine that, so as to build the virtual flexible equalizer bar under poised state.Using virtual Hookean spring determines the restraining force between sampled data value and virtual flexible beam, reaches these sampled values and virtual flexible beam Balance on mechanics, therefore, this balance can effectively be fitted the sampled data of sensor, obtain sampled data variation tendency more Plus it is smooth credible.
In actual application, it is preferred that threshold value could be arranged to three times of average forecasting error, i.e., when next sampling The error of the sampled value at quarter and virtual flexible beam predicted value is more than the three times of average forecasting error, then to illustrate sensor collection Abnormal signal, now just answers discarded sampling signal.The average forecasting error can be by all samplings on virtual flexible beam Error Absolute Value between the sampled value at moment and virtual flexible beam is averagely obtained.
Further, the step 3 of above-mentioned embodiment can also include, when spurious signal occurs in the sensor, use The signal value of virtual flexible beam prediction replaces sampled value, on the contrary then retain sampled value.
In reference picture 2, Fig. 3, a kind of embodiment of the invention, constraints structure is used as according to several sampled values Need to calculate displacement and the slope of virtual flexible beam when building virtual flexible beam.For purposes of illustration only, adopting with 4 sampling instants here It is illustrated exemplified by sample value.
Due to needing the sampled value of continuous 3 sampling instants before using sensor current time sampled value, therefore in journey In sequence firstly the need of open up a length be 3 queue, the sampled value for 3 moment before storage.In Fig. 2, p4It is current The sampled value at moment, p1To p3The sampled value of continuous 3 sampling instants before being respectively current time sampled value.If in the presence of one Virtual virtual flexible beam, its bending rigidity is EI, has a rigidity to be k in each sampling instantsVirtual linear spring Sampled value and virtual flexible beam are connected, so beam can twist and flexural deformation in the presence of virtual linear spring, use y tables Show.
In the state of the equilibrium, suffered shearing force is beam in the x-direction
(1)
Because the low order end of beam is free, therefore work as x>Shearing force when 3 on beam is zero, i.e.,
(2)
These active forces equally can produce moment of torsion to beam, therefore the moment of flexure that beam is subject to is
(3)
Because the rightmost side of beam is free, therefore in the position of x=3, the moment of flexure of beam should be zero, therefore can draw
(4)
Assuming that the bending rigidity of beam is EI, then the moment of flexure of beam and the relation of deformation are
(5)
Therefore the slope of beam can in the x-direction be integrated by the moment of flexure being distributed on beam obtains, i.e.,
(6)
Here c1It is a undetermined constant.The displacement of beam can be by obtaining to the further integration of above formula:
(7)
Here c2It is another undetermined constant.After the displacement for obtaining beam, the drawing that each virtual linear spring is produced to beam Power can just be calculated:
(8)
With reference to (2) formula and (4) formula, below equation group can be obtained with simultaneous:
(9)
Wherein,,.。
Solve equation, you can obtain unknown parameter f1~f4, and c1And c2.It can be obtained respectively from there through formula (6) and formula (7) The slope of each sampling instant and displacement on to beam.
For the virtual flexible beam that the sampled value based on more sampling instants is built as constraints, it can pass through The mode of recursion builds equation group (9).By observing A matrixes it can be found that the element of the matrix may be roughly divided into four pieces, a left side Upper 4 × 4 submatrix A11, the submatrix A of upper right 4 × 212, the submatrix A of lower-left 2 × 421, and the null matrix A of bottom right 2 × 222.For A11 For, its matrix element with 1/6 (k s /EI), 23/6(k s /EI), 33/6(k s /EI) ... form occur and arrange, A12With A21In element also change significantly rule.Similarly, the element in X and B also changes significantly rule.For example, for Sampling instant number be 5 as the virtual flexible beam of constraints for, A, X, B form directly can be write out by recursion:
,,
In specific scheme, sampling instant number is adjusted according to the type and purposes of sensor signal.Sampling instant number Mesh is more, and Lubricity is better, predicts more credible.But sampling instant is excessive, performance improves unobvious, can increase computing on the contrary Amount.When sampling sensor is that the gap sensor that the current sensor and factor data communications errors polluted by impulsive noise triggers is Example, the number of sampling instant can be taken between 4 ~ 15.
Solved by formula (9) after X, it is possible to which the oblique of each sampling instant on beam can be respectively obtained by formula (6) and formula (7) Rate and displacement.So as to build virtual flexible beam as constraints according to several sampled values, optimum fit curve is formed.
Further, in the present invention virtual flexible beam bending rigidityEIWith virtual linear spring ratek s It is adjustable, with It is adapted to the sensor signal of different frequency bands scope.EIWithk s Ratio it is bigger, the deformation of virtual flexible beam is smaller, and smoothing effect is also got over It is good, now it is adapted to the signal of narrower bandwidth;Otherwise the signal of broader bandwidth can be adapted to.
After virtual flexible beam is in poised state, next sampling instant can be predicted according to the extended position of virtual flexible beam Sampled value:For sampling instant is the virtual flexible beam of N number of constraints, the slope on beam at low order end sampling instant It is respectively with displacement
, (10)
; (11)
For real system as magnetic-levitation train, the signal that sensor is collected has a continuity, therefore for Future time instance h (h > N 1), the predicted value of signal is
(12)
For example, for 4 sampling instants shown in Fig. 3 are as the virtual flexible beam of constraints, in next sampling Carve (x=4) predicted value be
(13)
The following possible value of a period of time sensor signal can be predicted with virtual flexible beam, if in follow-up sampling In, the data of sensor collection and the difference of predicted value exceed the threshold value set, then may determine that the signal that sensor is collected It is probably false, is now accomplished by giving up these sampled values, and sampled value is replaced with predicted value.
In reference picture 4 and Fig. 5, a kind of embodiment, adjusted according to the type of sensor signal and different purposes Whole sampling time interval.Below using detection sensor as magnetic-levitation train gap sensor, when magnetic-levitation train gap sensor passes through Exemplified by track seam, come the selection of the time interval that further illustrates sampling instant.Pass through seam for correcting gap sensor When the spurious signal that produces for, the interval between sampling instant should accordingly increase.For example, for 20mm track seam and 30mm diameters gap sensor probe for, when magnetic-levitation train with 18km/h (5m/s) speed by the seam when, gap The time span of sensor probe output spurious signal is (0.02+0.03 × 2)/5=0.016s, as shown in Figure 4.Now adopt The time interval at sample moment can take 16 ~ 20ms, and this can be sampled by the sampled data to sensor and obtained.This Under situation, because virtual flexible beam is continuous, its value extended in estimation range is also continuous, therefore in these scopes Interior sampled value can be replaced by these predicted values.In general, " projection " being superimposed upon in normal signal is extended pre- Survey signal " trimming " to fall, as shown in Figure 5.For the higher situation of speed, sensor is also shorter by the time of seam, now Between be spaced and equally can effectively judge spurious signal and be corrected.By being connect to magnetic-levitation train gap sensor by track Spurious signal correction during seam.It can make magnetic-levitation train more steady when by track seam, reduce electromagnet and collide with rail The probability in road.
In another embodiment, when sensor is that the current sensor and factor data polluted by impulsive noise leads to Believe the gap sensor that error code triggers, burr typically only lasts for 1 ~ 2 controlling cycle to both sensors in short-term, when now sampling Interval between quarter can be taken as consistent with suspension controller controlling cycle.
A kind of magnetic-levitation train sensor signal processing method provided by the present invention is described in detail above.Herein In apply specific case the principle and embodiment of the present invention be set forth, the explanation of above example is only intended to side The core concept of the assistant solution present invention.It should be pointed out that for those skilled in the art, not departing from this hair On the premise of bright principle, some improvement and modification can also be carried out to the present invention, these are improved and modification also falls into power of the present invention In the protection domain that profit is required.

Claims (7)

1. a kind of magnetic-levitation train sensor signal processing method, it is characterised in that the described method comprises the following steps:
Step 1:Obtain sensor current time sampled value and before current time continuous some sampling instants sampled value, general Sampled value builds virtual flexible beam as constraints, forms optimum fit curve;
Wherein, the sampled value is as the virtual flexible beam of constraints structure:
Sampled value is connected on the longitudinal axis between virtual flexible beam by virtual linear spring, and virtual linear spring is according to sampled value Produce pulling force effect on virtual flexible beam with the position difference of virtual flexible beam so that virtual flexible beam be subjected to displacement with it is curved Song, multiple sampled values are acted on simultaneously, the stress and moment of flexure of beam is reached balance, are determined that deflection of beam is deformed, are built in balance Virtual flexible equalizer bar under state;
Step 2:The sampled value of next sampling instant is predicted according to the extended position of virtual flexible beam;
Step 3:Next sampling instant sampled value is obtained, if the sampled value of next sampling instant and the difference of predicted value exceed setting Threshold value, then it is assumed that occur spurious signal in sensor, it is on the contrary then think that sensor signal is normal.
2. magnetic-levitation train sensor signal processing method according to claim 1, it is characterised in that when the sensor goes out During existing spurious signal, sampled value is replaced with predicted value, it is on the contrary then retain sampled value.
3. magnetic-levitation train sensor signal processing method according to claim 1, it is characterised in that virtually flexible beam is anti- Curved rigidityEIWith virtual linear spring ratek s It is adjustable, to be adapted to the sensor signal of different frequency bands scope.
4. magnetic-levitation train sensor signal processing method according to claim 1, it is characterised in that the sampling instant Several types and purposes according to sensor is determined.
5. magnetic-levitation train sensor signal processing method according to claim 1, it is characterised in that the neighbouring sample moment it Between time interval be identical in once virtual flexible beam building process.
6. magnetic-levitation train sensor signal processing method according to claim 4, it is characterised in that the sampling instant it Between time interval according to the type and purposes of sensor determine.
7. magnetic-levitation train sensor signal processing method according to claim 1, it is characterised in that the prediction is next to adopt The sampling interval at sample moment, more than or equal to the sampling interval of sampling instant.
CN201510444147.4A 2015-07-27 2015-07-27 A kind of magnetic-levitation train sensor signal processing method Active CN104990717B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510444147.4A CN104990717B (en) 2015-07-27 2015-07-27 A kind of magnetic-levitation train sensor signal processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510444147.4A CN104990717B (en) 2015-07-27 2015-07-27 A kind of magnetic-levitation train sensor signal processing method

Publications (2)

Publication Number Publication Date
CN104990717A CN104990717A (en) 2015-10-21
CN104990717B true CN104990717B (en) 2017-10-20

Family

ID=54302559

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510444147.4A Active CN104990717B (en) 2015-07-27 2015-07-27 A kind of magnetic-levitation train sensor signal processing method

Country Status (1)

Country Link
CN (1) CN104990717B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105978848A (en) * 2015-12-04 2016-09-28 乐视致新电子科技(天津)有限公司 Processing method and device for collection of sensor data
CN108657014B (en) * 2018-07-11 2019-12-31 中国人民解放军国防科技大学 Method and device for processing position sensor signal of suspension system of magnetic-levitation train
CN113525098B (en) * 2020-04-20 2023-02-03 株洲中车时代电气股份有限公司 Suspension control method and device for magnetic suspension vehicle
CN113092143B (en) * 2021-04-30 2023-01-31 中车青岛四方机车车辆股份有限公司 Detection system for suspension system of maglev train
CN113819959B (en) * 2021-11-24 2022-02-08 中国空气动力研究与发展中心设备设计与测试技术研究所 Suspension system anomaly detection method based on Hailinge distance and correlation coefficient

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100430260C (en) * 2006-08-28 2008-11-05 中国人民解放军国防科学技术大学 Suspension method for controling module of electromagnetism type magnetic suspension train in normal conduction and low temperature
CN100497037C (en) * 2007-10-17 2009-06-10 中国人民解放军国防科学技术大学 Method and system for preventing magnet levitation train from adsorption for track
CN103522913B (en) * 2013-10-17 2015-10-28 中国人民解放军国防科学技术大学 For suspension control method and the device of EMS type low speed aerotrain
CN104477049B (en) * 2014-11-28 2016-08-17 中国人民解放军国防科学技术大学 Magnetic-levitation train based on virtual energy capture device-bridge self-excited vibration suppressing method
CN104477048B (en) * 2014-11-28 2016-06-08 中国人民解放军国防科学技术大学 A kind of electromagnetic type is often led the suspension control method of low-speed maglev train

Also Published As

Publication number Publication date
CN104990717A (en) 2015-10-21

Similar Documents

Publication Publication Date Title
CN104990717B (en) A kind of magnetic-levitation train sensor signal processing method
CN107745654B (en) Method and device for processing signals of relative positioning sensor of magnetic-levitation train
Lee et al. Analysis of dynamic interaction between catenary and pantograph with experimental verification and performance evaluation in new high-speed line
EP2409640A1 (en) Biological parameter monitoring method, computer program, and biological parameter monitoring device
US20170098127A1 (en) Measuring device, measuring system, measuring method, and program
CN110319990B (en) Bridge dynamic deflection monitoring method based on inclinometer optimized arrangement
CN104598753A (en) Bridge moving vehicle load recognition method based on Brakhage V method
CN111169476B (en) Motion trend prediction method and device, controller and automobile
CN114559988A (en) Method, device and system for measuring speed, mileage, station and motion state of train
CN113272659A (en) Embedded system for vibration detection and analysis
CN104457643A (en) Impulse noise filtering method and device for track geometry detection data
JP2009018641A (en) Vibration damping device for railroad vehicle
CN108268027B (en) Driving track optimization method and system
CN105222744A (en) For producing the apparatus and method of trigger pip and position measurement apparatus related to this in position measurement apparatus
JP5184306B2 (en) Estimating the height of overhead lines in electric railways
CN113655733A (en) EMC (electromagnetic compatibility) semi-physical simulation method for magnetic field of axle counter of rail transit vehicle
KR101040511B1 (en) Method for synchronizing positions of track irregularity data measured from railway track, and system for the same
CN111637991B (en) Steel rail temperature stress detection method and terminal equipment
CN109631800B (en) Method and device for detecting dynamic lifting amount of contact line
CN104102831A (en) Interpolation sampled value protection method
Kim An experimental study of the dynamic characteristics of the catenary-pantograph interface in high speed trains
CN104266590B (en) Simple linear positioning system
CN110763989B (en) Method for predicting service life of relay of urban rail transit vehicle
CN113212499A (en) Real-time speed measurement method and system at track seam crossing moment by utilizing gap sensor
CN110852006A (en) Method for calculating torque compensation of pantograph actuating motor based on wind pressure sensing

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