Next Article in Journal
Small Defect Detection Based on Local Structure Similarity for Magnetic Tile Surface
Next Article in Special Issue
Research on Spectrum Prediction Technology Based on B-LTF
Previous Article in Journal
Influence Analysis of SiC MOSFET’s Parasitic Capacitance on DAB Converter Output
Previous Article in Special Issue
Soft Fault Diagnosis of Analog Circuit Based on EEMD and Improved MF-DFA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Analog Circuit Soft Fault Diagnosis Method Based on Mathematical Morphology Fractal Dimension

1
School of Measurement-Control Technology and Communications Engineering, Harbin University of Science and Technology, Harbin 150080, China
2
Heilongjiang Network Space Research Center, Harbin 150090, China
3
Harbin Meteorological Bureau, Harbin 150028, China
4
Harbin Vocational College of Science and Technology, Harbin 150399, China
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(1), 184; https://doi.org/10.3390/electronics12010184
Submission received: 21 November 2022 / Revised: 27 December 2022 / Accepted: 28 December 2022 / Published: 30 December 2022
(This article belongs to the Special Issue Deep Learning in Image Processing and Pattern Recognition)

Abstract

:
It is difficult for traditional circuit-fault feature-extraction methods to accurately distinguish between nonlinear analog-circuit faults and analog-circuit faults with high fault rates and high diagnostic costs. To solve this problem, this paper proposes a method of mathematical morphology fractal dimension (VMD-MMFD) based on variational mode decomposition for soft-fault feature extraction in analog circuits. First, the signal is decomposed into variational modes to suppress the influence of environmental noise, and multiple high-dimensional eigenmode functions with different center frequencies are obtained. The fractal dimension of the signal feature information component IMF is calculated, and then, KPCA (Kernel Principal Component Analysis) is used to remove the overlapping and redundant parts of the data. The fault set obtained is used as the basis for judging the working state and the fault type of the circuit. The experimental results of the simulation circuits show that this method can be effectively used for circuit-fault diagnosis.

1. Introduction

In daily life and production, circuits can be divided into two types: analog and digital. Although analog circuits generally account for less than 20% of circuit components, 80% of faults are caused by analog circuits [1]. The fault types of analog circuits can be divided into soft faults and hard faults. A hard fault refers to catastrophic faults such as open circuits and short circuits in electronic circuits, which are easy to identify. A soft fault indicates that the value of the components has changed, and the deviation in value exceeds the allowable range [2]. When a soft fault problem occurs, the circuit can still work normally, but if the component is not replaced in time, the soft fault will be upgraded to a hard fault, which will cause significant damage to the entire circuit [3] and even endanger people’s lives and safety.
Many scholars have proposed many feature extraction methods for analog-circuit recognition. Common fault feature extraction methods of analog circuits include principal component analysis [4], factor analysis [5,6], and other linear extraction methods. These methods are relatively effective for linear circuits, but most circuits in daily life have nonlinear features. The above methods cannot reflect the non-stationary characteristics of the signal, resulting in the low separability of the extracted fault features, so there is a large classification error in fault-pattern recognition. Therefore, early soft fault research mainly introduced fuzzy algorithms, wavelet theory, and other means to determine the actual working conditions [7]. Although this algorithm improves the effect of fault diagnosis, some algorithms are seriously affected by the circuit state when analyzing fault characteristics, which makes the performance unstable. In order to solve this problem, in recent years, some scholars have combined mathematical morphology with wavelets and proposed a new nonlinear wavelet morphological wavelet. Reference [8] applied morphological wavelets to detect power disturbance and successfully identified bad conditions in the power transmission process. For example, Zhuang Ning [9] and others combined fractal dimension with EMD (empirical mode decomposition) to extract features of ECG signals to identify different emotional performances. Zheng Zhi [10] and others combined the LMD (local mean decomposition) decomposition method with generalized fractal dimensions to identify gear faults. However, actual measured signals are often accompanied by abnormal events, such as the noise impact and intermittent signals, which lead to modal aliasing and greatly affect the subsequent fault-diagnosis rate. At the same time, a common method to calculate the fractal dimension is the covering method (the box-counting method). This method has inevitable shortcomings. The method itself uses the regular grid division method, so the estimation of the fractal dimension is very unstable in some cases [11].
In view of the above problems, this paper proposes a method to calculate the fractal dimensions (VMD-MMFD-KPCA) of mathematical morphology based on variational mode decomposition (VMD). Using the variational mode decomposition (VMD) of the signal, the interference information in the signal is eliminated to the maximum extent possible to solve the problem of mode aliasing of traditional EMD decomposition. When using morphological fractal dimensions to calculate the dimensions of signals, the instability of traditional box-dimension estimation is avoided, including effectively distinguishing different fault types. On this basis, the KPCA dimension-reduction method is introduced to reduce the dimensions of high-dimensional data calculated by dimensions, eliminate the redundancy and duplication in the sample, and provide a data basis for subsequent fault diagnosis. Compared with traditional fault-diagnosis methods, it shows better performance of feature extraction and diagnosis. Finally, the effectiveness of this method is demonstrated by a simulation circuit.

2. Variational Modal Decomposition

VMD decomposition is a nonrecursive signal-decomposition algorithm. The intrinsic mode function (IMF) in different frequency bands is obtained by adaptively decomposing signals. In the process of solving the mode function, image continuation is used to avoid the endpoint effect in EMD (empirical mode decomposition) and other decomposition methods. VMD’s processing of nonlinear fault signals is helpful for extracting the characteristics of subsequent fault signals. The decomposed IMF has an independent center frequency and sparsity and can effectively avoid mode aliasing when the parameters are appropriate.
(1) In order to obtain a unidirectional spectrum, Hilbert transform is used to analyze and calculate each modal signal, and then the frequency shifting method is used to move the modal spectrum to the baseband:
δ ( t ) + j π t v k ( t ) e j ω k t
δ ( t ) is the Dirac function; v k and ω k are the first IMF component and its center frequency decomposed, respectively.
(2) The bandwidth is estimated by the square norm of the gradient, and the constraint expression is:
min { v k } , { ω k } k = 1 K t δ t + j π t v k ( t ) e j ω k t 2 2 s . t .                 k = 1 K v k ( t ) = f ( t )
t is the gradient calculation, and * is the convolution calculation symbol.
(3) In order to obtain the optimal solution more efficiently, the constrained problem is transformed into an unconstrained problem by using a Lagrangian operator τ ( t ) and penalty factor α . The expanded Lagrangian function expression is:
L ( { v k } , { ω k } , τ ) = α k = 1 K t δ t + j π t v k ( t ) e j ω k t 2 2                             + f ( t ) k = 1 K v k ( t ) 2 2 + τ ( t ) , f ( t ) k = 1 K v k ( t )
(4) The alternating direction multiplier method is used to iteratively update each modal component and the center frequency, and the saddle point of the unconstrained function is obtained, which is the best solution to the problem. The iterative update expression of v k , ω k , τ is as follows:
v ^ k n + 1 ( ω ) = f ^ ( ω ) i k v ^ i ( ω ) + τ ^ ( ω ) 2 1 + 2 α ( ω ω k ) 2
ω k n + 1 = 0 ω v ^ k ( ω ) 2 d ω 0 v ^ k ( ω ) 2 d ω
τ ^ n + 1 ( ω ) = τ ^ n ( ω ) + γ f ^ ( ω ) k = 1 K u ^ k n + 1 ω
In Equation (6) n is the number of iterations and γ is the noise tolerance parameter.
(5) Judgment of the iteration termination condition.
k = 1 K v ^ k n + 1 v ^ k n 2 2 k = 1 K v ^ k n 2 2 < ε
where ε is the judgment accuracy. When ε > 0, the iteration stops and the cycle ends. The step flow chart is shown in Figure 1:

3. Fractal Dimension Calculation Method Based on Mathematical Morphology

The key to estimating the fractal dimension is to measure the signal at different scales, which can be realized using multi-scale morphology [12]. In the process of gathering and covering the signal, the supremum function, i.e., the structural element g (n), is used for equivalent transformation, and the scale range ε = 1 , 2 , , ε max is analyzed. The algorithm includes two basic operators: expansion and corrosion. For a one-dimensional discrete time signal f ( n ) ( n = 0 , 1 , 2 , ) , the expansion and corrosion results of each scale are, respectively,
f Θ g ε ( n ) = ( ( f Θ g ) Θ g ) Θ g ε T i m e s
f g ε ( n ) = ( ( f g ) g ) g ε T i m e s
where Θ represents corrosion operation; represents expansion operation; g ε ( n ) represents the structural elements used in the scale ε ; f represents signal; g represents a structural element; and ε T i m e s is the number of expansion and corrosion. The coverage area of discrete signals at different scales is defined as
A g ( ε ) = n = 1 N ( f ε g ( n ) f Θ ε g ( n ) )
When ε approaches zero, A g ( ε ) satisfies:
lg A g ( ε ) ε 2 D M lg 1 ε + c
where D M is the Minkowski–Bouligand dimension of the signal, and c is a constant. The slope of the straight line obtained by a least-square fitting of the above equations lg A g ( ε ) ε 2 and lg 1 ε is D M , which is the final required f ( n ) fractal dimension.
Generally, the unit structure element g ( n ) is chosen to be { 0 , 0 , 0 } , because this structure not only ensures that the dimension estimation is not affected by the signal amplitude range but also reduces the computational complexity of the algorithm. In principle, the maximum analysis scale ε max is a positive integer less than N / 2 ( N is the number of discrete signal sampling points). When the data length is relatively large, appropriately reducing ε max can reduce the amount of calculation [13]. In this paper, the length is 256.

4. Kernel Principal Component Analysis

