# Feature Extraction of Double Pulse Metal Inert Gas Welding Based on Broadband Mode Decomposition and Locality Preserving Projection

Authors: Yanfeng Peng; Zucheng Wang; Kuanfang He; Yanfei Liu; Qingxian Li; Liangjiang Liu; Xianyu Zhu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

A novel adaptive signal decomposition algorithm, broadband mode decomposition (BMD), is proposed for analyzing nonstationary broadband signals. Unavoidable error will occur when applying former time-frequency methods to broadband signals, which is caused by Gibbs phenomenon and the calculation of extrema. To overcome that problem, BMD is proposed by searching in the associative dictionary that contains both broadband and narrowband signals. The procedure of the proposed method is as follows: First, the collected datasets are analyzed by BMD and the composite multiscale fuzzy entropies (CMFEs) of the obtained effective components are calculated. Then, locality preserving projection (LPP) is applied for further feature extraction. Analysis results show BMD is more effective when drawing broadband feature from noise and BMD is adaptive for the quality monitoring of DPMIG welding.


---

## Body

## 1. Introduction

Research on signal processing methods for extracting useful feature from noised metal inert gas (MIG) welding signals has been a research hotspot in welding quality estimation [1]. In the MIG welding procedure, the energy put into the welding line is mainly affected by the electrical inputs, which will immediately impact the quality. Thus, the feature embedded in the electrical data can reflect the quality of MIG welding [2]. Nevertheless, the high frequency electric data of double pulse metal inert gas (DPMIG) are modulated by the low frequency pulse frequencies and mixed with severe noise signals. Feature information which can reflect the welding quality usually is narrow-band signals, and the pulse signals such as square signals and sawtooth signals with “sharp corners” are broadband signals as their frequency bands are infinite [3]. Therefore, in order to extract narrowband feature which contains the quality information of welding from pulse signals and noise signals, a suitable time-frequency method should be applied.Many time-frequency analyses algorithms have been used for extracting feature from noise in engineering fields [3–5]. These adaptive time-frequency analysis algorithms can be classified into two categories: FT based algorithms and non-FT based algorithms [6–8]. The FT based algorithms including the variational mode decomposition (VMD) use the calculations of frequencies to fulfill the data analysis procedure [9]. VMD algorithm is an iteration procedure using Wiener filter for denoising. VMD was introduced by Dragomiretskiy and Zosso [10] and showed the effectiveness in analyzing nonstationary industrial signal compared with other algorithms [11–13]. Non-FT based algorithms including empirical mode decomposition (EMD) separate the original data into intrinsic mode functions (IMFs) through the calculation of the extrema [14–16]. Ensemble empirical mode decomposition (EEMD) was introduced by Wu et al. to increase the denoising ability of EMD by mixing various levels of noise in the data and calculating the average of the generated IMFs [17–19]. Researches have demonstrated the advantage and validity of EEMD and VMD nearly [11, 15]; the two methods are selected to compare with the BMD algorithm.The former time-frequency analysis methods are widely applied for extract feature from complex nonstationary signals. However, problems will appear when applying to broadband data as the bandwidth is infinite. Firstly, for FT based algorithms, the substance of the methods is multiple adaptive filters. Nevertheless, the Gibbs phenomenon indicates the high frequencies of broadband signals will attenuate or vanish after smoothing [20], which may generate a range of fluctuations at the mutant points of the decomposed results. Secondly, for non-FT based algorithms, the original data is decomposed into a series of “smooth” narrowband IMFs by computing the extrema envelops utilizing interpolating functions. Thus, an inevitable error will be generated when handling the broadband signals.For extracting useful information from complicated signals containing typical broadband modes and noise components, the broadband mode decomposition (BMD) algorithm is proposed. In BMD algorithm, an associative dictionary containing common broadband data, such as sawtooth, square, and narrowband data, is constructed. Then, the sparse solution is gained by seeking in the associative dictionary applying an optimization method. The BMD method utilizes a regulated smooth operator as the optimal object. The initial phase, the frequency, and the amplitude of common broadband data are used as the optimal arguments. The artificial chemical reaction optimization algorithm (ACROA) [21], which was introduced by Alatas, is employed for the solution of the optimal problems, and the sparse separation outcomes are gained during the optimal procedure. Compared with previous signal decomposition algorithms, BMD is suitable for extracting broadband features and narrowband features from noisy complex nonstationary signals.Broadband data including sawtooth and square data are used as the input voltage and current signals of DPMIG welding. As stated above, an inevitable error will be generated while dealing with broadband electric signals of inverter power sources when using previous adaptive signal decomposition signals. Thus, an algorithm for dealing with broadband data such as the introduced BMD algorithm is required. In this paper, the collected welding input electric datasets are applied for evaluating the validity of BMD.To evaluate the quality of DPMIG welding, suitable parameters should be calculated for the representation of features. Nevertheless, it is worth noting that not all the features are useful for evaluating the quality. Thus, suitable algorithms should be applied for merging the effective information from the mixed broadband and noisy modes, as well as reducing the feature dimension. Many developed feature extraction methods containing linear discriminant analysis (LDA), principal component analysis (PCA), Laplacian eigenmaps (LE), and locally linear embedding (LLE) [22–25] have been introduced to carry out the goal. Locality preserving projections (LPP) [26] was introduced by He et al. It can produce projective linear mapping which has the essential feature of the datasets in a low-dimensional space. LPP has been demonstrated to be suitable for complex data processing in a plenty of researches [27–29]. Therefore, it is selected for fusing the feature vectors generated from the obtained MIG welding electric datasets.The remaining contents are organized as below. In Section2, the existing disadvantages in analyzing the broadband data are given out. The iteration procedure of the BMD method is given in Section 3. In Section 4, the basic theories of LPP are stated. In Section 5, the collected experimental datasets of aluminum alloy DPMIG welding input voltage and current data are tested by the introduced algorithm. The conclusion is given out at last.


## 2. Existing Problems in Broadband Signal Analysis

### 2.1. VMD Method

The Fourier progression unfolds the data into a series of sinusoidal components. The square signal is taken as an illustration of broadband data; the Fourier progression is as follows:(1)squaret=∑i=1+∞12i−1sin2i−1t.As the bandwidth of a filter is restricted, the high frequencies of broadband signals will attenuate or vanish after filtering, which may induce a range of fluctuations at the mutant points of the decomposed results. When the bandwidth of the used filter is wider, the obtained filtering results are more accurate. This is named as Gibbs phenomenon [20].VMD decomposes a dataset into a series of band limited IMFs by using multiple Wiener filters. The filter parameters including the bandwidth and center frequency are gained by searching for the variational modes. As the substance of the methods is adaptive filters, it will be influenced by the Gibbs phenomenon when processing with broadband datasets. The iteration procedure for obtaining IMFs is as follows [10]:(1)
Calculate the mathematical formula ofIMFin by using the Hilbert transform:(2)δn+jπn∗IMFin,whereδ is the Dirac distribution and denotes convolution.(2)
Compute the center frequencyωi of IMFin; transform the formula in equation (2) to the frequency domain corresponding to ωi as follows:(3)δn+jπn∗IMFine−jwin.(3)
An optimization problem is constructed as equation (4); the optimal object is the L2 norm of the gradient value of the frequency domain formula. The separated results are gained by solving the optimization problem using the alternate direction multiplier method:(4)minIMFin,wi∑i∂tδn+jπn∗IMFine−jwin2S.T.xn=∑IMFin.


### 2.2. EEMD Method

By the continual extraction of the mean lower and upper envelopes applying the interpolating method, EMD separates the original data into a series of IMFs. The iteration procedure of EMD is as follows [8]:(1)
Sethin=xn, i=1, and j=1.(2)
Interpolate the extrema ofhin using the interpolating method to gain the lower envelope and upper envelope mun and mln; then calculate hin:(5)hi+1n=hin−mun+min2.(3)
Step (2) is repeated untilhi+1n accords with the condition of the IMFs [8]. The jth IMF is obtained:(6)IMFjn=hi+1n.(4)
Setrn=xn−IMFjn; set the residue rn to be the updated xn. Steps (1)–(3) are repeated until rn is monotonic data. The separated results are expressed as follows:(7)xn=∑j=1NIMFjn+rn.EMD is suitable for processing nonstationary datasets [14]. Nevertheless, mode mixture is one of the main disadvantages existing in EMD algorithm, which may induce serious error between the useful information when the data is strongly interfered by noise. Thus, to increase the mode mixing resisting and denoising ability, EEMD is introduced through the addition of various levels of noise in the data and calculating the average of the generated IMFs [15]. The additional noise datasets are irrelevant; thus the averaging procedure of EEMD can offset part of the noise.Non-FT based algorithms including EMD and EEMD apply the interpolating method to gain the envelope of the extrema. That will induce smooth upper and lower envelopes. Because the obtained IMFs are calculated from the envelopes, the smooth separation results will also be generated when conducting with broadband data.


## 2.1. VMD Method

The Fourier progression unfolds the data into a series of sinusoidal components. The square signal is taken as an illustration of broadband data; the Fourier progression is as follows:(1)squaret=∑i=1+∞12i−1sin2i−1t.As the bandwidth of a filter is restricted, the high frequencies of broadband signals will attenuate or vanish after filtering, which may induce a range of fluctuations at the mutant points of the decomposed results. When the bandwidth of the used filter is wider, the obtained filtering results are more accurate. This is named as Gibbs phenomenon [20].VMD decomposes a dataset into a series of band limited IMFs by using multiple Wiener filters. The filter parameters including the bandwidth and center frequency are gained by searching for the variational modes. As the substance of the methods is adaptive filters, it will be influenced by the Gibbs phenomenon when processing with broadband datasets. The iteration procedure for obtaining IMFs is as follows [10]:(1)
Calculate the mathematical formula ofIMFin by using the Hilbert transform:(2)δn+jπn∗IMFin,whereδ is the Dirac distribution and denotes convolution.(2)
Compute the center frequencyωi of IMFin; transform the formula in equation (2) to the frequency domain corresponding to ωi as follows:(3)δn+jπn∗IMFine−jwin.(3)
An optimization problem is constructed as equation (4); the optimal object is the L2 norm of the gradient value of the frequency domain formula. The separated results are gained by solving the optimization problem using the alternate direction multiplier method:(4)minIMFin,wi∑i∂tδn+jπn∗IMFine−jwin2S.T.xn=∑IMFin.


## 2.2. EEMD Method

By the continual extraction of the mean lower and upper envelopes applying the interpolating method, EMD separates the original data into a series of IMFs. The iteration procedure of EMD is as follows [8]:(1)
Sethin=xn, i=1, and j=1.(2)
Interpolate the extrema ofhin using the interpolating method to gain the lower envelope and upper envelope mun and mln; then calculate hin:(5)hi+1n=hin−mun+min2.(3)
Step (2) is repeated untilhi+1n accords with the condition of the IMFs [8]. The jth IMF is obtained:(6)IMFjn=hi+1n.(4)
Setrn=xn−IMFjn; set the residue rn to be the updated xn. Steps (1)–(3) are repeated until rn is monotonic data. The separated results are expressed as follows:(7)xn=∑j=1NIMFjn+rn.EMD is suitable for processing nonstationary datasets [14]. Nevertheless, mode mixture is one of the main disadvantages existing in EMD algorithm, which may induce serious error between the useful information when the data is strongly interfered by noise. Thus, to increase the mode mixing resisting and denoising ability, EEMD is introduced through the addition of various levels of noise in the data and calculating the average of the generated IMFs [15]. The additional noise datasets are irrelevant; thus the averaging procedure of EEMD can offset part of the noise.Non-FT based algorithms including EMD and EEMD apply the interpolating method to gain the envelope of the extrema. That will induce smooth upper and lower envelopes. Because the obtained IMFs are calculated from the envelopes, the smooth separation results will also be generated when conducting with broadband data.


## 3. BMD Algorithm

### 3.1. Dictionary

To decompose both the broadband modes and the AM-FM narrowband modes from noisy signals, an associative dictionary containing common broadband data, such as sawtooth, square, and narrowband data, is constructed at first, in which the narrowband dictionary was defined in [30]. The analysis results are gained by seeking in the 3 dictionaries defined as follows:(8)Dic1=A1squareω1n+θ1,D1,Dic2=A2sawtoochω2n+θ2,D2,Dic3=A3ncosω3n+θ3n.θ, ω, and A represent the original phase, the angular frequency, and the amplitude, respectively. In Dic3, the maximal angular frequency of A3n is restricted to be much less than ω3. θ3n is restricted to be a slowly changing signal to guarantee the elements in Dic3 are narrowband [30].


### 3.2. BMD Iteration Procedure

The BMD algorithm is introduced to solve the existing disadvantages of extracting broadband features when using previous signal decomposition algorithms, including the Gibbs phenomenon in FT based algorithms and the interpolating of extrema in non-FT based algorithms. In BMD, the associative dictionaries are established at first. Then, the sparse decomposition results are gained by seeking in the dictionaries. BMD utilizes a normalized smoothness objective function as the optimal object. The parameters of the typical broadband datasets are used as the optimal parameters. The ACROA is used for the solution of the optimal problems [21]. In the end, both the broadband IMFs and the narrowband IMFs can be separated from noisy signals by using the optimal procedure. The flow chart of BMD is illustrated in Figure 1. The main iteration steps are designed as follows:(1)
r0n is set to equal xn.(2)
Separate the direct-current (DC) signal:(9)IMF0n=ifftr^01,r1n=xn−IMF0n.r^0k is the FT of r0n; ifftr^01 is the IFT of r^01.(3)
i is set to equal 1; the optimal problem P1 is constructed as follows:(10)P1:MinimizeT1A1,ω1,θ1,D1,T2A2,ω2,θ2,D2,T3A3n,ω3n,θ3nS.T.xn=∑i=0NIMFin+resn,IMFij∈Dicj,Tj=D2IMFijn22+λD2xn−IMFijn22.Tj is a normalized smooth operator to gain the sparse solution of r0n [30], D2 is the quadratic differential form operator, and λ>0 is always set to equal 1.(4)
Solve L1 by using the ACROA to gain the optimalTjj=1,2,3, find the minimum Tj, and choose the optimal IMFni of Tj. The procedure of the method is introduced in as follows.(5)
Updaterin: ri+1n=rin−IMFin.(6)
If the terminal condition in equation (11) is achieved, terminate the iteration procedure; else, get back to step 3. ε is a preset parameter which is normally set to be 1×10−6:(11)IMFin−IMFi−1n22IMFi−1n22≤ε.Figure 1
The iteration procedure of BMD.


