CN106096771A - 支持向量机的参数选取方法和系统与负荷预测方法和系统 - Google Patents
支持向量机的参数选取方法和系统与负荷预测方法和系统 Download PDFInfo
- Publication number
- CN106096771A CN106096771A CN201610402774.6A CN201610402774A CN106096771A CN 106096771 A CN106096771 A CN 106096771A CN 201610402774 A CN201610402774 A CN 201610402774A CN 106096771 A CN106096771 A CN 106096771A
- Authority
- CN
- China
- Prior art keywords
- support vector
- vector machine
- load data
- variable
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012706 support-vector machine Methods 0.000 title claims abstract description 164
- 238000010187 selection method Methods 0.000 title claims abstract description 22
- 238000013277 forecasting method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 42
- 230000035772 mutation Effects 0.000 claims description 21
- 240000002853 Nelumbo nucifera Species 0.000 claims description 6
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 6
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 8
- 230000003044 adaptive effect Effects 0.000 abstract description 4
- 239000000428 dust Substances 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Genetics & Genomics (AREA)
- Primary Health Care (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种支持向量机的参数选取方法和系统以及负荷预测方法和系统,其是采用云遗传算法对支持向量机的参数进行了优化,根据正态云模型输出结果的随机性和稳定倾向性特点,随机性可以保持个体多样性从而避免搜索陷入局部极值,而稳定倾向性又可以很好地保护较优个体从而对全局最值进行自适应定位,有效提高搜索到最优解或满意解的速度,从而提高以选取的参数配置的支持向量机进行负荷预测的效率,而且在选取过程中获取了预测负荷数据与实际负荷数据之间的误差值以作为云遗传算法优化的效果,优化后使误差值将至最低,以此提高选取的参数的准确性。
Description
技术领域
本发明涉及电力系统技术领域,特别是涉及一种支持向量机的参数选取方法和系统与负荷预测方法和系统。
背景技术
科学的预测是正确决策的依据和保证。负荷预测是电力系统领域的一个传统研究问题,是指从已知的电力系统、经济、社会、气象等情况出发,通过对历史数据的分析和研究,探索事物之间的内在联系和发展变化规律,对负荷发展做出预先估计和推测。负荷预测是电力系统规划、计划、用电、调度等部门的基础工作,其重要性早已被人们所认识。
负荷预测本质上是对功率曲线进行拟合与回归,由于实时功率曲线受电力系统、经济、社会、气象等诸多因素影响,一般表现为复杂非线性特点,宜采用对复杂非线性特性具备较强学习能力的预测模型。
目前应用较多且比较成熟的预测方法,主要为支持向量机(Support VectorMachine,SVM)等。SVM同时考虑经验风险最小和结构风险最小,使模型具有较强的推广性,在小样本识别方面有较大优势,且SVM有严格的数学理论基础,其决策为全局最优。
SVM参数的选取策略目前尚没有统一的方法,SVM参数选取的优劣将直接影响模型的拟合和回归能力。现有技术中,较为常用的SVM参数优化算法包括网格搜索算法、粒子群算法等。利用这些算法虽然可以选取SVM参数,但搜索到最优解或满意解的速度太慢,依据选取的参数进行负荷预测的效率低。
发明内容
基于此,有必要针对现有的支持向量机搜索到最优解或满意解的速度太慢导致负荷预测效率低的问题,提供一种支持向量机的参数选取方法和系统与负荷预测方法和系统。
一种支持向量机的参数选取方法,包括以下步骤:
获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,支持向量机的各参数对应云遗传算法的各自变量;
根据对应的各参数的取值范围对各自变量均匀取值,形成初始种群,其中,初始种群中包括多个个体,个体包括取值后的各自变量;
针对每个个体,根据取值后的各自变量配置支持向量机,获取实时负荷数据和历史负荷数据以及实时负荷数据和历史负荷数据之间的相似度,从历史负荷数据中筛选出相似日的负荷数据,其中,相似日的负荷数据与实时负荷数据的相似度大于预设值;将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以训练支持向量机;将相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值;
选取初始种群中的部分个体,以选取的各个体对应的误差值的倒数作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
判断第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各自变量作为支持向量机的各参数;若否,则返回至针对每个个体,根据取值后的各自变量配置支持向量机的步骤。
一种支持向量机的参数选取系统,包括以下单元:
第一获取单元,用于获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,支持向量机的各参数对应云遗传算法的各自变量;
生成单元,用于根据对应的各参数的取值范围对各自变量均匀取值,形成初始种群,其中,初始种群中包括多个个体,个体包括取值后的各自变量;
第一训练单元,用于针对每个个体,根据取值后的各自变量配置支持向量机,获取实时负荷数据和历史负荷数据以及实时负荷数据和历史负荷数据之间的相似度,从历史负荷数据中筛选出相似日的负荷数据,其中,相似日的负荷数据与实时负荷数据的相似度大于预设值;将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以训练支持向量机;
第二获取单元,将相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值;
云操作单元,选取初始种群中的部分个体,以选取的各个体对应的误差值的倒数作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
判断单元,判断第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各自变量作为支持向量机的各参数;若否,则返回至第一训练单元。
根据上述本发明的支持向量机的参数选取方法和系统,其是采用云遗传算法对支持向量机的参数进行了优化,借助于云模型,利用云发生器代替传统的交叉、变异算子,根据正态云模型输出结果的随机性和稳定倾向性特点,随机性可以保持个体多样性从而避免搜索陷入局部极值,而稳定倾向性又可以很好地保护较优个体从而对全局最值进行自适应定位,有效提高搜索到最优解或满意解的速度,进而可以提高负荷预测的效率,而且在选取过程中获取了预测负荷数据与实际负荷数据之间的误差值以作为云遗传算法优化的效果,优化后使误差值将至最低,以此提高选取的参数的准确性。
一种负荷预测方法,包括以下步骤:
根据上述支持向量机的参数选取方法选取支持向量机的参数,根据选取的支持向量机的参数对支持向量机进行配置;
将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以再度训练支持向量机;
将实时负荷数据输入至再度训练后的支持向量机,获得对应实时负荷数据的日期之后的预设天数的负荷数据。
一种负荷预测系统,包括上述支持向量机的参数选取系统和以下单元:
配置单元,用于根据上述的支持向量机的参数选取系统选取的支持向量机的参数对支持向量机进行配置;
第二训练单元,将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以再度训练支持向量机;
第三获取单元,用于将实时负荷数据输入至再度训练后的支持向量机,获得对应实时负荷数据的日期之后的预设天数的负荷数据。
根据上述本发明的负荷预测方法和系统,其是采用云遗传算法对支持向量机的参数进行了优化,借助于云模型,利用云发生器代替传统的交叉、变异算子,根据正态云模型输出结果的随机性和稳定倾向性特点,随机性可以保持个体多样性从而避免搜索陷入局部极值,而稳定倾向性又可以很好地保护较优个体从而对全局最值进行自适应定位,有效提高搜索到最优解或满意解的速度,从而提高以选取的参数配置的支持向量机进行负荷预测的效率,而且在选取参数过程中,获取了训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值以作为云遗传算法优化的效果,优化后使误差值将至最低,以此提高选取的参数的准确性,从而降低支持向量机模型在实际预测时的预测误差,提高负荷预测的精度。
附图说明
图1为其中一个实施例的支持向量机的参数选取方法的流程示意图;
图2为其中一个实施例的负荷预测方法的流程示意图;
图3为其中一个实施例的负荷预测方法实际应用效果图;
图4为其中一个实施例的支持向量机的参数选取系统的结构示意图;
图5为其中一个实施例的负荷预测系统的结构示意图;
图6为其中一个实施例的负荷预测系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步的详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不限定本发明的保护范围。
参见图1所示,为本发明的支持向量机的参数选取方法的流程示意图。该实施例中的支持向量机的参数选取方法,包括以下步骤:
步骤S101:获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,支持向量机的各参数对应云遗传算法的各自变量;
在本步骤中,实质上是将支持向量机的各参数作为云遗传算法的各自变量。
步骤S102:根据对应的各参数的取值范围对各自变量均匀取值,形成初始种群,其中,初始种群中包括多个个体,个体包括取值后的各自变量;
在本步骤中,由于各参数都有对应的取值范围,也就是各自变量都有对应的取值范围,对各自变量进行取值可以形成多个个体。
步骤S103:针对每个个体,根据取值后的各自变量配置支持向量机,获取实时负荷数据和历史负荷数据以及实时负荷数据和历史负荷数据之间的相似度,从历史负荷数据中筛选出相似日的负荷数据,其中,相似日的负荷数据与实时负荷数据的相似度大于预设值;将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以训练支持向量机;将相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值;
在本步骤中,主要是根据个体中的参数值配置支持向量机,利用实时负荷数据和相似日的负荷数据对支持向量机进行训练,并获取误差值。
步骤S104:选取初始种群中的部分个体,以选取的各个体对应的误差值作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
步骤S105:判断第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各自变量作为支持向量机的各参数;若否,则返回至针对每个个体,根据取值后的各自变量配置支持向量机的步骤。
根据上述本发明的支持向量机的参数选取方法,其是采用云遗传算法对支持向量机的参数进行了优化,借助于云模型,利用云发生器代替传统的交叉、变异算子,根据正态云模型输出结果的随机性和稳定倾向性特点,随机性可以保持个体多样性从而避免搜索陷入局部极值,而稳定倾向性又可以很好地保护较优个体从而对全局最值进行自适应定位,有效提高搜索到最优解或满意解的速度,进而可以提高负荷预测的效率,而且在选取过程中获取了预测负荷数据与实际负荷数据之间的误差值以作为云遗传算法优化的效果,优化后使误差值将至最低,以此提高选取的参数的准确性。
在其中一个实施例中,对选取的各个体进行基于条件云的交叉操作的步骤包括以下步骤:
对于任何一组用于交叉的两个个体,对两个个体中的每个变量作以下运算:
根据Ex=Fmax、和计算第一特征数值,其中,Ex、En和He为第一特征数值,Fmax和Fmin分别为选取的各个体对应的适应度的最大值和最小值;
x0取交叉操作的两个个体对应的适应度中的较大值,将x0以及Ex、En和He输入至X条件云发生器,获得第一确定度;
根据和计算第二特征数值,其中,Ex’、En’和He’为第二特征数值,xf和xm分别为交叉操作的两个个体中对应的同一种变量,span为同一种变量的取值区间长度,Ff和Fm分别为两个个体对应的适应度,pop为种群大小;
将第一确定度以及Ex’、En’和He’输入至Y条件云发生器,获得与同一变量对应的第一子变量和第二子变量,用第一子变量替换xf,用第二子变量替换xm。
在本实施例中,交叉操作目的是利用Y条件云发生器对两个个体的信息进行交换,得到一对新个体,为维持种群大小不变,将新得到的一对个体取代原来的父母个体,使得种群中的个体之间的差异减小。
在其中一个实施例中,对第一个体组中的各个体进行基于基本云的变异操作的步骤包括以下步骤:
对于任何一个用于变异操作的个体,对该个体中的任一变量xp作以下运算:
根据和计算第三特征数值,其中,Ex”、En”和He”为第三特征数值,span(xp)表示变量xp的取值区间长度,pop为种群大小;
将Ex”、En”和He”输入至基本云发生器,获得与变量xp对应的第三子变量和第二确定度,生成服从[0,1]区间均匀分布的随机数,当第二确定度大于随机数时,用第三子变量替换变量xp。
在本实施例中,变异操作目的是利用基本云发生器对个体中的变量进行转换,将新的信息引入所选个体。变异操作的对象为经过交叉操作得到的第一个体组中的个体中的变量,在每个个体的变量完成变异后就得到新个体,将新个体取代原个体。
在其中一个实施例中,选取初始种群中的部分个体的步骤包括以下步骤:
基于轮盘赌方法的选择方式选取初始种群中的部分个体。
在本实施例中,轮盘赌选择是从种群中随机选择一些个体的方法,被选中的机率和它们的适应度值成比例,个体的适应度值愈高(也即对应的误差值愈小),被选中的概率也愈大,基于轮盘赌方法的选择方式选取的个体的适应度值较高,可以间接提高支持向量机参数选取的准确性。
根据上述支持向量机的参数选取方法,本发明还提供一种负荷预测方法,以下就本发明的负荷预测方法的实施例进行详细说明。
参见图2所示,为本发明的负荷预测方法的流程示意图。该实施例中的负荷预测方法包括以下步骤:
步骤S201:根据上述支持向量机的参数选取方法选取支持向量机的参数,根据选取的支持向量机的参数对支持向量机进行配置;
步骤S202:将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以再度训练支持向量机;
步骤S203:将实时负荷数据输入至再度训练后的支持向量机,获得对应实时负荷的日期之后的预设天数的负荷数据。
在本实施例中,采用云遗传算法对支持向量机的参数进行了优化,借助于云模型,利用云发生器代替传统的交叉、变异算子,根据正态云模型输出结果的随机性和稳定倾向性特点,随机性可以保持个体多样性从而避免搜索陷入局部极值,而稳定倾向性又可以很好地保护较优个体从而对全局最值进行自适应定位,有效提高搜索到最优解或满意解的速度,从而提高以选取的参数配置的支持向量机进行负荷预测的效率,而且在选取参数过程中,获取了训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值以作为云遗传算法优化的效果,优化后使误差值将至最低,以此提高选取的参数的准确性,从而降低支持向量机模型在实际预测时的预测误差,提高负荷预测的精度。
在其中一个实施例中,在再度训练支持向量机的步骤之后,还包括以下步骤:将相似日的负荷数据输入至再度训练后的支持向量机,获得再度训练后的支持向量机输出的预设天数的负荷数据,并获取再度训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值。
在本实施例中,可以利用历史数据来对再度训练后的支持向量机进行测试,获得再度训练后的支持向量机的预测误差值,以此来判断再度训练后的支持向量机进行负荷预测的精度。
在其中一个实施例中,负荷预测方法主要由两部分构成:
1)基于支持向量机和相似日的预测算法:采用支持向量机作为负荷预测的基础模型,核函数选择为径向基函数,并引入相似日选取作为数据预处理的主要部分;
径向基核函数如下式
其中,xi和xj为任意两个输入量,e为自然底数;参数σ为待定参数,其大小将在参数优化算法中进行调整;
2)对预测算法的参数取值的参数优化算法:为解决相似度阈值以及支持向量机参数的选取问题,选取云遗传算法对其进行优化,云遗传算法优化的效果以支持向量机模型对历史数据的拟合程度为判断依据。
预测算法包括以下步骤:
1、选定支持向量机的参数;
2、把实时数据和历史数据进行比对,在一定相似度阈值下筛选出一组相似日;
其中,实时数据和历史数据的相似度采用以下计算方法:
式中,Sim表示实时数据和历史数据的相似度,x表示实时功率数据构成的行向量,y表示历史功率数据构成的行向量,T为转置符号;
3、把该组相似日的功率数据作为输入,各相似日的后三天(天数可以根据需要设置)的功率数据作为输出,用以训练支持向量机模型;
4、取所得相似日的一部分(比如,取5个相似日)用以计算步骤3得到的支持向量机模型的拟合效果,其中,选取的相似日的功率数据作为输入,然后用模型输出的功率数据与各相似日的后三天的功率数据作误差分析,包括平均绝对百分比误差和均方根相对误差,相关公式如下:
一、平均绝对百分比误差εMAPE计算公式为
其中,n为用作误差分析的负荷值个数,Yi为i时刻预测负荷值,yi为i时刻实际负荷值,yN为被预测对象(风电/光伏/负荷)的额定功率;
二、均方根相对误差εRMSE计算公式为
各物理量含义与平均绝对百分比误差εMAPE计算公式中的相同;
三、总体误差值ε计算公式为
参数优化算法包括以下步骤:
1、取相似度阈值为0.9,支持向量机的参数作为云遗传算法的自变量,各变量取值范围皆为(0,100),云遗传算法的种群大小为100;
2、各个自变量根据各自的取值区间均匀取值,得到各个个体,形成初始种群;
3、把各个个体分别代入预测算法,得到训练好的支持向量机模型及其误差值,该误差值的倒数即为适应度值;
4、基于轮盘赌方法的选择操作,选出一部分个体;
5、对所选出的个体进行基于条件云的交叉操作;
6、对所选出的个体进行基于基本云的变异操作;
7、判断是否满足终止条件,若满足则输出最优个体,否则回到步骤3;
其中,满足以下任一规则即视为满足终止条件:
(1)最大迭代次数为100;
(2)最优个体的适应度在两次连续迭代过程中变化不超过10-8。
在实际进行预测时,将最优个体也即最优各参数输入预测算法,得到再度训练好的支持向量机模型,把实时数据输入,便可得到未来三天的预测功率数据。
条件云和基本云基于正态云模型。正态云模型是一个遵循正态分布规律、具有稳定倾向的随机数集,用期望值Ex、熵En、超熵He三个数值来表征。
参照图3,其为本发明的负荷预测方法对一个水泥厂的发电功率的实际预测效果,虚线为预测功率,实线为实际功率。经误差分析可得,其均方根相对误差为18.21%,平均绝对百分比误差为16.32%,具有较高的计算精度。在对同一预测对象进行多次计算,结果表明,云遗传算法优化的支持向量机模型平均在30次迭代后得到最优模型,收敛速度较快。
本发明涉及三种云发生器:
1)基本云发生器,包括以下步骤:
步骤1:输入Ex、En、He;
步骤2:生成期望值为En、方差为He的正态随机数En’;
步骤3:生成期望值为En’、方差为He的正态随机数x;
步骤4:计算确定度
步骤5:输出云滴(x,μ);
2)X条件云发生器,包括以下步骤:
步骤1:输入Ex、En、He,以及论域上特定的值x0;
步骤2:生成期望值为En、方差为He的正态随机数En’;
步骤3:计算确定度
步骤4:输出云滴(x0,μ);
3)Y条件云发生器,包括以下步骤:
步骤1:输入Ex、En、He,以及特定的确定度μ0;
步骤2:生成期望值为En、方差为He的正态随机数En’;
步骤3:计算x值
步骤5:输出云滴(x,μ0)。
参数优化算法中,基于条件云的交叉操作,对于任何一组用于交叉的两个个体,对个体的每个变量作以下运算:
1、特征数值Ex、En、He计算如下
Ex=Fmax
其中,Fmax和Fmin分别为适应度的最大和最小值;
2、x0取交叉操作的两个个体中较大的适应度,利用X条件云发生器生成一个云滴,并得到相应的确定度μ0;
3、特征数值Ex、En、He计算如下
其中,xf和xm分别为交叉操作的父个体和母个体的某一变量,span为该变量的取值区间长度,Ff和Fm则分别对应它们的适应度,pop为种群大小;
4、结合上述算得的确定度μ0,利用Y条件云发生器生成一对云滴,作为本交叉操作得到的一对儿女的相应变量。
参数优化算法中,基于基本云的变异操作,对于任何一个用于变异操作的个体,对任一变量xp作以下运算:
1、特征数值Ex、En、He计算如下
Ex=xp
其中span(xp)表示变量xp的取值区间长度,pop为种群大小;
2、利用基本云发生器生成一个云滴(x,μ);
3、生成服从[0,1]区间均匀分布的随机数rand,当μ>rand时,更新个体,也即取xp=x。
根据上述支持向量机的参数选取方法,本发明还提供一种支持向量机的参数选取系统,以下就本发明的支持向量机的参数选取系统的实施例进行详细说明。
参见图4所示,为本发明的支持向量机的参数选取系统的结构示意图。该实施例中的支持向量机的参数选取系统包括以下单元:
第一获取单元310,用于获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,支持向量机的各参数对应云遗传算法的各自变量;
生成单元320、用于根据对应的各参数的取值范围对各自变量均匀取值,形成初始种群,其中,初始种群中包括多个个体,个体包括取值后的各自变量;
第一训练单元330,用于针对每个个体,根据取值后的各自变量配置支持向量机,获取实时负荷数据和历史负荷数据以及实时负荷数据和历史负荷数据之间的相似度,从历史负荷数据中筛选出相似日的负荷数据,其中,相似日的负荷数据与实时负荷数据的相似度大于预设值;将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以训练支持向量机;
第二获取单元340,将相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值;
云操作单元350,选取初始种群中的部分个体,以选取的各个体对应的误差值的倒数作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
判断单元360,判断第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各自变量作为支持向量机的各参数;若否,则返回至第一训练单元330。
在其中一个实施例中,云操作单元350用于对于任何一组用于交叉的两个个体,对两个个体中的每个变量作以下运算:
根据Ex=Fmax、和计算第一特征数值,其中,Ex、En和He为第一特征数值,Fmax和Fmin分别为选取的各个体对应的适应度的最大值和最小值;
x0取交叉操作的两个个体对应的适应度中的较大值,将x0以及Ex、En和He输入至X条件云发生器,获得第一确定度;
根据和计算第二特征数值,其中,Ex’、En’和He’为第二特征数值,xf和xm分别为交叉操作的两个个体中对应的同一种变量,span为同一种变量的取值区间长度,Ff和Fm分别为两个个体对应的适应度,pop为种群大小;
将第一确定度以及Ex’、En’和He’输入至Y条件云发生器,获得与同一变量对应的第一子变量和第二子变量,用第一子变量替换xf,用第二子变量替换xm。
在其中一个实施例中,云操作单元用于对于任何一个用于变异操作的个体,对该个体中的任一变量xp作以下运算:
根据和计算第三特征数值,其中,Ex”、En”和He”为第三特征数值,span(xp)表示变量xp的取值区间长度,pop为种群大小;
将Ex”、En”和He”输入至基本云发生器,获得与变量xp对应的第三子变量和第二确定度,生成服从[0,1]区间均匀分布的随机数,当第二确定度大于随机数时,用第三子变量替换变量xp。
在其中一个实施例中,云操作单元350用于基于轮盘赌方法的选择方式选取初始种群中的部分个体。
本发明的支持向量机的参数选取系统与本发明的支持向量机的参数选取方法一一对应,在上述支持向量机的参数选取方法的实施例阐述的技术特征及其有益效果均适用于支持向量机的参数选取系统的实施例中。
根据上述负荷预测方法,本发明还提供一种负荷预测系统,以下就本发明的负荷预测系统的实施例进行详细说明。
参见图5所示,为本发明的负荷预测系统的结构示意图。该实施例中的负荷预测系统包括上述支持向量机的参数选取系统和以下单元:
配置单元410,用于根据上述的支持向量机的参数选取系统选取的支持向量机的参数对支持向量机进行配置;
第二训练单元420,将相似日的负荷数据作为支持向量机的输入,将对应相似日之后的预设天数的负荷数据作为支持向量机的输出,用以再度训练支持向量机;
第三获取单元430,用于将实时负荷数据输入至再度训练后的支持向量机,获得对应实时负荷数据的日期之后的预设天数的负荷数据。
在其中一个实施例中,如图6所示,负荷预测系统还包括第四获取单元440,用于将相似日的负荷数据输入至再度训练后的支持向量机,获得再度训练后的支持向量机输出的预设天数的负荷数据,并获取再度训练后的支持向量机输出的预设天数的负荷数据与相似日之后的预设天数的负荷数据之间的误差值。
本发明的负荷预测系统与本发明的负荷预测方法一一对应,在上述负荷预测方法的实施例阐述的技术特征及其有益效果均适用于负荷预测系统的实施例中。
在本发明中,“第一”、“第二”等序数词只是为了对所涉及的对象进行区分,并不是对对象本身进行限定。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种支持向量机的参数选取方法,其特征在于,包括以下步骤:
获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,所述支持向量机的各参数对应云遗传算法的各自变量;
根据对应的各参数的取值范围对各所述自变量均匀取值,形成初始种群,其中,所述初始种群中包括多个个体,所述个体包括取值后的各所述自变量;
针对每个个体,根据取值后的各所述自变量配置所述支持向量机,获取实时负荷数据和历史负荷数据以及所述实时负荷数据和所述历史负荷数据之间的相似度,从所述历史负荷数据中筛选出相似日的负荷数据,其中,所述相似日的负荷数据与所述实时负荷数据的相似度大于预设值;将所述相似日的负荷数据作为所述支持向量机的输入,将对应所述相似日之后的预设天数的负荷数据作为所述支持向量机的输出,用以训练所述支持向量机;将所述相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取所述训练后的支持向量机输出的预设天数的负荷数据与所述相似日之后的预设天数的负荷数据之间的误差值;
选取所述初始种群中的部分个体,以选取的各个体对应的误差值的倒数作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对所述第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
判断所述第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各所述自变量作为所述支持向量机的各参数;若否,则返回至所述针对每个个体,根据取值后的各所述自变量配置所述支持向量机的步骤。
2.根据权利要求1所述的支持向量机的参数选取方法,其特征在于,所述对选取的各个体进行基于条件云的交叉操作的步骤包括以下步骤:
对于任何一组用于交叉的两个个体,对两个个体中的每个变量作以下运算:
根据Ex=Fmax、和计算第一特征数值,其中,Ex、En和He为第一特征数值,Fmax和Fmin分别为选取的各个体对应的适应度的最大值和最小值;
x0取交叉操作的两个个体对应的适应度中的较大值,将x0以及Ex、En和He输入至X条件云发生器,获得第一确定度;
根据和计算第二特征数值,其中,Ex’、En’和He’为第二特征数值,xf和xm分别为交叉操作的两个个体中对应的同一种变量,span为所述同一种变量的取值区间长度,Ff和Fm分别为两个个体对应的适应度,pop为所述种群大小;
将所述第一确定度以及所述Ex’、En’和He’输入至Y条件云发生器,获得与所述同一变量对应的第一子变量和第二子变量,用所述第一子变量替换xf,用所述第二子变量替换xm。
3.根据权利要求1所述的支持向量机的参数选取方法,其特征在于,所述对所述第一个体组中的各个体进行基于基本云的变异操作的步骤包括以下步骤:
对于任何一个用于变异操作的个体,对该个体中的任一变量xp作以下运算:
根据Ex″=xp、和计算第三特征数值,其中,Ex”、En”和He”为第三特征数值,span(xp)表示变量xp的取值区间长度,pop为所述种群大小;
将所述Ex”、En”和He”输入至基本云发生器,获得与变量xp对应的第三子变量和第二确定度,生成服从[0,1]区间均匀分布的随机数,当所述第二确定度大于所述随机数时,用第三子变量替换所述变量xp。
4.根据权利要求1至3中任意一项所述的支持向量机的参数选取方法,其特征在于,所述选取所述初始种群中的部分个体的步骤包括以下步骤:
基于轮盘赌方法的选择方式选取所述初始种群中的部分个体。
5.一种负荷预测方法,其特征在于,包括以下步骤:
根据如权利要求1至4中任意一项所述的支持向量机的参数选取方法选取所述支持向量机的参数,根据选取的所述支持向量机的参数对所述支持向量机进行配置;
将所述相似日的负荷数据作为所述支持向量机的输入,将对应所述相似日之后的预设天数的负荷数据作为所述支持向量机的输出,用以再度训练所述支持向量机;
将所述实时负荷数据输入至所述再度训练后的支持向量机,获得对应所述实时负荷的日期之后的预设天数的负荷数据。
6.一种支持向量机的参数选取系统,其特征在于,包括以下单元:
第一获取单元,用于获取支持向量机的各参数的取值范围和云遗传算法的各自变量,其中,所述支持向量机的各参数对应云遗传算法的各自变量;
生成单元,用于根据对应的各参数的取值范围对各所述自变量均匀取值,形成初始种群,其中,所述初始种群中包括多个个体,所述个体包括取值后的各所述自变量;
第一训练单元,用于针对每个个体,根据取值后的各所述自变量配置所述支持向量机,获取实时负荷数据和历史负荷数据以及所述实时负荷数据和所述历史负荷数据之间的相似度,从所述历史负荷数据中筛选出相似日的负荷数据,其中,所述相似日的负荷数据与所述实时负荷数据的相似度大于预设值;将所述相似日的负荷数据作为所述支持向量机的输入,将对应所述相似日之后的预设天数的负荷数据作为所述支持向量机的输出,用以训练所述支持向量机;
第二获取单元,将所述相似日的负荷数据输入至训练后的支持向量机,获得训练后的支持向量机输出的预设天数的负荷数据,并获取所述训练后的支持向量机输出的预设天数的负荷数据与所述相似日之后的预设天数的负荷数据之间的误差值;
云操作单元,选取所述初始种群中的部分个体,以选取的各个体对应的误差值的倒数作为适应度值,根据选取的各个体对应的适应度值,对选取的各个体进行基于条件云的交叉操作,获得第一个体组;对所述第一个体组中的各个体进行基于基本云的变异操作,获得第二个体组;
判断单元,判断所述第二个体组中是否有个体满足终止条件,若是,则将该个体中取值后的各所述自变量作为所述支持向量机的各参数;若否,则返回至所述第一训练单元。
7.根据权利要求6所述的支持向量机的参数选取系统,其特征在于,所述云操作单元用于对于任何一组用于交叉的两个个体,对两个个体中的每个变量作以下运算:
根据Ex=Fmax、和计算第一特征数值,其中,Ex、En和He为第一特征数值,Fmax和Fmin分别为选取的各个体对应的适应度的最大值和最小值;
x0取交叉操作的两个个体对应的适应度中的较大值,将x0以及Ex、En和He输入至X条件云发生器,获得第一确定度;
根据和计算第二特征数值,其中,Ex’、En’和He’为第二特征数值,xf和xm分别为交叉操作的两个个体中对应的同一种变量,span为所述同一种变量的取值区间长度,Ff和Fm分别为两个个体对应的适应度,pop为所述种群大小;
将所述第一确定度以及所述Ex’、En’和He’输入至Y条件云发生器,获得与所述同一变量对应的第一子变量和第二子变量,用所述第一子变量替换xf,用所述第二子变量替换xm。
8.根据权利要求6所述的支持向量机的参数选取系统,其特征在于,所述云操作单元用于对于任何一个用于变异操作的个体,对该个体中的任一变量xp作以下运算:
根据和计算第三特征数值,其中,Ex”、En”和He”为第三特征数值,span(xp)表示变量xp的取值区间长度,pop为所述种群大小;
将所述Ex”、En”和He”输入至基本云发生器,获得与变量xp对应的第三子变量和第二确定度,生成服从[0,1]区间均匀分布的随机数,当所述第二确定度大于所述随机数时,用第三子变量替换所述变量xp。
9.根据权利要求6至8中任意一项所述的支持向量机的参数选取系统,其特征在于,所述云操作单元用于基于轮盘赌方法的选择方式选取所述初始种群中的部分个体。
10.一种负荷预测系统,其特征在于,包括如权利要求6至9中任意一项所述的支持向量机的参数选取系统和以下单元:
配置单元,用于根据如6至9中任意一项所述的支持向量机的参数选取系统选取的所述支持向量机的参数对所述支持向量机进行配置;
第二训练单元,将所述相似日的负荷数据作为所述支持向量机的输入,将对应所述相似日之后的预设天数的负荷数据作为所述支持向量机的输出,用以再度训练所述支持向量机;
第三获取单元,用于将所述实时负荷数据输入至所述再度训练后的支持向量机,获得对应所述实时负荷数据的日期之后的预设天数的负荷数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610402774.6A CN106096771A (zh) | 2016-06-07 | 2016-06-07 | 支持向量机的参数选取方法和系统与负荷预测方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610402774.6A CN106096771A (zh) | 2016-06-07 | 2016-06-07 | 支持向量机的参数选取方法和系统与负荷预测方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106096771A true CN106096771A (zh) | 2016-11-09 |
Family
ID=57227558
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610402774.6A Pending CN106096771A (zh) | 2016-06-07 | 2016-06-07 | 支持向量机的参数选取方法和系统与负荷预测方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106096771A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679478A (zh) * | 2017-09-27 | 2018-02-09 | 广州供电局有限公司 | 输电线路空间负荷状态的提取方法和系统 |
CN108171432A (zh) * | 2018-01-04 | 2018-06-15 | 南京大学 | 基于多维云模型-模糊支持向量机的生态风险评价方法 |
CN109118025A (zh) * | 2018-09-25 | 2019-01-01 | 新智数字科技有限公司 | 一种电力系统预测的方法和装置 |
CN109583520A (zh) * | 2018-12-27 | 2019-04-05 | 云南电网有限责任公司玉溪供电局 | 一种云模型与遗传算法优化支持向量机的状态评估方法 |
CN110082136A (zh) * | 2019-04-19 | 2019-08-02 | 电子科技大学 | 基于云遗传算法优化支持向量机的旋转机械故障诊断方法 |
CN110648248A (zh) * | 2019-09-05 | 2020-01-03 | 广东电网有限责任公司 | 一种发电站的控制方法、装置及设备 |
CN113361195A (zh) * | 2021-06-04 | 2021-09-07 | 国网浙江省电力有限公司杭州供电公司 | 一种基于lstm和相似日的配电网线路超短期负荷预测方法 |
-
2016
- 2016-06-07 CN CN201610402774.6A patent/CN106096771A/zh active Pending
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679478A (zh) * | 2017-09-27 | 2018-02-09 | 广州供电局有限公司 | 输电线路空间负荷状态的提取方法和系统 |
CN107679478B (zh) * | 2017-09-27 | 2020-04-14 | 广州供电局有限公司 | 输电线路空间负荷状态的提取方法和系统 |
CN108171432A (zh) * | 2018-01-04 | 2018-06-15 | 南京大学 | 基于多维云模型-模糊支持向量机的生态风险评价方法 |
WO2020063690A1 (zh) * | 2018-09-25 | 2020-04-02 | 新智数字科技有限公司 | 一种电力系统预测的方法和装置 |
CN109118025A (zh) * | 2018-09-25 | 2019-01-01 | 新智数字科技有限公司 | 一种电力系统预测的方法和装置 |
CN109583520A (zh) * | 2018-12-27 | 2019-04-05 | 云南电网有限责任公司玉溪供电局 | 一种云模型与遗传算法优化支持向量机的状态评估方法 |
CN109583520B (zh) * | 2018-12-27 | 2023-04-07 | 云南电网有限责任公司玉溪供电局 | 一种云模型与遗传算法优化支持向量机的状态评估方法 |
CN110082136A (zh) * | 2019-04-19 | 2019-08-02 | 电子科技大学 | 基于云遗传算法优化支持向量机的旋转机械故障诊断方法 |
CN110082136B (zh) * | 2019-04-19 | 2021-03-30 | 电子科技大学 | 基于云遗传算法优化支持向量机的旋转机械故障诊断方法 |
CN110648248A (zh) * | 2019-09-05 | 2020-01-03 | 广东电网有限责任公司 | 一种发电站的控制方法、装置及设备 |
CN110648248B (zh) * | 2019-09-05 | 2023-04-07 | 广东电网有限责任公司 | 一种发电站的控制方法、装置及设备 |
CN113361195A (zh) * | 2021-06-04 | 2021-09-07 | 国网浙江省电力有限公司杭州供电公司 | 一种基于lstm和相似日的配电网线路超短期负荷预测方法 |
CN113361195B (zh) * | 2021-06-04 | 2022-05-13 | 国网浙江省电力有限公司杭州供电公司 | 一种基于lstm和相似日的配电网线路超短期负荷预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106096771A (zh) | 支持向量机的参数选取方法和系统与负荷预测方法和系统 | |
Chang et al. | An improved neural network-based approach for short-term wind speed and power forecast | |
CN105631483B (zh) | 一种短期用电负荷预测方法及装置 | |
Lu et al. | US natural gas consumption prediction using an improved kernel-based nonlinear extension of the Arps decline model | |
Semero et al. | A PSO-ANFIS based hybrid approach for short term PV power prediction in microgrids | |
CN107704953A (zh) | Ewt分位数回归森林的短期风电功率概率密度预测方法 | |
CN109118025A (zh) | 一种电力系统预测的方法和装置 | |
Al-Ghandoor et al. | Electricity consumption in the industrial sector of Jordan: application of multivariate linear regression and adaptive neuro-fuzzy techniques | |
CN106295899B (zh) | 基于遗传算法与支持向量分位数回归的风电功率概率密度预测方法 | |
CN105844371A (zh) | 一种用电客户短期负荷需求预测方法及装置 | |
CN101383023A (zh) | 基于样本动态组织与温度补偿的神经网络短期电力负荷预测 | |
CN102102626A (zh) | 风电场短期功率预测方法 | |
CN103310286A (zh) | 一种具有时间序列特性的产品订单预测方法及装置 | |
CN104699890B (zh) | 一种城网配变重过载短期预警模型的建模方法 | |
CN106529732A (zh) | 基于神经网络与随机前沿分析的碳排放效率预测方法 | |
CN104463356A (zh) | 一种基于多维信息人工神经网络算法的光伏发电功率预测方法 | |
CN108647829A (zh) | 一种基于随机森林的梯级水电站群联合调度规则提取方法 | |
CN108596242A (zh) | 基于小波神经网络和支持向量机的电网气象负荷预测方法 | |
Papageorgiou et al. | Hybrid model for water demand prediction based on fuzzy cognitive maps and artificial neural networks | |
Huang et al. | A hybrid method for one-day ahead hourly forecasting of PV power output | |
CN103530700B (zh) | 城区配电网饱和负荷综合预测方法 | |
Dong et al. | A multifactorial short-term load forecasting model combined with periodic and non-periodic features-a case study of qingdao, China | |
Liu et al. | Exchange rates forecasting with least squares support vector machine | |
Tejaswi et al. | Machine learning algorithms-based solar power forecasting in smart cities | |
Jahan et al. | Intelligent system for power load forecasting in off-grid platform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161109 |
|
RJ01 | Rejection of invention patent application after publication |