CN103413053A - 一种基于内点法的电力***抗差状态估计方法 - Google Patents
一种基于内点法的电力***抗差状态估计方法 Download PDFInfo
- Publication number
- CN103413053A CN103413053A CN2013103676825A CN201310367682A CN103413053A CN 103413053 A CN103413053 A CN 103413053A CN 2013103676825 A CN2013103676825 A CN 2013103676825A CN 201310367682 A CN201310367682 A CN 201310367682A CN 103413053 A CN103413053 A CN 103413053A
- Authority
- CN
- China
- Prior art keywords
- delta
- dtri
- electric system
- state estimation
- variable
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于内点法的电力***抗差状态估计方法,属于电力***运行和控制技术领域。该方法包括:获取量测数据、拓扑分析以及计算电网参数;建立指数函数抗差状态估计模型;引入残差变量,对抗差状态估计模型进行等价转换,规范成内点法容易求解的形式;利用内点法对等价转换后的状态估计模型进行求解,并对海森矩阵进行数值近似。本发明通过对状态估计模型的等价转换,降低了电力***抗差状态估计的难度;通过对海森矩阵进行近似,大大提高了电力***抗差状态估计的计算效率;状态估计结果严格满足零注入等式约束;计算方法具有很强的鲁棒性。
Description
技术领域
本发明涉及一种基于内点法的电力***抗差状态估计方法,属于电力***运行和控制技术领域。
背景技术
电力***的状态估计就是以数据采集与监视控制***采集的实时量测为数据输入源,剔除其中的不良数据之后,结合电网模型,按照特定的状态估计模型,对电网的状态量进行最优估算。状态估计在能量管理***中处于极为重要的地位,是后续一系列高级应用的重要基石。
传统的电力***状态估计需要进行不良数据的辨识以及剔除,才能进行准确的估计。但现在希望能在进行状态估计的过程中自动降低不良数据所导致的估计偏差,因此可以采用抗差状态估计的模型。
在求解的算法方面,传统的牛顿拉夫逊方法存在收敛域小、对初值要求苛刻、无法处理不等式约束等缺点,在现场应用中,在重载的情况下,牛顿法有可能出现即便平启动也无法收敛的情况,同时也无法考虑节点负荷上下限等等的不等式约束,因此必须寻求更优的解法。而现代内点法是公认的优秀的算法。
发明内容
本发明的目的是提出一种基于内点法的电力***抗差状态估计方法,将现代内点法应用到抗差状态估计中,以解决实际工程中状态估计收敛域窄的问题。
本发明提出一种基于内点法的电力***抗差状态估计方法,包括以下步骤:
(1)从电力***的数据采集中心实时获取量测数据,量测数据包括节点电压幅值、节点有功功率、节点无功功率、支路有功功率和支路无功功率,并存储量测数据,根据量测数据进行拓扑分析,得到电力***结构图,根据量测数据计算电力***导纳矩阵Y;
(2)根据步骤(1)中量测数据、电力***结构图和电力***导纳矩阵Y,建立初始电力***抗差状态估计模型如下:
s.t. c(s)=O
其中,s为电力***状态量,s=[θ,V]T,mi(s)为电力***第i个量测方程的估计值,为上述量测数据中的第i个量测值,σ是核估计法的窗宽,窗宽σ的取值范围为:O<σ<1,c(s)=O表示零注入等式约束;
(3)对上述初始电力***抗差状态估计模型进行等价转换,得到电力***抗差状态估计模型如下:
s.t.c(s,t)=O
其中,mi(s)为电力***第i个量测方程的估计值,t为残差变量,残差变量中的第i个分量为:i=1,2,...,nmeas,nmeas为上述量测数据的个数,s为电力***状态量,s=[θ,V]T,夕是节点电压的相角,V是节点电压的幅值,T是矩阵转置;
(4)采用原-对偶内点法,对上述电力***抗差状态估计模型进行求解,具体实现过程如下:
(4-1)将上述电力***抗差状态估计模型改写为以下标准形式:
min f(x)
s.f h(x)=O
其中,x是电力***状态量s和残差变量t合并后的复合向量,x=[s,t]T,函数h(x)包含量测方程mi(s)的约束以及零注入约束,g(x)为电力***节点负荷与复合向量x之间的函数关系,为用户设定的电力***负荷上限,g为用户设定的电力***负荷下限;
(4-2)将上述标准形式中的不等式约束简化为等式约束,得到简化后的电力***抗差状态估计模型:
min f(x)
s.t h(x)=O
g(x)-1-g=0
(1,u)≥O
其中,(1,u)∈Rr为将不等式约束转换为等式约束的松弛向量,R为实数集,r为函数g(x)的维数;
(4-3)根据上述简化后的电力***抗差状态估计模型,构建一个增广拉格朗日函数如下:
(4-4)设置迭代初始值,将原始变量x中,电力***状态量s中的电压幅值初值设为1,电压相角初值设为0,将原始变量x中的残差变量t的初始值设为0,将对偶变量的初始值设为0,设置迭代次数的初始值k=0,设置最大迭代次数Kmax,Kmax的取值范围为:10~100,对于电力***中的零注入母线,设定残差变量的第i个分量
(4-5)对迭代次数k进行判断,若k<Kmax,则进行步骤(4-6),若k≥Kmax,则结束计算,并输出计算不收敛;
(4-6)根据下式计算简化后的电力***抗差状态估计模型的互补间隙CGap:
其中,l和u为松弛变量,z和w为对偶变量,r为函数g(x)的维数;
设定一个电力***抗差状态估计计算精度ε,一般可设为10-6,将互补间隙CGap与计算精度ε进行比较,若CGap<ε,则输出x,包括电力***状态量和残差变量,并结束计算,若CGap≥ε,则进行步骤(4-7);
(4-7)根据下式计算扰动因子μ:
其中τ为中心参数,取值范围为τ∈(0,1],r为函数g(x)的维数,CGap为上述互补间隙;
(4-8)按照以下卡罗需-库恩-塔克方程组,计算上述增广拉格朗日函数对原-对偶内点法中各原始变量和对偶变量的一阶偏导数:
(4-9)根据以下修正方程,计算得到在第k次迭代时原始变量x的修正量Δx、对偶变量y、z、w的修正量Δy、Δz、Δw和松弛变量1、u的修正量Δl、Δu:
其中:
设定一个转换阈值εswitch,一般设为计算精度ε的2~3倍即可,εswitch的取值范围为2×10-6~3×10-6,对原始-对偶互补间隙CGap进行判断,若CGap<εswitch,则由于x=[s,t]T,
其中nmeas为上述量测数据的个数,ti为残差变量的第i个分量,σ为核估计法的窗宽;
若CGap≥ε,则各元素采用下式的近似表达:
(4-11)根据步骤(4-10)的修正步长,更新上述原始变量和对偶变量:
使
(4-12)重复步骤(4-5)-步骤(4-12)。
本发明提出的基于内点法的电力***抗差状态估计方法,其优点是:
1、本发明状态估计方法对目标函数的海森矩阵进行数值近似之后,状态估计的计算效率较高,因而能够基本满足状态估计实时性的需求。
2、本发明状态估计方法中采用了内点法,因此收敛域广,即便在电力***运行负荷过大或出现故障等异常情况下,也能迅速得到状态估计的最终结果。
3、本发明状态估计方法,计算过程中严格满足零注入等式约束,并能考虑多种复杂的不等式约束,因此能够得到更加符合实际运行状态的电力***状态估计结果。
具体实施方式
本发明提出一种基于内点法的电力***抗差状态估计方法,包括以下步骤:
(1)从电力***的数据采集中心实时获取量测数据,量测数据包括节点电压幅值、节点有功功率、节点无功功率、支路有功功率和支路无功功率,并存储量测数据,根据量测数据进行拓扑分析,得到电力***结构图,根据量测数据计算电力***导纳矩阵Y;
(2)根据步骤(1)中量测数据、电力***结构图和电力***导纳矩阵Y,建立初始电力***抗差状态估计模型如下:
s.t. c(s)=0
其中,s为电力***状态量,s=[θ,V]T,mi(s)为电力***第i个量测方程的估计值,为上述量测数据中的第i个量测值,σ是核估计法的窗宽,窗宽σ的取值范围为:0<σ<1,c(s)=0表示零注入等式约束;
(3)对上述初始电力***抗差状态估计模型进行等价转换,得到电力***抗差状态估计模型如下:
s.t.c(s,t)=0
其中,mi(s)为电力***第i个量测方程的估计值,t为残差变量,残差变量中的第i个分量为:i=1,2,...,nmeas,nmeas为上述量测数据的个数,s为电力***状态量,s=[θ,V]T,θ是节点电压的相角,V是节点电压的幅值,T是矩阵转置;
(4)采用原-对偶内点法,对上述电力***抗差状态估计模型进行求解,具体实现过程如下:
(4-1)将上述电力***抗差状态估计模型改写为以下标准形式:
minf(x)
s.t h(x)=0
其中,x是电力***状态量s和残差变量t合并后的复合向量,x=[s,t]T,函数h(x)包含量测方程mi(s)的约束以及零注入约束,g(x)为电力***节点负荷与复合向量x之间的函数关系,为用户设定的电力***负荷上限,g为用户设定的电力***负荷下限;
(4-2)将上述标准形式中的不等式约束简化为等式约束,得到简化后的电力***抗差状态估计模型:
min f(x)
s.t h(x)=0
g(x)-l-g=0
(l,u)≥0
其中,(l,u)∈Rr为将不等式约束转换为等式约束的松弛向量,R为实数集,r为函数g(x)的维数;
(4-3)根据上述简化后的电力***抗差状态估计模型,构建一个增广拉格朗日函数如下:
(4-4)设置迭代初始值,将原始变量x中,电力***状态量s中的电压幅值初值设为1,电压相角初值设为0,将原始变量x中的残差变量t的初始值设为0,将对偶变量的初始值设为0,设置迭代次数初始值k=0,设置最大迭代次数Kmax,Kmax的取值范围为:10~100,对于电力***中的零注入母线,设定残差变量的第i个分量
(4-5)对迭代次数k进行判断,若k<Kmax,则进行步骤(4-6),若k≥Kmax,则结束计算,并输出计算不收敛;
(4-6)根据下式计算简化后的电力***抗差状态估计模型的互补间隙CGap:
其中,l和u为松弛变量,z和w为对偶变量,r为函数g(x)的维数;
设定一个电力***抗差状态估计计算精度ε,将互补间隙CGap与计算精度ε进行比较,若CGap<ε,则输出x,包括电力***状态量和残差变量,并结束计算,若CGap≥ε,则进行步骤(4-7);
(4-7)根据下式计算扰动因子μ:
其中τ为中心参数,取值范围为τ∈(0,1],r为函数g(x)的维数,CGap为上述互补间隙;
(4-8)按照以下卡罗需-库恩-塔克方程组,计算上述增广拉格朗日函数对原-对偶内点法中各原始变量和对偶变量的一阶偏导数:
(4-9)根据以下修正方程,计算得到在第k次迭代时原始变量x的修正量Δx、对偶变量y、z、w的修正量Δy、Δz、Δw和松弛变量l、u的修正量Δl、Δu:
其中:
设定一个转换阈值εswitch,εswitch的取值范围为2×10-6~3×10-6,对原始-对偶互补间隙CGap进行判断,若CGap<εswitch,则由于x=[s,t]T,在计算上述修正方程中的H(·)时,在中,各元素采用下式精确表达:
其中nmeas为上述量测数据的个数,ti为残差变量的第i个分量,σ为核估计法的窗宽;
(4-11)根据步骤(4-10)的修正步长,更新上述原始变量和对偶变量:
使
(4-12)重复步骤(4-5)-步骤(4-12)。
Claims (1)
1.一种基于内点法的电力***抗差状态估计方法,其特征在于,该方法包括以下步骤:
(1)从电力***的数据采集中心实时获取量测数据,量测数据包括节点电压幅值、节点有功功率、节点无功功率、支路有功功率和支路无功功率,并存储量测数据,根据量测数据进行拓扑分析,得到电力***结构图,根据量测数据计算电力***导纳矩阵Y;
(2)根据步骤(1)中量测数据、电力***结构图和电力***导纳矩阵Y,建立初始电力***抗差状态估计模型如下:
s.t. c(s)=0
其中,s为电力***状态量,s=[θ,V]T,mi(s)为电力***第i个量测方程的估计值,为上述量测数据中的第i个量测值,σ是核估计法的窗宽,窗宽σ的取值范围为:0<σ<1,c(s)=0表示零注入等式约束;
(3)对上述初始电力***抗差状态估计模型进行等价转换,得到电力***抗差状态估计模型如下:
s.t.c(s,t)=0
其中,mi(s)为电力***第i个量测方程的估计值,t为残差变量,残差变量t中的第i个分量为:i=1,2,...,nmeas,nmeas为上述量测数据的个数,s为电力***状态量,s=[θ,V]T,θ为节点电压的相角,V是节点电压的幅值,T是矩阵转置;
(4)采用原-对偶内点法,对上述电力***抗差状态估计模型进行求解,具体实现过程如下:
(4-1)将上述电力***抗差状态估计模型改写为以下标准形式:
min f(x)
s.t h(x)=0
其中,x是电力***状态量s和残差变量t合并后的复合向量,x=[s,t]T,函数h(x)包含量测方程mi(s)的约束以及零注入约束,g(x)为电力***节点负荷与复合向量x之间的函数关系,为用户设定的电力***负荷上限,g为用户设定的电力***负荷下限;
(4-2)将上述标准形式中的不等式约束简化为等式约束,得到简化后的电力***抗差状态估计模型:
min f(x)
s.t h(x)=0
g(x)-1-g=0
(l,u)≥0
其中,(l,u)∈Rr为将不等式约束转换为等式约束的松弛向量,R为实数集,r为函数g(x)的维数;
(4-3)根据上述简化后的电力***抗差状态估计模型,构建一个增广拉格朗日函数如下:
(4-4)设置迭代初始值,将原始变量x中,电力***状态量s中的电压幅值初值设为1,电压相角初值设为0,将原始变量x中的残差变量t的初始值设为0,将对偶变量的初始值设为0,设置迭代次数的初始值k=0,设置最大迭代次数Kmax,Kmax的取值范围为:10~100,对于电力***中的零注入母线,设定残差变量的第i个分量
(4-5)对迭代次数k进行判断,若k<Kmax,则进行步骤(4-6),若k≥Kmax,则结束计算,并输出计算不收敛;
(4-6)根据下式计算上述简化后的电力***抗差状态估计模型的互补间隙CGap:
其中,l和u为松弛变量,z和w为对偶变量,r为函数g(x)的维数;
设定一个电力***抗差状态估计计算精度ε,将互补间隙CGap与计算精度ε进行比较,若CGap<ε,则输出x,包括电力***状态量和残差变量,并结束计算,若CGap≥ε,则进行步骤(4-7);
(4-7)根据下式计算扰动因子μ:
其中τ为中心参数,取值范围为τ∈(0,1],r为函数g(x)的维数,CGap为上述互补间隙;
(4-8)按照以下卡罗需-库恩-塔克方程组,计算上述增广拉格朗日函数对原-对偶内点法中各原始变量和对偶变量的一阶偏导数:
(4-9)根据以下修正方程,计算得到在第k次迭代时原始变量x的修正量Δx、对偶变量y、z、w的修正量Δy、Δz、Δw和松弛变量l、u的修正量Δl、Δu:
其中:
设定一个转换阈值εswitch,εswitch的取值范围为2×10-6~3×10-6,对原始-对偶互补间隙CGap进行判断,若CGap<εswitch,则由于x=[s,t]T,在计算上述修正方程中的H(·)时,在中,各元素采用下式精确表达:
其中nmeas为上述量测数据的个数,ti为残差变量的第i个分量,σ为核估计法的窗宽;
若CGap≥ε,则各元素采用下式的近似表达:
(4-11)根据步骤(4-10)的修正步长,更新上述原始变量和对偶变量:
使 k=k+1
(4-12)重复步骤(4-5)~步骤(4-12)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310367682.5A CN103413053B (zh) | 2013-08-21 | 2013-08-21 | 一种基于内点法的电力***抗差状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310367682.5A CN103413053B (zh) | 2013-08-21 | 2013-08-21 | 一种基于内点法的电力***抗差状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103413053A true CN103413053A (zh) | 2013-11-27 |
CN103413053B CN103413053B (zh) | 2016-09-14 |
Family
ID=49606064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310367682.5A Active CN103413053B (zh) | 2013-08-21 | 2013-08-21 | 一种基于内点法的电力***抗差状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103413053B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103914613A (zh) * | 2014-03-04 | 2014-07-09 | 中国南方电网有限责任公司 | 电力***动态状态估计中异常情况的检测方法 |
CN104050355A (zh) * | 2014-03-25 | 2014-09-17 | 浙江工业大学 | 一种基于h∞滤波的机械臂外力估计方法 |
CN104102836A (zh) * | 2014-07-14 | 2014-10-15 | 国家电网公司 | 一种电力***快速抗差状态估计方法 |
CN104836230A (zh) * | 2015-06-02 | 2015-08-12 | 中国海洋石油总公司 | 一种海上油田群电网的状态估计方法 |
CN105322533A (zh) * | 2014-05-29 | 2016-02-10 | 河海大学 | 基于高斯-马尔科夫模型的自适应t型抗差状态估计方法 |
CN105552904A (zh) * | 2016-01-30 | 2016-05-04 | 清华大学 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
CN109274091A (zh) * | 2018-10-15 | 2019-01-25 | 同济大学 | 一种输配一体化并行状态估计方法 |
CN110428185A (zh) * | 2019-08-08 | 2019-11-08 | 河海大学 | 基于伪量测模型的电-热互联综合能源***抗差状态估计方法 |
CN112152198A (zh) * | 2020-08-12 | 2020-12-29 | 湖南大学 | 电力***低模型依赖性智能步长调整状态估计方法及*** |
CN114065118A (zh) * | 2021-11-02 | 2022-02-18 | 厦门大学 | 一种基于指数函数的电力***鲁棒状态估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599643A (zh) * | 2009-04-23 | 2009-12-09 | 清华大学 | 一种基于指数型目标函数的电力***抗差状态估计方法 |
-
2013
- 2013-08-21 CN CN201310367682.5A patent/CN103413053B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599643A (zh) * | 2009-04-23 | 2009-12-09 | 清华大学 | 一种基于指数型目标函数的电力***抗差状态估计方法 |
Non-Patent Citations (4)
Title |
---|
QIA DING等: ""Implementation of Interior Point Method Based Voltage/Reactive Power Optimization"", 《POWER ENGINEERING SOCIETY WINTER MEETING, 2000 IEEE 》 * |
Y-C.WU: ""Efficient two-level interior point method for optimal pumped hydrostorage scheduling exploiting the nonsparse matrix structure"", 《IEE PROC.-GENER. TRANSM. DISTRIB.》 * |
亓俊健等: ""电力***抗差状态估计研究综述"", 《电工电能新技术》 * |
卫志农等: ""含VSC-HVDC 的交直流***内点法最优潮流计算"", 《中国电机工程学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103914613A (zh) * | 2014-03-04 | 2014-07-09 | 中国南方电网有限责任公司 | 电力***动态状态估计中异常情况的检测方法 |
CN103914613B (zh) * | 2014-03-04 | 2017-02-01 | 中国南方电网有限责任公司 | 电力***动态状态估计中异常情况的检测方法 |
CN104050355B (zh) * | 2014-03-25 | 2017-06-23 | 浙江工业大学 | 一种基于h∞滤波的机械臂外力估计方法 |
CN104050355A (zh) * | 2014-03-25 | 2014-09-17 | 浙江工业大学 | 一种基于h∞滤波的机械臂外力估计方法 |
CN105322533A (zh) * | 2014-05-29 | 2016-02-10 | 河海大学 | 基于高斯-马尔科夫模型的自适应t型抗差状态估计方法 |
CN105322533B (zh) * | 2014-05-29 | 2017-11-03 | 河海大学 | 基于高斯‑马尔科夫模型的自适应t型抗差状态估计方法 |
CN104102836A (zh) * | 2014-07-14 | 2014-10-15 | 国家电网公司 | 一种电力***快速抗差状态估计方法 |
CN104836230A (zh) * | 2015-06-02 | 2015-08-12 | 中国海洋石油总公司 | 一种海上油田群电网的状态估计方法 |
CN105552904A (zh) * | 2016-01-30 | 2016-05-04 | 清华大学 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
CN109274091A (zh) * | 2018-10-15 | 2019-01-25 | 同济大学 | 一种输配一体化并行状态估计方法 |
CN110428185A (zh) * | 2019-08-08 | 2019-11-08 | 河海大学 | 基于伪量测模型的电-热互联综合能源***抗差状态估计方法 |
CN112152198A (zh) * | 2020-08-12 | 2020-12-29 | 湖南大学 | 电力***低模型依赖性智能步长调整状态估计方法及*** |
CN112152198B (zh) * | 2020-08-12 | 2022-06-21 | 湖南大学 | 电力***低模型依赖性智能步长调整状态估计方法及*** |
CN114065118A (zh) * | 2021-11-02 | 2022-02-18 | 厦门大学 | 一种基于指数函数的电力***鲁棒状态估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103413053B (zh) | 2016-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103413053A (zh) | 一种基于内点法的电力***抗差状态估计方法 | |
CN102801162B (zh) | 一种两阶段线性加权最小二乘电力***状态估计方法 | |
CN103208818B (zh) | 一种计及分布式电源的配网合环电流计算方法 | |
CN103248043B (zh) | 一种基于同步相角测量装置的电力***多区域分布式状态估计方法 | |
CN101599643A (zh) | 一种基于指数型目标函数的电力***抗差状态估计方法 | |
CN103454917B (zh) | 基于异步迭代模式的电力***分布式状态估计计算方法 | |
CN101964525B (zh) | 一种支持大规模电流量测的配电网状态估计方法 | |
CN102185308B (zh) | 一种考虑零注入量测等式约束的电力***状态估计方法 | |
CN102427229B (zh) | 基于修正牛顿法的带零注入约束的电力***状态估计方法 | |
CN103944165A (zh) | 一种大电网参数辨识估计方法 | |
CN102868157A (zh) | 一种基于最大指数绝对值目标函数的抗差状态估计方法 | |
CN105512502A (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
CN103474989A (zh) | 一种基于灵敏度分析的网络重构方法 | |
CN104037759B (zh) | 电力***注入转移分布因子的抗差最小二乘估计方法 | |
CN107749628A (zh) | 计及燃气机组无功调节及热电协调的多目标电压优化方法 | |
CN109494726A (zh) | 基于dlrnn神经网络的电力***稳定性在线评估方法 | |
CN102427227B (zh) | 考虑零注入约束的修正快速解耦电力***状态估计方法 | |
CN102280877B (zh) | 一种多量测断面的电力***不良支路参数辨识方法 | |
CN106099921A (zh) | 一种电力***时滞稳定裕度快速求解方法 | |
CN105305440B (zh) | 电力***状态的双曲余弦型最大指数绝对值抗差估计方法 | |
CN105977963B (zh) | 一种基于两阶段法的主配网一体化分布式状态估计方法 | |
CN106856334B (zh) | 一种考虑柔性直流控制特性的电力***状态估计方法 | |
CN109255541B (zh) | 一种基于最小二乘和一乘之和的配电网抗差状态估计方法 | |
CN103825270A (zh) | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 | |
CN110021931A (zh) | 一种计及模型不确定性的电力***辅助预测状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |