Research Article  Open Access
Application of Improved PSOBP Neural Network in Cold Load Forecasting of Mall AirConditioning
Abstract
A combination of JMP, PSOBP neural network, and Markov chain which aims at the low correlation between input and output data and the error of prediction model in the PSOBP neural network prediction model is proposed. First, the JMP data processing software is used to process the input data and eliminate the samples with low coupling degree. Then, obtaining the cooling load prediction results relies on the training from the PSOBP neural network. Finally, the final prediction results will be generated by eliminating the random errors using the Markov chain. The results show that the combination of the prediction methods has higher prediction accuracy and conforms to the change rule of the cooling load in shopping malls. Besides, the combination fits the actual application requirements as well.
1. Introduction
At present, the energy consumption of airconditioning occupies an increasing proportion of whole building energy consumption [1]. Its huge power consumption has increased the pressure on the grid. Some researchers have used ice storage airconditioning to solve this problem [2]. The operation of ice storage airconditioning requires a reasonable matching with the cooling capacity of the internal cooler and the ice trough in the peaks and troughs of the electricity price. Firstly, the cooling capacity of the building at every moment in the next day is predicted, and then the cooling capacity should be divided on the basis of the predicted results. Therefore, the dynamic prediction of the cooling load is the core of ice storage airconditioning.
For the airconditioning cooling load prediction, most researchers use the datadriven method [3], mainly including support vector machine (SVM), statistical regression, decision tree, genetic algorithm, and neural network algorithm. SVM is a common artificial intelligence method [4], which has the ability to transform nonlinear relations, but it takes too long to deal with data [5, 6]. Statistical regression is a simple and feasible prediction method, but its predictive power is lower than SVM, and it is difficult to choose the predictors of this method [7]. The decision tree method is a technique which is easy to understand and divides data into groups using a tree diagram, but its prediction results often deviate significantly from the actual situations [8], and it cannot address time series and nonlinear data very well [9]. Genetic algorithm is a powerful optimization algorithm for processing complex model problems [10]. When the input data is large and complex, the algorithm can obtain appropriate solutions according to objective functions or subjective judgments, but it also has disadvantages, such as nonunique results and the long calculating time [11, 12]. As a kind of prediction algorithms, neural network is widely used in various fields [13]. One of the neural networks BP (back propagation) is widely used for its powerful nonlinear mapping, selflearning, generalization, and fault tolerance, but the weaknesses of this method are localized minimization and slow convergence. Some researchers have used PSO (particle swarm optimization) [14] to optimize the initial threshold and weight of the BP neural network and to improve the convergence speed and accuracy of BP. However, due to the low matching degree between the input data and the output data and the errors under the process of combining model, the prediction results cannot meet the ideal requirements.
Based on the traditional artificial neural network, this paper shows the results from JMP data analysis software, which analyzes the correlation between input variables and output variables and eliminates irrelevant data. In addition, we used the Markov chain to correct the errors of the artificial neural network. The research shows that the improved model has higher prediction accuracy and convergence speed than those of the traditional artificial neural network prediction method, which is more suitable for the practical application of the project.
2. Thoughts on Dynamic Prediction Method of AirConditioning Cooling Load
The key and difficult point of the cooling load forecasting is to improve the prediction accuracy and timeliness. The traditional prediction of the airconditioning cooling load in shopping malls does not consider the marketrunning characteristic points, and it cannot distinguish the differences between weekends and business days [15]. Therefore, based on the characteristics of the mall, the data should be reanalyzed and integrated before forecasting. This specific prediction idea is shown in Figure 1.
According to the content above, JMP data analysis software is used to analyze the input data correlation to improve the neural network. The characteristic of the analysis of correlation can eliminate the input and output variables with low coupling degree. It can also fix the problems about slow training speed, and lowprecision vulnerability results from the low correlation between input and output variables. Therefore, the prediction accuracy of the neural network can be improved, and it can be applied to the dynamic prediction of cooling load of airconditioning in shopping mall buildings. In consideration of some relative errors in the predicted results of the combined model, the Markov chain is adopted to further modify the prediction results so that the accuracy can be improved.
3. Construction of Dynamic Prediction Model for Cooling Load of PSOBP Large Shopping Mall
3.1. The Prediction Model of the PSOBP Neural Network
The PSOBP neural network is used to predict the cooling load of airconditioning in shopping malls. In fact, the PSO algorithm is used to iteratively optimize the weight and threshold of the BP neural network, and then the weight and threshold with the best fitness are obtained for the prediction of the BP neural network.
In order to predict the cooling load of airconditioning in shopping malls, this paper adopts a threelayer network architecture. The commonly used input data for cold storage airconditioning cold load prediction include the outdoor dry bulb temperatures T and T−1 at the current and previous moments, respectively, and the solar radiation doses R, R1, and R2 at the current, previous, and the first two moments. At the current time, the outdoor humidity is H, the outdoor wind speed is W, and the previous time, the first two times, and the first four times of the cold load are C1, C2, and C4, respectively, with the current time cold load C as the output. In the input layer, the number of neurons is 10 and the output layer is the cold load at time T. The empirical formula (1) can be obtained for nodes of hidden layer of the neural network as 21:where A is the number of neurons in the hidden layer and B is the number of neurons in the input layer.(1)To figure out the fitness of the particles. The fitness value of the particles is determined by the predictionrelated data and the data obtained by the PSO algorithm, as shown in the following formula: where is the number of samples, is the observed value of sample , and is the predicted value of sample .(2)The fitness of the current particle is selected as the best fitness of the current particle and the best fitness of the historical optimum. The best historical fitness of the current particle and the best global fitness are selected as the historical optimal fitness of the current particle.(3)The velocities and positions of the particles are updated using formulas (3) and (4), respectively: where , , is the optimal step size of the individual particles, is the optimal step size of the particles, and are random numbers between 0 and 1, is the position of the particle , and is the speed of the particle .
Setting the maximum particle velocity interval as , when the particle velocity is higher than , the particle’s velocity is back to at this time.(4)Comparing the current iteration number with the maximum number of iterations : if , stop iterating; if , return to step (1).
3.2. The Improvement of the PSOBP Neural Network Based on JMP and Markov Chain
Although the prediction accuracy of the PSOBP neural network is better than that of some previous neural networks, it still cannot reach the actual ideal state of engineering, mainly because there are some deviations in the complexity of input data and combined prediction. Therefore, it is necessary to use the JMP data processing software [16] to process the input data and correct the prediction errors by the Markov chain.
Considering that some of the above data have low correlation between input data and output data, JMP data analysis software is used for correlation analysis and Pearson’s correlation significance test is conducted. [0, 0.3] is considered as weak correlation, [0.3, 0.5] as medium correlation, [0.5, 0.7] as strong correlation, and [0.7, 1.0] as strong correlation. The results of correlation analysis are shown in Table 1.

Combining with the analysis results of the foregoing and Table 1, we use the three input parameters R2, W, and C4 and good seven variables T, T1, H1, R, R1, 1, and CC2 as the input variables of the PSO and BP network prediction models and the parameter C as the output variable; combined with formula (1) for large public buildings, the cooling load prediction model structure is 7151, i.e., 7 input layer nodes, 15 hidden layer nodes, and 1 output layer node.
The Markov chain is a random process with no aftereffects [17]. It can derive the probability distribution of the next moment according to the time condition of a certain condition and obtain the state of the next moment regardless of other moments. The Markov chain is suitable for the correction of prediction problems with fluctuation [18].
The steps of Markov error correction for the predicted value of the above neural network model are as follows:(1)The error value matrix is calculated according to the predicted value matrix and the measured value matrix obtained by the PSOBP neural network, and the mean value and the standard deviation of obtained by formulas (5) and (6) are used to obtain the error state interval:(2)According to formula (7) and statistical method, the transition probability matrix is obtained according to the error state interval obtained in step (1): where is the probability transition matrix, is the error of state in which is the total number of occurrences, and is the error of state with being the number of transitions from state to state .(3)The state vector at step is obtained by using formula (8) to judge the state interval at step : where is the probability distribution under the initial condition and is the onestep probability transfer matrix, as shown in equation (6):
4. Simulation Analysis of Cooling Load Forecasting in Shopping Mall Buildings
4.1. Research and Analysis on the Mechanism of Cooling Load Change
According to the periodicity and similarity of the changing trend of cooling load, it is speculated that the changing law of cooling load in shopping malls should have similar characteristics. Therefore, taking a commercial center in Xi’an city as an example, the data of cooling load in a week and on a weekend are analyzed, respectively. According to the survey, the cooling load information from Monday to Friday in a week is similar, but it is not quite the same as that on weekends. To take the commercial center as the sample object, the cooling load variation chart, as shown in Figure 2, can be obtained to analyze the cooling load variation during the week and on weekends, and then use the line chart to represent it. It can be shown from Figure 3 that there is a difference in the variation of cooling load within a week and on a weekend. Therefore, in the subsequent experimental simulation, experiments are conducted on analyzing a weekday and a weekend, respectively.
4.2. Case Validation Analysis
According to the previous analysis, this experiment uses the actual data of a commercial center in Xi’an. The commercial center has ten floors, with a height of 40.6 meters and an area of 250,000 square meters. The building airconditioning area occupies about 187,600 square meters. The mall uses ice storage airconditioning for cooling in summer.
In this experiment, seven variables of outdoor dry bulb temperatures T and t − 1 at the current moment and the previous moment, solar radiation quantities R and r − 1 at the current moment and the previous moment, outdoor humidity H at the current moment, and cooling loads c − 1 and c − 2 at the previous moment and the first two moments are adopted as input variables of the PSOBP neural network, and cooling load C at the current moment is taken as the output variable of the neural network. Because of the characteristics of the commercial center, the experiment uses the input variables from 8:00 am to 10:00 pm in the first 50 days of the mall in June and July as training data and the data from July 26^{th} to July 29^{th} as the verification analysis data.
The test data is used for simulation to verify the prediction ability of the network. Figure 3 shows the predicted values of each model, where model 1 is the prediction of the weekday data without any distinctions, model 2 is the prediction of the weekly data but does not eliminate the input data with low correlation, and model 3 distinguishes and predicts the data on weekends and weekdays after eliminating the low correlation of input data. Figure 4 shows the prediction errors among the three models. According to Figures 3 and 4, eliminating the low correlation of input data could cause a better prediction.
It can be found from Figure 4 that the prediction error of model 3 is better than that of models 1 and 2, but it still has a peak error. Therefore, the Markov chain should be used to modify it and to improve the prediction accuracy of the model.
4.3. The Correction of the Markov Chain
The modeling process of the PSOBP neural network based on the Markov chain modification firstly divides the interval of the range according to the relative error obtained above. The method of interval division is the meanvariance method [19, 20]. Furthermore, the error is classified according to the divided interval, and the probability transition matrix is constructed according to the classification results to determine the state of the initial vector and the state of the predicted time period, and the modified value of the predicted time period is calculated according to the state of the predicted time period.
According to the relative error of the above analysis, the weekly error mean value and mean variance are obtained. The mean value of error over the weekend is , and the mean variance is going to be . After calculation, the state interval of Markov chain within a week and on the weekend is shown in Table 2.

According to the divided state interval, the error values of the PSOBP model within the prediction week and on the weekend are divided in turn, and the results are shown in Table 3.

According to the error state of Table 2, it is determined that the Markov chain onestep transition probability matrices and in the week and on the weekend are as follows:
According to the error state of Table 2, the state vector of each time period is determined, and the predicted value of the PSOBP neural network prediction model corrected by the Markov chain is solved according to the state transition vector and the probability transition matrix.
The corrected predicted value (model 4) is predicted by the aforementioned undivided weekly data (model 1), and the weekly data is divided, but the input data is not processed to perform prediction (model 2), which means the weekly data is distinguished and the weekend data is processed by comparing the input data to predict the predicted results (model 3). The results are shown in Figure 5.
(a)
(b)
(c)
(d)
RMSPE (root mean square percentage error) was used to evaluate the errors of the four prediction models, and the processing results are shown in Table 4:

The experimental results show that when the improved prediction of the PSOBP neural network combined with the Markov chain is performed, the maximum number of iterations is set to 2000 times, and the prediction model achieves a prediction accuracy of only 409 steps, which improves the validity of the prediction results. From the results of Figure 5 and Table 4, it can be seen that using the Markov chain to modify the PSOBP neural network figures out a higher prediction accuracy, and it can meet the requirements of practical engineering applications.
5. Conclusion
Aiming at the problem of cooling load prediction from ice storage airconditioning in malls, the paper proposes a prediction model to improve the PSOBP neural network and the method which uses the Markov chain to correct the error. The JMP data analysis software is used to analyze the correlation between input and output data. The PSOBP neural network is used to predict the cooling load of ice storage airconditioning, and the Markov chain is used to correct the prediction results. Applying the predictive model to the actual mall leads to the following conclusions:(1)The output cooling load of the ice storage airconditioning system in the shopping mall at time T is highly correlated with the 7 inputs as follows: the outdoor temperature at time T, the outdoor temperature at time t1, the solar radiation at time T, the solar radiation at time t1, the relative humidity at time T, the cooling load at time t1, and the cooling load at time t2.(2)The maximum relative error between the improved model and the unimproved model was reduced from 33.799% to 10.967%, and the root mean square error of the prediction model was reduced from 26.14% to 3.06%.(3)The prediction results show that the method we use can improve the prediction accuracy. Besides, it is more suitable for practical engineering application.
Nomenclature
:  Error matrix 
:  Measured value matrix 
:  Predictive matrix 
:  Probability transfer matrix 
:  State interval 
:  Total number of errors 
:  Error number 
A:  Number of hidden layer neurons 
B:  Number of neurons in the input layer 
C:  Cooling load 
l:  Particle step 
R:  Solar radiation 
T:  Outdoor temperature 
V:  Particle velocity 
W:  Outdoor wind speed 
H:  Outdoor humidity 
f:  Particle fitness 
X:  Particle position 
SVM:  Supportvector machine 
PSOBP:  Particle swarm optimizationback propagation 
BP:  Back propagation 
PSO:  Particle swarm optimization 
RMSPE:  Root mean square percentage error. 
Data Availability
The data used in this paper were from specific research projects, data collection, and a business center in Xi’an, and the data are safe and reliable. However, due to the need of later scientific research and the confidentiality of the market, the data cannot be disclosed for the time being.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Key R&D Projects in Shaanxi Province (2017zdcxlsf0302 and 2017ZDLSF165) and Shaanxi Provincial Education Department Industrialization Cultivation Project (17JF 016).
References
 L. Jingwan, W. Dong, and W.M. Zhu, “Application of neural network fuzzy predictive optimal control in VAV system,” Journal of System Simulation, vol. 22, no. 12, pp. 2909–2914, 2010. View at: Google Scholar
 J. Yu, Y. Wang, C. Xu, A. Zhao, and L. Yan, “Study on operation optimization of ice storage air conditioning system based on particle swarm optimization,” Journal of Xi‘an University of Architecture and Technology, vol. 50, no. 01, pp. 148–154, 2018. View at: Google Scholar
 Y. Wei, X. Zhang, Y. Shi et al., “A review of datadriven approaches for prediction and classification of building energy consumption,” Renewable and Sustainable Energy Reviews, vol. 82, pp. 1027–1047, 2018. View at: Publisher Site  Google Scholar
 V. Vapnik, S. E. Golowich, and A. Smola, “Support vector method for function approximation, regression estimation, and signal processing,” Advances in Neural Information Processing Systems, vol. 9, pp. 281–287, 1996. View at: Google Scholar
 H. X. Zhao and F. Magoulès, “Parallel support vector machines applied to the prediction of multiple buildings energy consumption,” Journal of Algorithms & Computational Technology, vol. 4, no. 2, pp. 231–249, 2010. View at: Publisher Site  Google Scholar
 Q. Li, Q. Meng, J. Cai, H. Yoshino, and A. Mochida, “Predicting hourly cooling load in the building: a comparison of support vector machine and different artificial neural networks,” Energy Conversion and Management, vol. 50, no. 1, pp. 90–96, 2009. View at: Publisher Site  Google Scholar
 L. G. Swan and V. I. Ugursal, “Modeling of enduse energy consumption in the residential sector: a review of modeling techniques,” Renewable and Sustainable Energy Reviews, vol. 13, no. 8, pp. 1819–1835, 2009. View at: Publisher Site  Google Scholar
 Z. Yu, F. Haghighat, B. C. M. Fung, and H. Yoshino, “A decision tree method for building energy demand modeling,” Energy and Buildings, vol. 42, no. 10, pp. 1637–1646, 2010. View at: Publisher Site  Google Scholar
 G. K. F. Tso and K. K. W. Yau, “Predicting electricity energy consumption: a comparison of regression analysis, decision tree and neural networks,” Energy, vol. 32, no. 9, pp. 1761–1768, 2007. View at: Publisher Site  Google Scholar
 H.G. Beyer, “Evolutionary algorithms in noisy environments: theoretical issues and guidelines for practice,” Computer Methods in Applied Mechanics and Engineering, vol. 186, no. 2–4, pp. 239–267, 2000. View at: Publisher Site  Google Scholar
 A. Azadeh, S. F. Ghaderi, S. Tarverdian, and M. Saberi, “Integration of artificial neural networks and genetic algorithm to predict electrical energy consumption,” Applied Mathematics and Computation, vol. 186, no. 2, pp. 1731–1741, 2007. View at: Publisher Site  Google Scholar
 K. Li and H. Su, “Forecasting building energy consumption with hybrid genetic algorithmhierarchical adaptive networkbased fuzzy inference system,” Energy and Buildings, vol. 42, no. 11, pp. 2070–2076, 2010. View at: Publisher Site  Google Scholar
 A. Setiawan, I. Koprinska, and V. G. Agelidis, “Very shortterm electricity load demand forecasting using support vector regression,” in Proceedings of the 2009 International Joint Conference on Neural Networks (IJCNN), pp. 2888–2894, Atlanta, GA, USA, June 2009. View at: Publisher Site  Google Scholar
 J. Huang and L. He, “Application of improved PSO—BP neural network in customer churn warning,” Procedia Computer Science, vol. 131, pp. 1238–1246, 2018. View at: Publisher Site  Google Scholar
 X. Zhou, Z. Fan, L. Guoqiang, and Y. Junwei, “Load prediction of air conditioning in shopping malls based on multiple nonlinear regression,” HVAC, vol. 48, no. 03, pp. 120–195, 2018. View at: Google Scholar
 X. Zhao, W. Zhao, and H. Chen, “Application of new method for pH determination and evaluation,” Journal of Pharmaceutical Analysis, vol. 38, no. 03, pp. 545–550, 2018. View at: Google Scholar
 Y. Liu, B. Liu, and A. Shi, “Modeling and control of network control system with interval random delay,” Journal of System Simulation, vol. 30, no. 02, pp. 654–663, 2018. View at: Google Scholar
 X. Lu, Y. Shen, J. Chen, and Z. Ji, “Study on ultrashort term wind power interval prediction considering wind power generation randomness,” Acta Solar Energy Sinica, vol. 38, no. 05, pp. 1307–1315, 2017. View at: Google Scholar
 Y. Huang, J. Peng, C. Li, L. Ding, and G. Sun, “Application of Markov theory in medium and longterm load forecasting,” Journal of Power System and Automation, vol. 23, no. 05, pp. 131–136, 2011. View at: Google Scholar
 Y. Jing, X. Zhang, and Y. Luo, “Water demand forecast of Qinghai province based on modified combination model,” Journal of Natural Resources, vol. 27, no. 06, pp. 1013–1021, 2012. View at: Google Scholar
Copyright
Copyright © 2019 JunQi Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.