The KPCA method is used to represent the nonlinear relationship between modeling data. It can effectively project linearly indivisible input data into a high-dimensional space that can be linearly separated, and then execute linear P C A in feature space H . Assume that the sample set is x 1 , x 2 , , x N R m , where: N is the number of samples and m is the number of variables. These samples are projected into the feature space H through a nonlinear mapping ϕ , which can be expressed as:
ϕ : x i R m ϕ ( x i ) R h
where h is the dimension of the feature space.
The dot product of vector ϕ ( x i ) and ϕ ( x j )   ( i   ,     j = 1   ,     2   ,       ,     N ) in feature space H is:
ϕ ( x i ) T ϕ ( x j ) = ϕ ( x i ) , ϕ ( x j ) = k ( x i , x j )
where: k is the kernel function.
The covariance matrix of samples in the high-dimensional feature space H is:
Q = 1 N 1 j = 1 N ϕ ( x j ) ϕ ( x j ) T
Similar to linear PCA, KPCA in the feature space is equivalent to solving the eigenvalue problem. Let the eigenvector be V k , the eigenvalue be μ k , and its characteristic equation be:
μ k V k = Q V k = 1 N 1 j = 1 N ϕ ( x j ) ϕ ( x j ) T V k =                         1 N 1 j = 1 N ϕ ( x j ) , V k ϕ ( x j )
In Equation (15), each eigenvector V k in the covariance matrix Q can be regarded as a linear combination of ϕ ( x 1 ) , ϕ ( x 2 ) , , ϕ ( x N ) [14], that is:
V k = i = 1 N v k , i ϕ ( x i ) ,                         k = 1 , 2 , , N
where v k and i are linear correlation coefficients.
Combining Equations (13) and (15), we can obtain
μ k i = 1 N v k , i K k , i = 1 N 1 i = 1 N v k , i j = 1 N K k , i K j i
Further simplification results in:
λ k v k = K v k
where λ k = ( N 1 ) μ k ; K R N × N is a kernel matrix, whose eigenvalues are λ 1 λ 2 λ N , and the corresponding eigenvectors are v 1 , v 2 , , v N . In the feature space H , the eigenvector V k of the covariance matrix Q in Equation (15) forms a matrix V f = V 1       V 2         V l         V l + 1         V N . One principal element (PC) is selected, and we obtain V ^ f = V 1       V 2         V l R N × l in the principal element space, since the eigenvector V k should meet the normal constraint in the feature space H , that is, V k , V k = 1 . Thus, the eigenvectors v 1 , v 2 , , v N of the kernel matrix K can be expressed as:
v k , v k = V k , V k λ k
The norm of 1 λ k is taken as the eigenvector, so we have
v k = 1 λ k v k
V ^ f in feature space H is expressed as
V ^ f = 1 λ 1 X T v l = X T V ^ Λ 1 / 2
where Λ = d i a g ( λ 1 , λ 2 , , λ l ) ; V ^ = [ v 1     v 2         v l   ] ; X = [ ϕ ( x 1 )         ϕ ( x 2 )         ϕ ( x N ) ] T .

5. Basic Theory of SVM

SVM can provide good generalization performance, and it has been widely used in small-sample classification [15]. As a classical nonlinear problem in the area of fault-location problems, SVM cannot find the classification hyperplane in the original sample space, so it is necessary to introduce a kernel function K ( x i , x j ) to map to a linearly separable high-dimensional space. In this paper, a Gaussian kernel function mapping transformation is used, and the expression is as follows:
K ( x i , x j ) = exp γ x i x j 2
where, x i x j 2 is the square of the distance between x i and x j , and r is the super parameter of the kernel function itself:
min ω , b , ζ 1 2 ω T ω + c i = 1 N ζ i s . t .             y i ( ω T ϕ ( x i ) ) 1 ζ i ,             i = 1 , 2 , N                       ζ i 0 ,             i = 1 , 2 , , N
where ζ i is the relaxation variable, c is the penalty factor, x i is the input i sample, and y i is the sample type.
The Lagrangian penalty factor α is introduced, and finally the dual type of soft interval Gaussian kernel SVM is obtained as follows:
min α 1 2 i = 1 N j = 1 N α i α j y i y j exp γ x i x j 2 i = 1 N α i s . t .       0 α i c ,             i = 1 , 2 , N ,                       i = 1 N α i y i = 0

6. Experiment and Simulation

The method based on VMD-MMFD is used for feature extraction and the fault diagnosis of fault circuits. The specific process is shown in Figure 2.
To verify the effectiveness of this method, simulation experiments and calculations were carried out through the software Multisim, Matlab, and Python. The Sallen–Key band-pass filter circuit was taken as an example to verify the analysis method of soft-fault diagnosis based on variational-mode decomposition and mathematical morphology fractal dimensions. The circuit diagram is shown in Figure 3.
The excitation signal f = 3sin (30) + 3sin (20) + noise is used as the input source for the fault diagnosis, the tolerance of resistance and capacitance of the circuit is 10%, and the value of components within this range is considered normal. In this experiment, the soft fault state is set when the nominal value fluctuates by 30%. In this experiment, 15-parameter fault combinations are configured. The fault of the amplifier is affected by many factors, so it is not discussed here. The fault mode setting is shown in Table 1, the original signal is shown in Figure 4, the noise signal is shown in Figure 5, and the mixed noise input signal is shown in Figure 6.
When VMD is used to decompose signals, the preset scale parameters and the second penalty factor are the main parameters that affect the decomposition accuracy. Therefore, for the VMD decomposition of measured signals, the reasonable selection of its parameters is the difficulty and key of this method. The center frequency of each order of IMF component obtained from the VMD decomposition of the signal is distributed from a low frequency to a high frequency. If the optimal preset scale parameter K is obtained, the center frequency of the last-order IMF component should be the maximum value, and the maximum center frequency value should remain stable. Based on the test and analysis of the VMD decomposition results of a large number of measured signals and reference [16], the second penalty factor α = 2000 is selected in this paper. Take fault type 3 and fault type 5 signals as examples of VMD decomposition, as shown in Figure 7 and Figure 8, the center frequency of the IMF component is obtained by decomposition, as shown in Figure 9 and Figure 10. It can be seen from the figures that when K = 4, the central frequency of the IMF component reaches the maximum and tends to be stable, the frequencies between modes do not overlap, and the impact of noise and mode aliasing is effectively suppressed. With the preset scale parameter K > 4, the center frequencies of IMF 5 and IMF 6 components become unstable and mode aliasing occurs. Therefore, the VMD decomposition of the signal is the best when k = 4.
In recent years, in the field of circuit fault diagnosis, the method of combining empirical mode decomposition with box-dimension estimation has been commonly used to extract features [17]. The experiment compares this traditional method with the model of feature extraction based on variational mode decomposition and mathematical morphology. Take four states of types 2–5 as examples, where q is the scale and the value range is [0, 10]. Each state corresponds to 10 fractal dimensions. With the increase in the scale q, fault states are depicted from different dimensions, Dq is the fractal dimension value under different scales q, and the spectrum analysis is carried out.
In Figure 11, the IMF1 component of Fault03 does not overlap with other components, and the IMF1 components of Fault02, Fault04, and Fault05 have aliasing and different degrees of crossing. With the increase in q, the aliasing and crossing become more and more serious, making them difficult to distinguish. In Figure 12, the IMF1 components of Fault02 and Fault04 have serious crossing phenomena, and it is difficult to distinguish the fault types. When the q value is small, the IMF1 of Fault03 and Fault05 have no aliasing phenomena. With the increase in q, the IMF1 of Fault04 has serious cross phenomena with the IMF1 of Fault03 and Fault05. In Figure 13, the IMF1 of Fault04 and Fault05 has seriously crossed, and it is impossible to distinguish the fault type. With the increase in the q value, the IMF1 of Fault03 and Fault04 has also mixed or even crossed phenomena. In Figure 14, when the q value is low, each fault type can be distinguished. With the increase in the q value, the IMF1 of Fault04 and Fault05 has aliasing and crossing phenomena, which affect the discrimination of the fault types. Later, with the increase in the scale of q, the IMF1 of Fault04 and Fault03-IMF1 has a certain degree of aliasing. The characteristic sample set of the IMF1 component dimension of the fault signal is shown in Table 2.
In Figure 15, the IMF2 of Fault04 and Fault05 has undergone serious aliasing and crossing, and the fault type of the IMF2 of Fault02 and Fault03 can be distinguished, in which aliasing and crossing have occurred. In Figure 16, the IMF2 of Fault02 and Fault03 has been separated, and the fault type is well distinguished, but the IMF2 of Fault04 and Fault05 shows aliasing with the increase in the q value, which is difficult to distinguish. In Fault02 and Fault04 in Figure 17, there is no aliasing, the effect of the fault type discrimination is good, and the IMF2 of Fault02 and Fault04 has serious aliasing and cross phenomena, so it is impossible to distinguish fault types. In Figure 18, when the scale of q is low, the state-differentiation effect of various fault types is good. With the increase in the q value, the IMF2 of Fault02, Fault03, Fault04, and Fault05 shows aliasing and crossing phenomena to different degrees. The characteristic sample set of the IMF2 component dimension of the fault signal is shown in Table 3.
In Figure 19, with the increase in the number of IMF components, the dimension values of several components of the IMF are no longer stable, because with the increase in the number of modal decomposition components, each component contains less and less original feature information. After the dimension calculation process, the fault types become difficult to distinguish. When the scale of q is small, Fault04 and Fault05 have been aliased and crossed to some extent. With the increase in the scale of the q value, the IMF3 of Fault02, Fault04, and Fault05 has been aliased and crossed to some extent. In Figure 20, the IMF3 aliasing of Fault02, Fault03, Fault04, and Fault05 is obviously crossed, which makes it difficult to distinguish the fault type*. In Figure 21, the IMF3 fault type * of Fault02 is not aliased with other components, and the IMF3 of Fault03, Fault04, and Fault05 is completely aliased, so the fault types cannot be distinguished. In Figure 22, the IMF3 of Fault03 and Fault05 is seriously aliased and crossed, which makes it impossible to classify fault types. The IMF3 of Fault04 and Fault02 is aliased and crossed with the increase in the scale of q, which makes it difficult to classify faults. The sample set of dimension characteristics of the IMF3 components of fault signals is shown in Table 4.
In Figure 23, with the increase in the number of IMF components, the dimension values of several IMF components fluctuate seriously, and the IMFs of Fault02, Fault03, Fault04, and Fault05 have all experienced serious aliasing and crossing. In Figure 24, the IMF4 in Fault02, Fault04, and Fault05 is aliased, so it is impossible to distinguish the fault type, and with the increase in the q value, the IMF3 in Fault03 is also aliased. In Figure 25, the IMF4 of Fault02, Fault04, and Fault05 shows serious aliasing and cross phenomena and even overlaps with the increase in the scale q, making it difficult to determine the fault type. Only the fault type of the IMF4 of Fault03 can be distinguished. In Figure 26, the IMF dimension maps of the four fault states fluctuate unstably, and the IMF4 of Fault02, Fault03, Fault04, and Fault05 is seriously aliased and crossed, making it impossible to judge the fault type. The characteristic sample set of the IMF4 component dimension of the fault signal is shown in Table 5.
The feature extraction of different states of signals is carried out through the fractal dimension. Although the number of data samples increases and the fault set becomes larger, the high-dimensional data contain a lot of miscellaneous and repetitive data, which seriously affect the accuracy of fault diagnosis. The feature set needs to be dimensionally reduced. In the literature [18], different kernel functions have different effects on data dimensionality reduction. Generally, the Gaussian radial basis function (Formula (25)) is selected to perform the kernel function PCA on data sets to achieve dimensionality reduction and fault classification in high-dimensional feature space.
k x , y = exp x y 2 2 σ 2
The cumulative contribution rate of the principal component calculated by the KPCA algorithm is shown in Table 6. It can be seen from Table 6 that the first three principal components contain almost all the information of the system, so the first three principal components are retained as the eigenvalues of the system; that is, the feature space is reduced from high-dimensional to three-dimensional. KPCA can be used to calculate the principal component contribution rate, effectively eliminating the samples containing too much interference information, which greatly improves the accuracy of subsequent diagnoses.
Therefore, considering the feature extraction method of the VMD-MMFD-KPCA, the results are shown in Figure 27, and various fault types can be distinguished in the following figure.
To further verify the effectiveness of the feature-extraction methods explained above, a support vector machine was used to complete fault diagnosis. Among them, the fault datasets that were used were all from the fault settings of various components of the simulation circuit. To achieve better experimental results, more than 50 signal samples were collected for each fault type, and more than 750 signal samples were collected for 15 fault types. Eighty percent of the samples from each fault type were randomly selected as the training sample, and the rest were selected as the test sample. The kernel function was introduced and determined in the previous theory, using Gaussian kernel function. After that, the data needed to be normalized. First, the data need to be filtered and sorted. For obviously abnormal data, we considered using the sample mean to replace it, so as to reduce the extent to which the abnormal sample data interferes with the results. To prevent the individual abnormal sample data volume from influencing the training results, the data were normalized according to Formula (26).
X i = X i X min X max X min
After 100 iterations, diagnosis is carried out in combination with different diagnosis models, as shown in Figure 28.
In the literature [19], there is a wavelet neural network (WNN) analog circuit soft fault diagnosis method based on kernel partial least square (KPLS) feature extraction, which uses support vector machine (KPLS-WNN-SVM) for fault classification. In the literature [20], a multifractal extraction method based on empirical mode decomposition is combined with support vector machine (EMD-MFD-SVM). The diagnostic rate of the diagnostic model combined with support vector machine (VMD-MMFD-KPCA-SVM) in this paper is shown in Table 7, which shows that the VMD-MMFD-KPCA-SVM method has the highest diagnostic accuracy.

7. Conclusions

This paper proposes a feature-extraction method based on VMD-MMFD-KPCA. Firstly, the fault signal is decomposed by VMD to obtain multiple modal IMF components, which greatly reduces the impact of modal aliasing. The appropriate IMF components were selected, and the fractal dimension of the mathematical morphology was calculated to obtain a large number of high-dimensional feature data sets. The data sets were reduced by KPCA to obtain low-dimensional data, which eliminated a large number of redundant and jumbled information in the high-dimensional data. The experimental results show that the feature-extraction method based on VMD-MMFD-KPCA has a good effect.

Author Contributions

Conceptualization, X.L.; methodology, X.L.; software, S.S.; validation, X.L., K.L., Q.W. and D.S.; writing—original draft preparation, C.Y., Z.L. and J.W.; writing—review and editing, C.Y.; project administration, X.L.; funding acquisition, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Heilongjiang Provincial Natural Science Foundation of China, grant number YQ2022F014.

Acknowledgments

The authors acknowledge Heilongjiang Provincial Natural Science Foundation of China (grant number YQ2022F014).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peng, W.; Yang, S. A new diagnosis approach for handling tolerance in analog and mixed-signal circuits by using fuzzy math. IEEE Trans. Circuits Syst. I Regul. Pap. 2005, 52, 2118–2127. [Google Scholar] [CrossRef]
  2. Zhu, Y.J.; Zhang, W.M.; Zhang, S.T.; Chen, J.J. Judgment and handling of “soft fault” of discrete component circuit. Chin. Sci. Technol. J. Database Eng. Technol. 2016, 9, 00300. [Google Scholar]
  3. Lu, X.M.; Zhao, H.; Wu, Q. Research on fractal method for soft fault diagnosis of nonlinear analog circuits. Model. Meas. Control. A 2017, 90, 58–73. [Google Scholar] [CrossRef]
  4. Han, H.T.; Ma, H.G.; Cao, J.F.; Zhang, J.L. Analog circuit fault diagnosis method based on nonlinear spectrum characteristics and kernel principal component analysis. J. Electr. Technol. 2012, 27, 248–2544. [Google Scholar]
  5. Wang, Y.H.; Lu, J.; Pan, G.Q.; Feng, J.C. Research on Feature Extraction for Diagnostics of Analog Circuit Based on Factor Analysis. Comput. Meas. Control 2014, 11, 25–27. [Google Scholar]
  6. Xiao, Y.; Feng, L. A novel neural-network approach of analog fault diagnosis based on kernel discriminant analysis and particle swarm optimization. Appl. Soft Comput. 2012, 12, 904–920. [Google Scholar] [CrossRef]
  7. Guo, Q.; Zhang, W.B.; Su, H.T. Analog circuit fault diagnosis based on WPT and FOAGRNN. Comput. Simul. 2020, 37, 355–359. [Google Scholar]
  8. Ji, T.Y.; Lu, Z.; Wu, Q.H. Detection of power disturbances using morphological gradient wavelet. Signal Process. 2008, 88, 255–267. [Google Scholar] [CrossRef]
  9. Zhuang, N.; Zeng, Y.; Tong, L.; Zhang, C.; Zhang, H.; Yan, B. Emotion Recognition from EEG Signals Using Multidimensional Information in EMD Domain. Biomed Res. Int. 2017, 2017, 317–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Zheng, Z.; Jiang, W.L.; Wang, Z.W.; Zhu, Y.; Yang, K. Gear fault diagnosis method based on local mean decomposition and generalized morphological fractal dimensions. Mech. Mach. Theory 2015, 91, 151–167. [Google Scholar] [CrossRef]
  11. Nayak, S.R.; Mishra, J. Fractal Dimension Based Generalized Box-Counting Technique with Application to Grayscale Images. Fractals 2020, 3, 8. [Google Scholar] [CrossRef]
  12. Zhang, K.; Cheng, J.S.; Yang, Y. A rolling bearing fault diagnosis method based on local mean decomposition and morphological fractal dimension. J. Vib. Shock. 2013, 32, 90–94. [Google Scholar]
  13. Wang, F.J.; Duan, S.L.; Yu, H.L.; Li, H.K. Diesel Engine Fault Diagnosis Based on EEMD and Morphological Fractal Dimension. J. Intern. Combust. Engine 2012, 30, 557–562. [Google Scholar]
  14. Yang, G.; Gu, X. Fault Diagnosis of Complex Chemical Processes Based on Enhanced Naive Bayesian Method. IEEE Trans. Instrum. Meas. 2020, 69, 4649–4658. [Google Scholar] [CrossRef]
  15. Gashteroodkhani, O.A.; Majidi, M.; Etezadi-Amoli, M.; Nematollahi, A.F.; Vahidi, B. A hybrid SVM-TT transform-based method for fault location in hybrid transmission lines with underground cables. Electr. Power Syst. Res. 2019, 170, 205–214. [Google Scholar] [CrossRef]
  16. Tang, G.J.; Wang, X.L. Variational modal decomposition method and its application in early fault diagnosis of rolling bearings. J. Vib. Eng. 2016, 29, 11. [Google Scholar]
  17. Han, D.Y.; Li, G.; Shi, P.M. Research on coupling fault diagnosis method of rotating machinery based on EMD and fractal box dimension. Vib. Shock. 2013, 32, 209–214. [Google Scholar]
  18. Li, C.Y.; Liu, G.D.; Tan, B. Comparative study of hyperspectral remote sensing data dimensionality reduction based on PCA and KPCA. Geospat. Inf. 2022, 20, 89–93+103. [Google Scholar]
  19. Cong, W.; Jing, B.; Yu, H.K. WNN soft fault diagnosis of analog circuit based on KPLS feature extraction. J. Cent. South Univ. (Sci. Technol.) 2014, 45, 1841–1846. [Google Scholar]
  20. Lu, X.M.; Zhao, H.; Wu, Q. Generalized Multi-fractal Method for Soft Fault Feature Extraction of Non-linear Analog Circuits Based on EMD. Rev. De La Fac. De Ing. 2016, 31, 192–203. [Google Scholar]
Figure 1. Variational Modal Decomposition Process.
Figure 1. Variational Modal Decomposition Process.
Electronics 12 00184 g001
Figure 2. VMD-MMFD Method Flow Chart.
Figure 2. VMD-MMFD Method Flow Chart.
Electronics 12 00184 g002
Figure 3. Sallen−Key band-pass filter circuit.
Figure 3. Sallen−Key band-pass filter circuit.
Electronics 12 00184 g003
Figure 4. Original Signal.
Figure 4. Original Signal.
Electronics 12 00184 g004
Figure 5. Noise Signal.
Figure 5. Noise Signal.
Electronics 12 00184 g005
Figure 6. Input signal of mixed noise.
Figure 6. Input signal of mixed noise.
Electronics 12 00184 g006
Figure 7. VMD Decomposition Sequence Diagram of Fault 3 Signal.
Figure 7. VMD Decomposition Sequence Diagram of Fault 3 Signal.
Electronics 12 00184 g007
Figure 8. VMD Decomposition Sequence Diagram of Fault 5 Signal.
Figure 8. VMD Decomposition Sequence Diagram of Fault 5 Signal.
Electronics 12 00184 g008
Figure 9. IMF Component Corresponding Spectrum of Fault 3.
Figure 9. IMF Component Corresponding Spectrum of Fault 3.
Electronics 12 00184 g009
Figure 10. IMF Component Corresponding Spectrum of Fault 5.
Figure 10. IMF Component Corresponding Spectrum of Fault 5.
Electronics 12 00184 g010
Figure 11. IMF1-EMD box fractal dimension map.
Figure 11. IMF1-EMD box fractal dimension map.
Electronics 12 00184 g011
Figure 12. IMF1-EMD mathematical morphology fractal dimension map.
Figure 12. IMF1-EMD mathematical morphology fractal dimension map.
Electronics 12 00184 g012
Figure 13. IMF1-VMD box fractal dimension map.
Figure 13. IMF1-VMD box fractal dimension map.
Electronics 12 00184 g013
Figure 14. IMF1-VMD mathematical morphology fractal dimension map.
Figure 14. IMF1-VMD mathematical morphology fractal dimension map.
Electronics 12 00184 g014
Figure 15. IMF2-EMD box fractal dimension map.
Figure 15. IMF2-EMD box fractal dimension map.
Electronics 12 00184 g015
Figure 16. IMF2-EMD mathematical morphology fractal dimension map.
Figure 16. IMF2-EMD mathematical morphology fractal dimension map.
Electronics 12 00184 g016
Figure 17. IMF2-VMD box fractal dimension map.
Figure 17. IMF2-VMD box fractal dimension map.
Electronics 12 00184 g017
Figure 18. IMF2-VMD mathematical morphology fractal dimension map.
Figure 18. IMF2-VMD mathematical morphology fractal dimension map.
Electronics 12 00184 g018
Figure 19. IMF3-EMD box fractal dimension map.
Figure 19. IMF3-EMD box fractal dimension map.
Electronics 12 00184 g019
Figure 20. IMF3-EMD mathematical morphology fractal dimension map.
Figure 20. IMF3-EMD mathematical morphology fractal dimension map.
Electronics 12 00184 g020
Figure 21. IMF3-VMD box fractal dimension map.
Figure 21. IMF3-VMD box fractal dimension map.
Electronics 12 00184 g021
Figure 22. IMF3-VMD mathematical morphology fractal dimension map.
Figure 22. IMF3-VMD mathematical morphology fractal dimension map.
Electronics 12 00184 g022
Figure 23. IMF4-EMD box fractal dimension map.
Figure 23. IMF4-EMD box fractal dimension map.
Electronics 12 00184 g023
Figure 24. IMF4-EMD mathematical morphology fractal dimension map.
Figure 24. IMF4-EMD mathematical morphology fractal dimension map.
Electronics 12 00184 g024
Figure 25. IMF4-VMD box fractal dimension map.
Figure 25. IMF4-VMD box fractal dimension map.
Electronics 12 00184 g025
Figure 26. IMF4-VMD mathematical morphology fractal dimension map.
Figure 26. IMF4-VMD mathematical morphology fractal dimension map.
Electronics 12 00184 g026
Figure 27. VMD-MMFD-KPCA feature extraction.
Figure 27. VMD-MMFD-KPCA feature extraction.
Electronics 12 00184 g027
Figure 28. Flow chart of SVM Classification.
Figure 28. Flow chart of SVM Classification.
Electronics 12 00184 g028
Table 1. Fault Type Description.
Table 1. Fault Type Description.
Type NumberTypeNominal ValueFault Zone
1Normal--
2R1  30%4 kΩ[4.4 kΩ, 5.2 kΩ]
3R1  30%4 kΩ[2.8 kΩ, 3.6 kΩ]
4R2 30%3 kΩ[3.3 kΩ, 3.9 kΩ]
5R2  30%3 kΩ[2.1 kΩ, 2.7 kΩ]
6R3  30%2 kΩ[2.2 kΩ, 2.6 kΩ]
7R3 30%2 kΩ[1.4 kΩ, 1.8 kΩ]
8R4 30%4 kΩ[4.4 kΩ, 5.2 kΩ]
9R4 30%4 kΩ[2.8 kΩ, 3.6 kΩ]
10R5 30%4 kΩ[4.4 kΩ, 5.2 kΩ]
11R5 30%4 kΩ[2.8 kΩ, 3.6 kΩ]
12C1 30%5 nF[5.5 nF, 6.5 nF]
13C1 30%5 nF[3.5 nF, 4.5 nF]
14C2 30%5 nF[5.5 nF, 6.5 nF]
15C2 30%5 nF[3.5 nF, 4.5 nF]
Table 2. Characteristic sample of the fault signal analysis dimension.
Table 2. Characteristic sample of the fault signal analysis dimension.
Fault TypeD1D2D3D4D5D6D7D8D9D10
IMF1-EMD box fractal dimensionFault 21.68701.67431.60611.61561.55671.48381.40951.49101.39741.4057
Fault 31.36021.37351.28531.29081.23031.26891.25951.24191.16551.2092
Fault 41.68541.65081.58351.60741.45561.41231.40731.33941.34621.1949
Fault 51.40711.42561.41331.45721.40971.40851.38321.41201.37961.4140
IMF1-EMD mathematical morphology fractal dimensionFault 21.65571.59951.66611.61241.50501.49231.48901.48201.47851.3868
Fault 31.29081.24971.25131.20511.24821.20101.17551.19801.14381.1610
Fault 41.58351.59811.53371.54021.54521.53721.41511.39501.40241.3920
Fault 51.42451.42881.42311.43051.38631.43541.35711.42901.32211.3840
IMF1-VMD box fractal dimensionFault 21.63771.66301.64471.60471.59881.54181.51821.54851.47241.4856
Fault 31.21761.21861.23421.20471.16361.16461.17091.09931.12361.0548
Fault 41.54391.51641.47291.44711.34721.34571.25161.19491.14851.0201
Fault 51.51521.40021.47891.46001.33851.29721.29551.29961.29571.2398
IMF1-VMD mathematical morphology fractal dimensionFault 21.75641.76291.74291.74801.70941.72061.65821.69471.56851.6582
Fault 31.35261.33091.34221.31641.29541.31681.32811.22941.27861.2202
Fault 41.67801.63331.55081.60141.54631.46241.39511.34721.32411.2779
Fault 51.46361.40811.42401.40961.38621.41581.40751.36441.34651.3457
Table 3. Characteristic sample of the fault signal analysis dimension.
Table 3. Characteristic sample of the fault signal analysis dimension.
Fault TypeD1D2D3D4D5D6D7D8D9D10
IMF2-EMD box fractal dimensionFault 21.67241.64091.59601.61041.58871.45921.50101.44421.38031.3475
Fault 31.20331.19101.27491.18321.24771.16371.13191.17491.09821.0771
Fault 41.57041.53661.53001.43051.43351.40651.39181.35011.33401.2630
Fault 51.47321.49911.47651.46391.37281.37071.34551.42091.34711.3389
IMF2-EMD mathematical morphology fractal dimensionFault 21.71911.69301.68061.77531.70651.64131.68321.65371.59671.6174
Fault 31.32621.33471.35541.30261.27841.30801.33051.22701.20341.1797
Fault 41.63671.59091.56301.54501.45991.45361.42991.37771.25981.2645
Fault 51.43371.44281.44551.39151.44171.39871.41691.38431.33391.3127
IMF2-VMD box fractal dimensionFault 21.66821.64301.64651.59381.62491.55141.52361.54581.52851.4540
Fault 31.41241.35321.31721.29331.34071.23311.20431.19301.14551.1187
Fault 41.53541.47611.46571.48081.43481.38261.33221.27791.23301.2255
Fault 51.34211.38181.32691.25321.22271.30781.23881.25531.15841.1871
IMF2-VMD mathematical morphology fractal dimensionFault 21.83361.76501.66971.66391.52061.45091.37641.23511.09541.0255
Fault 31.20661.20761.18191.17401.15431.16401.13721.09741.08671.0477
Fault 41.58621.58941.55841.48991.43711.36271.33941.26671.17441.0933
Fault 51.39141.40411.32281.29011.24951.18341.13111.09921.02160.9455
Table 4. Characteristic sample of the fault signal analysis dimension.
Table 4. Characteristic sample of the fault signal analysis dimension.
Fault TypeD1D2D3D4D5D6D7D8D9D10
IMF3-EMD box fractal dimensionFault 21.64251.58151.65681.55381.55911.51631.39471.39491.29071.2155
Fault 31.20931.27271.28971.19411.23411.19911.10681.15991.15441.0175
Fault 41.53981.47721.45631.47781.40281.34081.41951.34291.32431.2693
Fault 51.46921.53791.49321.37611.45321.28041.24491.18681.19601.1634
IMF3-EMD mathematical morphology fractal dimensionFault 21.52281.51091.47541.47981.50591.49501.44241.35191.34301.3664
Fault 31.37591.30591.36221.27731.27271.27931.24581.20061.11771.1418
Fault 41.58321.48831.50461.46101.35081.24761.19251.17731.08220.9293
Fault 51.40001.44091.39511.41751.43871.39771.34691.38751.29881.2890
IMF3-VMD box fractal dimensionFault 21.64951.68441.60461.52051.54951.46871.37871.38391.30421.2562
Fault 31.41911.30631.35251.22121.17641.12121.02170.84400.80910.7661
Fault 41.55371.49011.38371.40141.32901.20421.06991.08520.97330.8581
Fault 51.34491.37891.34861.29151.23831.14521.13821.07591.10571.0257
IMF3-VMD mathematical morphology fractal dimensionFault 21.77961.76421.74621.62091.49961.47731.36461.21681.03150.9121
Fault 31.31301.23901.25531.13221.11811.11541.04901.00630.88970.8612
Fault 41.58791.60341.60841.52441.45761.41531.31281.29011.18241.1364
Fault 51.36901.28181.21761.32541.28671.09910.96290.86770.98530.7289
Table 5. Characteristic sample of the fault signal analysis dimension.
Table 5. Characteristic sample of the fault signal analysis dimension.
Fault TypeD1D2D3D4D5D6D7D8D9D10
IMF4-EMD box fractal dimensionFault 21.70861.60091.67041.61451.66591.63131.57141.52831.42931.4673
Fault 31.36361.19811.33291.33071.20951.27521.25521.10511.18301.1496
Fault 41.63611.62151.52151.42401.51641.41661.46241.29871.39371.3317
Fault 51.45431.52681.50211.44391.32311.27081.39801.30521.16141.1513
IMF4-EMD mathematical morphology fractal dimensionFault 21.62371.57081.46931.53821.53811.54341.36511.35731.41111.3618
Fault 31.21921.19131.22631.16841.21601.09611.15271.09130.98661.0166
Fault 41.57511.46991.56101.51051.33411.22441.30161.13321.01310.9396
Fault 51.45141.44991.29241.42631.35371.36991.37631.21501.34151.2034
IMF4-VMD box fractal dimensionFault 21.73921.64741.59631.63301.48001.43251.58021.37411.27971.1849
Fault 31.18071.14131.11571.05960.90710.91120.78030.66630.61730.3559
Fault 41.66821.64141.49841.40611.52061.36891.26971.30841.12581.0662
Fault 51.39891.45841.28691.25341.20271.32111.27891.26851.10851.1401
IMF4-VMD mathematical morphology fractal dimensionFault 21.84311.75121.71051.62161.56551.55351.49401.39291.12520.9969
Fault 31.43501.32621.26001.27371.11331.15321.01510.88101.02700.7543
Fault 41.62081.65091.56601.56901.59821.35701.36411.18891.12231.1073
Fault 51.30011.50731.45251.38541.16311.32421.17711.05300.92130.9574
Table 6. Cumulative contribution rate.
Table 6. Cumulative contribution rate.
Cumulative Contribution RateKPCA1KPCA2KPCA3
Fault 2 (%)93.2695.5998.88
Fault 3 (%)95.3396.2199.11
Fault 4 (%)96.3498.88100
Fault 5 (%)97.7799.18100
Table 7. Fault Diagnosis Rate.
Table 7. Fault Diagnosis Rate.
Diagnostic ModelAverage Diagnostic Rate (%)
KPLS-WNN-SVM87.3
EMD-MFD-SVM91.9
VMD-MMFD-KPCA-SVM96.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lu, X.; Yang, C.; Wu, Q.; Wang, J.; Lu, Z.; Sun, S.; Liu, K.; Shao, D. Research on Analog Circuit Soft Fault Diagnosis Method Based on Mathematical Morphology Fractal Dimension. Electronics 2023, 12, 184. https://doi.org/10.3390/electronics12010184

AMA Style

Lu X, Yang C, Wu Q, Wang J, Lu Z, Sun S, Liu K, Shao D. Research on Analog Circuit Soft Fault Diagnosis Method Based on Mathematical Morphology Fractal Dimension. Electronics. 2023; 12(1):184. https://doi.org/10.3390/electronics12010184

Chicago/Turabian Style

Lu, Xinmiao, Cunfang Yang, Qiong Wu, Jiaxu Wang, Zihan Lu, Shuai Sun, Kaiyi Liu, and Dan Shao. 2023. "Research on Analog Circuit Soft Fault Diagnosis Method Based on Mathematical Morphology Fractal Dimension" Electronics 12, no. 1: 184. https://doi.org/10.3390/electronics12010184

APA Style

Lu, X., Yang, C., Wu, Q., Wang, J., Lu, Z., Sun, S., Liu, K., & Shao, D. (2023). Research on Analog Circuit Soft Fault Diagnosis Method Based on Mathematical Morphology Fractal Dimension. Electronics, 12(1), 184. https://doi.org/10.3390/electronics12010184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop