基于集合经验模态分解和LSTM的储能装置控制方法与流程
本发明属于电力负荷调度技术领域，涉及一种基于集合经验模态分解和lstm的储能装置控制方法。
背景技术：
作为电力系统经济调度中的一项重要内容，准确地负荷预测能够经济合理的安排电网内部发电机组的启停，保障电网的安稳运行，同时为电网调度计划安排、设备检修和电网改造扩建提供可靠的数据支撑。
近年来，随着用电领域的不断拓宽，用户数量的不断增加，加上新能源在电网中的渗透率越来越高，由于新能源出力具有间歇性和不确定性，使得电网负荷峰谷差逐步增大，供电可靠性降低，因此，在负荷预测的基础上，利用储能装置进行负荷的平移操作，基于目前电力市场分时电价的基础上可以使得峰谷差降低的同时获得经济效益。目前，负荷预测的方法也比较多，常用的有灰色理论法、专家系统法、支持向量机法和时间序列法等。针对负荷预测的不确定性，上述方法往往在预测精度上并不理想，随着神经网络的兴起，其强大的学习能力、自适应能力使得其在模式识别、智能机器人、自动控制等多领域取得了优异的成绩，由于负荷的不确定性，可以采用神经网路进行学习以提高负荷预测的精度。
目前，有人利用小波分析结合神经网络进行时间序列的预测，但是小波分析需要选取合适的母小波以及设置可行的分解层数，针对非线性、非平稳信号的分解的自适应效果差，预测精度仍然有待提高。
因此，研究一种预测精度高的储能装置控制方法具有十分重要的意义。
技术实现要素：
本发明的目的是为了克服上述现有技术中存在的问题，提供一种预测精度高的基于集合经验模态分解和lstm的储能装置控制方法具有十分重要的意义。
为了达到上述目的，本发明采用的技术方案为：
基于集合经验模态分解和lstm的储能装置控制方法，先基于集合经验模态分解和lstm对短期负荷进行预测，然后根据预测结果控制储能装置的充放电；
预测过程为：采集当前时间段以及距离当前时间段最近的n个时间段的历史短期负荷的数据构成一时间序列，对时间序列进行预处理得到多个子序列和残余分量，将各子序列和残余分量分别输入到各自对应的训练好的lstm模型中，一子序列对应一训练好的lstm模型，一残余分量对应一训练好的lstm模型，由训练好的lstm模型输出预测值后，对所有的预测值进行后处理得到下一时间段的未来短期负荷的数据；
lstm模型的训练过程即以连续n+2个时间段中前n+1个时间段的历史短期负荷的数据对应的子序列或残余分量为输入，以后1个时间段的历史短期负荷的数据对应的子序列或残余分量为理论输出，不断调整lstm模型的参数，直到模型损失函数值趋于收敛或达到最大迭代次数为止的过程，其中子序列或残余分量采用与预测过程相同的方法得到；
所述预处理是指先后进行归一化处理和集合经验模态分解，所述后处理是指先后进行重构和反归一化处理；
控制过程为：先读取当前时刻储能装置的荷电状态soc(t)，再根据下一时间段δt的预测结果估算δt内储能装置与外部的交换功率pex，再根据储能装置运行电压u计算δt内与外部的交换电流iex＝pex/u，当交换电流iex>0时，控制储能装置放电；反之，控制储能装置充电。
控制储能装置放电的具体过程为：
计算δt内储能装置最大可放电电流ioutmax，若ioutmax≥iex，标记储能装置能够充电，计算放电后储能装置的荷电状态soc(t+δt)作为参考；若ioutmax＜iex，提示缺电告警信息，计算δt内储能装置最大可放电电流ioutmax采用如下公式：
其中，soc(t-1)为放电初始时刻储能装置的荷电状态，t为储能装置可持续放电的时间，k和n'为取决于不同储能装置自身放电特性的常数，可以通过实验得到；
计算放电后储能装置的荷电状态采用如下公式：
无法持续供应负荷t＜t；
其中，t为在电流iout下的实际放电时间，cn为储能装置额定容量，soc(t)为放电完成后储能装置的荷电状态；
控制储能装置充电的具体过程为：
计算储能装置的可接受电流和储能装置在充电时间t内的充电容量cin，充电电流iin＜i0，储能装置以iin恒电流充电；否则估算储能装置在δt内的最大临界充电电流iinmax，并以iinmax为储能装置充电，同时提示电能过剩，需切除部分电源；
计算储能装置可接受电流采用如下公式：
其中，i0为开始充电时允许的最大初始电流值，η为储能装置的充电接受比，用来表征储能装置的充电接受特性；
计算储能装置在充电时间t内的充电容量cin的公式如下：
cin＝iintiin＜i0，t≥t；
其中，
最大临界充电电流iinmax的计算方法为：令放电仿真过程中的临界时刻t＝δt，即
本发明的基于集合经验模态分解和lstm的储能装置控制方法，集合经验模态分解是一种基于经验模态分解且针对非线性、非平稳信号能自适应分解信号的算法，不需要提前设定母小波，可根据数据自行确定分解层数，且克服了经验模态分解的模态混叠现象；考虑负荷数据在时间上具有相关性，因此可以采用适合于处理和预测时间序列中间间隔和延迟较长的重要事件的lstm模型进行负荷预测，单一利用lstm模型进行负荷预测没有考虑负荷数据的不平稳性，预测精度不够理想，因此可以结合集合经验模态分解和lstm进行负荷序列的预测。如将集合经验模态分解与svr(支持向量机)或深度置信网络等非记忆神经网络结合进行预测，其一般处理的是影响因素与预测结果之间的映射关系，在处理时序问题上时没有考虑历史数据之间的关系，而负荷相对具有时间上的联系，相对于集合经验模态分解结合记忆网络(即lstm模型)，其在时序预测上的准确率较低。将经验模态分解与神经网络结合进行预测，有时候会考虑将分解后的各子序列进行相关系数比较，选取有用的imf分量进行重构来得到一个与原始数据近似的新数据输入到单一网络模型进行训练预测，此方法虽然对预测的总体变化趋势不会产生太大影响，但是没有完全利用分解后的子序列，被抛弃的imf分量往往体现出随机性因素导致的数据波动，对于负荷预测，这种波动具有很好的参考价值，而且这种重构后的数据不一定具有平稳性，本发明通过对分解后的各子序列分别采用神经网络进行预测，再对各子序列预测结果进行重构可以提高预测的精度。
作为优选的技术方案：
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，n＝2。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，所述归一化处理的公式如下：
式中，xmax为归一化处理数据的最大值，xmin为归一化处理数据的最小值，x为归一化前的值，x′是归一化后的值。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，所述预处理还包括在归一化处理之前的异常值检测及修正处理；
异常值检测的过程为：将所有数据按大小顺序排列，记上四分位值为q1，即所有数据中只有1/4的数据大于q1，下四分位值为q2，即所有数据中只有1/4的数据小于q2，上界为(q1+1.5(q1-q2))，下界为(q2-1.5(q1-q2))，介于上下界之间的都是正常观测值，否则为异常值；
异常值修正的方法为：去除异常值后，将当天已有的负荷数据进行三次样条插值，三次样条插值(简称spline插值)是通过一系列形值点的一条光滑曲线，数学上通过求解三弯矩方程组得出曲线函数组的过程，用与异常值同时间段的三次样条插值上的负荷数据替代异常值。
三次样条插值建模：
设f(x)是区间[a,b]上的一个连续可微函数，在区间[a,b]上给定一组基点，由于负荷数据的因变量是日期，故可以用其在一天中的第几个时间段代替，当天去除异常点后的负荷数据个数为(n+1)：
a＝x0＜x1＜x2＜···＜xn＝b；
设函数s(x)满足条件：
(1)s(x)在每个子区间[xi-1,xi]上有表达式：si(x)＝aix3+bix2+cix+di；
(2)s(x)在区间[a,b]上有2阶连续导数；
根据以下已知条件取求解每个子区间的si(x)：
si(x)求出之后，可将异常值对应的时间段的
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，所述集合经验模态分解的步骤如下：
(1)将正态分布白噪声序列n(t)加到时间序列x(t)中；
x'(t)＝x(t)+n(t)；
(2)将添加正态分布白噪声序列后的时间序列x'(t)作为一个整体，然后进行经验模态分解，得到各imf分量；
(3)重复步骤(1)和(2)100次，每次加入新的正态分布白噪声序列；
(4)将每次得到的imf分量做集成平均处理。
本发明的集合经验模态分解法是在经验模态分解(emd)的基础上发展而来的，主要是将时间序列数据加上白噪声序列后整体进行经验模态分解处理，克服了后者在分解过程中出现混叠的现象，并且降低了后续过程中重构信号的偏差，前者能对信号降噪的同时更准确地复原信号的有用信息。
本发明的集合经验模态分解法将非平稳的负荷数据分解成频率不同的子序列以及残余分量，而且这些序列具有平稳性，平稳性意味着时序数据的均值、方差、协方差是与时间无关的常数，符合统计中的一致性，以便于更好的预测。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，本发明的lstm模型是一种时间递归神经网络，适合于处理和预测时间序列中间隔和延迟相较长的重要事件，其工作核心是通过遗忘门、输入门和输出门三个门来控制历史信息和当前信息的记忆和遗忘程度，lstm模型的遗忘门模型函数的表达式如下：
ft＝σ(wf[ht-1,xt]+bf)；
输入门模型函数的表达式如下：
it＝σ(wi·[ht-1,xt]+bi)；
输出门模型函数的表达式如下：
ot＝σ(wo·[ht-1,xt]+bo)；
最终的时间序列上的输出量的表达式如下：
ht＝ot*tanh(ct)；
上述各表达式中，wf和bf分别为遗忘门的权重系数矩阵和偏置项；wi和bi分别为输入门的权重系数矩阵和偏置项；wc和bc分别为细胞状态更新时的权重系数矩阵和偏置项；wo和bo分别为输出门的权重系数矩阵和偏置项；σ和tanh分别为sigmoid和双曲正切激活函数，ht-1表示的是上一个细胞的输出，xt表示当前细胞的输入，ft、it、ct、ot及ht分别为遗忘门的结果、输入门的结果、更新状态后的结果、输出门的结果和时间序列上的输出结果。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，lstm模型的训练步骤如下：
(1)确定lstm模型的结构(即模型的lstm层数及每层的神经元个数)、激活函数和损失函数，并初始化模型参数，令迭代次数s＝1，其中损失函数mae的表达式如下：
式中，yi为子序列或残余分量对应的预测值的期望值，
sigmoid激活函数表达式如下：
双曲正切激活函数表达式如下：
(2)前向计算每个神经元的输出值，即计算ft,it,ct,ot,ht五个向量的值；
(3)反向计算每个神经元的误差项值，lstm误差项的反向传播包括两个方向：一个是沿时间的反向传播，即从当前时间段t开始，计算每个时间段的误差项；另一个是将误差项向上一层传播；
(4)判断损失函数值是否收敛，同时判断是否达到最大迭代次数，如果有一个是，则停止更新；反之，则令迭代次数s＝s+1后，返回步骤(2)。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，重构采用基于遗传算法的bp算法，重构可以采用bp网络进行非线性拟合处理，但传统bp网络学习速度慢，容易陷入局部极小值，因此本发明采用遗传算法求解bp网络模型参数去改善bp网络模型收敛速度慢及容易陷入局部最优的缺点，步骤如下：
(1)bp网络建模，采用只含一层隐含层的bp网络，设置输入层神经元的个数为k，隐含层神经元的个数为m，输出层神经元的个数为1；
(2)利用遗传算法优化bp网络的权值和阈值，具体为：
(2.1)编码，对bp网络的权值和阈值进行二进制编码，bp网络中输入层有k个输入节点，隐含层有m个节点，输出层有1个输出节点，共有k*m+m*1个权值，m+1个阈值，所以编码长度l＝k*m+m+m+1，设置最大进化代数为r，以一组权值和阈值作为一个个体，随机生成m个个体作为初始群体p(0)，迭代次数λ＝0；
(2.2)计算适应度函数值，适应度函数的表达式如下：
式中，y′j为子序列或残余分量对应的预测值输入至bp网络后bp网络的输出值，yj为负荷数据归一化后的实际值，q为参与模型训练数据的个数；
(2.3)选择，根据个体的适应度函数值所量度的优劣程度决定它在下一代淘汰还是被遗传，用赌轮选择机制，令∑fj表示群体的适应度值之总和，fj表示群体中第j个个体的适应度值，它产生后代的能力刚好为其适应度值所占份额fj/∑fj；
(2.4)交叉，将交叉算子作用于群体；
(2.5)变异，将变异算子作用于群体，对群体中个体串的某些基因座上的基因值做变动；
(2.6)判断，若迭代次数λ＝r，则将进化过程中所得到的具有最大适应度值的个体进行解码得到bp网络的最优权值和阈值，并输出，终止计算；反之，则λ＝λ+1，返回步骤(2.2)；
(3)bp网络对lstm模型的所有预测值进行非线性重构。
传统的重构方法是将各预测结果简单相加，不能保证预测误差整体偏小，因此本发明采用了基于遗传算法的bp算法对各个预测结果进行重构，使得模型更加准确。此方法能够满足家庭负荷预测精度的需求，且实现方便。
如上所述的基于集合经验模态分解和lstm的储能装置控制方法，其特征在于，所述反归一化处理的公式如下：
z＝z'(xmax-xmin)+xmin；
式中，z′为反归一化处理前的值即重构后的数据值，z为反归一化处理后的值即为最终的负荷预测值。
有益效果：
(1)本发明的基于集合经验模态分解和lstm的储能装置控制方法，对于波动性、随机性以及不确定性较强的用电负荷有着较好的预测精度，能够更加合理地控制储能装置充放电；
(2)本发明的基于集合经验模态分解和lstm的储能装置控制方法，采用集合经验模态分解法将非平稳的负荷数据分解成频率不同的子序列以及残余分量，子序列以及残余分量的平稳性好，便于更好的预测；
(3)本发明的基于集合经验模态分解和lstm的储能装置控制方法，集合经验模态分解法是将时间序列数据加上白噪声序列后整体进行经验模态分解处理，克服了后者在分解过程中出现混叠的现象，并且降低了后续过程中重构信号的偏差，前者能对信号降噪的同时更准确地复原信号的有用信息；
(4)本发明的基于集合经验模态分解和lstm的储能装置控制方法，基于遗传算法的bp算法对各个预测结果进行重构而非传统重构方法中的将各预测结果简单相加，使得模型重构后的预测结果更加准确。
附图说明
图1为本发明的基于集合经验模态分解和lstm的储能装置控制方法的流程图；
图2为本发明对时间序列进行异常值检测的检测示意图；
图3为本发明对时间序列进行集合经验模态分解的流程图；
图4为本发明lstm模型的拓扑结构图；
图5为本发明优化bp网络得到基于遗传算法的bp网络的流程图。
具体实施方式
下面结合具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
基于集合经验模态分解和lstm的储能装置控制方法，如图1所示，步骤如下：
(1)训练lstm模型；
(1.1)采集连续n+2个时间段中前n+1个时间段的历史短期负荷的数据构成一时间序列，对时间序列进行预处理得到多个子序列和残余分量，n＝2，具体为：
(1.1.1)对时间序列进行异常值检测及修正处理，如图2所示；
异常值检测的过程为：将所有数据按大小顺序排列，记上四分位值为q1，即所有数据中只有1/4的数据大于q1，下四分位值为q2，即所有数据中只有1/4的数据小于q2，上界为(q1+1.5(q1-q2))，下界为(q2-1.5(q1-q2))，介于上下界之间的都是正常观测值，否则为异常值；
异常值修正的方法为：去除异常值后，将当天已有的负荷数据进行三次样条插值，三次样条插值(简称spline插值)是通过一系列形值点的一条光滑曲线，数学上通过求解三弯矩方程组得出曲线函数组的过程，用与异常值同时间段的三次样条插值上的负荷数据替代异常值；
三次样条插值建模：
设f(x)是区间[a,b]上的一个连续可微函数，在区间[a,b]上给定一组基点，由于负荷数据的因变量是日期，固可以用其在一天中的第几个时间段代替，当天去除异常点后的负荷数据个数为(n+1)：
a＝x0＜x1＜x2＜···＜xn＝b；
设函数s(x)满足条件：
a.s(x)在每个子区间[xi-1,xi]上有表达式：si(x)＝aix3+bix2+cix+di；
b.s(x)在区间[a,b]上有2阶连续导数；
根据以下已知条件取求解每个子区间的si(x)：
si(x)求出之后，可将异常值对应的时间段
(1.1.2)对修正后的时间序列进行归一化处理，归一化处理的公式如下：
式中，xmax为归一化处理数据的最大值，xmin为归一化处理数据的最小值，x为归一化前的值，x′是归一化后的值；
(1.1.3)对归一化处理处理后的数据进行集合经验模态分解(即eemd分解)，如图3所示，集合经验模态分解的步骤如下：
(1.1.3.1)将正态分布白噪声序列n(t)加到时间序列x(t)中，正态分布白噪声序列的标准差设置一般为0.01～0.4；
x'(t)＝x(t)+n(t)；
(1.1.3.2)将添加正态分布白噪声序列后的时间序列x'(t)作为一个整体，然后进行经验模态分解，得到各imf分量(imfs)；
(1.1.3.3)重复步骤(1.1.3.1)和(1.1.3.2)100次，每次加入新的正态分布白噪声序列；
(1.1.3.4)将每次得到的imf分量做集成平均处理；
(1.2)以连续n+2个时间段中前n+1个时间段的历史短期负荷的数据对应的子序列或残余分量为输入，以后1个时间段的历史短期负荷的数据对应的子序列或残余分量为理论输出，不断调整lstm模型的参数，直到模型损失函数值趋于收敛或达到最大迭代次数为止的过程，具体为：
(1.2.1)确定lstm模型的结构(即模型的lstm层数及每层的神经元个数)、激活函数和损失函数，并初始化模型参数，lstm模型的拓扑结构图如图4所示，lstm遗忘门模型函数的表达式如下：
ft＝σ(wf·[ht-1,xt]+bf)；
输入门模型函数的表达式如下：
it＝σ(wi·[ht-1,xt]+bi)；
输出门模型函数的表达式如下：
ot＝σ(wo·[ht-1,xt]+bo)；
令迭代次数s＝1，其中lstm模型的损失函数mae的表达式如下：
式中，yi为子序列或残余分量对应的预测值的期望值，
sigmoid激活函数表达式如下：
双曲正切激活函数表达式如下：
(1.2.2)前向计算每个神经元的输出值，即计算ft,it,ct,ot,ht五个向量的值；
(1.2.3)反向计算每个神经元的误差项值，lstm误差项的反向传播包括两个方向：一个是沿时间的反向传播，即从当前时间段t开始，计算每个时间段的误差项；另一个是将误差项向上一层传播；
(1.2.4)判断损失函数值是否收敛，同时判断是否达到最大迭代次数，如果有一个是，则停止更新；反之，则令迭代次数s＝s+1后，返回步骤(1.2.2)；
最终的时间序列上的输出量的表达式如下：
ht＝ot*tanh(ct)；
上述各表达式中，wf和bf分别为遗忘门的权重系数矩阵和偏置项；wi和bi分别为输入门的权重系数矩阵和偏置项；wc和bc分别为细胞状态更新时的权重系数矩阵和偏置项；wo和bo分别为输出门的权重系数矩阵和偏置项；σ和tanh分别为sigmoid和双曲正切激活函数，ht-1表示的是上一个细胞的输出，xt表示当前细胞的输入，ft、it、ct、ot及ht分别为遗忘门的结果、输入门的结果、更新状态后的结果、输出门的结果和时间序列上的输出结果；
(2)采集当前时间段以及距离当前时间段最近的n个时间段的历史短期负荷的数据构成一时间序列，对时间序列进行预处理得到多个子序列和残余分量，预处理步骤与步骤(1.1)相同；
(3)将预处理得到的各子序列和残余分量分别输入到各自对应的训练好的lstm模型中，一子序列对应一训练好的lstm模型，一残余分量对应一训练好的lstm模型，由训练好的lstm模型输出预测值；
(4)对预测值进行重构，重构采用基于遗传算法的bp算法，步骤如下：
s1：首先优化bp网络得到基于遗传算法的bp算法，具体流程如图5所示，
(4.1)bp网络建模，采用只含一层隐含层的bp网络，设置输入层神经元的个数为k，隐含层神经元的个数为m，输出层神经元的个数为1；
(4.2)利用遗传算法优化bp网络的权值和阈值，具体为：
(4.2.1)编码，对bp网络的权值和阈值进行二进制编码，bp网络中输入层有k个输入节点，隐含层有m个节点，输出层有1个输出节点，共有k*m+m*1个权值，m+1个阈值，所以编码长度l＝k*m+m+m+1，设置最大进化代数为r，随机生成m个个体作为初始群体p(0)，迭代次数λ＝0；
(4.2.2)计算适应度函数值，适应度函数的表达式如下：
式中，y′j为子序列或残余分量对应的预测值输入至bp网络后bp网络的输出值，yj为负荷数据归一化后的实际值，q为参与模型训练数据的个数；
(4.2.3)选择，根据个体的适应度函数值所量度的优劣程度决定它在下一代淘汰还是被遗传，用赌轮选择机制，令∑fj表示群体的适应度值之总和，fj表示群体中第j个个体的适应度值，它产生后代的能力刚好为其适应度值所占份额fj/∑fj；
(4.2.4)交叉，将交叉算子作用于群体；
(4.2.5)变异，将变异算子作用于群体，对群体中个体串的某些基因座上的基因值做变动；
(4.2.6)判断，若迭代次数λ＝r，则将进化过程中所得到的具有最大适应度值的个体进行解码得到bp网络的最优权值和阈值，并输出，终止计算；反之，则λ＝λ+1，返回步骤(4.2.2)；
s2：优化后的bp网络对lstm模型的所有预测值进行非线性重构；
(5)对重构后的数据进行反归一化处理，反归一化处理的公式如下：
z＝z'(xmax-xmin)+xmin；
式中，z′为反归一化处理前的值即重构后的数据值，z为反归一化处理后的值即为最终的负荷预测值；
(6)根据预测结果调整储能装置(如蓄电池)的运行状态，具体为：
先读取当前时刻储能装置的荷电状态soc(t)，再根据下一时间段δt的预测结果估算δt内储能装置与外部的交换功率pex，再根据储能装置运行电压u计算δt内与外部的交换电流iex＝pex/u，当交换电流iex>0时，控制储能装置放电；反之，控制储能装置充电，控制储能装置放电的具体过程为：
计算δt内储能装置最大可放电电流ioutmax，若ioutmax≥iex，标记储能装置能够充电，计算放电后储能装置的荷电状态soc(t+δt)作为参考；若ioutmax＜iex，提示缺电告警信息，计算δt内储能装置最大可放电电流ioutmax采用如下公式：
其中，soc(t-1)为放电初始时刻储能装置的荷电状态，t为储能装置可持续放电的时间，k和n'为取决于不同储能装置自身放电特性的常数，可以通过实验得到；
计算放电后储能装置的荷电状态采用如下公式：
无法持续供应负荷t＜t；
其中，t为在电流iout下的实际放电时间，cn为储能装置额定容量，soc(t)为放电完成后储能装置的荷电状态；
控制储能装置充电的具体过程为：
计算储能装置的可接受电流和储能装置在充电时间t内的充电容量cin，充电电流iin＜i0，储能装置以iin恒电流充电；否则估算储能装置在δt内的最大临界充电电流iinmax，并以iinmax为储能装置充电，同时提示电能过剩，需切除部分电源；
计算储能装置可接受电流采用如下公式：
其中，i0为开始充电时允许的最大初始电流值，η为储能装置的充电接受比，用来表征储能装置的充电接受特性；
计算储能装置在充电时间t内的充电容量cin的公式如下：
cin＝iintiin＜i0，t≥t；
其中，
最大临界充电电流iinmax的计算方法为：令放电仿真过程中的临界时刻t＝δt，即
经验证，本发明通过对不平稳的负荷数据进行集合经验模态分解处理转化为平稳性序列后进行长短期记忆网路(即lstm模型)训练预测并结合改进的bp神经网路(基于遗传算法的bp算法)进行重构，可以提高电力负荷预测的精度，而且泛化能力较强，且实现过程相对简单，具有很好的实用价值。
- 还没有人留言评论。精彩留言会获得点赞！