The Green Photovoltaic Industry Installed Capacity Forecast in China: Based on Grey Relation Analysis, Improved Signal Decomposition Method, and Artificial Bee Colony Algorithm

,


Introduction
In recent years, the development of new energy industries has gradually become a focus of countries all over the world. e rapid development of the economy has become more and more dependent on energy demands, not only in the pursuit of high efficiency, low cost, and stabile energy but also to improve the cleanliness and environmental friendliness of energy. Among many clean energy sources, solar energy is widely distributed with large energy values and has high development potential, representing new growth potential for the world's economy. According to the National Bureau of Statistics, in 2018, China's annual power generation reached 6.8 trillion kWh, of which thermal power generation accounted for 75.08%, while photovoltaic power generation was only 0.0539 trillion kWh, accounting for only 0.79%. China's photovoltaic power generation still has large room for improvement. In the "13th Five-Year Plan" for China's national power development, by 2020, the installed capacity of photovoltaic power generation is planned to reach 110 million kilowatts [1]. In 2017, the newly installed capacity of photovoltaic power generation in China was 53 million kilowatts, a year-to-year increase of 53.62% [2]. After continuous construction throughout the year, by the end of 2017, the total installed capacity of photovoltaic power generation in China reached 130 million kilowatts, and the "13th Five-Year Plan" was exceeded in advance.
From the perspective of PV installation, China has become the world's largest PV application market. Except for the occasional impact of the international market, the global PV growth rate has slowed. e remainder of the annual increase rate is above 100%, which indicates that the development center of the photovoltaic market is moving from Europe to emerging markets like China, Japan, and the United States. In 2018, China's photovoltaic power generation had a new installed capacity of 10.6 million kilowatts, accounting for about 20% of the new installed capacity of the global photovoltaic market, and continues to maintain a good development trend [3]. However, China's PV industry is greatly affected by the import and export trade situation, and the domestic PV market is still immature. erefore, it is reasonable to develop China's PV capacity and installed PV capacity, especially for the spread of centralized and distributed PV installed capacity [4]. Installation that is too rapid will lead to an increase in the phenomenon of "abandonment of light" and a serious overcapacity; a slow installation speed and insufficient installed capacity will lead China to miss the growth forecast of its economic development. erefore, a high-precision PV installed capacity forecasting model is proposed. Meanwhile, the innovations of this paper are as follows: (1) In order to accurately predict the cumulative installed capacity of photovoltaic power in China, this paper proposes a new support vector forecasting model based on improved empirical mode decomposition and an artificial bee colony algorithm to predict the cumulative installed capacity of photovoltaics in China. (2) Since the initial parameters of the traditional LSSVM model include the regularization parameter c and the width g of the radial basis function, the ABC (artificial bee colony) algorithm can solve the chattering and premature problems of LSSVM parameter optimization and determine the optimal weight and threshold. (3) e traditional LSSVM algorithm cannot fully determine the time-frequency characteristics of the time series data and thus greatly affects the performance of the forecasting model. en, the accuracy of the LSSVM cannot be ensured. e raw data are separated by CEEMD (complete ensemble empirical mode decomposition), which expands the magnitude of the data to achieve the data requirements of the machine algorithm. Based on the above reasons, this paper applies two different algorithms, CEEMD and ABC, into the LSSVM algorithm.
e main contents of this article are as follows: Section 2 introduces recent related literature; Section 3 introduces the mathematical principle of the signal decomposition model, including the artificial bee colony algorithm and the LSSVM algorithm, as well as the flow of the overall forecasting algorithm; Section 4 concludes the flow chart of China's cumulative PV installed capacity based on the CEEMD-ABC-LSSVM forecasting model; Section 5 applies the proposed CEEMD-ABC-LSSVM model to the prediction of cumulative PV installed capacity in China. By comparing the forecasting results of LSSVM, ABC-LSSVM, and EMD-ABC-LSSVM, the experimental results prove the practicality and accuracy of the proposed model in predicting China's cumulative PV installed capacity; Section 6 presents several forward-looking conclusions.

Literature Review
In recent years, many scholars have studied the photovoltaic industry's development system, which is mainly composed of two major parts-the industrial main body and the industrial environment [5]. Among them, the main body of the industry refers to the main production factors within the photovoltaic industry. e industrial environment mainly includes factors such as the level of social and economic development, the relevant policy system of the photovoltaic industry, market demand, and industrial capital investment, which affect the development of the industry. According to the system framework, the photovoltaic industry's development system is divided into two subsystems: the photovoltaic industry subsystem and the industrial development environment subsystem. Xiao et al. [6] argued that in the photovoltaic (PV) technology value chain, improving technology is an important factor in reducing PV costs. Investigating the policies of Germany and China proves that the deployment support and investment equations of PV part manufacturing plants affect the PV industry's important factors. In [7], Adaramola studied the technical and economic performance of an 80 kW solar photovoltaic gridconnected power generation system. In the case of solar photovoltaic power generation costs and global average solar radiation, the proportion of solar photovoltaic power generation and its average energy costs were obtained. Finally, it is economically feasible to develop solar photovoltaic systems for grid-connected power generation in northeastern Nigeria. Perez et al. [8] proposed a model predictive control (MPC) method to manage the energy generated by gridconnected photovoltaic (PV) power plants with energy storage in real time to optimize their economic benefits. Jafarzadeh et al. [9] believes that the randomness of solar energy resources is one of the factors affecting their largescale diffusion in the power system. e system for solar power plant modeling uses fuzzy logic and the predicted fuzzy system uses type-1 and interval-2 Takagi-Sugeno-Kang (TSK) fuzzy systems. e results of this example show that the type-2 TSK model provides the best performance index based on the type-2 prefactor and clear results.
PV installed capacity prediction is a complex nonlinear solution problem. So far, many scholars have proposed various prediction models, such as grey theory [10,11], multiple regression [12,13], and time series [14,15]. In recent years, various intelligent algorithms and signal decomposition models have also been widely used in installed capacity prediction [16][17][18][19][20][21][22][23][24][25]. Yuan et al. [10] defined the approximate real starting point for a steady-state growth trend sequence with random perturbations and calculated the generalized grey relation analysis (GRA) on this basis. e main reasons for the expansion of the photovoltaic industry were analyzed. It was found that the global installed capacity of photovoltaics is the main factor affecting the development of the photovoltaic industry. Shi et al. [11] proposed a periodic relationship model using the horizontal and vertical coordinates of the curve for the traditional grey correlation model. Ultimately, this model only refers to the wave period but is independent of it. e swing can show positive and negative correlations. Wolff et al. [12] evaluated multiple numerical weather prediction (NWP) parameters to improve their ability to improve photovoltaic power generation prediction capabilities. e importance of functionality is determined by a random forest algorithm. Furthermore, by using support vector regression, the random forest and linear regression models performed PV power predictions to test the resulting highest ranked features. Raghoebarsing and Reinders [13] outlined the current state of Suriname's energy sector and photovoltaic (PV) systems and examined the role that PV systems can play in the country's future energy transition. According to the government and nongovernment agencies that have planned photovoltaic projects, it is predicted that future PV systems will generate 0.8 TWh of electricity to meet the power demands of 2027. Sun et al. [14] proposed a new predictive model for multistep wind speed prediction (WSF). A twostage signal decomposition method combining set empirical mode decomposition (EEMD) and variational mode decomposition (VMD) was used to decompose empirical wind speed data. Further, the parameters in the proposed wavelet neural network (WNN) model were optimized to improve the prediction performance, which was realized by a hybrid backtracking search, and its effectiveness was illustrated. Pavlyshenko [15] studied the use of machine learning models in sales forecast analysis and proved that when there is a small amount of historical data in a specific sales time series, this effect can be used to make sales forecasts. e results prove that machine learning can improve sales time.
is shows the performance of predictive models of sequence prediction.
Signal decomposition and machine intelligence algorithms have been widely used in the field of hybrid prediction. Hu et al. [16] proposed a hybrid prediction method that includes EWT (empirical wavelet transform), CSA (coupled simulated annealing), and LSSVM (least squares support vector machine) to improve the accuracy of shortterm wind speed predictions. e proposed hybrid model has been shown to predict the average half-hour wind speed series obtained from a windmill farm in northwestern China. Wang et al. [17] proposed a variational mode decomposition (VMD) based on phase space refactoring (PSR) and the genetic algorithm. e results show that Wang et al.'s model is superior to all other comparison models, which proves that this model has superior multistep performance. Li et al. [18] developed four time series prediction techniques, including the metabolic grey model (MGM), the autoregressive integrated moving average (ARIMA), the grey model (GM)-ARIMA, and the nonlinear metabolic grey model (NMGM). A more accurate prediction of installed capacity can help to develop a better energy sustainability strategy. Tsaur and Lin [19] defined a solar building as a building-affiliated photovoltaic (BAPV) system and proposed a technology acceptance model (TAM) to predict and explain the public acceptance of BAPV. Wu et al. [20] considered the output radiation, cloud cover, and other weather conditions as the main influencing factors affecting PV production, considering the short-term predictions of operational plans, exchange sources, and peak load matching. Five types of prediction modules have been developed, including ARIMA, SVM, ANN, ANFIS, and a combined model using genetic algorithms. e prediction results show that the combined model has higher precision and efficiency.
LSSVM is a practical machine learning method that has been widely used in many fields. In order to ensure performance, it is crucial to choose the right technology to obtain the optimized hyper parameters of the LSSVM algorithm. Wu and Peng [21] proposed a bat algorithm-(BA-) optimized least squares support vector machine (LSSVM) hybrid model to improve the prediction of wind energy production. In order to effectively select the input of the LSSVM, the parameters in the LSSVM are optimized by the BA to verify the learning ability and generalization ability of the LSSVM. e results of the example show that the accuracy of the prediction model proposed in the paper can be improved by about 20%. Xiang et al. [22] proposed a decomposition-based artificial bee colony algorithm (ABC) to deal with multiobjective optimization problems (MAOPs). eir example shows that this new algorithm has good convergence and diversity. e quality of its solution is high, and its speed is fast. Zhao et al. [23] used the least squares support vector machine (LSSVM) method to establish the relationship between the reliability index and design variables and used the artificial bee colony (ABC) algorithm for reliability optimization. e numerical results show that this method not only satisfies the design accuracy but also improves optimization of the reliability parameters, such as the mechanical parameters, ground stress, and internal pressure of the rock mass. In [24], Xiao et al. proposed a fault diagnosis method based on heterogeneous information fusion in order to accurately determine fault type based on obtained information. e fusion feature was then introduced into the artificial bee colony algorithm-(ABC-) optimized least squares support vector machine (LSSVM). Simulation experiments show that the proposed method has higher diagnostic accuracy. Mustaffa and Yusof [25] used an enhanced artificial bee colony (EABC) to obtain the ideal values for the hyper parameters of the LSSVM, namely, the regularization parameters, the gamma parameters, and the kernel parameters, σ (2). en, the LSSVM was used as the forecasting model. is model was applied to the prediction of financial time series data, and the simulation results ensure the accuracy of parameter selection.

Principles of Forecasting Model
3.1. LSSVM. Given a set of training data samples (x i , y i ) N , among which x i ∈ R m is the m-dimensional data samples and y i ∈ R is a sample output, the LSSVM optimization algorithm is as follows [26]: where φ(x i ): R m ⟶ R mf is the mapping function from the original space to a high dimensional space, w ∈ R mf is the weight vector, e ∈ R is the tolerance, b is the offset, and c is the normalized coefficient. According to the objective function and constraints, we establish the Lagrange function as where a i is the LaGrange factor. For the KKT conditions, the L derivative can be obtained as After cancelling w and e, we obtain the following matrix equation: where I is the identity matrix, e optimal decision function is

Artificial Bee Colony Algorithm.
A swarm of bees can adapt to most living environments on the planet. We divided the bee colony into different types of work (such as collecting honey, breeding, and defense) to ensure a division of labor and cooperation. At the same time, there is a special "information transmission mechanism" between bees [27,28]. e bees communicate information to inform other bees around them and distribute work. e ABC algorithm is a swarm intelligence optimization algorithm designed by simulating the process of bees constantly searching for excellent honey sources. e ABC algorithm is a highly efficient and simple intelligent algorithm that is suitable for many function optimization problems. e honey source position in the algorithm represents a feasible solution for the target optimization function. e optimal honey source represents the optimal solution of the function optimization problem. e corresponding relationship between honey bee behavior and the optimization function is shown in Table 1.
In the ABC algorithm, the honey bee collecting process is divided into four stages: the initialization stage, the hired bee stage, observing the bee, and the detection bee stage [29].
(1) Initialization phase: the initial honey source location is randomly generated within the feasible solution domain. ese locations represent the initial source of honey (a feasible solution) and calculate their fitness values. (2) Hired bee stage: hired bees search the neighborhood at the current honey source location and record better honey source information. e recorded honey source information is then transmitted to the observation bee through the "swing dance." When the search for the honey source is completed and a better source of honey cannot be found, the hired bee gives up the honey source and turns into a detective bee.
(3) Observing the bee: this stage entails observing the bee waiting in the "swinging dance area" and receiving the honey source information transmitted by the employed bee. After receiving the information, the observation bee decides whether to choose to follow the hire bee to the honey source. If it chooses to follow, the observation bee turns into a hired bee. (4) Detection bee stage: if the employer is unable to find a better source after multiple searches, the hired bee will abandon the current honey source and convert into a detection bee. e scouting bee then begins to randomly search for a better source of honey. When a better source of honey is found, the scout is converted back to a hired bee.
During the entire search process, the bees constantly change between these three different roles until the end of the algorithm. us, during the operation of the algorithm, the three bees constantly change roles, and the possibility of the algorithm finding the global optimal solution is greatly improved. At the same time, the speed of the calculation method is also improved. e bees adjust their respective divisions of labor through the transformation of their roles and complete the transformation of their roles through the sharing of information.
e transformation diagram is shown in Figure 1. 4 Mathematical Problems in Engineering e ABC algorithm is a very simple and easy-to-calculate intelligent algorithm. e specific steps for its implementation are as follows [23].
(1) Initialize the population: when the algorithm starts running, the number of bee colonies is np and the number of hired bees is SN. Randomly generate NP honey sources, which are usually set as SN � NP/2. e feasible domain dimension of the algorithm is D. e maximum number of iterations of the algorithm is maxcycle, and the maximum number of stagnations is limit. Within the feasible domain of the objective function, we randomly generate the hired bees' SN honey source locations according to formula (6), which is recorded as X � X j i | i � 1, 2, . . . , SN; j � 1, 2, . . . , D}, and calculate the fitness of the honey source fitness according to formula (7): fitness � 1 1 + fitness , where rand(0, 1) refers to a random number representing a uniform distribution between (0, 1). D is the dimension of the optimized space, the value of k is selected from D 1, 2, . . . , SN { }, and k ≠ i, X j i ≠ X j k . maxcycle is used to specify the maximum number of iterations of the algorithm. limit is used to specify the maximum number of times a hired bee stays in a single source of honey.
(2) Hired bee search phase: this stage involves hiring a bee to conduct a honey source search to find the candidate a new honey source formula according to formula (8). e new honey source location will update according to formula (9): where φ ij is a random number between [− 1, 1], k ∈ (1, 2, . . . , SN) and X ij ≠ X kj . e old and new honey sources are selected by the greedy rule of formula (9); that is, the fitness values of the new and old honey sources is compared according to their fitness, and then a larger value of honey is taken. If the fitness value fitness of the new honey source V ij is better than the old honey source X ij , we use the new honey source V ij to replace the location of the old honey source, X ij . Otherwise, we save honey source location X ij , which is still the current honey source location. At the same time, the number of stagnating bees in the honey source is increased by one. (3) Observation bee search phase: at this stage, the observation bee selects the honey source according to the roulette method based on the honey source information obtained by the hired bee to perform the honey source search. If we choose to follow the hired bee, we search for the new honey source according to the method of step 2. e probability of the observation phase is as follows: (4) Detection bee search stage: each time the number of the stagnation of the honey source changes, a judgment is made. If this number reaches the threshold limit and the bee still does not find a better source, the bee gives up the honey source of the current search. e bee then turns into a scouting bee and begins a random search until a new, better source of honey is found. (5) Optimal solution judgment: in this stage, we judge whether the stop criterion is reached and then output the optimal honey source position, fitness value, and other parameters; otherwise, we continue to run step 2.

LSSVM Optimized by ABC.
Since the initial parameters of the traditional LSSVM model include the regularization parameter c and the width g of the radial basis function, the ABC algorithm can solve the chattering and premature problems of the LSSVM parameter optimization and determine the optimal weight and threshold. In the ABC algorithm, the fitness  Mathematical Problems in Engineering function in the ABC algorithm is calculated; the initial population size and the maximum evolution algebra maxgen are set, and the genetic operations, such as selection, improved crossover, and mutation, are performed on the individuals in the population. Finally, for the global excellent fitness, we use optimal fitness to obtain the regularization parameter c and the width g of the radial basis function. e flow chart of LSSVM optimized by ABC is shown in Figure 2.

Basic Principle of CEEMD.
In 1998, Huang et al. [27] proposed Empirical Mode Decomposition (EMD). is analysis method can decompose any complex signal according to different time-frequency scales and then generate a series of intrinsic mode function components (IMFs) from different frequencies (from a high frequency to a low frequency).
Huang et al. proposed three hypotheses: (1) Any complex signal consists of several M components (2) e number of extreme points of each IMF component is consistent with the number of zeros, and the envelope signal diagram is symmetrical around the time axis (3) Any signal can be decomposed into several M components at any time, and IMF full component superposition reconstruction can completely recover the signal.
In 2010, Yeh et al. [29] made some improvements to the EMD algorithm. By adding positive and negative paired auxiliary noise, the residual auxiliary noise in the reconstructed signal can be eliminated and the computational efficiency can also be improved. is method is called complete ensemble empirical mode decomposition (CEEMD). e steps of the CEEMD algorithm are as follows: (1) Firstly, fixed-intensity white noise is added into the original signal; then, the noise-added signal is subjected to EMD decomposition to obtain an IMF component, and n different kinds of white noise are continuously decomposed for corresponding times. en, the obtained components are averaged: where IMF 1 is the first-order of IMF component, n is the number of times the noise is added, e i indicates the i-th component of the EMD decomposition, ε indicates the proportion of noise added, and ω i is the white noise added. (2) Next, we calculate the residual by subtracting the first-order IMF component: (3) en, we decompose r i + εe 1 [ω i ], i � 1, 2, . . . , n and continue to iteratively decompose until the component needed to satisfy IMF 1 emerges. en, we obtain the overall average of the whole IMF 1 . In this way, the second-order of the original signal IMF 2 is obtained, which is en, take their overall average to obtain the k + 1th-order component of the original signal IMF k+1 , which is (5) Continue iterative screening and stop when the number of extreme values of the residual does not exceed two:

Mathematical Problems in Engineering
where r is the residual value and k is the number of IMF components. (6) e original signal recovered from each IMF component and the residual sequence can be expressed as Compared with EMD, CEEMD not only solves the modal aliasing problem perfectly but also has higher computational efficiency and complete signal decomposition. erefore, the CEEMD method is more suitable for filtering denoising.
After the signal is decomposed by CEEMD, several IMF components with high to low frequency are obtained from small to large, including IMF 1 , IMF 2 , . . ., IMF k , IMF k+1 . e high-frequency IMF component with a small order generally contains more noise, and the noise of the low-frequency IMF component with a large order has little effect. ere must be a component in the middle of the IMF component IMF k as the demarcation point. e components between IMF k+1 to IMF n refer to the signal dominant zone, and the other is the dominant zone of noise. In 2007, Boundraa and Cexus [30] proposed a method based on continuous mean square error to determine the demarcation point for IMF k . e discriminant formula of this method is as follows: where y k (t i ) is the selected signal reconstructed from the k-th layer to the residual, y k+1 (t i ) is the selected signal reconstructed from the k + 1-th layer to the residual, e k is the energy, and n is the signal length.

The Forecasting Model of CEEMD-ABC-LSSVM
Due to the late development of China's PV market, the historical data of the installed capacity are insufficient, and it is difficult to realize the data requirements of intelligent algorithms. erefore, it is very easy to create homogeneity among the influencing factors of the original data as well as different scale information characteristics for the input data. e traditional LSSVM algorithm cannot fully determine the time-frequency characteristics of the time series data, which greatly affects the performance of the forecasting model. us, the accuracy of the LSSVM cannot be ensured. e raw data are separated by CEEMD, which expands the magnitude of the data to achieve the data requirements of the machine algorithm. Based on the above reasons, this paper applies two different algorithms (CEEMD and ABC) into the LSSVM algorithm and proposes a Chinese PV installed capacity forecasting model called the CEEMD-ABC-LSSVM combined forecasting model. e forecasting steps of the overall algorithm are as follows: (1) Time series data decomposition. Use the CEEMD algorithm to decompose x(t) and then obtain n IMF components and one residual r n .

Data Processing.
Based on [1-3, 5-9], this paper selects the industry prosperity index (the annual average of the electric power production and supply industry), the cost of power generation (the average value of the photovoltaic ongrid electricity price for centralized and distributed general projects), and the domestic PV Output value (100 million), PV module production, foreign market demand, domestic market demand, international PV trade policy easing, PV tariff, GDP, population, total social electricity consumption,   Mathematical Problems in Engineering photovoltaic industry investment, photovoltaic power generation, macro policy support, and photovoltaic power generation technology as preliminary indicators affecting China's PV installed capacity. Among these indicators, relaxation of international PV trade policy, macro policy support, and photovoltaic power technology maturity are qualitative indicators. rough the Delphi method, experts score the data. Since China's PV installed capacity is a complex nonlinear problem with many influencing factors, the relationship between each index and PV installed capacity is vague. erefore, this paper uses the grey correlation method to define the distance between the influencing factors and the installed capacity of PV, as well as its initial impact. e effective collection of indicator sets can effectively improve the speed and accuracy of the forecasting model. In addition, based on the same data test set, Matlab R2014a is used for programming. is test was performed on a computer with an Intel Core i5-6200U, 8G memory, and Windows 10 Professional Edition as a support.

Normalization of Indicators.
If any data are missing for some years, in order to eliminate the dimensional difference between the data of different metrics and to ensure the availability of data, the data should first be rationally adjusted and supplemented. Normalize all the corrected data and use the Z-score data normalization method to standardize the N sets of data in the data sets of the M indicators: where Z mn is the normalized data, x mn is the corrected data, x n is the average of x mn , and σ is the standard deviation of x mn . Since the cost of photovoltaic power generation is a negative index, in order to achieve the forwardization of data, the inverse index forward processing method is adopted: where X mn refers to the data after forwardization, max is the maximum value of the normalized data, and min refers to the minimum value of the normalized data.

Grey Relational Analysis.
In this paper, the normalized data for 1990-2018, which provide the preliminary influencing indicator for China's PV installed capacity forecast, are used as the input for the GRA. en, the degree between each impact index and China's PV installed capacity can be calculated, as shown in Table 2.
Exclude the indicators below 0.5500. en, the domestic market demand, PV tariffs, macro policy support, PV industry investment, and photovoltaic power technology maturity are eliminated, and the main factors affecting China's PV installed capacity can be collected. e main factors include total social electricity consumption, GDP, power generation cost, population, photovoltaic power generation, PV module production, industry prosperity index, domestic PV output value, international PV trade policy easing, and foreign market demand. e normalized index values of the main factors affecting China's PV installed capacity (1990-2018) are shown in Table 3.

PV Installed Capacity Forecasting Based on the CEEMD-ABC-LSSVM Model.
e original historical load sequence is decomposed by CEEMD. Input the signal sequence of China's historical PV installed capacity and then load the CEEMD model with the whole forecasting algorithm. With this process, we obtain 4 IMFs and 1 residual. e decomposition results are shown in Figure 4.
It can be seen from Figure 4 that the time series of China's PV installed capacity has clear multiscale features, and the four IMF components present information for high and low variations with different fluctuation scales. Among them, the frequency of IMF 1 is relatively high, which reflects the random noise information of the original time series data. Meanwhile, the residual frequency is low and the change is stable. e trend information reflecting the timing of China's PV installed capacity reflects the overall change characteristics of China's PV installed capacity. e number of IMFS obtained by CEEMD decomposition is related to the characteristics of the original data time series itself. Based on the nature of these IMF components, we can understand the practical significance of CEEMD decomposition. e original Chinese PV installed capacity sequence and each of the IMF components decomposed from the original data sequence are independent of each other, and the respective IMF components are mutually orthogonal. e last item obtained by decomposition removes the remainder of the original sequence from all IMF components. Indeed, this item represents the average trend of the original Chinese PV installed capacity change.
e influential factor indicator set as well as each IMF component and residual are used as the inputs for the training set and the test set of the ABC-LSSVM model. In this way, the predictive model is trained. After many training cycles, the accuracy of the test set forecasting result increases significantly. e main adaptive and nonadaptive parameters of ABC-LSSVM multiple training cycles are shown in Table 4.
We reconstructed the forecasting series to obtain the forecasting results and relative errors of the components and residuals of the time series decomposition of China's PV installed capacity since 1990, as shown in Figure 5.
Based on an analysis of Figure 5, we find that the proposed CEEMD-ABC-LSSVM forecasting model has high forecasting accuracy for each IMF and residual prediction sequence of the PV installed capacity, and the     relative error of prediction is less than 0.8. Among them, the relative error of the residual signal timing sequence is less than 4%, while the residual frequency is low, and the change is stable. e trend information reflecting the timing of China's PV installed capacity reflects the overall variation characteristics of China's PV installed capacity. erefore, the CEEMD-ABC-LSSVM forecasting model proposed in this paper has been proven to be highly scientific and applicable.

Model Comparison Analysis.
In order to compare the superiority of the proposed model, this paper also uses data for the main factors affecting China's PV installed capacity in 2019-2040, such as LSSVM, ABC-LSSVM, EMD-ABC-LSSVM, and CEEMD-ABC-LSSVM. Algorithm training was performed on the four forecasting models, and the degree of fitting for the four forecasting models is shown in Figure 6.
In order to objectively compare the accuracy of a variety of models, common statistical indicators, including RMSE, R 2 , and MRE, are adopted, and the index calculation formula is as follows: where q i is a predicted value, q i is a sample mean, and n is a sample number. e calculation results for the four models are compared and shown in Table 5.
From this analysis, we can draw some conclusions, as follows: (1) According to the calculation results of the index, the accuracy of the signal decomposition forecasting algorithm is higher than that of the single algorithm (2) e index values of RMSE and MRE are lower than 2 and 50, respectively, indicating that signal decomposition plays an important role in the combined forecasting algorithm (3) e three indicators of the CEEMD-ABC-LSSVM combined forecasting model proposed in this paper are better than those of EMD-ABC-LSSVM, which indicates that the proposed forecasting model has better applicability to China's PV installed capacity forecasting problem Another important indicator of the machine intelligence algorithm for prediction is the operation speed of the algorithm. erefore, this paper compares the four algorithms of LSSVM, ABC-LSSVM, EMD-ABC-LSSVM, and CEEMD-ABC-LSSVM after 100 training cycles. e average operation time of each set is shown in Table 6.
rough this comparison, we can draw some conclusions based on the experimental results, as follows: (1) e ABC-LSSVM algorithm reduces the number of repetitive process executions of the LSSVM algorithm, which needs to select the optimal parameters and takes up fewer resources. us, the algorithm enables faster training speeds; because of its greater reliance on computation, this algorithm takes longer to perform calculations during the training process. is aspect makes ABC-LSSVM weaker than the other two forecasting models.
(2) However, the combined forecasting model of CEEMD-ABC-LSSVM proposed in this paper is similar to the average training time of EMD-ABC-LSSVM, and its training accuracy is much higher than that of the latter, which reflects the superiority of the proposed model.
(3) e accuracy comparison of the forecasting model has also proven that the forecasting model proposed in this paper will not fall into the local optimum. e superiority of the proposed combination algorithm has been proven based on both its accuracy and operation speed.

Prediction Results.
We applied GM (1, 1) to predict the total electricity consumption, GDP, power generation cost, population, photovoltaic power generation, PV module production, industry prosperity index, domestic PV output value, international PV trade policy easing, and foreign market demand from 2019 to 2035, which was used as the input data for the CEEMD-ABC-LSSVM forecasting model. Finally, we examined China's cumulative photovoltaic installed capacity from 2019 to 2035, as shown in Figure 7.

Conclusion
is article combined the CEEMD algorithm with the ABC-LSSVM algorithm to construct a CEEMD-ABC-LSSVM forecasting model for China's PV installed capacity. After selecting the main influential factors, we used data from 1990 to 2010 as the training set for the CEEMD-ABC-LSSVM, EMD-ABC-LSSVM, ABC-LSSVM, and LSSVM forecasting models and data from 2011 to 2018 as the test set. en, we compared the results from the forecasted results of all training and testing sets. Finally, we applied the CEEMD-ABC-LSSVM forecasting model to forecast China's future PV installed capacity in 2019-2035. We find that China's PV installed capacity will surpass 4000 GW around 2035. And, its installed capacity will increase year by year. Moreover, China's PV industry development will maintain a steady overall growth rate (compound annual growth rate). Our final conclusions and policy recommendations are as follows: (1) From 2019 to 2035, forecasts of China's photovoltaic industry show a gradual upward trend, except for growth rate. From the perspective of market certainty, although the government has a corresponding PV installed capacity index for every year, internal and external uncertainty factors are faced by the photovoltaic industry. At the same time, international PV market volatility also produces continuous market uncertainty. (2) China's PV installed capacity will be transferred to a scale effect stage around 2030. At that time, excessive PV industry investment will yield sustained GDP growth. In the future, China will continue to develop renewable energy technologies. Meanwhile, China must pursue better opportunities in the energy investment field and continue to reduce its pollutants and carbon emissions.
(3) China must pay more attention to the development of energy consumption terminal power to successfully realize its energy conservation and emission reduction responsibilities to the international community. As China's photovoltaic power generation industry gradually enters the era of affordable Internet access, the cost of photovoltaics on grid electricity will experience a gradual reduction. Indeed, China's PV industry still has significant room for future development. China's PV installed capacity will gradually become the main force among the global PV installed capacity.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
All the authors were involved in the data collection and preprocessing phase, model constructing, empirical research, result analysis and discussion, and manuscript preparation. All the authors have approved the submitted manuscript.