### 3.3. Solution for Optimization Problem P1

In P1, all three dictionaries should be searched by the ACROA. ForDic1 and Dic2, the ACROA is straightly used, as only the 4 optimal parameters containing θ,ω,AandD should be obtained. For Dic3, assume the dimension of xn is N×1; 3×N optimal parameters must be obtained in every iteration procedure. That will generate a mass of calculation. Thus, the optimal problem of seeking in Dic3 is transformed into obtaining the optimal filter parameters as follows:(1)
Construct a bandpass filterχkλ, λ=ω,ωb,ωc (shown in Figure 2) is the parameter vector of the bandpass filter, and ω, ωb, and ωc determine the decay rate, the bandwidth, and the center frequency of the filter, respectively:(12)χkλ=sinωk−ωc+ωb+π2ω,ωc−ωb−π2ω≤k<ωc−ωb,1,ωc−ωb≤k≤ωc+ωb,cosωk−ωc−ωb,ωc+ωb<k≤ωc+ωb+π2ω,0,else.(2)
Solve the optimization problem P2 by the ACROA:(13)P2:MinimizeT3=D2ifftχkλr^ik22+λD2rin−ifftχkλr^ik22S.T.IMFi3n=ifftχkλr^ik.r^ik is the Fourier transformation of rin and ifft denotes the inverse Fourier transformation.(3)
IMFi3n is calculated by using the optimal parameter λ0:(14)IMFi3n=ifftχkλ0r^ik.Figure 2
The bandpass filter.To reduce the calculated quantity of BMD algorithm, the optimization problem P1 in time domain is translated to frequency domain. The representative solution of P1 is constructed by applying a bandpass filter. The decay rateω, the bandwidth ωb, and the center frequency ωc of the filter are optimized during the iteration process. And the optimal objective function T3 is used to constrain the obtained results to be effective components. Therefore, the final generated component is obtained by filtering the original signal using the optimized bandpass filter.


### 3.4. Simulation Analysis

Definesquareωn,rate as the square signals, the angular frequency is denoted as ω, and the duty ratio is denoted as rate. The composite signal x(n) of equation (15) contains a square component, an AM-FM sine component, and a noised signal n(t) with the signal noise ratio (SNR) of 3; the sampling rate is 12 kHz. The time domain curves of x(n) and its elements are illustrated in Figure 3:Figure 3
The time-domain waveforms ofx(n). (a) Mixed signal x(n). (b) Square component. (c) AM-FM sinusoidal component. (d) Noise component.
(a)(b)(c)(d)(15)xn=square20πn,0.5+2cos300πne−2t2+nt.For the sake of comparison,x(n) is separated using BMD, EEMD, and VMD, respectively. Figure 4 shows the decomposed results of x(n) by applying VMD. A series of apparent disturbances are shown in the figure, and the broadband feature of the square data is lacking, which is because the intrinsic theory of VMD is adaptive filtering; it will be influenced by Gibbs phenomenon. Figure 5 illustrates the separated results by applying EEMD. It can be seen that the obtained IMF1 corresponding to the square signal seems like a sine signal. Apparently, the obtained useful information is not correct. This is because EEMD uses an interpolating algorithm to gain the extrema envelopes, which will induce smooth IMFs. As the envelopes are calculated to obtain the IMFs, the eventual separated results of EEMD will turn into smooth IMFs when analyzing broadband datasets. Figure 6 illustrates the separated IMFs by using BMD. Therefore, the generated IMF1 is close to the real square component and the mutant feature is maintained.Figure 4
The decomposition results ofx(n) generated by VMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)Figure 5
The decomposition results ofx(n) generated by EEMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)Figure 6
The decomposition results ofx(n) generated by BMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)For the purpose of comparison, the accuracy parameters of the separated IMFs obtained by the 3 algorithms are illustrated in Table1, and the accuracy parameters containing the error of energy Ei, the coefficient of correlation ri, and the time T are applied [30]. Table 1 illustrates the IMFs obtained by BMD are closer to the square component in xn. To compare the computation time, the separating procedure is conducted on identical computer. The analysis results indicate the BMD algorithm requires more computational time compared with EEMD and VMD because of the complicated optimal procedure.Table 1
Evaluation parameters ofx(n) of the generated components.
Methodr1E1r2E2T(s)VMD0.96150.06860.97650.04644.7231EEMD0.80140.79830.90640.17838.9837BMD0.99880.00200.97820.013836.2351To analyze the performance of BMD in mode mixture phenomenon, a simulation signaly(n) shown in equation (16) is decomposed by VMD, EEMD, and BMD, respectively. To test the robust performance of BMD in different noise intensity, the SNRs of y(n) are set to be −10, −5, and 0, respectively. The sampling rate is 12 kHz. The time domain curves of y(n) and its elements are illustrated in Figure 7:(16)xn=square40πn,0.5+nt.Figure 7
The time-domain waveforms of simulation signaly(n) and its components. (a–c) The mixed signal y(n) as SNR = −10, −5, and 0, respectively. (d) The square component of y(n). (e–g) The noise components of the signals shown in (a), (b), and (c), respectively.
(a)(b)(c)(d)(e)(f)(g)Figures8 and 9 show the decomposed results of y(n) by applying VMD and EEMD. In Figure 8, it can be seen that, due to the low SNR, VMD treats the square signal as narrow band signal. Moreover, apparent disturbances are shown in Figure 8(a) when the SNR equals −10, which is caused by mode mixture. In Figure 9 of the EEMD generated decomposition results, the effective components are still close to sine signals due to the calculation of envelops. Meanwhile, when SNR equals −10 and −5, there are obvious disturbances mixed by the strong noise. Figure 10 shows the generated decomposition results of BMD. In Figures 10(a) and 10(b), BMD treats the square component as narrowband signal, which is caused by the low SNRs. Otherwise, the generated components are mixed by other modes, where it can be inferred that when a decomposed component of BMD belongs to the narrowband signal dictionary Dic3, there may still exist mode mixing problem when the noise is strong enough. However, the decomposed results of BMD more approach the real square component when the noise intensity varies, which indicates BMD is more robust to noise. For further analysis, the parameters of the separated IMFs obtained by the 3 algorithms are illustrated in Table 2; it can be easily concluded that although BMD costs more calculation time than the other 2 methods, it is still more robust to noise and behaves better in mode mixing.Figure 8
The decomposed results generated by VMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Figure 9
The decomposed results generated by EEMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Figure 10
The decomposed results generated by BMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Table 2
Evaluation parameters ofy(n) of the generated components.
MethodSNRr1E1T(s)VMD−100.77330.42384.1871−50.84320.35763.232400.88970.26265.7304EEMD−100.61741.49087.4127−50.70330.84176.926400.82030.57297.9327BMD−100.90410.129240.5328−50.93770.060435.986300.97210.011438.4571


## 3.1. Dictionary

To decompose both the broadband modes and the AM-FM narrowband modes from noisy signals, an associative dictionary containing common broadband data, such as sawtooth, square, and narrowband data, is constructed at first, in which the narrowband dictionary was defined in [30]. The analysis results are gained by seeking in the 3 dictionaries defined as follows:(8)Dic1=A1squareω1n+θ1,D1,Dic2=A2sawtoochω2n+θ2,D2,Dic3=A3ncosω3n+θ3n.θ, ω, and A represent the original phase, the angular frequency, and the amplitude, respectively. In Dic3, the maximal angular frequency of A3n is restricted to be much less than ω3. θ3n is restricted to be a slowly changing signal to guarantee the elements in Dic3 are narrowband [30].


## 3.2. BMD Iteration Procedure

The BMD algorithm is introduced to solve the existing disadvantages of extracting broadband features when using previous signal decomposition algorithms, including the Gibbs phenomenon in FT based algorithms and the interpolating of extrema in non-FT based algorithms. In BMD, the associative dictionaries are established at first. Then, the sparse decomposition results are gained by seeking in the dictionaries. BMD utilizes a normalized smoothness objective function as the optimal object. The parameters of the typical broadband datasets are used as the optimal parameters. The ACROA is used for the solution of the optimal problems [21]. In the end, both the broadband IMFs and the narrowband IMFs can be separated from noisy signals by using the optimal procedure. The flow chart of BMD is illustrated in Figure 1. The main iteration steps are designed as follows:(1)
r0n is set to equal xn.(2)
Separate the direct-current (DC) signal:(9)IMF0n=ifftr^01,r1n=xn−IMF0n.r^0k is the FT of r0n; ifftr^01 is the IFT of r^01.(3)
i is set to equal 1; the optimal problem P1 is constructed as follows:(10)P1:MinimizeT1A1,ω1,θ1,D1,T2A2,ω2,θ2,D2,T3A3n,ω3n,θ3nS.T.xn=∑i=0NIMFin+resn,IMFij∈Dicj,Tj=D2IMFijn22+λD2xn−IMFijn22.Tj is a normalized smooth operator to gain the sparse solution of r0n [30], D2 is the quadratic differential form operator, and λ>0 is always set to equal 1.(4)
Solve L1 by using the ACROA to gain the optimalTjj=1,2,3, find the minimum Tj, and choose the optimal IMFni of Tj. The procedure of the method is introduced in as follows.(5)
Updaterin: ri+1n=rin−IMFin.(6)
If the terminal condition in equation (11) is achieved, terminate the iteration procedure; else, get back to step 3. ε is a preset parameter which is normally set to be 1×10−6:(11)IMFin−IMFi−1n22IMFi−1n22≤ε.Figure 1
The iteration procedure of BMD.


## 3.3. Solution for Optimization Problem P1

In P1, all three dictionaries should be searched by the ACROA. ForDic1 and Dic2, the ACROA is straightly used, as only the 4 optimal parameters containing θ,ω,AandD should be obtained. For Dic3, assume the dimension of xn is N×1; 3×N optimal parameters must be obtained in every iteration procedure. That will generate a mass of calculation. Thus, the optimal problem of seeking in Dic3 is transformed into obtaining the optimal filter parameters as follows:(1)
Construct a bandpass filterχkλ, λ=ω,ωb,ωc (shown in Figure 2) is the parameter vector of the bandpass filter, and ω, ωb, and ωc determine the decay rate, the bandwidth, and the center frequency of the filter, respectively:(12)χkλ=sinωk−ωc+ωb+π2ω,ωc−ωb−π2ω≤k<ωc−ωb,1,ωc−ωb≤k≤ωc+ωb,cosωk−ωc−ωb,ωc+ωb<k≤ωc+ωb+π2ω,0,else.(2)
Solve the optimization problem P2 by the ACROA:(13)P2:MinimizeT3=D2ifftχkλr^ik22+λD2rin−ifftχkλr^ik22S.T.IMFi3n=ifftχkλr^ik.r^ik is the Fourier transformation of rin and ifft denotes the inverse Fourier transformation.(3)
IMFi3n is calculated by using the optimal parameter λ0:(14)IMFi3n=ifftχkλ0r^ik.Figure 2
The bandpass filter.To reduce the calculated quantity of BMD algorithm, the optimization problem P1 in time domain is translated to frequency domain. The representative solution of P1 is constructed by applying a bandpass filter. The decay rateω, the bandwidth ωb, and the center frequency ωc of the filter are optimized during the iteration process. And the optimal objective function T3 is used to constrain the obtained results to be effective components. Therefore, the final generated component is obtained by filtering the original signal using the optimized bandpass filter.


## 3.4. Simulation Analysis

Definesquareωn,rate as the square signals, the angular frequency is denoted as ω, and the duty ratio is denoted as rate. The composite signal x(n) of equation (15) contains a square component, an AM-FM sine component, and a noised signal n(t) with the signal noise ratio (SNR) of 3; the sampling rate is 12 kHz. The time domain curves of x(n) and its elements are illustrated in Figure 3:Figure 3
The time-domain waveforms ofx(n). (a) Mixed signal x(n). (b) Square component. (c) AM-FM sinusoidal component. (d) Noise component.
(a)(b)(c)(d)(15)xn=square20πn,0.5+2cos300πne−2t2+nt.For the sake of comparison,x(n) is separated using BMD, EEMD, and VMD, respectively. Figure 4 shows the decomposed results of x(n) by applying VMD. A series of apparent disturbances are shown in the figure, and the broadband feature of the square data is lacking, which is because the intrinsic theory of VMD is adaptive filtering; it will be influenced by Gibbs phenomenon. Figure 5 illustrates the separated results by applying EEMD. It can be seen that the obtained IMF1 corresponding to the square signal seems like a sine signal. Apparently, the obtained useful information is not correct. This is because EEMD uses an interpolating algorithm to gain the extrema envelopes, which will induce smooth IMFs. As the envelopes are calculated to obtain the IMFs, the eventual separated results of EEMD will turn into smooth IMFs when analyzing broadband datasets. Figure 6 illustrates the separated IMFs by using BMD. Therefore, the generated IMF1 is close to the real square component and the mutant feature is maintained.Figure 4
The decomposition results ofx(n) generated by VMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)Figure 5
The decomposition results ofx(n) generated by EEMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)Figure 6
The decomposition results ofx(n) generated by BMD. (a) The IMF corresponding to the square component. (b) The IMF corresponding to the AM-FM sinusoidal component.
(a)(b)(c)For the purpose of comparison, the accuracy parameters of the separated IMFs obtained by the 3 algorithms are illustrated in Table1, and the accuracy parameters containing the error of energy Ei, the coefficient of correlation ri, and the time T are applied [30]. Table 1 illustrates the IMFs obtained by BMD are closer to the square component in xn. To compare the computation time, the separating procedure is conducted on identical computer. The analysis results indicate the BMD algorithm requires more computational time compared with EEMD and VMD because of the complicated optimal procedure.Table 1
Evaluation parameters ofx(n) of the generated components.
Methodr1E1r2E2T(s)VMD0.96150.06860.97650.04644.7231EEMD0.80140.79830.90640.17838.9837BMD0.99880.00200.97820.013836.2351To analyze the performance of BMD in mode mixture phenomenon, a simulation signaly(n) shown in equation (16) is decomposed by VMD, EEMD, and BMD, respectively. To test the robust performance of BMD in different noise intensity, the SNRs of y(n) are set to be −10, −5, and 0, respectively. The sampling rate is 12 kHz. The time domain curves of y(n) and its elements are illustrated in Figure 7:(16)xn=square40πn,0.5+nt.Figure 7
The time-domain waveforms of simulation signaly(n) and its components. (a–c) The mixed signal y(n) as SNR = −10, −5, and 0, respectively. (d) The square component of y(n). (e–g) The noise components of the signals shown in (a), (b), and (c), respectively.
(a)(b)(c)(d)(e)(f)(g)Figures8 and 9 show the decomposed results of y(n) by applying VMD and EEMD. In Figure 8, it can be seen that, due to the low SNR, VMD treats the square signal as narrow band signal. Moreover, apparent disturbances are shown in Figure 8(a) when the SNR equals −10, which is caused by mode mixture. In Figure 9 of the EEMD generated decomposition results, the effective components are still close to sine signals due to the calculation of envelops. Meanwhile, when SNR equals −10 and −5, there are obvious disturbances mixed by the strong noise. Figure 10 shows the generated decomposition results of BMD. In Figures 10(a) and 10(b), BMD treats the square component as narrowband signal, which is caused by the low SNRs. Otherwise, the generated components are mixed by other modes, where it can be inferred that when a decomposed component of BMD belongs to the narrowband signal dictionary Dic3, there may still exist mode mixing problem when the noise is strong enough. However, the decomposed results of BMD more approach the real square component when the noise intensity varies, which indicates BMD is more robust to noise. For further analysis, the parameters of the separated IMFs obtained by the 3 algorithms are illustrated in Table 2; it can be easily concluded that although BMD costs more calculation time than the other 2 methods, it is still more robust to noise and behaves better in mode mixing.Figure 8
The decomposed results generated by VMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Figure 9
The decomposed results generated by EEMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Figure 10
The decomposed results generated by BMD. (a–c) The generated effective component ofy(n) as SNR = −10, −5, and 0, respectively. The red dotted line is the real component and the black solid line is the generated effective component.
(a)(b)(c)Table 2
Evaluation parameters ofy(n) of the generated components.
MethodSNRr1E1T(s)VMD−100.77330.42384.1871−50.84320.35763.232400.88970.26265.7304EEMD−100.61741.49087.4127−50.70330.84176.926400.82030.57297.9327BMD−100.90410.129240.5328−50.93770.060435.986300.97210.011438.4571


## 4. CMFE and LPP

### 4.1. CMFE

The DPMIG welding quality is mainly determined by the stability of the injected energy, which can be reflected by the nonlinearity degree of the input electrical signals. One of the most effective ways to measure the nonlinearity degree of nonstationary signals is entropy. When the entropy value is bigger, the nonlinearity degree of the collected DPMIG welding electrical signals is higher, which means the welding quality is worse. In the real physical world, the boundaries between different classes may be ambiguous, and it is difficult to determine whether an input pattern totally belongs to a given class or the other. Therefore, fuzzy entropy (FE) was proposed by applying the exponential function and removing the mean of matching template. The fuzzy entropy estimation procedure can be described as follows:(1)
A normalized time seriesxi:1≤i≤N with length N is considered for the sake of convenience. For given m, n, and r, a vector set sequence Xim,i=1,2,…,N−m+1 is constructed in the form as(17)Xim=xi,xi+1,…,xi+m−1−x0i,i=1,2,…,N−m+1,whereXim represents m consecutive u values, commencing with the ith point, and is obtained by removing their baseline:(18)x0i=m−1∑k=0m−1xi+k.(2)
For a certain vectorXim, the distance dijm between Xim and Xjm is defined as the maximum absolute difference of the corresponding scalar components, namely,(19)dijm=dXim,Xjm=maxk∈0,m−1ui+k−u0i−uj+k−u0j.(3)
Calculate the similarity degreeDijm of Xim and Xjm through the fuzzy function μdijm,n,r:(20)Dijm=μdijm,n,r=e−ln2dijm/rn.(4)
Define functionφmn,r as(21)φmn,r=1N−m∑i=1N−m1N−m−1∑j=1j≠iN−mDijm.(5)
Similarly, form + 1, repeat steps (1)–(4) and the function φm+1n,r can be obtained as(22)φm+1n,r=1N−m∑i=1N−m1N−m−1∑j=1j≠iN−mDijm+1.(6)
Finally, the FE of the original sequence is defined as the negative natural logarithm of the deviation ofφmn,r from φm+1n,r:(23)FuzzyEnm,n,r=limN⟶∞lnφmn,r−lnφm+1n,r.When the length of datasetsN is finite, FuzzyEnm,n,r can be estimated by the statistic as(24)FuzzyEnm,n,r,N=lnφmn,r−lnφm+1n,r.However, FE still lacked a clear statistical and physical significance for short-term time series or coarse-grained series. To alleviate the shortening of the time series in the coarse graining of MFE, CMFE applies the mean multiscale fuzzy entropy (MFE) of different coarse graining time series with the same scale factor as the extracted feature. If the CMFE of a given signal is larger, the mean MFE with a specific scale of different coarse graining time series of the signal is larger, which means that the given signal contains more mode information and is more complicated. To verify the stability of welding quantitatively and intuitively, the CMFE of the collected inverter power source electrical signals is calculated. A larger CMFE means a larger uncertainty in the inverter power source output energy, which will lead to a worse welding quantitatively. The calculation process of CMFE is as follows:(1)
Suppose a signal functionxn, set the number of data points as N, and calculate ykτ=yk,1τ,yk,2τ,…,yk,jτ of different coarse-grained time series with scale factor as τ. ykτ is the kth coarse-grained time series of xn with the scale factor of τ:(25)yk,jτ=1τ∑i=j−1τ+kjτ+k−1xi,1≤j≤Nτ,1≤k≤τ.(2)
Secondly, the results of fuzzy entropy of each coarse-grained time seriesykτ obtained by calculation were averaged to obtain the coarse-grained time series CMFE of the scale factor τ.(26)CMFExn,τ=1τ∑k=1τFuzzyEnykτ.


### 4.2. LPP Algorithm

After the evaluating parameters are calculated from the electric datasets, the feature with the dimension ofM is obtained. Nevertheless, the obtained features usually contain redundant information and are high dimensional. Therefore, suitable algorithm should be used for the useful information extraction and the feature dimension reduction. PCA is a widely applied algorithm for the reduction of dimension [28, 29], which means increasing the global feature and lack of the ability for maintaining the local structure. Thus, LPP algorithm which can extract the local feature of the given datasets is used for further feature extraction in this article.LPP maintains the local feature through the extraction of the approximation of the given dataset.F=f1,f2,…,fNi is set to be the feature vectors. Ni and J are the amount of the feature vectors and the data samples. LPP seeks a projection matrix to reduce the dimensions of the feature vectors into F=f1′,f2′,…,fNi′; J′J′≪J is the dimension of the fitted feature vector F′. For choosing a map like that, the object is constructed as follows [28, 29]:(27)min∑i,j=1Nifi−fj2Si,j.The parameterS is obtained by the nearest neighborhood and induces a punishment when the neighbor area fi is far away from fj. The formula of S is as follows:(28)Si,j=exp−fi−fj22σ2,fiorfjis between theknearest neighbors offjorfi,respectively,0,k represents the “locality” degree of the local area and σ shows the major kernel value. The objective functions are optimized to make sure fj is near fi; thus, fj′ is also near fi′. Therefore, the local feature is maintained. Set fi′=aTfi; then, the optimal solution is altered as follows:(29)Jla=12∑i,j=1Nifi′−fj′2Si,j=12∑i,j=1NiaTfi−aTfj2Si,j=aTFD−SFTa=aTFLFTa.A Laplacian matrixL=D−S is constructed. The bigger Dii obtains a more important fi′. Furthermore, the restriction constraint is set up as follows:(30)F′TDF′=1⟹aFDFTa=1.The optimal problem is constructed as follows:(31)minaaFLFTasubject toaFDFTa=1.The constrained optimal object function can be transformed into an eigenvalue question by using Lagrange multiplier algorithm. The transform matrix is gained as follows:(32)FLFTa=λFDFTa,where λ represents the weight value. The dimension reduction formula is as follows:(33)F′=ATFA=a1,a2,…,alT,where A is a matrix with the dimension of J×J′. The input vector F is transformed before a dimension reduced vector F′ is gained. The local useful information is maintained by F′ and the feature is extracted. The optimal linear local structure is searched in the matrix. Therefore, LPP is adaptive for noisy signals and isolated feature values as the local structure is maintained.PCA is applied to extract the subspace of maximal difference. Therefore, it can maintain the global information. However, the useful features of the local structure of experimental signals like electric datasets are more meaningful. In addition, LPP has been identified to be adaptive for dealing with the electric datasets. At the same time, researchers have made comparisons between PCA and LPP by using various industrial datasets; the analysis results have demonstrated that PCA performed better in the useful information extraction of features. Therefore, LPP algorithm, which can maintain the useful information, is used in the introduced BMD algorithm for extracting effective feature better. In this paper, the CMFEs of the first two decomposed components are applied as the calculated parameter of the analyzed signals; the number of the generated parameters is determined by the number of the scale factorτ. In the experimental analysis part, τ is set to 1, 2, 3, and 4 as we find the parameters become more chaos when τ is bigger than 4. Therefore, there are 8 generated CMFE parameters with different τ. Then, the CMFE parameters are fused into one parameter vector by LPP to merge the effective information.


## 4.1. CMFE

The DPMIG welding quality is mainly determined by the stability of the injected energy, which can be reflected by the nonlinearity degree of the input electrical signals. One of the most effective ways to measure the nonlinearity degree of nonstationary signals is entropy. When the entropy value is bigger, the nonlinearity degree of the collected DPMIG welding electrical signals is higher, which means the welding quality is worse. In the real physical world, the boundaries between different classes may be ambiguous, and it is difficult to determine whether an input pattern totally belongs to a given class or the other. Therefore, fuzzy entropy (FE) was proposed by applying the exponential function and removing the mean of matching template. The fuzzy entropy estimation procedure can be described as follows:(1)
A normalized time seriesxi:1≤i≤N with length N is considered for the sake of convenience. For given m, n, and r, a vector set sequence Xim,i=1,2,…,N−m+1 is constructed in the form as(17)Xim=xi,xi+1,…,xi+m−1−x0i,i=1,2,…,N−m+1,whereXim represents m consecutive u values, commencing with the ith point, and is obtained by removing their baseline:(18)x0i=m−1∑k=0m−1xi+k.(2)
For a certain vectorXim, the distance dijm between Xim and Xjm is defined as the maximum absolute difference of the corresponding scalar components, namely,(19)dijm=dXim,Xjm=maxk∈0,m−1ui+k−u0i−uj+k−u0j.(3)
Calculate the similarity degreeDijm of Xim and Xjm through the fuzzy function μdijm,n,r:(20)Dijm=μdijm,n,r=e−ln2dijm/rn.(4)
Define functionφmn,r as(21)φmn,r=1N−m∑i=1N−m1N−m−1∑j=1j≠iN−mDijm.(5)
Similarly, form + 1, repeat steps (1)–(4) and the function φm+1n,r can be obtained as(22)φm+1n,r=1N−m∑i=1N−m1N−m−1∑j=1j≠iN−mDijm+1.(6)
Finally, the FE of the original sequence is defined as the negative natural logarithm of the deviation ofφmn,r from φm+1n,r:(23)FuzzyEnm,n,r=limN⟶∞lnφmn,r−lnφm+1n,r.When the length of datasetsN is finite, FuzzyEnm,n,r can be estimated by the statistic as(24)FuzzyEnm,n,r,N=lnφmn,r−lnφm+1n,r.However, FE still lacked a clear statistical and physical significance for short-term time series or coarse-grained series. To alleviate the shortening of the time series in the coarse graining of MFE, CMFE applies the mean multiscale fuzzy entropy (MFE) of different coarse graining time series with the same scale factor as the extracted feature. If the CMFE of a given signal is larger, the mean MFE with a specific scale of different coarse graining time series of the signal is larger, which means that the given signal contains more mode information and is more complicated. To verify the stability of welding quantitatively and intuitively, the CMFE of the collected inverter power source electrical signals is calculated. A larger CMFE means a larger uncertainty in the inverter power source output energy, which will lead to a worse welding quantitatively. The calculation process of CMFE is as follows:(1)
Suppose a signal functionxn, set the number of data points as N, and calculate ykτ=yk,1τ,yk,2τ,…,yk,jτ of different coarse-grained time series with scale factor as τ. ykτ is the kth coarse-grained time series of xn with the scale factor of τ:(25)yk,jτ=1τ∑i=j−1τ+kjτ+k−1xi,1≤j≤Nτ,1≤k≤τ.(2)
Secondly, the results of fuzzy entropy of each coarse-grained time seriesykτ obtained by calculation were averaged to obtain the coarse-grained time series CMFE of the scale factor τ.(26)CMFExn,τ=1τ∑k=1τFuzzyEnykτ.


## 4.2. LPP Algorithm

After the evaluating parameters are calculated from the electric datasets, the feature with the dimension ofM is obtained. Nevertheless, the obtained features usually contain redundant information and are high dimensional. Therefore, suitable algorithm should be used for the useful information extraction and the feature dimension reduction. PCA is a widely applied algorithm for the reduction of dimension [28, 29], which means increasing the global feature and lack of the ability for maintaining the local structure. Thus, LPP algorithm which can extract the local feature of the given datasets is used for further feature extraction in this article.LPP maintains the local feature through the extraction of the approximation of the given dataset.F=f1,f2,…,fNi is set to be the feature vectors. Ni and J are the amount of the feature vectors and the data samples. LPP seeks a projection matrix to reduce the dimensions of the feature vectors into F=f1′,f2′,…,fNi′; J′J′≪J is the dimension of the fitted feature vector F′. For choosing a map like that, the object is constructed as follows [28, 29]:(27)min∑i,j=1Nifi−fj2Si,j.The parameterS is obtained by the nearest neighborhood and induces a punishment when the neighbor area fi is far away from fj. The formula of S is as follows:(28)Si,j=exp−fi−fj22σ2,fiorfjis between theknearest neighbors offjorfi,respectively,0,k represents the “locality” degree of the local area and σ shows the major kernel value. The objective functions are optimized to make sure fj is near fi; thus, fj′ is also near fi′. Therefore, the local feature is maintained. Set fi′=aTfi; then, the optimal solution is altered as follows:(29)Jla=12∑i,j=1Nifi′−fj′2Si,j=12∑i,j=1NiaTfi−aTfj2Si,j=aTFD−SFTa=aTFLFTa.A Laplacian matrixL=D−S is constructed. The bigger Dii obtains a more important fi′. Furthermore, the restriction constraint is set up as follows:(30)F′TDF′=1⟹aFDFTa=1.The optimal problem is constructed as follows:(31)minaaFLFTasubject toaFDFTa=1.The constrained optimal object function can be transformed into an eigenvalue question by using Lagrange multiplier algorithm. The transform matrix is gained as follows:(32)FLFTa=λFDFTa,where λ represents the weight value. The dimension reduction formula is as follows:(33)F′=ATFA=a1,a2,…,alT,where A is a matrix with the dimension of J×J′. The input vector F is transformed before a dimension reduced vector F′ is gained. The local useful information is maintained by F′ and the feature is extracted. The optimal linear local structure is searched in the matrix. Therefore, LPP is adaptive for noisy signals and isolated feature values as the local structure is maintained.PCA is applied to extract the subspace of maximal difference. Therefore, it can maintain the global information. However, the useful features of the local structure of experimental signals like electric datasets are more meaningful. In addition, LPP has been identified to be adaptive for dealing with the electric datasets. At the same time, researchers have made comparisons between PCA and LPP by using various industrial datasets; the analysis results have demonstrated that PCA performed better in the useful information extraction of features. Therefore, LPP algorithm, which can maintain the useful information, is used in the introduced BMD algorithm for extracting effective feature better. In this paper, the CMFEs of the first two decomposed components are applied as the calculated parameter of the analyzed signals; the number of the generated parameters is determined by the number of the scale factorτ. In the experimental analysis part, τ is set to 1, 2, 3, and 4 as we find the parameters become more chaos when τ is bigger than 4. Therefore, there are 8 generated CMFE parameters with different τ. Then, the CMFE parameters are fused into one parameter vector by LPP to merge the effective information.


## 5. Experimental Analysis

### 5.1. Data Collection

The experimental electric datasets are collected from the built aluminum alloy DPMIG welding experimental platform illustrated in Figure11. The experimental platform applies a DPMIG welding machine and a weld robot to proceed with the welding procedure. The aluminum welding machine model is DPMIG-5, which is responsible for setting welding process parameters such as welding peak voltage, welding base value voltage, welding peak current, welding base value current, duty cycle, welding frequency, and the starting and stopping welding experiments. The robot includes the welding robot ABB IRB1410 and the robot controller to control welding speed and walking direction. The data acquisition system consists of a sensor box, an acquisition box, a computer, and a communication cable, which is used to store welding arc data with high speed and accuracy. The wire feeding mechanism includes a wire feeding disc and a wire feeding controller, which can provide two kinds of wire feeding modes for selection: manual wire feeding and automatic wire feeding. In the welding process experiment of argon gas storage tank, inert gas is provided during the transportation of welding wire to isolate the air outside the welding area and prevent oxidation.Figure 11
The aluminum alloy DPMIG welding test bed. ①–⑧ are the welding robot, the wire feeder, the acquisition box, the computer, the welding machine, the electrical signal sensor box, the argon gas storage tank, and the inverter power source.The scantling of the aluminum alloy is 0.25 m × 0.2 m × 0.005 m. The argon flow rate is 0.2 L/s, and the semidiameter of the welding stick is 0.6 mm. The datasets are obtained by using electrical transducers. To test the validity of BMD in various operating conditions, 4 experiments with various current frequency, duty cycle, and welding speed are applied. A 48-second data collection procedure is conducted in each experiment, and the sample frequency is set to be 12 kHz. Table3 shows the experimental parameters. Figure 12 shows the captured electrical datasets of the 4 experiments. The welding appearances are shown in Figure 13. The welding appearance of dataset 3 illustrated in Figure 13(c) has apparent incomplete-penetration and fusion caused by the arc breaking of the input broadband current data.Table 3
Parameters of the welding experiment.
DatasetWelding speed (mm/s)Inputting currentPeak-to-peak value (A)Duty cycle (%)Pulse frequency (Hz)115050221.55050235505024560505Figure 12
The collected inverter power source output current signals. (a–d) The time-domain waveform of datasets 1–4, respectively.
(a)(b)(c)(d)Figure 13
The appearance of the weld shapes. (a–d) Photos of the weld shapes corresponding to datasets 1–4, respectively.
(a)(b)(c)(d)


### 5.2. Experimental Signal Decomposition

One-second current data in dataset 3 is chosen for the comparison of VMD, EEMD, and BMD, and the time domain waveform of the signal is illustrated in Figure14. The original signal is decomposed into 4 IMFs and a residual component. The feature corresponding to the quality of welding is extracted from the input broadband and noisy disturbance. Thus, the useful components, the input broadband square component, and the noisy components should be separated. The generated results of VMD are illustrated in Figure 15. The results indicate that VMD is unable to differentiate the useful components and broadband square components, which implies it lacks the ability of separating broadband and narrowband components. The generated results of VMD are illustrated in Figure 16. The IMFs generated by EEMD are quite distinct from the input square signal as it handles the broadband square signals to be sinusoidal signals that will generate serious mode confusion. The generated results of BMD are illustrated in Figure 17. The IMF1 of BMD is similar to the input square signal with the parameters of Table 3, and the IMF2 maintains the mutant feature of the quality of DPMIG welding.Figure 14
The selected experimental signal of dataset 3.Figure 15
Decomposition results generated by VMD.
(a)(b)(c)(d)(e)Figure 16
Decomposition results generated by EEMD.
(a)(b)(c)(d)(e)Figure 17
Decomposition results generated by BMD.
(a)(b)(c)(d)(e)


### 5.3. LPP Feature Extraction

Each of the four datasets is divided into 48 equal samples; each sample contains one-second signal. All the samples are decomposed by VMD, BMD, and EEMD, respectively. Then, the CMFEs of the generated effective components are calculated to verify the stability of welding quantitatively and intuitively (the scale factorτ is set to 1, 2, 3, 4). Then, the features are fused by LPP, the fused feature vectors are shown in Figures 18–20, and the mean CMFEs in each experiment are shown in Table 4. It can be seen in Figure 20 and Table 4 that the CMFEs generated by applying BMD become larger when the welding speed increases, which means that the welding quality is worse. This is because the electrical energy becomes more dispersed with the increase in the welding speed. In addition, the larger pulse frequency matches better with the welding speed, and the arc energy becomes more stable. Thus, the CMFEs of dataset 3 are larger than that of dataset 4. However, the CMFEs obtained by BMD and EEMD are confusing, which further indicates the extracted feature is disturbed by other modes from the broadband signal and noise. To compare the computation time, the mean time cost parameters of the analysis procedure of the 4 datasets by applying the three adaptive decomposition methods combined with CMFE and LPP are given out in Table 4. The analysis results show that it still costs more time when applying BMD method to the experimental signals.Figure 18
The LPP-fused CMFEs of the VMD-generated results.Figure 19
The LPP-fused CMFEs of the EEMD-generated results.Figure 20
The LPP-fused CMFEs of the BMD-generated results.Table 4
The mean fused CMFEs and the mean time cost of the generated results.
MethodCMFET(s)Dataset 1Dataset 2Dataset 3Dataset 4The mean time cost of the 4 datasetsVMD0.11070.13950.18210.1536213.1275EEMD0.11450.13410.19370.1547410.2391BMD0.09970.11910.16540.13221315.5237The simulated and experimental analyzing results indicate BMD exceeds the other two methods in both the precision and denoising performance. Furthermore, in the experimental analysis, BMD performed better when combined with LPP in the quality feature extraction of aluminum alloy DPMIG welding.


