[go: up one dir, main page]

CN116720103B - Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density - Google Patents

Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density

Info

Publication number
CN116720103B
CN116720103B CN202310611976.1A CN202310611976A CN116720103B CN 116720103 B CN116720103 B CN 116720103B CN 202310611976 A CN202310611976 A CN 202310611976A CN 116720103 B CN116720103 B CN 116720103B
Authority
CN
China
Prior art keywords
fault
cloud
under
probability density
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.)
Active
Application number
CN202310611976.1A
Other languages
Chinese (zh)
Other versions
CN116720103A (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.)
Ningbo Lidou Intelligent Technology Co ltd
Original Assignee
Ningbo Lidou Intelligent Technology Co ltd
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 Ningbo Lidou Intelligent Technology Co ltd filed Critical Ningbo Lidou Intelligent Technology Co ltd
Priority to CN202310611976.1A priority Critical patent/CN116720103B/en
Publication of CN116720103A publication Critical patent/CN116720103A/en
Application granted granted Critical
Publication of CN116720103B publication Critical patent/CN116720103B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/28Testing of electronic circuits, e.g. by signal tracer
    • G01R31/2851Testing of integrated circuits [IC]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computational Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Probability & Statistics with Applications (AREA)
  • Algebra (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

The invention discloses a large-scale circuit fault diagnosis method and system of a Bayesian model based on Epanechnikov kernel density estimation of cloud theory, which belong to the field of large-scale analog circuit fault diagnosis, wherein fault information of a modularized integrated circuit after circuit tearing is collected, digital characteristics of clouds under each fault are generated by a reverse generator, a plurality of cloud drops under each fault mode are obtained by a forward cloud generator, membership values of each characteristic and corresponding fault state under a sample set are calculated by using a standard cloud model, cloud drops under each fault characteristic are determined according to a membership matrix, a fault sample set is expanded, and posterior probability density values under each fault category are calculated by using the Bayesian model based on Epanechnikov kernel density estimation, so that a fault diagnosis result is obtained. The method combines the cloud theory and the Bayesian kernel density, effectively solves the problems of single state information, less fault sample set and the like after the network is torn, integrates data information and priori knowledge, and greatly improves the fault accuracy.

Description

Cloud theory and nuclear density-based large-scale circuit fault diagnosis method and system
Technical Field
The invention belongs to the technical field of large-scale analog circuit fault diagnosis, and particularly relates to a cloud theory-based Bayesian model large-scale circuit fault diagnosis method and system for Epanechnikov kernel density estimation.
Background
Since the development of integrated circuits has been carried out in small, medium and large scale stages, and in the very large scale stage, industrial production has urgent requirements for large scale circuit testing, and the diagnosis of large scale circuit faults has not been fully mature in theory and method.
The accurate module-level fault positioning and clear identification detection result is a topic which needs to be solved by electronic engineering urgently, is a key step of practical application in theory, has huge number of devices contained in topology and continuity of each device parameter in an effective range for a large-scale analog circuit, so that the working condition of the circuit is more complicated, the state information of a large-scale circuit network is single after tearing, a simulation sample is difficult to obtain, the number of fault sample sets formed by tearing node voltages in each fault mode is small, accurate diagnosis of all faults is difficult to carry out, and the traditional fault analysis method is difficult to adapt to the current requirements.
Disclosure of Invention
Aiming at the defects or improvement demands of the prior art, the invention provides a Bayesian model large-scale circuit fault diagnosis method and system based on the Epanechnikov kernel density estimation of cloud theory, so as to improve the diagnosis accuracy and diagnosis efficiency of circuit diagnosis.
To achieve the above object, according to one aspect of the present invention, there is provided a method for diagnosing a large-scale circuit fault of a bayesian model of an Epanechnikov kernel density estimation based on cloud theory, including:
Establishing a large-scale circuit schematic diagram, selecting the voltage of a modularized tearing node as a diagnosis signal, carrying out coding classification on fault types according to different faults of components, and collecting the voltage of the tearing node under different faults to form an initial sample set;
Generating digital characteristic expectations, entropy and super entropy of cloud under each fault by using an inverse generator of a cloud model according to an initial sample set formed by voltages of tearing nodes under different faults;
according to the digital characteristic expectation, entropy and super entropy of the cloud under each fault, a plurality of cloud drops under each fault are obtained by utilizing a forward cloud generator so as to restore a plurality of measurable node voltage values;
Calculating each voltage value under the initial sample set and the membership value of the corresponding fault state by using a standard cloud model to form a membership matrix, and determining cloud drops under each fault voltage value according to the membership matrix of the fault so as to expand the initial sample set;
Dividing the expanded initial sample set into a plurality of folds for cross verification, and calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation;
And classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result.
In some alternative embodiments, the selecting the voltage of the modular tear node as the diagnostic signal comprises:
and carrying out modularized tearing on the large-scale circuit, tearing the limited accessible node into modularization on the premise of the limited accessible node, and selecting the voltage of the torn node in each fault mode as a diagnosis signal to obtain a fault characteristic vector.
In some alternative embodiments, the initial sample set of voltages of tearing nodes at different faults uses an inverse generator of a cloud model to generate digital feature expectations, entropy and super entropy of the cloud at each fault, comprising:
Note the initial sample set x ij=[xi1,xi2,…,xin of circuit faults (i=1, 2,..n; j=1, 2,..m.), a j-th class fault that characterizes the i-th fault, n is the dimension of the feature vector in the fault mode, m is the total number of fault modes, then Obtaining the mean value of the initial sample set of the j-th type faultFrom the following components A first order central moment K of the initial sample set of the j-th class of faults is obtained, from the following componentsObtaining the variance S of an initial sample set of the j-th type fault;
since the certainty is unknown, the digital characteristic expectations E x, the entropy E n and the super entropy H e of the cloud model obtained by the inverse cloud generator CG -1 without certainty are as follows:
in some alternative embodiments, the obtaining a plurality of cloud droplets per fault using a forward cloud generator includes:
On the basis of the cloud { E x,En,He } of each type of faults, setting up the number Num of cloud drops, utilizing a forward cloud model CG (E x,En,He, num), firstly, for the j type of faults, generating a model with E n as an expectation, A normal random number E' ni that is the variance;
using the generated random number, with E x as the expectation, Generating a normal random number x k for variance and calculating the x k membershipWherein x k of u k with certainty represents a cloud droplet in the domain;
Repeating the steps until the set cloud drop number Num is generated, and obtaining a cloud sample set drop (x k,uk) of the j-th type fault.
In some alternative embodiments, the calculating the membership value of each voltage value and the corresponding fault state in the initial sample set by using the standard cloud model forms a membership matrix, including:
Describing membership degrees of fault feature vectors corresponding to fault types through a cloud model, substituting standard cloud parameters corresponding to the j-th type fault of the i-th feature vector of an initial fault sample x ij into the standard cloud model formed under each fault Obtaining the corresponding membership of x ij, wherein u ij is the membership of the j-th fault of the i-th feature vector of the initial fault sample, E nij' is the normal random number generated by the entropy E nij and the super entropy H eij 2 of the j-th fault of the i-th feature vector, and E xij is the digital feature expectation of the j-th fault of the i-th feature vector;
After the membership degree of various faults of the initial fault sample under all feature vectors is obtained, membership degree moment is obtained: m represents the total number of failure modes;
the fuzzy vector is formed by fuzzy comprehensive evaluation uncertainty index Obtaining the membership degree of uncertainty, wherein n is a fuzzy comprehensive evaluation vector (u i1,ui2,…,uin), i=1, 2, & n, beta is the maximum membership degree in (u i1,ui2,…,uin), and gamma is the second maximum membership degree in (u i1,ui2,…,uin);
calculating the uncertain membership degree of each fault under the membership degree to form an uncertain membership degree fuzzy characteristic vector And (3) representing the uncertain membership degree of the ith eigenvector to obtain a membership degree matrix of m rows and n+1 columns:
In some alternative embodiments, the dividing the expanded initial sample set into a plurality of fold cross-validations, calculating posterior probability density values for each fault class in the test set using a bayesian model of an Epanechnikov kernel density estimation, includes:
Cross-verifying the expanded initial sample set into a training set and a testing set;
In the Bayesian model, a likelihood function of the Epanechnikov kernel density estimation can be obtained through a plurality of independent samples distributed in the test set, posterior distribution of the likelihood function of the Epanechnikov kernel density estimation is obtained according to the Bayesian theorem, and then posterior probability density values under each fault category are obtained.
In some alternative embodiments, the classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain the diagnosis result includes:
Classifying fault data by using an Epanechnikov kernel density estimation Bayesian model, regarding voltage data points of faults to be detected as a vector, and estimating probability density functions of each fault type by using an Epanechnikov kernel density estimation method;
And projecting the voltage data points of the faults to be detected onto the probability density function of each fault class, calculating the probability density value under each class, and classifying the voltage data points of the faults to be detected into the class with the highest probability density value.
In some optional embodiments, the projecting the voltage data points of the fault to be detected onto the probability density function of each fault class, calculating the probability density value under each class, and classifying the voltage data points of the fault to be detected into the class with the highest probability density value includes:
Taking the fault voltage data point to be detected as a vector, wherein the vector comprises a plurality of features, subtracting the posterior mean value of the class from the features in the vector, dividing each element obtained by the posterior square difference of the class, substituting the obtained element into an Epanechnikov kernel function as an independent variable, and taking the average value of the obtained new vector to obtain the probability density value of the fault sample to be detected under the class;
calculating the reciprocal of the square root of the determinant of the posterior variance of the class, multiplying the reciprocal by the obtained probability density value to normalize the probability value;
and calculating the posterior probability density value of the fault voltage data point to be detected under each fault category according to the Bayesian theorem, and classifying the fault voltage data point to be detected into the category with the highest probability density value.
According to another aspect of the present invention, there is provided a system for diagnosing a large-scale circuit fault of a bayesian model of an Epanechnikov kernel density estimation based on a cloud theory, including:
The circuit acquisition module to be diagnosed is used for establishing a large-scale circuit schematic diagram, selecting the voltage of the modularized tearing node as a diagnosis signal, carrying out coding classification on fault types according to different faults of components and parts, and collecting the voltage of the tearing node under different faults to form an initial sample set;
The cloud theory calculation module is used for generating digital characteristic expectations, entropies and super entropies of clouds under each fault by using an initial sample set formed by voltages of tearing nodes under different faults and using a reverse generator of a cloud model, and obtaining a plurality of cloud drops under each fault by using the forward cloud generator according to the digital characteristic expectations, the entropies and the super entropies of the clouds under each fault so as to restore a plurality of measurable node voltage values;
The membership degree expansion sample module is used for calculating each voltage value under the initial sample set and the membership degree value corresponding to the fault state to form a membership degree matrix by utilizing the standard cloud model, and determining cloud drops under each fault voltage value according to the membership degree matrix of the fault so as to expand the initial sample set;
The fault diagnosis module is used for dividing the expanded initial sample set into a plurality of fold cross-validation, calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation, and classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result.
According to another aspect of the present invention there is provided a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of any of the methods described above.
In general, the above technical solutions conceived by the present invention, compared with the prior art, enable the following beneficial effects to be obtained:
The Bayesian model of the Epanechnikov kernel density estimation based on the cloud theory is used for large-scale circuit fault diagnosis, cloud drops generated by combining a forward cloud model and a reverse cloud model in the cloud theory are used for expanding a fault sample set, the generation of a sample based on the cloud fuses the randomness and ambiguity uncertainty theory of the cloud model, characteristic data deviating from normal distribution is filtered, the model has certain anti-noise capability, the Bayesian model of the Epanechnikov kernel function does not need to make any assumption on the distribution of target variables, and the model based on the circuit fault diagnosis obtained on the basis can be better suitable for different data sets and problems, so that the accuracy and the diagnosis efficiency of circuit fault diagnosis are effectively improved.
Drawings
FIG. 1 is a schematic flow chart of a circuit fault diagnosis method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a video amplifying circuit according to an embodiment of the present invention;
FIG. 3 is a tearing diagram of a video amplifying circuit according to an embodiment of the present invention;
FIG. 4 is a fault cloud sample generation process provided by an embodiment of the present invention;
FIG. 5 is a graph comparison of cloud models under normal and fault conditions provided by an embodiment of the present invention;
FIG. 6 is a confusion matrix of fault diagnosis results provided by an embodiment of the present invention;
Fig. 7 is a schematic structural diagram of a circuit fault diagnosis device according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
Referring to fig. 1, a schematic flow chart of a method provided by an embodiment of the present invention is shown, and the method shown in fig. 1 includes the following steps:
(1) Establishing a simulation model of the video amplifying circuit, wherein a simulation circuit diagram is shown in fig. 2, carrying out modularized division on the video amplifying circuit under the condition that node tearing is met, and collecting voltage state information of torn nodes to obtain fault characteristic vectors, wherein tolerance of resistance in the circuit is 5%, considering that a semiconductor device and a passive device are easier to generate faults, and according to practical experience, the circuit is provided with 9 fault modes including normal states, as shown in table 1;
TABLE 1 failure states and categories
Fault labels Fault state Fault coding
F0 Normal state 0
F1 V 1 open circuit 1
F2 R 1 short circuit 2
F3 Q 4 BE extremely short circuit 3
F4 Q 2 CE extremely short circuit 4
F5 R 9 short circuit 5
F6 V 2 open circuit 6
F7 V 3 open circuit 7
F8 Q 10 CE extremely short circuit 8
In step (1), the large-scale circuit is subjected to modularized tearing, and one network needs to meet the following conditions that 1, the torn node is an accessible node and comprises a public node, 2, no topological relation and no coupling between parameters exist among all sub-networks after tearing, 3, all sub-networks are minimum sub-networks relative to the torn node, all sub-networks are only used for enabling the node to be torn, namely, the sub-network scale is minimum, the fewer the number of faults is combined, and the higher the diagnosis rate is.
(2) Sending an initial sample set formed by the voltage of the tearing node under each fault mode to a reverse cloud generator to obtain cloud { E x,En,He };
(3) A plurality of cloud drops under each fault mode are obtained by utilizing a forward cloud generator algorithm to restore a plurality of measurable node voltage values;
(4) Calculating membership values of each voltage characteristic and corresponding fault state under the initial sample set by using a standard cloud model, and determining cloud drops under each fault voltage characteristic according to a membership matrix of an original fault to expand the initial sample set;
(5) Dividing the expanded initial sample set into a plurality of folds for cross verification, and calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation;
(6) And classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result.
In this embodiment, the above step (2) may be implemented by:
The cloud sample generation needs a reverse cloud and forward cloud model algorithm at the same time, and for each circuit fault, the cloud can be used for describing the uncertain conversion between the qualitative concept and the fixed value of the fault mode so as to reflect the uncertainty of things or human knowledge concepts in the natural world, namely ambiguity and randomness, which not only give explanation from random theory and fuzzy set theory, but also reflect the relevance between the ambiguity and randomness, so that the mapping between quantification and qualitative is formed, as shown in fig. 4, and the generation process of the fault cloud sample is realized. Elements in a sample set in a fault mode can be represented by cloud drops, the sampled voltage is converted from a fixed value to a qualitative concept by an inverse cloud generator, and the qualitative concept is represented by a digital feature { E x,En,He }, wherein the digital feature is expected to be E x, entropy E n and super-entropy H e.
The algorithm of the one-dimensional reverse cloud generator is realized as follows:
Inputting a quantitative value of an initial fault sample set;
outputting an expected value E x, entropy E n and super entropy H e of the qualitative concept A represented by cloud drops;
note the initial sample set x ij=[xi1,xi2,…,xin of circuit faults (i=1, 2,., n; j=1, 2, m), calculating a sample mean Absolute central moment of first-order sampleSample variance
In the domain space, a cloud is the point that best represents the qualitative concept, which is expected to be the central value in the domain space, which is expected
Entropy is determined by randomness and ambiguity of qualitative concepts, and represents a measurable granularity of a qualitative concept, E n is a measure of randomness of the qualitative concept, reflects the discrete degree of the cloud, reflects the margin of the qualitative concept, reflects the range of values of the cloud which can be accepted by the qualitative concept in the domain space, is a measure of ambiguity of the qualitative concept, and generally, the larger the entropy is, the larger the range of values of the cloud which can be accepted by the qualitative concept is, the more the qualitative concept is ambiguous, and reflects the relativity between the randomness and the ambiguity, and the entropy can be obtained by the average value of samples
Super entropy is a measure of uncertainty of entropy, revealing the cohesiveness of uncertainty of all points of language values in the domain space and the association of ambiguity and randomness, indirectly reflecting the thickness of cloud, and is obtained by sample variance and super entropy
In this embodiment, the above step (3) may be implemented by:
The forward cloud generator is a mapping from qualitative concepts to quantitative representations thereof, and generates cloud droplets according to the digital features { E x,En,He } of the cloud obtained above, each cloud droplet being a specific implementation of the concept.
The algorithm of the one-dimensional forward cloud generator is realized as follows:
Three digital characteristic values E x,En,He representing the qualitative concept A and a cloud drop number Num are input;
Outputting quantitative values of Num cloud drops, wherein each cloud drop represents the certainty degree of the concept A;
the generation is expected to take E n as a matter of course, A normal random number E' ni that is the variance;
using the generated random number, with E x as the expectation, Generating a normal random number x k for variance and calculating the membership of x k Wherein x k of u k with certainty represents a cloud droplet in the domain;
Repeating the steps until a set number of cloud drops Num is generated, and obtaining a cloud sample set drop (x k,uk).
In this example, the number of droplets may be num=10000.
As shown in fig. 5, the cloud model and the sample distribution situation with 10000 cloud droplets are generated by the fault labels F0, F3, F5 and F7 respectively.
From comparison, the expectation, the certainty (membership), and the randomness (discrete degree) are seen from the high to low limit of the cloud model analysis. The range of the cloud drops accepted by each fault in the domain space and the aggregation degree of the cloud drops are greatly different, and the normal state and the fault state are greatly different.
In this embodiment, the above step (4) may be implemented by:
in order to avoid excessively long feature vectors when a fault sample set is expanded, the embodiment of the invention selects a cloud membership matrix under the feature nodes to obtain cloud droplets with larger occurrence probability and certainty as the feature vectors so as to facilitate the next step of fault mode identification.
In order to generate a cloud sample set without increasing data dimension, describing membership degrees of each fault feature vector corresponding to each fault type through a cloud model, calculating the membership degree of the j-th fault of the i-th feature vector by utilizing a standard cloud model formed under each fault, and substituting standard cloud parameters corresponding to the j-th fault of the i-th feature of an initial fault sample x ij into the cloud modelAnd obtaining a corresponding membership degree, wherein u ij is the membership degree of the j-th class fault of the i-th feature vector of the initial fault sample, and E nij' is a normal random number generated by the entropy E nij and the super entropy H eij 2 of the j-th class fault of the i-th feature vector.
Meanwhile, the uncertainty causes of data errors and external interference in the fault diagnosis process are considered, a fuzzy vector is formed by fuzzy comprehensive evaluation uncertainty indexes, and the membership calculation formula of the uncertainty is as follows: n is the fuzzy comprehensive evaluation vector (u i1,ui2,…,uin), i=1, 2,..n the number of membership degrees, β is the maximum membership degree of (u i1,ui2,…,uin), and γ is the second largest membership degree of (u i1,ui2,…,uin). Calculating the uncertain membership degree of each fault under the membership degree to form an uncertain membership degree fuzzy characteristic vector as shown in the formula In the middle ofAn uncertain membership representing the ith feature. By integrating the above formula, the membership matrix of m rows and n+1 columns can be obtained again.
Taking the influence of external interference into consideration, collecting 200 cloud drops under each membership degree u+/-0.01, and finally obtaining 11 multiplied by 9 multiplied by 200 fault samples with 11 fault characteristic attributes and 9 fault categories.
In this embodiment, the above step (5) may be implemented by:
the sample set is cross-validated into a training set and a test set.
The epannechnikov kernel density estimation is a non-parametric probability density estimation method that can be used to estimate arbitrary distributions, while the bayesian model is a bayesian-based probability model that can be used to estimate posterior distributions of unknown parameters. In the bayesian model, unknown parameters are regarded as random variables, posterior distribution of the parameters is updated through observation data, and n independent samples X 1,X2,…Xn with the same distribution are obtained, so that the likelihood function expression of the Epanechnikov kernel density estimation can be obtained as follows:
Wherein h is a bandwidth parameter, determines the smoothness of a kernel function, and obtains posterior distribution of pi (f|X 1,…,Xn)∞L(f|X1,…,Xn) pi (f) according to Bayesian theorem;
In this bayesian model, it is assumed that the variable y is generated from an unknown probability distribution, and then an prior probability distribution is constructed using an Epanechnikov kernel function, which represents the density distribution of y throughout the input space. And combining the observed data x with the prior probability distribution through Bayesian rules to obtain posterior probability distribution. Finally, the posterior probability distribution can be used to predict the target variable value for the new data point. Unlike traditional parametric regression models, the bayesian model based on the Epanechnikov kernel function does not require any assumptions on the distribution of the target variables, and thus can better accommodate different datasets and problems.
In this embodiment, step (6) may be implemented by:
The method comprises classifying fault data by using a Bayesian model of Epanechnikov kernel density estimation, regarding each fault voltage data point as a vector, and estimating a probability density function of each fault category by using an Epanechnikov kernel density estimation method. The data points to be classified are then projected onto the probability density function of each fault class, and their probability density values under each class are calculated. Finally, the data points are classified into the category with the highest probability density value. Specifically, the present invention relates to a method for manufacturing a semiconductor device.
Regarding the fault data point to be tested as a vector, the fault data point to be tested comprises a plurality of characteristics, subtracting the characteristics in the vector from the posterior mean value of the class, dividing each element obtained by the posterior variance of the class as an independent variable, substituting the independent variable into an Epanechnikov kernel function, obtaining a new vector average value, namely a probability density value of a sample to be tested under the class, calculating the reciprocal of the square root of the determinant of the posterior variance of the class, and multiplying the probability density value obtained in the last step to normalize the probability value. And calculating the posterior probability density value of the data point to be classified under each fault category according to the Bayesian theorem, and classifying the data point to be classified into the category with the highest probability density value.
The confusion matrix of the fault classification results is shown in fig. 6, and the diagnosis accuracy is 99.4%.
According to the large-scale circuit fault diagnosis method of the Bayesian model based on the Epanechnikov kernel density estimation of the cloud theory, disclosed by the invention, a fault sample set is expanded by utilizing cloud drops generated by combining a forward cloud model and a reverse cloud model in the cloud theory, the situation that the number of samples after network tearing is small or a simulation sample is difficult to obtain under a specific condition is effectively solved, the Bayesian model based on the Epanechnikov kernel density estimation is used for carrying out pattern recognition, the probability density function of Gaussian distribution to estimate class conditions is solved, and the diagnosis degree of fault diagnosis is improved due to the fact that the probability density function is too fit caused by the small number.
As shown in fig. 7, in another embodiment of the present invention, there is further provided a system for diagnosing a large-scale circuit fault of a bayesian model of an Epanechnikov kernel density estimation based on a cloud theory, including:
The circuit acquisition module to be diagnosed is used for establishing a large-scale circuit schematic diagram, selecting the voltage of the modularized tearing node as a diagnosis signal, carrying out coding classification on fault types according to different faults of components and parts, and collecting the voltage of the tearing node under different faults to form an initial sample set;
The cloud theory calculation module is used for generating digital characteristic expectations, entropies and super entropies of clouds under each fault by using an initial sample set formed by voltages of tearing nodes under different faults and using a reverse generator of a cloud model, and obtaining a plurality of cloud drops under each fault by using the forward cloud generator according to the digital characteristic expectations, the entropies and the super entropies of the clouds under each fault so as to restore a plurality of measurable node voltage values;
The membership degree expansion sample module is used for calculating each voltage value under the initial sample set and the membership degree value corresponding to the fault state to form a membership degree matrix by utilizing the standard cloud model, and determining cloud drops under each fault voltage value according to the membership degree matrix of the fault so as to expand the initial sample set;
The fault diagnosis module is used for dividing the expanded initial sample set into a plurality of fold cross-validation, calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation, and classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result.
The specific implementation of each module may refer to the description of the above method embodiment, and this embodiment will not be repeated.
It should be noted that each step/component described in the present application may be split into more steps/components, or two or more steps/components or part of operations of the steps/components may be combined into new steps/components, according to the implementation needs, to achieve the object of the present application.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.

Claims (7)

1. A cloud theory-based method for diagnosing a large-scale circuit fault of a bayesian model of an Epanechnikov kernel density estimation, comprising the steps of:
Establishing a large-scale circuit schematic diagram, selecting the voltage of a modularized tearing node as a diagnosis signal, carrying out coding classification on fault types according to different faults of components, and collecting the voltage of the tearing node under different faults to form an initial sample set;
Generating digital characteristic expectations, entropy and super entropy of cloud under each fault by using an inverse generator of a cloud model according to an initial sample set formed by voltages of tearing nodes under different faults;
Obtaining a plurality of cloud drops under each fault by utilizing digital characteristic expectation, entropy and super entropy of the cloud under each fault by utilizing a forward cloud generator to restore a plurality of measurable node voltage values, describing randomness and ambiguity of the fault mode and the relevance between the two by utilizing the cloud for each circuit fault mode, wherein elements in a fault sample set under the fault mode can be regarded as formed cloud drops;
Calculating each voltage value under the initial sample set and the membership value of the corresponding fault state by using a standard cloud model to form a membership matrix, and determining cloud drops under each fault voltage value according to the membership matrix of the fault so as to expand the initial sample set;
Dividing the expanded initial sample set into a plurality of folds for cross verification, and calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation;
classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result;
dividing the expanded initial sample set into a plurality of folds for cross-validation, calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation, wherein the method comprises the following steps:
Cross-verifying the expanded initial sample set into a training set and a testing set;
In a Bayesian model, a likelihood function of Epanechnikov kernel density estimation can be obtained through a plurality of independent samples distributed in the test set, posterior distribution of the likelihood function of the Epanechnikov kernel density estimation is obtained according to the Bayesian theorem, and then posterior probability density values under each fault category are obtained;
The step of classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result comprises the following steps:
Classifying fault data by using an Epanechnikov kernel density estimation Bayesian model, regarding voltage data points of faults to be detected as a vector, and estimating probability density functions of each fault type by using an Epanechnikov kernel density estimation method;
Projecting voltage data points of faults to be detected onto probability density functions of each fault class, calculating probability density values under each class, and classifying the voltage data points of the faults to be detected into the class with the highest probability density values;
Projecting the voltage data points of the faults to be detected onto probability density functions of each fault class, calculating probability density values under each class, and classifying the voltage data points of the faults to be detected into the class with the highest probability density values, wherein the method comprises the following steps:
Taking the fault voltage data point to be detected as a vector, wherein the vector comprises a plurality of features, subtracting the posterior mean value of the class from the features in the vector, dividing each element obtained by the posterior square difference of the class, and substituting the posterior square difference into an Epanechnikov kernel function as an independent variable to obtain a new vector average value, and obtaining the probability density value of the fault sample to be detected under the class;
calculating the reciprocal of the square root of the determinant of the posterior variance of the class, multiplying the reciprocal by the obtained probability density value to normalize the probability value;
and calculating the posterior probability density value of the fault voltage data point to be detected under each fault category according to the Bayesian theorem, and classifying the fault voltage data point to be detected into the category with the highest probability density value.
2. The fault diagnosis method according to claim 1, wherein selecting the voltage of the modular tear node as the diagnosis signal comprises:
and carrying out modularized tearing on the large-scale circuit, tearing the limited accessible node into modularization on the premise of the limited accessible node, and selecting the voltage of the torn node in each fault mode as a diagnosis signal to obtain a fault characteristic vector.
3. The fault diagnosis method according to claim 1 or 2, wherein the initial sample set of voltages of tearing nodes at different faults uses a reverse generator of a cloud model to generate digital characteristic expectations, entropy and super entropy of the cloud at each fault, comprising:
Note that the initial sample set x ij=[xi1,xi2,…,xin of circuit faults (i=1, 2,..n; j=1, 2,..m) represents the j-th fault of the i-th fault feature, n is the dimension of the feature vector in the fault mode, and m is the total number of fault modes, and then the fault is represented by the following formula Obtaining the mean value of the initial sample set of the j-th type faultFrom the following components A first order central moment K of the initial sample set of the j-th class of faults is obtained, from the following componentsObtaining the variance S of an initial sample set of the j-th type fault;
since the certainty is unknown, the digital characteristic expectations E x, the entropy E n and the super entropy H e of the cloud model obtained by the inverse cloud generator CG -1 without certainty are as follows:
4. a fault diagnosis method according to claim 3, wherein said obtaining a plurality of cloud droplets for each fault using a forward cloud generator comprises:
On the basis of the cloud { E x,En,He } of each type of faults, setting up the number Num of cloud drops, utilizing a forward cloud model CG (E x,En,He, num), firstly, for the j type of faults, generating a model with E n as an expectation, A normal random number E' ni that is the variance;
Using the generated random number, generating a random number expected by E x, Generating a normal random number x k for variance and calculating the x k membershipWherein x k of u k with certainty represents a cloud droplet in the domain;
Repeating the steps until the set cloud drop number Num is generated, and obtaining a cloud sample set drop (x k,uk) of the j-th type fault.
5. The fault diagnosis method according to claim 4, wherein calculating the membership matrix of each voltage value and the membership value of the corresponding fault state in the initial sample set using the standard cloud model comprises:
Describing membership degrees of fault feature vectors corresponding to fault types through a cloud model, substituting standard cloud parameters corresponding to the j-th type fault of the i-th feature vector of an initial fault sample x ij into the standard cloud model formed under each fault Obtaining the corresponding membership of x ij, wherein u ij is the membership of the j-th fault of the i-th feature vector of the initial fault sample, E nij' is the normal random number generated by the entropy E nij and the super entropy H eij 2 of the j-th fault of the i-th feature vector, and E xij is the digital feature expectation of the j-th fault of the i-th feature vector;
After the membership degree of various faults of the initial fault sample under all feature vectors is obtained, membership degree moment is obtained: m represents the total number of failure modes;
the fuzzy vector is formed by fuzzy comprehensive evaluation uncertainty index Obtaining the membership degree of uncertainty, wherein n is a fuzzy comprehensive evaluation vector (u i1,ui2,…,uin), i=1, 2, & n, beta is the maximum membership degree in (u i1,ui2,…,uin), and gamma is the second maximum membership degree in (u i1,ui2,…,uin);
calculating the uncertain membership degree of each fault under the membership degree to form an uncertain membership degree fuzzy characteristic vector And (3) representing the uncertain membership degree of the ith eigenvector to obtain a membership degree matrix of m rows and n+1 columns:
6. A cloud theory-based, bayesian model based system for diagnosing faults in a large scale of circuit for estimating the density of an Epanechnikov kernel, comprising:
The circuit acquisition module to be diagnosed is used for establishing a large-scale circuit schematic diagram, selecting the voltage of the modularized tearing node as a diagnosis signal, carrying out coding classification on fault types according to different faults of components and parts, and collecting the voltage of the tearing node under different faults to form an initial sample set;
The cloud theory calculation module is used for generating digital characteristic expectations, entropies and super entropies of clouds under each fault by using an initial sample set formed by voltages of tearing nodes under different faults and using a reverse generator of a cloud model, and obtaining a plurality of cloud drops under each fault by using the forward cloud generator according to the digital characteristic expectations, the entropies and the super entropies of the clouds under each fault so as to restore a plurality of measurable node voltage values;
The membership degree expansion sample module is used for calculating each voltage value under the initial sample set and the membership degree value corresponding to the fault state to form a membership degree matrix by utilizing the standard cloud model, and determining cloud drops under each fault voltage value according to the membership degree matrix of the fault so as to expand the initial sample set;
The fault diagnosis module is used for dividing the expanded initial sample set into a plurality of fold cross-validation, calculating posterior probability density values under each fault category in the test set by using a Bayesian model of Epanechnikov kernel density estimation, and classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result;
The method for calculating the posterior probability density value under each fault category in the test set by using the Bayesian model of the Epanechnikov kernel density estimation comprises the following steps:
Cross-verifying the expanded initial sample set into a training set and a testing set;
In a Bayesian model, a likelihood function of Epanechnikov kernel density estimation can be obtained through a plurality of independent samples distributed in the test set, posterior distribution of the likelihood function of the Epanechnikov kernel density estimation is obtained according to the Bayesian theorem, and then posterior probability density values under each fault category are obtained;
The step of classifying the newly acquired test sample into the category with the highest posterior probability density value to obtain a diagnosis result comprises the following steps:
Classifying fault data by using an Epanechnikov kernel density estimation Bayesian model, regarding voltage data points of faults to be detected as a vector, and estimating probability density functions of each fault type by using an Epanechnikov kernel density estimation method;
Projecting voltage data points of faults to be detected onto probability density functions of each fault class, calculating probability density values under each class, and classifying the voltage data points of the faults to be detected into the class with the highest probability density values;
Projecting the voltage data points of the faults to be detected onto probability density functions of each fault class, calculating probability density values under each class, and classifying the voltage data points of the faults to be detected into the class with the highest probability density values, wherein the method comprises the following steps:
Taking the fault voltage data point to be detected as a vector, wherein the vector comprises a plurality of features, subtracting the posterior mean value of the class from the features in the vector, dividing each element obtained by the posterior square difference of the class, and substituting the posterior square difference into an Epanechnikov kernel function as an independent variable to obtain a new vector average value, and obtaining the probability density value of the fault sample to be detected under the class;
calculating the reciprocal of the square root of the determinant of the posterior variance of the class, multiplying the reciprocal by the obtained probability density value to normalize the probability value;
and calculating the posterior probability density value of the fault voltage data point to be detected under each fault category according to the Bayesian theorem, and classifying the fault voltage data point to be detected into the category with the highest probability density value.
7. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 5.
CN202310611976.1A 2023-05-26 2023-05-26 Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density Active CN116720103B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310611976.1A CN116720103B (en) 2023-05-26 2023-05-26 Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310611976.1A CN116720103B (en) 2023-05-26 2023-05-26 Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density

Publications (2)

Publication Number Publication Date
CN116720103A CN116720103A (en) 2023-09-08
CN116720103B true CN116720103B (en) 2025-08-26

Family

ID=87872509

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310611976.1A Active CN116720103B (en) 2023-05-26 2023-05-26 Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density

Country Status (1)

Country Link
CN (1) CN116720103B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118485451B (en) * 2024-05-24 2024-11-22 海南天地神农科技有限公司 Agricultural product traceability method and system based on source data verification mechanism

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2464364A1 (en) * 2001-10-17 2003-04-24 Commonwealth Scientific And Industrial Research Organisation Method and apparatus for identifying diagnostic components of a system
CN107505133A (en) * 2017-08-10 2017-12-22 滁州学院 The probability intelligent diagnosing method of rolling bearing fault based on adaptive M RVM

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7539597B2 (en) * 2001-04-10 2009-05-26 Smartsignal Corporation Diagnostic systems and methods for predictive condition monitoring
CN106597260B (en) * 2016-12-29 2020-04-03 合肥工业大学 Analog circuit fault diagnosis method based on continuous wavelet analysis and ELM network
CN108051660A (en) * 2017-10-31 2018-05-18 华北电力大学(保定) A kind of transformer fault combined diagnosis method for establishing model and diagnostic method
CN109345032B (en) * 2018-10-30 2021-09-21 南京航空航天大学 Particle filter multi-crack-propagation prediction method based on dynamic crack number
US10901786B2 (en) * 2019-03-09 2021-01-26 International Business Machines Corporation Tool for identifying sources of operating system jitter from power consumption measurements
CN114036857B (en) * 2021-11-26 2024-06-14 天津大学 Flexible riser risk evaluation method based on quantitative analysis and probabilistic reasoning

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2464364A1 (en) * 2001-10-17 2003-04-24 Commonwealth Scientific And Industrial Research Organisation Method and apparatus for identifying diagnostic components of a system
CN107505133A (en) * 2017-08-10 2017-12-22 滁州学院 The probability intelligent diagnosing method of rolling bearing fault based on adaptive M RVM

Also Published As

Publication number Publication date
CN116720103A (en) 2023-09-08

Similar Documents

Publication Publication Date Title
Nguyen et al. Spatial-temporal recurrent graph neural networks for fault diagnostics in power distribution systems
CN109145516B (en) A fault identification method for analog circuits based on an improved extreme learning machine
CN113484738A (en) Circuit fault diagnosis method based on multi-feature information fusion
CN114048819B (en) Power distribution network topology identification method based on attention mechanism and convolutional neural network
CN111523778A (en) Power grid operation safety assessment method based on particle swarm algorithm and gradient lifting tree
CN113111618B (en) A Fault Diagnosis Method for Analog Circuits Based on Improved Empirical Wavelet Transform
CN108919059A (en) A kind of electric network failure diagnosis method, apparatus, equipment and readable storage medium storing program for executing
CN112561303B (en) Power system dynamic analysis method based on ensemble learning and grid topology changes
CN113435575B (en) Transient stability evaluation method of gate graph neural network based on unbalanced data
CN116720103B (en) Large-scale circuit fault diagnosis method and system based on cloud theory and kernel density
CN108595884A (en) Power system transient stability appraisal procedure and device
CN109374997A (en) Detection and evaluation method of power quality disturbance in hybrid power system based on VMD initialization S transform
CN112183606A (en) Power system fault identification and classification method and system based on C4.5 algorithm
CN111273126A (en) Power distribution network topology rapid sensing method
CN117493841A (en) Power distribution network line fault prediction method based on BiGRU-Attention model
CN117195117A (en) Fault scene reduction method, device, equipment and storage medium
Sahoo et al. Model order estimation for low-frequency oscillations in power systems by an advanced K-mean clustering approach
CN115128410A (en) A TPA-LSTM-based Partial Discharge Fault Mode Recognition Method for DC Cables
Banerjee et al. Data-driven static load model parameter estimation with confidence factor
CN115841179A (en) Power system situation sensing method based on graph digital twins
CN118095050A (en) Power system transient stability assessment method considering sample imbalance based on artificial intelligence
Liang et al. Advanced voltage abnormality detection in real-vehicle battery systems using self-organizing map neural networks and adaptive threshold
Nayak et al. Analysis of wind characteristics using ARMA & Weibull distribution
Zhu et al. Research on circuit fault diagnosis method based on multi-feature information fusion
CN111368933A (en) Power distribution network transient process fault classification method and system based on Softmax regression

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant