CN100485342C - Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault - Google Patents

Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault Download PDF

Info

Publication number
CN100485342C
CN100485342C CNB2006100418821A CN200610041882A CN100485342C CN 100485342 C CN100485342 C CN 100485342C CN B2006100418821 A CNB2006100418821 A CN B2006100418821A CN 200610041882 A CN200610041882 A CN 200610041882A CN 100485342 C CN100485342 C CN 100485342C
Authority
CN
China
Prior art keywords
feature
vector machine
integrated
intelligent
supporting vector
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.)
Expired - Fee Related
Application number
CNB2006100418821A
Other languages
Chinese (zh)
Other versions
CN1811367A (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.)
Shenji Group Kunming Machine Tool Co., Ltd.
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CNB2006100418821A priority Critical patent/CN100485342C/en
Publication of CN1811367A publication Critical patent/CN1811367A/en
Application granted granted Critical
Publication of CN100485342C publication Critical patent/CN100485342C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The present invention discloses an integrated support vector machine mixed intelligent diagnosis method of machine failure. Said method includes the following steps: respectively adopting and lifting small wave packet, utilizing frequency band and empirical mode decomposition process to decompose vibration signal according to the eigen mode component and extract time domain statistical character of decomposed signal to form total characteristic set; providing characteristic distance evaluation technique and characteristic evaluation index; utilizing said characteristic evaluation index to select most sensitive characteristic set from the total characteristic set; using said most sensitive characteristic set as diagnosis characteristics and creating integrated support vector machine mixed intelligent diagnosis model so as to implement intelligent diagnosis of machine failure state.

Description

The integrated supporting vector machine mixed intelligent diagnosing method of mechanical fault
Technical field
The invention belongs to mechanical equipment fault intelligent diagnostics field, be specifically related to a kind of integrated supporting vector machine mixed intelligent diagnosing method of mechanical fault.
Background technology
In order to break away from the problem that mechanical fault diagnosis is too dependent on professional and technical personnel and diagnostician, efficient, the reliable intelligent diagnostics of canbe used on line, in recent years, people are applied to the intelligent trouble diagnosis of plant equipment with artificial intelligence technologys such as fuzzy theory, expert system, neural network and cluster analyses, have obtained certain effect in practice.But it is found that in actual applications these technology no all roses, for example: the fuzzy fault diagnosis often need manually be determined subordinate function and fuzzy relation matrix by priori, but in fact, there are many difficulties in subordinate function and the fuzzy relation matrix that obtains to conform to the equipment actual conditions; Also be faced with many difficulties in the expert system, as lack the expression of effective fault diagnosis knowledge, experimental knowledge is obtained difficulty etc., and in addition, the operating personnel of expert system also need to possess higher level, is difficult for popularizing; The master sample of some need be provided during neural network classification, the master sample that obtains all faults of equipment is quite to be difficult for, in addition, neural network all needs raw data is carried out feature extraction in actual applications, if do not do feature extraction, and directly participate in calculating with raw data, the structure of network will be very huge, can't satisfy the requirement of real time on-line monitoring diagnosis; Cluster analysis can not need the failure criterion sample, it forms a plurality of classes with data sample, data sample difference between the requirement class should be big as much as possible, and the difference between the data sample should be as much as possible little in the same class, promptly satisfy the principle of " similarity between minimized class; similarity in the maximized class ", therefore, undesirable to the cluster effect of early stage Weak fault or combined failure.
In engineering practice, fault diagnosis for large complicated plant equipment, often measurand is many, the data volume of analyzing and processing is big, carry out the fault sample famine of intelligent diagnostics model training, and plant equipment is the system of a complexity, time variation, randomness, many-sided factor such as ambiguity makes fuzzy theory, expert system, traditional intelligent diagnosing method such as neural network and cluster analysis is difficult to the malfunction of equipment operation is made identification accurately and effectively, particularly for potential, early stage Weak fault lacks effective recognition and diagnostic means.In order to solve the existing problem in intelligent trouble diagnosis field, urgently need research and introduce new theory and technology, new, intelligent diagnosis technology and method are efficiently proposed.
Support vector machine is a kind of novel intelligence learning machine that development in recent years is got up, and is suitable for the mechanical disorder pattern-recognition of few sample, and has obtained certain achievement in fault diagnosis.Because the high complexity of space-time, the execution of support vector machine usually is a kind of approximate treatment yet in actual applications; Simultaneously, the selection problem of kernel function parameter makes support vector machine that study very easily take place or owes to learn phenomenon, directly influences it and promotes performance.
Integrated supporting vector machine can improve the classification performance of single support vector machine effectively, yet, when very faint or fault sample is not true to type when the feature of initial failure, still may produce wrong diagnostic result.So need improve the quality of diagnostic message by the modern signal processing method in the link of feature extraction, thereby improve the accuracy and the reliability of fault diagnosis.Select suitable modern signal analytical approach, the fault that machine essential fault, initial failure etc. are difficult to indicate all may be discerned, for the intelligent diagnostics of early stage incipient fault provides foundation.Integrated use artificial intelligence and systematic analytic method advantage is separately learnt from other's strong points to offset one's weaknesses, and has complementary advantages, can improve the rapidity and the accuracy of diagnosis effectively, reduce misdiagnosis rate and rate of missed diagnosis, determine correct time and position that fault takes place, estimate its order of severity and development trend.
Summary of the invention
The objective of the invention is to overcome above-mentioned the deficiencies in the prior art, a kind of integrated supporting vector machine mixed intelligent diagnosing method of mechanical fault is provided, this method precision height, cost are low, simple and reliable, be convenient to use in the engineering practice, by Lifting Wavelet packet transform and these two kinds of modern signal processing technology of empirical mode decomposition, characteristic distance assessment technology and integrated supporting vector machine are effectively merged, realize intelligent diagnostics to the failure operation state.
Technical scheme of the present invention is to solve like this: undertaken by following step:
1) adopts the Lifting Wavelet bag by the eigenmodes component vibration signal to be decomposed respectively, extract the time domain statistical nature of decomposed signal, constitute all feature sets by frequency range and empirical mode decomposition;
2) propose characteristic distance assessment technology and feature evaluation index, from all feature sets, choose sensitive features collection according to feature evaluation index size;
The sensitive features collection that 3) will obtain is set up the integrated supporting vector machine mixed intelligent diagnostic model as diagnostic characteristic, realizes the intelligent diagnostics to equipment running status.
Said characteristic distance assessment technology and feature evaluation index are as follows:
At first, for an original vibration signal X (t), calculate 1 time domain statistical nature F Time, 1 time domain statistical nature is meant one or more in average, variance, root amplitude, effective value, peak value, measure of skewness, kurtosis, waveform index, peak value index, nargin index and the pulse index; Utilize the Lifting Wavelet packet transform that X (t) is decomposed the m layer, the m value is 2,3 or 4, obtains 2 mIndividual decomposed signal calculates 1 statistical nature respectively to each decomposed signal, obtains 1 * 2 altogether mFeature F behind the individual Lifting Wavelet packet transform WpAdopt empirical mode decomposition that X (t) is decomposed and obtain preceding n eigenmodes component, the n value is 4 to 8 integer, calculates the relative energy feature F of n eigenmodes component EmdThereby, constitute all feature set F Total=[F TimeF WpF Emd];
Then, with all feature set F TotalThe class interval S bWith in the class apart from S wRatio be set at apart from evaluation index J A
At last, according to feature evaluation index J ASize, from all feature set F TotalThe middle feature evaluation index J that selects greater than given threshold value ρ APairing feature is made as ρ with the pairing threshold value of maximum classification accuracy, thereby constitutes sensitive features collection F Sensitivity
Described integrated supporting vector machine mixed intelligent diagnostic model method for building up is as follows:
At first utilize packed algorithm, be called for short the Bagging algorithm, from training sample, generate T training sample subclass TR at random t(t=1,2 ..., T), the T value is 10 to 100 integer, utilizes T multi-category support vector machines to these subclass TR then tTrain, obtain T sub-classifier f tIf T sub-classifier f tAverage error in classification be
Figure C200610041882D0007110820QIETU
, will f = 1 / E ^ Be the fitness function of genetic algorithm, to T sub-classifier f tIntegrated result be optimized, obtain the weight vector w of the optimization weights of sorter, thereby constitute integrated supporting vector machine, with sensitive features collection F greater than predetermined threshold value λ=1/T SensitivityAs diagnostic characteristic, integrated supporting vector machine is trained, can set up the hybrid intelligent diagnostic model, utilize the hybrid intelligent diagnostic model to realize the intelligent diagnostics of mechanical fault.
Because the present invention has realized the feature extraction based on the modern signal analytical technology, the selection and the mixing of integrated intelligent sorting technique on algorithm of sensitive features collection, the present invention has the following significant advantage that is different from classic method:
1) carries out characteristic distance assessment and can remove the uncorrelated or redundant information that primitive character is concentrated effectively, choose sensitive features collection, thereby improve the classification accuracy and the operation efficiency of sorter;
2) set up the integrated supporting vector machine mixed intelligent model and carry out intelligent classification, can carry out few sample mode identification well, improved the classification performance and the anti-noise ability of single sorter effectively, for diagnosis accurately and effectively provides assurance;
3) whole process has realized feature extraction, feature selecting and the pattern-recognition mixing on algorithm, for the intelligent diagnostics of mechanical fault provides effective practical new technology.
Description of drawings
Fig. 1 is the integrated supporting vector machine mixed intelligent diagnostic flow chart of mechanical fault;
Fig. 2 is that certain electric locomotive traveling portion wheel is to structural representation;
Fig. 3 be sensitive features collection select apart from evaluation graph;
Fig. 4 is under the different feature evaluation metrics-thresholds, the classification performance comparison diagram of integrated supporting vector machine and single support vector machine and Bagging algorithm;
Fig. 5 is under the different noise content, the classification performance comparison diagram of integrated supporting vector machine and single support vector machine and Bagging algorithm.
Embodiment
Accompanying drawing is specific embodiments of the invention;
Below in conjunction with accompanying drawing content of the present invention is described in further detail:
1) the vibratory output signal to equipment carries out Lifting Wavelet packet transform and empirical mode decomposition respectively, extract the time domain statistical nature of each different frequency range component, constitute all feature sets, all feature sets that obtain are carried out the characteristic distance assessment, choose sensitive features collection;
The sensitive features collection that 2) will obtain is trained integrated supporting vector machine as diagnostic characteristic, and the hybrid intelligent diagnostic model that trains can carry out intelligent classification to equipment running status.
According to the hybrid intelligent diagnostic structure flow process of foregoing invention content and Fig. 1, at first, structure primitive character collection.
For an original vibration signal X (t), extract 11 time domain statistical nature F Time(average, variance, root amplitude, effective value, peak value, measure of skewness, kurtosis, waveform index, peak value index, nargin index and pulse index).
For original vibration signal X (t), data length is L, and the j yardstick is with approximation signal { s J+1(k) } subdivision is strange, even sample sequence { s j(2k+1) } and { s j(2k) }, adopt the interpolation subdividing principle to predict strange sample, then detail signal d with N even sample j(k) be
d j ( k ) = s j + 1 ( 2 k + 1 ) - Σ m = 1 N p ( m ) s j + 1 ( 2 m + k - N ) - - - ( 1 )
P (m) is a predictive coefficient in the formula (1), order: P=[p (1) ..., p (N)] T, employing formula (1) is tried to achieve
Figure C200610041882D0008110645QIETU
Individual detail signal d j(k) upgrade even sample { s j(2k) }, approximation signal s then j(k) be
s j ( k ) = s j + 1 ( 2 k ) - Σ m = 1 N ~ u ( m ) d j ( m + k - N ~ / 2 - 1 ) - - - ( 2 )
U (m) is a update coefficients in the formula (2), order: U = [ u ( 1 ) , · · · , u ( N ) ~ ] T 。In formula (1) and formula (2), choose different N and
Figure C200610041882D00084
The scaling function and the wavelet function of different vanishing moments will be obtained having.
With N and
Figure C200610041882D00085
The small echo with impact signal characteristic of structure is a basis function, adopts to promote strategy to wavelets Subspace W jDecompose, thereby obtain subspace behind the wavelet package transforms
Figure C200610041882D00086
Signal be X j={ x J, n, l, j, n, 1 ∈ Z}, x J, n, lThe 1st data for the n frequency band of j yardstick.
To original signal X (t), utilize the Lifting Wavelet packet transform to decompose 3 layers, obtain 8 sections wavelet packet frequency band coefficients, to every section coefficient each extract with the original signal processing in 11 identical statistical natures, obtain 88 feature F behind the Lifting Wavelet packet transform altogether Wp
Empirical mode decomposition can be decomposed into arbitrary signal several eigenmodes components and a remainder.So-called eigenmodes component is exactly function or the signal that satisfies 2 conditions: 1. in whole data sequence, the quantity of extreme point must equate with the quantity of crossing at 0 or differ one at most.2. in any point, the average of the lower envelope line that coenvelope line that the local maximum point of data sequence is determined and local minizing point determine is 0, and promptly signal is about the local symmetry of time shaft.The essence of empirical mode decomposition is the process of a screening, and through after a series of decomposition, time series X (t) can be expressed as n eigenmodes component f i(t) and a remainder r n(t) sum, promptly
x ( t ) = Σ i = 1 n f i ( t ) + r n ( t ) - - - ( 3 )
N the eigenmodes component f that obtains in the formula (3) i(t), its frequency is arranged from big to small, f 1(t) contained frequency is the highest, f n(t) contained frequency is minimum, remainder r n(t) be the monotonic sequence of a non-concussion.
To original signal X (t), utilize empirical mode decomposition to obtain the first six eigenmodes component, extract the relative energy feature F of 6 eigenmodes components Emd, obtain all feature set F Total=[F TimeF WpF Emd].
Then, to primitive character collection F TotalCarry out feature selecting.
Suppose c mode class ω 1, ω 2..., ω cThe associating set of eigenvectors be { q (i, k), i=1,2 ..., c; K=1,2 ..., N i, q wherein (i, k)Be k feature among the ω i, N iBe ω iThe number of middle proper vector.Feature selecting can be divided into three steps:
The first step: calculate ω iMean distance in the class between all proper vectors is as follows
S i = 1 2 N i Σ j = 1 N j 1 N i - 1 Σ k = 1 N i | q ( i , j ) - q ( i , k ) | - - - ( 4 )
To S i(i=1,2 ..., c) ask and obtain that distance is in the average class after average
S w = 1 c Σ i = 1 c S i - - - ( 5 )
Second step: the between class distance that calculates c mode class is as follows
S b = 1 c Σ i = 1 c | μ ( i ) - μ | - - - ( 6 )
Wherein: μ ( i ) = 1 N i Σ k = 1 N i q ( i , k ) Be ω iIn the average of all features, μ = 1 c Σ i = 1 c 1 N i Σ k = 1 N i q ( i , k ) Population mean for c mode class sample.
The 3rd step: the ratio J of distance in definition class spacing and the class ABe the distance evaluation index
J A = S b S w - - - ( 7 )
As can be seen, distance and big average between class distance just have good separability in the little average class from the definition of (7) formula, therefore select the J greater than certain threshold value APairing feature constitutes sensitive features collection F Sensitivity, the pairing threshold value of maximum classification accuracy is made as ρ.
At last, set up the integrated supporting vector machine mixed intelligent diagnostic model, to sensitive features collection F SensitivityCarry out intelligent classification.
(Support Vector Machines SVMs) is based on a kind of small sample algorithm for pattern recognition of Statistical Learning Theory and structural risk minimization to support vector machine.In the integrated supporting vector machine algorithm, at first utilize the Bagging algorithm from training sample, to generate T training sample subclass TR at random t(t=1,2 ..., T), utilize T multi-category support vector machines then to these subclass TR tTrain, obtain T sub-classifier f tIn the present invention, establish T sub-classifier f tAverage error in classification be
Figure C200610041882D0007110820QIETU
, will f = 1 / E ^ As the fitness function of genetic algorithm, to T sub-classifier f tIntegrated result be optimized, obtain the weight vector w of the optimization weights of sorter, thereby constitute integrated supporting vector machine greater than predetermined threshold value λ=1/T.Its algorithm is as shown in table 1.
Table 1 integrated supporting vector machine algorithm
Figure C200610041882D00106
Figure C200610041882D00111
The supporting vector machine mixed intelligent diagnostic model that utilization trains can be diagnosed equipment running status, thereby can obtain the type of fault.
With reference to shown in Figure 1, be the integrated supporting vector machine mixed intelligent diagnostic flow chart of mechanical fault, the vibration signal of plant equipment picks up through sensor and data acquisition system (DAS), and after the pre-service of signal pre-processing system, can obtain original time domain vibration signal; 88 statistical natures after 11 temporal signatures, Lifting Wavelet bag that calculate original time domain vibration signal decompose and 6 relative energy features after the empirical mode decomposition constitute all feature sets; Utilize the characteristic distance assessment technology that all feature sets are carried out feature selecting, constitute sensitive features collection; Sensitive features set pair integrated supporting vector machine with training sample is trained, and makes up the integrated supporting vector machine mixed intelligent diagnostic model, utilizes this model to carry out the real-time intelligent diagnosis to plant equipment, obtains the hybrid intelligent diagnostic result at last.
With reference to shown in Figure 2, for certain electric locomotive traveling portion wheel to structural representation, each is taken turns being made up of two bearings 1 and 5, two wheels 2 and 4 and axletrees 3, degree of will speed up sensor 6 is installed in and carries out vibration signal above the bearing seat and pick up.
With reference to shown in Figure 3, for sensitive features collection select apart from evaluation graph, according to preset threshold ρ, thus can determine greater than this threshold value apart from evaluation index J APairing feature is sensitive features collection, and wherein, horizontal ordinate is the feature sequence number, and ordinate is apart from evaluation index J A
With reference to shown in Figure 4, under different feature evaluation metrics-thresholds, integrated vector machine and the single classification performance comparison diagram of holding vector machine and Bagging algorithm held, by with classic method more as can be seen, under different threshold values, the present invention has good classifying quality, and the feature of selecting when ρ=20 is sensitive features collection, wherein, horizontal ordinate is threshold value ρ, and ordinate is a classification accuracy.
With reference to shown in Figure 5, under different noise content, the classification performance comparison diagram of integrated SVMs and single SVMs and Bagging algorithm, by with classic method more as can be seen, under different noise content, the present invention has very high diagnosis efficiency height and anti-noise ability, wherein, horizontal ordinate is a noise content number percent, and ordinate is a classification accuracy.
Embodiment:
This embodiment has provided the specific implementation process of the present invention in engineering practice, simultaneous verification should the invention validity.
Certain passenger-carrying version electric locomotive traveling portion by six accessory whorls to forming.Wheel is to being made up of an axletree and two wheels, and every accessory whorl is to combining with axle box again.Wheel to structure as shown in Figure 2.Carry out signals collecting by the acceleration transducer that is installed in axletree bearing top, sample frequency is 12.8KHz, and data length is 8192 points.
Get this electric locomotive bearing each 36 groups of four kinds of operating modes (normal, outer ring fault, rolling body fault, outer ring and rolling body combined failure) vibration datas down, wherein 22 groups as training data, other 14 groups as test sample book.Concentrate all feature set F that extract from the training data of electric locomotive bear vibration Total=[F TimeF WpF Emd] form by 105 features, these features apart from the evaluation index value as shown in Figure 3.
In order to verify the classification performance of integrated SVMs method, also analyze simultaneously to single SVMs sorting algorithm with based on the Bagging algorithm classification algorithm of SVMs.Get the number T=20 of integrated SVMs herein, can select the predetermined threshold value λ=1/T=0.05 among the integrated SVMs.
The classification results of single SVMs, Bagging algorithm and integrated SVMs is as shown in table 2.Wherein nicety of grading is that each test repeats 10 times average result.
SVMs, the Bagging of the different threshold value correspondences of table 2 and the nicety of grading of integrated Bagging are relatively
Figure C200610041882D00131
As shown in table 2, along with the increase of threshold value, number of features reduces gradually.In training process, the nicety of grading of single SVMs is all greater than 98.95%, and in ρ=0.5,5 or obtained maximal value (100%) at 20 o'clock; Yet the nicety of grading of the training sample of Bagging algorithm and integrated SVMs all is 100%.It can also be seen that in addition in integrated SVMs, the integrated number (9.3) of average SVMs is approximately half of overall number (20), therefore compares with the Bagging algorithm, can select the test duration of integrated SVMs to reduce significantly (dropping to 0.02 second from 0.12 second).
At different feature selecting threshold values, the test sample book classification results of single SVMs, Bagging algorithm and integrated SVMs as shown in Figure 4.
As can be seen from Figure 4, for the test result of single SVMs, when not carrying out feature selecting, its classification accuracy only is 63.33%; Classification accuracy improves along with the increase of threshold value, is up to 97.05%, this moment threshold value ρ=20, this feature that shows that ρ=20 o'clock are selected is sensitive features collection.Can find that the classification performance of integrated SVMs is better than other two kinds of methods, when ρ=20, integrated SVMs obtains best classifying quality (100%).
Further verify the extensive performance of integrated SVMs, in test data, add the random noise of different content.Be without loss of generality, the number T=20 that gets integrated SVMs in Bagging algorithm and integrated SVMs experimentizes, threshold value ρ=20.The classification accuracy of different noise content as shown in Figure 5.As can be seen from the figure, along with the increase of noise content, for integrated SVMs, when noise content less than 20% the time, its nicety of grading is all greater than 90%.When noise content was between 0 and 40%, the nicety of grading of integrated SVMs all was higher than the Bagging algorithm.Further increase along with noise content, the nicety of grading of the nicety of grading of integrated SVMs and Bagging algorithm is suitable, all be about 70%, but in integrated SVMs, the integrated number (9.3) of average SVMs is approximately half of integrated number (20) of Bagging algorithm, can select the test duration of integrated SVMs to reduce significantly (dropping to 0.02 second from 0.12 second).These illustrate that all integrated SVMs counting yield height, anti-noise ability are strong.
From the analysis of Fig. 4 and Fig. 5 as can be known, the classification performance of integrated SVMs is better than single SVMs and Bagging algorithm, fault diagnosis for the locomotive wheel set bearing, when the threshold setting in the characteristic distance assessment technology is ρ=20, its classification accuracy is 100%, this explanation integrated supporting vector machine mixed intelligent diagnosing method can be diagnosed out outer ring fault, rolling body fault and outer ring and the rolling body combined failure of bearing exactly, has good intelligent diagnostics ability.

Claims (1)

1. the integrated supporting vector machine mixed intelligent diagnosing method of a mechanical fault comprises:
1) adopts the Lifting Wavelet bag by the eigenmodes component original vibration signal to be decomposed respectively, extract the time domain statistical nature of decomposed signal and original vibration signal, constitute all feature sets by frequency range and empirical mode decomposition;
It is characterized in that:
2) propose ratio with distance in the class spacing of all feature sets and the class as the distance evaluation index, from all feature sets, choose sensitive features collection according to distance evaluation index size;
At first, for an original vibration signal X (t), calculate l time domain statistical nature F Time, l time domain statistical nature is meant one or more in average, variance, root amplitude, effective value, peak value, measure of skewness, kurtosis, waveform index, peak value index, nargin index and the pulse index; Utilize the Lifting Wavelet packet transform that X (t) is decomposed the m layer, the m value is 2,3 or 4, obtains 2 mIndividual decomposed signal calculates l statistical nature respectively to each decomposed signal, obtains l * 2 altogether mFeature F behind the individual Lifting Wavelet packet transform WpAdopt empirical mode decomposition that X (t) is decomposed and obtain preceding n eigenmodes component, the n value is 4-8, calculates the relative energy feature F of n eigenmodes component EmdThereby, constitute all feature set F Total=[F TimeF WpF Emd];
Then, with all feature set F TotalThe class interval S bWith in the class apart from S wRatio be set at apart from evaluation index J A
At last, according to distance evaluation index J ASize, from all feature set F TotalIn select greater than given threshold value ρ apart from evaluation index J APairing feature is made as ρ with the pairing threshold value of maximum classification accuracy, thereby constitutes sensitive features collection F Sensitivity
The sensitive features collection that 3) will obtain is set up the integrated supporting vector machine mixed intelligent diagnostic model based on packed algorithm and genetic algorithm as diagnostic characteristic, realizes the intelligent diagnostics to equipment running status;
At first utilize packed algorithm from training sample, to generate T training sample subclass TR at random t(t=1,2 ..., T), the T value is 10-100, utilizes T multi-category support vector machines to these subclass TR then tTrain, obtain T sub-classifier f tIf T sub-classifier f tAverage error in classification be
Figure C200610041882C00031
Will f = 1 / E ^ As the fitness function of genetic algorithm, to T sub-classifier f tIntegrated result be optimized, obtain the weight vector w of the optimization weights of T sub-classifier, thereby constitute integrated supporting vector machine, with sensitive features collection F greater than predetermined threshold value λ=1/T SensitivityAs diagnostic characteristic, integrated supporting vector machine is trained, can set up the hybrid intelligent diagnostic model, utilize the hybrid intelligent diagnostic model to realize the intelligent diagnostics of mechanical fault.
CNB2006100418821A 2006-03-03 2006-03-03 Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault Expired - Fee Related CN100485342C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100418821A CN100485342C (en) 2006-03-03 2006-03-03 Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100418821A CN100485342C (en) 2006-03-03 2006-03-03 Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault

Publications (2)

Publication Number Publication Date
CN1811367A CN1811367A (en) 2006-08-02
CN100485342C true CN100485342C (en) 2009-05-06

Family

ID=36844417

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100418821A Expired - Fee Related CN100485342C (en) 2006-03-03 2006-03-03 Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault

Country Status (1)

Country Link
CN (1) CN100485342C (en)

Families Citing this family (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101750210B (en) * 2009-12-24 2012-05-30 重庆大学 Fault diagnosis method based on OLPP feature reduction
CN101832857B (en) * 2010-03-31 2012-06-27 桂林电子科技大学 Method for diagnosing spindle bearing fault of high speed flash frequency observation device of electron beam welder
CN102179722A (en) * 2010-12-20 2011-09-14 西安瑞特快速制造工程研究有限公司 Estimation method for operational reliability of numerical control machine based on proportional hazards model
CN102539159A (en) * 2010-12-24 2012-07-04 中国船舶研究设计中心 Fault diagnosis method for valve mechanism of diesel engine
CN102175449B (en) * 2011-03-18 2012-07-18 天津工业大学 Blade fault diagnostic method based on strain energy response of wind-driven generator
CN102620928A (en) * 2012-03-02 2012-08-01 燕山大学 Wind-power gear box fault diagnosis method based on wavelet medium-soft threshold and electronic-magnetic diaphragm (EMD)
CN102736027A (en) * 2012-07-18 2012-10-17 南京因泰莱配电自动化设备有限公司 Circuit breaker failure diagnosis method based on circuit breaker dynamic property test instrument
CN102829974B (en) * 2012-08-07 2015-01-07 北京交通大学 LMD (local mean decomposition) and PCA (principal component analysis) based rolling bearing state identification method
CN103115789B (en) * 2013-01-17 2015-01-21 西安交通大学 Second generation small-wave support vector machine assessment method for damage and remaining life of metal structure
CN103163864B (en) * 2013-03-07 2014-11-19 电子科技大学 Method for optimizing mechanical equipment state estimation
CN103217280B (en) * 2013-03-18 2015-10-28 西安交通大学 The multivariable support vector machine prediction method of aero-engine rotor residual life
CN103674511B (en) * 2013-03-18 2016-06-22 北京航空航天大学 A kind of mechanical wear part Performance Evaluation based on EMD-SVD and MTS and Forecasting Methodology
CN103207952A (en) * 2013-04-23 2013-07-17 华北科技学院 Slope displacement prediction method
CN103542929A (en) * 2013-05-20 2014-01-29 西北工业大学 Method for extracting signal features of bearings on basis of wavelet packet energy matrixes
CN103335840A (en) * 2013-07-02 2013-10-02 中煤科工集团西安研究院 Intelligent diagnosis method for faults of mining drilling machine gearbox
CN103323228A (en) * 2013-07-02 2013-09-25 中煤科工集团西安研究院 Mining drill fault intelligent identification method
FR3011408A1 (en) * 2013-09-30 2015-04-03 Orange RE-SAMPLING AN AUDIO SIGNAL FOR LOW DELAY CODING / DECODING
CN103471870B (en) * 2013-10-08 2016-01-20 中国矿业大学(北京) Based on the characteristic index of two-dimensional projections to the susceptibility of coal mine equipment state recognition and cluster Journal of Sex Research method
CN103528820B (en) * 2013-10-10 2015-11-18 辽宁大学 A kind of Fault Diagnosis of Roller Bearings based on distance evaluation factor potential-energy function
CN104655423B (en) * 2013-11-19 2017-09-15 北京交通大学 A kind of Fault Diagnosis of Roller Bearings based on time-frequency domain multi-dimensional vibration Fusion Features
CN103821750B (en) * 2014-03-05 2016-04-06 北京工业大学 A kind of ventilator stall based on electric current and surge monitoring and diagnostic method
CN103868692B (en) * 2014-03-18 2016-02-10 电子科技大学 Based on the rotary machinery fault diagnosis method of Density Estimator and K-L divergence
CN105134619B (en) * 2015-09-28 2017-03-01 北京航空航天大学 A kind of fault diagnosis based on wavelet energy, manifold dimension-reducing and dynamic time warping and health evaluating method
CN105784364A (en) * 2016-03-02 2016-07-20 西安交通大学 Bearing fault diagnosis method based on total experience mode decomposition and fractal box dimensions
CN105784340B (en) * 2016-03-03 2018-03-23 电子科技大学 Valve fault diagnosis method based on Hybrid Intelligent Technology
CN106092574B (en) * 2016-05-30 2018-04-17 西安工业大学 Based on the Method for Bearing Fault Diagnosis for improving EMD decomposition and sensitive features selection
CN106198014B (en) * 2016-06-29 2018-05-29 潍坊学院 A kind of envelope Analysis Method based on empirical mode decomposition and spectrum kurtosis
CN106198000A (en) * 2016-07-11 2016-12-07 太原理工大学 A kind of rocker arm of coal mining machine gear failure diagnosing method
CN106441946B (en) * 2016-11-15 2018-10-19 重庆工商大学 Vehicle hydraulic damper fault recognition method based on vibration signal and system
CN106841844A (en) * 2016-12-14 2017-06-13 广西大学 A kind of syntonic type electric locomotive power supply system of train method for diagnosing faults
CN107378779A (en) * 2017-06-29 2017-11-24 西安交通大学 A kind of grinding trembling detection method based on axis feeding motor current signal
CN107727395B (en) * 2017-07-21 2019-12-03 中国矿业大学 A kind of Method for Bearing Fault Diagnosis based on full variation and uncompensation distance assessment
US10991381B2 (en) 2018-04-09 2021-04-27 Well Checked Systems International LLC System and method for machine learning predictive maintenance through auditory detection on natural gas compressors
CN109238728B (en) * 2018-09-12 2020-06-02 温州大学 Method and system for diagnosing faults of parts on vehicle engine
CN109323860A (en) * 2018-10-31 2019-02-12 广东石油化工学院 A kind of rotating machinery gearbox fault data set optimization method
CN110030187B (en) * 2019-03-20 2021-05-14 杭州哲达智慧能源技术研究院有限公司 Pump health diagnosis method based on support vector machine
CN109827777B (en) * 2019-04-01 2020-12-18 哈尔滨理工大学 Rolling bearing fault prediction method based on partial least square method extreme learning machine
CN110132598B (en) * 2019-05-13 2020-10-09 中国矿业大学 Fault noise diagnosis algorithm for rolling bearing of rotating equipment
CN110208377A (en) * 2019-06-19 2019-09-06 南京邮电大学 A kind of more characteristic parameters damage degree assessment method based on Lamb wave
CN110595751B (en) * 2019-09-19 2020-12-18 华东理工大学 Early fault characteristic wavelet reconstruction method guided by Gini index and application thereof
CN111337256B (en) * 2020-03-27 2020-12-29 西安交通大学 Method for diagnosing fault depth local migration of rolling bearing weighted by domain asymmetry factor
CN112651426A (en) * 2020-12-01 2021-04-13 华北电力大学(保定) Fault diagnosis method for rolling bearing of wind turbine generator
CN112945556B (en) * 2021-01-26 2022-07-05 大连海事大学 Bearing fault diagnosis method based on wavelet packet decomposition and optimal rejection classification strategy
CN114894460B (en) * 2022-05-10 2023-03-28 同济大学 IMU-based method for monitoring and evaluating damage state of anti-seismic support and hanger

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
ART-KOHONEN neural network for fault diagnosis ofrotatingmachinery. B.S.Yang.Mechanical Systems Signal Processing,Vol.18 . 2004
ART-KOHONEN neural network for fault diagnosis ofrotatingmachinery. B.S.Yang.Mechanical Systems Signal Processing,Vol.18. 2004 *
Ensembling Neural Networks:Many Could Be Better Than All. Zhi-Hua Zhou.Artificial Intelligence,Vol.137 No.1-2. 2002
Ensembling Neural Networks:Many Could Be Better Than All. Zhi-Hua Zhou.Artificial Intelligence,Vol.137 No.1-2. 2002 *
一种基于SVM和EMD的齿轮故障诊断方法. 于德介等.机械工程学报,第41卷第1期. 2005
基于bagging支持向量机集成的入侵检测研究. 谷雨等.微电子学与计算机,第22卷第5期. 2005
基于bagging支持向量机集成的入侵检测研究. 谷雨等.微电子学与计算机,第22卷第5期. 2005 *
基于bagging的选择性聚类集成. 唐伟等.软件学报,第16卷第4期. 2005
基于bagging的选择性聚类集成. 唐伟等.软件学报,第16卷第4期. 2005 *
经验模式分解模糊特征提取的支持向量机混合诊断模型. 胡桥等.西安交通大学学报,第39卷第3期. 2005

Also Published As

Publication number Publication date
CN1811367A (en) 2006-08-02

Similar Documents

Publication Publication Date Title
CN100485342C (en) Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault
CN102944418B (en) Wind turbine generator group blade fault diagnosis method
CN108181107B (en) The Wind turbines bearing mechanical method for diagnosing faults of meter and more class objects
CN109635677B (en) Compound fault diagnosis method and device based on multi-label classification convolutional neural network
CN103728551B (en) A kind of analog-circuit fault diagnosis method based on cascade integrated classifier
Zhang et al. Research on bearing fault diagnosis of wind turbine gearbox based on 1DCNN-PSO-SVM
CN102944416B (en) Multi-sensor signal fusion technology-based fault diagnosis method for wind turbine blades
CN110161343B (en) Non-invasive real-time dynamic monitoring method for external powered device of intelligent train
CN108827605A (en) A kind of mechanical breakdown characteristic automatic extraction method based on improvement sparseness filtering
CN106197999A (en) A kind of planetary gear method for diagnosing faults
CN104819846A (en) Rolling bearing sound signal fault diagnosis method based on short-time Fourier transform and sparse laminated automatic encoder
CN108168682A (en) A kind of GIL On-line Faults monitoring system based on vibration signal support vector machines
CN113255078A (en) Bearing fault detection method and device under unbalanced sample condition
CN105760839A (en) Bearing fault diagnosis method based on multi-feature manifold learning and support vector machine
CN206504869U (en) A kind of rolling bearing fault diagnosis device
CN106198000A (en) A kind of rocker arm of coal mining machine gear failure diagnosing method
CN106596116A (en) Vibration fault diagnosis method of wind generating set
CN106092578A (en) A kind of machine tool mainshaft bearing confined state online test method based on wavelet packet and support vector machine
CN103033359A (en) Wind turbine generator main gear fault diagnosis method of multi-feature multilevel comprehensive judgment
CN103335617A (en) Vibrational signal based railway track geometric deformation detection method
CN106596149A (en) Method for monitoring and diagnosing flat wheel state of rail vehicle based on VMD
CN107909118A (en) A kind of power distribution network operating mode recording sorting technique based on deep neural network
CN103398843A (en) Classification method of planetary gear box sun gear faults based on multi-class relevance vector machine
CN113158984B (en) Bearing fault diagnosis method based on complex Morlet wavelet and lightweight convolution network
CN105678343A (en) Adaptive-weighted-group-sparse-representation-based diagnosis method for noise abnormity of hydroelectric generating set

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: SHENJI GROUP KUMING TOOL CO., LTD.

Free format text: FORMER OWNER: XI'AN COMMUNICATION UNIV.

Effective date: 20091023

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20091023

Address after: No 23 Dam Road, Yunnan, Kunming

Patentee after: Shenji Group Kunming Machine Tool Co., Ltd.

Address before: No. 28, Xianning Road, Xi'an, Shaanxi

Patentee before: Xi'an Jiaotong University

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090506

Termination date: 20170303