CN102520341A - Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm - Google Patents

Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm Download PDF

Info

Publication number
CN102520341A
CN102520341A CN2011103965629A CN201110396562A CN102520341A CN 102520341 A CN102520341 A CN 102520341A CN 2011103965629 A CN2011103965629 A CN 2011103965629A CN 201110396562 A CN201110396562 A CN 201110396562A CN 102520341 A CN102520341 A CN 102520341A
Authority
CN
China
Prior art keywords
fault
test
sample
value
class
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.)
Granted
Application number
CN2011103965629A
Other languages
Chinese (zh)
Other versions
CN102520341B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN 201110396562 priority Critical patent/CN102520341B/en
Publication of CN102520341A publication Critical patent/CN102520341A/en
Application granted granted Critical
Publication of CN102520341B publication Critical patent/CN102520341B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Test And Diagnosis Of Digital Computers (AREA)
  • Tests Of Electronic Circuits (AREA)

Abstract

The invention discloses an analog circuit fault diagnosis method based on a Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm, which comprises the following steps of: carrying out fault diagnosis by adopting a kernelized fuzzy C-means clustering algorithm and firstly judging whether a new fault exits in a test sample, if YES, a diagnostic model of a new fault sample is trained to join a diagnosis system, or else, the fault positioning is carried out on the test sample according to a Bayes fault classification standard. In the invention, the wavelet transform pretreatment is carried out on the fault sample, and the multi-feature fusion is carried out on the wavelet coefficient energy value and the wavelet coefficient fractal dimension value of the sample to extract fault characteristics; and the frequency of an optimal measurable node and/or a test signal is selected through taking a maximum class inter-class distance as a basis. Compared with the prior art, the analog circuit fault diagnosis method realizes that the analog circuit fault diagnosis method, the new fault of an analog circuit can be effectively diagnosed, and the diagnosis accuracy can be improved.

Description

Analog circuit fault diagnosis method based on Bayes-KFCM algorithm
Technical Field
The invention relates to a fault diagnosis method for an analog circuit, in particular to a fault diagnosis method for an analog circuit based on a Bayes-KFCM algorithm.
Background
The test and diagnosis of the analog circuit mainly aim at the functional test of the circuit, the current sequence division in the test process according to the simulation can be divided into simulation before test and simulation after test, research documents and patent data in recent years show that the intelligent fault diagnosis method is the focus and key point of the current research, the method belongs to the simulation before test, and the intelligent diagnosis method can partially solve the problems of fuzziness, uncertainty and the like of the analog circuit fault. The most common and widely applied intelligent diagnosis method is a neural network method, and in recent years, an analog circuit fault diagnosis method based on SVM has been greatly developed and becomes an important branch of the intelligent diagnosis method. However, both of the two conventional intelligent methods belong to the instructor learning algorithm, that is, the existing generic labels of the fault samples must be known when the diagnostic model is trained, but when a new fault occurs, neither of the two fault diagnosis algorithms can effectively diagnose the new fault. In addition, when a circuit generates a new fault class, the diagnosis system based on the instructor learning algorithm needs to retrain all sample sets, and the diagnosis efficiency is low. Therefore, in order to effectively diagnose new faults and also diagnose known fault classes, a fault classification algorithm based on the learning of no instructor is a suitable intelligent diagnosis algorithm.
In addition, the selection and extraction of fault characteristics are a key technology in fault diagnosis of the analog circuit, and because components of the analog circuit have tolerance characteristics, the output response of the circuit has randomness in a certain range, and the output of the circuit often presents nonlinear characteristics. Therefore, the fault characteristics of different faults of the analog circuit may have a cross or overlapping phenomenon, which becomes a fuzzy fault and cannot be correctly diagnosed and positioned. In order to solve the problem of identification of fuzzy fault characteristics, a signal processing and dimension reduction method is adopted to extract and select the characteristics of a fault response signal in the fault diagnosis of an analog circuit. The most common signal processing method is wavelet transform, which calculates the energy entropy of wavelet coefficients as fault features, and has been proven to be effective in extracting fault features. However, the overlapping and crossing phenomena of fuzzy fault features cannot be completely solved by only adopting a single fault feature, and if a plurality of kinds of feature information are selected to reflect the fault features from different sides, certain information complementarity exists. Therefore, the method for realizing the fusion of different characteristics by means of the information fusion technology is an effective way for extracting the fault characteristics of the analog circuit.
Disclosure of Invention
The invention aims to solve the technical problems that the diagnosis efficiency is low and new faults cannot be effectively diagnosed due to the fact that the existing fault diagnosis method adopts a teacher learning algorithm, and provides an analog circuit fault diagnosis method based on a Bayes-KFCM algorithm.
A fault diagnosis method for an analog circuit based on a Bayes-KFCM algorithm comprises the following steps:
a, selecting an optimal measurable node and a test signal frequency of a circuit to be tested;
b, inputting test signals to the circuit to be tested, simulating various typical fault states, and collecting voltage output values of the optimal measurable nodes to obtain fault data serving as training data; under the same test signal and measurable node, collecting data of the test circuit in the actual working state as test data;
c, respectively extracting the characteristics of the fault data and the test data, and denoising to generate a training sample set and a test sample set;
step D, training the fault diagnosis model by using the training sample set, and performing fault diagnosis on the test sample set by using the trained fault diagnosis model;
the step D specifically comprises the following steps:
step D1, training the fault diagnosis model by using a KFCM clustering algorithm, which specifically comprises the following steps: mapping the training sample set to a high-dimensional space through a kernel function; then by blurringCClustering by using a mean method, stopping the algorithm when the ratio of the number of correctly clustered samples to all the clustered samples is greater than or equal to a preset threshold value, finishing the training, taking the trained clustering model as a diagnosis model, and simultaneously obtaining the clustering centers of various training samples and the distance value of the training sample with the largest distance from the clustering center in each training sample
Figure 2011103965629100002DEST_PATH_IMAGE002
Wherein
Figure 2011103965629100002DEST_PATH_IMAGE004
nClass number of training sample;
step D2, mapping the test sample to a high-dimensional space through a kernel function, and calculating the distance from the test sample to various cluster centers in the high-dimensional space
Figure 2011103965629100002DEST_PATH_IMAGE006
Wherein
Figure 548275DEST_PATH_IMAGE004
(ii) a When in use
Figure 2011103965629100002DEST_PATH_IMAGE008
If so, the test sample is a new fault sample, clustering the new fault sample by adopting a KFCM clustering algorithm, and adding a new clustering model into the diagnosis system; otherwise, fault location is carried out on the test sample by using a Bayes fault classification criterion, wherein the Bayes fault classification criterion is as follows:
Figure 2011103965629100002DEST_PATH_IMAGE010
in the formula,
Figure 2011103965629100002DEST_PATH_IMAGE012
is the firstiThe number of training samples of a class,
Figure 2011103965629100002DEST_PATH_IMAGE014
is the number of all the training samples,is a test specimenxFrom the first to the secondiThe distance between the cluster centers of the class training samples,
Figure 2011103965629100002DEST_PATH_IMAGE018
is all thatiThe mean value of the distance between the class training sample and the cluster center thereof, Bayes fault classification and representation test samplexBelong to have the largest
Figure 2011103965629100002DEST_PATH_IMAGE020
Fault class of value.
Further, the invention adopts a method of multi-fault feature fusion and selection during feature extraction to overcome the problem that the overlapping and crossing phenomena of fuzzy fault features cannot be thoroughly solved by a single fault feature, which is specifically as follows:
the feature extraction in the step C is specifically carried out according to the following method:
step C1, carrying out multilayer wavelet decomposition on the acquired voltage value to decompose the voltage value into detail coefficients and approximation coefficients;
step C2, calculating the energy entropy of each layer of detail coefficients and approximation coefficients, and using a vector formed by the energy entropy of the multilayer wavelet coefficients as a first feature representation of the voltage signal;
step C3, calculating the fractal dimension value of each layer of detail coefficient and approximation coefficient, and using the vector formed by the fractal dimension values of the multi-layer wavelet coefficients as a second feature representation of the voltage signal;
and step C4, fusing the two characteristics by adopting a linear summation method, wherein the linear summation fusion formula is expressed as follows:
Figure 2011103965629100002DEST_PATH_IMAGE022
in the formula,a feature vector representing the fusion is represented by,
Figure 2011103965629100002DEST_PATH_IMAGE026
and
Figure 2011103965629100002DEST_PATH_IMAGE028
respectively representing an energy entropy value and a fractal dimension value which are obtained by calculating a layer of wavelet decomposition coefficient D of the signal,
Figure 2011103965629100002DEST_PATH_IMAGE030
and
Figure 2011103965629100002DEST_PATH_IMAGE032
respectively representing the wavelet coefficient energy value and the wavelet coefficient fractal dimension value in the fusionA weight, and
step C5, calculating fusion characteristic vector
Figure 672874DEST_PATH_IMAGE024
The total correlation magnitude of each feature and the rest features in the database are sorted from high to low, and the cumulative contribution value of the correlation of the selected features is larger than that of the correlation of the other featuresK% of the total amount ofhAnd (4) carrying out feature dimension reduction on each feature,Kthe value range of (1) is (0, 100); sample setting
Figure 2011103965629100002DEST_PATH_IMAGE036
Therein containPSamples, each sample having a feature dimension ofWThen sample
Figure 327977DEST_PATH_IMAGE036
The characteristic correlation coefficient vector is obtained according to the following method:
first, a sample is calculated
Figure 97088DEST_PATH_IMAGE036
The sample mean vector of (d) is:
Figure 2011103965629100002DEST_PATH_IMAGE038
then, the samples are calculated
Figure 786827DEST_PATH_IMAGE036
The covariance matrix of (1), wherein the character is represented by
Figure 2011103965629100002DEST_PATH_IMAGE040
Represents each element in the covariance matrix:
Figure 2011103965629100002DEST_PATH_IMAGE042
Figure 2011103965629100002DEST_PATH_IMAGE044
then, according to the matrixCalculating a correlation matrix:
Figure 2011103965629100002DEST_PATH_IMAGE048
finally, the rest of each feature pair is calculated respectivelyW-1 sum of correlation coefficients of features:
Figure 2011103965629100002DEST_PATH_IMAGE050
the feature dimension is thenWOf (2) a sample
Figure 45507DEST_PATH_IMAGE036
The feature correlation coefficient vector of (a) is:
Figure 2011103965629100002DEST_PATH_IMAGE052
wherein the feature-dependent cumulative contribution value
Figure 2011103965629100002DEST_PATH_IMAGE054
Calculated according to the following formula:
Figure 2011103965629100002DEST_PATH_IMAGE056
furthermore, in order to make the selection of the optimal measurable node and/or the test signal frequency more representative and thus improve the fault diagnosis, the invention also adopts the maximum intra-class inter-class distance as the selection basis to select the optimal measurable node and/or the test signal frequency, which specifically comprises the following steps:
the test signal frequency is selected by:
a1, obtaining an amplitude-frequency response curve of a circuit to be tested;
a2, selecting an inflection point on an amplitude-frequency response curve and frequencies near the inflection point as a frequency set to be selected;
step A3, simulating some typical test faults manually, collecting the response voltage value of the circuit to be tested under the excitation of all the frequency to be tested at the output end of the circuit as a fault sample value, calculating the intra-class inter-class distance of different test fault class samples, and selecting the frequency to be tested with the maximum intra-class inter-class distance of the test fault class as the test frequency.
The optimal measurable node is selected by the following method:
step A4, taking all testable test nodes in the circuit to be tested as test nodes to be selected, manually simulating some typical test faults, taking selected test signals as excitation sources to be loaded to the circuit to be tested, and collecting voltage values of all test faults on all test nodes to be selected as fault sample values;
step A5, calculating the inter-class distance of the test fault class sample in each test node to be selected, and selecting the front point with the maximum inter-class distance of all fault classesMA plurality of test nodes, each of which is connected to a test node,Mthe test nodes are preset integers which are smaller than the total number of the test nodes to be selected.
Compared with the prior art, the analog circuit fault diagnosis method has the following beneficial effects:
(1) the method has the advantages that two characteristics of the energy entropy and the fractal dimension after wavelet transformation are extracted simultaneously, characteristic fusion is carried out through a linear summation method, the characteristics of fault signals are extracted from various side faces, certain complementarity is achieved, and the one-sidedness of single characteristics is compensated.
(2) Selecting the top with the accumulated contribution value of the feature correlation larger than 90 percent by calculating the correlation size of each feature component in the feature vector and sorting the correlation sizeshThe characteristic can effectively remove redundant and irrelevant characteristics in the characteristic vector, reduce the dimension of the fault characteristic vector and improve the efficiency of a diagnosis algorithm.
(3) A kernel fuzzy C-means clustering (KFCM) algorithm based on instructor-free learning is adopted as a fault diagnosis algorithm, so that the limitation that the traditional instructor-based learning algorithm needs a fault class label is overcome, and a new fault class can be effectively diagnosed. In addition, in the face of new fault classes, the traditional intelligent diagnosis system needs to retrain all sample sets, diagnosis efficiency is low, and only a diagnosis model of the new fault classes needs to be trained by adopting a guide-free learning algorithm and is added into the diagnosis system.
(4) A fault classification criterion is constructed according to the characteristics of a KFCM algorithm and a Bayes decision theory, and the criterion can quickly and accurately judge and position faults.
Drawings
FIG. 1 is a flow chart of an analog circuit fault diagnostic method of the present invention;
FIG. 2 is a flow chart of a method of multiple fault signature fusion and selection in the method of the present invention.
Detailed Description
The technical scheme of the invention is explained in detail in the following with the accompanying drawings:
the analog circuit fault diagnosis method of the present invention, as shown in fig. 1, includes the steps of:
and A, selecting the optimal measurable node and the test signal frequency of the circuit to be tested.
In order to enable the selection of the optimal measurable node and/or the optimal frequency of the test signal to be more representative and improve the fault diagnosis, the method adopts the maximum intra-class inter-class distance as a selection basis to select the optimal measurable node and/or the optimal frequency of the test signal. Specifically, step a specifically includes:
a1, obtaining an amplitude-frequency response curve of a circuit to be tested;
a2, selecting an inflection point on an amplitude-frequency response curve and frequencies near the inflection point as a frequency set to be selected;
a3, simulating some typical test faults manually, collecting the response voltage values of the circuit to be tested under the excitation of all the to-be-selected frequencies at the output end of the circuit as fault sample values, calculating the intra-class inter-class distances of different test fault class samples, and selecting the to-be-selected frequency with the maximum intra-class inter-class distance of the test fault class as the test frequency; wherein, the inter-class distance
Figure DEST_PATH_IMAGE058
The calculation of (a) is the prior art, and the calculation formula is as follows:
Figure 2011103965629100002DEST_PATH_IMAGE060
wherein
Figure DEST_PATH_IMAGE064
Figure DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE068
In the formula (I), whereincAs the number of the categories,n i is composed ofiThe number of samples of a class,P i is the firstiThe prior probability of a class sample,are respectively asiThe feature vector of a class is then calculated,
Figure DEST_PATH_IMAGE072
is shown asiThe mean vector of the class sample set is,mrepresents the sample lumped mean vector of all classes,
Figure DEST_PATH_IMAGE074
referred to as the inter-class dispersion matrix,
Figure DEST_PATH_IMAGE076
the intra-class dispersion matrix is generally considered that the inter-class dispersion is as large as possible, and the intra-class dispersion is as small as possible, so that classification is facilitated;
step A4, taking all testable test nodes in the circuit to be tested as test nodes to be selected, manually simulating some typical test faults, taking selected test signals as excitation sources to be loaded to the circuit to be tested, and collecting voltage values of all test faults on all test nodes to be selected as fault sample values;
step A5, calculating the inter-class distance of the test fault class sample in each test node to be selected, and selecting the front point with the maximum inter-class distance of all fault classesMA plurality of test nodes, each of which is connected to a test node,Mto prepareFirstly, setting an integer less than the total number of the test nodes to be selected.
B, inputting test signals to the circuit to be tested, simulating various typical fault states, and collecting voltage output values of the optimal measurable nodes to obtain fault data serving as training data; under the same test signal and measurable node, data of the test circuit in the actual working state is collected as test data.
And C, respectively extracting the characteristics of the fault data and the test data, denoising, and generating a training sample set and a test sample set.
The invention adopts a method of multi-fault feature fusion and selection during feature extraction to overcome the problem that the overlapping and crossing phenomena of fuzzy fault features cannot be thoroughly solved by a single fault feature, and specifically as shown in fig. 2, the method comprises the following steps:
step C1, carrying out multilayer wavelet decomposition on the acquired voltage value to decompose the voltage value into detail coefficients and approximation coefficients;
step C2, calculating the energy entropy of each layer of detail coefficients and approximation coefficients, and using a vector formed by the energy entropy of the multilayer wavelet coefficients as a first feature representation of the voltage signal;
step C3, calculating the fractal dimension value of each layer of detail coefficient and approximation coefficient, and using the vector formed by the fractal dimension values of the multi-layer wavelet coefficients as a second feature representation of the voltage signal;
and step C4, fusing the two characteristics by adopting a linear summation method, wherein the linear summation fusion formula is expressed as follows:
Figure 371928DEST_PATH_IMAGE022
in the formula,
Figure 882412DEST_PATH_IMAGE024
a feature vector representing the fusion is represented by,
Figure 613608DEST_PATH_IMAGE026
and
Figure 420021DEST_PATH_IMAGE028
respectively representing an energy entropy value and a fractal dimension value which are obtained by calculating a layer of wavelet decomposition coefficient D of the signal,
Figure 236667DEST_PATH_IMAGE030
and
Figure 285264DEST_PATH_IMAGE032
respectively represent the weight of wavelet coefficient energy value and wavelet coefficient fractal dimension value in fusion, and
Figure 870966DEST_PATH_IMAGE034
(ii) a In this embodiment, the weight
Figure 113860DEST_PATH_IMAGE030
And
Figure 417802DEST_PATH_IMAGE032
all values of (A) are 0.5;
step C5, calculating fusion characteristic vector
Figure 4510DEST_PATH_IMAGE024
The total correlation magnitude of each feature and the rest features in the database are sorted from high to low, and the cumulative contribution value of the correlation of the selected features is larger than that of the correlation of the other featuresK% of the total amount ofhAnd (4) carrying out feature dimension reduction on each feature,Kis (0, 100), in the present embodiment, the top of the feature correlation cumulative contribution value greater than 90% is selectedhA feature; sample setting
Figure 710298DEST_PATH_IMAGE036
Therein containPSamples, each sample having a feature dimension ofWThen sample
Figure 858513DEST_PATH_IMAGE036
The characteristic correlation coefficient vector is obtained according to the following method:
first, a sample is calculated
Figure 649752DEST_PATH_IMAGE036
The sample mean vector of (d) is:
Figure 40151DEST_PATH_IMAGE038
then, the samples are calculated
Figure 600445DEST_PATH_IMAGE036
The covariance matrix of (1), wherein the character is represented by
Figure 247458DEST_PATH_IMAGE040
Represents each element in the covariance matrix:
Figure 198097DEST_PATH_IMAGE042
Figure 392187DEST_PATH_IMAGE044
then, according to the matrix
Figure 806987DEST_PATH_IMAGE046
Calculating a correlation matrix:
Figure 624902DEST_PATH_IMAGE048
finally, the rest of each feature pair is calculated respectivelyW-1 sum of correlation coefficients of features:
Figure 374421DEST_PATH_IMAGE050
the feature dimension is thenWOf (2) a sample
Figure 857355DEST_PATH_IMAGE036
The feature correlation coefficient vector of (a) is:
Figure 142974DEST_PATH_IMAGE052
wherein the feature-dependent cumulative contribution value
Figure 318740DEST_PATH_IMAGE054
Calculated according to the following formula:
Figure 555555DEST_PATH_IMAGE056
in this embodiment, the denoising is performed by using a conventional wavelet transform soft threshold method.
And D, training the fault diagnosis model by using the training sample set, and performing fault diagnosis on the test sample set by using the trained fault diagnosis model.
The invention provides a new algorithm combining a KFCM algorithm and a Bayes classification criterion, which is called Bayes-KFCM algorithm in the invention.
KFCM is a guiding-free learning algorithm, original characteristic space data are mapped to a high-dimensional space through a kernel function, and then clustering is carried out through a fuzzy C mean value method. Suppose there are training sampleskClass failure, then the number of clusters iskAnd taking a certain preset accuracy as a condition for stopping the KFCM clustering algorithm, wherein the accuracy is calculated by the following formula:
Figure DEST_PATH_IMAGE078
the objective function of the KFCM algorithm is:
Figure DEST_PATH_IMAGE080
the constraint conditions are as follows:
Figure DEST_PATH_IMAGE082
Figure DEST_PATH_IMAGE086
Figure DEST_PATH_IMAGE088
wherein,
Figure DEST_PATH_IMAGE090
Figure DEST_PATH_IMAGE092
in the formula
Figure DEST_PATH_IMAGE094
Is a membership matrix which is a matrix of membership,the degree of membership is represented by,is a weighted index of the number of bits in the bit stream,v j is the center of the cluster in the input space,cis the number of categories of the cluster and,the non-linear mapping is carried out,Kis a nuclear letterAnd (4) counting.
According to the Bayes optimal decision principle, a new sample should be classified into the category with the optimal posterior probability:
Figure DEST_PATH_IMAGE102
wherein the prior probability
Figure DEST_PATH_IMAGE104
Class correlation density
Figure DEST_PATH_IMAGE106
Can be determined by a pseudo-density function
Figure DEST_PATH_IMAGE108
So as to obtain the compound with the characteristics of,
Figure DEST_PATH_IMAGE110
hence posterior probability
Figure 2011103965629100002DEST_PATH_IMAGE112
Equal to:
Figure DEST_PATH_IMAGE114
wherein
Figure 2011103965629100002DEST_PATH_IMAGE116
Is usually a constant, wherein
Figure 987322DEST_PATH_IMAGE012
Is the firstiThe number of training samples of a class,is the number of all the training samples,is a test specimenxFrom the first to the secondiThe distance between the centers of the fault-like clusters,
Figure 932647DEST_PATH_IMAGE018
is all thatiThe summed average of class training samples and their cluster center distances. The classification criterion based on Bayes' optimal decision principle described above can be expressed in the form:
Figure 2011103965629100002DEST_PATH_IMAGE118
specifically, step B specifically includes:
step D1, training the fault diagnosis model by using a KFCM clustering algorithm, which specifically comprises the following steps: mapping the training sample set to a high-dimensional space through a kernel function; then by blurringCClustering by using a mean method, stopping the algorithm when the ratio of the number of correctly clustered samples to all the clustered samples is greater than or equal to a preset threshold value, finishing the training, taking the trained clustering model as a diagnosis model, and simultaneously obtaining the clustering centers of various training samples and the distance value of the training sample with the largest distance from the clustering center in each training sample
Figure 396864DEST_PATH_IMAGE002
Wherein
Figure 391496DEST_PATH_IMAGE004
nClass number of training sample;
step D2, mapping the test sample to a high-dimensional space through a kernel function, and calculating the distance from the test sample to various cluster centers in the high-dimensional space
Figure 643486DEST_PATH_IMAGE006
Wherein
Figure 854893DEST_PATH_IMAGE004
(ii) a When in use
Figure 483321DEST_PATH_IMAGE008
If so, the test sample is a new fault sample, clustering the new fault sample by adopting a KFCM clustering algorithm, and adding a new clustering model into the diagnosis system; otherwise, fault location is carried out on the test sample by using a Bayes fault classification criterion, wherein the Bayes fault classification criterion is as follows:
Figure 598038DEST_PATH_IMAGE010
in the formula,is the firstiThe number of training samples of a class,
Figure 719633DEST_PATH_IMAGE014
is the number of all the training samples,
Figure 151751DEST_PATH_IMAGE016
is a test specimenxFrom the first to the secondiThe distance between the cluster centers of the class training samples,
Figure 229298DEST_PATH_IMAGE018
is all thatiThe mean value of the distance between the class training sample and the cluster center thereof, Bayes fault classification and representation test samplexBelong to have the largest
Figure 88669DEST_PATH_IMAGE020
Fault class of value.

Claims (7)

1. A fault diagnosis method for an analog circuit based on a Bayes-KFCM algorithm comprises the following steps:
a, selecting an optimal measurable node and a test signal frequency of a circuit to be tested;
b, inputting test signals to the circuit to be tested, simulating various typical fault states, and collecting voltage output values of the optimal measurable nodes to obtain fault data serving as training data; under the same test signal and measurable node, collecting data of the test circuit in the actual working state as test data;
c, respectively extracting the characteristics of the fault data and the test data, and denoising to generate a training sample set and a test sample set;
step D, training the fault diagnosis model by using the training sample set, and performing fault diagnosis on the test sample set by using the trained fault diagnosis model;
the method is characterized in that the step D specifically comprises the following steps:
step D1, training the fault diagnosis model by using a KFCM clustering algorithm, which specifically comprises the following steps: mapping the training sample set to a high-dimensional space through a kernel function; then by blurringCClustering by using a mean method, stopping the algorithm when the ratio of the number of correctly clustered samples to all the clustered samples is greater than or equal to a preset threshold value, finishing the training, taking the trained clustering model as a diagnosis model, and simultaneously obtaining the clustering centers of various training samples and the distance value of the training sample with the largest distance from the clustering center in each training sample
Figure 2011103965629100001DEST_PATH_IMAGE002
Wherein
Figure 2011103965629100001DEST_PATH_IMAGE004
nClass number of training sample;
step D2, mapping the test sample to a high-dimensional space through a kernel function, and calculating the distance from the test sample to various cluster centers in the high-dimensional space
Figure 2011103965629100001DEST_PATH_IMAGE006
Wherein
Figure 791804DEST_PATH_IMAGE004
(ii) a When in use
Figure 2011103965629100001DEST_PATH_IMAGE008
If so, the test sample is newThe fault sample is clustered by adopting a KFCM clustering algorithm on the new fault sample, and a new clustering model is added into the diagnosis system; otherwise, fault location is carried out on the test sample by using a Bayes fault classification criterion, wherein the Bayes fault classification criterion is as follows:
Figure 2011103965629100001DEST_PATH_IMAGE010
in the formula,
Figure 2011103965629100001DEST_PATH_IMAGE012
is the firstiThe number of training samples of a class,
Figure 2011103965629100001DEST_PATH_IMAGE014
is the number of all the training samples,
Figure 2011103965629100001DEST_PATH_IMAGE016
is a test specimenxFrom the first to the secondiThe distance between the cluster centers of the class training samples,
Figure 2011103965629100001DEST_PATH_IMAGE018
is all thatiThe mean value of the distance between the class training sample and the cluster center thereof, Bayes fault classification and representation test samplexBelong to have the largest
Figure 2011103965629100001DEST_PATH_IMAGE020
Fault class of value.
2. The analog circuit fault diagnosis method based on the Bayes-KFCM algorithm according to claim 1, wherein the value of the threshold in step D1 is 90%.
3. The analog circuit fault diagnosis method based on the Bayes-KFCM algorithm as claimed in claim 1, wherein said feature extraction in step C is specifically according to the following method:
step C1, carrying out multilayer wavelet decomposition on the acquired voltage value to decompose the voltage value into detail coefficients and approximation coefficients;
step C2, calculating the energy entropy of each layer of detail coefficients and approximation coefficients, and using a vector formed by the energy entropy of the multilayer wavelet coefficients as a first feature representation of the voltage signal;
step C3, calculating the fractal dimension value of each layer of detail coefficient and approximation coefficient, and using the vector formed by the fractal dimension values of the multi-layer wavelet coefficients as a second feature representation of the voltage signal;
and step C4, fusing the two characteristics by adopting a linear summation method, wherein the linear summation fusion formula is expressed as follows:
Figure 2011103965629100001DEST_PATH_IMAGE022
in the formula,a feature vector representing the fusion is represented by,
Figure 2011103965629100001DEST_PATH_IMAGE026
and
Figure 2011103965629100001DEST_PATH_IMAGE028
respectively representing an energy entropy value and a fractal dimension value which are obtained by calculating a layer of wavelet decomposition coefficient D of the signal,and
Figure 2011103965629100001DEST_PATH_IMAGE032
respectively represent the weight of wavelet coefficient energy value and wavelet coefficient fractal dimension value in fusion, and
Figure DEST_PATH_IMAGE034
step C5, calculating fusion characteristic vector
Figure 585536DEST_PATH_IMAGE024
The total correlation magnitude of each feature and the rest features in the database are sorted from high to low, and the cumulative contribution value of the correlation of the selected features is larger than that of the correlation of the other featuresK% of the total amount ofhAnd (4) carrying out feature dimension reduction on each feature,Kthe value range of (1) is (0, 100); sample setting
Figure DEST_PATH_IMAGE036
Therein containPSamples, each sample having a feature dimension ofWThen sampleThe characteristic correlation coefficient vector is obtained according to the following method:
first, a sample is calculated
Figure 813441DEST_PATH_IMAGE036
The sample mean vector of (d) is:
Figure DEST_PATH_IMAGE038
then, the samples are calculated
Figure 847256DEST_PATH_IMAGE036
The covariance matrix of (1), wherein the character is represented by
Figure DEST_PATH_IMAGE040
Represents each element in the covariance matrix:
Figure DEST_PATH_IMAGE042
Figure DEST_PATH_IMAGE044
then, according to the matrix
Figure DEST_PATH_IMAGE046
Calculating a correlation matrix:
Figure DEST_PATH_IMAGE048
finally, the rest of each feature pair is calculated respectivelyW-1 sum of correlation coefficients of features:
Figure DEST_PATH_IMAGE050
the feature dimension is thenWOf (2) a sample
Figure 18516DEST_PATH_IMAGE036
The feature correlation coefficient vector of (a) is:
Figure DEST_PATH_IMAGE052
wherein the feature-dependent cumulative contribution valueCalculated according to the following formula:
Figure DEST_PATH_IMAGE056
4. the analog circuit fault diagnosis method based on Bayes-KFCM algorithm as claimed in claim 3, wherein the weight occupied by wavelet coefficient energy value and wavelet coefficient fractal dimension value in feature fusion
Figure 132971DEST_PATH_IMAGE030
And
Figure 334145DEST_PATH_IMAGE032
all values of (A) are 0.5.
5. The analog circuit fault diagnosis method based on Bayes-KFCM algorithm as claimed in claim 3,Kis at a value of 90.
6. The analog circuit fault diagnosis method based on the Bayes-KFCM algorithm according to any of claims 1 to 5, characterized in that the test signal frequency is selected by:
a1, obtaining an amplitude-frequency response curve of a circuit to be tested;
a2, selecting an inflection point on an amplitude-frequency response curve and frequencies near the inflection point as a frequency set to be selected;
step A3, simulating some typical test faults manually, collecting the response voltage value of the circuit to be tested under the excitation of all the frequency to be tested at the output end of the circuit as a fault sample value, calculating the intra-class inter-class distance of different test fault class samples, and selecting the frequency to be tested with the maximum intra-class inter-class distance of the test fault class as the test frequency.
7. A method of fault diagnosis for an analog circuit based on the Bayes-KFCM algorithm according to any of claims 1-5 wherein the optimal measurable node is selected by:
step A4, taking all testable test nodes in the circuit to be tested as test nodes to be selected, manually simulating some typical test faults, taking selected test signals as excitation sources to be loaded to the circuit to be tested, and collecting voltage values of all test faults on all test nodes to be selected as fault sample values;
step A5, calculating the inter-class distance of the test fault class samples in each test node to be selected, and selecting all fault classesFront of maximum distance between classesMA plurality of test nodes, each of which is connected to a test node,Mthe test nodes are preset integers which are smaller than the total number of the test nodes to be selected.
CN 201110396562 2011-12-05 2011-12-05 Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm Expired - Fee Related CN102520341B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110396562 CN102520341B (en) 2011-12-05 2011-12-05 Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110396562 CN102520341B (en) 2011-12-05 2011-12-05 Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm

Publications (2)

Publication Number Publication Date
CN102520341A true CN102520341A (en) 2012-06-27
CN102520341B CN102520341B (en) 2013-12-18

Family

ID=46291327

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110396562 Expired - Fee Related CN102520341B (en) 2011-12-05 2011-12-05 Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm

Country Status (1)

Country Link
CN (1) CN102520341B (en)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103198217A (en) * 2013-03-26 2013-07-10 X·Q·李 Fault detection method and system
CN104090227A (en) * 2014-07-16 2014-10-08 电子科技大学 Measuring point selection method in analog circuit fault diagnosis
CN104502834A (en) * 2014-12-30 2015-04-08 南京农业大学 OFRF (output frequency response function) modeling based nonlinear analog circuit fault parameter identification method
CN104573740A (en) * 2014-12-22 2015-04-29 山东鲁能软件技术有限公司 SVM classification model-based equipment fault diagnosing method
CN104635146A (en) * 2015-02-06 2015-05-20 南京农业大学 Analog circuit fault diagnosis method based on random sinusoidal signal test and HMM (Hidden Markov Model)
CN104849650A (en) * 2015-05-19 2015-08-19 重庆大学 Analog circuit fault diagnosis method based on improvement
CN105044562A (en) * 2015-08-26 2015-11-11 云南电网有限责任公司电力科学研究院 Distribution network fault positioning method based on Bayes formula algorithm
CN105652182A (en) * 2015-12-28 2016-06-08 北京航天测控技术有限公司 Circuit board fault positioning system and circuit board fault positioning method based on circuit network and graph search
WO2017024692A1 (en) * 2015-08-07 2017-02-16 合肥工业大学 Method of simulating analog circuit diagnostic fault using single measurement node
CN106444468A (en) * 2016-05-03 2017-02-22 重庆青年职业技术学院 Information machine adapter detection system and method
CN106896306A (en) * 2017-04-26 2017-06-27 国网上海市电力公司 A kind of GIS oscillatory surges pressure test signal antinoise method
CN107220670A (en) * 2017-05-27 2017-09-29 重庆大学 Supervised Artifical Taste system features extracting method is had based on wavelet transform
CN107843779A (en) * 2017-09-26 2018-03-27 国网河南省电力公司民权县供电公司 A kind of Power System Fault Record classifying and analyzing method and system based on fuzzy clustering
CN105699882B (en) * 2016-01-22 2018-05-29 电子科技大学 Analog-circuit fault diagnosis method based on vibration test technology
CN108470462A (en) * 2018-05-18 2018-08-31 上海会为智能技术有限公司 A kind of signal lamp failure detection method, device, equipment and medium
CN109726770A (en) * 2019-02-19 2019-05-07 金陵科技学院 A kind of analog circuit fault testing and diagnosing method
CN109885018A (en) * 2019-01-18 2019-06-14 广东寰球智能科技有限公司 A kind of management method of the expert think tank system towards mental package equipment
CN110579709A (en) * 2019-08-30 2019-12-17 西南交通大学 fault diagnosis method for proton exchange membrane fuel cell for tramcar
CN111307453A (en) * 2020-03-20 2020-06-19 朗斯顿科技(北京)有限公司 Transmission system fault diagnosis method based on multi-information fusion
CN111523594A (en) * 2020-04-23 2020-08-11 湖州师范学院 Improved KNN fault classification method based on LDA-KMEDOIDS
CN113159088A (en) * 2021-01-07 2021-07-23 中国地质大学(武汉) Fault monitoring and diagnosis method based on multi-feature fusion and width learning
CN113409284A (en) * 2021-06-28 2021-09-17 北京百度网讯科技有限公司 Circuit board fault detection method, device, equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614787A (en) * 2009-07-07 2009-12-30 南京航空航天大学 Analogical Electronics method for diagnosing faults based on M-ary textural classification device
CN101871994A (en) * 2010-06-11 2010-10-27 南京航空航天大学 Method for diagnosing faults of analog circuit of multi-fractional order information fusion
CN101907681A (en) * 2010-07-15 2010-12-08 南京航空航天大学 Analog circuit dynamic online failure diagnosing method based on GSD-SVDD

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614787A (en) * 2009-07-07 2009-12-30 南京航空航天大学 Analogical Electronics method for diagnosing faults based on M-ary textural classification device
CN101871994A (en) * 2010-06-11 2010-10-27 南京航空航天大学 Method for diagnosing faults of analog circuit of multi-fractional order information fusion
CN101907681A (en) * 2010-07-15 2010-12-08 南京航空航天大学 Analog circuit dynamic online failure diagnosing method based on GSD-SVDD

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103198217B (en) * 2013-03-26 2016-06-22 X·Q·李 A kind of fault detection method and system
CN103198217A (en) * 2013-03-26 2013-07-10 X·Q·李 Fault detection method and system
CN104090227A (en) * 2014-07-16 2014-10-08 电子科技大学 Measuring point selection method in analog circuit fault diagnosis
CN104090227B (en) * 2014-07-16 2016-10-26 电子科技大学 A kind of measuring point system of selection in analog circuit fault diagnosing
CN104573740B (en) * 2014-12-22 2018-07-03 山东鲁能软件技术有限公司 A kind of equipment fault diagnosis method based on svm classifier model
CN104573740A (en) * 2014-12-22 2015-04-29 山东鲁能软件技术有限公司 SVM classification model-based equipment fault diagnosing method
CN104502834A (en) * 2014-12-30 2015-04-08 南京农业大学 OFRF (output frequency response function) modeling based nonlinear analog circuit fault parameter identification method
CN104635146B (en) * 2015-02-06 2017-04-19 南京农业大学 Analog circuit fault diagnosis method based on random sinusoidal signal test and HMM (Hidden Markov Model)
CN104635146A (en) * 2015-02-06 2015-05-20 南京农业大学 Analog circuit fault diagnosis method based on random sinusoidal signal test and HMM (Hidden Markov Model)
CN104849650A (en) * 2015-05-19 2015-08-19 重庆大学 Analog circuit fault diagnosis method based on improvement
CN104849650B (en) * 2015-05-19 2018-03-02 重庆大学 One kind is based on improved analog-circuit fault diagnosis method
US10539613B2 (en) 2015-08-07 2020-01-21 Hefei University Of Technology Analog circuit fault diagnosis method using single testable node
WO2017024692A1 (en) * 2015-08-07 2017-02-16 合肥工业大学 Method of simulating analog circuit diagnostic fault using single measurement node
CN105044562A (en) * 2015-08-26 2015-11-11 云南电网有限责任公司电力科学研究院 Distribution network fault positioning method based on Bayes formula algorithm
CN105652182B (en) * 2015-12-28 2018-10-02 北京航天测控技术有限公司 A kind of board failure positioning system and method based on circuit network and graph search
CN105652182A (en) * 2015-12-28 2016-06-08 北京航天测控技术有限公司 Circuit board fault positioning system and circuit board fault positioning method based on circuit network and graph search
CN105699882B (en) * 2016-01-22 2018-05-29 电子科技大学 Analog-circuit fault diagnosis method based on vibration test technology
CN106444468A (en) * 2016-05-03 2017-02-22 重庆青年职业技术学院 Information machine adapter detection system and method
CN106444468B (en) * 2016-05-03 2023-09-05 重庆青年职业技术学院 Information machine adapter detection system and method
CN106896306B (en) * 2017-04-26 2019-04-30 国网上海市电力公司 A kind of GIS oscillatory surge pressure test signal antinoise method
CN106896306A (en) * 2017-04-26 2017-06-27 国网上海市电力公司 A kind of GIS oscillatory surges pressure test signal antinoise method
CN107220670A (en) * 2017-05-27 2017-09-29 重庆大学 Supervised Artifical Taste system features extracting method is had based on wavelet transform
CN107220670B (en) * 2017-05-27 2020-07-14 重庆大学 Method for extracting characteristics of supervised artificial taste system based on discrete wavelet transform
CN107843779A (en) * 2017-09-26 2018-03-27 国网河南省电力公司民权县供电公司 A kind of Power System Fault Record classifying and analyzing method and system based on fuzzy clustering
CN108470462A (en) * 2018-05-18 2018-08-31 上海会为智能技术有限公司 A kind of signal lamp failure detection method, device, equipment and medium
CN109885018A (en) * 2019-01-18 2019-06-14 广东寰球智能科技有限公司 A kind of management method of the expert think tank system towards mental package equipment
CN109885018B (en) * 2019-01-18 2020-10-09 广东寰球智能科技有限公司 Management method of expert intelligent library system for metal packaging equipment
CN109726770A (en) * 2019-02-19 2019-05-07 金陵科技学院 A kind of analog circuit fault testing and diagnosing method
CN110579709A (en) * 2019-08-30 2019-12-17 西南交通大学 fault diagnosis method for proton exchange membrane fuel cell for tramcar
CN110579709B (en) * 2019-08-30 2021-04-13 西南交通大学 Fault diagnosis method for proton exchange membrane fuel cell for tramcar
CN111307453B (en) * 2020-03-20 2021-11-12 朗斯顿科技(北京)有限公司 Transmission system fault diagnosis method based on multi-information fusion
CN111307453A (en) * 2020-03-20 2020-06-19 朗斯顿科技(北京)有限公司 Transmission system fault diagnosis method based on multi-information fusion
CN111523594A (en) * 2020-04-23 2020-08-11 湖州师范学院 Improved KNN fault classification method based on LDA-KMEDOIDS
CN111523594B (en) * 2020-04-23 2024-08-02 湖州师范学院 Improved KNN fault classification method based on LDA-KMEDOIDS
CN113159088B (en) * 2021-01-07 2022-07-15 中国地质大学(武汉) Fault monitoring and diagnosis method based on multi-feature fusion and width learning
CN113159088A (en) * 2021-01-07 2021-07-23 中国地质大学(武汉) Fault monitoring and diagnosis method based on multi-feature fusion and width learning
CN113409284A (en) * 2021-06-28 2021-09-17 北京百度网讯科技有限公司 Circuit board fault detection method, device, equipment and storage medium
CN113409284B (en) * 2021-06-28 2023-10-20 北京百度网讯科技有限公司 Circuit board fault detection method, device, equipment and storage medium

Also Published As

Publication number Publication date
CN102520341B (en) 2013-12-18

Similar Documents

Publication Publication Date Title
CN102520341B (en) Analog circuit fault diagnosis method based on Bayes-KFCM (Kernelized Fuzzy C-Means) algorithm
CN110287983B (en) Single-classifier anomaly detection method based on maximum correlation entropy deep neural network
Ishida et al. Deep learning-based wafer-map failure pattern recognition framework
CN103728551B (en) A kind of analog-circuit fault diagnosis method based on cascade integrated classifier
CN108737406A (en) A kind of detection method and system of abnormal flow data
CN106124175A (en) A kind of compressor valve method for diagnosing faults based on Bayesian network
CN109389171B (en) Medical image classification method based on multi-granularity convolution noise reduction automatic encoder technology
CN105223475B (en) Shelf depreciation chromatogram characteristic algorithm for pattern recognition based on Gaussian parameter fitting
CN109165604A (en) The recognition methods of non-intrusion type load and its test macro based on coorinated training
CN108333468B (en) The recognition methods of bad data and device under a kind of active power distribution network
CN108647707B (en) Probabilistic neural network creation method, failure diagnosis method and apparatus, and storage medium
CN105606914A (en) IWO-ELM-based Aviation power converter fault diagnosis method
CN102435910A (en) Power electronic circuit health monitoring method based on support vector classification
CN113516228B (en) Network anomaly detection method based on deep neural network
CN109214444B (en) Game anti-addiction determination system and method based on twin neural network and GMM
CN108171119B (en) SAR image change detection method based on residual error network
CN113657491A (en) Neural network design method for signal modulation type recognition
CN112861066B (en) Machine learning and FFT (fast Fourier transform) -based blind source separation information source number parallel estimation method
CN115600088A (en) Distribution transformer fault diagnosis method based on vibration signals
CN111310719B (en) Unknown radiation source individual identification and detection method
CN116400168A (en) Power grid fault diagnosis method and system based on depth feature clustering
CN110458071B (en) DWT-DFPA-GBDT-based optical fiber vibration signal feature extraction and classification method
CN116527346A (en) Threat node perception method based on deep learning graph neural network theory
CN112085062A (en) Wavelet neural network-based abnormal energy consumption positioning method
CN104573728B (en) A kind of texture classifying method based on ExtremeLearningMachine

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
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: 20131218

Termination date: 20161205