Next Article in Journal
Advanced Machine Learning Functionalities in the Medium Voltage Distributed Monitoring System QuEEN: A Macro-Regional Voltage Dips Severity Analysis
Previous Article in Journal
Harmonic Profile Enhancement of Grid Connected Fuel Cell through Cascaded H-Bridge Multi-Level Inverter and Improved Squirrel Search Optimization Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Photovoltaic Power Quality Analysis Based on the Modulation Broadband Mode Decomposition Algorithm

1
Hunan Provincial Key Laboratory of Health Maintenance for Mechanical Equipment, Hunan University of Science and Technology, Xiangtan 411201, China
2
College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082, China
3
National Innovation Center of Advanced Rail Transit Equipment, Zhuzhou 412000, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(23), 7948; https://doi.org/10.3390/en14237948
Submission received: 4 November 2021 / Revised: 17 November 2021 / Accepted: 24 November 2021 / Published: 27 November 2021

Abstract

:
The Broadband Mode Decomposition (BMD) method was previously proposed to solve the Gibbs phenomenon that occurs during photovoltaic signal decomposition; its main idea is to build a dictionary which contains signal features, and to search in the dictionary to solve the problem. However, BMD has some shortcomings; especially if the relative bandwidth of the decomposed signal is not small enough, it may treat a square wave signal as several narrowband signals, resulting in a deviation in the decomposition effect. In order to solve the problem of relative bandwidth, the original signal is multiplied by a high-frequency, single-frequency signal, and the wideband signal is processed as an approximate wideband signal. This is the modulation broadband mode decomposition algorithm (MBMD) proposed in this article. In order to further identify and classify the disturbances in the photovoltaic direct current (DC) signal, the experiment uses composite multi-scale fuzzy entropy (CMFE) to calculate the components after MBMD decomposition, and then uses the calculated value in combination with the back propagation (BP) neural network algorithm. Simulation and experimental signals verify that the method can effectively extract the characteristics of the square wave component in the DC signal, and can successfully identify various disturbance signals in the photovoltaic DC signal.

1. Introduction

The application of photovoltaic technology is becoming more and more widespread in daily life, and the generation industry of photovoltaic power is at a stage of rapid development. In photovoltaic systems, including DC and alternating current (AC) systems, there have been considerable economic benefits from the wide use of various related equipment such as frequency converters, inverters, rectifiers, and charging piles. However, in the power grid, many power quality problems such as harmonic, noise, and distortion issues will be produced due to the imbalance, impact, and nonlinearity of the load in the circuit. The use of electrical appliances and the stability of the power supply network will be seriously affected by these power quality problems [1,2,3]. Reference [4] proposes an evaluation scheme for renewable solar energy, which is aimed at the risk–benefit assessment of historical and traditional buildings. In terms of the photovoltaic direct current, the identification and classification of disturbance in photovoltaic DC signals is the key and basis for power quality analysis, evaluation, and management [5,6], therefore the classification and identification of power quality disturbance has high research value. Reference [5] proposed a comprehensive solution, which is a real-time analysis of power quality disturbances, and reference [6] proposed a new type of framework for the recognition of complex power quality disturbances; this framework is based on the multiple fusion convolutional neural network (MFCNN), which could automatically extract and fuse signal features from multiple sources.
In the engineering field, in order to extract features from signals that contain noise, many time–frequency analysis methods [7,8,9,10] are used. Reference [7] proposed a new method of emitter identification, which is based on spectral features and variational mode decomposition. In order to monitor family activities, reference [8] proposed a feature extraction method of acoustic signals, which is based on the non-negative matrix factorization (NMF) classifier. Reference [9] is a very classical algorithm, namely the empirical mode decomposition method (EMD), which was proposed by N. E. Huang et al. This method and the subsequent improved algorithm played a key part in signal feature extraction and signal decomposition. In order to extract time–frequency information features from vibration signals, wavelet packet transform was applied in reference [10]. These adaptive time–frequency analysis methods for signal analysis can be roughly divided into two categories according to whether they are based on the Fourier transform. For example, the variational mode decomposition (VMD) algorithm is based on the Fourier transform. All of these methods realize signal decomposition by computing in the frequency domain. The essence of the VMD algorithm is an adaptive filter [11], and the method was proposed by D. Zosso and K. Dragomiretskiy. VMD shows good performance [12,13,14] compared with previous methods, especially in the analysis of complex, non-stationary signals, which has been proved in many references. Reference [12] proposed a novel method of online chatter identification for milling processes. In order to make a prediction of short-term wind speed, reference [13] combined linear and nonlinear prediction models and proposed a hybrid method based on variational mode decomposition. To detect the blade imbalance fault in the marine current turbine (MCT), reference [14] proposed a novel VMD denoising method. The original signal was decomposed into several intrinsic mode functions, which was obtained by calculating the envelope of the extreme point; this is the essence of a method not based on the Fourier transform [15], such as the ensemble empirical mode decomposition (EEMD), the local mode decomposition (LMD), and the empirical mode decomposition (EMD) methods. The EEMD algorithm was proposed by Z. Wu et al., and compared with EMD, EEMD improves the anti-noise ability of the EMD method by averaging the decomposition results and mixing different levels of white noise into the original signal [16]. The VMD and EEMD methods are superior and effective [17,18], which can be shown by recent studies. Hence, the EEMD method was selected as a representative of the methods that are based on the non-Fourier transform, while the VMD method was selected as a representative of the methods that are based on the Fourier transform; these two algorithms were compared with the algorithm proposed in this paper.
Due to the incorporation and disconnection of linear loads in the circuit, an interference signal similar to a square wave signal is generated in the original photovoltaic DC signal, and this type of signal is generally sharp. Since signals with sharp corners, such as sawtooth signals and square wave signals, have frequency bands that are infinite and can be considered as wideband signals, the existing time–frequency analysis methods are expected show the Gibbs phenomenon when processing this kind of signal.
First, the essence of the algorithm based on the Fourier transform is multiscale adaptive filtering. However, after filtering, a series of interference may be caused at the breakpoint of the decomposition result; the characteristics of the sharp corner of the broadband signal disappearing or attenuating is known as the Gibbs phenomenon [19]. For broadband signals, the use of an improved Hilbert–Huang transform (HHT) and the energy entropy algorithm [20] to extract the features of alternating current (AC) signal square waves was proposed by K. He et al. Second, the essence of the methods that are based on the non-Fourier transform is the use of the interpolation function to calculate the envelope of the extreme point, and the division of the original signal into a number of intrinsic mode functions (IMFs) with a “smooth” narrow band. Therefore, when dealing with broadband components, errors will inevitably occur. Consequently, Y. Peng et al. proposed the Broadband Mode Decomposition (BMD) algorithm [21,22]. The innovative point of the BMD algorithm is to construct an association dictionary, which contains common broadband data such as narrowband data, square wave data, and sawtooth data. The sparse solution is then obtained by using the optimization method to search in the association dictionary. Compared with the previous signal decomposition algorithms, the BMD algorithm takes the adjustment differential operator as the optimal object; it can also extract signal features from complex signals, including narrowband signal features and broadband signal features. However, when applied to a wideband signal, the BMD algorithm may treat it as several narrowband components because its relative bandwidth of is not sufficiently small. Therefore, to denoise photovoltaic DC signals, a broadband mode decomposition (MBMD) method is proposed in this paper, which is based on the modulation differential operator. In order to make the decomposition result more accurate on the basis of the BMD algorithm, a high-frequency, single-frequency signal is multiplied by the MBMD so that the relative bandwidth of the effective wideband signal is converted into one that is sufficiently small and far lower than 1, and the wideband signal is then regarded as an approximate wideband signal.
Recently, J. Zheng et al. proposed the composite multiscale fuzzy entropy (CMFE) algorithm, which is based on the multiscale fuzzy entropy (MFE) algorithm. It has been proved to be suitable for processing non-stationary signals in engineering [23,24]. Therefore, it is applied to construct the feature vectors of the MBMD decomposition results. Then, the BP neural network algorithm is used to identify and classify the photovoltaic DC current disturbance.
The remaining sections of the paper are as follows. Section 2 introduces the disturbance state in the photovoltaic DC current and establishes a photovoltaic DC mathematical model. Section 3 introduces the problems in the identification of photovoltaic DC disturbances. Section 4 describes the details of the MBMD algorithm. Section 5 explains the other algorithms used. Section 6 describes the simulation analysis of the feature extraction of the photovoltaic signal disturbance model. Section 7 analyzes the collected experimental data sets of photovoltaic direct current signals. Section 8 provides the conclusion. All the data in this paper come from the data collected by the photovoltaic experimental platform in the literature [25].

2. Photovoltaic Electrical Signal Model

According to the analysis of the measured signals, the mathematical model of the DC signal is initially established, which consists of the DC signal, a harmonic signal superimposed on the DC, a distortion signal caused by the impact load, a noise signal caused by external interference, electromagnetic interference, etc.

2.1. Harmonic Signal

In normal operation, the current of the DC system contains abundant harmonic signals. The harmonic signal refers to the AC component superimposed on the DC. Therefore, the established harmonic signal model is as follows:
U 1 ( t ) = k U 0 ( t ) i = 1 n cos ( 2 π i ω t ) 0 < k < 0.1
where, U 0 is the amplitude of the DC current, k is the amplitude coefficient of harmonics, and n is the maximum number of harmonics.

2.2. Distortion Signal

Due to the impact, nonlinearity, and imbalance of the load, various disturbance components in the photovoltaic DC system increase. These disturbance components are distortion signals. For the photovoltaic DC system, there are two kinds of disturbances that have the greatest impact on photovoltaic power quality: current mutation and current pulse. The sudden change of current refers to the steep rise and fall of the current caused when the load switch in the system is opened and closed. For linear loads such as light bulbs, when the load switch is opened or closed, the current change is mostly straight up and down; for nonlinear loads such as motors, when the load switch is closed or opened, the current changes are characterized by slopes that rise and fall. This is determined by the characteristics of the nonlinear load.
Therefore, the linear load current mutation signal model is established as follows:
U 2 ( t ) = U 0 ( t ) ± A [ U 0 ( t t 1 ) U 0 ( t t 2 ) ]
where A is the amplitude coefficient of disturbance, t 1 is the initial time of disturbance, and t 2 is the end time of disturbance.
The nonlinear load current mutation signal model is the following:
U 3 ( t ) = U 0 ( t ) + [ k 1 × ( t t 1 + b ) 2 + a 1 ] + [ k 2 × ( t t 2 + b ) 2 + a 2 ] + A [ U 0 ( t t 1 b ) U 0 ( t t 2 + b ) ]
where k 1 and k 2 are the slope of the signal, a 1 and a 2 are the amplitudes of signal change, b is the duration of the signal rise or fall, t 1 is the initial time of disturbance, and t 2 is the end time of disturbance.
A current pulse is produced when faults such as a grounding fault or an automatic restart occur in the operation of the photovoltaic DC system, and refers to the current generally rising quickly and then falling quickly with a large amplitude change, but with a relatively short duration, generally 10–20 ms.
Therefore, the current pulse signal model is established as:
U 4 ( t ) = U 0 ( t ) + B × [ sin 2 π ω 2 ( t t 1 ) sin 2 π ω 2 ( t t 2 ) ]
where B is the amplitude coefficient of disturbance, t 1 is the initial time of disturbance, t 2 is the end time of disturbance, and ω 2 is the frequency of the pulse signal.

3. Problems in the Identification of Photovoltaic DC Disturbances

3.1. Gibbs Phenomenon

The Fourier series expands a signal into the sum of sinusoidal signals. Taking square wave signals as an example, the Fourier expansion of broadband signals is as follows:
s q u a r e ( t ) = i = 1 + 1 2 i 1 sin ( 2 i 1 ) t
Because the bandwidth of the filter is limited, the Gibbs phenomenon will appear after filtering. The Gibbs phenomenon is a phenomenon that may cause a series of disturbances at the breakpoint of the decomposition result, because the high frequency part of the broadband signal will attenuate or disappear after the filter is filtered. The wider the filter bandwidth, the more accurate the filter result is. Set s q u a r e ( ω t , r a t e ) to be the square signal, where rate is the duty cycle and ω is the angular frequency. In Equation (6), the mixed signal x 1 ( t ) includes a noise signal with a signal-to-noise ratio (SNR) of 10, a simple square signal, and an amplitude modulation and frequency modulation (AM–FM) sinusoidal signal. The time domain waveform of x 1 ( t ) and its components are shown in Figure 1.
x 1 ( t ) = s q u a r e ( 20 π t , 50 ) + 0.5 cos ( 10 π t + t 2 ) sin ( 300 π t ) + n ( t )
In order to obtain square wave components, the Butterworth bandpass filter is used to filter x 1 ( t ) , and the result after filtering is shown in Figure 2. Through experiments, the optimized filter parameters are determined (passband cutoff frequency is 48 hz, stopband cutoff frequency is 90 hz). As can be seen in Figure 2a, after filtering, there is an inevitable disturbance in the result due to the Gibbs phenomenon.

3.2. VMD Algorithm

By applying several Wiener filters, the VMD algorithm can adaptively divide the original signal into multiple components with limited frequency bands. By optimizing the variational model, the center frequency and bandwidth of the filter are obtained. The iterative process of the i-th IMF has the following three steps:
(1)
Calculate the analytical formula of I M F i ( t ) through the Hilbert transform.
[ δ ( t ) + j π t ] × I M F i ( t )
(2)
Estimate the center frequency ω i of the analyzed signal I M F i ( t ) , and convert the spectrum of the analyzed signal into the frequency baseband corresponding to ω i .
[ δ ( t ) + j π t ] × I M F i ( t ) e j w i t
(3)
Take the L2 norm of the analytic signal gradient as the optimization objective function, and take its bandwidth as the optimization parameter. Obtain the decomposition results by the alternating direction multiplier method (ADMM) in the optimization process. Describe the optimization problem as follows:
min { I M F i ( t ) } ,   { w i }   i   t [ δ ( t ) + j π t ] × I M F i ( t ) e j w i t 2 , S .   T .   x ( t ) = I M F i ( t )
When analyzing broadband signals, the VMD algorithm will be affected by the Gibbs phenomenon because its essence is the multi-scale adaptive Wiener filtering. The decomposition result of x 1 ( t ) in Equation (6) by VMD is shown in Figure 3. The curve has obvious disturbances, and the sudden information loss of the square wave signal is more serious.

3.3. EEMD Algorithm

EMD uses an interpolation algorithm to iteratively extract the mean value of the upper and lower envelopes and decompose the original signal into several IMFs.
(1)
Set h i ( t ) = x ( t ) , j = 1 and i = 1 .
(2)
Use the cubic spline function to interpolate the poles of h i ( t ) to obtain the lower envelope m l ( t ) and the upper envelope m u ( t ) , and then update h i ( t ) :
h i + 1 ( t ) = h i ( t ) m u ( t ) + m i ( t ) 2
(3)
Repeat step 2 until h i + 1 ( t ) meets the requirements defined by IMF; then, extract the j-th component:
I M F j ( t ) = h i + 1 ( t )
(4)
Set r ( t ) = x ( t ) I M F j ( t ) , update x ( t ) with the remaining r ( t ) , and then repeat the previous three steps until r ( t ) becomes a monotone function. The decomposition result can be expressed in the following form:
x ( t ) = j = 1 N I M F j ( t ) + r ( t )
In non-stationary signal processing, the EMD algorithm shows good performance. However, mode mixing is a major shortcoming of EMD; if the analyzed signal is a strong noise signal, this shortcoming may cause severe errors. Therefore, in order to improve the ability to resist mode mixing and noise interference, the EEMD algorithm is proposed. In contrast, the difference is that the EEMD algorithm averages the decomposition results and adds different levels of white noise to the original signal. After the averaging process of the EEMD algorithm, part of the noise in the original signal will be canceled out because the added white noise signals are uncorrelated.
The essence of the methods that are not based on the Fourier transform is to obtain the envelope of the extreme points by using interpolation functions to generate smooth envelopes such as the EMD algorithm and the EEMD algorithm. When dealing with broadband components, the final decomposition results will also become smooth because the IMFs are obtained by calculating the envelopes. The decomposed results corresponding to x 1 ( t ) in Equation (6) by applying EEMD are shown in Figure 4. It can be clearly seen that the IMF component obtained by EEMD is close to a sinusoidal signal; it is evident that the characteristics of the extracted square wave signal have changed.

4. MBMD Algorithm

4.1. Narrowband Signal

The definition of the narrowband signal is shown in Equation (13):
x n a r r o w ( t ) = A ( t ) cos ( ω t + ϕ ( t ) )
where ϕ ( t ) is function, and this function changes slowly. A ( t ) is limited by the frequency band, ω is its central frequency, which is far less than ω . The relative bandwidth is 2 ω / ω , and it is far lower than 1.
In Equation (13), the regulated difference operator of the BMD algorithm can effectively distinguish narrowband signals from noise because wideband or narrowband signals are always smoother than noise, and the operator D ( 2 ) can represent the smoothness of nonstationary signals. To constrain the generated components in order for them to be smoother than the original signal, the regulated differential operator is constructed in the BMD algorithm. The operator can almost be applied to broadband signals, although it is suitable for narrowband signals.

4.2. Broadband Signal

For broadband signals, the Fourier series of sawtooth signals is as follows, while the Fourier series of square wave signals is given by Equation (5).
s a w t o o t h ( t ) = i = 1 + ( 1 ) ( i 1 ) ( 2 i 1 ) 2 sin ( 2 i 1 ) t
According to Equations (5) and (14), as the frequency increases, the magnitude of the sinusoidal component decreases. Therefore, the general form of the broadband signals can be constructed as follows:
x b r o a d ( t ) = i = 1 + A i sin i ω t + θ i ( t )
when ω , A 0 , and A i decrease along with the increasement of ω , θ ( t ) changes slowly.
The frequency of the wideband signal is a multiple of the narrowband signal shown in Equation (15), but even when the frequency amplitude decreases rapidly, the broadband signal can be distinguished from the noise because the wideband signal is still smoother than the noise signal. However, because the modulation differential operator was originally designed for narrowband signals, there are problems with broadband signals. For x n a r r o w ( t ) , the first Fourier series A cos ( ω t + ϕ ( t ) ) multiplied by an envelope A ( t ) / A should be its theoretically generated component. As for the narrowband signal, there is only one dominant frequency ω .
However, there are multiple principal frequencies in x b r o a d ( t ) which have different amplitudes such as A i 1 sin ( i + 1 ) ω t + θ i 1 ( t ) and A i sin i ω t + θ i ( t ) . Among the sinusoidal components, the frequency interval is ω , but the center frequency of the broadband signal is also ω . The definition of relative bandwidth is the value obtained by dividing the bandwidth by the center frequency, which is ω / ω = 1 . In Section 4.1, the relative bandwidth of narrowband signals has been proven to be much less than 1. The BMD algorithm may regard A i sin i ω t + θ i ( t ) and A i 1 sin ( i 1 ) ω t + θ i 1 ( t ) as two different components because the relative bandwidth is not small enough, and then the algorithm may divide the broadband signal into several narrowband components when applied to a broadband signal with a strong noise interference.

4.3. Modulated Differential Operator for Broadband Signals

First, the broadband component is multiplied by the high-frequency sinusoidal signal as follows:
x b r o a d ( t ) = x b r o a d ( t ) sin ( ω t ) = i = 1 + A i sin i ω t + θ i ( t )   sin ( ω t )
where sin ( ω t ) is a high-frequency, single-frequency signal with ω = M ω , and M is a preset positive integer, which can usually be set to 5–10. x b r o a d ( t ) is a modulated signal, and it can then be converted to the following form:
x b r o a d ( t ) = 1 2 i = 1 + A i cos ( ω i ω ) t θ i ( t )   + 1 2 i = 1 + A i cos ( ω + i ω ) t + θ i ( t )
The center frequency of the broadband signal is moved from ω to ω ω and ω + ω . The relative bandwidths are much lower than 1, which are ω / ( ω ω ) and ω / ( ω + ω ) , respectively. After modulation, the relative bandwidth of a narrowband signal can be easily proved to also be much lower than 1. Therefore, the modulation difference operator can be constructed as follows:
T j = D ( 2 ) I M F i j ( t ) 2 2 + λ D ( 2 ) x ( t ) sin ( ω t ) I M F i j ( t ) 2 2
D ( 2 ) I M F i j ( t ) 2 2 is minimized to ensure the smoothness of the decomposed IMFs, and D ( 2 ) x ( t ) sin ( ω t ) I M F i j ( t ) 2 2 is minimized to rapidly reduce the smoothness of the residual to obtain the sparsest decomposition result of the original signal.

4.4. MBMD Algorithm

The main process of MBMD algorithm is outlined below, which is based on the modulation difference operator constructed in Equation (18). Figure 5 shows the flow chart of MBMD obtaining IMF components.
(1)
r 0 ( t ) is set to equal x ( t ) .
(2)
Separate direct current signals from the original signal.
I M F 0 ( t ) = i f f t [ r ^ 0 ( 1 ) ]
r 1 ( t ) = x ( t ) I M F 0 ( t )
where i is set to equal 1, r ^ 0 ( k ) is the FT of r 0 ( t ) , i f f t [ r ^ 0 ( 1 ) ] is the IFT of r ^ 0 ( 1 ) .
(3)
Construct the optimal problem P1 as follows:
P 1 :   M i n i m i z e   T 1 ( A 1 , ω 1 , θ 1 , D 1 ) , T 2 ( A 2 , ω 2 , θ 2 , D 2 ) , T 3 [ A 3 ( t ) , ω 3 ( t ) , θ 3 ( t ) ] S .   T .   x ( t ) = i = 0 N I M F i ( t ) + r e s ( t ) ,   I M F i j D i c j
T j = D ( 2 ) I M F i j ( t ) 2 2 + λ D ( 2 ) x ( t ) sin ( ω t ) I M F i j ( t ) 2 2
where D ( 2 ) is the quadratic differential form operator, T j is the modulated differential operator defined in Equation (18), λ > 0 is always set to equal 1. ω = M ω , ω equals the maximum value for the frequency of r i ( n ) . M can be generally set to 5–10.
(4)
Use ACROA to solve L1 to obtain the optimal value T j ( j = 1 , 2 , 3 ) , find the minimum T j and choose the optimal I M F i ( t ) of T j . Figure 5 illustrates the flow of this method.
(5)
Update r i ( t ) : r i + 1 ( t ) = r i ( t ) I M F i ( t ) after getting the best IMF.
(6)
Terminate the iteration process if the termination condition of Equation (23) is met; otherwise, return to step 3 and continue to separate the signal until the termination condition is met.
I M F i ( t ) I M F i 1 ( t ) 2 2 / I M F i 1 ( t ) 2 2 ε
Figure 5. Flow chart of MBMD to obtain IMF.
Figure 5. Flow chart of MBMD to obtain IMF.
Energies 14 07948 g005
The MBMD method is used to decompose the photovoltaic DC current signal, while the DC signal, pulse signal, distortion and harmonic signals and the noise signal are separated and extracted.

5. Introduction of Other Algorithms

CMFE Algorithm

In order to alleviate the shortening of the time series in the multi-scale fuzzy entropy (MFE) coarse-grained time series, CMFE applies the average MFE of different coarse-grained time series, which has the same proportion of factors as the extracted features. If the CMFE of a given signal is larger, it means that the given signal is more complex and contains more pattern information. The smaller the fluctuation range of the CMFE, the better the decomposition effect of an algorithm is and the more accurate the extracted signal features become. The calculation process of the CMFE algorithm is as follows:
(1)
Calculate a different coarse-grained time series y k ( τ ) = { y k , 1 ( τ ) , y k , 2 ( τ ) , K , y k , j ( τ ) } with a scale factor τ , for a given signal x ( n ) with N data points.
y k , j τ = 1 τ i = ( j 1 ) τ + k j τ + k 1 x i ,   1 j N / τ , 1 k τ
(2)
Calculate the fuzzy entropy of each coarse-grained time series y k τ , and then obtain the CMFE with the scale factor τ by averaging all the fuzzy entropies.
C M F E ( x ( n ) , τ ) = 1 τ k = 1 τ F u z z y E n ( y k ( τ ) )

6. Simulation Analysis

The simulation signal of a photovoltaic DC current disturbance is established, as shown in Figure 6. The second harmonic with an amplitude of 0.02 A and a frequency of 100 Hz is set in the signal, where disturbance A is the disturbance generated by linear load access or connection; disturbance B is the disturbance caused by nonlinear load access or connection; disturbance C is the current pulse disturbance, and n 1 ( t ) is white noise with an SNR of 25. The signal model formula is as follows:
f ( t ) = 2 + 0.02 cos ( 200 π t ) + A + B + B 1 [ 2 × ( t 2.5 ) 2 + 0.5 ] + B 2 [ 2 × ( t 4 ) 2 ]   + B 3 [ 2 × ( t 5.5 ) 2 0.5 ] + B 4 [ 2 × ( t 7 ) 2 ] + C sin ( 50 π t ) + n 1 ( t ) A = 2 0.2 t 0.7 1.5 1.0 t 1.7 0 t = e l s e   ,   B = 0.5 2.5 t 3.5 0.5 5.5 t 6.5 0 t = e l s e   B 1 = 1 2 t 2.5 0 t = e l s e   ,   B 2 = 1 3.5 t 4 0 t = e l s e B 3 = 1 5 t 5.5 0 t = e l s e ,   B 4 = 1 6.5 t 7 0 t = e l s e ,   C = 2.5 9 t 9.2 0 t = e l s e
When analyzing broadband signals, the VMD algorithm is still affected by the Gibbs phenomenon because its essence is the multi-scale adaptive Wiener filtering. Figure 7 shows the decomposition results of the photovoltaic DC simulation electrical signals by VMD.
It can be seen in Figure 7a that the mutation information of the square wave signal in the IMF1 component decomposed by VMD is lost, the Gibbs phenomenon appears, and the missing characteristic information appears separately in the remaining IMF components. Figure 7d is the residual component.
To obtain the envelope of the extreme points, the EMD method and the EEMD method typically use interpolation functions, which are non-Fourier transform-based methods, thereby generating a smooth envelope. When dealing with the broadband component, since the IMF is obtained by calculating the envelope, the final decomposition results are also smoothed. Figure 8 shows the decomposition result of the photovoltaic current simulation signal using EEMD. It can be seen that the obtained IMF1 and IMF2 are components with partial signal characteristics, the IMF3 component is close to a sinusoidal signal, and IMF4 is close to the DC component. Evidently, the extracted signal characteristics have changed.
It can be seen in Figure 9 that the MBMD method is used to decompose the DC current signal, which can separate and extract the DC signal, pulse signal, distortion and harmonic signals, and the noise signal.
Table 1 lists the accuracy parameters of the separated IMF obtained by the three algorithms, including the accuracy parameters of the correlation coefficient r 1 , the energy error E 1 , and the time T to make a more accurate comparison.
In order to compare the effect of signal decomposition under different signal-to-noise ratios, two additional simulation models are constructed here. Except for noise, the other components remain unchanged. White noise n 2 ( t ) with an SNR of 10 and white noise n 3 ( t ) with an SNR of 5 were added, respectively. The two signal models constructed are shown in Figure 10 and Figure 11.
Figure 12 and Figure 13 show the VMD decomposition results of these two signal models. It can be seen that with the decrease of SNR, the components of the VMD containing noise increase, especially IMF2 and IMF3. Due to the increase of noise, the decomposed components also contain part of the noise signals, so that the characteristics of the square wave signal are not obvious. Figure 14 and Figure 15 are the decomposition results of EEMD. Some components obtained by EEMD are relatively smooth, and the increased noise components are reflected in IMF1. Figure 16 and Figure 17 shows the decomposition results of MBMD. Table 2 and Table 3 show the relevant parameters (energy error E 1 , correlation coefficient r 1 , and time T), respectively. It can be seen that the decomposition effect of MBMD is the best under different SNRs.

7. Experimental Analysis

The flow chart of experimental analysis is roughly shown in Figure 18:

7.1. Data Collection

The experimental photovoltaic data were collected from the experimental platform for small photovoltaic measurement and the control system. As shown in Figure 19, the instrument used in the experiment has two 18 V, 50 W solar panels of 1480 mm × 680 mm, and the solar panels store electrical energy in batteries. The rated voltage of the battery as a stable power source is 12 V. The controller of this photovoltaic system functions as a combiner box and a controller at the same time, which can control the charging and discharging of the battery and solar panels. The inverter parameters are 12 V, 1000 W. The linear loads used are 36 W AC bulbs and 50 W DC bulbs. The nonlinear loads used are 8 W AC motors, 28 W AC motors, and 20.4 W DC motors. The signal measurement module contains four sensors; among them, the DC sensor uses a Hall sensor, while the rest are ordinary current and voltage sensors. These signals are collected by our self-made signal acquisition software, a specific collection software on the computer. The sampling time is 20 s, and the signal acquisition frequency is 2 kHz. In the experiment, the switching sequence of each load was set first, and 50 sets of signals, including the no-load state, were obtained. After collecting the signals, the 15th set of signals was selected as the experimental signals.

7.2. Experimental Analysis of Feature Extraction

Select the 15th set of experimental signals. Figure 20 shows the signal collected through the photovoltaic experimental platform; for convenience of comparison, Figure 21a, Figure 22a and Figure 23a are all current signals collected in the experiment. It can be seen in Figure 20 that the transient rise and fall of the current in the original experimental signal caused by the switch of the DC bulb is similar to a regular rectangular signal, which is a broadband signal. Moreover, as the DC motors are turned off, two sets of current pulse signals are generated.
The 15th set of collected signals is decomposed by the VMD algorithm. It can be seen in Figure 21b that after the current distortion generated by the DC load switch is decomposed, the mutation information of the Chinese wave signal of the IMF1 component of the VMD decomposition is lost because of the Gibbs phenomenon; the mutation information and the current pulse signal are thus divided into a set as narrowband frequency signals.
The EEMD algorithm is used to decompose the collected 15th set of signals. It can be seen in Figure 22, among the IMF components obtained by the EEMD decomposition method, that I M F 3 is the closest to the original current signal. When processing broadband components, the final decomposition result will also become smooth because the IMF is obtained by calculating the envelope.
The MBMD algorithm was used to decompose the collected signals of set 15, as shown in Figure 23. The MBMD algorithm divides the experimental signals into: Figure 23b, DC component; Figure 23c, broadband component; Figure 23d, current distortion component (narrowband signal); Figure 23e, second harmonic component; Figure 23f, harmonic and noise component. Since the second harmonics of the photovoltaic DC signal account for a large part of the signal, and the third harmonics and other high harmonics have very little content, the third harmonics and other high harmonics along with the noise components are divided together. It can be seen that the MBMD algorithm has a good decomposition effect on the feature information of the broadband signal, and is suitable for extracting effective features.
In order to compare the decomposition effect of the three algorithms of MBMD, VMD, and EEMD on the photovoltaic DC experimental signal more closely, CMFE was used to calculate the three algorithms. In addition, 40 sets of four common photovoltaic DC disturbance signals were collected in the experiment. The MBMD, VMD, and EEMD algorithms were used to decompose each disturbance, and the IMF1 components obtained by each of the three algorithms were extracted by CMFE. To unify the standard, the scale factor of the composite multi-scale fuzzy entropy (CMFE) was set to 4, and a graph with a scale factor of 3 was used for comparison. The CMFE comparison chart of the three algorithms under the four disturbance states is presented in Figure 24, Figure 25, Figure 26 and Figure 27.
Because harmonic disturbance always exists in the photovoltaic DC circuit, the disturbance states measured by the four experiments all contain harmonics. From the CMFE diagram above, it can be seen that among the various resolution methods, the MBMD method has the smallest CMFE fluctuation after decomposition, followed by VMD, and EEMD has the largest fluctuation; therefore, the MBMD decomposition has the best effect. In order to effectively and accurately distinguish the type of photovoltaic DC signal disturbance, the BP neural network algorithm was used to extract and classify the results after CMFE extraction.

7.3. Disturbance Identification

At this point, 40 sets of photovoltaic DC signal samples of the following four states were selected from the collected experimental signals. Because harmonics always exists in the photovoltaic DC circuit, and although it is also a kind of disturbance, it is marked as state 1; the other three states have a disturbance in addition to harmonics. State 1 is the harmonics; State 2 is the current distortion of the DC linear load plus harmonics; State 3 is the current distortion of the DC non-linear load plus harmonics; State 4 is the current pulse plus harmonics.
Comparing the CMFE decomposed by MBMD from Figure 25 and Figure 26, it can be seen that the magnitude of the CMFE decomposed by the two disturbance states is very close. In order to prevent the BP neural network algorithm from mixing these two types of signals into one state, the first three IMF components decomposed by different algorithms are calculated by CMFE, and the CMFE of the three IMF components is used as the standard to classify the disturbance state. To illustrate that the CMFE of the three IMF components can be classified, we used the MBMD decomposition algorithm as an example: We took 40 sets of sample signals for each state (160 sets of signals in total) and decomposed the signals in the four disturbance states in the samples. Subsequently, we took the first three IMF components obtained by decomposition, calculated the respective CMFE values, and also calculated the value range as well as the mean and range of the CMFE. Table 4 shows the CMFE values of the three IMF components obtained by MBMD decomposition. It can be seen that four states can be distinguished by comparing the values of IMF1, IMF2, and IMF3.
According to the data in Table 4, the CMFE value range of IMF1 in State 2 is 0.0049–0.0180, and the CMFE value range of IMF1 in State 3 is 0.0126–0.0305, with an intersection of 0.0126–0.0180. If only the CMFE value of IMF1 is used for classification, the disturbance categories within this intersection range cannot be correctly classified, hence three IMF components are used for classification. The CMFE value of IMF2 in State 2 ranges from 0.0021 to 0.0043, while that of IMF2 in State 3 ranges from 0.0311 to 0.2566. Therefore, even if IMF1 in States 2 and 3 in the intersection cannot be classified successfully, it can be classified by IMF2. By observing and comparing IMF1, IMF2, and IMF3 in the four states in the table, it can be concluded that the four states can be distinguished by these value classes. For comparison, the parameters of VMD and EEMD are calculated, as shown in Table 5 and Table 6, respectively.
By observing and comparing the parameters of the three tables, it can be seen that the variance of CMFE calculated by MBMD is the smallest, followed by VMD, and the variance of EEMD is the largest. For the characteristic value of a state, it must be that the smaller the fluctuation, the smaller the variance is, which is better; in addition, the four figures in the previous section, Figure 24, Figure 25, Figure 26 and Figure 27, also showed the superiority of MBMD in extracting characteristic values. Therefore, the CMFE values of the three IMF components obtained by MBMD were used as samples, and the BP neural network algorithm was used to identify their disturbance categories.
The BP neural network is a feedforward network; it is multi-layered, and it is trained through error back propagation. In order to minimize the error sum of squares of the network, it continuously adjusts the weight and threshold of the network through backpropagation and uses the fastest descent method. We took 30 sets of samples in each state as training samples, and 10 sets as test samples. The BP neural network algorithm was used to classify the test samples. The BP classification diagrams of the three signal decomposition algorithms are shown in Figure 28, Figure 29 and Figure 30. The classification accuracy rate of each of the four state classifications and the parameters of the total classification accuracy rate are shown in Table 7.
It can be seen in Figure 28 that the predicted classification results of the 40 test samples are highly consistent with the actual categories, while Figure 29 and Figure 30 still illustrate a large number of classification errors. Table 7 lists the specific classification parameters. It can be seen that the MBMD algorithm only misclassified one case in the classification of State 1, and the classification of other states was 100% successful, with a total accuracy of 97.50%; on the other hand, the VMD and EEMD algorithms had an unsatisfactory classification effect of only about 70%. It can be seen that the classification effect of these two algorithms for State 1, State 3, and State 4 is still good, both of which are above 80%. However, when State 2 is classified, a large number of errors occur, with only 20% and 10% accuracy, respectively. The main reason is that State 2 is a harmonic plus square wave signal, and the two algorithms are not ideal in extracting the features of the other wave signal due to the defects caused by the Gibbs phenomenon or extreme point interpolation. The experimental results show that the BP neural network combined with the MBMD algorithm can effectively identify the disturbance types in photovoltaic DC signals.

8. Conclusions

The photovoltaic DC current signal and its disturbance were taken as the research objects, and a series of algorithm experiments on feature extraction and classification was conducted. A new method of identifying the disturbance state of the photovoltaic DC current is consequently proposed.
(1)
The categories and causes of photovoltaic DC current signal and its disturbance are proposed to be studied, and the corresponding mathematical model is then constructed. This mathematical model can be used for the simulation of photovoltaic DC signals and other papers about the photovoltaic DC current.
(2)
The MBMD method is proposed to solve the problem of broadband signal feature extraction due to the Gibbs phenomenon in the photovoltaic DC signal disturbance analysis. In past time–frequency methods such as EEMD, VMD, and other defects caused by extreme point interpolation or the Gibbs phenomenon, it is proposed that the MBMD method can be avoided by searching in the dictionary. The experimental and simulation results verify that the method can effectively extract the characteristics of broadband components in photovoltaic DC signals. This algorithm can denoise, especially for DC current signals containing square wave signals, and the denoising effect is better than other algorithms.
(3)
Combined with the MBMD algorithm, CMFE, and the BP neural network algorithm, a set of recognition methods for photovoltaic DC signal disturbance are proposed. The results show that the recognition accuracy of this method is high, and it can effectively identify the disturbance in the photovoltaic DC signal. It can be used as a special method to detect and identify disturbance types in photovoltaic DC signals.
(4)
When analyzing complex signals composed of wideband components, the application of MBMD can be extended to signal processing in other fields. However, because the optimization process is more complicated, MBMD takes more time than EEMD and VMD. Therefore, in the future, we will study some fast optimization methods for MBMD and some other improved algorithms.

Author Contributions

Algorithm proposal and innovation, Z.W. and Y.P.; simulation experiment and data analysis, Z.W. and H.G.; experimental supervisors, Y.L. (Yanfei Liu) and Y.G.; writing—original draft preparation, Z.W.; equipment debugging and data management, Y.L. (Yi Liu), S.L., and C.F.; writing—review and editing, Z.W., Y.P. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China under Grant No. 51805161, the Research Project of Hunan Provincial Department of Education under Grant No. 19C0769, the Hunan Provincial Natural Science Foundation of China under Grant Nos. 2021JJ30260 and 2018JJ3182, the Scientific Research Fund of Hunan Provincial Education Department under Grant No. 19B191, and the Hunan Innovative Province Construction Special Foundation under Grant No. 2019RS2062.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dash, P.K.; Satapathy, P.; Nayak, P.; Sahani, M. Islanding and power quality disturbance monitoring in microgrid using adaptive cross variational mode decomposition and reduced kernel ridge regression. Int. Trans. Electr. Energy Syst. 2020, 30, e12364. [Google Scholar] [CrossRef]
  2. Leite, D.; Decker, L.; Santana, M.; Souza, P. EGFC: Evolving gaussian fuzzy classifier from never-ending semi-supervised data streams—With application to power quality disturbance detection and classification. In Proceedings of the 2020 IEEE International Conference on Fuzzy Systems (FUZZ-IEEE), Glasgow, UK, 19–24 July 2020; pp. 1–9. [Google Scholar] [CrossRef]
  3. Mishra, M. Power quality disturbance detection and classification using signal processing and soft computing techniques: A comprehensive review. Int. Trans. Electr. Energy Syst. 2019, 29, e12008. [Google Scholar] [CrossRef] [Green Version]
  4. López, C.S.P.; Lucchi, E.; Leonardi, E.; Durante, A.; Schmidt, A.; Curtis, R. Risk-benefit assessment scheme for renewable solar solutions in traditional and historic buildings. Sustainability 2021, 13, 5246. [Google Scholar] [CrossRef]
  5. Chen, Z.; Li, M.; Ji, T.; Wu, Q. Real-time recognition of power quality disturbance-based deep belief network using embedded parallel computing platform. IEEJ Trans. Electr. Electron. Eng. 2020, 15, 519–526. [Google Scholar] [CrossRef]
  6. Qiu, W.; Tang, Q.; Liu, J.; Yao, W. An automatic identification framework for complex power quality disturbances based on multi-fusion convolutional neural network. IEEE Trans. Ind. Inf. 2020, 16, 3233–3241. [Google Scholar] [CrossRef]
  7. Satija, U.; Trivedi, N.; Biswal, G.; Ramkumar, B. Specific emitter identification based on variational mode decomposition and spectral features in single hop and relaying scenarios. IEEE Trans. Inf. Forensics Secur. 2019, 14, 581–591. [Google Scholar] [CrossRef]
  8. Lee, S.; Pang, H. Feature extraction based on the non-negative matrix factorization of convolutional neural networks for monitoring domestic activity with acoustic signals. IEEE Access 2020, 8, 122384–122395. [Google Scholar] [CrossRef]
  9. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  10. Yen, G.Y.; Lin, K.C. Wavelet packet feature extraction for vibration monitoring. IEEE Trans. Ind. Electron. 2000, 47, 650–667. [Google Scholar] [CrossRef] [Green Version]
  11. Dragomiretskiy, K.; Zosso, D. Variational mode decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  12. Yang, K.; Wang, G.; Dong, Y.; Zhang, Q.; Shang, L. Early chatter identification based on an optimized variational mode decomposition. Mech. Syst. Signal Process. 2019, 115, 238–254. [Google Scholar] [CrossRef]
  13. Ali, M.; Khan, A.; Rehman, N.U. Hybrid multiscale wind speed forecasting based on variational mode decomposition. Int. Trans. Electr. Energy Syst. 2018, 28, e2466. [Google Scholar] [CrossRef]
  14. Wei, J.; Xie, T.; Wang, T. A VMD denoising-based imbalance fault detection method for marine current turbine. In Proceedings of the IECON 2020—46th Annual Conference of the IEEE Industrial Electronics Society, Marina Bay Sands Expo and Convention Centre, Singapore, 18–21 October 2020; pp. 1–6. [Google Scholar] [CrossRef]
  15. Sweeney, K.T.; McLoone, S.F.; Ward, T.E. The use of ensemble empirical mode decomposition with canonical correlation analysis as a novel artifact removal technique. IEEE Trans. Biomed. Eng. 2013, 60, 97–105. [Google Scholar] [CrossRef] [Green Version]
  16. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–59. [Google Scholar] [CrossRef]
  17. Yi, C.; Wang, D.; Fan, W.; Tsui, K.; Lin, J. EEMD-based steady-state indexes and their applications to condition monitoring and fault diagnosis of railway axle bearings. Sensors 2018, 18, 704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Hao, Y.; Du, Z.; Xing, Z. Leakage signal analysis of urban gas pipeline based on improved variational mode decomposition. Int. J. Pattern Recognit. Artif. Intell. 2020, 34, 1–16. [Google Scholar] [CrossRef]
  19. Ridder, F.D.; Pintelon, R.; Schoukens, J.; Verheyden, A. Reduction of the gibbs phenomenon applied on nonharmonic time base distortions. IEEE Trans. Instrum. Meas. 2005, 54, 1118–1125. [Google Scholar] [CrossRef]
  20. He, K.; Zhang, Z.; Xiao, S.; Li, X. Feature extraction of ac square wave saw arc characteristics using improved hilbert–huang transformation and energy entropy. Measurement 2013, 46, 1385–1392. [Google Scholar] [CrossRef]
  21. Peng, Y.; Li, Z.; He, K.; Liu, Y.; Lu, Q.; Li, Q.; Liu, L.; Luo, R. Broadband mode decomposition and its application to the quality evaluation of welding inverter power source signals. IEEE Trans. Ind. Electron. 2020, 67, 9734–9746. [Google Scholar] [CrossRef]
  22. Peng, Y.; Li, Z.; He, K.; Liu, Y.; Lu, Q.; Li, Q.; Liu, L.; Luo, R. Quality monitoring of aluminum alloy DPMIG welding based on broadband mode decomposition and MMC-FCH. Measurement 2020, 158, 107683–107694. [Google Scholar] [CrossRef]
  23. Zheng, J.; Jiang, Z.; Pan, H. Sigmoid-based refined composite multiscale fuzzy entropy and t-SNE based fault diagnosis approach for rolling bearing. Measurement 2018, 129, 332–342. [Google Scholar] [CrossRef]
  24. Zheng, J.; Pan, H.; Cheng, J. Rolling bearing fault detection and diagnosis based on composite multiscale fuzzy entropy and ensemble support vector machines. Mech. Syst. Signal Process. 2017, 85, 746–759. [Google Scholar] [CrossRef]
  25. Wang, Z.; Peng, Y.; Liu, Y.; He, K.; Wu, Z.; Li, B. The calculation method of pv direct current energy based on modulated broadband mode decomposition and compound simpson integral algorithm. IEEE Access 2021, 9, 51403–51415. [Google Scholar] [CrossRef]
Figure 1. Time domain signals. (a) Mixed signals x 1 ( t ) ; (b) Square wave component; (c) AM–FM sine component; (d) Noise composition.
Figure 1. Time domain signals. (a) Mixed signals x 1 ( t ) ; (b) Square wave component; (c) AM–FM sine component; (d) Noise composition.
Energies 14 07948 g001
Figure 2. Filtering result of x 1 ( t ) . (a) Square wave component; (b) Residues.
Figure 2. Filtering result of x 1 ( t ) . (a) Square wave component; (b) Residues.
Energies 14 07948 g002
Figure 3. Decomposition results of VMD for x 1 ( t ) . (a) Square wave components corresponding to IMF; (b) IMF corresponding to the AM–FM sine component; (c) Noise component.
Figure 3. Decomposition results of VMD for x 1 ( t ) . (a) Square wave components corresponding to IMF; (b) IMF corresponding to the AM–FM sine component; (c) Noise component.
Energies 14 07948 g003
Figure 4. Decomposition results x 1 ( t ) by EEMD. (a) Square wave components corresponding to IMF; (b) IMF corresponding to the AM–FM sine component. (c) Noise component.
Figure 4. Decomposition results x 1 ( t ) by EEMD. (a) Square wave components corresponding to IMF; (b) IMF corresponding to the AM–FM sine component. (c) Noise component.
Energies 14 07948 g004
Figure 6. Photovoltaic DC simulation signal diagram.
Figure 6. Photovoltaic DC simulation signal diagram.
Energies 14 07948 g006
Figure 7. VMD simulation signal decomposition results. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Figure 7. VMD simulation signal decomposition results. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Energies 14 07948 g007
Figure 8. EEMD simulation signal decomposition results. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Figure 8. EEMD simulation signal decomposition results. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Energies 14 07948 g008
Figure 9. Simulation signal decomposition diagram of MBMD. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Figure 9. Simulation signal decomposition diagram of MBMD. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Energies 14 07948 g009
Figure 10. Photovoltaic signal simulation diagram under the noise of SNR 10.
Figure 10. Photovoltaic signal simulation diagram under the noise of SNR 10.
Energies 14 07948 g010
Figure 11. Photovoltaic signal simulation diagram under the noise of SNR 5.
Figure 11. Photovoltaic signal simulation diagram under the noise of SNR 5.
Energies 14 07948 g011
Figure 12. VMD decomposition diagram under the noise of SNR 10. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Figure 12. VMD decomposition diagram under the noise of SNR 10. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Energies 14 07948 g012
Figure 13. VMD decomposition diagram under the noise of SNR 5. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Figure 13. VMD decomposition diagram under the noise of SNR 5. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) Residue.
Energies 14 07948 g013
Figure 14. EEMD decomposition diagram under the noise of SNR 10. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Figure 14. EEMD decomposition diagram under the noise of SNR 10. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Energies 14 07948 g014
Figure 15. EEMD decomposition diagram under the noise of SNR 5. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Figure 15. EEMD decomposition diagram under the noise of SNR 5. (a) IMF1 component; (b) IMF2 component; (c) IMF3 component; (d) IMF4 component; (e) Residue.
Energies 14 07948 g015
Figure 16. MBMD decomposition diagram under the noise of SNR 10. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Figure 16. MBMD decomposition diagram under the noise of SNR 10. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Energies 14 07948 g016
Figure 17. MBMD decomposition diagram under the noise of SNR 5. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Figure 17. MBMD decomposition diagram under the noise of SNR 5. (a) DC signal; (b) Wideband signal; (c) Narrowband signal; (d) Harmonic signal; (e) Noise signal; (f) A 0.1 s signal extract from (d).
Energies 14 07948 g017
Figure 18. Flow chart of experimental analysis.
Figure 18. Flow chart of experimental analysis.
Energies 14 07948 g018
Figure 19. Experimental platform for small photovoltaic measurement and the control system. A: Measuring system: ① Analog-to-digital conversion module; ② Signal measurement module; ③ Computer terminal signal acquisition software; B: Photovoltaic system: ④ Photovoltaic system controller; ⑤ Inverter; ⑥ Storage battery (12 V); C: DC loads: ⑦ DC motor (20.4 W); ⑧ DC bulb (50 W); D: AC loads: ⑨ AC motor (28 W); ⑩ AC motor (8 W); ⑪ AC bulb (36 W); E: Switches: ⑫ Load switches; ⑬ Master switch.
Figure 19. Experimental platform for small photovoltaic measurement and the control system. A: Measuring system: ① Analog-to-digital conversion module; ② Signal measurement module; ③ Computer terminal signal acquisition software; B: Photovoltaic system: ④ Photovoltaic system controller; ⑤ Inverter; ⑥ Storage battery (12 V); C: DC loads: ⑦ DC motor (20.4 W); ⑧ DC bulb (50 W); D: AC loads: ⑨ AC motor (28 W); ⑩ AC motor (8 W); ⑪ AC bulb (36 W); E: Switches: ⑫ Load switches; ⑬ Master switch.
Energies 14 07948 g019
Figure 20. Experimental signal diagram.
Figure 20. Experimental signal diagram.
Energies 14 07948 g020
Figure 21. VMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) IMF1 component; (c) IMF2 component; (d) IMF3 component; (e) Residue.
Figure 21. VMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) IMF1 component; (c) IMF2 component; (d) IMF3 component; (e) Residue.
Energies 14 07948 g021
Figure 22. EEMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) IMF1 component; (c) IMF2 component; (d) IMF3 component; (e) IMF4 component; (f) Residue.
Figure 22. EEMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) IMF1 component; (c) IMF2 component; (d) IMF3 component; (e) IMF4 component; (f) Residue.
Energies 14 07948 g022
Figure 23. MBMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) DC signal; (c) Broadband signal; (d) Narrowband signal; (e) Harmonic signal; (f) Noise signal; (g) A 0.1 s signal extract from (d).
Figure 23. MBMD decomposition diagram of the experimental signal. (a) Experimental current signal; (b) DC signal; (c) Broadband signal; (d) Narrowband signal; (e) Harmonic signal; (f) Noise signal; (g) A 0.1 s signal extract from (d).
Energies 14 07948 g023
Figure 24. MBMD, VMD, EEMD comparison diagram of the harmonic signal.
Figure 24. MBMD, VMD, EEMD comparison diagram of the harmonic signal.
Energies 14 07948 g024
Figure 25. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the square wave signal.
Figure 25. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the square wave signal.
Energies 14 07948 g025
Figure 26. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the narrowband signal.
Figure 26. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the narrowband signal.
Energies 14 07948 g026
Figure 27. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the pulse signal.
Figure 27. MBMD, VMD, EEMD comparison diagram of the harmonic signal plus the pulse signal.
Energies 14 07948 g027
Figure 28. The BP classification diagram of the value of the components decomposed by MBMD calculated by CMFE.
Figure 28. The BP classification diagram of the value of the components decomposed by MBMD calculated by CMFE.
Energies 14 07948 g028
Figure 29. The BP classification diagram of the value of the components decomposed by VMD calculated by CMFE.
Figure 29. The BP classification diagram of the value of the components decomposed by VMD calculated by CMFE.
Energies 14 07948 g029
Figure 30. The BP classification diagram of the value of the components decomposed by EEMD calculated by CMFE.
Figure 30. The BP classification diagram of the value of the components decomposed by EEMD calculated by CMFE.
Energies 14 07948 g030
Table 1. Evaluation parameters of the three algorithms (SNR 25).
Table 1. Evaluation parameters of the three algorithms (SNR 25).
Methodr1E1T(s)
VMD0.95430.01594.7134
EEMD0.80520.58359.3256
MBMD0.99810.001534.3517
Table 2. Evaluation parameters of the three algorithms (SNR 10).
Table 2. Evaluation parameters of the three algorithms (SNR 10).
Methodr1E1T(s)
VMD0.93140.03144.9525
EEMD0.78370.74939.7264
MBMD0.99790.001735.5524
Table 3. Evaluation parameters of the three algorithms (SNR 5).
Table 3. Evaluation parameters of the three algorithms (SNR 5).
Methodr1E1T(s)
VMD0.91970.05375.2497
EEMD0.75290.86229.9633
MBMD0.99770.001936.1146
Table 4. CMFE evaluation parameters of MBMD in four states.
Table 4. CMFE evaluation parameters of MBMD in four states.
StateIMFCMFE of MBMD
Value RangeMean ValueVariance
State 1IMF11.8781–3.53632.61170.1686
IMF20.5583–0.72610.64840.0023
IMF31.0725–1.27321.13430.0027
State 2IMF10.0049–0.01800.01031.5732 × 10−5
IMF20.0021–0.00430.01075.5222 × 10−5
IMF30.1245–0.32330.17440.0019
State 3IMF10.0126–0.03050.02135.0052 × 10−5
IMF20.0311–0.25660.07840.0040
IMF31.1450–1.54701.26870.0109
State 4IMF10.1281–0.23850.17170.0010
IMF20.3203–0.66470.58490.0068
IMF30.0598–0.59880.14300.0078
Table 5. CMFE evaluation parameters of VMD in four states.
Table 5. CMFE evaluation parameters of VMD in four states.
StateIMFCMFE of VMD
Value RangeMean ValueVariance
State 1IMF10.9901–2.67821.70520.2026
IMF20.5742–2.62460.73390.1657
IMF30.9360–15.66474.474914.7161
State 2IMF10.0028–0.04910.01217.3699 × 10−5
IMF20.0056–1.81550.86550.6079
IMF30.0179–3.05921.77960.3520
State 3IMF10.0353–0.11840.05355.1415 × 10−4
IMF21.6848–3.25172.44040.2222
IMF30.0306–3.92501.36621.3549
State 4IMF10.2456–0.70480.41280.0074
IMF20.0048–2.66491.39390.5016
IMF30.0179–3.05921.77960.3520
Table 6. CMFE evaluation parameters of EEMD in four states.
Table 6. CMFE evaluation parameters of EEMD in four states.
StateIMFCMFE of EEMD
Value RangeMean ValueVariance
State 1IMF11.8885–7.72242.41291.0112
IMF21.8442–4.48972.49840.3383
IMF30.9122–19.35054.263916.1200
State 2IMF10.0014–0.08020.02523.5337 × 10−4
IMF20.0016–2.59090.92280.7383
IMF30.0133–4.24641.76780.4866
State 3IMF10.0408–0.18390.10880.0012
IMF20.0345–4.13192.35240.5414
IMF30.0411–4.47431.46651.7213
State 4IMF10.0029–0.74390.41280.0168
IMF20.0024–3.64051.41890.7393
IMF35.0444–9.45366.86230.5020
Table 7. BP classification parameters of the three algorithms.
Table 7. BP classification parameters of the three algorithms.
StateParameterAlgorithm
MBMDVMDEEMD
State 1Error point112
Accuracy90%90%80%
State 2Error point089
Accuracy100%20%10%
State 3Error point011
Accuracy100%90%90%
State 4Error point012
Accuracy100%90%80%
TotalError point11114
Accuracy97.50%72.50%65.00%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Peng, Y.; Liu, Y.; Guo, Y.; Liu, Y.; Geng, H.; Li, S.; Fan, C. Photovoltaic Power Quality Analysis Based on the Modulation Broadband Mode Decomposition Algorithm. Energies 2021, 14, 7948. https://doi.org/10.3390/en14237948

AMA Style

Wang Z, Peng Y, Liu Y, Guo Y, Liu Y, Geng H, Li S, Fan C. Photovoltaic Power Quality Analysis Based on the Modulation Broadband Mode Decomposition Algorithm. Energies. 2021; 14(23):7948. https://doi.org/10.3390/en14237948

Chicago/Turabian Style

Wang, Zucheng, Yanfeng Peng, Yanfei Liu, Yong Guo, Yi Liu, Hongyan Geng, Sai Li, and Chao Fan. 2021. "Photovoltaic Power Quality Analysis Based on the Modulation Broadband Mode Decomposition Algorithm" Energies 14, no. 23: 7948. https://doi.org/10.3390/en14237948

APA Style

Wang, Z., Peng, Y., Liu, Y., Guo, Y., Liu, Y., Geng, H., Li, S., & Fan, C. (2021). Photovoltaic Power Quality Analysis Based on the Modulation Broadband Mode Decomposition Algorithm. Energies, 14(23), 7948. https://doi.org/10.3390/en14237948

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

Article Metrics

Back to TopTop