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

. A novel adaptive signal decomposition algorithm, broadband mode decomposition (BMD), is proposed for analyzing non-stationary 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 eﬀective components are calculated. Then, locality preserving projection (LPP) is applied for further feature extraction. Analysis results show BMD is more eﬀective when drawing broadband feature from noise and BMD is adaptive for the quality monitoring of DPMIG welding.


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. us, 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]. erefore, 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][4][5].
ese adaptive time-frequency analysis algorithms can be classified into two categories: FT based algorithms and non-FT based algorithms [6][7][8].
e 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][12][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][15][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][18][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. e 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. us, 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. en, the sparse solution is gained by seeking in the associative dictionary applying an optimization method. e BMD method utilizes a regulated smooth operator as the optimal object. e initial phase, the frequency, and the amplitude of common broadband data are used as the optimal arguments. e 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.
us, 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. us, 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][23][24][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 lowdimensional space. LPP has been demonstrated to be suitable for complex data processing in a plenty of researches [27][28][29]. erefore, it is selected for fusing the feature vectors generated from the obtained MIG welding electric datasets. e remaining contents are organized as below. In Section 2, the existing disadvantages in analyzing the broadband data are given out. e 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. e conclusion is given out at last.

Existing Problems in Broadband
Signal Analysis

VMD Method.
e Fourier progression unfolds the data into a series of sinusoidal components. e square signal is taken as an illustration of broadband data; the Fourier progression is as follows: 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.
is is named as Gibbs phenomenon [20].
VMD decomposes a dataset into a series of band limited IMFs by using multiple Wiener filters. e 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. e iteration procedure for obtaining IMFs is as follows [10]: (1) Calculate the mathematical formula of IMF i (n) by using the Hilbert transform: where δ is the Dirac distribution and denotes convolution.

Mathematical Problems in Engineering
(2) Compute the center frequency ω i of IMF i (n); transform the formula in equation (2) to the frequency domain corresponding to ω i as follows: (3) An optimization problem is constructed as equation (4); the optimal object is the L 2 norm of the gradient value of the frequency domain formula. e separated results are gained by solving the optimization problem using the alternate direction multiplier method:

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. e iteration procedure of EMD is as follows [8]: (1) Set h i (n) � x(n), i � 1, and j � 1.
(2) Interpolate the extrema of h i (n) using the interpolating method to gain the lower envelope and upper envelope m u (n) and m l (n); then calculate h i (n): (3) Step (2) is repeated until h i+1 (n) accords with the condition of the IMFs [8]. e jth IMF is obtained: (4) Set r(n) � x(n) − IMF j (n); set the residue r(n) to be the updated x(n). Steps (1)-(3) are repeated until r(n) is monotonic data. e separated results are expressed as follows: 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. us, 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]. e 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. at 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.

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]. e analysis results are gained by seeking in the 3 dictionaries defined as follows: θ, ω, and A represent the original phase, the angular frequency, and the amplitude, respectively. In Di c 3 , the maximal angular frequency of A 3 (n) is restricted to be much less than ω 3 . θ 3 (n) is restricted to be a slowly changing signal to guarantee the elements in Dic 3 are narrowband [30].

BMD Iteration Procedure.
e 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. en, the sparse decomposition results are gained by seeking in the dictionaries. BMD utilizes a normalized smoothness objective function as the optimal object. e parameters of the typical broadband datasets are used as the optimal parameters. e 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. e flow chart of BMD is illustrated in Figure 1. e main iteration steps are designed as follows: (1) r 0 (n) is set to equal x(n).
(2) Separate the direct-current (DC) signal: r 0 (k) is the FT of r 0 (n); ifft[r 0 (1)] is the IFT of r 0 (1). (3) i is set to equal 1; the optimal problem P1 is constructed as follows: Mathematical Problems in Engineering 3 T j is a normalized smooth operator to gain the sparse solution of r 0 (n) [30], D (2) is the quadratic differential form operator, and λ > 0 is always set to equal 1.

No
Extract the DC component

Continue or end the iteration process
Obtain the minimum value T j between the optimized T 1 , (4) Solve L1 by using the ACROA to gain the optimal T j(j�1, 2,3) , find the minimum T j , and choose the optimal IMF i (n) of T j . e procedure of the method is introduced in as follows.
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 :

Solution for Optimization Problem P1.
In P1, all three dictionaries should be searched by the ACROA. For Dic 1 and Dic 2 , the ACROA is straightly used, as only the 4 optimal parameters containing θ, ω, A and D should be obtained. For Dic 3 , assume the dimension of x(n) is N × 1; 3 × N optimal parameters must be obtained in every iteration procedure. at will generate a mass of calculation. us, the optimal problem of seeking in Dic 3 is transformed into obtaining the optimal filter parameters as follows: (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: (2) Solve the optimization problem P2 by the ACROA: r i (k) is the Fourier transformation of r i (n) and ifft denotes the inverse Fourier transformation.
is calculated by using the optimal parameter λ 0 : To reduce the calculated quantity of BMD algorithm, the optimization problem P1 in time domain is translated to frequency domain. e representative solution of P1 is constructed by applying a bandpass filter. e decay rate ω, the bandwidth ω b , and the center frequency ω c of the filter are optimized during the iteration process. And the optimal objective function T 3 is used to constrain the obtained results to be effective components. erefore, the final generated component is obtained by filtering the original signal using the optimized bandpass filter.

Simulation Analysis.
Define square(ωn, rate) as the square signals, the angular frequency is denoted as ω, and the duty ratio is denoted as rate. e 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. e time domain curves of x(n) and its elements are illustrated in Figure 3: 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 IMF 1 corresponding to the square signal seems like a sine signal. Apparently, the obtained useful information is not correct. is 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. erefore, the generated IMF 1 is close to the real square component and the mutant feature is maintained.
For the purpose of comparison, the accuracy parameters of the separated IMFs obtained by the 3 algorithms are illustrated in Table 1, and the accuracy parameters containing the error of energy E i , the coefficient of correlation r i , and the time T are applied [30]. Table 1 illustrates the IMFs obtained by BMD are closer to the square component in x(n). To compare the computation time, the separating procedure is conducted on identical computer. e analysis results indicate the BMD algorithm requires more Mathematical Problems in Engineering computational time compared with EEMD and VMD because of the complicated optimal procedure. To analyze the performance of BMD in mode mixture phenomenon, a simulation signal y(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. e sampling rate is 12 kHz. e time domain curves of y(n) and its elements are illustrated in Figure 7: x(n) � square(40πn, 0.5) + n(t).
(16) Figures 8 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 Dic 3 , 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.

CMFE.
e 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. erefore, fuzzy entropy (FE) was proposed by applying the exponential function and removing the mean of matching template. e fuzzy entropy estimation procedure can be described as follows: (1) A normalized time series x(i): 1 ≤ i ≤ N { } with length N is considered for the sake of convenience. For given m, n, and r, a vector set sequence X m i , i � 1, 2, . . . , N − m + 1 is constructed in the form as where X m i represents m consecutive u values, commencing with the ith point, and is obtained by removing their baseline: (2) For a certain vector X m i , the distance d m ij between X m i and X m j is defined as the maximum absolute difference of the corresponding scalar components, namely, Similarly, for m + 1, repeat steps (1)-(4) and the function φ m+1 (n, r) can be obtained as When the length of datasets N is finite, FuzzyEn(m, n, r) can be estimated by the statistic as FuzzyEn(m, n, r, N) � ln φ m (n, r) − ln φ m+1 (n, r). (24) 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. e calculation process of CMFE is as follows: (1) Suppose a signal function x(n), set the number of data points as N, and calculate y (τ) k � y (τ) k,1 , y (τ) k,2 , . . . , y (τ) k,j of different coarse-grained time series with scale factor as τ. y (τ) k is the kth coarsegrained time series of x(n) with the scale factor of τ:   (2) Secondly, the results of fuzzy entropy of each coarsegrained time series y (τ) k obtained by calculation were averaged to obtain the coarse-grained time series CMFE of the scale factor τ.

LPP Algorithm.
After the evaluating parameters are calculated from the electric datasets, the feature with the dimension of M is obtained. Nevertheless, the obtained features usually contain redundant information and are high dimensional. erefore, 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. us, 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 � f 1 , f 2 , . . . , f N i } is set to be the feature vectors. N i 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 is the dimension of the fitted feature vector F ′ . For choosing a map like that, the object is constructed as follows [28,29]: e parameter S is obtained by the nearest neighborhood and induces a punishment when the neighbor area f i is far away from f j . e formula of S is as follows: k represents the "locality" degree of the local area and σ shows the major kernel value. e objective functions are optimized to make sure f j is near f i ; thus, f j ′ is also near f i ′ . erefore, the local feature is maintained. Set f i ′ � a T f i ; then, the optimal solution is altered as follows: A Laplacian matrix L � D − S is constructed. e bigger D ii obtains a more important f i ′ . Furthermore, the restriction constraint is set up as follows: e optimal problem is constructed as follows: min a aFLF T a subject to aFDF T a � 1.
(31) e constrained optimal object function can be transformed into an eigenvalue question by using Lagrange multiplier algorithm. e transform matrix is gained as follows: where λ represents the weight value. e dimension reduction formula is as follows: where A is a matrix with the dimension of J × J ′ . e input vector F is transformed before a dimension reduced vector F ′ is gained. e local useful information is maintained by F ′ and the feature is extracted. e optimal linear local structure is searched in the matrix. erefore, 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. erefore, 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. erefore, 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. erefore, there are 8 generated CMFE parameters with different τ. en, the CMFE parameters are fused into one parameter vector by LPP to merge the effective information.

Data Collection.
e experimental electric datasets are collected from the built aluminum alloy DPMIG welding experimental platform illustrated in Figure 11. e experimental platform applies a DPMIG welding machine and a weld robot to proceed with the welding procedure. e 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. e robot includes the welding robot ABB IRB1410 and the robot controller to control welding speed and walking direction. e 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. e 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.
e scantling of the aluminum alloy is 0.25 m × 0.2 m × 0.005 m. e argon flow rate is 0.2 L/s, and the semidiameter of the welding stick is 0.6 mm. e 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. Table 3 shows the experimental parameters. Figure 12 shows the captured electrical datasets of the 4 experiments. e welding appearances are shown in Figure 13. e 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.

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 Figure 14.
e original signal is decomposed into 4 IMFs and a residual component. e feature corresponding to the quality of welding is extracted from the input broadband and noisy disturbance. us, the useful components, the input broadband square component, and the noisy components should be separated. e generated results of VMD are illustrated in Figure 15. e 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. e generated results of VMD are illustrated in Figure 16. e 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. e generated results of BMD are illustrated in Figure 17. e IMF 1 of BMD is similar to the input square signal with the parameters of Table 3, and the IMF 2 maintains the mutant feature of the quality of DPMIG welding.

LPP Feature Extraction.
Each of the four datasets is divided into 48 equal samples; each sample contains onesecond signal. All the samples are decomposed by VMD, BMD, and EEMD, respectively. en, 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). en, 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. is 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. us, 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. e analysis results show that it still costs more time when applying BMD method to the experimental signals. e 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.  Figure 11: e 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.

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. e 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.
e 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. e 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. us, fast iteration algorithms of BMD should be studied in the future.

Data Availability
All data included in this study are available upon request to the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.