CN105426583B - 一种基于同步的同质传感器融合处理方法 - Google Patents

一种基于同步的同质传感器融合处理方法 Download PDF

Info

Publication number
CN105426583B
CN105426583B CN201510742237.1A CN201510742237A CN105426583B CN 105426583 B CN105426583 B CN 105426583B CN 201510742237 A CN201510742237 A CN 201510742237A CN 105426583 B CN105426583 B CN 105426583B
Authority
CN
China
Prior art keywords
data
estimation
period
sensor
sample
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
Application number
CN201510742237.1A
Other languages
English (en)
Other versions
CN105426583A (zh
Inventor
曾伟
叶远誉
江峰
范瑞祥
曹蓓
黄明山
舒志猛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
STATE GRID JIANGXI ELECTRIC POWER Co
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Henan Xuji Instrument Co Ltd
Original Assignee
STATE GRID JIANGXI ELECTRIC POWER Co
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Henan Xuji Instrument Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by STATE GRID JIANGXI ELECTRIC POWER Co, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd, Henan Xuji Instrument Co Ltd filed Critical STATE GRID JIANGXI ELECTRIC POWER Co
Priority to CN201510742237.1A priority Critical patent/CN105426583B/zh
Publication of CN105426583A publication Critical patent/CN105426583A/zh
Application granted granted Critical
Publication of CN105426583B publication Critical patent/CN105426583B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

一种基于同步的同质传感器融合处理方法，中央处理器以T秒为采样周期，同质传感器定时对被测信号进行采样和量化并获得关于时间的数据序列，将采样数据按周期进行截取，获得m个周期序列为：，每个周期样本中包含了N个数据点，即，其中。从每个周期序列各时刻点中抽取一个点值组成单一传感器多次测量的一维数据序列，构造周期数据序列，即，每组一维数据序列中包括了m个数据点，即。对得到的数据样本进行基于极大似然估计和最小二乘法估计的优化组合预处理，用于指导电力系统周期采样。本发明将不同类传感器进行无差别融合，得到等周期数据样本，以消减不同传感器差异和信号采集过程的随机误差，该方法适用于电力系统数据预处理，改善数据质量。

Description

一种基于同步的同质传感器融合处理方法
技术领域
本发明涉及一种基于同步的同质传感器融合处理方法，属电力设备传感器测量技术领域。
背景技术
传统的用电设备能实时地感知电网的有效供电能力是有序用电的基础,电网能实时获取用电设备用电真实状态是智能电网的前提。传统的供用电信息采集往往由单一的传感器来完成，即使采用多个(种)传感器也多是分时使用，所以是从多个侧面孤立地反映电网的信息。随着技术的进步，这些测量数据需要融合处理，即指利用多个传感器的输出推断出一个有效的信息。
电力系统所使用的传感器受到生产厂家生产环境、采用技术的不同。不同厂家的仪器甚至同一厂家的不同批次的仪器测得同一对象的某个生理数据都有差异，尤其在幅度上差别明显，因此对海量数据对比、自动分析等造成困难。
发明内容
本发明的目的是，针对电力设备传感器采集信息存在的问题，本发明公开了一种基于同步的同质传感器融合处理方法。
实现本发明的技术方案是，一种基于同步的同质传感器融合处理方法，将电力系统不同类传感器进行无差别融合，进行优化组合预处理，得到等周期数据样本；以消减不同传感器差异和信号采集过程中的随机误差，改善数据质量。
本发明提供了一种基于同步的同质传感器融合处理方法，中央处理器以T秒为采样周期，同质传感器(传感器观测的是同一物理现象)定时对被测信号进行采样和量化并获得关于时间的数据序列X(N)，将采样数据按周期进行截取，获得m个周期序列为：Y1(t),Y2(t),…Ym(t)，每个周期样本中包含了N个数据点，即Yi＝[Xi(1),Xi(2),…Xi(N)]，其中i＝1,2,…,m。从每个周期序列各时刻点中抽取一个点值组成单一传感器多次测量的一维数据序列，构造按周期同步构建的周期数据序列，即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]。对得到的数据样本Y′t进行基于极大似然估计和最小二乘法估计的优化组合预处理，用于指导电力系统周期采样数据指导。
一种基于同步的同质传感器融合处理方法，包括以下步骤：
步骤一：中央处理器以T(秒)为采样周期，同质传感器定时对系统被测信号进行采样和量化，并得到相同采样频率下的数据样本X(N)；
步骤二：将采样数据按周期进行截取，获得m个周期序列为：Y1(t),Y2(t),…Ym(t)，每个周期样本中包含了N个数据点，即Yi＝[Xi(1),Xi(2),…Xi(N)]，其中i＝1,2,…,m；
步骤三：构造按周期同步构建的周期数据序列，即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]；
步骤四：对数据样本Y′t进行基于极大似然估计和最小二乘法估计的优化组合预处理。
步骤一中，所述同质传感器定时对系统被测信号进行采样和量化的具体步骤为：
所述中央处理器采样周期为T(秒)，包含在系统内部的m个传感器定时对系统内部的被测信号进行采样和量化，进而得到相同采样频率下的数据样本Xi(N)，其中i＝1,2,…,m；所述m个传感器可以不同时、不同位置，但是被检测或采集信号的特征相同。
步骤三中，所述构造按周期同步构建的周期数据序列的具体步骤为：
将采集到的电力系统数据在按周期进行截取后，将每个周期数据中各个时刻点数据视为同步后的同一对象的多次测量结果，即各个周期对应时刻点数据为一个样本，多个周期的数据中相对应时刻的每个点值即构成一个数据序列，从而可视为单一传感器多次测量的一维数据序列，即构造了按周期同步构建的周期数据序列，也即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]。
步骤四中，所述对数据样本Y′t进行基于极大似然估计的优化组合预处理的步骤为：
(a)假设在给定的时刻，待测环境特征为X，传感器的值为Y,则该传感器的测量模型为：Y＝f(X)+V，其中，V是符合高斯分布的噪声项；所谓数据融合就是由N个传感器得到测量值Y1、Y2、…、YN,并按某种估计准则从这些测量值中得到特征参数X的最优估计；
(b)寻找合适的准则函数，即当X被估计为X(Y)时所产生的损失最小的函数；取损失函数为均匀损失：
(c)在损失函数L的基础上，定义相应估计风险的函数R：
其中，p(x)、p(x|y)表示概率分 布；
(d)取风险最小为估计准则，即
其中，为求导操作；可以得到符合式(1)的最优估计(最大后验估计)为：
(e)在具有N个传感器的系统中，相应的信息融合可以看作是在观测值Y1、Y2、…、YN下，值X具有最大后验的估计为
(f)取(a)中f(X)为X的线性函数，即Yi＝AiXi+Vi，Ai为权重系数，Vi为噪声；p(Y|X)服从高斯噪声分布，即(Y|X)～N(AiX,Vi)，N(·)表示高斯分布函数；X满足高斯噪声分布其中为均值，VX为方差；
采用“模糊先验”的概念，即对所有可能参数X均采用p(x)＝1,得到：
其中，p(·)表示概率，∏(·)表示联合分布律，此时最大后验估计即简化为极大似然估计，相应的融合计算公式为：
其中，Vf为X|Y分布的协方差，(·)T表示矩阵转置操作，(·)-1表示求逆矩阵操作；
所述对数据样本Yi'进行基于最小二乘法估计的优化组合预处理按以下列步骤进行：
对最大似然估计信息融合的计算以后,再考虑最小二乘估计,相应的融合计算公式为：
其中，为参数估计值。
本发明的有益效果是，本发明将不同类传感器进行无差别融合，进行优化组合预处理，得到等周期数据样本，以消减不同传感器差异和信号采集过程中的随机误差，该方法适用于电力系统数据预处理，改善数据质量。
附图说明
图1是本发明实施例一种基于同步的同质传感器融合处理方法的流程图。
具体实施方式
本发明的具体实施方式如图1所示。
本实施例一种基于同步的同质传感器融合处理方法，其主要包括以下步骤：
步骤A：中央处理器以T(秒)为采样周期，同质传感器定时对系统被测信号进行采样和量化，并得到相同采样频率下的数据样本X(N)；
步骤B：将采样数据按周期进行截取，获得m个周期序列为：Y1(t),Y2(t),…Ym(t)，每个周期样本中包含了N个数据点，即Yi＝[Xi(1),Xi(2),…Xi(N)]，其中i＝1,2,…,m；
步骤C：构造按周期同步构建的周期数据序列，即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]；
步骤D：对数据样本Y′t进行基于极大似然估计和最小二乘法估计的优化组合预处理。
本实施例所属步骤A：中央处理器以T(秒)为采样周期，同质传感器定时对系统被测信号进行采样和量化，并得到相同采样频率下的数据样本X(N)的具体步骤为：
中央处理器采样周期为T(秒)，包含在系统内部的m个传感器(可以不同时、不同位置，但是被检测或采集信号的特征相同)定时对系统内部的被测信号进行采样和量化，进而得到相同采样频率下的数据样本Xi(N)，其中i＝1,2,…,m。
本实施例所属步骤C：构造按周期同步构建的周期数据序列，即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]的具体步骤为：
将采集到的电力系统数据在按周期进行截取后，将每个周期数据中各个时刻点数据视为同步后的同一对象的多次测量结果，即各个周期对应时刻点数据为一个样本，多个周期的数据中相对应时刻的每个点值即构成一个数据序列，从而可视为单一传感器多次测量的一维数据序列，即构造了按周期同步构建的周期数据序列，也即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]。
本实施例所属步骤D：对数据样本Y′t进行基于极大似然估计和最小二乘法估计的优化组合预处理的具体步骤按照两种方案进行：
方案一：对数据样本Yt'进行基于极大似然估计的优化组合预处理以下列步骤进行：
S101：假设在给定的时刻，待测环境特征为X，传感器的值为Y,则该传感器的测量模型为：Y＝f(X)+V其中，V是符合高斯分布的噪声项。所谓数据融合就是由N个传感器得到测量值Y1、Y2、…、YN,并按某种估计准则从这些测量值中得到特征参数X的最优估计。
S102：寻找合适的准则函数，即当X被估计为X(Y)时所产生的损失最小。取损失函数为均匀损失：
S103：在损失函数L的基础上，定义相应估计风险的函数R：
其中，p(x)、p(x|y)表示概率分 布；
S104：取风险最小为估计准则，即
其中，为求导操作，可以得到符合式(1)的最优估计(最大后验估计)为:
S105：在具有N个传感器的系统中，相应的信息融合可以看作是在观测值Y1、Y2、…、YN下，值X具有最大后验的估计为
S106：取S101中f(X)为X的线性函数，即Yi＝AiXi+Vi，Ai为权重系数，Vi为噪声。p(Y|X)服从高斯噪声分布，即(Y|X)～N(AiX,Vi)，N(·)表示高斯分布函数；X满足高斯噪声分布其中为均值，VX为方差。采用“模糊先验”的概念，即对所有可能参数X均采用p(x)＝1,得到：
其中，p(·)表示概率，∏(·)表示联合分布律，此时最大后验估计即简化为极大似然估计，相应的融合计算公式为：
其中，Vf为X|Y分布的协方差，(·)T表示矩阵转置操作，(·)-1表示求逆矩阵操作。
方案二：对数据样本Y′t进行基于最小二乘法估计的优化组合预处理以下列步骤进行：
S201：依据方案一中关于对最大似然估计信息融合的计算以后,再考虑最小二乘估计，相应的融合计算公式为：
其中，为参数估计值。
本发明实施例公开了一种基于同步的同质传感器融合处理的方法，中央处理器以T秒为采样周期，同质传感器(传感器观测的是同一物理现象)定时对被测信号进行采样和量化并获得关于时间的数据序列X(N)，将采样数据按周期进行截取，获得m个周期序列为：Y1(t),Y2(t),…Ym(t)，每个周期样本中包含了N个数据点，即Yi＝[Xi(1),Xi(2),…Xi(N)]，其中i＝1,2,…,m。从每个周期序列各时刻点中抽取一个点值组成单一传感器多次测量的一维数据序列，构造按周期同步构建的周期数据序列，即Y′1(t),Y′2(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Y′i＝[X1(1),X2(2),…Xm(N)]。对得到的数据样本Y′t进行基于极大似然估计和最小二乘法估计的优化组合预处理，用于指导电力系统周期采样数据指导。本发明将不同类传感器进行无差别融合，得到等周期数据样本，以消减不同传感器差异和信号采集过程中的随机误差，该方法适用于电力系统数据预处理，改善数据质量。

Claims (3)

1.一种基于同步的同质传感器融合处理方法，其特征在于，所述方法将电力系统不同类传感器进行无差别融合，进行优化组合预处理，得到等周期数据样本；以消减不同传感器差异和信号采集过程中的随机误差，改善数据质量；
所述方法包括以下步骤：
步骤一：中央处理器以T秒为采样周期，同质传感器定时对系统被测信号进行采样和量化，并得到相同采样频率下的数据样本X(N)；
步骤二：将采样数据按周期进行截取，获得m个周期序列为：Y1(t),Y2(t),…Ym(t)，每个周期样本中包含了N个数据点，即Yi＝[Xi(1),Xi(2),…Xi(N)]，其中i＝1,2,…,m；
步骤三：构造按周期同步构建的周期数据序列，即Y1′(t),Y2′(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Yi′＝[X1(1),X2(2),…Xm(N)]；
步骤四：对数据样本Yi'进行基于极大似然估计和最小二乘法估计的优化组合预处理；
所述同质传感器定时对系统被测信号进行采样和量化的步骤如下：
所述中央处理器采样周期为T秒，包含在系统内部的m个传感器定时对系统内部的被测信号进行采样和量化，进而得到相同采样频率下的数据样本Xi(N)，其中i＝1,2,…,m；所述m个传感器可以不同时、不同位置，但是被检测或采集信号的特征相同；
所述构造按周期同步构建的周期数据序列的具体为：
将采集到的电力系统数据在按周期进行截取后，将每个周期数据中各个时刻点数据视为同步后的同一对象的多次测量结果，即各个周期对应时刻点数据为一个样本，多个周期的数据中相对应时刻的每个点值即构成一个数据序列，从而可视为单一传感器多次测量的一维数据序列，即构造了按周期同步构建的周期数据序列，也即Y1′(t),Y2′(t),…Y′N(t)，每组一维数据序列中包括了m个数据点，即Yi′＝[X1(1),X2(2),…Xm(N)]。
2.根据权利要求1所述一种基于同步的同质传感器融合处理方法，其特征在于，所述对数据样本Yi'进行基于极大似然估计的优化组合预处理的步骤为：
(a)假设在给定的时刻，待测环境特征为X，传感器的值为Y,则该传感器的测量模型为：Y＝f(X)+V，其中，V是符合高斯分布的噪声项；所谓数据融合就是由N个传感器得到测量值Y1、Y2、…、YN,并按某种估计准则从这些测量值中得到特征参数X的最优估计；
(b)寻找合适的准则函数，即当X被估计为X(Y)时所产生的损失最小的函数；取损失函数为均匀损失：
(c)在损失函数L的基础上，定义相应估计风险的函数R：
其中，p(x)、p(x|y)表示概率分布；
(d)取风险最小为估计准则，即
其中，为求导操作；可以得到符合式(1)的最大后验估计为：
(e)在具有N个传感器的系统中，相应的信息融合可以看作是在观测值Y1、Y2、…、YN下，值X具有最大后验的估计为
(f)取(a)中f(X)为X的线性函数，即Yi＝AiXi+Vi，Ai为权重系数，Vi为噪声；p(Y|X)服从高斯噪声分布，即(Y|X)～N(AiX,Vi)，N(·)表示高斯分布函数；X满足高斯噪声分布其中为均值，VX为方差；
采用“模糊先验”的概念，即对所有可能参数X均采用p(x)＝1,得到：
其中，p(·)表示概率，П(·)表示联合分布律，此时最大后验估计即简化为极大似然估计，相应的融合计算公式为：
其中，Vf为X|Y分布的协方差，(·)T表示矩阵转置操作，(·)-1表示求逆矩阵操作。
3.根据权利要求1所述一种基于同步的同质传感器融合处理方法，其特征在于，所述对数据样本Yi'进行基于最小二乘法估计的优化组合预处理按以下列步骤进行：
对最大似然估计信息融合的计算以后,再考虑最小二乘估计,相应的融合计算公式为：
其中，为参数估计值。
CN201510742237.1A 2015-11-03 2015-11-03 一种基于同步的同质传感器融合处理方法 Active CN105426583B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510742237.1A CN105426583B (zh) 2015-11-03 2015-11-03 一种基于同步的同质传感器融合处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510742237.1A CN105426583B (zh) 2015-11-03 2015-11-03 一种基于同步的同质传感器融合处理方法

Publications (2)

Publication Number Publication Date
CN105426583A CN105426583A (zh) 2016-03-23
CN105426583B true CN105426583B (zh) 2019-11-26

Family

ID=55504792

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510742237.1A Active CN105426583B (zh) 2015-11-03 2015-11-03 一种基于同步的同质传感器融合处理方法

Country Status (1)

Country Link
CN (1) CN105426583B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105954695B (zh) * 2016-04-20 2020-04-17 国家电网公司 一种基于同步的同质传感器突变参数识别方法与装置
CN106019353B (zh) * 2016-07-29 2020-01-17 山东省济南生态环境监测中心 一种环境辐射检测方法
JP6348163B2 (ja) * 2016-12-15 2018-06-27 ファナック株式会社 制御装置及び制御システム
CN108334932B (zh) * 2017-11-27 2022-03-29 中科观世（北京）科技有限公司 基于目标信号特征的频率区分方法
US11287530B2 (en) * 2019-09-05 2022-03-29 ThorDrive Co., Ltd Data processing system and method for fusion of multiple heterogeneous sensors
CN114739448A (zh) * 2022-03-09 2022-07-12 中铁第四勘察设计院集团有限公司 一种数据处理方法、系统、装置和存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565778A (zh) * 2011-12-12 2012-07-11 中国科学院遥感应用研究所 一种自动提取伪不变特征的遥感图像相对辐射校正方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9450277B2 (en) * 2013-09-13 2016-09-20 Recycling Coordinators, Inc. Systems for recycling volatile battery sources

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565778A (zh) * 2011-12-12 2012-07-11 中国科学院遥感应用研究所 一种自动提取伪不变特征的遥感图像相对辐射校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于协整的多传感器测试仿真模型的研究;黎亮等;《系统仿真学报》;20140708;第26卷(第7期);第1559-1565页 *

Also Published As

Publication number Publication date
CN105426583A (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
CN105426583B (zh) 一种基于同步的同质传感器融合处理方法
KR101958674B1 (ko) 시퀀스 재귀 필터링 3차원 변분(3d-var) 기반의 실측 해양 환경 데이터 동화방법
CN109685314B (zh) 一种基于长短期记忆网络的非侵入负荷分解方法和系统
CN107247259B (zh) 基于神经网络的k分布海杂波形状参数估计方法
CN110048416B (zh) S-g滤波和自适应mp算法的低频振荡模态辨识方法
CN106419893A (zh) 睡眠状态检测方法和装置
CN105954695B (zh) 一种基于同步的同质传感器突变参数识别方法与装置
CN103020478A (zh) 一种海洋水色遥感产品真实性检验的方法
CN109815855B (zh) 一种基于机器学习的电子设备自动测试方法及系统
RU2016142190A (ru) Коррекция посторонних эхосигналов epi
CN103675544A (zh) 基于优化算法的电力系统故障信号检测与波形识别方法
CN103278867A (zh) 一种基于经验模态分解的自动气象站数据质量控制方法
CN107192411A (zh) 一种航天器多遥测参数分析窗口划分方法
CN115718906A (zh) 一种多能源系统多源异构数据融合方法及系统
CN108931776A (zh) 一种高精度的匹配场定位方法
CN104320144A (zh) 稀疏度自适应信号重构方法
CN107525686B (zh) 一种航天器多遥测参数分析窗口划分方法
CN104002019A (zh) 基于电弧混沌特性优化算法的焊接材料工艺性评价方法
CN102750547B (zh) 一种基于压缩传感的水果大小分级方法
CN107884752A (zh) 一种基于压缩感知的外辐射源雷达对目标检测方法
CN108170951A (zh) 基于采样数据时空匹配示踪试验的纵向离散系数确定方法
KR101608155B1 (ko) 임계치 파형 작성 장치
CN107219347B (zh) 一种用于获取土壤湿度的传感器布置优化方法
CN109670243B (zh) 一种基于勒贝格空间模型的寿命预测方法
CN104111272A (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
GR01 Patent grant
GR01 Patent grant