<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJSST</journal-id><journal-title-group><journal-title>Open Journal of Safety Science and Technology</journal-title></journal-title-group><issn pub-type="epub">2162-5999</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojsst.2015.52008</article-id><article-id pub-id-type="publisher-id">OJSST-57342</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject><subject> Social Sciences&amp;Humanities</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Research of Railway Line State Detection Signal Processing Method Based on EMD
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>en-Fa</surname><given-names>Zhu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hui-Zhen</surname><given-names>Ma</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiao-Dong</surname><given-names>Chai</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Shu-Bin</surname><given-names>Zhen</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Library, Shanghai University of Engineering Science, Shanghai, China</addr-line></aff><aff id="aff1"><addr-line>College of Urban Railway Transportation, Shanghai University of Engineering Science, Shanghai, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>zhuwenfa1986@163.com(EZ)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>15</day><month>05</month><year>2015</year></pub-date><volume>05</volume><issue>02</issue><fpage>63</fpage><lpage>68</lpage><history><date date-type="received"><day>3</day>	<month>April</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>18</month>	<year>June</year>	</date><date date-type="accepted"><day>24</day>	<month>June</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper, an EMD de-noising algorithm is proposed based on the statistical feature of random noise, which can eliminate the noise impaction digital integrator generated by the collected railway line state detection signals using strap-down inertial technology. Firstly, the first IMF component of the noise-dominant modes treated by the process “random sort-sum-average-reconstruc-tion”, the signal-to-noise ratio is improved while the noise power is weakened in this process. Then the signal-to-noise cut-off can be determined according to the characters of noise autocorrelation function. Finally, the global threshold could be selected by the noise-dominant mode component, so as to realize the function of filtering. The simulation and validation based on the collected railway line acceleration data using the EMD de-noising algorithm show that the noise in railway line state acceleration detection signals can be effectively eliminated using this method.
 
</p></abstract><kwd-group><kwd>Track Status</kwd><kwd> EMD</kwd><kwd> Autocorrelation Function</kwd><kwd> Statistical Properties</kwd><kwd> Filtering</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In the use of Strap-down Inertial Navigation System [<xref ref-type="bibr" rid="scirp.57342-ref1">1</xref>] to detection track state, the collected acceleration signal should be have integrals two times to get the movement of the vector. The acceleration signal from the inertial measurement unit on the scene contains a certain amount of noise which is easy to cause the integrator saturation, so the acceleration signal effect of the reduction noise infects the accuracy of track long wave irregularity detection.</p></sec><sec id="s2"><title>2. Hilbert-Huang Transform</title><p>Hilbert Huang transform [<xref ref-type="bibr" rid="scirp.57342-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.57342-ref3">3</xref>] is a signal processing method proposed by Dr. Hung of the America NASA. it has obvious advantages in dealing with nonlinear and non-stationary signal. HHT mainly contains the following two parts: Empirical Mode Decomposition (EMD) and Hilbert transform, the empirical mode decomposition (EMD) is core of the HHT. According to the fluctuation of real size step by step signal in the signal decomposition (for signal smoothing processing), EMD get a series of different characteristics of the scale of IMF (Intrinsic Mode Function, the intrinsic mode function). The IMF should meet the following two conditions:</p><p>1) the IMF, an number of extreme points and zero points should be equal or the maximum difference of 1;</p><p>2) signal locally symmetric on the time axis. the any point of signal, the mean of the envelope curve point determined by the signal local maxima and the envelope point determination and local minimum value is 0.</p></sec><sec id="s3"><title>3. The EMD Denoising Algorithm Based on the Statistical Characteristic of Noise</title><sec id="s3_1"><title>3.1. Characteristics of the Random Noise Statistic</title><p>For discrete signal which length is N, the power is calculated as:</p><disp-formula id="scirp.57342-formula500"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1480126x5.png"  xlink:type="simple"/></disp-formula><p>If the original signal amplitude of each element was kept invariant, its position will be ordering up. According to the type 1, signal through random sort, its power remains unchanged. The first random noise of IMF component by EMD decomposition can still keep the random properties of approximation. However, power is obtained the first IMF component random sort the number which is after “random sort of cumulative average reconstruction” is inversely proportional, namely power decreases with the increase of random sort. as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. Therefore, the noise energy can be reduced through the signal with noise is decomposed imfi “EMD of random scheduling and cumulative average”.</p></sec><sec id="s3_2"><title>3.2. The Statistical Properties of the Random Noise from the Autocorrelation Function</title><p>The random signal of the autocorrelation function is defined as:</p><disp-formula id="scirp.57342-formula501"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1480126x7.png"  xlink:type="simple"/></disp-formula><p>Here the use of normalized autocorrelation function to represent the the correlation degree of random signal in different time values. That is:</p><disp-formula id="scirp.57342-formula502"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1480126x8.png"  xlink:type="simple"/></disp-formula><p>In the formula, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x10.png" xlink:type="simple"/></inline-formula></p><p>Although the normalized autocorrelation function of the random noise and signal get maximum value at zero, but the point of the time difference outside zero have following difference:</p><p>The random noise from the autocorrelation function have the maximum value at zero point. at the other point, the autocorrelation function decays rapidly decrease small on the contrary. but the general signal autocorrelation function does not have such a characteristic. According to the difference of random noise and general signal autocorrelation function can accurately judge the signal-to-noise cutoff point between the dominant mode of K.</p></sec><sec id="s3_3"><title>3.3. Denoising Algorithm of Constructing EMD</title><p>According to the two points of the statistical characteristics to random noise the above, the construction of a new denoising method-EMD denoising method based on the noise statistical characteristics [<xref ref-type="bibr" rid="scirp.57342-ref4">4</xref>] -[<xref ref-type="bibr" rid="scirp.57342-ref6">6</xref>] , the specific steps are as follows in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>1) After the signal with noise is decomposed by EMD, the N of intrinsic mode components in IMF were got. the decomposed the last trend component is recorded as the first N IMF.</p><p>2) Record <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x12.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x13.png" xlink:type="simple"/></inline-formula></p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Algorithm flow chart</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x14.png"/></fig><p>3) A new component is obtained after the random ordering, and the cumulative,</p><disp-formula id="scirp.57342-formula503"><graphic  xlink:href="http://html.scirp.org/file/6-1480126x15.png"  xlink:type="simple"/></disp-formula><p>4) Repeat step 3 R times, a new power weakened noise dominant mode was got after calculated accumulation</p><p>average. record <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x17.png" xlink:type="simple"/></inline-formula></p><p>5) Get a new SNR signal with noise is improved from the reconstruction. Record <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x19.png" xlink:type="simple"/></inline-formula></p><p>6) Received signal after inhibition will be regarded as the noise repeated S times steps after 1 - 5.</p><p>7) The EMD decomposition, calculate the autocorrelation function of N imf component, according to the characteristics of the autocorrelation function obtained signal-to-noise leading the dividing point between modes K</p><p>(8) Global threshold processing method selection to the dominant modal components of noise record</p><disp-formula id="scirp.57342-formula504"><graphic  xlink:href="http://html.scirp.org/file/6-1480126x21.png"  xlink:type="simple"/></disp-formula><p>In the formula, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x23.png" xlink:type="simple"/></inline-formula>as the components of the threshold, L is the length of the signal, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1480126x24.png" xlink:type="simple"/></inline-formula>is the standard deviation.</p><p>9) Obtain the signal after denoising, get the obtained and reconstructed.</p></sec><sec id="s3_4"><title>3.4. Simulation Analysis</title><p>Signal denoising by using the statistical characteristic of noise in EMD based de-noising method respectively, which is obtained by the superposition of signal after Gauss white noise after the EMD decomposition, take the first IMF component is “random scheduling and cumulative average” R times, then with the other IMF component reconstruction, continue to repeat S times (), the noise power is weakened, improved SNR signal were got, as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>On the EMD decomposition of IMF, calculation of the various components of the autocorrelation function x corr (imf), in <xref ref-type="fig" rid="fig3">Figure 3</xref>. According to the random noise of the autocorrelation function to the statistical proper-</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Signal to noise ratio improvement after SNR</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x25.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Each component of the normalized autocorrelation function</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x26.png"/></fig><p>ties, the signal-to-noise cutoff point can be determined. Therefore the global threshold selection method for denoising is obtained, after reconstruction, the signal after denoising was got, as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p></sec></sec><sec id="s4"><title>4. Experimental Verification</title><p>In the experiment to test the car loaded with inertial measurement unit through a simulation of line section, collecting sensitive acceleration signal to trolley moving along X, Y, Z axis. First of all, using the method of mean filter DC quantity contained in the micro accelerometer signal. Using EMD de-noising method based on the statistical characteristic of noise, the acceleration signal is on denoising processing after eliminating direct flow. And the waveform and frequency spectrum of signals before and after respectively were got as shown in Figures 5-7. Above, denoising method based on noise statistical characteristics can remove random noise effectively from the acceleration signal.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> After denoising</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x27.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Waveform and frequency spectrum of X axis acceleration signal before and after denoising</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x28.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Waveform and frequency spectrum of Y axis acceleration signal before and after denoising</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x29.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Waveform and frequency spectrum of Z axis acceleration signal before and after denoising</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1480126x30.png"/></fig></sec><sec id="s5"><title>5. Conclusion</title><p>This article is on the correlation function of studying the statistical properties of random noise. EMD filtering method is proposed based on the statistical characteristics. Using this method the acceleration signal exported an inertial measurement unit to deal with denoising. The acceleration signal denoised is integral operated. Through the attitude matrix the motion of the information is converted from the carrier coordinates system to the geographic coordinate system. The car trajectory can be got [<xref ref-type="bibr" rid="scirp.57342-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.57342-ref8">8</xref>] . Experiments show: EMD filtering method based on the statistical characteristics of the noise can remove the noise signal effectively. The requirements of detection accuracy can be met better.</p></sec><sec id="s6"><title>Acknowledgements</title><p>The project is jointly supported by the Shanghai Tertiary Education Specialized Fund for Planning to Support Young Teacher’s Trainings (ZZGJD13072), National Natural Science Foundation of China (51478258), and the Shanghai Graduate Education Innovation Project in Layout and Construction Project (13sc002).</p></sec></body><back><ref-list><title>References</title><ref id="scirp.57342-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Zheng, S.B., Lin, J.H. and Lin, G.B. (2007) Implementation of Detecting Maglev Track Long Wave Irregularity Based on Inertial Measurement Principle. Journal of Electronic Measurement and Instrument, 21, 61-65.</mixed-citation></ref><ref id="scirp.57342-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Huang, N.E., Shen, Z., Long, S.R., Wu, M.C., Shih, H.H., Zheng, Q., et al. (1998) The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-Stationary Time Series Analysis. Proceedings of the Royal Society of London A, 454, 56-78.</mixed-citation></ref><ref id="scirp.57342-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Rilling, G. and Flandfin, P. (2008) One or Two Frequencies? The Empirical Mode Decomposition Answers. IEEE Transactions on Signal Processing, 56, 85-96. http://dx.doi.org/10.1109/TSP.2007.906771</mixed-citation></ref><ref id="scirp.57342-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Boudraa, A.-O. and Cexus, J.-C. (2007) EMD-Based Signal Filtering. IEEE Transactions on Instrumentation and Measurement, 56, 2196-2202. http://dx.doi.org/10.1109/TIM.2007.907967</mixed-citation></ref><ref id="scirp.57342-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Sun, J.-B., Zhu, Y.-L., Liu, L.-Q. and Huang, J.-C. (2010) Study on De-Noising Method Based on EMD Used in Insulator Leakage Current. Journal of North China Electric Power University, 37, 1-5, 22.</mixed-citation></ref><ref id="scirp.57342-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Wang, Z. (2010) Research on EMD Algorithm and Its Application in Signal Denoising. Harbin Engineering University, 34-36.</mixed-citation></ref><ref id="scirp.57342-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, W.F., Chai, X.D., Zheng, S.B., Li, L.M. and Luo, Y.J. (2012) Developing IMU SignaI Acquisition and Processing System Based on LabVIEW. Computer Measurement &amp; Control, 20, 1697-1698, 1706.</mixed-citation></ref><ref id="scirp.57342-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Liu, X.-C., Chai, X.D., Zheng, S.-B., Zhu, W.-F. and Jiang, Q.-Z. (2013) De-noising Method for Acceleration Signal Measured by Inertial Measurement Unit. Instrument Technique and Sensor, 10, 97-99.</mixed-citation></ref></ref-list></back></article>