CN109916625B - Single-channel gear box multi-fault separation dual-core micro-processing system - Google Patents

Single-channel gear box multi-fault separation dual-core micro-processing system Download PDF

Info

Publication number
CN109916625B
CN109916625B CN201910232499.1A CN201910232499A CN109916625B CN 109916625 B CN109916625 B CN 109916625B CN 201910232499 A CN201910232499 A CN 201910232499A CN 109916625 B CN109916625 B CN 109916625B
Authority
CN
China
Prior art keywords
signal
reconstruction
module
channel
separation
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
CN201910232499.1A
Other languages
Chinese (zh)
Other versions
CN109916625A (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.)
Changsha University of Science and Technology
Original Assignee
Changsha University of Science and 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 Changsha University of Science and Technology filed Critical Changsha University of Science and Technology
Priority to CN201910232499.1A priority Critical patent/CN109916625B/en
Publication of CN109916625A publication Critical patent/CN109916625A/en
Application granted granted Critical
Publication of CN109916625B publication Critical patent/CN109916625B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention relates to a single-channel gearbox multi-fault separation dual-core micro-processing system, and belongs to the technical field of gearbox fault diagnosis and signal processing. The invention relates to a dual-core micro-processing system of a DSP and an ARM, wherein a DSP core comprises a phase space reconstruction parameter estimation module, a phase space reconstruction module, a reconstruction parameter adjustment module, a reconstruction signal separation module and a Fourier transform module, and the ARM core comprises a time-frequency domain spectrogram drawing module; the invention adopts the time delay reconstruction method for expansion, thereby avoiding the human intervention of the reconstructed signal by taking experience as the judgment standard; the influence of noise and redundancy in the gearbox on a reconstructed signal is reduced by adopting a principal component analysis method; and estimating a source signal probability density function by utilizing a finite support sample kernel function and FastICA fusion algorithm, further obtaining a nonlinear function according with the statistical characteristics of the source signal, and finally realizing effective separation of the signals. The method greatly reduces the influence of improper selection of the reconstruction parameters and the separation nonlinear function on the multi-fault separation of the gearbox.

Description

Single-channel gear box multi-fault separation dual-core micro-processing system
Technical Field
The invention relates to the technical field of gearbox fault diagnosis and signal processing analysis, in particular to a single-channel gearbox multi-fault separation dual-core micro-processing system.
Background
The gear box is an important mechanism for transmitting power of mechanical equipment, and is one of important research objects in mechanical failure diagnosis. Therefore, the research on the fault diagnosis technology of the gearbox is significant for guaranteeing the safe operation of mechanical equipment. The difficulty of fault detection is complex faults, and the estimation of the number of complex faults is the difficulty of complex fault diagnosis. In order to understand the operating state of the gearbox, the acceleration sensor typically picks up its vibration signal. However, the sensor can only be installed on the outer surface of the gear box, which causes the transmission path of the vibration signal collected by the sensor to be too long and the vibration signal to be a composite of a plurality of vibration signals, including redundant vibration of each component, and besides the influence of redundant vibration of the gear box, the influence of different interference sources of various electronic noises from the collection channel and large noise of the surrounding environment also exists, which causes the number of the actually collected sources of each channel to be often more than the number of the sensor, which increases the difficulty of diagnosis, so that a plurality of sensors are generally installed in a laboratory to collect signals. However, in practical engineering problems, due to cost problems and environmental problems, multiple sensors may not be installed, so that the single-channel blind source separation technical scheme for diagnosing the gearbox mixed fault becomes a research hotspot in recent years.
Most blind source signal separation methods are only suitable for positive or over-definite situations, so that aiming at the special underdetermined situation of a single channel, a virtual multi-channel idea is provided, the underdetermined situation is converted into the positive or over-definite situation by using a channel expansion technology, for example, patent documents with application publication number CN 109029973A and publication number 2018, 12 and 18 disclose a method for realizing mixed fault diagnosis of a single-channel gearbox, signals are mapped to a plurality of Inherent Modal Functions (IMF) by adopting Empirical Mode Decomposition (EMD), IMF components are selected by using a kurtosis criterion and a correlation coefficient criterion and then are reconstructed with source single-channel signals to form multi-channel signals, and the single channel is expanded to be multi-channel. However, EMD has poor anti-aliasing effects and it fails when the source signal is not a natural mode. In order to improve the method, a single-channel signal is decomposed by using Ensemble Empirical Mode Decomposition (EEMD), for example, patent documents with an authorization publication number of CN 107192553B and a publication date of 2018, 3, month and 2 disclose a gearbox composite fault diagnosis method based on blind source separation, white noise is added to a vibration signal after noise reduction to carry out EEMD signal decomposition, and the defect of poor saw resistance effect is overcome. However, the method for expanding the multichannel signal does not have a standard evaluation index, and still depends on human judgment, and needs to be based on prior information about the source signal. In gear fault detection, the C-C method of the phase space reconstruction technology can map a mixed signal source to a multi-dimensional space without human intervention. For example, patent documents with an authorization publication number of CN 106513879B and a publication date of 2019, 1 and 15 disclose a spark discharge state identification and detection method based on a chaos theory, and a reconstruction parameter is obtained by using a C-C method to reconstruct a phase space of a discharge state, so that expansion of a human intervention signal is avoided. Theoretically the reconstructed signal is equivalent to the source signal, but if the embedding dimensions and delay time are chosen correctly, an inappropriate choice of reconstruction parameters may lead to reconstruction errors.
In addition, signals of expansion reconstruction need to be separated, and FastICA is a fast algorithm of self-adaptive multi-channel blind source separation developed based on an independent analysis (ICA) algorithm and is widely applied to blind source separation. For example, patent document No. CN 103575523B, published as 2015, 12, 75523, and 9, discloses a rotating machine fault diagnosis method based on FastICA-spectral kurtosis-envelope spectrum analysis, and a FastICA method based on negative entropy maximization is adopted to decouple and separate measured multi-channel acceleration signals, so that the acquired mixed vibration signals are effectively separated. However, there is a wrong choice for the selection of the nonlinear function in the FastICA algorithm, which causes the FastICA algorithm to fall into a local extreme and thus deteriorates the separation performance, so that the selection of the nonlinear function is very important but difficult.
Therefore, the effect of improper selection of reconstruction parameters and the correct selection of the nonlinear function are effectively reduced, and the separation of multi-source faults in the gearbox is important.
Disclosure of Invention
In order to reduce the influence of improper signal reconstruction parameter selection on the reconstructed signal, and simultaneously avoid the influence of FastICA on the separation effect due to wrong selection of a nonlinear function. The invention provides a double-core micro-processing system for multi-fault separation of a single-channel gear box, which can analyze and process acceleration vibration signals of a single channel and effectively separate fault signals in mixed signals.
In order to achieve the technical purpose, the technical scheme of the invention is realized as follows:
a single-channel gearbox multi-fault separation dual-core micro-processing system is composed of a DSP and an ARM core, and comprises a phase space reconstruction parameter estimation module, a phase space reconstruction module, a reconstruction parameter adjustment module, a reconstruction signal separation module, a Fourier transform module and a time-frequency domain spectrogram drawing module.
The DSP and the ARM of the dual-core micro-processing system are clear in division of labor for data analysis and processing, and the division of labor is as follows:
the DSP core receives external input signal and determines the number of processesAccording to the object, let x ═ (1,2, … N) denote a time series of length N; firstly, a phase space reconstruction parameter estimation module solves the estimation value of the phase space reconstruction parameter, the estimation algorithm adopts a representative algorithm C-C method that the embedded dimension and the time delay are mutually dependent, and only a time window tau is ensuredwThe method comprises the steps of obtaining a single-channel signal, obtaining a multi-channel signal by using a finite support sample kernel function, obtaining a multi-channel signal by using a finite FasICA fusion algorithm, obtaining a multi-source time-frequency domain signal, obtaining a multi-channel signal by using a finite support sample kernel function, obtaining a multi-source time-frequency domain signal by using a finite support sample kernel function, obtaining a multi-channel signal by using a multi-channel signal, obtaining a multi-source time-frequency domain signal by using a multi-channel signal by using a time delay method, obtaining a group of matrix X (N- (m-1) tau-1) × m by using a phase space reconstruction module according to estimated phase space reconstruction parameters m and tau, adjusting the reconstruction parameter by using a principal component analysis method, reducing the influence of noise and redundancy in the expanded multi-channel signal by finding the principal component in X, and then mapping the signal to a new multi-.
And the ARM core draws the time frequency signal of the separated signal finally processed by the DSP and then is in butt joint with peripheral communication equipment.
For errors existing in the selection of parameters of the phase space reconstruction method, no phase space reconstruction parameter estimation algorithm can correctly estimate reconstruction parameters so far, so that the phase space reconstructed by the reconstruction parameters obtained according to the existing phase space reconstruction parameter estimation algorithm has reconstruction errors, and the most important reconstruction errors are that the reconstruction signals contain noise and redundant information. The principal component analysis method can take the reconstructed multi-dimensional signals as an original signal matrix, and map the reconstructed signals to a new multi-dimensional space by searching the principal components of the signals, so that noise and redundancy are reduced.
The specific processing process of the adjusting module of the phase space reconstruction parameter in the DSP kernel of the dual-core micro-processing system comprises the following steps:
s41: solving the covariance matrix cov (X) of the reconstructed signal and the covariance matrix cov (X) of the phase space matrix as shown in the following formula
Figure GDA0002538765390000031
In Cov (X), Cov (x) in the diagonalii) Is the variance of each dimension itself. Cov (x)ij) (i ≠ j) is the covariance of the i dimension and the j dimension. If τ is greater than its optimum value, noise is contained in the reconstructed phase space, i.e., Cov (x)ii) Will be smaller. If m is greater than its optimum, the reconstruction size is too large and Cov (x)ij) Will be larger. If Cov (x)ii) Large and Cov (x)ij) (i ≠ j) is less than the set value, the process is finished, if not, the process is switched to the next step
S42 calculating the eigenvalue lambda ∈ RmFeature vector P ∈ Rm×m
S43, removing the eigenvalue with the eigenvalue only accounting for 15 percent of the sum of the eigenvalues, namely removing the noise and redundant components, leaving the fault components of the gearbox, and finally changing the new eigenvector into P' ∈ Rm×e(e<m)
S44, retrieving a new reconstructed signal y (y ∈ R) according to the new eigenvalue and eigenvector((N-(m-1)τ-1)×e)) Thus, by eliminating the smallest eigenvalues and corresponding eigenvectors, the redundant dimension is reduced, enabling the characterization of a true real system to be revealed.
In order to separate a source signal from a mixed signal, an accurate probability density function G (y) needs to be known, a specific nonlinear function is used for replacing G (y) for FsatICA, however, under the condition that a source signal such as a gear box is unknown, the selection of the nonlinear function becomes the largest influence factor for separating the mixed signal, in order to solve the problem, a FastICA fusion algorithm based on a limited support sample kernel function is adopted by a reconstructed signal separation module, the probability density function of the reconstructed signal is estimated by the algorithm, the nonlinear function which accords with the statistical characteristic of the source signal is obtained, negative entropy is taken as a target function, and the advantages of the FastICA algorithm batch processing calculation method are combined, so that blind separation of the mixed signal is rapidly and accurately realized.
The specific processing process of the reconstruction signal separation module in the DSP kernel of the dual-core micro-processing system comprises the following steps:
s51: centering the reconstructed signal y to make the mean value of the reconstructed signal y become 0;
s52: whitening the signal to remove the correlation of the data;
s53: selecting the number m' of signals to be estimated, and setting iteration number p ← 1;
s54: selecting an initial weight vector (random) Wn
S55: sorting the reconstructed signals in a non-decreasing order to obtain y ═ y { (y)1,y2,…yNIs obtained by
Figure GDA0002538765390000041
Estimating a probability density function of a kernel function of a finite support sample, wherein
Figure GDA0002538765390000042
For M finite support samples yi,yi+1,…yi+MAnd M denotes the length of the window, i.e. M denotes the length of the window
Figure GDA0002538765390000043
Deriving a lowest order kernel function by satisfying boundary conditions of limited support samples and basic conditions of a probability density function
Figure GDA0002538765390000044
In the formula, mui=(yi+M-yi) (ii) obtaining G according to the conditions of the kernel functioni=15μi -5And/16, finally obtaining the derivative of the probability density p (y) of random y:
Figure GDA0002538765390000045
thereby estimating a probability density function, solving a nonlinear function g, and solving the solutions of n optimization problems based on the principle of' separating one by one
Figure GDA0002538765390000046
The constraint conditions are as follows:
Figure GDA0002538765390000047
according to Kuhn Tucker conditions and | W | 1, may be passed through f (W) ═ E { zG (W)Tz) } + β W is 0, the optimal solution is solved, β is a constant, and a Newton iterative formula is adopted to solve the solution;
s56: order to
Figure GDA0002538765390000048
Wherein g' (. cndot.) represents the first derivative of g (-);
s57: reduction of Newton iterations to a reduced number of newton iterations through a series of reduction processes
Figure GDA0002538765390000049
Figure GDA00025387653900000410
Wn+1←Wn+1/‖Wn+1‖;
S58: if W isnIf not, returning to the step S56, otherwise, continuing;
s59: let n equal to n +1, if n ≦ m', return to step S54.
The method has the technical effects that the advantages of each core in the double-core are utilized on hardware, the DSP core is utilized to carry out phase space reconstruction parameter estimation, phase space reconstruction, reconstruction parameter adjustment, reconstruction signal separation and Fourier transformation, and the ARM core is utilized to carry out signal separation time-frequency domain spectrogram drawing. Aiming at the problems existing in the traditional single-channel blind source separation, the method for adjusting the parameters of the reconstructed signal by adopting a principal component analysis method is provided on software, so that the influence of noise and redundant signals in the reconstructed signal is reduced, finally, the source signal which is recovered due to the inaccuracy of a selected nonlinear function can be unsatisfactory when a traditional wider FastICA separation algorithm is used, and the FastICA fusion algorithm based on a limited support sample kernel function is adopted, so that the problem is solved, and finally, a high-speed and effective processing system can be provided for a diagnosis object with a complex environment and running condition, namely a gearbox.
Drawings
FIG. 1 is a block diagram of the system modules of the present invention;
FIG. 2 is a flow chart of a method of the present invention;
FIG. 3 is a flow chart of a reconstruction parameter adjustment module according to the present invention;
fig. 4 is a flow chart of the reconstructed signal separation module of the present invention.
Detailed Description
For a more clear description of the invention, reference is now made to the following detailed description and accompanying drawings.
Referring to fig. 1, the dual-core microprocessor system is composed of a DSP and an ARM core, and includes a phase space reconstruction parameter estimation module, a phase space reconstruction module, a reconstruction parameter adjustment module, a reconstruction signal separation module, and a time-frequency domain spectrogram drawing module. The phase space reconstruction parameter estimation module, the phase space reconstruction module, the reconstruction parameter adjustment module, the reconstruction signal separation module and the Fourier transform module belong to a DSP kernel operation module, and the time-frequency domain spectrogram drawing module belongs to an ARM module.
Referring to fig. 2, in the signal processing flow chart of the processing system, the DSP core receives an external input signal, establishes a processing data object, and sets x to (1,2, … N) to indicate that a time series with a length of N is a processing object. Firstly, a phase space reconstruction parameter estimation module in a DSP kernel obtains an estimation value of a phase space reconstruction parameter according to an acceleration vibration signal of a single-channel gear box, the estimation algorithm adopts a representative algorithm C-C method that an embedded dimension and time delay are interdependent, and only a time window tau is ensuredw(m-1) × τ is constant, where m is the embedding dimension and τ is the lag time; the phase space reconstruction module reconstructs the parameter m andtau, performing multichannel expansion on the single-channel signal by adopting a time delay method, and reconstructing a phase space to obtain a group of matrixes X (N- (m-1) × tau-1) × m;
referring to fig. 3, a flow chart of the reconstruction parameter adjusting module, the reconstruction parameter adjusting module adjusts the reconstruction parameters of the signals reconstructed according to the preliminarily estimated reconstruction parameters by using a principal component analysis method, and solves a covariance matrix cov (x) of the reconstructed signals, and a covariance matrix cov (x) of the phase space matrix is as shown in the following formula:
Figure GDA0002538765390000051
since if τ is greater than its optimal value, noise is contained in the reconstructed phase space, i.e., Cov (x)ii) Will be smaller. If m is greater than its optimum, the reconstruction size is too large and Cov (x)ij) Will be larger, so look at Cov (x)ii) And Cov (x)ij) (i ≠ j) size, Cov (x)ii) Large and Cov (x)ij) (i ≠ j) is less than the set value, the process is ended, if not, the process is switched to the next step, and then the characteristic value lambda ∈ R of the matrix is calculatedmFeature vector P ∈ Rm×mRemoving the characteristic value with the characteristic value only accounting for 15% of the total sum of the characteristic values, namely removing noise and redundant components, leaving the characteristic value of the suspected fault component of the gear box, and finally changing the new characteristic vector into P' ∈ Rm×e(e<m) recovering a new reconstructed signal y (y ∈ R) from the new eigenvalues and eigenvectors((N -(m-1)τ-1)×e)) Therefore, the redundant dimension is reduced by eliminating the minimum characteristic value and the corresponding characteristic vector, and the reconstruction parameter adjusting module completes the adjustment of the reconstruction parameter of the reconstruction phase space signal.
The above process completes the process of expanding from single channel to multiple channels and converting from under-determined condition to multi-source mixed signal of positive or over-determined condition, and then the operation is to carry out blind source separation on a multi-channel signal.
Referring to fig. 4, a detailed flowchart of the reconstructed signal separation module, the reconstructed signal separation module applies a FastICA fusion algorithm based on a finite support sample kernel to the adjusted new multi-dimensional phase space signal.
S51: centering the reconstructed signal y to make the mean value of the reconstructed signal y become 0;
s52: whitening the signal to remove the correlation of the data;
s53: selecting the number m' of signals to be estimated, and setting iteration number p ← 1;
s54: selecting an initial weight vector (random) WnProviding service for subsequent Newton iteration;
s55: sorting the pre-processing signals y in a non-decreasing order to obtain y ═ y1,y2,…yNIs obtained by
Figure GDA0002538765390000061
Estimating a probability density function of a kernel function of a finite support sample, wherein
Figure GDA0002538765390000062
For M finite support samples yi,yi+1,…yi+MAnd M denotes the length of the window, i.e. M denotes the length of the window
Figure GDA0002538765390000063
Deriving a lowest order kernel function by satisfying boundary conditions of limited support samples and basic conditions of a probability density function
Figure GDA0002538765390000064
In the formula, mui=(yi+M-yi) (ii) obtaining G according to the conditions of the kernel functioni=15μi -5And/16, finally obtaining the derivative of the probability density p (y) of random y:
Figure GDA0002538765390000065
thereby estimating the probability densityDegree function, calculating nonlinear function g, and solving the solutions of n optimization problems based on the principle of' separating one by one
Figure GDA0002538765390000066
The constraint condition is
Figure GDA0002538765390000067
According to the Kuhn-Tucker condition and | W | 1, can be passed through f (W) ═ E { zG (W)Tz) } + β W is 0, the optimal solution is solved, β is a constant, and a Newton iterative formula is adopted to solve the solution;
s56: order to
Figure GDA0002538765390000068
Wherein g' (. cndot.) represents the first derivative of g (-);
s57: reduction of Newton iterations to a reduced number of newton iterations through a series of reduction processes
Figure GDA0002538765390000069
Figure GDA00025387653900000610
Wn+1←Wn+1/‖Wn+1‖;
S58: judgment of WnConvergence of if WnIf the divergence is not the case, returning to the step S56, otherwise, continuing the next step;
s59: let n equal to n +1, if n ≦ m', return to step S54.
Finally, the fault source signal is separated, and N is { N ═ N1,N2…NLAnd L is the estimated number of fault sources }. And the Fourier transform module performs Fourier transform on the separated time domain signal to further obtain a corresponding frequency domain signal.
And the ARM core draws the time frequency signal of the separated signal finally processed by the DSP and then is in butt joint with peripheral communication equipment.

