# Fault Diagnosis of Nonlinear Analog Circuit Based on Generalized Frequency Response Function and LSSVM Classifier Fusion

Authors: Jialiang Zhang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

For fault diagnosis of nonlinear analog circuit, a novel method based on generalized frequency response function (GFRF) and least square support vector machine (LSSVM) classifier fusion is presented. The sinusoidal signal is used as the input of analog circuit, and then, the generalized frequency response functions are estimated directly by the time-domain formulations. The discrete Fourier transform of measurement data is avoided. After obtaining the generalized frequency response functions, the amplitudes of the GFRFs are chosen as the fault feature parameters. A classifier fusion algorithm based on least square support vector machine (LSSVM) is used for fault identification. Two LSSVM multifault classifiers with different kernel functions are constructed as subclassifiers. Fault diagnosis experiments of resistor-capacitance (RC) circuit and Sallen Key filter are carried out, respectively. The results show that the estimated GFRFs of the circuit are accurate, and the fault diagnosis method can get high recognition rate.


---

## Body

## 1. Introduction

Analog circuits are widely used in various electronic systems and play a very important role. When the fault occurred, they will affect the performance of the relevant electronic systems and even make the systems completely invalid. Therefore, it is of great significance to research the fault diagnosis of analog circuits. In order to study the problem of fault diagnosis of analog circuits, a variety of diagnosis methods have been proposed [1–3]. A wavelet neural network diagnosis model is presented to improve the diagnosis capability of soft fault for analog circuit, and the adaptive unscented Kalman filter is used to optimize the parameters of wavelet neural network [4]. The fuzzy classifier is constructed to diagnose the faults of analog electronic circuit, and the simulation before test approach is used to obtain signature parameters by analyzing the frequency response [5]. By using kurtosis and entropy as feature information, the diagnosis of analog circuit is studied based on artificial neural network [6]. A modular fault diagnosis system for analog electronic circuits is designed based on wavelet transform and neural network [7]. The features are obtained by wavelet fractal analysis and kernel principal component analysis, and novel linear ridgelet network approach is used for analog circuit fault diagnosis [8]. The relationship between the input and output of analog circuit is very complex, which is often nonlinear. There are still many problems in fault diagnosis of nonlinear analog circuits, especially in feature information acquisition and fault identification.Volterra series is an important mathematical tool for theoretical analysis and practical application of nonlinear systems. Volterra frequency-domain kernel, also known as generalized frequency response function, is an extension of frequency response function of linear system in nonlinear system. It provides a description method of frequency-domain characteristics of nonlinear system, which can reflect the inherent characteristics of the system [9–14]. The generalized random signal is used as the input, and the first two-order GFRFs of the nonlinear system are solved according to the input spectrum and the output spectrum [15]. The measurement of Volterra frequency-domain kernel for nonlinear dynamic 3×3 multi-input multi-output system is studied, and the characteristic of Volterra frequency-domain kernel is analyzed [16]. In [17], a new complex orthogonal least square algorithm is used to truncate the Volterra series and estimate the generalized frequency response function. A novel method for high-order frequency response function measurement of nonlinear system is presented based on cross-correlation technology, and the system input is sinusoidal signal [18]. For the prediction of nonlinear aeroelastic vibration response and flutter, a frequency component identification method based on correlation analysis is used to estimate the GFRFs [19]. In [20], a two-step strategy is used to estimate the generalized frequency response function of the nonlinear controlled plant in the closed-loop state. For the above estimation methods, the input and output data need to be Fourier-transformed in the process of solving GFRFs. An estimation method based on time-domain data is introduced to calculate GFRFs of single-input single-output nonlinear systems [21]. It does not need Fourier transform for time-domain data and reduces the computational complexity. The feature information in frequency domain is obtained by the generalized frequency response function, which can be used for fault diagnosis of nonlinear system. In [22], a new fault diagnosis method for nonlinear analog circuit based on simplified estimation of GFRFs and hybrid-kernel least squares support vector machine is proposed.After obtaining the frequency-domain characteristic information of nonlinear system by using the generalized frequency response function, it is very important to identify the fault accurately according to the data. As an effective machine learning method, support vector machine can be used to study classification problems [23, 24]. It can realize small sample learning and has good generalization ability. It has been used in fault diagnosis of many systems [25–28]. In [29], for the real-time fault diagnosis of high-voltage circuit breakers, fault detection is carried out by adaptive KPCA, and then fault recognition is completed by SVM classifier. The mechanical fault diagnosis of high voltage circuit breaker is studied by combining genetic algorithm with support vector machine [30]. The wavelet theory is used to obtain feature information, and then the support vector machine is used to detect bearing fault of asynchronous motor [31]. The state of asynchronous motor is monitored based on Hilbert-Park transform and support vector machine [32]. By combining support vector machine and multisensor data fusion method, the fault diagnosis of the motor is studied [33]. In [34], the statistical features are classified by support vector machine, and the fault diagnosis of the hydraulic brake system is carried out. In [35], a modified support vector machine classifier is used to analyze the fault of high voltage circuit breaker based on the waveform of coil current and travel waveforms. A two-stage scheme based on support vector machine is presented for the fault detection and classification of electrical pitch drives in offshore wind farms [36]. The least squares support vector machine is proposed by changing the risk function of support vector machine [37]. The training of LSSVM only needs to solve a system of linear equations, which reduces the computational complexity.In this paper, a new fault diagnosis method of nonlinear analog circuit is proposed based on generalized frequency response function and LSSVM classifier fusion algorithm. In order to simplify calculation, the sinusoidal signal is used as the input, and the time-domain measurements are used to directly estimate the GFRFs. After obtaining the generalized frequency response functions, their amplitudes of GFRFs are chosen as the feature parameters of analog circuits. Two multifault classifiers are constructed as subclassifier using least square support vector machine, where the kernel functions are polynomial function and Gaussian radial basis function, respectively. By using a classifier fusion algorithm based on the binary LSSVM classification function of the subclassifiers, the diagnosis result can be derived. Two experiments for the fault diagnosis of RC circuit and Sallen Key filter verified the effectiveness of the proposed method.


## 2. Simplified Estimation of GFRFs Using Time-Domain Measurements

The Volterra series is an extension of Taylor series and is widely used to describe and analyze nonlinear systems. A single-input single-output nonlinear system can be expressed as a Volterra series:(1)yt=∑n=1∞ynt,(2)ynt=∫−∞∞∫−∞∞⋯∫−∞∞hnτ1,τ2,…,τn∏i=1nut−τidτ1,…,dτn,where ut∈R is the input, yt∈R is the output, and hnτ1,τ2,…,τn is the nth order Volterra kernel.The multidimensional Fourier transform of the nth order Volterra kernelhnτ1,τ2,…,τn is given by(3)Hnjω1,jω2,…,jωn=∫−∞∞∫−∞∞⋯∫−∞∞hnτ1,τ2,…,τne−jω1τ1+ω2τ2+⋯+ωnτndτ1,…,dτn,where Hnjω1,jω2,…,jωn represents the nth-order Volterra frequency-domain kernel, also known as the nth-order generalized frequency response function, which is a mapping form of the nth-order Volterra kernel in the frequency domain.The frequency-domain output of a single-input single-output nonlinear system can be expressed as a generalized frequency response function:(4)Yjω=∑n=1∞Ynjω,Ynjω=12πn−1∫−∞∞⋯∫−∞∞Hnjω−jω1−⋯−jωn−1,jω1,…,jωn−1⋅Ujω−jω1−⋯−jωn−1Ujω1,…,Ujωn−1dω1,…,dωn−1,where Yjω is the frequency-domain output of the system, Ynjω is the nth order frequency-domain output of the system, and U⋅ is the input spectrum.The discrete form of the single-input single-output nonlinear system described by equations (1) and (2) can be expressed as(5)yk=∑n=1∞ynk,ynk=∑n=1∞∑−∞∞⋯∑−∞∞hnτ1,τ2,⋯,τn∏i=1nut−τi,k∈Z.If a weakly nonlinear system can be described by the first-order Volterra series, the frequency-domain form of the discrete Volterra model of the system is given by(6)Yω=H1ωUω+∑p+q=ωH2p,qUpUq,where H1ω is the first-order generalized frequency response function, H2p,q is the second order generalized frequency response function, and U⋅ is the input spectrum.Let the input be a sinusoidal signalut=Acosωt, and then, the output of the nonlinear system can be expressed as [21](7)yk=AReH1ωejωk+2A22ReH2ω,ωej2ωk+2A22ReH2ω,−ω+ek,where Re⋅ is the real part of a complex number, and ek is the zero mean Gaussian white noise.Define(8)H1ω=R1ω+jI1ω,(9)H2ω,ω=R2ω,ω+jI2ω,ω,(10)H2ω,−ω=R0,where R1ω and I1ω represent the real part and imaginary part of H1ω, respectively, R2ω,ω and I2ω,ω represent the real part and imaginary part of H2ω,ω, respectively, and R0 is a constant.From the equations (7)∼(10), the following can be obtained:(11)yk=AR1ωcoskω+AI1ω−sinkω+2A22R2ω,ωcos2kω−2A22I2ω,ωsin2kω+2A22R0+ek.Fork=1,2,…,N, equation (11) can be written as the matrix form:(12)Y=Xθ+E,where Y=y1,y2,…,yNT, θ=R′,I′,R2,I2,R0T, E=e1,e2,…,eNT, and(13)X=Acosω−Asinω2A22cos2ω−2A22sin2ω2A22Acos2ω−Asin2ω2A22cos4ω−2A22sin4ω2A22⋮AcosNω−AsinNω2A22cos2Nω−2A22sin2Nω2A22.Equation (12) can be solved by using standard least squares principle:(14)θ^=XTX−1XTY,where θ^=R^1,I^1,R^2,I^2,R^0T.The generalized frequency response function can reflect the essential characteristics of nonlinear systems. For a nonlinear system, the GFRF is different in different fault modes. Therefore, fault diagnosis of nonlinear analog circuits can be realized by using GFRF to obtain characteristic information. After estimating the generalized frequency response functions of the nonlinear analog circuit by equation (14), their amplitudes are chosen as the feature parameters for fault diagnosis in this paper.


## 3. Multifault Classifier Design Based on LSSVM Classifier Fusion

Support vector machine is a machine learning method developed on the basis of statistical learning theory. It uses the principle of structural risk minimization to improve the generalization ability, which can solve the problems of small samples and local minima and can effectively classify data. However, the support vector machine needs to solve the quadratic programming problem so that the calculation is very complicated. Least squares support vector machine is a deformed support vector machine algorithm. The training of LSSVM only requires solving a system of linear equations, which is easy to solve and has high training efficiency.Let the training sample data beS1:xi,yii=1M, where xi∈Rp is the input space sample, yi∈−1,1 is the class label, and M is the sample size. The binary classification problem of least squares support vector machine can be described as(15)min12ωΤω+12C∑i=1Mξi2s.t.yiωΤφxi+b=1−ξi,i=1,2,…,M,where ω is the weight vector of the classification hyperplane, C>0 is the penalty factor, ξi is the relaxation variable, φ⋅ is the nonlinear mapping, and b is the classification threshold.Define the Lagrange function as(16)Lω,b,ξi,αi=12ωΤω+12C∑i=1Mξi2−∑i=1MαiyiωΤφxi+b−1+ξi,where αi≥0 is the Lagrange multiplier.Letting the partial derivatives of the Lagrange function forω, b, ξi and αi be zero, the constraint optimization problem of LSSVM can be transformed into a system of linear equations:(17)0YTYH⋅bα=01,where Y=y1,y2,⋯,yMT, H=Ω+c−1⋅I, Ω is the M-dimensional symmetric matrix, Ωi,j=yiyjKxi,xj, K⋅,⋅ is the kernel function, I is the M-dimensional unit matrix, α=α1,α2,⋯,αMT is the Lagrange multiplier vector, and 1=1,1,…,1T.The decision function is given by(18)yx=sign∑i=1MαiyiKx,xi+b,where x is the sample vector to be classified.For the classification of nonlinear data, the SVM classifier maps the feature data to a high-dimensional feature space using a nonlinear function and then completes the classification. Therefore, when using the least squares support vector machine to construct a classifier, the choice of kernel function is very important and will directly affect the classification result. Commonly used kernel functions include polynomial kernel function and Gaussian radial basis kernel function. When traditional LSSVM deals with the multiclassification problem, one type of function is selected as the kernel function, and one classifier is constructed for classification. Because different kernel functions have different characteristics, their classification performance for the same feature data is also different.In order to improve the fault recognition rate, a LSSVM classifier fusion method is adopted to identify faults in nonlinear analog circuits. Polynomial kernel function and Gaussian radial basis kernel function are used to construct two LSSVM classifiers as subclassifiers. Then, the two subclassifiers are fused to obtain the diagnosis result. The polynomial kernel function is a global kernel function, and the Gaussian radial basis kernel function is a local kernel function. Therefore, this LSSVM fusion algorithm can take into account the global characteristics and local characteristics of the feature data obtained through the generalized frequency response functions of the nonlinear analog circuits.The structures of the designed LSSVM subclassifiers are “one-against-one,” and the voting method is used for fault decisions of subclassifiers. The structure of nonlinear circuit fault diagnosis based on LSSVM fusion is shown in Figure1, where N is the number of categories of data to be classified.Figure 1
The structure of nonlinear circuit fault diagnosis based on LSSVM fusion.The fusion strategy of support vector machine has a great impact on the accuracy of the classifier. Therefore, choosing an appropriate fusion algorithm is very important for fault identification.According to decision function (18) of the binary SVM, define(19)y′x=αiyiKx,xi+b.The functiony′x represents the distance from sample x to the optimal hyperplane, which can indirectly reflect the degree of support for decision results. In this paper, the weight coefficients are obtained on the basis of the function y′x of each binary SVM in the subclassifier, and then the weighted fusion is performed to obtain the fault identification result.Assume that the fault modes of the nonlinear analog circuit aref1,f2,…,fN, N is the number of modes, the obtained GFRFs amplitude data of the nonlinear analog circuit is x, and the decision results of the two LSSVM subclassifiers are D1 and D2. Define the weight of the LSSVM subclassifier as(20)wj=∑k=1Vj∑i=1MαjkiyjkiKjkx,xjki+bjk,where j=1,2, M is the training sample capacity for each binary LSSVM classifier in subclassifier j, Vj is the number of binary LSSVM classifiers with a decision result of Dj, αjki and yjki are the Lagrangian multiplier and class label of the kth binary LSSVM classifier in subclassifier j, and Kjkx,xjki and bjk are the kernel function and the classification threshold of the K-th binary LSSVM classifier in the subclassifier j.By fusing the decision results of the two subclassifiers, the fault identification result is given by(21)D=argmaxfl∈Q∑j=12wjδj,l,(22)δj,l=1,Dj=fl,0,Dj≠fl,where Q⊆f1,f2,…,fN is the set of fault modes formed by the decision results of the two LSSVM subclassifiers.The result of subclassifier is given by(23)Dj=argmaxfk∑m=1pAjkm,k=1,2,…,N,(24)Ajkm=1Djm=fk0Djm≠fk,m=1,2,…,p,where Djm is the decision result of the mth binary LSSVM classifier of subclassifier j, and p=NN−1/2.After obtaining the frequency-domain feature information of the nonlinear analog circuit based on the generalized frequency response functions, the fusion algorithm of least square support vector machine is used for fault diagnosis. The fault diagnosis steps are as follows:(i)
Step 1. Input the frequency-domain feature vector x of the nonlinear analog circuit into subclassifier 1 and subclassifier 2, respectively, and solve the value of the αjkiyjkiKjkx,xjki+bjk, k=1,…,NN−1/2 of the binary LSSVM classifier in the subclassifiers(ii)
Step 2. For each LSSVM subclassifier, calculate the decision result by equation (18), and then make a voting based on all of the results(iii)
Step 3. According to the classification function of the binary LSSVMs in the subclassifiers, calculate the weights of the LSSVM subclassifiers by using equation (20)(iv)
Step 4. According to the decision results D1 and D2 of the subclassifier, and the corresponding weights w1 and w2, calculate the fusion result by equations (21)∼(24)Because the fault diagnosis method of nonlinear analog circuits based on LSSVM fusion uses the classification function of binary LSSVM to obtain the weight of the subclassifier, the classification reliability of the two subclassifiers with different types of kernel functions is fully considered.


## 4. Experiment and Analysis

Two experiments of nonlinear analog circuits are used to verify the effectiveness of the diagnosis method.


### 4.1. Resistor-Capacitance (RC) Circuit

A resistor-capacitance circuit is shown in Figure2, where us is the supply voltage, R1 and R2 are the resistances, and C is the capacitance. The relationship between voltage and current of the nonlinear resistance R2 is i2=a1u+a2u2, where a1=1×10−3 and a2=5×10−4.Figure 2
A resistor-capacitance circuit.According to Kirchhoff’s voltage law and current law, we can obtain(25)uint=R1i1t+ut,i1t=i2t+i3t,where i3t=Cdut/dt.Let the power supply voltageuint be the input and let the voltage ut on R2 be the output. The nonlinear differential equation can be obtained as(26)dutdt+b1ut+b2u2t=b3uint,where b1=a1R1+1/R1C, b2=a2/C, and b3=1/R1C.The nominal values of components in the analog circuit system areC=10−4F and R1=103Ω. The simulation experiment has been performed using MATLAB R2014a. The CPU clock speed of the computer is 2.3 GHz, and the main memory is 8 GB.Let the input signal beuint=cos2πt. The sampling frequency is 100 Hz, and sampling length is 10 s. The output spectrum of the nonlinear circuit is shown in Figure 3. There are significant peaks at frequencies 1 Hz and 2 Hz, which are 0.48 and 0.03, respectively. The output of analog circuit has obvious second harmonics so that the circuit has obvious nonlinear characteristics. Because the harmonic amplitudes above the third harmonic are approximately equal to 0, the analog circuit can be described by the first two-order generalized frequency response functions.Figure 3
The output spectrum of analog circuit.For the nonlinear circuit system described by Equation (26), the actual value of GFRFs can be derived [38]:(27)H1jω=b3jω+b1,H2jω1,jω2=−b2H1jω1H1jω2jω1+ω2+b1.When the output noise is Gaussian white noise with signal-to-noise ratios (SNR) of 35 dB, 30 dB, and 25 dB, respectively, the simplified estimation method is used to obtain the generalized frequency response functions of the analog circuit. The output of the nonlinear circuit is shown in Figure4. As the SNR gradually decreases, the deformation of the output becomes more and more obvious. The estimation results are shown in Table 1, where ω=2π. Under different signal-to-noise ratios, the GFRFs of nonlinear analog circuit can be accurately estimated. As the signal-to-noise ratio decreases, the estimation error gradually increases.Figure 4
The output of analog circuit: (a) SNR = 35 dB, (b) SNR = 30 dB, (c) SNR = 25 dB.
(a)(b)(c)Table 1
The estimation results of GFRFs.
GFRFsActual value [38]Estimated valueSNR = 35 dBSNR = 30 dBSNR = 25 dBH1ω0.455–0.143i0.462–0.150i0.459–0.152i0.459–0.153iH2ω,ω−0.019 + 0.044i−0.019 + 0.048i−0.020 + 0.048i−0.023 + 0.050iH2ω,−ω−0.057−0.059−0.060−0.060Assume that, under normal conditions, the variation ranges ofC and R1 are within 5%. If the actual value of C or R1 differs from the nominal value by 50%, the component will fail. The analog circuit contains five fault modes, and the fault description is shown in Table 2.Table 2
The fault description of analog circuit.
Fault numberFault typeParameter valueF0NormalC∈0.95×10−4F∼1.05×10−4F, R1∈0.95×103Ω∼1.05×103ΩF1C increasedC∈1.45×10−4F∼1.55×10−4F, R1∈0.95×103Ω∼1.05×103ΩF2C decreasedC∈0.45×10−4F∼0.55×10−4F, R1∈0.95×103Ω∼1.05×103ΩF3R1 increasedC∈0.95×10−4F∼1.05×10−4F, R1∈1.45×103Ω∼1.55×103ΩF4R1 decreasedC∈0.95×10−4F∼1.05×10−4F, R1∈0.45×103Ω∼0.55×103ΩThe Monte Carlo method is used for fault diagnosis simulation experiment of the analog circuit. The output noise is Gaussian white noise, and the SNR is 25 dB. 200 sets of input and output sample data of the analog circuit are collected for each fault. The amplitudes of the GFRFsH1ω, H2ω,ω and H2ω,−ω are solved by the proposed simplified estimation method. The amplitude distributions of the five fault types are shown in Figure 5.Figure 5
The amplitude distributions of the five fault types.The 100 sets of sample data are used as training samples, and the remaining 100 sets are used as test samples. The fault diagnosis of the analog circuit is carried out by the LSSVM fusion algorithm and the traditional single LSSVM multiclassification algorithm, respectively. The order of the polynomial function is 2, and the width of the Gaussian radial basis function is4×10−3. The fault diagnosis steps of the LSSVM fusion algorithm are shown in Section 3, and the results are obtained using Step 1 to Step 4.The fault recognition rate is shown in Table3. The recognition rates of the LSSVM fusion algorithm are higher than those of the traditional single LSSVM classifiers.Table 3
The fault recognition rate of RC circuit.
Fault numberFault typeRecognition rateLSSVM fusion algorithmSingle LSSVM (polynomial kernel function)Single LSSVM (Gaussian radial basis kernel function)F0Normal939186F1C increased999895F2C decreased909088F3R1 increased10010096F4R1 decreased100100100


### 4.2. Sallen Key Filter

A Sallen Key filter circuit is shown in Figure6, where R1=1kΩ, R2=3kΩ, R3=2kΩ, R4=4kΩm, R5=4kΩ, C1=5nF, C2=5nF, and the input voltage is u=cos2000πtV.Figure 6
A Sallen Key filter circuit.Assuming that the actual value ofR4 or R5 differs from the nominal value by 50%, the component will fail. There are five fault modes, normal (F0), R4 increased by 50% (F1), R4 decreased by 50% (F2), R5 increased by 50% (F3), and R5 decreased by 50% (F4). The Monte Carlo method is used for fault diagnosis experiment of the Sallen Key filter circuit, and the tolerances of the circuit components are 5%. The input and output data of the circuit are obtained using Multisim10. The diagnosis algorithm has been performed using MATLAB R2014a.The output noise is Gaussian white noise, and the SNR is 40 dB. The sampling frequency is 105 Hz, and sampling length is 10 ms. The output of the Sallen Key filter in normal state is shown in Figure7.Figure 7
The output of Sallen Key filter.200 sets of input and output sample data of the analog circuit are collected for each fault. The amplitudes of the GFRFsH1ω, H2ω,ω and H2ω,−ω are solved by the proposed simplified estimation method, where ω=2000π. The amplitude distributions of the five fault types of the Sallen Key filter are shown in Figure 8.Figure 8
The amplitude distributions of the five fault types of Sallen Key filter.The 100 sets of sample data are used as training samples, and the remaining 100 sets are used as test samples. The fault diagnosis of the Sallen Key filter is carried out by the LSSVM fusion algorithm. The order of the polynomial function is 3, and the width of the Gaussian radial basis function is2×10−3. In addition, the single LSSVM and multikernel LSSVM [39] multiclassification algorithm are used for comparative experiments.The fault recognition rate is shown in Table4. The recognition rates of the LSSVM fusion algorithm are significantly higher than those of the single LSSVM and multikernel LSSVM classifiers.Table 4
The fault recognition rate of Sallen Key filter.
Fault numberFault typeRecognition rateLSSVM fusion algorithm (%)Single LSSVM (polynomial kernel function) (%)Single LSSVM (Gaussian radial basis kernel function) (%)Multikernel LSSVM (%)F0Normal94858890F1R4 increased81697173F2R4 decreased96878988F3R5 increased97868891F4R5 decreased83727576Therefore, the proposed diagnosis method for nonlinear analog circuit based on the generalized frequency response function and LSSVM fusion algorithm has good diagnostic performance.


## 4.1. Resistor-Capacitance (RC) Circuit

A resistor-capacitance circuit is shown in Figure2, where us is the supply voltage, R1 and R2 are the resistances, and C is the capacitance. The relationship between voltage and current of the nonlinear resistance R2 is i2=a1u+a2u2, where a1=1×10−3 and a2=5×10−4.Figure 2
A resistor-capacitance circuit.According to Kirchhoff’s voltage law and current law, we can obtain(25)uint=R1i1t+ut,i1t=i2t+i3t,where i3t=Cdut/dt.Let the power supply voltageuint be the input and let the voltage ut on R2 be the output. The nonlinear differential equation can be obtained as(26)dutdt+b1ut+b2u2t=b3uint,where b1=a1R1+1/R1C, b2=a2/C, and b3=1/R1C.The nominal values of components in the analog circuit system areC=10−4F and R1=103Ω. The simulation experiment has been performed using MATLAB R2014a. The CPU clock speed of the computer is 2.3 GHz, and the main memory is 8 GB.Let the input signal beuint=cos2πt. The sampling frequency is 100 Hz, and sampling length is 10 s. The output spectrum of the nonlinear circuit is shown in Figure 3. There are significant peaks at frequencies 1 Hz and 2 Hz, which are 0.48 and 0.03, respectively. The output of analog circuit has obvious second harmonics so that the circuit has obvious nonlinear characteristics. Because the harmonic amplitudes above the third harmonic are approximately equal to 0, the analog circuit can be described by the first two-order generalized frequency response functions.Figure 3
The output spectrum of analog circuit.For the nonlinear circuit system described by Equation (26), the actual value of GFRFs can be derived [38]:(27)H1jω=b3jω+b1,H2jω1,jω2=−b2H1jω1H1jω2jω1+ω2+b1.When the output noise is Gaussian white noise with signal-to-noise ratios (SNR) of 35 dB, 30 dB, and 25 dB, respectively, the simplified estimation method is used to obtain the generalized frequency response functions of the analog circuit. The output of the nonlinear circuit is shown in Figure4. As the SNR gradually decreases, the deformation of the output becomes more and more obvious. The estimation results are shown in Table 1, where ω=2π. Under different signal-to-noise ratios, the GFRFs of nonlinear analog circuit can be accurately estimated. As the signal-to-noise ratio decreases, the estimation error gradually increases.Figure 4
The output of analog circuit: (a) SNR = 35 dB, (b) SNR = 30 dB, (c) SNR = 25 dB.
(a)(b)(c)Table 1
The estimation results of GFRFs.
GFRFsActual value [38]Estimated valueSNR = 35 dBSNR = 30 dBSNR = 25 dBH1ω0.455–0.143i0.462–0.150i0.459–0.152i0.459–0.153iH2ω,ω−0.019 + 0.044i−0.019 + 0.048i−0.020 + 0.048i−0.023 + 0.050iH2ω,−ω−0.057−0.059−0.060−0.060Assume that, under normal conditions, the variation ranges ofC and R1 are within 5%. If the actual value of C or R1 differs from the nominal value by 50%, the component will fail. The analog circuit contains five fault modes, and the fault description is shown in Table 2.Table 2
The fault description of analog circuit.
Fault numberFault typeParameter valueF0NormalC∈0.95×10−4F∼1.05×10−4F, R1∈0.95×103Ω∼1.05×103ΩF1C increasedC∈1.45×10−4F∼1.55×10−4F, R1∈0.95×103Ω∼1.05×103ΩF2C decreasedC∈0.45×10−4F∼0.55×10−4F, R1∈0.95×103Ω∼1.05×103ΩF3R1 increasedC∈0.95×10−4F∼1.05×10−4F, R1∈1.45×103Ω∼1.55×103ΩF4R1 decreasedC∈0.95×10−4F∼1.05×10−4F, R1∈0.45×103Ω∼0.55×103ΩThe Monte Carlo method is used for fault diagnosis simulation experiment of the analog circuit. The output noise is Gaussian white noise, and the SNR is 25 dB. 200 sets of input and output sample data of the analog circuit are collected for each fault. The amplitudes of the GFRFsH1ω, H2ω,ω and H2ω,−ω are solved by the proposed simplified estimation method. The amplitude distributions of the five fault types are shown in Figure 5.Figure 5
The amplitude distributions of the five fault types.The 100 sets of sample data are used as training samples, and the remaining 100 sets are used as test samples. The fault diagnosis of the analog circuit is carried out by the LSSVM fusion algorithm and the traditional single LSSVM multiclassification algorithm, respectively. The order of the polynomial function is 2, and the width of the Gaussian radial basis function is4×10−3. The fault diagnosis steps of the LSSVM fusion algorithm are shown in Section 3, and the results are obtained using Step 1 to Step 4.The fault recognition rate is shown in Table3. The recognition rates of the LSSVM fusion algorithm are higher than those of the traditional single LSSVM classifiers.Table 3
The fault recognition rate of RC circuit.
Fault numberFault typeRecognition rateLSSVM fusion algorithmSingle LSSVM (polynomial kernel function)Single LSSVM (Gaussian radial basis kernel function)F0Normal939186F1C increased999895F2C decreased909088F3R1 increased10010096F4R1 decreased100100100


## 4.2. Sallen Key Filter

A Sallen Key filter circuit is shown in Figure6, where R1=1kΩ, R2=3kΩ, R3=2kΩ, R4=4kΩm, R5=4kΩ, C1=5nF, C2=5nF, and the input voltage is u=cos2000πtV.Figure 6
A Sallen Key filter circuit.Assuming that the actual value ofR4 or R5 differs from the nominal value by 50%, the component will fail. There are five fault modes, normal (F0), R4 increased by 50% (F1), R4 decreased by 50% (F2), R5 increased by 50% (F3), and R5 decreased by 50% (F4). The Monte Carlo method is used for fault diagnosis experiment of the Sallen Key filter circuit, and the tolerances of the circuit components are 5%. The input and output data of the circuit are obtained using Multisim10. The diagnosis algorithm has been performed using MATLAB R2014a.The output noise is Gaussian white noise, and the SNR is 40 dB. The sampling frequency is 105 Hz, and sampling length is 10 ms. The output of the Sallen Key filter in normal state is shown in Figure7.Figure 7
The output of Sallen Key filter.200 sets of input and output sample data of the analog circuit are collected for each fault. The amplitudes of the GFRFsH1ω, H2ω,ω and H2ω,−ω are solved by the proposed simplified estimation method, where ω=2000π. The amplitude distributions of the five fault types of the Sallen Key filter are shown in Figure 8.Figure 8
The amplitude distributions of the five fault types of Sallen Key filter.The 100 sets of sample data are used as training samples, and the remaining 100 sets are used as test samples. The fault diagnosis of the Sallen Key filter is carried out by the LSSVM fusion algorithm. The order of the polynomial function is 3, and the width of the Gaussian radial basis function is2×10−3. In addition, the single LSSVM and multikernel LSSVM [39] multiclassification algorithm are used for comparative experiments.The fault recognition rate is shown in Table4. The recognition rates of the LSSVM fusion algorithm are significantly higher than those of the single LSSVM and multikernel LSSVM classifiers.Table 4
The fault recognition rate of Sallen Key filter.
Fault numberFault typeRecognition rateLSSVM fusion algorithm (%)Single LSSVM (polynomial kernel function) (%)Single LSSVM (Gaussian radial basis kernel function) (%)Multikernel LSSVM (%)F0Normal94858890F1R4 increased81697173F2R4 decreased96878988F3R5 increased97868891F4R5 decreased83727576Therefore, the proposed diagnosis method for nonlinear analog circuit based on the generalized frequency response function and LSSVM fusion algorithm has good diagnostic performance.


## 5. Conclusions

According to generalized frequency response function and LSSVM classifier fusion, a fault diagnosis method for nonlinear analog circuits is presented in this paper. The sinusoidal signal is used as input, and the first two-order generalized frequency response functions of the circuit system are directly estimated by the time-domain measurements. The amplitudes of the generalized frequency response functions are used as the feature parameters. In order to improve the accuracy of diagnosis, a classifier fusion algorithm based on LSSVM is presented for fault identification. The polynomial kernel function and Gaussian radial basis function are used to construct two different LSSVM classifiers, and then the decision fusion is carried out by the binary LSSVMs of these two classifiers. Two fault diagnosis experiments for RC circuit and Sallen Key filter are used to verify the effectiveness of the proposed method. The results show that the method can achieve good fault diagnosis performance.For fault diagnosis of complex circuits in real time applications, the diagnosis can be carried out by the proposed method only when the input signal is sinusoidal. Furthermore, the fault diagnosis of analog circuit under the excitation of multitone signal will be studied.



---

*Source: 8274570-2020-09-22.xml*