## 5.1. Data Collection

The experimental electric datasets are collected from the built aluminum alloy DPMIG welding experimental platform illustrated in Figure11. The experimental platform applies a DPMIG welding machine and a weld robot to proceed with the welding procedure. The aluminum welding machine model is DPMIG-5, which is responsible for setting welding process parameters such as welding peak voltage, welding base value voltage, welding peak current, welding base value current, duty cycle, welding frequency, and the starting and stopping welding experiments. The robot includes the welding robot ABB IRB1410 and the robot controller to control welding speed and walking direction. The data acquisition system consists of a sensor box, an acquisition box, a computer, and a communication cable, which is used to store welding arc data with high speed and accuracy. The wire feeding mechanism includes a wire feeding disc and a wire feeding controller, which can provide two kinds of wire feeding modes for selection: manual wire feeding and automatic wire feeding. In the welding process experiment of argon gas storage tank, inert gas is provided during the transportation of welding wire to isolate the air outside the welding area and prevent oxidation.Figure 11
The aluminum alloy DPMIG welding test bed. ①–⑧ are the welding robot, the wire feeder, the acquisition box, the computer, the welding machine, the electrical signal sensor box, the argon gas storage tank, and the inverter power source.The scantling of the aluminum alloy is 0.25 m × 0.2 m × 0.005 m. The argon flow rate is 0.2 L/s, and the semidiameter of the welding stick is 0.6 mm. The datasets are obtained by using electrical transducers. To test the validity of BMD in various operating conditions, 4 experiments with various current frequency, duty cycle, and welding speed are applied. A 48-second data collection procedure is conducted in each experiment, and the sample frequency is set to be 12 kHz. Table3 shows the experimental parameters. Figure 12 shows the captured electrical datasets of the 4 experiments. The welding appearances are shown in Figure 13. The welding appearance of dataset 3 illustrated in Figure 13(c) has apparent incomplete-penetration and fusion caused by the arc breaking of the input broadband current data.Table 3
Parameters of the welding experiment.
DatasetWelding speed (mm/s)Inputting currentPeak-to-peak value (A)Duty cycle (%)Pulse frequency (Hz)115050221.55050235505024560505Figure 12
The collected inverter power source output current signals. (a–d) The time-domain waveform of datasets 1–4, respectively.
(a)(b)(c)(d)Figure 13
The appearance of the weld shapes. (a–d) Photos of the weld shapes corresponding to datasets 1–4, respectively.
(a)(b)(c)(d)


## 5.2. Experimental Signal Decomposition

One-second current data in dataset 3 is chosen for the comparison of VMD, EEMD, and BMD, and the time domain waveform of the signal is illustrated in Figure14. The original signal is decomposed into 4 IMFs and a residual component. The feature corresponding to the quality of welding is extracted from the input broadband and noisy disturbance. Thus, the useful components, the input broadband square component, and the noisy components should be separated. The generated results of VMD are illustrated in Figure 15. The results indicate that VMD is unable to differentiate the useful components and broadband square components, which implies it lacks the ability of separating broadband and narrowband components. The generated results of VMD are illustrated in Figure 16. The IMFs generated by EEMD are quite distinct from the input square signal as it handles the broadband square signals to be sinusoidal signals that will generate serious mode confusion. The generated results of BMD are illustrated in Figure 17. The IMF1 of BMD is similar to the input square signal with the parameters of Table 3, and the IMF2 maintains the mutant feature of the quality of DPMIG welding.Figure 14
The selected experimental signal of dataset 3.Figure 15
Decomposition results generated by VMD.
(a)(b)(c)(d)(e)Figure 16
Decomposition results generated by EEMD.
(a)(b)(c)(d)(e)Figure 17
Decomposition results generated by BMD.
(a)(b)(c)(d)(e)


## 5.3. LPP Feature Extraction

Each of the four datasets is divided into 48 equal samples; each sample contains one-second signal. All the samples are decomposed by VMD, BMD, and EEMD, respectively. Then, the CMFEs of the generated effective components are calculated to verify the stability of welding quantitatively and intuitively (the scale factorτ is set to 1, 2, 3, 4). Then, the features are fused by LPP, the fused feature vectors are shown in Figures 18–20, and the mean CMFEs in each experiment are shown in Table 4. It can be seen in Figure 20 and Table 4 that the CMFEs generated by applying BMD become larger when the welding speed increases, which means that the welding quality is worse. This is because the electrical energy becomes more dispersed with the increase in the welding speed. In addition, the larger pulse frequency matches better with the welding speed, and the arc energy becomes more stable. Thus, the CMFEs of dataset 3 are larger than that of dataset 4. However, the CMFEs obtained by BMD and EEMD are confusing, which further indicates the extracted feature is disturbed by other modes from the broadband signal and noise. To compare the computation time, the mean time cost parameters of the analysis procedure of the 4 datasets by applying the three adaptive decomposition methods combined with CMFE and LPP are given out in Table 4. The analysis results show that it still costs more time when applying BMD method to the experimental signals.Figure 18
The LPP-fused CMFEs of the VMD-generated results.Figure 19
The LPP-fused CMFEs of the EEMD-generated results.Figure 20
The LPP-fused CMFEs of the BMD-generated results.Table 4
The mean fused CMFEs and the mean time cost of the generated results.
MethodCMFET(s)Dataset 1Dataset 2Dataset 3Dataset 4The mean time cost of the 4 datasetsVMD0.11070.13950.18210.1536213.1275EEMD0.11450.13410.19370.1547410.2391BMD0.09970.11910.16540.13221315.5237The simulated and experimental analyzing results indicate BMD exceeds the other two methods in both the precision and denoising performance. Furthermore, in the experimental analysis, BMD performed better when combined with LPP in the quality feature extraction of aluminum alloy DPMIG welding.


## 6. Conclusions

Combined with LPP, the BMD algorithm is introduced for the quality feature extraction of nonstationary broadband electrical signals of aluminum alloy double pulse metal inert gas (DPMIG) welding, and monitoring the influence of the electrical signals on the arc stability and welding appearance. The disadvantages of previous time-frequency algorithms in FT based and non-FT based algorithms induced by the calculation of extrema and Gibbs phenomenon can be alleviated by seeking in the dictionary using ACROA. The superiority of BMD is tested by simulated and experimental data. Analysis results indicate that the proposed method can precisely draw the broadband useful features and evaluate the quality of aluminum alloy DPMIG welding. The proposed method is tested by aluminum alloy DPMIG welding electrical signals. BMD algorithm can also be applied to analyze other broadband industrial data. Nevertheless, the calculation of BMD is still relatively large because of the optimization procedure. Thus, fast iteration algorithms of BMD should be studied in the future.



---

*Source: 7576034-2020-07-18.xml*