# Multifractal Analysis for Soft Fault Feature Extraction of Nonlinear Analog Circuits

Authors: Xinmiao Lu; Hong Zhao; Haijun Lin; Qiong Wu
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

Aiming at the nonstationarity and nonlinearity of soft fault signals of nonlinear analog circuits, the use of multifractal detrended fluctuation analysis can effectively reveal the dynamic behavior hidden in multiscale nonstationary signals. This paper adopts a new method that uses multifractal detrended fluctuation analysis to calculate the multifractal singularity spectrum of soft fault signals of nonlinear analog circuits. Moreover, this method endows the parameters of the spectrum with definite physical meanings including width, maximum singular index, minimum singular index, and corresponding singularity index of the extreme point. Therefore, this method can be applied to characterize the internal dynamic mechanism of the soft fault signals of nonlinear analog circuits, making it suitable for the feature extraction of fault circuits. All multifractal feature parameters can be organized into a feature set, which will be then input to a support vector machine, and fault detection for the nonlinear analog circuit can be conducted via the support vector machine.


---

## Body

## 1. Introduction

The testing and fault diagnosis of nonlinear analog circuits have always been a hot research topic in the field of circuit testing. Although many important results have been achieved, the study of nonlinear systems is still in the stage of approximation process and numerical calculation. System modeling analysis, diagnosis theory, and its practicality all await further study. The nature of fault diagnosis of nonlinear analog circuits is pattern recognition, and the key to diagnosis is to construct fault features that reflect the nature of the tested circuit. Different approaches can be used to describe nonlinear analog circuits, including the Volterra series method, Wiener series method, and wavelet analysis-based feature extraction method of recent years.The Volterra series can more accurately establish common fault models of nonlinear analog circuits. However, owing to the problem of dimensionality involved in using the Volterra series, in practice, only the low-order finite term can be taken. Furthermore, for some nonanalytic nonlinear systems, Volterra series expansion cannot be used. In addition, the terms in the Volterra series are not orthogonal to each other, making it difficult to extract the fault features of the circuit to be diagnosed, which reduces the effect of the fault diagnosis [1–3].Lin et al. [4] solved this problem using the Wiener series. Gaussian white noise excitation is applied to a nonlinear system with unknown characteristic parameters, and an orthogonal functional series expansion is calculated for its response in order to obtain a set of functions that characterize system features. A time-correlation function is used to estimate the correlation function of each order to solve the Wiener kernels of the system, and a BP neural network is used to realize fault diagnosis. However, this method must be improved to obtain the Wiener kernels of the tested circuit more accurately and more efficiently [5].The feature extraction theory based on wavelet analysis is the most utilized theory in the fault diagnosis of analog circuits. For the fault diagnosis of analog circuits, the biggest advantage of wavelet analysis and artificial neural networks is that the time-domain localized features are excellent, high frequencies and low frequencies are separated, the system is insensitive to circuit noise, and the noise removal ability is satisfactory, especially when there are many nonlinear components in the tested analog circuit [6–8]. When the ambient noise is particularly high, the frequency separation feature of wavelet analysis can function effectively, greatly reducing the difficulty of extracting circuit fault features. However, this method also requires a great amount of experimental data in order to form an accurate classification of the fault phenomena [9–11].Chaotic dynamics and fractal theory developed in recent years provide a good direction for the study of nonlinear systems. Fractal theory is not only at the leading edge and an important branch of nonlinear science, but also an emerging transdisciplinary science. The use of the fractal method can accurately extract characteristics of the nonlinear signal. The majority of researchers use a single fractal to represent the overall irregular performance of signals. Zhang et al. [12] used the wavelet fractal method to perform a fault diagnosis of a three-phase bridge rectifier circuit. Mao et al. [13] also achieved satisfactory results by using analog circuit fault diagnosis based on a support vector machine (SVM) and fractal characteristics.Peng et al. [14] proposed a detrended fluctuation analysis (DFA) that analyzes power-law functions; this analysis can effectively characterize the long-range correlation of the nonlinear time series of complex dynamic systems. However, the DFA method uses only a second-order wave function to analyze the scaling behavior of a sequence of events, obtaining only a single scaling exponent. Therefore, this method is only suitable for nonstationary time series with single scaling behavior (also known as a single fractal), which lacks a description of local singular characteristics. “Multifractal” not only describes the irregularity of signals as a whole but also more finely reflects the description of the local characteristics of the signal. Multifractal depicts the local scaling characteristics of the fractal subsets of different scales and scaling exponents distributed over the subsets.Kantelhardt et al. [15] not only extended a one-dimensional fractal DFA but also proposed an MF-DFA method by combining the one-dimensional fractal DFA with the standard partition function-based multifractal formalism. This method can effectively eliminate the interference trends and fluctuation scale of time series measurement, thus effectively overcoming the defects in DFA methods when analyzing multiscale multifractal nonstationary time series. In addition, this method can more accurately estimate the multifractal spectrum.This paper introduces the MF-DFA method into nonlinear analog circuits and applies the MF-DFA method to analyze the collected signals for fault diagnosis based on the characteristic parameters of the multifractal spectrum of the fault signals.


## 2. Estimating Multifractal Spectrum Using the MF-DFA Method

Assume that the time series to be detected has a signal lengthN{xt,t=1,2,…,N}. The DFA method is conducted as follows.Step 1.
Calculate the average value of the original time series. Sum the difference of the terms between the given time seriesxt and the average value, producing a new contour sequence y(i):(1)yi=∑t=1ixt-x-i=1,2,…,N,where x-=(1/N)∑t=1Nxt and x- is the average value of the time sequence {xt,t=1,2,…,N}.Step 2.
Divide the new sequencey(i) obtained from the transform into m nonoverlapping subintervals of equal length s so that m=int(N/s), which is a rounding process for N/s. Because N may not be a multiple of s, to avoid the loss of information from the original time series, the new sequence y(i) is also divided in the same way, so that all figures in series y(i) are in the calculation. Finally, the 2m subintervals of equal set length s can be obtained.Step 3.
Remove the fluctuation trend of each small section using the least squares method. Assume thatyv(i) is the mth order of the least squares fitting polynomial of the v(v=1,2,…,2m)th small interval. The function is the kth-order polynomial fitting of the vth subinterval for s observations, and the fitting equation obtained is as follows:(2)yvi=a0+a1i+a2i2+⋯+akiki=1,2,…,s;k=1,2,…
In the formula,ak is the fitting polynomial, and k is the order of the fitting polynomial.Step 4.
Eliminate the trend within each subinterval, and calculate the mean square errorFv2s.
Whenv=1,2,…,Ns, there is(3)Fv2s=1s∑i=1sYv-1s+i-y0i2.
Whenv=Ns+1,Ns+2,…,2Ns, there is(4)Fv2s=1s∑i=1sYN-v-Nss+i-yvi2.Step 5.
For the2Ns detrended intervals, take the average value of Fv2(s) to obtain the qth wave function: (5)Fqs=12N∑v=12NsFv2sq/21/qq≠0,Fqs=exp14N∑v=12Nsln⁡Fv2sq=0.
Whenq=2, MF-DFA degenerates into DFA; different q values describe the impact on Fq(s) from different degrees of fluctuation. The mean square error Fv2(s) can eliminate the trend of each small section using the polynomial fitting of each subinterval, which is more beneficial for the singularity recognition of local fluctuation.Step 6.
By analyzing the relationship in the double logarithmic plotFq(s)~sh(q), the scaling exponent of the wave function h(q) can be determined, so a power-law relationship exists:(6)Fqs∝shq.
Eachs corresponds to a wave function Fq(s). The wave function value Fq(s) is an increasing function of s. Draw the relationship plot of lg⁡Fq(s) with respect to lg⁡s, and the slope obtained from the linear fitting of the least squares method is h(q), also known as the generalized Hurst exponent. When h(q) is constant, the time series appears as a single fractal. When h(q)=0.5, the time series appears either without a correlation or with a short correlation. When h(q) has a nonlinear relationship with q, the time series appears to have multifractal characteristics.Step 7.
Between the generalized Hurst exponenth(q) from the discrete time series and the mass index τ(q) from classic multifractal theory, there exists(7)τq=qhq-1.
In this way, the relationship between the generalized Hurst exponenth(q) and the scaling index τ(q) from the partition function-based method can be obtained. With the Legendre transform from statistical physics, (8)αq=dτqdq,fα=qαq-τq,the relationship between the multifractal spectral index α, f(α), and the generalized Hurst exponent h(q) is obtained:(9)αq=hq+qdhqdq,fα=qαq-hq+1.


## 3. Feature Extraction of Multifractal Spectrum

The multifractal spectrum obtained from the MF-DFA method is a set of parameters that can finely portray the dynamic behavior of a multifractal time series. This paper will adopt six types of commonly used multifractal spectrum parameters as the fault characteristic values of soft fault signals of a nonlinear analog circuit. The parameters of multifractal singularity spectrumα~f(α) are as follows:(1) The width of the multifractal singularity spectrumΔα=αmax-αmin: αmin corresponds to the minimum probability subset, αmax corresponds to the maximum probability subset, and Δα represents the multifractal strength of a time series, a greater value of which would mean a more uneven distribution of probability measure of the time series and a more intense change in the signals. Δα can describe the degree of uneven distribution of circuit detection signals.(2) The dimensional differenceΔf=f(αmin)-f(αmax): f(αmin) represents the fractal dimension of the maximum fluctuation point set, and f(αmax) represents the fractal dimension of the minimum fluctuation point set. Δf reflects the frequency ratio of the maximum fluctuation to the minimum fluctuation. The value of Δf reflects the proportion of the number of most stable subsets and most volatile subsets in the fluctuation in the circuit monitor signals.


## 4.  Experimental Simulation and Analysis

The MF-DFA method can be used to evaluate the fault feature extraction based on six parameters of the multifractal spectrum. To assess the performance of MF-DFA for soft fault signal analysis of nonlinear analog circuits, an experimental fault diagnosis simulation of chaotic circuits previously described in the literature [16] is conducted. The structure and parameters of the circuit are shown in Figure 1.Figure 1
Chaotic circuits system.In this experiment, the selected excitation frequency is 0.155159 Hz, the amplitude is 0.7414148 V, and the two-dimensional phase plane of node voltagesU1 and U2 at this time is shown in Figure 2.Figure 2
Two-dimensional phase plane of node voltagesU1 and U2.The normal tolerance of resistor components and capacitor components is 5% of their nominal value. In this test, a drift of plus or minus 30% and 50% of the nominal value is set as a fault condition. Only the signals at output nodeU1 are collected and analyzed, while the other measureable nodes are not included in this discussion. The fault mode setting is shown in Table 1. A total of 20 types of parametric fault combinations are tested in the experiment, in which the fault state corresponding to each fault uses a representation of “(element 1 AND element 2) ± 30%,” meaning that each element in parentheses has a drift of either positive or negative 30%. For instance, “(R1 AND R2) + 30%” indicates a fault type in which resistors R1 and R2 each have a 30% drift from their nominal values. “ROAi” indicates equivalent input impedance of amplifier. PSpice 9.1 software is used to simulate the circuits in a variety of fault states and normal states 50 times each and to collect 2000 voltage signals in the [0,30]-s time period from output node U1 as the original features of each fault state.Table 1
Fault mode setting.
Fault codes
Fault states
Fault codes
Fault states
0
Normal
10
R
11
+
30
%
1
R
1
+
30
%
11
C
1
+
30
%
2
R
1
-
30
%
12
C
1
-
30
%
3
R
2
+
30
%
13
C
2
+
30
%
4
R
2
-
30
%
14
C
2
-
30
%
5
R
4
-
30
%
15
R
5
-
50
%
6
R
6
+
30
%
16
R
9
+
50
%
7
R
7
-
30
%
17
R
O
A
1
+
50
%
8
R
8
+
30
%
18
R
O
A
3
-
50
%
9
R
10
-
30
%
19
R
O
A
5
+
50
%In this analysis, signals from faults 12 and 16 are used as examples. The signal time series waveform collected at nodeU1 is shown in Figure 3. As can be seen from the figure, the time series waveforms from fault 12 and fault 16 are fairly close to each other, and it is difficult to distinguish fault 12 and fault 16 from the time series plot. The MF-DFA is used to process these signals.Figure 3
Time series waveforms from fault 12 and fault 16.Figure4 shows the function diagram of h(q)~q of two types of faults. It can be derived from the figure that in both types of faults there is a nonlinear relationship between h(q) and q, indicating that both faults have multifractal characteristics. With the q value fixed, the relationship between the detrended fluctuation logarithmic function of two faults and different logarithmic scale lg⁡Fq(s)~lg⁡s is shown in Figures 5 and 6.Figure 4
h
(
q
)~q of fault 12 and fault 16.Figure 5
lg
⁡
F
q
(
s
)~lgs of fault 12.Figure 6
lg
⁡
F
q
(
s
)~lgs of fault 16.The Legendre transform can be used to obtain the numerical solution ofα and f(α). The multifractal spectrum f(α) of two types of faults is shown in Figure 7. Six parameter values [αmin,αmax,Δα,f(αmin),f(αmax), and Δf] of fault 12 and fault 16 can be obtained from the figure. Using these six commonly used multifractal spectrum parameters as fault characteristic values of soft fault signals of the nonlinear analog circuit, the features of circuits under different states can be easily distinguished.Figure 7
f
(
α
)~α of fault 12 and fault 16.Using the MF-DFA for collected data from 20 types of faults, different multifractal singularity spectrum graphs can be obtained (hereby omitted), thereby extracting different fault characteristic parameters from different multifractal singularity spectrum graphs. Eventually, this constitutes 20 groups of fault characteristic parameters, as shown in Table2, where different fault conditions have different multifractal characteristic parameters. Using the MF-DFA, we can effectively extract the fault features to distinguish different states.Table 2
Multifractal characteristic parameters of fault signals.
Fault codes
α
m
i
n
α
m
a
x
Δ
α
f
(
α
m
i
n
)
f
(
α
m
a
x
)
Δ
f
0
0.651
1.726
1.075
0.227
0.019
0.208
1
0.617
1.533
0.916
0.351
−0.058
0.409
2
0.669
1.835
1.166
0.272
−0.121
0.393
3
0.626
1.525
0.899
0.245
−0.083
0.328
4
0.630
1.517
0.887
0.236
−0.117
0.353
5
0.704
2.138
1.434
0.265
−0.126
0.391
6
0.688
1.841
1.153
0.263
−0.113
0.376
7
0.713
1.926
1.213
0.300
0.076
0.224
8
0.666
1.701
1.035
0.473
−0.143
0.616
9
0.683
1.714
1.031
0.384
−0.161
0.545
10
0.721
1.644
0.923
0.359
−0.157
0.516
11
0.740
1.528
0.788
0.333
0.101
0.232
12
0.642
1.532
0.890
0.363
0.056
0.307
13
0.631
1.927
1.296
0.244
−0.153
0.397
14
0.680
1.707
1.027
0.259
−0.150
0.409
15
0.704
1.817
1.113
0.323
0.097
0.226
16
0.667
1.423
0.756
0.310
0.116
0.194
17
0.637
1.662
1.025
0.214
0.017
0.197
18
0.671
1.804
1.133
0.229
0.025
0.204
19
0.634
1.701
1.067
0.210
0.015
0.195To verify the effectiveness of the approach, the extracted multifractal feature set was inputted into the support vector machine for the experiment of fault classification. Among the 50 feature samples for each fault mode, 30 feature samples were randomly selected to train the support vector machine; the rest of the feature samples were used to test the support vector machine. The accuracy of the fault diagnosis results was as high as 93%, thus proving the effectiveness of the multifractal analysis for soft fault feature extraction of nonlinear analog circuits.


## 5. Conclusion

The fractal theory can effectively analyze the inherent nature and characteristics of nonlinear signals and quantitatively describe the irregularities of fractal characteristics and signals of nonlinear systems. The soft fault signals of an analog circuit are nonstationary and nonlinear. This paper adopts MF-DFA to analyze the soft fault signals of a nonlinear analog circuit. Based on the multifractal characteristic parameters of different fault modes, this method is verified by a simulation experiment using a chaotic circuit system. The multifractal feature parameters are inputted to the support vector machine for fault detection. The results confirmed that the method can effectively extract fault features and provide favorable conditions for the realization of analog circuit fault detection.



---

*Source: 7305702-2016-05-10.xml*