CN101614554A - 连续物理量测量装置及方法 - Google Patents

连续物理量测量装置及方法 Download PDF

Info

Publication number
CN101614554A
CN101614554A CN200910158370A CN200910158370A CN101614554A CN 101614554 A CN101614554 A CN 101614554A CN 200910158370 A CN200910158370 A CN 200910158370A CN 200910158370 A CN200910158370 A CN 200910158370A CN 101614554 A CN101614554 A CN 101614554A
Authority
CN
China
Prior art keywords
steady
state process
value
sampled value
state
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
Application number
CN200910158370A
Other languages
English (en)
Other versions
CN101614554B (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.)
Sanchuan Power Equipment Co ltd
Original Assignee
Baoding Sanchuan Electric 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 Baoding Sanchuan Electric Co Ltd filed Critical Baoding Sanchuan Electric Co Ltd
Priority to CN2009101583707A priority Critical patent/CN101614554B/zh
Publication of CN101614554A publication Critical patent/CN101614554A/zh
Application granted granted Critical
Publication of CN101614554B publication Critical patent/CN101614554B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开了一种连续物理量测量装置及方法,包括用来输出时基和时标的定时单元;采样单元,在所述时基控制下以满足奈奎斯特采样定理的时间间隔对连续物理量进行采样,并输出采样值;判断单元,判断采样值是否进入或者离开稳态过程;计算单元,在稳态过程中对所述采样单元的采样值计算稳态值;输出单元,根据判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及所述稳态值。本发明可以在不产生混迭的同时减少输出数据量。

Description

连续物理量测量装置及方法
技术领域
本发明涉及连续物理量的测量和记录技术,尤其是涉及热工、电力系统中的温度、压力、流量、电压、电流、功率、相角等的数字测量和记录。
背景技术
连续物理量的值在相当长一段时间维持基本不变或变化不大,称为稳态过程;两个稳态过程之间的短时间的过渡称为暂态过程。对于稳态过程和暂态过程组成的过程物理量采样,按照奈奎斯特(也称作香农)采样定理,采样间隔要小于其暂态过程的最小时间常数τmin的一半。按照这样的采样间隔采集数据会产生海量记录数据。在遥测情况下,更要占用大量通信资源传送这些数据。因此不便记录和传送。
为了减少记录和传送的数据量,现有的办法是人为规定一个采样时间间隔T,只记录或传送该间隔整倍数点上的数据,如图1所示。显然,T不满足采样定理,必定造成混迭,致使记录或传送的数据不能正确反映物理过程的状态和变化,数据错误而无法使用。
如果取T时间段内的测量平均值替代瞬时值,也不能解决采集数据的错误,除了产生混迭外,还会产生额外误差,使本来积等于两数相乘的三个量,其平均值不再存在积的关系。例如,变量A、B、C,其中,A=B·C,B=B+Δb,C=C+Δc,A≠B·C。
发明内容
本发明的目的旨在至少解决现有技术中的上述问题之一。
为此,本发明的实施例提出一种既减少输出或传送的数据量,又不产生混迭的、数据测量准确的连续物理量测量装置和方法。
根据本发明的一个方面,本发明的实施例提出了一种连续物理量测量装置,包括用来输出时基和时标的定时单元;采样单元,在所述时基控制下以满足奈奎斯特采样定理的时间间隔对连续物理量进行采样,并输出采样值;判断单元,判断采样值是否进入或者离开稳态过程;计算单元,在稳态过程中对所述采样单元的采样值计算稳态值;输出单元,根据判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及所述稳态值。
根据本发明进一步的实施例,所述判断单元根据统计法的t分布或者滤波器输出进行判断。
根据本发明进一步的实施例,还包括预处理单元,用于对所述采样单元输出的采样值进行标度变换、坏数据去除、重抽样和/或有效值计算。
根据本发明进一步的实施例,所述计算单元根据均值算法或低通滤波算法计算所述稳态值。
根据本发明进一步的实施例,还包括记录单元,所述记录单元对所述稳态过程的开始时间、结束时间以及在所述结束时间对应的稳态值进行记录。
根据本发明的另一方面,本发明的实施例提出一种连续物理量测量方法,所述测量方法包括以下步骤:在定时器输出时基的控制下以满足奈奎斯特采样定理的时间间隔对连续物理量进行采样,并输出采样值xk和时标;判断采样值xk是否进入或者离开稳态过程;以及根据判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及在稳态过程中对采样值计算获得的稳态值X。
根据本发明进一步的实施例,所述稳态值X根据均值算法或低通滤波算法计算获得。
根据本发明进一步的实施例,所述判断步骤包括:计算采样值xk对应的均值xk和方差
Figure G2009101583707D00021
判断
Figure G2009101583707D00022
是否服从t分布;若服从,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;若不服从,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。所述t分布判断可以依据以下公式: | x k - x ‾ k | ≤ k · t α / 2 ( k - 1 ) · s ^ k , 其中α为风险系数,一般取值0.05左右。
根据本发明进一步的实施例,所述判断步骤包括:计算采样值xk对应的均值xk和方差
Figure G2009101583707D00031
判断是否满足 | x k - x ‾ k | ≤ A · s ^ k , 其中A为给定值,A在3~10之间;若满足,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;若不满足,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。
根据本发明进一步的实施例,所述判断步骤包括:计算采样值xk对应的均值xk;判断是否满足|xk-xk|≤δ·xn,其中δ为给定值,一般在2%到10%之间,xn为所述连续物理量对应的额定值;若满足,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;若不满足,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。
根据本发明进一步的实施例,所述判断步骤包括:对采样值xk进行滤波,分别得到采样值xk的滤波分量;分别判断所述滤波分量是否超出对应设定的上限值;以及若均未超出,判断采样值xk未离开稳态过程;若存在一个分量超出,判断采样值xk离开稳态过程。
根据本发明进一步的实施例,所述判断步骤包括:对采样值xk进行滤波,分别得到采样值xk的滤波分量;分别判断所述滤波分量是否超出对应设定的下限值;以及若存在一个分量超出,判断采样值xk未进入稳态过程;若均未超出,判断采样值xk进入稳态过程。
根据本发明进一步的实施例,还包括对采样值xk进行标度变换、坏数据剔除、重抽样和/或有效值计算的预处理步骤。
根据本发明进一步的实施例,还包括对所述稳态过程的开始时间、结束时间以及在所述结束时间对应的稳态值进行记录的步骤。
本发明可以将输出或传送的数据大幅度减少,同时可以避免数据出现混迭,提高数据测量的准确性,保证了信号不失真地被输出或记录下来。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为现有连续物理量的采样输出示意图;
图2为本发明实施例的连续物理量测量装置结构示意图;
图3为本发明实施例的稳态过程中数据及时标的测量示意图;
图4显示了本发明实施例的连续物理量测量方法流程图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
现在参考图2,图2为本发明实施例的连续物理量测量装置结构示意图。如图2中所示,测量装置包括定时单元20、采样单元12、判断单元14、计算单元16和输出单元18。定时单元20输出时基(定时的基准信号)和时标(一般可以理解为年月日时分秒),用于采样单元12的定时。采样单元12可以在定时单元20的输出时基控制下,按照时间间隔Δt对连续物理量进行采样,其中Δt满足奈奎斯特采样定理。
为了简化正弦量的计算,按照等相位Δφ采样(如锁相环输出控制的采样),Δφ与Δt存在关系,所以,认为是上述采样的变种,而这种变种也是普遍公知的。
采样单元12将不断获得的采样值xk输出给判断单元14,判断单元14则判断接收到的每个采样值xk是否进入或者离开稳态过程。
在一个实施例中,判断单元14可以根据统计法的t分布判断采样值xk是否进入稳态过程或者离开稳态过程。
如果为平稳随机过程,由于随机因素的影响,落在3σ内概率在99.7%之内;如果x服从正态分布,则
Figure G2009101583707D00041
服从t分布。(参考浙江大学的《概率论与数理统计》高等教育版和茆诗书的《概率论与数理统计》中国统计出版社第2版)。在测量中,递推X和
Figure G2009101583707D00042
对于第k次,就能够判断数据xk是否仍属于稳态还是不属于。
在本发明的实施例中,判断单元14对于采样单元12输出的采样值xk,可以递推得到其对应的均值xk和方差
Figure G2009101583707D00051
因此通过判断是否服从t分布,从而可以判断采样值xk是否离开稳态过程或者进入稳态过程。
具体来说,若判断
Figure G2009101583707D00053
服从t分布,则采样值xk未离开当前稳态过程或者采样值xk进入稳态过程;若不服从,判断采样值xk离开当前稳态过程或者采样值xk未进入稳态过程。
当然,判断单元14也可以采用高斯分布等其他统计法对采样值xk的稳态状态进行判断,这里不局限于上述具体实施例。
在一个实施例中,判断单元14可以根据滤波器输出判断采样值是否进入或者离开稳态过程。这里,滤波器可以是任意1阶或更高阶(例如2阶、3阶等)的滤波器,例如n阶的卡尔曼滤波器。
判断单元14通过对采样值xk进行滤波,从而可以得到采样值xk对应的滤波分量。显然,对于不同阶的滤波算法,得到的滤波分量个数也不同。这里,本发明不局限于具体个数的滤波分量,任意合适的滤波算法均可落在本发明的保护范围内。
判断单元14分别判断所得的滤波分量是否超出对应设定的上限值,若均未超出,判断采样值xk未离开稳态过程;若存在一个分量超出,判断采样值xk离开稳态过程。
类似地,通过分别判断上述滤波分量是否超出对应设定的下限值,则可以判断采样值xk是否进入稳态过程。具体来说,若存在一个分量超出,判断采样值xk未进入稳态过程;若均未超出,判断采样值xk进入稳态过程。
当然,本发明可以采用其他合适方法来判断采样值xk是否进入或者离开稳态过程。
这样,判断单元14判断新数据是否仍处于相对稳定的数值附近,还是出现了超出稳态过程的突变。
在判断采样值未离开稳态过程时,判断单元14将采样值xk输入给计算单元16。计算单元16用于在稳态过程中,对采样单元12的采样值进行稳态值计算。对于连续物理量,在稳态过程能够以一个值X,例如平均值代表其稳态值大小。对于离开稳态过程的采样值,判断单元14可以不输入到计算单元16中。
计算单元16可以根据均值算法或低通滤波算法计算稳态过程对应输出采样值的稳态值。稳态值X通过稳态时间段内的均值算法或低通滤波算法获得,以避免随机干扰、提高测量的精度。
在一个实施例中,稳态值X对应的均值算法公式例如表示为 x ‾ k = 1 k [ ( k - 1 ) · x ‾ k - 1 + x k ) ] , 其中xk表示稳态过程获得的采样值,k表示这些采样值的个数,xk即为根据这些采样值计算得到的稳态值X。k=1时,则x1=x1
在一个实施例中,稳态值X可以根据低通滤波器的递推算法计算获得,例如利用公式表示为xk=α·xk-1+(1-α)·xk的一阶低通滤波算法,其中α为设定的常数,与物理量的时间常数有关。
当然,本领域技术人员显然可知,本发明计算单元16还可以采用加权均值法或更高阶(例如2阶、3阶等)滤波器来获得稳态值X。本发明不局限于上述具体实施例。
对于上述稳态值X计算的递推算法,计算单元16利用新的采样值xk对当前稳态值X进行修正,得到对应稳态过程的新稳态值X。这里,对X修正与之前X的计算相同,利用当前稳态过程时间段内的采样数据获得,只不过这里进一步包括了新数据xk。输出单元18则根据判断单元14的判断结果,并结合定时单元20输出的时标,输出每个稳态过程对应的开始和结束时间,以及每个稳态过程中对应的稳态值X。
具体来说,对于连续物理量,假定判断单元14判断其以往位于稳态过程,计算单元14以一个值X代表该稳态过程中对应已采样数据的稳态值大小。在采样单元12以时间间隔采样之后,获得一个新采样数据xk;判断单元14判断xk是否仍在稳态过程,是,则由计算单元16利用xk修正X,继续稳态过程;不是,则由输出单元18输出该稳态过程的结束时间以及该稳态过程的稳态值X。
另外,如果以往不在稳态过程,对于新数据xk判断单元14判断其是否进入到一个稳态过程;若没有进入,则继续当前状态;否则,由输出单元18输出稳态过程的开始时标,并由计算单元16将xk作为该稳态过程的初始稳态数据x1,建立该稳态过程对应的新的稳态值X;如此继续。
在判断单元14利用统计法的t分布进行判断的实施例中,判断单元14所需的均值xk的算法,可以与计算单元16计算稳态值X的算法相同。这样,判断单元14可以将采样单元12的采样值xk输出给计算单元16,由计算单元16递推计算得到对应的均值和方差,并按图示虚线箭头返给判断单元14。
若判断单元14根据计算单元16返回的数值结合具体t分布判断采样值xk仍在稳态过程或进入稳态过程,则将计算单元16计算的均值作为该稳态过程的当前稳态值,并输出稳态值X给输出单元18。
对于输出单元18,其可以根据判断单元14的判断,仅输出在稳态过程结束时对应的稳态值X,也可以根据远程需要,将稳态过程中不同时间对应的稳态值X输出。
在一个实施例中,本发明的测量装置还包括可以记录单元(图中未显示),记录单元对每个稳态过程的开始时间、结束时间以及在结束时间对应的稳态值进行记录。
例如对于图3所示采样的连续物理量,按照上述方式,输出单元18对其稳态过程输出的数据可以是:
T0,{X1},T1,T2,{X2},T3,...,
其中,{X1}、{X2}、...为各个稳态过程对应的稳态值,如上文所示,每个稳态过程输出的稳态值可以是不同时间对应的稳态值集合。T0(t=0)、T1、T2、...为表示对应稳态过程开始时间和结束时间的时标。
对于记录单元,则其可以记录数据如下:
T0,X1,T1,T2,X2,T3,...,记为{Dw}-稳态记录。这里,X1、X2是对应每个稳态过程结束时间的稳态值。
一般地,以上输出和记录的稳态数据比现有的间隔T时间输出一次数据获得的数据量少(多个T都在一个稳态过程)且避免了混迭错误。并且,它们真实反映了系统在该时间段上的特性,从而可作为后续系统计算的无误输入。
由于稳态时间较暂态时间要长得多,所以输出或传送的数据大幅度减少。另外,又保证了信号不失真地被记录下来。当然,由记录的数据也能够复现原信号。
在一个实施例中,测量装置还可以包括预处理单元(图中未显示),预处理单元可以对采样单元12输出的采样值进行预处理,然后直接或重抽样输出到计算单元14。一般地,预处理包括标度变换,还可以包括坏数据剔除,而对于正弦量,例如交流电压、交流电流、交流功率等,还可以包括有效值的计算,依据应用情况,有效值还可以重抽样。在这种情况下,预处理单元输出的是重抽样后的有效值。
现在参考图4,该图显示了本发明实施例的连续物理量测量方法流程图。
首先,可以在定时器输出时基的控制下,按照预定时间间隔Δt对连续物理量进行采样,并分别输出采样值和时标(步骤102),其中Δt满足奈奎斯特采样定理。然后,判断采样值是否进入或者离开稳态过程(步骤104),根据对应的判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及在稳态过程中对采样值计算获得的稳态值(步骤106)。定时器在输出每个时基时会相应地给出对应的时标,因此根据时基可确定每个采样值对应的时标。
在一个实施例中,步骤104可以根据统计法的t分布判断采样值是否进入稳态过程或者离开稳态过程。
对于统计法的t分布:在方差分析中,如果为平稳随机过程,由于随机因素的影响,落在3σ内概率在99.7%之内。现实中,若σ未知,可以用
Figure G2009101583707D00081
代替,且
Figure G2009101583707D00082
服从t分布,在测量中,递推X和
Figure G2009101583707D00083
对于第k次采样,就能够判断数据xk是否仍属于稳态还是不属于。
在本发明的实施例中,对于采样输出的采样值xk,可以递推得到其对应的均值xk和方差
Figure G2009101583707D00084
例如,根据均值公式: x ‾ k = 1 k [ ( k - 1 ) · x ‾ k - 1 + x k ) ] ,
以及方差公式: s ^ k 2 = k - 2 k - 1 s ^ k - 1 2 + ( x ‾ k - x ‾ k - 1 ) 2 + 1 k - 1 ( x ‾ k - x k ) 2 获得对应的均值和方差,其中k为稳态过程或者离开稳态过程后对应的采样值的数量。
因此通过判断
Figure G2009101583707D00087
是否服从t分布,从而可以判断采样值xk是否离开稳态过程或者进入稳态过程。
具体来说,若判断
Figure G2009101583707D00088
服从t分布,则采样值xk未离开当前稳态过程或者采样值xk进入稳态过程;若不服从,判断采样值xk离开当前稳态过程或者采样值xk未进入稳态过程。
如上文所述,
Figure G2009101583707D00091
在一个实施例中,t分布判断可以依照以下判据公式: | x k - x ‾ k | ≤ k · t α / 2 ( k - 1 ) · s ^ k , 其中α为给定的风险系数,α值在0.05左右。若满足上述公式,则表示xk服从t分布,否则表示不服从。
1)对于采样和预处理输出的xk
2)若k=1:则x1=x1 s ^ 1 2 = 0 ;
否则,按上述公式计算xk
Figure G2009101583707D00094
3)若 | x k - x ‾ k | ≤ k · t α / 2 ( k - 1 ) · s ^ k , 则表示为稳态过程。即表示xk为稳态过程或者进入稳态过程,且稳态值X=xk
对于xk为稳态过程的,令k=k+1,并根据步骤2和步骤3继续下一个采样值的判断;对于xk为进入稳态过程的,令k=1,将xk作为该稳态过程的初始稳态数据x1,并输出时标,然后重复上述步骤2和步骤3。
否则,表示从稳态过程离开或者未进入稳态过程,令k=1,对于从稳态过程离开的,还要输出时标,重复上述步骤2和步骤3。
在一个实施例中,t分布的判拒也可以简化成 | x k - x ‾ k | ≤ A · s ^ k , 其中A为给定系数,在3~10之间。或者,更进一步可以简化成|xk-xk|≤Δ,Δ=δ·xn其中δ为给定常数,根据精度要求一般在2%~10%之间取值,优选为5%左右,xn为被测连续物理量对应的额定值。
并且在第3)步的判断中,若对于初始测量(t=0)的采样值x1,若判断为稳态过程,则记录该稳态过程的开始时间。
当然,这里也可以采用高斯分布等其他统计法对采样值xk的稳态状态进行判断,本发明不局限于上述具体实施例。
在一个实施例中,可以根据滤波器输出判断采样值是否进入或者离开稳态过程。这里,滤波器可以是任意1阶或更高阶(例如2阶、3阶等)的滤波器,例如n阶卡尔曼滤波器。
通过对采样值xk进行滤波,从而可以得到采样值xk对应的滤波分量。显然,对于不同阶的滤波器,得到的滤波分量个数也不同。这里,本发明不局限于具体个数的滤波分量,任意合适的滤波算法均可落在本发明的保护范围内。
分别判断所得的滤波分量是否超出对应设定的上限值,若均未超出,判断采样值xk未离开稳态过程;若存在一个分量超出,判断采样值xk离开稳态过程。
类似地,通过分别判断上述滤波分量是否超出对应设定的下限值,则可以判断采样值xk是否进入稳态过程。具体来说,若存在一个分量超出,判断采样值xk未进入稳态过程;若均未超出,判断采样值xk进入稳态过程。
当然,本发明可以采用其他合适方法来判断采样值xk是否进入或者离开稳态过程。
这样,可以判断新数据是否仍处于相对稳定的数值附近,还是出现了超出稳态过程的突变。
在判断采样值未离开稳态过程时,将采样值xk进行稳态值计算。对于连续物理量,在稳态过程能够以一个值X,例如平均值代表其稳态大小。
这里,可以根据均值算法或低通滤波算法计算稳态过程对应输出采样值的稳态值。稳态值X通过稳态时间段内的均值算法或低通滤波算法获得,以避免随机干扰、提高测量的精度。
在一个实施例中,稳态值X对应的均值算法公式例如表示为 x ‾ k = 1 k [ ( k - 1 ) · x ‾ k - 1 + x k ) ] , 其中xk表示稳态过程获得的采样值,k表示这些采样值的个数,xk即为根据这些采样值计算得到的稳态值X。k=1时,则x1=x1
在一个实施例中,稳态值X可以根据低通滤波器的递推算法计算获得,例如利用公式表示为xk=α·xk-1+(1-α)·xk的一阶低通滤波算法,其中α为设定的常数。
当然,本领域技术人员显然可知,本发明还可以采用更高阶(例如2阶、3阶等)滤波算法或者其他均值算法来获得稳态值X。本发明不局限于上述具体实施例。
在利用统计法的t分布进行判断的实施例中,所需的均值xk可以与稳态值X的计算方法相同。
对于上述稳态值X计算的递推算法,利用新的采样值xk对当前稳态值X进行修正,得到对应稳态过程的新稳态值X。
根据判断结果,并结合定时器输出的时标,输出每个稳态过程对应的开始和结束时间,以及每个稳态过程中对应的稳态值X。
具体来说,对于连续物理量,假定判断其以往位于稳态过程,则计算一个值X代表该稳态过程中对应已采样数据的稳态值大小。在以时间间隔采样之后,获得一个新采样数据xk;判断xk是否仍在稳态过程,是,则利用xk修正X,继续稳态过程;不是,则输出该稳态过程的结束时间以及该稳态过程的稳态值X。
另外,如果以往不在稳态过程,对于新数据xk判断其是否进入到一个稳态过程;若没有进入,则继续当前状态;否则,输出稳态过程的开始时标,并将xk作为稳态数据,建立该稳态过程对应的新的稳态值X;如此继续。
对于稳态值输出,可以根据判断结果,仅输出在稳态过程结束时对应的稳态值X,也可以根据远程需要,将稳态过程中不同时间对应的稳态值X输出。
在一个实施例中,本发明还包括可以记录步骤,从而对每个稳态过程的开始时间、结束时间以及在结束时对应的稳态值进行记录。
以上输出和记录的稳态数据比现有的间隔T时间输出一次数据获得的数据量少(多个T都在一个稳态过程)且避免了混迭错误。并且,它们真实反映了系统在该时间段上的特性,从而可作为后续系统计算的无误输入。
由于稳态时间较暂态时间要长得多,所以输出或传送的数据大幅度减少。另外,又保证了信号不失真地被记录下来。当然,由记录的数据也能够复现原信号。
在一个实施例中,还可以包括预处理步骤,用于对采样出的采样值进行预处理,然后直接或重抽样输出用于后续判断和计算。一般地,预处理包括标度变换,还可以包括坏数据剔除,而对于正弦量的测量物理量,还可以包括有效值的计算。
需要指出的是,对于轮询系统,例如在规定的时间间隔t(例如电力调度自动化系统中t=3秒)或规定的时刻(例如用电监视系统中每个整点开始的每刻钟),需要采样终端输出一次数据的情况下,本发明的测量装置和方法可以不直接输出稳态过程开始和结束的时标,即,时标不明确输出而是暗含在轮询中。当输出的数据在稳态过程时,仅输出对应所需的稳态值X,而如果处于暂态过程,则输出一个采样数值如xk,并标志该数据为非稳态数据。这种变化仍落在本发明的范围内。
本发明可以应用于热工、电力系统中的温度、压力、流量、电压、电流、功率和相角等测量和记录中,能够真实、清晰地反映系统参数的变化情况,并且显著减少记录和输出数据的数量。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同限定。

Claims (17)

1.一种连续物理量测量装置,其特征在于,所述测量装置包括:
定时单元,所述定时单元输出时基和时标;
采样单元,在所述时基控制下以满足奈奎斯特采样定理的时间间隔对连续物理量进行采样,并输出采样值;
判断单元,判断采样值是否进入或者离开稳态过程;
计算单元,在稳态过程中对所述采样单元的采样值计算稳态值;
输出单元,根据判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及所述稳态值。
2.如权利要求1所述的测量装置,其特征在于,所述判断单元根据统计法的t分布或者滤波器输出进行判断。
3.如权利要求1所述的测量装置,其特征在于,还包括预处理单元,用于对所述采样单元输出的采样值进行标度变换、坏数据去除、重抽样和/或有效值计算。
4.如权利要求1所述的测量装置,其特征在于,所述计算单元根据均值算法或低通滤波算法计算所述稳态值。
5.如权利要求1所述的测量装置,其特征在于,还包括记录单元,所述记录单元对所述稳态过程的开始时间、结束时间以及在所述结束时间对应的稳态值进行记录。
6.一种连续物理量测量方法,其特征在于,所述测量方法包括以下步骤:
在定时器输出时基的控制下以满足奈奎斯特采样定理的时间间隔对连续物理量进行采样,并输出采样值xk和时标;
判断采样值xk是否进入或者离开稳态过程;以及
根据判断结果以及所述时标,输出稳态过程的开始时间和结束时间以及在稳态过程中对采样值计算获得的稳态值X。
7.如权利要求6所述的测量方法,其特征在于,所述稳态值X根据均值算法或低通滤波算法计算获得。
8.如权利要求6所述的测量方法,其特征在于,所述判断步骤包括:
计算采样值xk对应的均值xk和方差
Figure A2009101583700003C1
判断是否服从t分布;
若服从,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;
若不服从,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。
9.如权利要求8所述的测量方法,其特征在于,所述t分布判断依据以下公式:
| x k - x ‾ k | ≤ k · t α / 2 ( k - 1 ) · s ^ k , 其中α为风险系数。
10.如权利要求6所述的测量方法,其特征在于,所述判断步骤包括:
计算采样值xk对应的均值xk和方差
Figure A2009101583700003C4
判断是否满足 | x k - x ‾ k | ≤ A · s ^ k , 其中A为给定值;
若满足,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;
若不满足,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。
11.如权利要求10所述的测量方法,其特征在于,所述A位于3~10之间。
12.如权利要求6所述的测量方法,其特征在于,所述判断步骤包括:
计算采样值xk对应的均值xk
判断是否满足|xk-xk|≤δ·xn,其中δ为给定值,xn为所述连续物理量对应的额定值;
若满足,判断采样值xk未离开稳态过程或者判断采样值xk进入稳态过程;
若不满足,判断采样值xk离开稳态过程或者判断采样值xk未进入稳态过程。
13.如权利要求12所述的测量方法,其特征在于,所述δ位于2%到10%之间。
14.如权利要求6所述的测量方法,其特征在于,所述判断步骤包括:
对采样值xk进行滤波,分别得到采样值xk的滤波分量;
分别判断所述滤波分量是否超出对应设定的上限值;以及
若均未超出,判断采样值xk未离开稳态过程;
若存在一个分量超出,判断采样值xk离开稳态过程。
15.如权利要求6所述的测量方法,其特征在于,所述判断步骤包括:
对采样值xk进行滤波,分别得到采样值xk的滤波分量;
分别判断所述滤波分量是否超出对应设定的下限值;以及
若存在一个分量超出,判断采样值xk未进入稳态过程;
若均未超出,判断采样值xk进入稳态过程。
16.如权利要求6所述的测量方法,其特征在于,还包括对采样值xk进行标度变换、坏数据剔除、重抽样和/或有效值计算的预处理步骤。
17.如权利要求6所述的测量方法,其特征在于,还包括对所述稳态过程的开始时间、结束时间以及在所述结束时间对应的稳态值进行记录的步骤。
CN2009101583707A 2009-07-08 2009-07-08 连续物理量测量装置及方法 Expired - Fee Related CN101614554B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101583707A CN101614554B (zh) 2009-07-08 2009-07-08 连续物理量测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101583707A CN101614554B (zh) 2009-07-08 2009-07-08 连续物理量测量装置及方法

Publications (2)

Publication Number Publication Date
CN101614554A true CN101614554A (zh) 2009-12-30
CN101614554B CN101614554B (zh) 2011-07-20

Family

ID=41494331

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101583707A Expired - Fee Related CN101614554B (zh) 2009-07-08 2009-07-08 连续物理量测量装置及方法

Country Status (1)

Country Link
CN (1) CN101614554B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011003347A1 (zh) * 2009-07-08 2011-01-13 保定市三川电气有限责任公司 连续物理量测量装置及方法
CN102331535A (zh) * 2011-06-09 2012-01-25 郝玉山 交流电物理量测量和数据采集装置和方法
CN102393214A (zh) * 2011-06-09 2012-03-28 郝玉山 连续物理量数据采集方法和装置
CN102914329A (zh) * 2012-10-23 2013-02-06 保定市三川电气有限责任公司 物理量数字测量或遥测的稳态处理方法及系统
CN102928014A (zh) * 2012-10-23 2013-02-13 保定市三川电气有限责任公司 电力系统数字测量或遥测处理的方法及装置
CN103217581A (zh) * 2013-03-11 2013-07-24 河北省电力公司 基于稳态遥测技术实现线路参数辨识的方法和系统
WO2014063635A1 (zh) * 2012-10-23 2014-05-01 Hao Yushan 数据采样方法与系统及其在参数辨识中的应用方法与系统
CN104457911A (zh) * 2014-12-26 2015-03-25 上海神开石油化工装备股份有限公司 一种液位检测系统及液位检测方法
CN104483836A (zh) * 2014-11-17 2015-04-01 国家电网公司 一种遥测数据的在线稳态处理方法
CN107589390A (zh) * 2017-10-25 2018-01-16 国家电网公司 一种直流电子式互感器阶跃响应测试系统显控实现方法
CN110646709A (zh) * 2019-09-29 2020-01-03 南瑞集团有限公司 一种适用于次同步振荡监控装置的数据汇集方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6122959A (en) * 1998-01-14 2000-09-26 Instrumented Sensor Technology, Inc. Method and apparatus for recording physical variables of transient acceleration events
CN1185496C (zh) * 2002-09-23 2005-01-19 清华大学 一种同步发电机空载电势相量的纯电气测量方法
CN101350518A (zh) * 2008-09-05 2009-01-21 国电南瑞科技股份有限公司 数字化保护测控装置采样值处理方法
CN101403778A (zh) * 2008-10-23 2009-04-08 泰豪科技股份有限公司 一种自适应流水线式交流参数测量方法

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011003347A1 (zh) * 2009-07-08 2011-01-13 保定市三川电气有限责任公司 连续物理量测量装置及方法
CN102331535B (zh) * 2011-06-09 2014-06-04 郝玉山 交流电物理量测量和数据采集装置和方法
CN102331535A (zh) * 2011-06-09 2012-01-25 郝玉山 交流电物理量测量和数据采集装置和方法
CN102393214A (zh) * 2011-06-09 2012-03-28 郝玉山 连续物理量数据采集方法和装置
WO2012167749A1 (zh) * 2011-06-09 2012-12-13 Hao Yushan 连续物理量数据采集方法和装置
US9143245B2 (en) 2011-06-09 2015-09-22 Yushan Hao Methods and devices for acquiring continuous physical signals
EA027159B1 (ru) * 2011-06-09 2017-06-30 Юйшань Хао Способы и устройства для получения значений непрерывных физических величин
EP2720002A1 (en) * 2011-06-09 2014-04-16 Yushan Hao Continuous physical quantity data collection method and device
EP2720002A4 (en) * 2011-06-09 2014-11-26 Yushan Hao METHOD AND DEVICE FOR CONTINUOUS DATA COLLECTION ON PHYSICAL QUANTITIES
CN102914329A (zh) * 2012-10-23 2013-02-06 保定市三川电气有限责任公司 物理量数字测量或遥测的稳态处理方法及系统
WO2014063635A1 (zh) * 2012-10-23 2014-05-01 Hao Yushan 数据采样方法与系统及其在参数辨识中的应用方法与系统
WO2014063633A1 (zh) * 2012-10-23 2014-05-01 Hao Yushan 物理量数字测量或遥测的稳态处理方法及系统
CN102928014A (zh) * 2012-10-23 2013-02-13 保定市三川电气有限责任公司 电力系统数字测量或遥测处理的方法及装置
US9438449B2 (en) 2012-10-23 2016-09-06 Yushan Hao Data sampling method and system, and application method and system thereof in parameter identification
CN102928014B (zh) * 2012-10-23 2015-05-13 保定市三川电气有限责任公司 电力系统数字测量或遥测处理的方法及装置
CN102914329B (zh) * 2012-10-23 2015-07-08 保定市三川电气有限责任公司 物理量数字测量或遥测的稳态处理方法及系统
CN103217581A (zh) * 2013-03-11 2013-07-24 河北省电力公司 基于稳态遥测技术实现线路参数辨识的方法和系统
CN103217581B (zh) * 2013-03-11 2016-03-30 河北省电力公司 基于稳态遥测技术实现线路参数辨识的方法和系统
CN104483836A (zh) * 2014-11-17 2015-04-01 国家电网公司 一种遥测数据的在线稳态处理方法
CN104483836B (zh) * 2014-11-17 2017-08-29 国家电网公司 一种遥测数据的在线稳态处理方法
CN104457911A (zh) * 2014-12-26 2015-03-25 上海神开石油化工装备股份有限公司 一种液位检测系统及液位检测方法
CN104457911B (zh) * 2014-12-26 2018-01-02 上海神开石油化工装备股份有限公司 一种液位检测系统及液位检测方法
CN107589390A (zh) * 2017-10-25 2018-01-16 国家电网公司 一种直流电子式互感器阶跃响应测试系统显控实现方法
CN110646709A (zh) * 2019-09-29 2020-01-03 南瑞集团有限公司 一种适用于次同步振荡监控装置的数据汇集方法及系统
CN110646709B (zh) * 2019-09-29 2022-08-23 南瑞集团有限公司 一种适用于次同步振荡监控装置的数据汇集方法及系统

Also Published As

Publication number Publication date
CN101614554B (zh) 2011-07-20

Similar Documents

Publication Publication Date Title
CN101614554B (zh) 连续物理量测量装置及方法
CN101614555B (zh) 连续物理量测量装置及方法
CN101876693B (zh) 一种基于电能计量芯片的终端校准系统
Koch et al. Modeling and control of aggregated heterogeneous thermostatically controlled loads for ancillary services
US4218737A (en) Revenue metering system for power companies
CN101655545B (zh) 电能表现场校验方法
EP3026773B1 (en) System and method for modelling load in an electrical power network
CN108693496A (zh) 一种基于参量退化方程的智能电能表误差预估方法
CN105656204A (zh) 一种智能变电站遥信闭环测试方法
CN101479613B (zh) 瞬时确定交流电网上的信号畸变率的方法和相关装置
CN104537233A (zh) 一种基于核密度估计的配电网伪量测生成方法
Azad et al. Grid serving deployment of smart meter data in the context of distribution grid automation
US10366403B2 (en) Distributed forecasting and pricing system
CN103884910A (zh) 一种适用于频率偏移的电力系统相量计算方法
CN102313834B (zh) 一种电力数据的转换系统、方法与装置
CN110133539A (zh) Ups的健康度预测方法、系统和计算机可读存储介质
CN112557999A (zh) 一种便于校验误差的多用户电能表和误差校验方法
RU2717378C1 (ru) Устройство и способ для режимного ограничения электропотребления объектов инфраструктурного типа
Lu et al. Exponentially weighed moving average control chart for gamma distribution with type I censoring
CN114690628A (zh) 伪雅克比矩阵参数自整定的mimo紧格式无模型控制方法
Razo-Hernandez et al. IRIG-B decoder based on FPGA for synchronization in PMUs by considering different input formats
Khurram et al. Real-time grid and der co-simulation platform for validating large-scale der control schemes
CN116908537B (zh) 一种电流电压频率计算电路和方法
CN117039893B (zh) 配电网状态确定方法、装置及电子设备
CN205120262U (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
ASS Succession or assignment of patent right

Owner name: SANCHUAN ELECTRIC CO., LTD.

Free format text: FORMER OWNER: SANCHUAN ELECTRIC CO LTD, BAODING CITY

Effective date: 20110909

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20110909

Address after: High tech Zone Hebei city Baoding province show street 071051 No. 295 high tech Innovation Service Center C block 3 layer

Patentee after: Sanchuan Power Equipment Co.,Ltd.

Address before: 071051 building C, Pioneer Center, hi tech Zone, Hebei, Baoding

Patentee before: Baoding Sanchuan Electric Co.,Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110720

Termination date: 20210708