<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">CMES</journal-id>
<journal-id journal-id-type="nlm-ta">CMES</journal-id>
<journal-id journal-id-type="publisher-id">CMES</journal-id>
<journal-title-group>
<journal-title>Computer Modeling in Engineering &#x0026; Sciences</journal-title>
</journal-title-group>
<issn pub-type="epub">1526-1506</issn>
<issn pub-type="ppub">1526-1492</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">25396</article-id>
<article-id pub-id-type="doi">10.32604/cmes.2023.025396</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Research on Short-Term Load Forecasting of Distribution Stations Based on the Clustering Improvement Fuzzy Time Series Algorithm</article-title>
<alt-title alt-title-type="left-running-head">Research on Short-Term Load Forecasting of Distribution Stations Based on the Clustering Improvement Fuzzy Time Series Algorithm</alt-title>
<alt-title alt-title-type="right-running-head">Research on Short-Term Load Forecasting of Distribution Stations Based on the Clustering Improvement Fuzzy Time Series Algorithm</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Gu</surname><given-names>Jipeng</given-names>
</name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Zhang</surname><given-names>Weijie</given-names>
</name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-3" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Zhang</surname><given-names>Youbing</given-names>
</name><xref ref-type="aff" rid="aff-1">1</xref><email>youbingzhang@zjut.edu.cn</email></contrib>
<contrib id="author-4" contrib-type="author">
<name name-style="western"><surname>Wang</surname><given-names>Binjie</given-names>
</name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-5" contrib-type="author">
<name name-style="western"><surname>Lou</surname><given-names>Wei</given-names>
</name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-6" contrib-type="author">
<name name-style="western"><surname>Ye</surname><given-names>Mingkang</given-names>
</name><xref ref-type="aff" rid="aff-3">3</xref></contrib>
<contrib id="author-7" contrib-type="author">
<name name-style="western"><surname>Wang</surname><given-names>Linhai</given-names>
</name><xref ref-type="aff" rid="aff-3">3</xref></contrib>
<contrib id="author-8" contrib-type="author">
<name name-style="western"><surname>Liu</surname><given-names>Tao</given-names>
</name><xref ref-type="aff" rid="aff-4">4</xref></contrib>
<aff id="aff-1"><label>1</label><institution>College of Information Engineering, Zhejiang University of Technology</institution>, <addr-line>Hangzhou, 310023</addr-line>, <country>China</country></aff>
<aff id="aff-2"><label>2</label><institution>Electric Power Research Institute of State Grid Anhui Electric Power Co., Ltd.,</institution> <addr-line>Hefei, 230601</addr-line>, <country>China</country></aff>
<aff id="aff-3"><label>3</label><institution>State Grid Zhejiang Electric Power Co., Ltd., Wencheng County Power Supply Company</institution>, <addr-line>Wenzhou, 325000</addr-line>, <country>China</country></aff>
<aff id="aff-4"><label>4</label><institution>Zhejiang Tusheng Power Transmission and Transformation Engineering Co., Ltd., Tusheng Technology Branch</institution>, <addr-line>Wenzhou, 325000</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Youbing Zhang. Email: <email>youbingzhang@zjut.edu.cn</email></corresp>
</author-notes>
<pub-date date-type="collection" publication-format="electronic"><year>2023</year></pub-date>
<pub-date date-type="pub" publication-format="electronic">
<day>2</day><month>3</month><year>2023</year>
</pub-date>
<volume>136</volume>
<issue>3</issue>
<fpage>2221</fpage>
<lpage>2236</lpage>
<history>
<date date-type="received">
<day>09</day><month>7</month><year>2022</year>
</date>
<date date-type="accepted">
<day>03</day><month>11</month><year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2023 Gu et al.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Gu et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_CMES_25396.pdf"></self-uri>
<abstract>
<p>An improved fuzzy time series algorithm based on clustering is designed in this paper. The algorithm is successfully applied to short-term load forecasting in the distribution stations. Firstly, the K-means clustering method is used to cluster the data, and the midpoint of two adjacent clustering centers is taken as the dividing point of domain division. On this basis, the data is fuzzed to form a fuzzy time series. Secondly, a high-order fuzzy relation with multiple antecedents is established according to the main measurement indexes of power load, which is used to predict the short-term trend change of load in the distribution stations. Matlab/Simulink simulation results show that the load forecasting errors of the typical fuzzy time series on the time scale of one day and one week are [&#x2212;50, 20] and [&#x2212;50, 30], while the load forecasting errors of the improved fuzzy time series on the time scale of one day and one week are [&#x2212;20, 15] and [&#x2212;20, 25]. It shows that the fuzzy time series algorithm improved by clustering improves the prediction accuracy and can effectively predict the short-term load trend of distribution stations.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Short-term load forecasting</kwd>
<kwd>fuzzy time series</kwd>
<kwd>K-means clustering</kwd>
<kwd>distribution stations</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label><title>Introduction</title>
<p>Power demand forecasting plays an important role in modern power system research. Accurate prediction of power load in different periods in the future will improve the management level of the power system [<xref ref-type="bibr" rid="ref-1">1</xref>,<xref ref-type="bibr" rid="ref-2">2</xref>]. According to the forecasting time, load forecasting can be divided into three main forms: short term, long term and super long term. Short-term load forecasting is the prediction of power load in the next few minutes to a week. Accurate short-term load forecasting will help to formulate a reasonable power production plan. It also can avoid excessive waste of power resources and improve the economic benefits of the power system [<xref ref-type="bibr" rid="ref-3">3</xref>&#x2013;<xref ref-type="bibr" rid="ref-5">5</xref>]. Short-term load forecasting methods mainly fall into two categories: statistical methods and machine learning methods. Statistical methods mainly include regression analysis, time series, Markov chain, etc. Machine learning methods mainly include support vector machines (SVM), artificial neural networks, etc. [<xref ref-type="bibr" rid="ref-6">6</xref>&#x2013;<xref ref-type="bibr" rid="ref-8">8</xref>].</p>
<p>In the study of statistical methods, short-term load forecasting was performed by using an improved Gaussian process regression model with multi-core covariance, and the interval prediction results at a certain confidence level were obtained [<xref ref-type="bibr" rid="ref-9">9</xref>]. The gray time series modified by Markov was used to predict the power load trend, and the prediction accuracy was higher than that of single algorithm models such as time series and Markov chain [<xref ref-type="bibr" rid="ref-10">10</xref>]. A blind Kalman filter algorithm was proposed for short-term load forecasting, which has great advantages in load profile analysis and peak load forecasting by predicting unknown matrix and state alternation estimation [<xref ref-type="bibr" rid="ref-11">11</xref>]. The mixed random forest algorithm and mean generating functions were used to form a hybrid short-term load forecasting model, which has better prediction accuracy for peak and valley conditions with the large change of load data [<xref ref-type="bibr" rid="ref-12">12</xref>]. An adaptive hybrid fractal model was proposed for power system short-term load forecasting, composed of composite linear fractal difference function, iterative learning and optimization algorithm, and has higher accuracy than commonly used time series methods [<xref ref-type="bibr" rid="ref-13">13</xref>]. A load forecasting method based on the combination of multiple phase space reconstruction (MPSR) and support vector regression (SVR) was proposed, which takes into account the coupling relationship between multiple energy loads and has high prediction efficiency and accuracy [<xref ref-type="bibr" rid="ref-14">14</xref>]. A seasonal autoregressive integrated moving average and variance-covariance prediction method was proposed, which considered the interaction of multiple performance indicators and had high prediction accuracy&#x00A0;[<xref ref-type="bibr" rid="ref-15">15</xref>].</p>
<p>In the study of machine learning methods, the joint learning method based on recurrent neural networks was used to predict short-term load changes, which has a good prediction effect [<xref ref-type="bibr" rid="ref-16">16</xref>]. The prediction interval of the neural network was used to construct a lower bound estimation method, which can quantify the potential uncertainty factors related to prediction and predict the load change trend with high accuracy in a short time [<xref ref-type="bibr" rid="ref-17">17</xref>]. The two-stage attention mechanism based on the long-term and short-term memory (LSTM) neural network was introduced for the probability prediction of short-term regional load, and the prediction model has higher accuracy and generalization ability [<xref ref-type="bibr" rid="ref-18">18</xref>]. A time convolutional neural network that integrated channel and time attention mechanism was proposed for short-term load forecasting of the power system, which effectively expressed the nonlinear relationship between meteorological factors and power load [<xref ref-type="bibr" rid="ref-19">19</xref>]. In [<xref ref-type="bibr" rid="ref-20">20</xref>], a holographic integrated forecasting method for short-term power load was proposed, which integrates multi-category and multi-state load information into four levels (data set, sampling space, prediction model and decision), which can comprehensively integrate information throughout the whole life cycle of the forecasting process and greatly improve the effect of short-term load forecasting. In [<xref ref-type="bibr" rid="ref-21">21</xref>], Box-Cox conversion processing and parameter fitting of Copula model were carried out for loaded load data, and a data-driven deep confidence network was proposed to predict the hourly load of the power system. In [<xref ref-type="bibr" rid="ref-22">22</xref>], nonlinear exogenous recurrent neural network (NARX), Elman neural network and autoregressive moving average (ARMA) were used for short-term load forecasting, and it was found that the average absolute percentage errors of NARX, Elman and ARMA were 5.53%, 3.42% and 10.28%, respectively. In [<xref ref-type="bibr" rid="ref-23">23</xref>], a personalized federated learning method was proposed, with high prediction accuracy in individual consumer load forecasting. In addition to the above literature, An et al. [<xref ref-type="bibr" rid="ref-24">24</xref>&#x2013;<xref ref-type="bibr" rid="ref-27">27</xref>] also adopted different methods to study load forecasting of the power system.</p>
<p>The statistical method is simple in principle and fast in the calculation, but it has limited ability to deal with nonlinear variables. While the machine learning method can approach the nonlinear function relationship with arbitrary precision in principle, it is difficult to mine the timing characteristics between data [<xref ref-type="bibr" rid="ref-28">28</xref>,<xref ref-type="bibr" rid="ref-29">29</xref>]. How to integrate intelligent algorithms into typical time series prediction methods to improve the processing ability of its nonlinear function is worthy of in-depth research. Time series analysis includes two parts: time series modeling and time series forecasting. Modeling is the rational cognition of the internal development law of things, while forecasting is the specific performance of future development trend based on the model. The classical time series analysis method can deal with most realistic problems, but it cannot deal with the imprecise, incomplete or fuzzy realistic problems. For example, the change of power load is often expressed as &#x201C;sudden increase, sudden decrease, steady&#x201D; and other language variables. Although these phenomena can be described with accurate numerical values, it is very difficult to obtain historical data, similar to the above vague or incomplete data are more. In this case, fuzzy time series expressed by linguistic variables can be used for prediction.</p>
<p>At present, the theoretical research of fuzzy time series mainly focuses on the reasonable division of fuzzy interval. Fuzzy interval has a great influence on the calculation process and prediction accuracy of the model, which is the basis of establishing fuzzy time series prediction model. In [<xref ref-type="bibr" rid="ref-30">30</xref>], the concept of fuzzy time series was put forward, which was a pioneer in the theory and application of fuzzy time series. In [<xref ref-type="bibr" rid="ref-31">31</xref>,<xref ref-type="bibr" rid="ref-32">32</xref>], the maximum and minimum values of the sample data were rounded as the domain of the model, and then the fixed interval length was selected to divide the domain equally. The membership function setting of this method was simple and the calculation speed was fast, but the fuzzy set corresponding to the interval was not accurate, which led to the low prediction accuracy. In [<xref ref-type="bibr" rid="ref-33">33</xref>,<xref ref-type="bibr" rid="ref-34">34</xref>], the distribution characteristics of sample data were used to divide the domain. The interval with dense sample data was reduced, and the interval with sparse sample data was expanded. The interval division based on statistical characteristics was more reasonable, and the prediction accuracy of the model was improved to a certain extent. In [<xref ref-type="bibr" rid="ref-35">35</xref>&#x2013;<xref ref-type="bibr" rid="ref-37">37</xref>], neural networks and optimization algorithms were used to divide the domain, aiming to find the optimal number and length of intervals. The prediction accuracy of this method was greatly improved, but the divided interval was difficult to be interpreted in natural language, which weakened the advantage of fuzzy theory in the application of time series prediction. Using the clustering algorithm to cluster the sample data, and taking the clustering results as the basis for domain division, each interval represents a clearer practical significance, and the prediction effect is also very accurate, so this kind of algorithm is very meaningful.</p>
<p>Power load forecasting is the basis of power system operation management and real-time regulation. It is an important link in the reform of the electric power system and the transformation of energy structure to constantly improve the power load forecasting technology and seek a more accurate load forecasting model. It is helpful for decision makers to make reasonable power grid dispatching plan and maintain the safe and economic operation of the power grid. Combined with the scenario of short-term load forecasting in the distribution station, this paper focuses on the prediction performance of fuzzy time series and its improved algorithm. The innovation points of this paper are as follows:
<list list-type="order">
<list-item>
<p>The typical fuzzy time series can predict the time series data containing fuzzy information or incomplete information, and the improved fuzzy time series still has this feature.</p></list-item>
<list-item>
<p>The typical fuzzy time series adopts the principle of equal division of the domain. While the improved fuzzy time series divides the domain according to the probability distribution characteristics of the data, and the division of the domain is more reasonable.</p></list-item>
<list-item>
<p>Fuzzy time series represented by high-order fuzzy relations can better fit the nonlinear relationship between data. This paper presents a concrete design method for the third order fuzzy time series.</p></list-item>
<list-item>
<p>The prediction algorithm proposed in this paper can be applied to other fields, such as traffic flow and weather, etc. The prediction accuracy of the algorithm is high, and the algorithm is universal.</p></list-item>
</list></p>
</sec>
<sec id="s2">
<label>2</label><title>K-Means Clustering and Fuzzy Time Series</title>
<sec id="s2_1">
<label>2.1</label><title>K-Means Clustering Method</title>
<p>Clustering is the process of classifying and organizing samples with similar characteristics in a data set. As one of the most famous partition clustering algorithms, K-means clustering is widely used due to its simplicity and efficiency. It is a clustering analysis algorithm based on iterative solution process [<xref ref-type="bibr" rid="ref-38">38</xref>,<xref ref-type="bibr" rid="ref-39">39</xref>]. The algorithm steps are as follows:
<list list-type="simple">
<list-item><label>1)</label><p><italic>k</italic> samples are randomly selected from the data set containing <italic>n</italic> samples as the initial clustering center <italic>C</italic><sub><italic>i</italic></sub> (<italic>i</italic> &#x003D; 1, 2, &#x2026;, <italic>k</italic>).</p></list-item>
<list-item><label>2)</label><p>According to the initial value of each cluster center, the absolute distance <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula> from each sample <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> to the cluster center <italic>C</italic><sub><italic>i</italic></sub> was calculated. Then the samples are divided into class clusters <italic>L</italic><sub><italic>i</italic></sub> (<italic>i</italic> &#x003D; 1, 2, &#x2026;, <italic>N</italic><sub><italic>i</italic></sub>) according to the shortest distance, <italic>N</italic><sub><italic>i</italic></sub> is the number of samples contained in each cluster.</p></list-item>
<list-item><label>3)</label><p>Calculating the mean of each cluster <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2208;</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munder><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, which is used as the clustering center after the algorithm is updated.</p></list-item>
<list-item><label>4)</label><p>Repeating Steps 2 and 3 until all samples cannot be redistributed.</p></list-item>
</list></p>
</sec>
<sec id="s2_2">
<label>2.2</label><title>Fuzzy Time Series</title>
<p>Fuzzy time series analysis refers to the use of fuzzy mathematics theory to study the deep nonlinear relationship contained in the time series containing fuzzy or uncertain information. The relevant definitions involved are as follows [<xref ref-type="bibr" rid="ref-40">40</xref>&#x2013;<xref ref-type="bibr" rid="ref-42">42</xref>]:</p>
<p>Definition 1: Let U be the domain of time series, and divide the domain into <italic>n</italic> ordered subintervals, i.e., U &#x003D; {<italic>u</italic><sub>1</sub>, <italic>u</italic><sub>2</sub>, &#x2026;, <italic>u</italic><sub>n</sub>}, then the fuzzy set <italic>A</italic> defined on domain U can be expressed as:</p>
<p><disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mi>A</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /></mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /></mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /></mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>where <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mo>&#x22C5;</mml:mo><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is membership function, <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the membership of <italic>u</italic><sub><italic>i</italic></sub> to <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:mi>A</mml:mi></mml:math></inline-formula>, <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x03F5;</mml:mi></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>]</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi></mml:math></inline-formula>.</p>
<p>Definition 2: Let <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:mi>X</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> be a subset of the real number set <italic>R</italic>, and <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> be a fuzzy set defined on the domain <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><mml:mi>X</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>. If <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>is a set composed of <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>, then <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is a fuzzy time series defined on <inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><mml:mi>X</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>.</p>
<p>Definition 3: If <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, the fuzzy logic relationship (FLR) can be expressed as</p>
<p><inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">&#x2192;</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the left components of FLR, <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the right components of FLR, and <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><mml:mi>R</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is the fuzzy logic relationship from <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> to <inline-formula id="ieqn-21"><mml:math id="mml-ieqn-21"><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>, then the first-order model can be expressed as:</p>
<p><disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>F</mml:mi><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>&#x2218;</mml:mo></mml:mrow></mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-22"><mml:math id="mml-ieqn-22"><mml:msup><mml:mi></mml:mi><mml:mrow><mml:mo>&#x2218;</mml:mo></mml:mrow></mml:msup></mml:math></inline-formula> represents compositional operation. Similarly, the <italic>n</italic>-order model of <inline-formula id="ieqn-23"><mml:math id="mml-ieqn-23"><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> can be obtained, that is, the relationship between <inline-formula id="ieqn-24"><mml:math id="mml-ieqn-24"><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-25"><mml:math id="mml-ieqn-25"><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>n</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, which is expressed as follows:</p>
<p><disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>F</mml:mi><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>&#x2218;</mml:mo></mml:mrow></mml:msup><mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Definition 4: If the left components of multiple fuzzy logic relationships are the same, but the right components are different, these fuzzy logic relationships can be combined into a relationship, that is, a fuzzy logic relationship group.</p>
<p><disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">&#x2192;</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo></mml:math></disp-formula></p>
</sec>
</sec>
<sec id="s3">
<label>3</label><title>Fuzzy Time Series Prediction Model</title>
<sec id="s3_1">
<label>3.1</label><title>Typical Fuzzy Time Series Prediction (FTS)</title>
<p>Typical fuzzy time series prediction includes the following steps:</p>
<p>(1) Based on historical data, the domain is determined. Let <inline-formula id="ieqn-26"><mml:math id="mml-ieqn-26"><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> be a time series, then the domain U can be defined as:</p>
<p><disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:mrow><mml:mtext>U</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-27"><mml:math id="mml-ieqn-27"><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-28"><mml:math id="mml-ieqn-28"><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> are the minimum and maximum values in <inline-formula id="ieqn-29"><mml:math id="mml-ieqn-29"><mml:mi>X</mml:mi></mml:math></inline-formula>, and <inline-formula id="ieqn-30"><mml:math id="mml-ieqn-30"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-31"><mml:math id="mml-ieqn-31"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> are two customized real numbers.</p>
<p>(2) Dividing the domain equally to form multiple numerical subintervals. It is worth noting that too large or too small subinterval range will have a certain impact on the final prediction results, and its range must be selected reasonably.</p>
<p>(3) Fuzzy sets are defined to fuzzify the data. Suppose that there are <italic>m</italic> subintervals formed after the domain is divided equally, then the corresponding fuzzy set is defined as <inline-formula id="ieqn-32"><mml:math id="mml-ieqn-32"><mml:mi>A</mml:mi><mml:mo>=</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula>, and <inline-formula id="ieqn-33"><mml:math id="mml-ieqn-33"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the language variable of <inline-formula id="ieqn-34"><mml:math id="mml-ieqn-34"><mml:mi>A</mml:mi></mml:math></inline-formula>. Fuzzify the historical data according to the triangular membership function, and the process is as follows:</p>
<p><disp-formula id="ueqn-6"><mml:math id="mml-ueqn-6" display="block"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>11</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
<p><disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>21</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>22</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
<p><inline-formula id="ieqn-35"><mml:math id="mml-ieqn-35"><mml:mo>&#x22EF;</mml:mo></mml:math></inline-formula></p>
<p><disp-formula id="ueqn-7"><mml:math id="mml-ueqn-7" display="block"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
<p>(4) Establish fuzzy logic relation and determine fuzzy relation matrix. If <inline-formula id="ieqn-36"><mml:math id="mml-ieqn-36"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the fuzzy membership of data <italic>i</italic> and <inline-formula id="ieqn-37"><mml:math id="mml-ieqn-37"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the fuzzy membership of data <italic>j</italic>, then the fuzzy logic relationship can be expressed as <inline-formula id="ieqn-38"><mml:math id="mml-ieqn-38"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">&#x2192;</mml:mo><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. Thus, the fuzzy logic relation <inline-formula id="ieqn-39"><mml:math id="mml-ieqn-39"><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> between the values of each fuzzy time series can be written, and the fuzzy relation matrix can be obtained by combining the fuzzy logic relations <inline-formula id="ieqn-40"><mml:math id="mml-ieqn-40"><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x222A;</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>&#x22EF;</mml:mo></mml:math></inline-formula>.</p>
<p>(5) Defuzzification and prediction. According to the established fuzzy relation, the time series data can be predicted. During the prediction, the fuzzy quantity needs to be de fuzzified to obtain an accurate numerical quantity. The defuzzification method adopts the weighted average method, as shown in <xref ref-type="disp-formula" rid="eqn-7">Eq. (7)</xref>.</p>
<p><disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:mrow><mml:mtext>Y</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mo>&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>p</mml:mi></mml:math></disp-formula>where <inline-formula id="ieqn-41"><mml:math id="mml-ieqn-41"><mml:mi>p</mml:mi></mml:math></inline-formula> is the number of fuzzy relations corresponding to the fuzzy set <inline-formula id="ieqn-42"><mml:math id="mml-ieqn-42"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> obtained by fuzzifying the data <inline-formula id="ieqn-43"><mml:math id="mml-ieqn-43"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, <inline-formula id="ieqn-44"><mml:math id="mml-ieqn-44"><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the weight, and <inline-formula id="ieqn-45"><mml:math id="mml-ieqn-45"><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the intermediate value of the interval <inline-formula id="ieqn-46"><mml:math id="mml-ieqn-46"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> corresponding to the fuzzy set <inline-formula id="ieqn-47"><mml:math id="mml-ieqn-47"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>.</p>
</sec>
<sec id="s3_2">
<label>3.2</label><title>Improved Fuzzy Time Series Prediction (IFTS)</title>
<p>In the typical fuzzy time series prediction method, the equal division principle is used to divide the domain. However, from the probability distribution characteristics of the data, it can be seen that the data distribution is generally uneven. At this time, the equal division principle is unreasonable to divide the domain. Therefore, this paper considers the cluster centers as the basis for domain division, and the resulting subinterval will be more reasonable and have data similarity. The principles of domain division are as follows:</p>
<p><disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mtext>&#xA0;</mml:mtext><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mtext>&#xA0;</mml:mtext><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mtable columnalign="center center" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-48"><mml:math id="mml-ieqn-48"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is each subinterval, <inline-formula id="ieqn-49"><mml:math id="mml-ieqn-49"><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is the average value of cluster centers, <inline-formula id="ieqn-50"><mml:math id="mml-ieqn-50"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is the minimum difference of cluster centers, and the calculation formulas of <inline-formula id="ieqn-51"><mml:math id="mml-ieqn-51"><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-52"><mml:math id="mml-ieqn-52"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> are as follows:</p>
<p><disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:msub><mml:mi>l</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /></mml:mrow><mml:mn>2</mml:mn></mml:math></disp-formula>where <inline-formula id="ieqn-53"><mml:math id="mml-ieqn-53"><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is cluster center.</p>
<p>In addition, the typical fuzzy time series prediction method often uses the first-order fuzzy relationship, but the higher-order fuzzy relationship can better reflect the nonlinear function relationship between data. Therefore, this paper adopts the third-order fuzzy relation, that is, the data at t &#x2212; 2, t&#x00A0;&#x2212;&#x00A0;1 and t time are used to predict the data at t &#x002B; 1 time. The prediction process is as follows:</p>
<p><disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle scriptlevel="0"><mml:mrow><mml:mo maxsize="1.2em" minsize="1.2em">|</mml:mo></mml:mrow></mml:mstyle><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mstyle scriptlevel="0"><mml:mrow><mml:mo maxsize="1.2em" minsize="1.2em">|</mml:mo></mml:mrow></mml:mstyle></mml:math></disp-formula></p>
<p><disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:math></disp-formula></p>
<p><disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:mi>Z</mml:mi><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:math></disp-formula></p>
<p><disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:mi>Y</mml:mi><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>P</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msubsup><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup><mml:mo>]</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>S</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-54"><mml:math id="mml-ieqn-54"><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, <inline-formula id="ieqn-55"><mml:math id="mml-ieqn-55"><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, <inline-formula id="ieqn-56"><mml:math id="mml-ieqn-56"><mml:mi>Z</mml:mi><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, <inline-formula id="ieqn-57"><mml:math id="mml-ieqn-57"><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-58"><mml:math id="mml-ieqn-58"><mml:mi>Y</mml:mi><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> are the set calculation variables, <inline-formula id="ieqn-59"><mml:math id="mml-ieqn-59"><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> represents the data at time <italic>i</italic>, <inline-formula id="ieqn-60"><mml:math id="mml-ieqn-60"><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> represents the predicted value at time <italic>i</italic> &#x002B; 1, <italic>P</italic>, <italic>S</italic> are temporary variables for data accumulation, with values as shown in <xref ref-type="table" rid="table-1">Table 1</xref>, <inline-formula id="ieqn-61"><mml:math id="mml-ieqn-61"><mml:mrow><mml:mo>[</mml:mo><mml:msubsup><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> is the interval <inline-formula id="ieqn-62"><mml:math id="mml-ieqn-62"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> with membership of 1 in the fuzzy set <inline-formula id="ieqn-63"><mml:math id="mml-ieqn-63"><mml:msub><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, and <inline-formula id="ieqn-64"><mml:math id="mml-ieqn-64"><mml:mi>M</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msubsup><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> is the median value of interval <inline-formula id="ieqn-65"><mml:math id="mml-ieqn-65"><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>.</p>
<table-wrap id="table-1">
<label>Table 1</label>
<caption><title>Values of P, S</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th rowspan="2">Variable</th>
<th align="center" colspan="2">In [<inline-formula id="ieqn-66"><mml:math id="mml-ieqn-66"><mml:msubsup><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>]</th>
<th align="center" colspan="2">Not in [<inline-formula id="ieqn-67"><mml:math id="mml-ieqn-67"><mml:msubsup><mml:mi>A</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>]</th>
</tr>
<tr>
<th><italic>P</italic></th>
<th><italic>S</italic></th>
<th><italic>P</italic></th>
<th><italic>S</italic></th>
</tr>
</thead>
<tbody>
<tr>
<td><inline-formula id="ieqn-68"><mml:math id="mml-ieqn-68"><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><inline-formula id="ieqn-69"><mml:math id="mml-ieqn-69"><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td>1</td>
<td>0</td>
<td>0</td>
</tr>
<tr>
<td><inline-formula id="ieqn-70"><mml:math id="mml-ieqn-70"><mml:mi>Z</mml:mi><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>P</italic> &#x002B; <inline-formula id="ieqn-71"><mml:math id="mml-ieqn-71"><mml:mi>Z</mml:mi><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>S</italic> &#x002B; 1</td>
<td><italic>P</italic></td>
<td><italic>S</italic></td>
</tr>
<tr>
<td><inline-formula id="ieqn-72"><mml:math id="mml-ieqn-72"><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>P</italic> &#x002B; <inline-formula id="ieqn-73"><mml:math id="mml-ieqn-73"><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>S</italic> &#x002B; 1</td>
<td><italic>P</italic></td>
<td><italic>S</italic></td>
</tr>
<tr>
<td><inline-formula id="ieqn-74"><mml:math id="mml-ieqn-74"><mml:mi>Y</mml:mi><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>P</italic> &#x002B; <inline-formula id="ieqn-75"><mml:math id="mml-ieqn-75"><mml:mi>Y</mml:mi><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td><italic>S</italic> &#x002B; 1</td>
<td><italic>P</italic></td>
<td><italic>S</italic></td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3">
<label>3.3</label><title>Evaluation of Prediction Accuracy</title>
<p>The prediction accuracy of time series is often tested by residual size. In this paper, four indexes including mean square error (MSE), root mean square error (RMSE), mean relative error (MRE) and efficiency coefficient NS are selected to evaluate the prediction accuracy of the model, which can be expressed by the following formulas:</p>
<p><disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:mrow><mml:mtext>MSE</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>y</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mi>n</mml:mi></mml:mfrac></mml:math></disp-formula></p>
<p><disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mrow><mml:mtext>RMSE</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:msqrt><mml:mn>1</mml:mn><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:msubsup><mml:mrow><mml:mo>&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>y</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:msqrt></mml:math></disp-formula></p>
<p><disp-formula id="eqn-21"><label>(21)</label><mml:math id="mml-eqn-21" display="block"><mml:mrow><mml:mtext>MRE</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>n</mml:mi></mml:mfrac><mml:msubsup><mml:mrow><mml:mo>&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mo>|</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>y</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mfrac><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mn>100</mml:mn><mml:mi mathvariant="normal">&#x0025;</mml:mi></mml:math></disp-formula></p>
<p><disp-formula id="eqn-22"><label>(22)</label><mml:math id="mml-eqn-22" display="block"><mml:mrow><mml:mtext>NS</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>y</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>y</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:math></disp-formula>where <inline-formula id="ieqn-76"><mml:math id="mml-ieqn-76"><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the measured value, <inline-formula id="ieqn-77"><mml:math id="mml-ieqn-77"><mml:msub><mml:mrow><mml:mover><mml:mi>y</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the predicted value, <inline-formula id="ieqn-78"><mml:math id="mml-ieqn-78"><mml:mover><mml:mi>y</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover></mml:math></inline-formula> is the average value of the measured value, and <italic>n</italic> is the numbers of data. When EMS &#x003D; 0, RMSE &#x003D; 0, MARE &#x003D; 0, NS &#x003D; 1, the model fitting effect is the best.</p>
<p>In addition, it is pointed out in [<xref ref-type="bibr" rid="ref-43">43</xref>] that the modified Diebold-Mariano test statistic (MDM) can be used to evaluate whether the loss function between different models is statistically significant. If the loss function error of each model is statistically different, the prediction ability of each model is obviously different.</p>
<p>Assuming that the sequence formed from the measured values is {<inline-formula id="ieqn-79"><mml:math id="mml-ieqn-79"><mml:msub><mml:mrow><mml:mtext>f</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:math></inline-formula>} and the sequence formed from the predicted values is {<inline-formula id="ieqn-80"><mml:math id="mml-ieqn-80"><mml:msub><mml:mrow><mml:mover><mml:mrow><mml:mtext>f</mml:mtext></mml:mrow><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:math></inline-formula>}, the prediction error of the model can be expressed as 
<xref ref-type="disp-formula" rid="eqn-23">Eq. (23)</xref>.</p>
<p><disp-formula id="eqn-23"><label>(23)</label><mml:math id="mml-eqn-23" display="block"><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>f</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>The loss function at time <italic>t</italic> can be expressed as <xref ref-type="disp-formula" rid="eqn-24">Eq. (24)</xref>.</p>
<p><disp-formula id="eqn-24"><label>(24)</label><mml:math id="mml-eqn-24" display="block"><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mo>=</mml:mo><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>f</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>L</italic> is the sign of the loss function.</p>
<p>The prediction results of different models are different, resulting in different prediction errors and loss function values. When the performance of the two prediction models is compared, the H<sub>0</sub> hypothesis theory can be used.</p>
<p><disp-formula id="eqn-25"><label>(25)</label><mml:math id="mml-eqn-25" display="block"><mml:msub><mml:mrow><mml:mtext>H</mml:mtext></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>:</mml:mo></mml:mrow><mml:mtext>&#xA0;</mml:mtext><mml:mi>E</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:math></disp-formula></p>
<p><disp-formula id="eqn-26"><label>(26)</label><mml:math id="mml-eqn-26" display="block"><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>e</italic><sub>1<italic>t</italic></sub> is the prediction error of model 1, <italic>e</italic><sub>2<italic>t</italic></sub> is the prediction error of model 2. It is worth noting that the two models satisfying the H<sub>0</sub> hypothesis will have the same predictive performance.</p>
<p>MDM is an important indicator to evaluate the predictive ability of a model and can be expressed as follows:</p>
<p><disp-formula id="eqn-27"><label>(27)</label><mml:math id="mml-eqn-27" display="block"><mml:mrow><mml:mtext>MDM</mml:mtext></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mover><mml:mi>y</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mrow><mml:mover><mml:mi>w</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow></mml:mfrac><mml:msqrt><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mi>h</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mrow><mml:mi>h</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>h</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>n</mml:mi></mml:mfrac></mml:msqrt></mml:math></disp-formula></p>
<p><disp-formula id="eqn-28"><label>(28)</label><mml:math id="mml-eqn-28" display="block"><mml:msup><mml:mrow><mml:mover><mml:mi>w</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>r</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mrow><mml:mo>&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>h</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mover><mml:mi>r</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><disp-formula id="eqn-29"><label>(29)</label><mml:math id="mml-eqn-29" display="block"><mml:msub><mml:mrow><mml:mover><mml:mi>r</mml:mi><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>n</mml:mi></mml:mfrac><mml:msubsup><mml:mrow><mml:mo>&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mi>j</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>y</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>y</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>h</italic> is the number of steps before prediction.</p>
</sec>
</sec>
<sec id="s4">
<label>4</label><title>Application Case Analysis</title>
<sec id="s4_1">
<label>4.1</label><title>Analysis of Statistical Characteristic of Power Load Measured Value</title>
<p>The power load data in this paper comes from a distribution station of a city in Zhejiang Province, China. It is the real load data obtained based on the intelligent meter advanced measurement system. Sampling is conducted every 15 min to obtain 35,040 load data in 2021. According to the above data, the correlation statistical analysis of power load can be carried out by Matlab/Simulink software and its programming. <xref ref-type="fig" rid="fig-1">Fig. 1a</xref> shows the change trend of power load in 2021, <xref ref-type="fig" rid="fig-1">Fig. 1b</xref> shows the partial diagram of power load at a certain stage in 2021, and <xref ref-type="fig" rid="fig-2">Fig. 2</xref> shows the frequency distribution histogram of power load in different intervals.</p>
<fig id="fig-1">
<label>Figure 1</label>
<caption><title>Variation trend of power load in 2021</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-1.tif"/>
</fig><fig id="fig-2">
<label>Figure 2</label>
<caption><title>Frequency distribution histogram of power load</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-2.tif"/>
</fig>
<p>As can be seen from <xref ref-type="fig" rid="fig-1">Fig. 1a</xref>, the change of power load shows strong nonlinearity and seasonality, that is, the power load in summer is significantly higher than that in other seasons, which is basically consistent with the power load in Southern China. It can be seen from <xref ref-type="fig" rid="fig-1">Fig. 1b</xref> that the short-term change of power load has obvious periodicity. The periodicity of short-term load can be used to predict the short-term load change in the future, so as to formulate a reasonable power grid dispatching plan.</p>

<p>It can be seen from <xref ref-type="fig" rid="fig-2">Fig. 2</xref> that the power load frequency in different intervals is different, and the smaller the interval, the better the probability density curve fitting, but the interval has certain limitations. This also shows that it is very necessary to divide the fuzzy interval reasonably when using fuzzy time series to predict power load.</p>
</sec>
<sec id="s4_2">
<label>4.2</label><title>Analysis of Power Load Forecasting Results</title>
<p>The simulation platform uses Matlab/Simulink software, and the programming language uses m language. Based on this, the relevant algorithm programs are written in this paper. The load forecasting results of one day and one week are observed respectively, and the short-term power load forecasting results under different algorithms are compared and analyzed. <xref ref-type="fig" rid="fig-3">Fig. 3</xref> shows the prediction results and errors under the typical fuzzy time series in a day, <xref ref-type="fig" rid="fig-4">Fig. 4</xref> shows the prediction results and errors under the clustering improved fuzzy time series in a day. <xref ref-type="fig" rid="fig-5">Fig. 5</xref> shows the prediction results and errors under the typical fuzzy time series in a week, <xref ref-type="fig" rid="fig-6">Fig. 6</xref> shows the prediction results and errors under the clustering improved fuzzy time series in a week. <xref ref-type="table" rid="table-2">Table 2</xref> shows the MSE, RMSE, MRE and NS of the two algorithms at different sampling times. The loss functions used in the MDM test follow the metric functions defined earlier in this paper, namely mean square error (MSE), root mean square error (RMSE), and mean relative error (MRE). MDM statistic values of the two algorithms under different sampling times are shown in <xref ref-type="table" rid="table-3">Table 3</xref>.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption><title>FTS prediction results and errors in a day</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-3.tif"/>
</fig><fig id="fig-4">
<label>Figure 4</label>
<caption><title>IFTS prediction results and errors in a day</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-4.tif"/>
</fig><fig id="fig-5">
<label>Figure 5</label>
<caption><title>FTS prediction results and errors in a week</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-5.tif"/>
</fig><fig id="fig-6">
<label>Figure 6</label>
<caption><title>IFTS prediction results and errors in a week</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-6.tif"/>
</fig><table-wrap id="table-2">
<label>Table 2</label>
<caption><title>Evaluation index value under two methods</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th rowspan="2">Index</th>
<th align="center" colspan="2">Day</th>
<th align="center" colspan="2">Week</th>
</tr>
<tr>
<th>FTS</th>
<th>IFTS</th>
<th>FTS</th>
<th>IFTS</th>
</tr>
</thead>
<tbody>
<tr>
<td><italic>MSE</italic></td>
<td>97.778</td>
<td>33.598</td>
<td>76.218</td>
<td>32.102</td>
</tr>
<tr>
<td><italic>RMSE</italic></td>
<td>9.882</td>
<td>5.796</td>
<td>8.730</td>
<td>5.665</td>
</tr>
<tr>
<td><italic>MRE</italic></td>
<td>0.036</td>
<td>0.027</td>
<td>0.043</td>
<td>0.033</td>
</tr>
<tr>
<td><italic>NS</italic></td>
<td>0.991</td>
<td>0.995</td>
<td>0.978</td>
<td>0.990</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-3">
<label>Table 3</label>
<caption><title>MDM statistic values under two methods</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Index</th>
<th>MDM-MSE</th>
<th>MDM-RMSE</th>
<th>MDM-MRE</th>
</tr>
</thead>
<tbody>
<tr>
<td>Day</td>
<td>&#x2212;2.6966&#x002A;&#x002A;&#x002A;</td>
<td>&#x2212;3.8436&#x002A;&#x002A;&#x002A;</td>
<td>&#x2212;3.3300&#x002A;&#x002A;&#x002A;</td>
</tr>
<tr>
<td>Week</td>
<td>&#x2212;5.8875&#x002A;&#x002A;&#x002A;</td>
<td>&#x2212;7.1704&#x002A;&#x002A;&#x002A;</td>
<td>&#x2212;6.9300&#x002A;&#x002A;&#x002A;</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>It can be seen from <xref ref-type="fig" rid="fig-3">Fig. 3</xref> that the daily power load prediction error of typical fuzzy time series is between &#x2212;50 and 20. As can be seen from <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, the prediction error of daily power load in clustering improved fuzzy time series is between &#x2212;20 and 15. It shows that the prediction algorithm designed in this paper is reasonable and can improve the prediction accuracy to a certain extent.</p>

<p>It can be seen from <xref ref-type="fig" rid="fig-5">Fig. 5</xref> that the weekly power load prediction error of typical fuzzy time series is between &#x2212;50 and 30. As can be seen from <xref ref-type="fig" rid="fig-6">Fig. 6</xref>, the prediction error of weekly power load in clustering improved fuzzy time series is between &#x2212;20 and 25. Both algorithms can effectively predict the trend of power load change on a one-week time scale, but the prediction accuracy of the clustering improved fuzzy time series is higher than that of the typical fuzzy time series, indicating that the proposed algorithm in this paper is reasonable and effective.</p>
<p>It can be seen from <xref ref-type="table" rid="table-2">Table 2</xref> that the evaluation indexes <italic>MSE</italic>, <italic>RMSE</italic> and <italic>MRE</italic> of the fuzzy time series prediction method improved by clustering are less than those of the typical fuzzy time series prediction method. The evaluation index <italic>NS</italic> of the fuzzy time series prediction method improved by clustering is better than the typical fuzzy time series prediction method, which further reflects the accuracy of the prediction results.</p>

<p>In <xref ref-type="table" rid="table-3">Table 3</xref>, &#x002A;, &#x002A;&#x002A; and &#x002A;&#x002A;&#x002A; indicate that the null hypothesis is rejected at the significance level of 10%, 5% and 1%, respectively, that is, &#x002A; represents <italic>p</italic> &#x003C; 0.1, &#x002A;&#x002A; represents <italic>p</italic> &#x003C; 0.05, &#x002A;&#x002A;&#x002A; represents <italic>p</italic>&#x00A0;&#x003C;&#x00A0;0.01.</p>

<p>As can be seen from <xref ref-type="table" rid="table-3">Table 3</xref>, at the time scale of one day and one week, various MDM statistical values between FTS and IFTS models do not accept the null hypothesis at the significance level of 1%. It indicates that the prediction accuracy of FTS and IFTS models is significantly different. Furthermore, the IFTS model can improve the prediction accuracy, but there is still much room for improvement. In the subsequent study, other clustering algorithms will be considered to further improve the prediction performance of fuzzy time series.</p>

</sec>
<sec id="s4_3">
<label>4.3</label><title>Comparison with Other Load Forecasting Methods</title>
<p>Typical time series forecasting methods include moving average model (MA), autoregressive model (AR), autoregressive moving average model (ARMA) and so on. Fuzzy time series prediction method is a kind of combined form, and its performance is generally better than that of a single time series prediction method, so it is not compared with this kind of algorithm.</p>
<p>There are two parts that affect the prediction performance of fuzzy time series: 1. The division criterion of the domain; 2. Defuzzification method. In [<xref ref-type="bibr" rid="ref-44">44</xref>&#x2013;<xref ref-type="bibr" rid="ref-46">46</xref>], the central method is used as the main method of defuzzification, so the comparison experiment with this kind of algorithm will be more practical significance. <xref ref-type="fig" rid="fig-7">Figs. 7</xref> and <xref ref-type="fig" rid="fig-8">8</xref> show the load forecasting results and errors of fuzzy time series based on the central defuzzification method (CD-FTS) on the time scale of one day and one week.</p>
<fig id="fig-7">
<label>Figure 7</label>
<caption><title>CD-FTS prediction results and errors in a day</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-7.tif"/>
</fig><fig id="fig-8">
<label>Figure 8</label>
<caption><title>CD-FTS prediction results and errors in a week</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMES_25396-fig-8.tif"/>
</fig>
<p>It can be seen from <xref ref-type="fig" rid="fig-7">Figs. 7</xref> and <xref ref-type="fig" rid="fig-8">8</xref> that the CD-FTS method can effectively predict the short-term variation trend of load, and the prediction error range on the time scale of one day and one week is [&#x2212;80, 40]. The prediction error of the CD-FTS is much larger than that of the IFTS, which further demonstrates the superiority of the proposed algorithm in this paper.</p>

</sec>
</sec>
<sec id="s5">
<label>5</label><title>Conclusion</title>
<p>This paper presents a fuzzy time series prediction method based on K-means clustering and applies it to the actual case of short-term power load forecasting in distribution stations. The following conclusions are obtained through Matlab/Simulink simulation case analysis: (1) The fuzzy time series prediction method improved by clustering has a more reasonable division of the fuzzy theory domain. (2) The fuzzy time series expressed by higher-order fuzzy relation can better fit the nonlinear relation between data. (3) The load forecasting evaluation indexes MSE/RMSE/MRE/NS of typical fuzzy time series on the time scale of one day are 97.778/9.882/0.036/0.991. The corresponding values of the improved fuzzy time series are 33.598/5.796/0027/0.995. The MSE/RMSE/MRE value of the typical fuzzy time series is greater than the corresponding value of the improved fuzzy time series, while the NS value of the typical fuzzy time series is less than the corresponding value of the improved fuzzy time series. It shows that the improved fuzzy time series method can improve the forecasting accuracy of power load. (4) Similarly, the load forecasting evaluation indexes MSE/RMSE/MRE of the typical fuzzy time series is larger than the corresponding value of the improved fuzzy time series in the one-week time scale, while the NS value is small. It further shows that the design of the improved fuzzy time series algorithm is reasonable and effective. (5) The null hypothesis is not accepted for various MDM statistic values between fuzzy time series and improved fuzzy time series at the 1% level of significance. The MDM values of the two algorithms are different, indicating that their prediction performance is different.</p>
<p>In addition, the prediction algorithm proposed in this paper can be applied to prediction research in other fields, and it has certain universality. However, the factors affecting power load change are not single, such as temperature, humidity and other environmental factors. Whether the clustering improved fuzzy time series prediction method proposed in this paper is still accurate in power load prediction under multiple influencing factors needs to be further verified in future research.</p>
</sec>
</body>
<back>
<sec><title>Funding Statement</title>
<p>This work was supported by the <funding-source>National Natural Science Foundation of China</funding-source> under Grant <award-id>51777193</award-id>.</p>
</sec>
<sec sec-type="COI-statement"><title>Conflicts of Interest</title>
<p>The authors declare that they have no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear"><title>References</title>
<ref id="ref-1"><label>1.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ali</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Adnan</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Tariq</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Optimum control strategies for short term load forecasting in smart grids</article-title>. <source>International Journal of Electrical Power &#x0026; Energy Systems</source><italic>,</italic> <volume>113</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>792</fpage>&#x2013;<lpage>806</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2019.06.010</pub-id></mixed-citation></ref>
<ref id="ref-2"><label>2.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lu</surname>, <given-names>Y. T.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>G. C.</given-names></string-name>, <string-name><surname>Huang</surname>, <given-names>S. Q.</given-names></string-name></person-group> (<year>2022</year>). <article-title>A short-term load forecasting model based on mixup and transfer learning</article-title>. <source>Electric Power Systems Research</source><italic>,</italic> <volume>207</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.107837</pub-id></mixed-citation></ref>
<ref id="ref-3"><label>3.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ilyas</surname>, <given-names>O.</given-names></string-name>, <string-name><surname>Berat</surname>, <given-names>E. S.</given-names></string-name>, <string-name><surname>Harun</surname>, <given-names>O.</given-names></string-name></person-group> (<year>2021</year>). <article-title>A combined deep learning application for short term load forecasting</article-title>. <source>Alexandria Engineering Journal</source><italic>,</italic> <volume>60</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>3807</fpage>&#x2013;<lpage>3818</lpage>. <pub-id pub-id-type="doi">10.1016/j.aej.2021.02.050</pub-id></mixed-citation></ref>
<ref id="ref-4"><label>4.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yang</surname>, <given-names>X. D.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>Z. Y.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>Y. B.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Resilience-oriented co-deployment of remote-controlled switches and soft open point in distribution networks</article-title>. <source>IEEE Transactions on Power Systems</source><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2022.3176024</pub-id></mixed-citation></ref>
<ref id="ref-5"><label>5.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name>Jalali, S. M. J.</string-name>, <string-name><surname>Ahmadian</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Khosravi</surname>, <given-names>A.</given-names></string-name></person-group> (<year>2021</year>). <article-title>A novel evolutionary-based deep convolutional neural network model for intelligent load forecasting</article-title>. <source>IEEE Transactions on Industrial Informatics</source><italic>,</italic> <volume>17</volume><italic>(</italic><issue>12</issue><italic>),</italic> <fpage>8243</fpage>&#x2013;<lpage>8253</lpage>. <pub-id pub-id-type="doi">10.1109/TII.2021.3065718</pub-id></mixed-citation></ref>
<ref id="ref-6"><label>6.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lin</surname>, <given-names>W. X.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Benoit</surname>, <given-names>B.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Spatial-temporal residential short-term load forecasting via graph neural networks</article-title>. <source>IEEE Transactions on Smart Grid</source><italic>,</italic> <volume>12</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>5373</fpage>&#x2013;<lpage>5384</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2021.3093515</pub-id></mixed-citation></ref>
<ref id="ref-7"><label>7.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yang</surname>, <given-names>X. D.</given-names></string-name>, <string-name><surname>Xu</surname>, <given-names>C. B.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>Y. B.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Real-time coordinated scheduling for ADNs with soft open points and charging stations</article-title>. <source>IEEE Transactions on Power Systems</source><italic>,</italic> <volume>36</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>5486</fpage>&#x2013;<lpage>5499</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2021.3070036</pub-id></mixed-citation></ref>
<ref id="ref-8"><label>8.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yang</surname>, <given-names>X. D.</given-names></string-name>, <string-name><surname>Xu</surname>, <given-names>C. B.</given-names></string-name>, <string-name><surname>He</surname>, <given-names>H. B.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Flexibility provisions in active distribution networks with uncertainties</article-title>. <source>IEEE Transactions on Sustainable Energy</source><italic>,</italic> <volume>12</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>553</fpage>&#x2013;<lpage>567</lpage>.</mixed-citation></ref>
<ref id="ref-9"><label>9.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zong</surname>, <given-names>W. T.</given-names></string-name>, <string-name><surname>Wei</surname>, <given-names>Z. N.</given-names></string-name>, <string-name><surname>Sun</surname>, <given-names>G. Q.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Short-term load interval prediction based on improved gaussian process regression model</article-title>. <source>Proceeding of the CSU-EPSA</source><italic>,</italic> <volume>29</volume><italic>(</italic><issue>8</issue><italic>),</italic> <fpage>22</fpage>&#x2013;<lpage>28</lpage>.</mixed-citation></ref>
<ref id="ref-10"><label>10.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lu</surname>, <given-names>X. S.</given-names></string-name>, <string-name><surname>Pan</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>K.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>K.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Markov modified grey-time series electric load forecasting method</article-title>. <source>Automation Technology and Application</source><italic>,</italic> <volume>41</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>132</fpage>&#x2013;<lpage>136</lpage>.</mixed-citation></ref>
<ref id="ref-11"><label>11.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Shalini</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Angshul</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Victor</surname>, <given-names>E.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Blind Kalman filtering for short-term load forecasting</article-title>. <source>IEEE Transactions on Power System</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>4916</fpage>&#x2013;<lpage>4919</lpage>.</mixed-citation></ref>
<ref id="ref-12"><label>12.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Fan</surname>, <given-names>G. F.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>L. Z.</given-names></string-name>, <string-name><surname>Yu</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Applications of random forest in multivariable response surface for short-term load forecasting</article-title>. <source>International Journal of Electrical Power and Energy Systems</source><italic>,</italic> <volume>139</volume><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2022.108073</pub-id></mixed-citation></ref>
<ref id="ref-13"><label>13.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname>, <given-names>X. L.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>J.</given-names></string-name></person-group> (<year>2022</year>). <article-title>An adaptive hybrid fractal model for short-term load forecasting in power systems</article-title>. <source>Electric Power Systems Research</source><italic>,</italic> <volume>207</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.107858</pub-id></mixed-citation></ref>
<ref id="ref-14"><label>14.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname>, <given-names>H. M.</given-names></string-name>, <string-name><surname>Tang</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Pu</surname>, <given-names>Y.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Short-term load forecasting of multi-energy in integrated energy system based on multivariate phase space reconstruction and support vector regression mode</article-title>. <source>Electric Power Systems Research</source><italic>,</italic> <volume>210</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.108066</pub-id></mixed-citation></ref>
<ref id="ref-15"><label>15.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Cui</surname>, <given-names>M. J.</given-names></string-name>, <string-name><surname>Ke</surname>, <given-names>D. P.</given-names></string-name>, <string-name><surname>Gan</surname>, <given-names>D.</given-names></string-name></person-group> (<year>2015</year>). <article-title>Statistical scenarios forecasting method for wind power ramp events using modified neural networks</article-title>. <source>Journal of Modern Power Systems and Clean Energy</source><italic>,</italic> <volume>3</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>371</fpage>&#x2013;<lpage>380</lpage>. <pub-id pub-id-type="doi">10.1007/s40565-015-0138-7</pub-id></mixed-citation></ref>
<ref id="ref-16"><label>16.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Navid</surname>, <given-names>F. M.</given-names></string-name>, <string-name><surname>Katarina</surname>, <given-names>G.</given-names></string-name>, <string-name><surname>Syed</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Distributed load forecasting using smart meter data: Federated learning with recurrent neural networks</article-title>. <source>International Journal of Electrical Power and Energy Systems</source><italic>,</italic> <volume>137</volume><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>12</lpage>.</mixed-citation></ref>
<ref id="ref-17"><label>17.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Hao</surname>, <given-names>Q.</given-names></string-name>, <string-name><surname>Srinivasan</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Khosravi</surname>, <given-names>A.</given-names></string-name></person-group> (<year>2014</year>). <article-title>Short-term load and wind power forecasting using neural network-based prediction intervals</article-title>. <source>IEEE Transactions on Neural Networks and Learning System</source><italic>,</italic> <volume>25</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>303</fpage>&#x2013;<lpage>315</lpage>. <pub-id pub-id-type="doi">10.1109/TNNLS.2013.2276053</pub-id>; <pub-id pub-id-type="pmid">24807030</pub-id></mixed-citation></ref>
<ref id="ref-18"><label>18.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lin</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Ma</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Zhu</surname>, <given-names>J. G.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Short-term load forecasting based on LSTM networks considering attention mechanism</article-title>. <source>International Journal of Electrical Power and Energy Systems</source><italic>,</italic> <volume>137</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.107818</pub-id></mixed-citation></ref>
<ref id="ref-19"><label>19.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tang</surname>, <given-names>X. L.</given-names></string-name>, <string-name><surname>Chen</surname>, <given-names>H. X.</given-names></string-name>, <string-name><surname>Xiang</surname>, <given-names>W. H.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Short-term load forecasting using channel and temporal attention based temporal convolutional network</article-title>. <source>Electric Power Systems Research</source><italic>,</italic> <volume>205</volume><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2021.107761</pub-id></mixed-citation></ref>
<ref id="ref-20"><label>20.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhou</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Jin</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Holographic ensemble forecasting method for short-term power load</article-title>. <source>IEEE Transactions on Smart Grid</source><italic>,</italic> <volume>10</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>425</fpage>&#x2013;<lpage>434</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2017.2743015</pub-id></mixed-citation></ref>
<ref id="ref-21"><label>21.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ouyang</surname>, <given-names>T. H.</given-names></string-name>, <string-name><surname>He</surname>, <given-names>Y. S.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>H. J.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Modeling and forecasting short-term power load with copula model and deep belief network</article-title>. <source>IEEE Transactions on Emerging Topics in Computational Intelligence</source><italic>,</italic> <volume>3</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>127</fpage>&#x2013;<lpage>136</lpage>. <pub-id pub-id-type="doi">10.1109/TETCI.2018.2880511</pub-id></mixed-citation></ref>
<ref id="ref-22"><label>22.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Alhmoud</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Nawafleh</surname>, <given-names>Q.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Short-term load forecasting for Jordan power system based on NARX-ELMAN neural network and ARMA model</article-title>. <source>IEEE Canadian Journal of Electrical and Computer Engineering</source><italic>,</italic> <volume>44</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>356</fpage>&#x2013;<lpage>363</lpage>. <pub-id pub-id-type="doi">10.1109/ICJECE.2021.3076124</pub-id></mixed-citation></ref>
<ref id="ref-23"><label>23.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Gao</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Hug</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Personalized federated learning for individual consumer load forecasting</article-title>. <source>CSEE Journal of Power and Energy Systems</source>. <pub-id pub-id-type="doi">10.17775/CSEEJPES.2021.07350</pub-id></mixed-citation></ref>
<ref id="ref-24"><label>24.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>An</surname>, <given-names>Y. F.</given-names></string-name>, <string-name><surname>Zhai</surname>, <given-names>X. Q.</given-names></string-name></person-group> (<year>2022</year>). <article-title>SVR-DEA model of carbon tax pricing for China&#x2019;s thermal power industry</article-title>. <source>Science of the Total Environment</source><italic>,</italic> <volume>734</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2020.139438</pub-id>; <pub-id pub-id-type="pmid">32460083</pub-id></mixed-citation></ref>
<ref id="ref-25"><label>25.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Kazem</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Sharifi</surname>, <given-names>E.</given-names></string-name>, <string-name><surname>Hussain</surname>, <given-names>F. K.</given-names></string-name></person-group> (<year>2013</year>). <article-title>Support vector regression with chaos-based firefly algorithm for stock market price forecasting</article-title>. <source>Applied Soft Computing</source><italic>,</italic> <volume>13</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>947</fpage>&#x2013;<lpage>958</lpage>. <pub-id pub-id-type="doi">10.1016/j.asoc.2012.09.024</pub-id></mixed-citation></ref>
<ref id="ref-26"><label>26.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lin</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Lin</surname>, <given-names>Z. X.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Forecasting the realized volatility of stock price index: A hybrid model integrating CEEMDAN and LSTM</article-title>. <source>Expert Systems with Applications</source><italic>,</italic> <volume>206</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>117736</fpage>. <pub-id pub-id-type="doi">10.1016/j.eswa.2022.117736</pub-id></mixed-citation></ref>
<ref id="ref-27"><label>27.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liang</surname>, <given-names>Y. H.</given-names></string-name>, <string-name><surname>Lin</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Lu</surname>, <given-names>Q.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Forecasting gold price using a novel hybrid model with ICEEMDAN and LSTM-CNN-CBAM</article-title>. <source>Expert Systems with Applications</source><italic>,</italic> <volume>206</volume><italic>(</italic><issue>10</issue><italic>),</italic> <fpage>117847</fpage>. <pub-id pub-id-type="doi">10.1016/j.eswa.2022.117847</pub-id></mixed-citation></ref>
<ref id="ref-28"><label>28.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname>, <given-names>L. H.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>B.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Energy market prediction with novel long short-term memory network: Case study of energy futures index volatility</article-title>. <source>Energy</source><italic>,</italic> <volume>211</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>118634</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2020.118634</pub-id></mixed-citation></ref>
<ref id="ref-29"><label>29.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lin</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Liao</surname>, <given-names>Q. D.</given-names></string-name>, <string-name><surname>Lin</surname>, <given-names>Z. X.</given-names></string-name></person-group> (<year>2022</year>). <article-title>A novel hybrid model integrating modified ensemble empirical mode decomposition and LSTM neural network for multi-step precious metal prices prediction</article-title>. <source>Resources Policy</source><italic>,</italic> <volume>78</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>102884</fpage>. <pub-id pub-id-type="doi">10.1016/j.resourpol.2022.102884</pub-id></mixed-citation></ref>
<ref id="ref-30"><label>30.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Song</surname>, <given-names>Q.</given-names></string-name>, <string-name><surname>Chissom</surname>, <given-names>B. S.</given-names></string-name></person-group> (<year>1993</year>). <article-title>Forecasting enrollments with fuzzy time series&#x2014;Part I</article-title>. <source>Fuzzy Sets and Systems</source><italic>,</italic> <volume>54</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1016/0165-0114(93)90355-L</pub-id></mixed-citation></ref>
<ref id="ref-31"><label>31.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname>, <given-names>S. M.</given-names></string-name></person-group> (<year>1996</year>). <article-title>Forecasting enrollments based on fuzzy time series</article-title>. <source>Fuzzy Sets and Systems</source><italic>,</italic> <volume>81</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>311</fpage>&#x2013;<lpage>319</lpage>. <pub-id pub-id-type="doi">10.1016/0165-0114(95)00220-0</pub-id></mixed-citation></ref>
<ref id="ref-32"><label>32.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lee</surname>, <given-names>M. H.</given-names></string-name>, <string-name><surname>Efendi</surname>, <given-names>R.</given-names></string-name>, <string-name><surname>Ismail</surname>, <given-names>Z.</given-names></string-name></person-group> (<year>2009</year>). <article-title>Modified weighted for enrollment forecasting based on fuzzy time series</article-title>. <source>MATEMATIKA: Malaysian Journal of Industrial and Applied Mathematics</source><italic>,</italic> <volume>25</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>67</fpage>&#x2013;<lpage>78</lpage>.</mixed-citation></ref>
<ref id="ref-33"><label>33.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Huarng</surname>, <given-names>K. H.</given-names></string-name></person-group> (<year>2006</year>). <article-title>Ratio-based lengths of intervals to improve fuzzy time series forecasting</article-title>. <source>IEEE Transactions on Systems, Man, and Cybernetics&#x2013;Part B: Cybernetics</source><italic>,</italic> <volume>36</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>328</fpage>&#x2013;<lpage>340</lpage>. <pub-id pub-id-type="doi">10.1109/TSMCB.2005.857093</pub-id>; <pub-id pub-id-type="pmid">16602593</pub-id></mixed-citation></ref>
<ref id="ref-34"><label>34.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tahseen</surname>, <given-names>A. J.</given-names></string-name>, <string-name><surname>Syed</surname>, <given-names>M. A. B.</given-names></string-name></person-group> (<year>2008</year>). <article-title>A refined fuzzy time series model for stock market forecasting</article-title>. <source>Physica A: Statistical Mechanics and its Applications</source><italic>,</italic> <volume>387</volume><italic>(</italic><issue>12</issue><italic>),</italic> <fpage>2857</fpage>&#x2013;<lpage>2862</lpage>. <pub-id pub-id-type="doi">10.1016/j.physa.2008.01.099</pub-id></mixed-citation></ref>
<ref id="ref-35"><label>35.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Aladag</surname>, <given-names>C. H.</given-names></string-name>, <string-name><surname>Yolcu</surname>, <given-names>U.</given-names></string-name>, <string-name><surname>Egrioglu</surname>, <given-names>E.</given-names></string-name></person-group> (<year>2010</year>). <article-title>A high order fuzzy time series forecasting model based on adaptive expectation and artificial neural network</article-title>. <source>Mathematics and Computers in Simulation</source><italic>,</italic> <volume>81</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>875</fpage>&#x2013;<lpage>882</lpage>. <pub-id pub-id-type="doi">10.1016/j.matcom.2010.09.011</pub-id></mixed-citation></ref>
<ref id="ref-36"><label>36.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Egrioglu</surname>, <given-names>E.</given-names></string-name></person-group> (<year>2010</year>). <article-title>Finding an optimal interval length in high order fuzzy time series</article-title>. <source>Expert Systems with Applications</source><italic>,</italic> <volume>37</volume><italic>(</italic><issue>7</issue><italic>),</italic> <fpage>5052</fpage>&#x2013;<lpage>5055</lpage>. <pub-id pub-id-type="doi">10.1016/j.eswa.2009.12.006</pub-id></mixed-citation></ref>
<ref id="ref-37"><label>37.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Egrioglu</surname>, <given-names>E.</given-names></string-name></person-group> (<year>2011</year>). <article-title>A new approach based on the optimization of the length of intervals in fuzzy time series</article-title>. <source>Journal of Intelligent &#x0026; Fuzzy Systems</source><italic>,</italic> <volume>22</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>15</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.3233/IFS-2010-0470</pub-id></mixed-citation></ref>
<ref id="ref-38"><label>38.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>R.</given-names></string-name>, <string-name><surname>Gao</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>J. L.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Electric vehicle charging demand forecasting method based on clustering analysis</article-title>. <source>Power System Protection and Control</source><italic>,</italic> <volume>48</volume><italic>(</italic><issue>16</issue><italic>),</italic> <fpage>37</fpage>&#x2013;<lpage>44</lpage>.</mixed-citation></ref>
<ref id="ref-39"><label>39.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yan</surname>, <given-names>L. P.</given-names></string-name>, <string-name><surname>Hong</surname>, <given-names>W. C.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Evaluation and forecasting of wind energy investment risk along the belt and road based on a novel hybrid intelligent model</article-title>. <source>Computer Modeling in Engineering &#x0026; Sciences</source><italic>,</italic> <volume>128</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>1069</fpage>&#x2013;<lpage>1102</lpage>. <pub-id pub-id-type="doi">10.32604/cmes.2021.016499</pub-id></mixed-citation></ref>
<ref id="ref-40"><label>40.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname>, <given-names>G.</given-names></string-name>, <string-name><surname>Ding</surname>, <given-names>H. L.</given-names></string-name></person-group> (<year>2018</year>). <article-title>Stabilization algorithm of fuzzy time series based on principal component analysis</article-title>. <source>Control and Decision</source><italic>,</italic> <volume>33</volume><italic>(</italic><issue>9</issue><italic>),</italic> <fpage>1643</fpage>&#x2013;<lpage>1648</lpage>.</mixed-citation></ref>
<ref id="ref-41"><label>41.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Xian</surname>, <given-names>S. D.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>T. J.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Fuzzy time series prediction model based on improved wolf pack algorithm</article-title>. <source>Control Theory &#x0026; Applications</source><italic>,</italic> <volume>37</volume><italic>(</italic><issue>7</issue><italic>),</italic> <fpage>1638</fpage>&#x2013;<lpage>1643</lpage>.</mixed-citation></ref>
<ref id="ref-42"><label>42.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Zhu</surname>, <given-names>H.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Fuzzy segmentation of multivariate time series with KPCA and G-G clustering</article-title>. <source>Control and Decision</source><italic>,</italic> <volume>36</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>115</fpage>&#x2013;<lpage>124</lpage>.</mixed-citation></ref>
<ref id="ref-43"><label>43.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Lin</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Lu</surname>, <given-names>Q.</given-names></string-name>, <string-name><surname>Tan</surname>, <given-names>B.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Forecasting energy prices using a novel hybrid model with variational mode decomposition</article-title>. <source>Energy</source><italic>,</italic> <volume>246</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>123366</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2022.123366</pub-id></mixed-citation></ref>
<ref id="ref-44"><label>44.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zeng</surname>, <given-names>D. L.</given-names></string-name>, <string-name><surname>Lu</surname>, <given-names>J. Y.</given-names></string-name>, <string-name><surname>Zheng</surname>, <given-names>Y. F.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Combined fuzzy time series prediction method for fault prediction of EML pulse capacitors</article-title>. <source>IEEE Transactions on Plasma Science</source><italic>,</italic> <volume>49</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>905</fpage>&#x2013;<lpage>913</lpage>. <pub-id pub-id-type="doi">10.1109/TPS.2020.3029840</pub-id></mixed-citation></ref>
<ref id="ref-45"><label>45.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ahmed</surname>, <given-names>T. S.</given-names></string-name>, <string-name><surname>Patrick</surname>, <given-names>J. N.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Heuristic hidden Markov model for fuzzy time series forecasting</article-title>. <source>International Journal of Intelligent Systems Technologies and Applications</source><italic>,</italic> <volume>20</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>146</fpage>&#x2013;<lpage>166</lpage>. <pub-id pub-id-type="doi">10.1504/IJISTA.2021.119030</pub-id></mixed-citation></ref>
<ref id="ref-46"><label>46.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yousif</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Mahmod</surname>, <given-names>O.</given-names></string-name>, <string-name><surname>Akram</surname>, <given-names>A. A.</given-names></string-name></person-group> (<year>2021</year>). <article-title>A novel stochastic fuzzy time series forecasting model based on a new partition method</article-title>. <source>IEEE Access</source><italic>,</italic> <volume>9</volume><italic>,</italic> <fpage>80236</fpage>&#x2013;<lpage>80252</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2021.3084048</pub-id></mixed-citation></ref>
</ref-list>
</back>
</article>