CN105205502A - 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 - Google Patents
一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 Download PDFInfo
- Publication number
- CN105205502A CN105205502A CN201510733268.0A CN201510733268A CN105205502A CN 105205502 A CN105205502 A CN 105205502A CN 201510733268 A CN201510733268 A CN 201510733268A CN 105205502 A CN105205502 A CN 105205502A
- Authority
- CN
- China
- Prior art keywords
- load
- matrix
- probability
- sigma
- markov
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,包括:找出电压降落的时间点,在电压降落的时间点对应的扰动时刻,进行负荷动特性提取与分类;判断负荷类别间的变化是否具有马氏性,将所有数据按时间平均分段,对每段数据基于最大似然的思想建立马氏链的概率转移矩阵;判断数字特征是否改变,如果未改变,则转入步骤五,否则依据矩阵相应的数字特征对时间段的负荷数据进行聚类,对有改变了数据特征的各时间段负荷数据求概率转移矩阵;利用马尔可夫蒙特卡罗模拟,描述负荷变化情况;利用隐马尔科夫模型处理反映负荷类别转化的序列。本方法对马尔柯夫链蒙特卡罗模拟进行了改进,有效降低了矩阵迭代后陷入稳态的可能性。
Description
技术领域
本发明涉及一种基于马尔柯夫蒙特卡罗的负荷特性分类综合方法。
背景技术
负荷建模是电力系统建模中的一个基础和关键问题。建立能够反映负荷特性准确的负荷模型一直是挑战性的难点问题。负荷建模的最大困难在于负荷的随机性、时变性,它包括负荷大小的改变和负荷组成成分的改变,虽然如此,但是负荷特性存在一定的规律性。
为了在把握规律的基础上,解决负荷构成的随机性、时变性问题,对负荷动特性进行分类与综合。负荷动特性分类与综合是将不同时间采集的同一变电站的动态负荷特性数据中负荷构成成分相似的归为一类,并通过综合的方法建立每一类的负荷模型;且综合得到的模型能反映该类负荷的动特性特征,同时具有对不同强度扰动的内插外推能力。
发明内容
为解决现有技术存在的不足,本发明公开了一种针对电力负荷随机性、时变性的分类特点,提供了一种精确度高的针对电力负荷随机性、时变性的分类综合方法,克服了现有负荷建模方法无法考虑负荷的随机性、时变性问题。
为实现上述目的,本发明的具体方案如下:
一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,包括以下步骤:
步骤一:找出电压降落的时间点,在电压降落的时间点对应的扰动时刻,进行负荷动特性提取与分类;
步骤二:判断负荷类别间的变化是否具有马氏性,如果是,则进行步骤三,否则,结束;
步骤三:将所有数据按时间平均分段,对每段数据基于最大似然的思想建立马氏链的概率转移矩阵;
步骤四:判断数字特征即概率转移矩阵的稳态分布作为一个向量是否改变,当向量的中心距超过规定门阈值时,认为数字特征已经改变,依据概率转移矩阵相应的数字特征对时间段的负荷数据进行聚类,对有改变了数据特征的各时间段负荷数据求概率转移矩阵,然后转入步骤五,如果数字特征未改变,则直接转入步骤五;
步骤五:利用概率转移矩阵生成累积状态转移率矩阵,用马尔柯夫蒙特卡罗仿真(MCMC)对负荷状态在整个时段内的变化进行内插与外推,得到完整的负荷数据序列;
步骤六:利用隐马尔科夫模型处理步骤五的结果即反映负荷类别转化的序列,得到更为精确的负荷类别变化数据。
进一步的,所述步骤四中,数字特征取为各类别的稳态分布:
式中,P为概率转移矩阵,pij为概率转移矩阵的第i行第j列的元素,函数的返回值为矩阵的任意行矢量,π*表示第一步求取的最优稳态分布,求J的最小值。
对于这样一个非线性优化问题,应用遗传算法,容易得到近似极值,其中稳态分布的存在性可转化为最小化Pm的列矢量与列矢量均值的差的二范数。
进一步的,所述步骤四中,依据矩阵相应的数字特征对时间段的负荷数据进行聚类:设x为聚类向量,将目标函数设为样本到所属类别中心的距离,采用模糊C均值聚类(fuzzyc-meansalgorithm,FCM),用迭代的方法求取最小目标函数。
FCM算法把样本空间x={x1,x2,…,xn}分为c类(2≤c≤n),任一样本点xi不会严格划分为某一类。FCM用模糊划分,使用隶属度uij(0≤uij≤1)来确定样本点xi属于第j(0≤j≤c)类的程度。如样本空间x的一个模糊子集所对应的隶属度矩阵是一个模糊隶属度矩阵,用U={Uij}表示。隶属度矩阵U具有如下性质:
式中,uij为隶属度。
FCM算法就是在式(1)的约束条件下使目标函数J最小化,即:
式中,m∈[1,∞)是模糊加权系数;cj是c类中第j类的聚类中心;是xi到cj的欧拉距离。
其步骤如下:
(1)给定聚类数c,模糊加权系数m,迭代停止阈值ε;设迭代次数k=0,最大迭代次数kmax;用值在[0,1]间的随机数初始化隶属矩阵U,使其满足式(1)中的约束条件;
(2)由式(2)计算聚类中心
式中,uij为隶属度,xi为样本;
(3)由聚类中心cj更新隶属度矩阵U(k+1),即
式中,dij为样本和各类中心的距离,上标(k+1)表示第k次迭代的相应数值;
(4)给定收敛的判别精度ε>0,若||U(k+1)-U(k)||≤ε,停止迭代;否则置k=k+1,并返回步骤(2);
(5)得到x的一个最优模糊C划分隶属度矩阵U={Uij}和聚类中心c={cij}。
进一步的,所述步骤五中,在进行内插和外推时,用概率转移矩阵生成累积状态转移率矩阵:
Pcum中元素的取值如下:
其中,pij为概率转移矩阵P的第i行第j列的元素。换言之,经过本步骤的处理,累积状态转移率矩阵的Pcum,il为P的所有第i行前j列元素之和。
进一步的,所述步骤五中,利用马尔柯夫蒙特卡罗仿真MCMC法生成给定时间长度负荷类别变化时间序列(给定时间上限)步骤如下:
1)随机产生一个在区间[0,1]内的数u;
2)将u与Pcum的第p行元素(元素行号与当前所处状态号相同)进行比较,假设u落在累积状态转移率矩阵Pcum的Pcum,pq到Pcum,p(q+1)这个范围内,则认为负荷状态的序列的下一时刻状态为q;
3)如果生成的序列已经满足时间长度的要求,则转步骤1)继续,否则,令当前状态变为q,转步骤1)步继续,直到达到时间要求。
进一步的,所述步骤五中,上述步骤在实现时,每次模拟只对相邻故障时间间隔内的负荷状态进行模拟。
进一步的,所述步骤五中,因为故障时间点随机出现,各时间点之间的时间间隔不统一,相邻两个时间点之间可能还有其他的负荷状态,所以要统一步长为常数△t,累积状态转移率矩阵的处理相应也要进行改进:
转移概率矩阵中,P12表示经过一个步长的时间后,从1状态转化为2状态的概率,也可以表示,现在为状态2,一个步长的时间之前的状态为1的概率,所以
式中,P12,k为已知△t时刻,负荷状态为1,则第(k+1)△t时刻,负荷处于第2种状态的概率。
P12,kt=P21,k为已知第(k+1)△t时刻,负荷状态为1,则△t时刻,负荷处于第2种状态的概率。
问题可以细化为假设已知△t、10△t两个时刻的负荷状态分别为a与b,已知马尔柯夫链P,采用蒙特卡洛仿真内插确定[△t,10△t]时间段内各时间点负荷状态。以确定5△t时刻的负荷情况为例,P4、(P5)T分别为首、末两端与5△t时刻负荷状态的转换概率分布。
根据马尔可夫性,△t时刻的情况与5△t时刻有关,而5△t时刻的情况与△t时刻有关,即一个马尔柯夫过程的正向和反向都具有马尔柯夫性,所以,5△t时刻的情况应该结合△t、10△t两个时刻的负荷状态进行确定。
对P4进行处理,令
Pcum1中最后一列元素为1,元素的取值如下:
对(P5)T进行处理,令
Pcum2中元素的取值如下:
再将Pcum2各行归一化。
从△t、10△t两个时刻对5△t时刻的情况进行蒙特卡罗仿真,将Pcum=(Pcum1+Pcum2)/2,然后利用前文所述的利用MCMC法生成给定时间长度负荷类别变化时间序列的处理方法进行仿真。
进一步的,所述步骤六中,隐马尔科夫模型HMM是一个双重随机过程,它有两个组成部分:
a马尔可夫链----描述类别的转移,用转移概率描述;
b一般随机过程----描述状态与观察序列间的关系,用观察值概率描述。
两个随机过程的关系如图1所示,显而易见,两者之间是嵌套的关系。
一个隐马尔科夫模型HMM常用五元数组来描述:λ=(Q,O,A,B,π)。其中:状态的有限集合Q={q1,q2,…,qn,},n为马氏链的状态数,即对应本文中的各种负荷类别;观察值的有限集合O={o1,o2,…,om},m为每种类别对应的可能的观察值数目;状态转移概率矩阵A;观察值概率矩阵B;初始概率分布向量π。
用Baum-Welch算法,基于HMM模型识别负荷特性的变化。步骤如下:
1)建立初始模型即待训练模型λ0;
2)基于λ0以及观察值序列O,训练新模型λ;
3)如果logP(O|λ)-logP(O|λ0)<D,说明训练已经达到预期效果,算法结束,式中,P(O|λ)表示在实际模型为λ的情况下,观察序列为O的概率,D为设定门阈值;
4)否则,令λ0=λ,继续步骤2)工作。
训练过程前,需要定义两个辅助变量。第一个辅助变量为ξt(i,j),表示给定模型λ和观察序列条件下,从i到j的转移概率:
式中,αt(i)=P(o1,o2,…,ot,xt=i),aij=P(xt+1=j|xt=i),bj(xt+1)=P(ot+1|xt+1=j),βt+1(j)=P(ot+1,ot+2,…,oT|qt=i,λ)。
第二个辅助变量为为t时刻处于状态Si的概率,为整个过程中从状态Si转出次数的预期,为从Si跳转到Sj次数的预期。
训练过程:为当t=1时,处于qi的概率。
分别为状态转移概率矩阵A;观察值概率矩阵B中的元素的估计值。当收敛或达到预置次数时,即可得到状态转移概率矩阵A作为本发明的最终结果。
本发明的有益效果:
1.马氏链的概率转移矩阵更能体现负荷特性各类别之间的转化规律。
2.可精确确定负荷类别的故障时间点是随机出现的,各时间点之间的时间间隔不统一,相邻两个时间点之间可能还有其他的负荷类别,基于最大似然思想建立马氏链的概率转移矩阵,体现负荷类别转化时的真实情况,误差较小。
3.数字特征是对马氏链特性的体现,当各数字特征都明显改变的时候,反映的正是马氏链的变化。通过数字特征对时间段进行聚类,针对不同类的时段分别建立概率矩阵,有效减小了误差。
4.本方法对马尔柯夫链蒙特卡罗模拟进行了改进,有效降低了矩阵迭代后陷入稳态的可能性。
5.本发明之电力负荷时变性分类综合方法适用于针对同一变电站在一段时间内负荷变化情况的描述,能够满足误差要求,极具工程实用价值。
附图说明
图1两个随机过程的关系示意图。
图2本发明针对电力负荷时变性的分类综合方法流程方框图。
具体实施方式:
下面结合附图对本发明进行详细说明:
如图2所示,执行步骤01,开始;
其次,执行步骤02,计算统计量此统计量服从自由度为(n-1)2的χ2分布。
设mij表示X(0)(t)从状态i经过一步转移到状态j的频数,并将转移频数矩阵的各列之和除以各行各列的总和,所得到的值记为P'j:
记
则统计量
给定置信度α,查表可得到值,如果则认为该序列符合马氏性,否则不是马氏链。
再次,执行步骤03,采取基于PSO的非线性规划求取概率转移矩阵。
本文采取基于粒子群模拟退火算法(particleswarmoptimizationsimulatedannealingalgorithm,PSOSA)的非线性规划求取概率转移矩阵。其步骤如下
1)约束条件为当i与k都是1到N之间的整数时,
其中,Pik为转移矩阵各元素。
2)确定适应度函数。适应度函数的输入量是马尔可夫矩阵,适应度函数的输出量是概率之和。某相邻两个故障时刻的负荷类别分别为甲乙,若按理论计算,即计及步长地按马氏链的规律转化,得到相应矩阵,将此矩阵中表明从甲到乙的概率计入计数器,以此类推,将各次转化的概率叠加,即得到适应度函数的输出量,问题也转化为求此适应度函数的最大值。
其中,m为检测时间点数,ti为第i个检测点,si为第i个检测点处的负荷类别,为中第si行第si+1列元素。
3)初始化粒子群,对于每一个粒子计算目标函数值,确定自身最优位置pi和全体最优位置pg。
4)执行一次SA搜索,在温度T下从初始位置x0开始,对解产生随机扰动生成新解并用Metropolis准则判定是否接受,重复进行Lk(马尔科夫链长度)次直到得到一个新解x'.执行退火操作降低温度,Tk+1=K*Tk。
5)执行一次PSO搜索,粒子群向个体最优与全局最优移动,并更新个体最优与全局最优。
其中,w为惯性因子,x为当前粒子所在位置,v为例子飞翔速度。
6)若f(x')<f(pg),随机在全体粒子中选取一个粒子i令xi=x',返回5);若f(x')>f(pg),则令x0=pg,返回4);若f(x')=f(pg)且模拟退火算法和粒子群算法都满足终止条件,结束计算。
最终得到x的最优值,即概率转移矩阵
然后,执行步骤04,求各类别的稳态分布
式中,P为概率转移矩阵,pij为概率转移矩阵的第i行第j列的元素,函数的返回值为矩阵的任意行矢量,π*表示第一步求取的最优稳态分布。求J的最小值。
对于这样一个非线性优化问题,应用遗传算法,容易得到近似极值,其中稳态分布的存在性可转化为最小化Pm的列矢量与列矢量均值的差的二范数。
依据矩阵相应的数字特征对时间段的负荷数据进行聚类:设x为聚类向量,将目标函数设为样本到所属类别中心的距离,采用模糊C均值聚类(fuzzyc-meansalgorithm,FCM),用迭代的方法求取最小目标函数。
FCM算法把样本空间x={x1,x2,…,xn}分为c类(2≤c≤n),任一样本点xi不会严格划分为某一类。FCM用模糊划分,使用隶属度uij(0≤uij≤1)来确定样本点xi属于第j(0≤j≤c)类的程度。如样本空间x的一个模糊子集所对应的隶属度矩阵是一个模糊隶属度矩阵,用U={Uij}表示。隶属度矩阵U具有如下性质:
式中,uij为隶属度。
FCM算法就是在式(1)的约束条件下使目标函数J最小化,即:
式中,m∈[1,∞)是模糊加权系数;cj是c类中第j类的聚类中心;是xi到cj的欧拉距离。
其步骤如下:
(1)给定聚类数c,模糊加权系数m,迭代停止阈值ε;设迭代次数k=0,最大迭代次数kmax;用值在[0,1]间的随机数初始化隶属矩阵U,使其满足式(1)中的约束条件;
(2)由式(2)计算聚类中心
式中,uij为隶属度,xi为样本。
(3)由聚类中心cj更新隶属度矩阵U(k+1),即
式中,dij为样本和各类中心的距离,上标(k+1)表示第k次迭代的相应数值,。
(4)给定收敛的判别精度ε>0,若||U(k+1)-U(k)||≤ε,停止迭代;否则置k=k+1,并返回步骤(2);
(5)得到x的一个最优模糊C划分隶属度矩阵U={Uij}和聚类中心c={cij}。
接着,执行步骤05,对系统进行内插和外推。
用概率转移矩阵生成累积状态转移率矩阵。
Pcum中元素的取值如下:
利用MCMC法生成给定时间长度负荷类别变化时间序列(给定时间上限)步骤如下:
1)随机产生一个在区间[0,1]内的数u;
2)将u与Pcum的第p行元素(元素行号与当前所处状态号相同)进行比较,假设u落在累积状态转移率矩阵Pcum的Pcum,pq到Pcum,p(q+1)这个范围内,则认为负荷状态的序列的下一时刻状态为q;
3)如果生成的序列已经满足时间长度的要求,则转步骤1)继续,否则,令当前状态变为q,转步骤1)步继续,直到达到时间要求;
4)每次模拟只对相邻故障时间间隔内的负荷状态进行模拟。
但是,本文所解决的问题有其特殊性,因为故障时间点随机出现,各时间点之间的时间间隔不统一,相邻两个时间点之间可能还有其他的负荷状态,所以要统一步长为常数△t,累积状态转移率矩阵的处理相应也要进行改进处理。
转移概率矩阵中,P12表示经过一个步长的时间后,从1状态转化为2状态的概率,也可以表示,现在为状态2,一个步长的时间之前的状态为1的概率,所以
式中,P12,k为已知△t时刻,负荷状态为1,则第(k+1)△t时刻,负荷处于第2种状态的概率。:
P12,kt=P21,k为已知第(k+1)△t时刻,负荷状态为1,则△t时刻,负荷处于第2种状态的概率。。
问题可以细化为假设已知△t、10△t两个时刻的负荷状态分别为a与b,已知马尔柯夫链P,采用蒙特卡洛仿真内插确定[△t,10△t]时间段内各时间点负荷状态。以确定5△t时刻的负荷情况为例,P4、(P5)T分别为首、末两端与5△t时刻负荷状态的转换概率分布。
根据马尔可夫性,△t时刻的情况与5△t时刻有关,而5△t时刻的情况与△t时刻有关,即一个马尔柯夫过程的正向和反向都具有马尔柯夫性,所以,5△t时刻的情况应该结合△t、10△t两个时刻的负荷状态进行确定。
对P4进行处理,令
Pcum1中最后一列元素为1,元素的取值如下:
对(P5)T进行处理,令
Pcum2中元素的取值如下:
再将Pcum2各行归一化。
从△t、10△t两个时刻对5△t时刻的情况进行蒙特卡罗仿真,将Pcum=(Pcum1+Pcum2)/2,然后利用前文所述的利用MCMC法生成给定时间长度负荷类别变化时间序列的处理方法进行仿真。
最后,执行步骤06,用Baum-Welch算法,基于HMM模型识别负荷特性的变化。步骤如下:
1)初始模型(待训练模型)λ0,
2)基于λ0以及观察值序列O,训练新模型λ;
3)如果logP(O|λ)-logP(O|λ0)<D,说明训练已经达到预期效果,算法结束,式中,P(O|λ)表示在实际模型为λ的情况下,观察序列为O的概率,D为设定门阈值;
4)否则,令λ0=λ,继续第2步工作。
给定模型λ和观察序列条件下,从i到j的转移概率定义为ξt(i,j)。
式(6-1)中,αt(i)=P(o1,o2,…,ot,xt=i),aij=P(xt+1=j|xt=i),bj(xt+1)=P(ot+1|xt+1=j),βt+1(j)=P(ot+1,ot+2,…,oT|qt=i,λ)。
第二个辅助变量为为t时刻处于状态Si的概率,为整个过程中从状态Si转出次数的预期,为从Si跳转到Sj次数的预期。
训练过程:为当t=1时,处于qi的概率。
与分别为状态转移概率矩阵A;观察值概率矩阵B中的元素的估计值。当收敛或达到预置次数时,即可得到状态转移概率矩阵A作为本发明的最终结果。
其中,步骤3、步骤4和步骤5中,建立概率转移矩阵并以此进行内插与外推均在MATLAB仿真平台中完成。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,包括以下步骤:
步骤一:找出电压降落的时间点,在电压降落的时间点对应的扰动时刻,进行负荷动特性提取与分类;
步骤二:判断负荷类别间的变化是否具有马氏性,如果是,则进行步骤三,否则,结束;
步骤三:将所有数据按时间平均分段,对每段数据基于最大似然的思想建立马氏链的概率转移矩阵;
步骤四:判断数字特征即概率转移矩阵的稳态分布作为一个向量是否改变,当向量的中心距超过规定门阈值时,认为数字特征已经改变,依据概率转移矩阵相应的数字特征对时间段的负荷数据进行聚类,对有改变了数据特征的各时间段负荷数据求概率转移矩阵,然后转入步骤五,如果数字特征未改变,则直接转入步骤五;
步骤五:利用概率转移矩阵生成累积状态转移率矩阵,用马尔柯夫蒙特卡罗仿真对负荷状态在整个时段内的变化进行内插与外推,得到完整的负荷数据序列;
步骤六:利用隐马尔科夫模型处理步骤五的结果即反映负荷类别转化的序列,得到更为精确的负荷类别变化数据。
2.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤四中,数字特征取为各类别的稳态分布:
式中,P为概率转移矩阵,pij为概率转移矩阵的第i行第j列的元素,函数的返回值为矩阵的任意行矢量,π*表示第一步求取的最优稳态分布,求J的最小值。
3.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤四中,依据矩阵相应的数字特征对时间段的负荷数据进行聚类:设x为聚类向量,将目标函数设为样本到所属类别中心的距离,采用模糊C均值聚类,用迭代的方法求取最小目标函数。
4.如权利要求3所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,用迭代的方法求取最小目标函数:
FCM算法把样本空间x={x1,x2,…,xn}分为c类,2≤c≤n,任一样本点xi不会严格划分为某一类,FCM用模糊划分,使用隶属度uij(0≤uij≤1)来确定样本点xi属于第j(0≤j≤c)类的程度。如样本空间x的一个模糊子集所对应的隶属度矩阵是一个模糊隶属度矩阵,用U={Uij}表示,隶属度矩阵U具有如下性质:
式中,uij为隶属度;
FCM算法就是在式(1)的约束条件下使目标函数J最小化,即:
式中,m∈[1,∞)是模糊加权系数;cj是c类中第j类的聚类中心;是xi到cj的欧拉距离。
5.如权利要求4所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,FCM算法就是在式(1)的约束条件下使目标函数J最小化,其步骤如下:
(1)给定聚类数c,模糊加权系数m,迭代停止阈值ε;设迭代次数k=0,最大迭代次数kmax;用值在[0,1]间的随机数初始化隶属矩阵U,使其满足式(1)中的约束条件;
(2)由式(2)计算聚类中心
式中,uij为隶属度,xi为样本;
(3)由聚类中心cj更新隶属度矩阵U(k+1),即
式中,dij为样本和各类中心的距离,上标(k+1)表示第k次迭代的相应数值;
(4)给定收敛的判别精度ε>0,若||U(k+1)-U(k)||≤ε,停止迭代;否则置k=k+1,并返回步骤(2);
(5)得到x的一个最优模糊C划分隶属度矩阵U={Uij}和聚类中心c={cij}。
6.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤五中,在进行内插和外推时,用概率转移矩阵生成累积状态转移率矩阵:
Pcum中元素的取值如下:
其中,pij为概率转移矩阵P的第i行第j列的元素,经过本步骤的处理,累积状态转移率矩阵的Pcum,il为P的所有第i行前j列元素之和。
7.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤五中,利用马尔柯夫蒙特卡罗仿真MCMC法生成给定时间长度负荷类别变化时间序列步骤如下:
1)随机产生一个在区间[0,1]内的数u;
2)将u与Pcum的第p行元素进行比较,假设u落在累积状态转移率矩阵Pcum的Pcum,pq到Pcum,p(q+1)这个范围内,则认为负荷状态的序列的下一时刻状态为q;
3)如果生成的序列已经满足时间长度的要求,则转步骤1)继续,否则,令当前状态变为q,转步骤1)继续,直到达到时间要求。
8.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤五中,因为故障时间点随机出现,各时间点之间的时间间隔不统一,相邻两个时间点之间可能还有其他的负荷状态,所以要统一步长为常数△t,累积状态转移率矩阵的处理相应也要进行改进:
转移概率矩阵中,P12表示经过一个步长的时间后,从1状态转化为2状态的概率,也可以表示,现在为状态2,一个步长的时间之前的状态为1的概率,所以
式中,P12,k为已知△t时刻,负荷状态为1,则第(k+1)△t时刻,负荷处于第2种状态的概率。
P12,kt=P21,k为已知第(k+1)△t时刻,负荷状态为1,则△t时刻,负荷处于第2种状态的概率。
9.如权利要求1所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,所述步骤六中,隐马尔科夫模型HMM是一个双重随机过程,它有两个组成部分:
a马尔可夫链----描述类别的转移,用转移概率描述;
b一般随机过程----描述状态与观察序列间的关系,用观察值概率描述;
一个隐马尔科夫模型HMM常用五元数组来描述:λ=(Q,O,A,B,π),其中:状态的有限集合Q={q1,q2,…,qn,},n为马氏链的状态数;观察值的有限集合O={o1,o2,…,om},m为每种类别对应的可能的观察值数目;状态转移概率矩阵A;观察值概率矩阵B;初始概率分布向量π。
10.如权利要求9所述的一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法,其特征是,用Baum-Welch算法,基于HMM模型识别负荷特性的变化,步骤如下:
1)建立初始模型即待训练模型λ0;
2)基于λ0以及观察值序列O,训练新模型λ;
3)如果logP(O|λ)-logP(O|λ0)<D,说明训练已经达到预期效果,算法结束,式中,P(O|λ)表示在实际模型为λ的情况下,观察序列为O的概率,D为设定门阈值;
4)否则,令λ0=λ,继续步骤2)工作;
训练过程前,需要定义两个辅助变量,第一个辅助变量为ξt(i,j),表示给定模型λ和观察序列条件下,从i到j的转移概率:
式中,αt(i)=P(o1,o2,…,ot,xt=i),aij=P(xt+1=j|xt=i),bj(xt+1)=P(ot+1|xt+1=j),βt+1(j)=P(ot+1,ot+2,…,oT|qt=i,λ);
第二个辅助变量为为t时刻处于状态Si的概率,为整个过程中从状态Si转出次数的预期,为从Si跳转到Sj次数的预期;
训练过程:为当t=1时,处于qi的概率;
分别为状态转移概率矩阵A、观察值概率矩阵B中的元素的估计值,当收敛或达到预置次数时,即可得到状态转移概率矩阵A作为最终结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510733268.0A CN105205502B (zh) | 2015-10-30 | 2015-10-30 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510733268.0A CN105205502B (zh) | 2015-10-30 | 2015-10-30 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105205502A true CN105205502A (zh) | 2015-12-30 |
CN105205502B CN105205502B (zh) | 2019-01-01 |
Family
ID=54953174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510733268.0A Active CN105205502B (zh) | 2015-10-30 | 2015-10-30 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105205502B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107392323A (zh) * | 2017-05-27 | 2017-11-24 | 中国电子科技集团公司第三十六研究所 | 一种两阶段衰退与冲击竞争风险的综合预防维修方法 |
CN107480817A (zh) * | 2017-08-09 | 2017-12-15 | 浙江群力电气有限公司 | 一种电压波动指标测量方法及系统 |
CN107834538A (zh) * | 2017-10-12 | 2018-03-23 | 浙江大学 | 基于序优化和马氏链的多状态电力系统冗余优化方法 |
CN108288039A (zh) * | 2018-01-25 | 2018-07-17 | 浙江群力电气有限公司 | 电压暂升和暂降的电能质量扰动识别方法、装置及设备 |
CN108830499A (zh) * | 2018-06-26 | 2018-11-16 | 广东石油化工学院 | 一种非侵入式负荷分解方法及系统 |
CN108964276A (zh) * | 2018-07-16 | 2018-12-07 | 河南理工大学 | 支持自动需求响应的用电管控终端和系统及负荷辨识方法 |
CN109301877A (zh) * | 2018-09-13 | 2019-02-01 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
CN109344459A (zh) * | 2018-09-13 | 2019-02-15 | 北京应用物理与计算数学研究所 | 一种蒙特卡罗程序的可视建模与转换方法及系统 |
CN110137951A (zh) * | 2019-05-16 | 2019-08-16 | 广东电力交易中心有限责任公司 | 基于节点电价的市场分区方法及装置 |
CN111309973A (zh) * | 2020-01-21 | 2020-06-19 | 杭州安脉盛智能技术有限公司 | 基于改进马尔可夫模型和改进k最近邻的缺失值填补方法 |
CN113408155A (zh) * | 2021-08-03 | 2021-09-17 | 中国人民解放军海军航空大学青岛校区 | 战时航材需求预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390117A (zh) * | 2013-08-08 | 2013-11-13 | 山东大学 | 一种电力负荷动特性的特征提取方法 |
US20130325782A1 (en) * | 2012-05-31 | 2013-12-05 | Nec Corporation | Latent variable model estimation apparatus, and method |
CN104218620A (zh) * | 2014-09-26 | 2014-12-17 | 国家电网公司 | 基于伪序贯蒙特卡洛模拟的主动配电网可靠性分析方法 |
CN104734171A (zh) * | 2015-04-16 | 2015-06-24 | 合肥工业大学 | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 |
-
2015
- 2015-10-30 CN CN201510733268.0A patent/CN105205502B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130325782A1 (en) * | 2012-05-31 | 2013-12-05 | Nec Corporation | Latent variable model estimation apparatus, and method |
CN103390117A (zh) * | 2013-08-08 | 2013-11-13 | 山东大学 | 一种电力负荷动特性的特征提取方法 |
CN104218620A (zh) * | 2014-09-26 | 2014-12-17 | 国家电网公司 | 基于伪序贯蒙特卡洛模拟的主动配电网可靠性分析方法 |
CN104734171A (zh) * | 2015-04-16 | 2015-06-24 | 合肥工业大学 | 一种用于配电网可靠性评估的电动汽车充电站建模方法及其应用 |
Non-Patent Citations (2)
Title |
---|
TONG WU等: "Markov Chain Monte Carlo Method for the Modeling of Wind Power Time Series", 《IEEE》 * |
于鹏 等: "含风电功率时域特性的风电功率序列建模方法", 《中国电机工程学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107392323B (zh) * | 2017-05-27 | 2020-05-08 | 中国电子科技集团公司第三十六研究所 | 一种两阶段衰退与冲击竞争风险的综合预防维修方法 |
CN107392323A (zh) * | 2017-05-27 | 2017-11-24 | 中国电子科技集团公司第三十六研究所 | 一种两阶段衰退与冲击竞争风险的综合预防维修方法 |
CN107480817A (zh) * | 2017-08-09 | 2017-12-15 | 浙江群力电气有限公司 | 一种电压波动指标测量方法及系统 |
CN107834538A (zh) * | 2017-10-12 | 2018-03-23 | 浙江大学 | 基于序优化和马氏链的多状态电力系统冗余优化方法 |
CN108288039A (zh) * | 2018-01-25 | 2018-07-17 | 浙江群力电气有限公司 | 电压暂升和暂降的电能质量扰动识别方法、装置及设备 |
CN108830499A (zh) * | 2018-06-26 | 2018-11-16 | 广东石油化工学院 | 一种非侵入式负荷分解方法及系统 |
CN108830499B (zh) * | 2018-06-26 | 2021-10-15 | 广东石油化工学院 | 一种非侵入式负荷分解方法及系统 |
CN108964276A (zh) * | 2018-07-16 | 2018-12-07 | 河南理工大学 | 支持自动需求响应的用电管控终端和系统及负荷辨识方法 |
CN109301877B (zh) * | 2018-09-13 | 2021-08-17 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
CN109344459A (zh) * | 2018-09-13 | 2019-02-15 | 北京应用物理与计算数学研究所 | 一种蒙特卡罗程序的可视建模与转换方法及系统 |
CN109301877A (zh) * | 2018-09-13 | 2019-02-01 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
CN109344459B (zh) * | 2018-09-13 | 2022-11-29 | 北京应用物理与计算数学研究所 | 一种蒙特卡罗程序的可视建模与转换方法及系统 |
CN110137951A (zh) * | 2019-05-16 | 2019-08-16 | 广东电力交易中心有限责任公司 | 基于节点电价的市场分区方法及装置 |
CN110137951B (zh) * | 2019-05-16 | 2021-06-08 | 广东电力交易中心有限责任公司 | 基于节点电价的市场分区方法及装置 |
CN111309973A (zh) * | 2020-01-21 | 2020-06-19 | 杭州安脉盛智能技术有限公司 | 基于改进马尔可夫模型和改进k最近邻的缺失值填补方法 |
CN111309973B (zh) * | 2020-01-21 | 2024-01-05 | 杭州安脉盛智能技术有限公司 | 基于改进马尔可夫模型和改进k最近邻的缺失值填补方法 |
CN113408155A (zh) * | 2021-08-03 | 2021-09-17 | 中国人民解放军海军航空大学青岛校区 | 战时航材需求预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105205502B (zh) | 2019-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105205502A (zh) | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 | |
Liu et al. | A novel hybrid methodology for short-term wind power forecasting based on adaptive neuro-fuzzy inference system | |
Tian et al. | Multi-step short-term wind speed prediction based on integrated multi-model fusion | |
CN111881505B (zh) | 一种基于ga-rbf算法的既有建筑多目标优化改造决策方法 | |
CN107730044A (zh) | 一种可再生能源发电和负荷的混合预测方法 | |
CN109117951A (zh) | 基于bp神经网络的概率潮流在线计算方法 | |
CN108694467A (zh) | 一种对配电网线损率进行预测的方法及系统 | |
CN110175386A (zh) | 变电站电气设备温度预测方法 | |
CN103676881A (zh) | 一种半导体生产线动态瓶颈分析方法 | |
CN103488869A (zh) | 一种最小二乘支持向量机的风力发电短期负荷预测方法 | |
CN104636985A (zh) | 一种改进bp神经网络的输电线路无线电干扰预测方法 | |
CN104133372A (zh) | 基于模糊神经网络的室温控制算法 | |
CN105260786A (zh) | 一种电力推进系统仿真可信度评估模型综合优化方法 | |
CN102768701A (zh) | 基于量子遗传算法的高压开关柜绝缘子电场优化方法 | |
CN109255726A (zh) | 一种混合智能技术的超短期风功率预测方法 | |
CN104217258A (zh) | 一种电力负荷条件密度预测方法 | |
CN103927460A (zh) | 一种基于rbf的风电场短期风速预测方法 | |
CN106096723A (zh) | 一种基于混合神经网络算法的用于复杂工业产品性能评估方法 | |
Guoqiang et al. | Study of RBF neural network based on PSO algorithm in nonlinear system identification | |
CN104134011A (zh) | 一种配电网小水电接入最优接纳能力的计算方法 | |
CN109840335A (zh) | 基于强化t-s模糊神经网络的径向锻造凹坑预测优化方法 | |
CN108898245A (zh) | 一种用于风电机组零部件的需求预测方法 | |
CN115829258A (zh) | 基于多项式混沌近似动态规划的电力系统经济调度方法 | |
CN105207233A (zh) | 基于Metropolis-Hastings与PSO相结合的无功优化方法 | |
El-Dahshan et al. | Artificial neural network and genetic algorithm hybrid technique for nucleus–nucleus collisions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |