CN102570979B - Iterative Teager energy operator demodulation method and system - Google Patents

Abstract

The invention discloses an iterative Teager energy operator (TEO) demodulation method and an iterative TEO demodulation system. The demodulation method comprises the following steps of: acquiring an original signal; performing energy operator demodulation to obtain a signal energy function and a signal differential energy function of the signal, and performing energy separation to obtain an amplitude envelope and instantaneous frequency; performing low-pass filtration to obtain an instantaneous amplitude envelope and the instantaneous frequency of the signal; subtracting the calculated instantaneous amplitude envelope from the amplitude envelope to obtain a residual signal component; and comparing the energy of the residual signal component with the energy of the original signal to obtain an energy difference, continuing the energy operator demodulation if the energy difference is greater than a preset value, otherwise finishing the energy operator demodulation. The iterative TEO demodulation can be implemented through a hardware platform consisting of a digital circuit and an analogue circuit or through a software algorithm on the basis of combination of an energy demodulation algorithm, zero-phase low-pass filtration and energy judgment, and the finally obtained demodulation signal component can be stored and displayed by a selected personal computer (PC) according to the type of the signal or can be displayed in real time by an analogue oscilloscope.

Description

A kind of iteration Teager energy operator demodulation method and system

Technical field

The present invention relates to a kind of method and system of signal processing, particularly a kind of iteration Teager energy operator demodulation method and system that is applicable to AM-FM signal.

Background technology

In the time there is local fault in the rotating machinery element in mechanical system, Mechanical Fault Vibration Signals often shows as modulation format, as rolling bearing fault vibration signal, gearbox fault vibration signal, Rotor Rubbing Fault vibration signal, in the envelope of these modulation signals and phase place, comprised abundant fault message, therefore demodulation analysis is the effective method of one of extracting fault signature and tracing trouble type.At present, in mechanical fault diagnosis field, conventional demodulation method has: Hilbert conversion demodulation, energy operator separation algorithm (Energy operator separation algorithm, EOSA), high pass absolute value are analyzed demodulation, detection filter demodulation, square solution mediation cyclostationary analysis etc.But these methods are only applicable to the AM-FM signal of simple component, and there is some limitation.For example, Hilbert converts at the modulation signal two ends that demodulate and the middle part that has sudden change will produce modulation, causes that demodulating error increases; The result curve that energy operator demodulation method obtains is not very smooth; High pass absolute value is analyzed the broad sense detection filter demodulation analysis methods such as demodulation, detection filter demodulation, square demodulation and is being taken absolute value or detection process may produce mixing effect, thereby in demodulation spectra, occurs false frequency content.Because most of Mechanical Fault Vibration Signals are all multicomponent AM-FM signals, therefore traditional method is first by bandpass filtering, it to be decomposed into the AM-FM signal of simple component, and then carries out demodulation in order to upper method.But the centre frequency of band pass filter and mainly dependence experience of bandwidth are selected, this will bring demodulating error, so that can not correctly extract fault signature.For this problem, propose Hilbert-Huang conversion, many components demodulation method such as TEO demodulation method, iteration Hilbert conversion based on empirical mode decomposition (EMD), and be applied in the failure diagnosis of rotating machinery.Wherein, Hilbert-Huang conversion and the TEO demodulation method based on EMD have good adaptivity, but all have the problems such as mode aliasing, boundary effect; Though and iteration Hilbert conversion can partly overcome the above problems, adaptivity is not strong and frequency demodulation precision is not high.

Therefore be badly in need of the strong and high signal iterative demodulation method and system of frequency demodulation precision of a kind of adaptivity.

Summary of the invention

In view of this, technical problem to be solved by this invention is to provide a kind of adaptivity by force and the high signal iterative demodulation method and system, particularly frequency demodulation precision aspect of frequency demodulation precision, is obviously better than the Hilbert conversion demodulation of iteration.

One of object of the present invention is to propose a kind of iteration Teager energy operator demodulation method; Two of object of the present invention is to propose a kind of iteration Teager energy operator demodulating system.

One of object of the present invention is achieved through the following technical solutions:

A kind of iteration Teager energy operator demodulation method provided by the invention, comprises the following steps:

S1: gather original signal;

S2: signal is carried out to energy operator demodulation, obtain signal energy function and the signal differentiation energy function of this signal;

S3: signal energy function and signal differentiation energy function are carried out to energy separation, obtain amplitude envelope and instantaneous frequency;

S4: amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtain instantaneous amplitude envelope and the instantaneous frequency of signal;

S5: amplitude envelope is deducted to the instantaneous amplitude envelope of obtaining, obtain residual signal component;

S6: residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation;

S7: in the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

Further, further comprising the steps of:

S8: demonstration and storage obtain instantaneous amplitude, instantaneous frequency and the residual signal component of signal each time.

Further, in described step S2, signal is carried out also comprising for primary signal and differentiating before energy operator demodulation, obtains differential signal,

Further, the energy separation in described step S3, specifically comprises the following steps:

S31: primary signal is carried out to energy operator computing as follows, obtain signal energy function:

Signal energy operator function expression formula is:

Wherein, ψ (x) represents signal energy function, and x (t) represents primary signal, represent single order differential signal, represent second-order differential;

S32: differential signal is carried out to energy operator computing as follows, obtain signal differentiation energy function:

Differential signal energy operator function expression is:

Wherein, represent signal differentiation energy function; represent three rank differential signals;

S33: signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, and filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum;

S34: the instantaneous amplitude envelope and the instantaneous frequency that obtain signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Wherein, | a (t) | represent the instantaneous amplitude of signal, ω (t) represents instantaneous frequency.

Further, in described step S4, amplitude envelope and instantaneous frequency are carried out to low-pass filtering employing zero phase low-pass filtering.

Two of object of the present invention is achieved through the following technical solutions:

Iteration Teager energy operator demodulating system provided by the invention, comprises signal acquisition module, energy operator demodulation module, energy separation module, low-pass filtering module, residual signal component module and energy judging module,

Described signal acquisition module, for gathering primary signal;

Described energy operator demodulation module, for signal is carried out to energy operator demodulation, obtains signal energy function and the signal differentiation energy function of this signal;

Described energy separation module, for signal energy function and signal differentiation energy function are carried out to energy separation, obtains amplitude envelope and instantaneous frequency;

Described low-pass filtering module, for amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtains instantaneous amplitude envelope and the instantaneous frequency of signal by following formula;

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Described residual signal component module, for amplitude envelope being deducted to the instantaneous amplitude envelope of obtaining, obtains residual signal component;

Described energy judging module, for residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation; In the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

Further, also comprise display module, for showing and storing the instantaneous amplitude, instantaneous frequency and the residual signal component that obtain each time signal.

Further, described energy operator demodulation module also comprises signal energy operator operation module, signal differentiation energy operator computing module, low-pass filtering module;

Described signal energy operator operation module, for primary signal is carried out to energy operator computing as follows, obtains signal energy function: energy operator function expression is:

Further, described energy operator demodulation module also comprises signal differentiation module, for primary signal is differentiated, and exports differential signal;

Described signal differentiation energy operator computing module, for differential signal is carried out to energy operator computing as follows, obtains signal differentiation energy function: differential signal energy operator function expression is:

Described low-pass filtering module, for signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum.

Further, described energy separation module also comprises instantaneous amplitude envelope computing module and instantaneous frequency computing module,

Described instantaneous amplitude envelope computing module, for obtain the instantaneous amplitude envelope of signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},$

Described instantaneous frequency computing module, for obtain the instantaneous frequency of signal by following formula:

$\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Further, also comprise zero phase low-pass filtering module, described zero phase low-pass filtering module, for carrying out low-pass filtering to instantaneous amplitude envelope and instantaneous frequency.

The invention has the advantages that: the present invention adopts the many components demodulation method based on Teager energy operator (TEO), the energy operator that the differential of the primary signal collecting and signal is realized by analog circuit or software calculates, obtain signal energy function and the signal differentiation energy function of this signal, respectively the signal energy function obtaining and signal differentiation energy function are carried out to low-pass filtering, filtering high frequency error component.Then obtain amplitude envelope and the instantaneous frequency of many component signals by energy separation; Carry out zero phase low-pass filtering for the amplitude envelope and the instantaneous frequency that produce, so just obtain instantaneous amplitude and the instantaneous frequency of demodulation for the first time; Amplitude envelope is deducted to a certain component instantaneous amplitude of obtaining and obtain remaining signal, it is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, transfer to the input of energy operator algorithm using this signal as iteration signal, again carry out demodulation.In the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete demodulating process, be called the demodulation of iteration Teager energy operator.Wherein, residual signal is called residue signal component; The signal amplitude envelope and the instantaneous frequency that obtain each time, and residue signal component can gather to PC and be shown by capture card, if analog signal also can directly access analog oscilloscope and show in real time.

The present invention is based on energy demodulating algorithm, zero phase low-pass filtering combines with energy judgement, realizes the demodulation of iteration Teager energy operator.The hardware platform that its realization can be built by digital circuit and analog circuit is realized, and also can realize software algorithm by programming.The restituted signal component finally obtaining can basis signal the preservation of type selecting PC and demonstration or analog oscilloscope show in real time.

Other advantage of the present invention, target and feature will be set forth to a certain extent in the following description, and to a certain extent, based on will be apparent to those skilled in the art to investigating below, or can be instructed from the practice of the present invention.The objects and other advantages of the present invention can be passed through specification below, claims, and in accompanying drawing, specifically noted structure realizes and obtains.

Brief description of the drawings

In order to make the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing, the present invention is described in further detail, wherein:

The flow chart of the iteration Teager energy operator demodulation method that Fig. 1 provides for the embodiment of the present invention;

The circuit diagram of realizing iteration Teager energy operator demodulation method that Fig. 2 provides for the embodiment of the present invention;

The signal that Fig. 3 provides for the embodiment of the present invention many component amplitudes time diagram that superposes;

The signal separation amplitude time diagram that Fig. 4 provides for the embodiment of the present invention;

The signal cross frequence time diagram that Fig. 5 provides for the embodiment of the present invention.

Embodiment

Below with reference to accompanying drawing, the preferred embodiments of the present invention are described in detail; Should be appreciated that preferred embodiment is only for the present invention is described, instead of in order to limit the scope of the invention.

The flow chart of the iteration Teager energy operator demodulation method that Fig. 1 provides for the embodiment of the present invention, the circuit diagram of realizing iteration Teager energy operator demodulation method that Fig. 2 provides for the embodiment of the present invention, the signal that Fig. 3 provides for the embodiment of the present invention many component amplitudes frequency diagram that superposes, the signal separation amplitude frequency diagram that Fig. 4 provides for the embodiment of the present invention, the signal cross frequence time diagram that Fig. 5 provides for the embodiment of the present invention, as shown in the figure: a kind of iteration Teager energy operator demodulation method provided by the invention, comprises the following steps:

S1: gather original signal;

S2: primary signal is differentiated, obtain differential signal, then signal is carried out to energy operator demodulation, obtain signal energy function and the signal differentiation energy function of this signal;

S3: signal energy function and signal differentiation energy function are carried out to energy separation, obtain amplitude envelope and instantaneous frequency, specifically comprise the following steps:

S31: primary signal is carried out to energy operator computing as follows, obtain signal energy function:

Signal energy operator function expression formula is:

Wherein, ψ (x) represents signal energy function, and x (t) represents primary signal, represent single order differential signal, represent second-order differential;

S32: differential signal is carried out to energy operator computing as follows, obtain signal differentiation energy function:

Differential signal energy operator function expression is:

Wherein, represent signal differentiation energy function; represent three rank differential signals;

S33: signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, and filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum;

S34: the instantaneous amplitude envelope and the instantaneous frequency that obtain signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Wherein, | a (t) | represent the instantaneous amplitude of signal, ω (t) represents instantaneous frequency.

S4: amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtain instantaneous amplitude envelope and the instantaneous frequency of signal; In the embodiment of the present invention, amplitude envelope and instantaneous frequency are carried out to low-pass filtering employing zero phase low-pass filtering.

S5: amplitude envelope is deducted to the instantaneous amplitude envelope of obtaining, obtain residual signal component;

S6: residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation;

S7: in the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

S8: demonstration and storage obtain instantaneous amplitude, instantaneous frequency and the residual signal component of signal each time.

The embodiment of the present invention also provides a kind of iteration Teager energy operator demodulating system, comprises signal acquisition module, energy operator demodulation module, energy separation module, low-pass filtering module, residual signal component module, energy judging module, zero phase low-pass filtering module and display module;

Described signal acquisition module, for gathering primary signal;

Described energy operator demodulation module, for signal is carried out to energy operator demodulation, obtains signal energy function and the signal differentiation energy function of this signal;

Described energy separation module, for signal energy function and signal differentiation energy function are carried out to energy separation, obtains amplitude envelope and instantaneous frequency;

Described low-pass filtering module, for amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtains instantaneous amplitude envelope and the instantaneous frequency of signal by following formula;

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Described residual signal component module, for amplitude envelope being deducted to the instantaneous amplitude envelope of obtaining, obtains residual signal component;

Described energy judging module, for residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation; In the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

Described zero phase low-pass filtering module, for instantaneous amplitude envelope and instantaneous frequency are carried out to low-pass filtering,

Described display module, for showing and storing the instantaneous amplitude, instantaneous frequency and the residual signal component that obtain each time signal.

Described energy operator demodulation module also comprises signal energy operator operation module, signal differentiation energy operator computing module, low-pass filtering module;

Described signal energy operator operation module, for primary signal is carried out to energy operator computing as follows, obtains signal energy function: energy operator function expression is:

Described energy operator demodulation module also comprises signal differentiation module, for primary signal is differentiated, and exports differential signal;

Described signal differentiation energy operator computing module, for differential signal is carried out to energy operator computing as follows, obtains signal differentiation energy function: differential signal energy operator function expression is:

Described low-pass filtering module, for signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum.

Described energy separation module also comprises instantaneous amplitude envelope computing module and instantaneous frequency computing module,

Described instantaneous amplitude envelope computing module, for obtain the instantaneous amplitude envelope of signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},$

Described instantaneous frequency computing module, for obtain the instantaneous frequency of signal by following formula:

$\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Elaborate the method that the invention provides the iteration Teager energy operator demodulation that embodiment provides below, be implemented as follows:

Be analog signal x (t) or digital signal x (n) for primary signal, first it differentiated obtain differential signal or then primary signal and differential signal are carried out respectively to energy operator ψ computing, obtain respectively signal energy function and signal differentiation energy function ψ (x) and wherein energy operator function expression is for AF-AM signal arbitrarily:

X (t)=a (t) cos[φ (t)], can obtain:

$\mathrm{\ψ}\left(a\right)={\left[\stackrel{\·}{a}\left(t\right)\right]}^{2}a\left(t\right)\stackrel{\·\·}{a}\left(t\right),$

For AF-AM signal arbitrarily, the variation of modulation signal changes much slow than carrier wave.

Wherein, a (t) and for the variation of carrier wave, be slowly, what therefore can be similar to is considered as constant it, so just has: ψ (a) ≈ 0, $\stackrel{\·\·}{\mathrm{\φ}}\left(t\right)\≈0.$ So obtain $\mathrm{\ψ}\left(x\right)\≈{\left[a\left(t\right)\stackrel{\·}{\mathrm{\φ}}\left(t\right)\right]}^2=a^2\left(t\right){\mathrm{\ω}}^2\left(t\right),$ Can obtain equally $\mathrm{\ψ}\left(\stackrel{\·}{x}\right)\≈a^2\left(t\right){\mathrm{\ω}}^4\left(t\right).$

There is error in above-mentioned simplification process:

$\begin{array}{c}E=a^2\left(t\right)\stackrel{\·\·}{\mathrm{\φ}}\left(t\right)\sin \left[2\mathrm{\φ}\left(t\right)\right]/2+{\cos }^2\left[\mathrm{\φ}\left(t\right)\right]\mathrm{\ψ}\left(a\right)\\ =\frac{1}{2}\mathrm{\ψ}\left[a\left(t\right)\right]+\frac{1}{2}{a}^2\left(t\right)\stackrel{\·\·}{\mathrm{\φ}}\left(t\right)\sin \left[2\mathrm{\φ}\left(t\right)\right]+\mathrm{\ψ}\left[a\left(t\right)\right]\cos \left[2\mathrm{\φ}\left(t\right)\right]\end{array},$

Wherein low frequency aberration component is:

$E_L=\frac{1}{2}\mathrm{\ψ}\left[a\left(t\right)\right],$

High frequency error component is:

$E_H=\frac{1}{2}{a}^2\left(t\right)\stackrel{\·\·}{\mathrm{\φ}}\left(t\right)\sin \left[2\mathrm{\φ}\left(t\right)\right]+\mathrm{\ψ}\left[a\left(t\right)\right]\cos \left[2\mathrm{\φ}\left(t\right)\right].$

Analyzing and producing error is mainly due to high frequency error component E _h, therefore adopting low-pass filtering, filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum.

Signal energy function and signal differentiation energy function are carried out to energy separation, obtain instantaneous amplitude and the instantaneous frequency of signal,

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

For a better understanding of the present invention, below, be described by specific embodiment.

The specific implementation process that the invention provides embodiment is as follows:

1. first generate two simple component real signal x by matlab program ₁and x (t) ₂(t), these two signals are all AM-FM signals, and their expression formula is as follows:

x ₁(t)＝[1+0.5cos(2π*20t)]sin[2π*500t+cos(2π*60t)]

x ₂(t)＝[0.2+0.2cos(2π*10t)]sin[2π*600t+cos(2π*60t)]。

2. by the stack of two signals, simulate multicomponent AM-FM signal x (t)=x ₁(t)+x ₂(t), as shown in Figure 3.

3. adopt matlab software programming iteration Teager energy operator demodulation method, isolate x ₁and x (t) ₂(t) amplitude A ₁and A (t) ₂, and their frequency f (t) ₁and f (t) ₂(t).The amplitude figure wherein demodulating as 4 and Fig. 5 as shown in.

4. effect analysis: in original AM-FM signal, theoretical instantaneous amplitude is A ₁(t)=1+0.5cos (2 π * 20t), A ₂(t)=0.2+0.2cos (2 π * 10t), instantaneous frequency f ₁(t)=500-60sin (2 π 60t), f ₂(t)=600-60sin (2 π 60t).As can be seen from Figure 4 and Figure 5, the A of the acquisition of process iteration Teager energy operator demodulation method ₁(t), A ₂(t), f ₁and f (t) ₂(t) very approaching with theoretical value.The more important thing is, the method can be carried out adaptive demodulation analysis, and manual intervention is few.

The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, obviously, those skilled in the art can carry out various changes and modification and not depart from the spirit and scope of the present invention the present invention.Like this, if these amendments of the present invention and within modification belongs to the scope of the claims in the present invention and equivalent technologies thereof, the present invention is also intended to comprise these changes and modification interior.

Claims (10)

1. an iteration Teager energy operator demodulation method, is characterized in that: comprise the following steps:

S1: gather original signal;

S2: signal is carried out to energy operator demodulation, obtain signal energy function and the signal differentiation energy function of this signal;

S3: signal energy function and signal differentiation energy function are carried out to energy separation, obtain amplitude envelope and instantaneous frequency;

S4: amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtain instantaneous amplitude envelope and the instantaneous frequency of signal;

S5: amplitude envelope is deducted to the instantaneous amplitude envelope of obtaining, obtain residual signal component;

S6: residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation;

S7: in the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

2. iteration Teager energy operator demodulation method according to claim 1, is characterized in that: further comprising the steps of:

S8: demonstration and storage obtain instantaneous amplitude, instantaneous frequency and the residual signal component of signal each time.

3. iteration Teager energy operator demodulation method according to claim 1, is characterized in that: in described step S2, signal is carried out also comprising for primary signal and differentiating before energy operator demodulation, obtain differential signal.

4. iteration Teager energy operator demodulation method according to claim 1, is characterized in that: the energy separation in described step S3, specifically comprises the following steps:

S31: primary signal is carried out to energy operator computing as follows, obtain signal energy function:

Signal energy operator function expression formula is:

Wherein, ψ (x) represents signal energy function, and x (t) represents primary signal, represent single order differential signal, represent second-order differential;

S32: differential signal is carried out to energy operator computing as follows, obtain signal differentiation energy function:

Differential signal energy operator function expression is:

Wherein, represent signal differentiation energy function; represent three rank differential signals;

S33: signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, and filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum;

S34: the instantaneous amplitude envelope and the instantaneous frequency that obtain signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Wherein, | a (t) | represent the instantaneous amplitude of signal, ω (t) represents instantaneous frequency.

5. iteration Teager energy operator demodulation method according to claim 1, is characterized in that: in described step S4, amplitude envelope and instantaneous frequency are carried out to low-pass filtering employing zero phase low-pass filtering.

6. an iteration Teager energy operator demodulating system, is characterized in that: comprise signal acquisition module, energy operator demodulation module, energy separation module, low-pass filtering module, residual signal component module and energy judging module,

Described signal acquisition module, for gathering primary signal;

Described energy operator demodulation module, for signal is carried out to energy operator demodulation, obtains signal energy function and the signal differentiation energy function of this signal;

Described energy separation module, for signal energy function and signal differentiation energy function are carried out to energy separation, obtains amplitude envelope and instantaneous frequency;

Described low-pass filtering module, for amplitude envelope and instantaneous frequency are carried out to low-pass filtering, obtains instantaneous amplitude envelope and the instantaneous frequency of signal by following formula;

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}},$

Described residual signal component module, for amplitude envelope being deducted to the instantaneous amplitude envelope of obtaining, obtains residual signal component;

Described energy judging module, for residual signal component is carried out to energy judgement, if the energy of residual signal component and the energy ratio of original signal are greater than predetermined value, return in S2 this residual signal component is carried out to energy operator demodulation; In the time that the energy of residual signal component and the energy ratio of original signal are less than predetermined value, complete iteration Teager energy operator demodulating process.

7. iteration Teager energy operator demodulating system according to claim 6, is characterized in that: also comprise display module, for showing and storing the instantaneous amplitude, instantaneous frequency and the residual signal component that obtain each time signal.

8. iteration Teager energy operator demodulating system according to claim 6, is characterized in that: described energy operator demodulation module also comprises signal differentiation module, for primary signal is differentiated, and exports differential signal.

9. iteration Teager energy operator demodulating system according to claim 6, it is characterized in that: described energy separation module also comprises signal energy separation module, signal differentiation energy module, low-pass filtering module, instantaneous amplitude envelope module and instantaneous frequency module, described signal energy separation module, for primary signal is carried out to energy operator computing as follows, obtain signal energy function: energy operator function expression is:

Described signal differentiation energy separation module, for differential signal is carried out to energy operator computing as follows, obtains signal differentiation energy function: differential signal energy operator function expression is:

Described low-pass filtering module, for signal energy function and signal differentiation energy function are adopted to zero phase low-pass filtering, filtering high frequency error component, obtains signal energy function and the signal differentiation energy function of error minimum;

Instantaneous amplitude envelope module, for obtain the instantaneous amplitude envelope of signal by following formula:

$\left|a\left(t\right)\right|\≈\frac{\mathrm{\ψ}\left(x\right)}{\sqrt{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}},$

Described instantaneous frequency module, for obtain the instantaneous frequency of signal by following formula:

$\mathrm{\ω}\left(t\right)\≈\sqrt{\frac{\mathrm{\ψ}\left(\stackrel{\·}{x}\right)}{\mathrm{\ψ}\left(x\right)}}.$

10. iteration Teager energy operator demodulating system according to claim 6, is characterized in that: also comprise zero phase low-pass filtering module, described zero phase low-pass filtering module, for carrying out low-pass filtering to amplitude envelope and instantaneous frequency.