Claims (3)

1. A single-channel gearbox multi-fault separation dual-core micro-processing system is characterized in that the micro-processing system consists of a DSP and an ARM dual-core, and the dual-core micro-processing system comprises a phase space reconstruction parameter estimation module, a phase space reconstruction module, a reconstruction parameter adjustment module, a reconstruction signal separation module, a Fourier transform module and a time-frequency domain spectrogram drawing module;
a DSP core in the single-channel gearbox multi-fault separation dual-core microprocessing system completes signal input, phase space reconstruction parameter estimation, phase space reconstruction, reconstruction parameter adjustment, reconstruction signal separation and Fourier transformation, and an ARM core completes drawing of a time-frequency domain spectrogram and communication work of a peripheral interface;
the phase space reconstruction module reconstructs signals by a time delay method, a single-channel signal X ═ (1,2, … N) is set as a time sequence with the length of N, a group of matrixes X ═ (N- (m-1) × tau-1) × m are obtained after the phase space is reconstructed, wherein m is an embedding dimension, tau is lag time, and a C-C method is adopted for the reconstruction parameter estimation module, and a time window tau needs to be ensuredw(m-1) × τ is constant;
the reconstruction parameter adjusting module is mainly used for adjusting a reconstruction signal X by a principal component analysis method, reducing the influence of noise and redundancy in the expanded multi-channel signal by removing a characteristic value and a corresponding characteristic vector of which the characteristic value accounts for 15% of the total sum of the characteristic values in the X, and mapping the signal to a new multi-dimensional space y ∈ R according to the new characteristic value and the new characteristic vector((N -(m-1)τ-1)×e)Wherein e is<m;
The reconstruction signal separation module adopts a FastICA fusion algorithm based on the limited support sample kernel function, and estimates the probability density function of the limited support sample kernel function of the reconstruction signal from the reconstruction signal
Figure FDA0002554174790000011
Figure FDA0002554174790000012
Wherein phi (-) is the lowest order kernel function, the lowest order kernel function phi (y) satisfying the boundary condition of the finite support sample and the basic condition of the probability density function
Figure FDA0002554174790000013
In the formula, mui=(yi+M-yi) (ii) obtaining G according to the conditions of the kernel functioni=15μi -5Further, a derivative H (y) of the probability density of the reconstructed signal
Figure FDA0002554174790000014
Calculating a source signal probability density function to finally obtain a nonlinear function based on the statistical characteristics of the source signals;
the single-channel gearbox multi-fault separation dual-core micro-processing system comprises the following processing processes:
s1: vibration signals of the single-channel multi-fault-source gearbox are input into the single-channel gearbox multi-fault separation dual-core micro-processing system;
s2: estimating reconstruction parameters m and tau of the single-channel multi-fault gearbox vibration signal to be analyzed by a phase space reconstruction parameter estimation module based on a C-C algorithm;
s3: expanding the single-channel signal of the gear box to be analyzed into a multi-channel signal X (N- (m-1) × τ -1) × m by a phase space reconstruction module based on a time delay method according to the reconstruction parameters obtained in the step S2;
s4: a reconstruction parameter adjusting module based on the principal component analysis method maps the reconstruction signal to a new multidimensional space to adjust the reconstruction parameters and output the multidimensional signal again by searching the principal component of the reconstruction signal X obtained in the step S3;
s5: a reconstructed signal separation module based on a finite support sample kernel function algorithm and a FastICA fusion algorithm estimates a probability density function to obtain a nonlinear function based on statistical characteristics of a reconstructed signal, and the signal separation is realized on a newly obtained reconstructed signal S4 by taking negative entropy as a target function and combining with FastICA algorithm batch processing calculation;
s6: the Fourier transform module performs Fourier transform on the separated time domain signal so as to obtain a corresponding frequency spectrum;
and S7, the time-frequency spectrogram drawing module draws a time-frequency spectrogram of the time-frequency signal obtained in the S6, and further positions and quantifies the fault of the gearbox according to the time-frequency spectrogram.
2. The single-channel gearbox multi-fault separation dual-core micro-processing system according to claim 1, wherein the step S4 comprises:
s41: solving the covariance matrix Cov (X) of the reconstructed signal if Cov (x)ii) Large and Cov (x)ij) If (i is not equal to j) is less than the set value, ending, otherwise, switching to the next step;
s42 calculating the eigenvalue lambda ∈ RmFeature vector P ∈ Rm×m
S43, removing the eigenvalue with the eigenvalue accounting for 15% of the sum of the eigenvalues, and changing the new eigenvalue into P' ∈ Rm×eWherein e is<m;
S44, recovering new reconstructed signal y, y ∈ R according to the new eigenvalue and eigenvector((N-(m-1)τ-1)×e)
3. The single-channel gearbox multi-fault separation dual-core micro-processing system according to claim 1, wherein the step S5 comprises:
s51: centering the reconstructed signal y to make the mean value of the reconstructed signal y become 0;
s52: whitening the signal to remove the correlation of the data;
s53: selecting the number m' of signals to be estimated, and setting iteration number p ← 1;
s54: randomly selecting an initial weight vector Wn
S55: estimating a probability density function by using a finite support sample kernel function algorithm, and solving a nonlinear function g;
s56: order to
Figure FDA0002554174790000021
Whereing' (. cndot.) represents the first derivative of g (-);
s57: according to
Figure FDA0002554174790000022
Wn+1←Wn+1/‖Wn+1‖;
S58: if W isnIf not, returning to S56, otherwise, continuing;
s59: let n equal to n +1, if n ≦ m', return to step S54.
CN201910232499.1A 2019-03-26 2019-03-26 Single-channel gear box multi-fault separation dual-core micro-processing system Active CN109916625B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910232499.1A CN109916625B (en) 2019-03-26 2019-03-26 Single-channel gear box multi-fault separation dual-core micro-processing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910232499.1A CN109916625B (en) 2019-03-26 2019-03-26 Single-channel gear box multi-fault separation dual-core micro-processing system

Publications (2)

Publication Number Publication Date
CN109916625A CN109916625A (en) 2019-06-21
CN109916625B true CN109916625B (en) 2020-08-18

Family

ID=66966791

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910232499.1A Active CN109916625B (en) 2019-03-26 2019-03-26 Single-channel gear box multi-fault separation dual-core micro-processing system

Country Status (1)

Country Link
CN (1) CN109916625B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110378065B (en) * 2019-07-29 2022-08-30 河海大学 Large-deformation plate nonlinear vibration signal phase space reconstruction method based on constant parameters
CN110608884B (en) * 2019-08-08 2021-01-05 桂林电子科技大学 Rolling bearing state diagnosis method based on self-attention neural network
CN111999607B (en) * 2020-07-20 2021-08-20 中国南方电网有限责任公司超高压输电公司广州局 Method and device for separating partial discharge narrow-band interference blind source under single-channel signal

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010061133A (en) * 1999-12-28 2001-07-07 이계안 A power transmit error estimating apparatus between final reduction gear and differential ring gear for ff cars
CN103868694A (en) * 2014-03-26 2014-06-18 东南大学 Embedded variable-rotation-speed bearing fault diagnosis device
CN104122086A (en) * 2014-08-11 2014-10-29 长沙理工大学 Kriging model based wind turbine gearbox fault diagnosis method
CN108956131A (en) * 2018-07-30 2018-12-07 西安英特迈思信息科技有限公司 A kind of express locomotive EEF bogie gear-box intellectual monitoring unit
CN109029973A (en) * 2018-06-22 2018-12-18 石家庄铁道大学 The method for realizing the diagnosis of single channel gear-box mixed fault

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010061133A (en) * 1999-12-28 2001-07-07 이계안 A power transmit error estimating apparatus between final reduction gear and differential ring gear for ff cars
CN103868694A (en) * 2014-03-26 2014-06-18 东南大学 Embedded variable-rotation-speed bearing fault diagnosis device
CN104122086A (en) * 2014-08-11 2014-10-29 长沙理工大学 Kriging model based wind turbine gearbox fault diagnosis method
CN109029973A (en) * 2018-06-22 2018-12-18 石家庄铁道大学 The method for realizing the diagnosis of single channel gear-box mixed fault
CN108956131A (en) * 2018-07-30 2018-12-07 西安英特迈思信息科技有限公司 A kind of express locomotive EEF bogie gear-box intellectual monitoring unit

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于小波包和独立分量分析的微弱多源故障声发射信号分离;王向红 等;《上海交通大学学报》;20160531;第50卷(第5期);第757-763页 *

Also Published As

Publication number Publication date
CN109916625A (en) 2019-06-21

Similar Documents

Publication Publication Date Title
CN109916625B (en) Single-channel gear box multi-fault separation dual-core micro-processing system
CN105068048B (en) Distributed microphone array sound localization method based on spatial sparsity
CN111274525B (en) Tensor data recovery method based on multi-linear augmented Lagrange multiplier method
CN117034043B (en) Intelligent building comprehensive energy consumption monitoring method and system based on multi-energy Internet of things
CN115013298B (en) Real-time performance online monitoring system and monitoring method of sewage pump
CN115235612B (en) Intelligent fault diagnosis system and method for servo motor
CN112084847B (en) Hyperspectral image denoising method based on multichannel truncated nuclear norm and total variation regularization
CN112906158A (en) Mechanical fault diagnosis method based on multi-sensor multivariate data fusion
CN111458746A (en) Tunnel microseismic waveform arrival time picking method based on U-Net neural network
CN111010356A (en) Underwater acoustic communication signal modulation mode identification method based on support vector machine
CN113935249A (en) Upper-layer ocean thermal structure inversion method based on compression and excitation network
CN115311531A (en) Ground penetrating radar underground cavity target automatic detection method of RefineDet network model
CN111693937A (en) Near-field signal source positioning method based on sparse reconstruction and without gridding
CN117009916A (en) Actuator fault diagnosis method based on multi-sensor information fusion and transfer learning
CN114722923A (en) Light electromechanical equipment fault diagnosis method
CN116908684A (en) Motor fault prediction method and device, electronic equipment and storage medium
CN117113243A (en) Photovoltaic equipment abnormality detection method
CN111428772A (en) Photovoltaic system depth anomaly detection method based on k-nearest neighbor adaptive voting
CN116363536A (en) Unmanned aerial vehicle inspection data-based power grid infrastructure equipment defect archiving method
CN113869289B (en) Multi-channel ship radiation noise feature extraction method based on entropy
CN114866110A (en) Frequency hopping signal parameter estimation method based on combination of elastic network model and generalized ADMM
CN109856673B (en) High-resolution Radon transformation data separation technology based on dominant frequency iterative weighting
CN113311364A (en) Permanent magnet synchronous motor inverter open-circuit fault diagnosis method based on multi-core SVM
Zhang et al. Research on Electronic Circuit Fault Diagnosis Method Based on SWT and DCNN-ELM
CN117711423B (en) Mixed underwater sound signal separation method and system combining auditory scene analysis and deep learning

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