CN106936628B - 一种计及传感器故障的分数阶网络系统状态估计方法 - Google Patents
一种计及传感器故障的分数阶网络系统状态估计方法 Download PDFInfo
- Publication number
- CN106936628B CN106936628B CN201710082556.3A CN201710082556A CN106936628B CN 106936628 B CN106936628 B CN 106936628B CN 201710082556 A CN201710082556 A CN 201710082556A CN 106936628 B CN106936628 B CN 106936628B
- Authority
- CN
- China
- Prior art keywords
- fractional order
- moment
- network system
- follows
- sensor fault
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/142—Network analysis or design using statistical or mathematical methods
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/06—Management of faults, events, alarms or notifications
- H04L41/0631—Management of faults, events, alarms or notifications using root cause analysis; using analysis of correlation between notifications, alarms or events based on decision criteria, e.g. hierarchy, tree or time analysis
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Algebra (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种计及传感器故障的分数阶网络系统状态估计方法,用于分析分数阶网络系统因传感器故障而引起的数据丢包情况下的状态估计问题。该方法的具体步骤如下:首先对传感器故障情况下,发生数据随机丢包的系统模型进行了分析,建立了考虑数据随机丢包情况下的分数阶网络系统模型。接着,以传统的分数阶扩展卡尔曼滤波状态估计方法为基础,结合传感器故障情况下的分数阶网络系统模型,设计出了改进分数阶扩展卡尔曼状态估计方法。本方法适用于传感故障引起数据随机丢包情况下的分数阶网络系统状态估计问题,且流程简单易于实现。
Description
技术领域
本发明涉及一种计及传感器故障的分数阶网络系统状态估计方法,属于网络系统分析与控制技术领域。
背景技术
网络系统的分析与控制对于保证网络系统安全稳定的运行具有重要的意义。近年来,随着传感器技术的发展,网络系统的实时在线监测与控制成为了众多研究人员所关注的焦点。在现有的研究中,借助于传感器所获取的实时量测信息,通过设计动态的状态估计器,是实现网络系统实时分析与控制的主要途径。
一般情况下,现场数据通过传感器进行量测,然后通过信息传输通道传到控制中心,但是需要注意的是,传感器所量测的信息并不总是真实的,因为其会受到外界的干扰,以及信号的衰减,甚至传感器故障的影响。所以,在进行网络系统的动态估计器时必须计及量测信号发生丢包的情况。
分数阶网络系统由于可以更加精确的描述系统的结构,近年来被广泛应用于各个领域,如运用在电力系统网络中,可以更加精确的对电力系统中的节点电压和电流进行预测和估计。但是,在现有的分数阶网络系统研究中,计及传感器失败所引起的数据丢包现象主要集中于线性的分数阶网络,而对非线性分数阶网络计及传感器故障的分析和研究,国内外鲜有相关报道。为了进一步拓展分数阶网络的应用,本发明设计了计及传感器故障下的非线性分数阶网络系统的状态估计方法,并从理论上给予了证明。最后,实际的分数阶网络系统算例测试验证了本发明方法的有效性和实用性。
发明内容
发明目的:针对现有技术中存在的问题,本发明提供一种计及传感器失效的非线性分数阶网络状态估计方法。
技术方案:一种计及传感器故障的分数阶网络系统状态估计方法,包括如下部分:
1)计及传感器故障的分数阶网络系统建模
对于传感器故障下,系统量测数据发生随机丢包的离散非线性分数阶网络系统,其状态方程xk+1和输出方程yk分别为:
Δγxk+1=f(xk)+wk
yk=Γkh(xk)+vk
式中:xk+1表示k+1时刻的状态矢量,yk表示k时刻的输出矢量,f(·)和h(·)对应于两个可用泰勒级数展开的非线性函数,wk和vk分别为k时刻的系统噪声值和量测噪声值,二者相互独立无关,满足的协方差矩阵分别为Qk和Rk,式中γj和Γk计算公式如下
式中n≥0是分数阶阶次,j≥0代表不同时刻,是符合伯努利分布的二进制标量,其取值为0或1;期望和方差分别为πi,πi(1-πi),即满足(其中P(·)表示某件事发生的概率)
在建立传感器网络失败所引起量测信号丢包的模型之后,则可以通过如下方法对量测信号数据丢包情况下的非线性分数阶网络系统进行状态估计。
2)初始化k时刻的估计初始值和估计误差协方差Pk,估计时刻最大值N;
式中E(·)表示对某变量进行求期望运算,(·)T表示求矩阵转置。
3)计算k时刻的系统函数的雅克比矩阵,计算公式如下
式中表示求函数在变量处的偏导。
4)计算k时刻的反馈增益矩阵Kk,计算公式如下
式中[·]-1为矩阵逆运算符,是阿达玛乘积算子,Υ式中计算公式如下
5)计算k+1时刻的状态估计协方差矩阵Pk+1,计算公式如下
6)计算k+1时刻的状态估计值计算公式如下
7)若k+1<N则进行下一时刻的迭代估计;反之,则结束迭代,输出估计结果。
8)算法的证明过程如下
证明:估计误差ek+1可以表示为
基于泰勒级数展开可得到:
所以估计误差ek+1可以近似等价为:
对k+1时刻的状态估计协方差矩阵Pk+1推导化简可得到:
式中简化上式中的一项为:
并定义如下变量:
则k+1时刻的状态估计协方差矩阵Pk+1可进一步简化为:
通过完成观测增益的平方项,则:
通过联立上述两式,可得到:
进而可求得滤波反馈增益:
当且仅当估计误差协方差矩阵Pk+1取得最小值,此时协方差矩阵Pk+1为
所以k+1时刻的状态估计值可以通过如下公式求取
附图说明
图1为本发明实施例的方法流程图;
图2为实施例的量测信号图;
图3为实施例运用本发明和传统方法的状态估计结果图,其中(a)为本发明方法的状态估计结果图,(b)为传统方法的状态估计结果图;
图4为实施例采用本发明和传统方法的估计误差图,其中(a)为本发明方法的估计误差图,(b)为传统方法的估计误差图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,计及传感器故障的分数阶网络系统状态估计方法,方法在计算机中是依次按照如下步骤实现的:
1)初始化k时刻的估计初始值和估计误差协方差Pk,估计时刻最大值N;
2)计算k时刻的系统函数的雅克比矩阵,计算公式如下
3)计算k时刻的反馈增益矩阵Kk,计算公式如下
4)计算k+1时刻的状态估计协方差矩阵Pk+1,计算公式如下
5)计算k+1时刻的状态估计值计算公式如下
6)若k+1<N则进行下一时刻的迭代估计;反之,则结束迭代,输出估计结果。
为了验证本发明方法的有效性,下面介绍本发明的一个实施例,考虑如下非线性分数阶网络系统
Δγxk+1=sin(xk)+wk
式中分数阶阶次n=0.9,因传感器失效所引起的量测信号丢包率为系统噪声wk和量测噪声vk所满足的协方差矩阵分别为
在运用本发明方法对实施例非线性分数阶网络进行状态估计时,状态估计的初始值x0=0.9;初始状态估计误差协方差矩阵为P0=1,最大迭代估计时刻N=100。
分别运用本发明计及传感器故障的分数阶网络系统状态估计方法,以及传统分数阶卡尔曼滤波状态估计方法对实施例分数阶网络系统进行变量估计,不同算法的状态估计结果如图3所示,状态估计误差如图4所示。
综合图3和图4所示的测试结果,可以得出如下结论:由于传感器失效会引起量测信号丢失,所以在对分数阶网络系统进行状态估计器设计时必须计及量测信号丢失的情形。
Claims (1)
1.一种计及传感器故障的分数阶网络系统状态估计方法,其特征在于,包括如下步骤:
1)计及传感器故障的分数阶网络系统建模
对于传感器故障下,系统量测数据发生随机丢包的离散非线性分数阶网络系统,其状态方程xk+1和输出方程yk分别为:
Δγxk+1=f(xk)+wk
yk=Γkh(xk)+vk
式中:xk+1表示k+1时刻的状态矢量,yk表示k时刻的输出矢量,f(·)和h(·)对应于两个可用泰勒级数展开的非线性函数,wk和vk分别为k时刻的系统噪声值和量测噪声值,二者相互独立无关,满足的协方差矩阵分别为Qk和Rk,式中γj和Γk计算公式如下
式中n≥0是分数阶阶次,j≥0代表不同时刻,是符合伯努利分布的二进制标量,其取值为0或1;期望和方差分别为πi,πi(1-πi),满足P(·)表示某件事发生的概率
2)初始化k时刻的估计初始值和估计误差协方差Pk,估计时刻最大值N;
3)计算k时刻的系统函数的雅克比矩阵,计算公式如下
4)计算k时刻的反馈增益矩阵Kk,计算公式如下
5)计算k+1时刻的状态估计协方差矩阵Pk+1,计算公式如下
6)计算k+1时刻的状态估计值计算公式如下
7)若k+1<N则进行下一时刻的迭代估计;反之,则结束迭代,输出估计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710082556.3A CN106936628B (zh) | 2017-02-16 | 2017-02-16 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710082556.3A CN106936628B (zh) | 2017-02-16 | 2017-02-16 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106936628A CN106936628A (zh) | 2017-07-07 |
CN106936628B true CN106936628B (zh) | 2019-10-18 |
Family
ID=59424186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710082556.3A Active CN106936628B (zh) | 2017-02-16 | 2017-02-16 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106936628B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107977489A (zh) * | 2017-11-08 | 2018-05-01 | 南京邮电大学 | 一种复杂网络的保性能状态估计器的设计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101666825A (zh) * | 2009-08-28 | 2010-03-10 | 西南交通大学 | 一种基于电压稳定性局部指标的电网电压稳定在线监测方法 |
CN103248043A (zh) * | 2013-04-28 | 2013-08-14 | 西南交通大学 | 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法 |
CN104462015A (zh) * | 2014-11-26 | 2015-03-25 | 河海大学 | 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法 |
CN104794101A (zh) * | 2015-04-08 | 2015-07-22 | 河海大学 | 一种分数阶非线性系统状态估计方法 |
CN104820788A (zh) * | 2015-05-15 | 2015-08-05 | 河海大学 | 计及Lévy噪声的分数阶扩展卡尔曼滤波方法 |
CN105958470A (zh) * | 2014-10-20 | 2016-09-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN106355029A (zh) * | 2016-09-14 | 2017-01-25 | 国网天津市电力公司 | 一种配电系统关键量测解耦快速辨识的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120022713A1 (en) * | 2010-01-14 | 2012-01-26 | Deaver Sr Brian J | Power Flow Simulation System, Method and Device |
US9627886B2 (en) * | 2012-03-27 | 2017-04-18 | Mitsubishi Electric Research Laboratoriies, Inc. | State estimation for power system using hybrid measurements |
-
2017
- 2017-02-16 CN CN201710082556.3A patent/CN106936628B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101666825A (zh) * | 2009-08-28 | 2010-03-10 | 西南交通大学 | 一种基于电压稳定性局部指标的电网电压稳定在线监测方法 |
CN103248043A (zh) * | 2013-04-28 | 2013-08-14 | 西南交通大学 | 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法 |
CN105958470A (zh) * | 2014-10-20 | 2016-09-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN104462015A (zh) * | 2014-11-26 | 2015-03-25 | 河海大学 | 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法 |
CN104794101A (zh) * | 2015-04-08 | 2015-07-22 | 河海大学 | 一种分数阶非线性系统状态估计方法 |
CN104820788A (zh) * | 2015-05-15 | 2015-08-05 | 河海大学 | 计及Lévy噪声的分数阶扩展卡尔曼滤波方法 |
CN106355029A (zh) * | 2016-09-14 | 2017-01-25 | 国网天津市电力公司 | 一种配电系统关键量测解耦快速辨识的方法 |
Non-Patent Citations (4)
Title |
---|
"计及PMU量测的分步线性状态估计模型";厉超,卫志农,倪明;《电网技术》;20140820;第38卷(第6期);全文 * |
"计及负荷静态特性的电力系统静态状态估计算法";陆子刚,卫志农,孙国强;《河海大学学报(自然科学版)》;20140217(第1期);全文 * |
"计及零注入约束的电力系统动态状态估计";杨韵;《广东电力》;20160922;第29卷(第9期);全文 * |
Esmaeil Ghahremani;Innocent Kamwa."Local and Wide-Area PMU-Based Decentralized Dynamic State Estimation in Multi-Machine Power Systems".《IEEE Transactions on Power Systems》.2016,第31卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN106936628A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Variational Bayesian approach for ARX systems with missing observations and varying time-delays | |
CN107590317B (zh) | 一种计及模型参数不确定性的发电机动态估计方法 | |
Sayed | Diffusion adaptation over networks | |
US9627886B2 (en) | State estimation for power system using hybrid measurements | |
CN108319570B (zh) | 一种异步多传感器空时偏差联合估计与补偿方法及装置 | |
CN108448585B (zh) | 一种基于数据驱动的电网潮流方程线性化求解方法 | |
Chen et al. | Impulse response estimation with binary measurements: A regularized FIR model approach | |
Moreau et al. | Estimation of power spectral density from laser Doppler data via linear interpolation and deconvolution | |
CN106972949B (zh) | 一种基于自适应补偿技术的分数阶网络系统状态估计方法 | |
CN104076332A (zh) | 一种雷达均匀线性阵列幅度和相位的估计方法 | |
Fan | Data fusion-based distributed prony analysis | |
CN106707234B (zh) | 一种联合时延差与角度测量的传感器网络目标定位方法 | |
Qi et al. | Estimation of distribution function for control valve stiction estimation | |
CN106936628B (zh) | 一种计及传感器故障的分数阶网络系统状态估计方法 | |
CN114239796A (zh) | 一种基于扩展卡尔曼滤波的电力系统状态估计方法 | |
Zhang et al. | Robust adaptive Unscented Kalman Filter with gross error detection and identification for power system forecasting-aided state estimation | |
Mårtensson et al. | Covariance analysis in SISO linear systems identification | |
CN106878076B (zh) | 计及数据丢包和增益扰动的分数阶网络系统状态估计方法 | |
CN109638811B (zh) | 基于模型等值的配电网电压功率灵敏度鲁棒估计方法 | |
CN106202926A (zh) | 基于多节点协同探测的空间系统偏差配准优化方法 | |
Khan et al. | Tracking of wireless mobile nodes in he presence of unknown path-loss characteristics | |
Biswas et al. | A model-agnostic method for PMU data recovery using optimal singular value thresholding | |
CN116088303A (zh) | 一种不确定复杂动态网络状态时变递推估计方法 | |
Eappen et al. | Semiparametric algorithm for processing MST radar data | |
Galvão et al. | Identification of fractional-order transfer functions using exponentially modulated signals with arbitrary excitation waveforms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |