首页 > 专利 > 杭州电子科技大学 > 基于长短时信息融合的逆变器直流回路电容故障检测方法专利详情

基于长短时信息融合的逆变器直流回路电容故障检测方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-03-22
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2019-07-23
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2021-01-01
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-03-22
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910222023.X 申请日 2019-03-22
公开/公告号 CN109946542B 公开/公告日 2021-01-01
授权日 2021-01-01 预估到期日 2039-03-22
申请年 2019年 公开/公告年 2021年
缴费截止日
分类号 G01R31/00G01R27/26 主分类号 G01R31/00
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 0
引用专利数量 0 被引证专利数量 0
非专利引证
引用专利 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 徐晓滨、俞卓辰、高海波、侯平智、胡燕祝、黄大荣、韩德强 第一发明人 徐晓滨
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 7
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江千克知识产权代理有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
周希良
摘要
本发明涉及一种基于长短时信息融合的逆变器直流回路电容故障检测方法。本发明先将采集的逆变器直流回路电压瞬时值以一定周期计算得到每一周期的电压瞬时值峰‑峰值,按一定间隔选取阈值,利用连续的Sigmoid隶属度函数将处理后的瞬时电压峰‑峰值的信息转化成短时报警证据;利用长短时信息融合公式对短时报警证据进行融合得到全局报警证据,进行报警决策,误报率和漏报率的平方和最小所对应的阈值为最优阈值;在线获取数据,利用训练得到的最优阈值得到短时报警证据,利用长短时信息融合公式进行融合得到长时报警证据,进行报警决策。本发明通过长短时信息融合得到长时报警证据,可以给出比单个短时报警证据更为准确的故障检测结果。
  • 摘要附图
    基于长短时信息融合的逆变器直流回路电容故障检测方法
  • 说明书附图:图1
    基于长短时信息融合的逆变器直流回路电容故障检测方法
  • 说明书附图:图2
    基于长短时信息融合的逆变器直流回路电容故障检测方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2021-01-01 授权
2 2019-07-23 实质审查的生效 IPC(主分类): G01R 31/00 专利申请号: 201910222023.X 申请日: 2019.03.22
3 2019-06-28 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.基于长短时信息融合的逆变器直流回路电容故障检测方法,其特征在于该方法包括以下步骤:
(1)设定逆变器直流回路电容故障报警器的辨识框架为Θ={NA,A},其中NA表示电容器容量处于正常状态,A表示电容器容量衰减到正常状态的x%,亦即,异常状态,50≤x≤
80;
(2)在逆变器运行时,以一定的采样频率采集逆变器直流回路电压瞬时值VDC,以一定的周期计算VDC的峰-峰值Vpp,得到样本序列Vpp(t),t=0,1,2,3,…,是Vpp在第t个周期的值,单位为V,定义Votp为最优阈值,记max(Vpp)和min(Vpp)分别是Vpp的最大值和最小值,那么Votp∈[min(y),max(y)];利用Sigmoid隶属度函数构建正常状态NA、异常状态A的模糊隶属度函数μA(Vpp(t))、μNA(Vpp(t)),如式(1)-(2)所示:
其中参数a∈[0.1,100];
(3)将Vpp(t)分别代入步骤(2)关于异常状态A、正常状态NA的模糊隶属度函数μNA(Vpp(t))、μA(Vpp(t))中,即可得到短时报警信度mt(NA)、mt(A),分别如下(3)-(4)式所示:
mt(NA)=μNA(Vpp(t))                      (3)
mt(A)=μA(Vpp(t))                         (4)
那么,在每个周期,能够获得由上两个信度获得向量形式的短时报警证据mt=(mt(NA),mt(A)),其中mt(A)和mt(NA)表示Vpp(t)支持“正常状态NA”和“异常状态A”这两个事件发生的程度;
(4)在步骤(3)的基础上得到了每一周期的短时报警证据之后,利用长短时信息融合公式得到融合后的全局报警证据,记为m1:t=(m1:t(A),m1:t(NA)),具体步骤如下:
(4-1)当1≤t≤2时,全局报警证据等于该周期的短时报警证据即m1:t=(mt(A),mt(NA));
(4-2)当t≥3时,求取m1:t=(m1:t(A),m1:t(NA)),具体用于计算的长短时信息融合公式如下
这里m1:t-2(A)是通过递归融合计算得到的长时报警证据,包含了0到t-2周期所有的报警信息,对于mt(A|D)的求解,若mt(NA)>mt(A),则D=NA,那么mt(A|NA)=0;若mt(A)≥mt(NA),则D=A,那么mt(A|A)=1;同理可得mt-1(A|D),mt-2(A|D),其中τt和νt表示的是融合权重,λt和γt表示条件信度权重;
对于融合权重和条件信度权重的求取,假定正常状态N A的理想证据是I1=(mI1(A),mI1(NA))=(1,0),异常状态A的理想证据是I2=(mI2(A),mI2(NA))=(0,1),那么式(5)中τt、νt、λt和γt的具体确定步骤如下:
(4-3)λt和γt的确定具体步骤如下
(4-3-1)如果长时报警证据m1:t-2(A)≥m1:t-2(NA),计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的距离为
其中参数P∈[100,200];
(4-3-2)根据式(6)和式(7)可以计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的相似度为
得到λt和γt为
(4-3-3)如果长时报警证据m1:t-2(A)其中参数P∈[100,200];
(4-3-4)根据式(6)和式(7)计算单步短时报警证据mt-1、mt和发出警报的理想证据I2之间的相似度为
得到λt和γt为
(4-4)根据步骤(4-3-4)得到的λt和γt,将mt-1、mt融合为一个证据mt-1:t为mt-1:t=λtmt-1+γtmt                   (16)
(4-5)τt和νt的确定具体步骤如下:
(4-5-1)求得mt-1:t和m1:t-2的信息熵如式(17)-(18)所示
Ht-1:t=-(mt-1:t(A)logb(mt-1:t(A))+mt-1:t(NA)logb(mt-1:t(NA)))    (17)H1:t-2=-(m1:t-2(A)logb(m1:t-2(A))+m1:t-2(NA)logb(m1:t-2(NA)))     (18)其中b=2;
(4-5-2)τt和νt可由下式(19)确定
(5)根据步骤(4)得到全局报警证据m1:t=(m1:t(A),m1:t(NA)),对于全局报警证据,如果m1:t(A)≥m1:t(NA),则发出警报,否则,不警报;
(6)在分别获取W>2000个Vpp在正常和异常下的训练样本值后,最小化这W个样本的误报率和漏报率的平方和,求得最优阈值Votp。
说明书

技术领域

[0001] 本发明涉及一种基于长短时信息融合的逆变器直流回路电容故障检测方法,属于电力电子设备故障诊断领域。

背景技术

[0002] 逆变器作为船舶电力推进系统的核心部件,直流回路电容器的故障是逆变器的常见故障之一。电容器故障一般包括短路、漏电和电容量衰减等。在实际系统运行过程中,电容器电容量衰减必然会发生,但由于电容器的具体工作环境不同,其衰减的程度和特征并没有规律,具有不确定性。电容器在现阶段主要有电解电容器、陶瓷电容器和薄膜电容器等。其中,电解电容器作为逆变器常用电容器,它的寿命一般很短,通常少于逆变器的一般服役时间,随着逆变器工作时间的变长,电容器内部电解液的损耗增多,电解液不断蒸发变稠,最终会导致电容逐渐老化,具体表现在电容量的下降。此时逆变器仍旧能够工作,但逆变器在带载运行时会经常出现过流、欠压等现象,增加设备的安全隐患。此类故障属于微小故障,电力保护系统不容易被激活,但若忽视此故障长期运行必将影响设备的使用寿命,因此电容器老化的微小故障在线诊断是必不可少的。

发明内容

[0003] 本发明的目的是提出一种基于长短时信息融合的逆变器直流回路电容故障检测方法,采集逆变器直流回路电压瞬时值,按一定周期计算电压瞬时值的峰-峰值,将其用Sigmoid隶属度阈值函数转换为相应的短时报警证据,通过长短时信息融合公式对每一短时报警证据进行融合,得到可靠的长时报警证据,利用长时报警证据进行决策。由于长时报警证据包含了短时和历史的电容电压瞬时值变化信息,利用其做出的故障报警决策比任何单个时段振动信息给出的决策结果更为可靠和准确。
[0004] 本发明包括以下步骤:
[0005] (1)设定逆变器直流回路电容故障报警器的辨识框架为Θ={NA,A},其中A表示电容器容量处于正常状态,NA表示电容器容量衰减到正常状态的x%,亦即,异常状态,50≤x≤80。
[0006] (2)在逆变器运行时,以一定的采样频率采集逆变器直流回路电压瞬时值VDC,以一定的周期计算VDC的峰-峰值Vpp,得到样本序列Vpp(t),t=0,1,2,3,…,是Vpp在第t个周期的值,单位为V,定义Votp为最优阈值,记max(Vpp)和min(Vpp)分别是Vpp的最大值和最小值,那么Votp∈[min(y),max(y)];利用Sigmoid隶属度函数构建正常状态A、异常状态NA的模糊隶属度函数μA(Vpp(t))、μNA(Vpp(t)),如式(1)-(2)所示:
[0007]
[0008]
[0009] 其中参数a∈[0.1,100]。
[0010] (3)将Vpp(t)分别带入步骤(2)关于正常状态NA、异常状态A的模糊隶属度函数μNA(Vpp(t))、μA(Vpp(t))中,即可得到短时报警信度mt(NA)、mt(A),分别如下(3)-(4)式所示:
[0011] mt(NA)=μNA(Vpp(t))                      (3)
[0012] mt(A)=μA(Vpp(t))                         (4)
[0013] 那么,在每个周期,能够获得由上两个信度获得向量形式的短时报警证据mt=(mt(NA),mt(A)),其中mt(A)和mt(NA)表示Vpp(t)支持“报警A”和“正常NA”这两个事件发生的程度。
[0014] (4)在步骤(3)的基础上得到了每一周期的短时报警证据之后,利用长短时信息融合公式得到融合后的全局报警证据,记为m1:t=(m1:t(A),m1:t(NA)),具体步骤如下:
[0015] (4-1)当1≤t≤2时,因在之前没有相关的证据信息,所以t=1或2时的全局报警证据等于该周期的短时报警证据即m1:t=(mt(A),mt(NA));
[0016] (4-2)当t≥3时,求取m1:t=(m1:t(A),m1:t(NA)),具体用于计算的长短时信息融合公式如下
[0017]
[0018] 这里m1:t-2(A)是通过递归融合计算得到的长时报警证据,包含了0到t-2周期所有的报警信息,对于mt(A|D)的求解,若mt(NA)>mt(A),则D=NA,那么mt(A|NA)=0;若mt(A)≥mt(NA),则D=A,那么mt(A|A)=1;同理可得mt-1(A|D),mt-2(A|D),其中τt和νt表示的是融合权重,λt和γt表示条件信度权重;
[0019] 对于融合权重和条件信度权重的求取,假定发出警报(A)的理想证据是I1=(mI1(A),mI1(NA))=(1,0),未报警(NA)的理想证据是I2=(mI2(A),mI2(NA))=(0,1),那么式(5)中τt、νt、λt和γt的具体确定步骤如下:
[0020] (4-3)λt和γt的确定具体步骤如下
[0021] (4-3-1)如果长时报警证据m1:t-2(A)≥m1:t-2(NA),计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的距离为
[0022]
[0023]
[0024] 其中参数P∈[100,200];
[0025] (4-3-2)根据式(6)和式(7)可以计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的相似度为
[0026]
[0027]
[0028] 可以得到λt和γt为
[0029]
[0030] (4-3-3)如果长时报警证据m1:t-2(A)
[0031]
[0032]
[0033] 其中参数P∈[100,200];
[0034] (4-3-4)根据式(6)和式(7)可以计算单步短时报警证据mt-1、mt和发出警报的理想证据I2之间的相似度为
[0035]
[0036]
[0037] 可以得到λt和γt为
[0038]
[0039] (4-4)根据步骤(4-3-4)得到的λt和γt,将mt-1、mt融合为一个证据mt-1:t为[0040] mt-1:t=λtmt-1+γtmt                   (16)
[0041] (4-5)τt和νt的确定具体步骤如下:
[0042] (4-5-1)求得mt-1:t和m1:t-2的信息熵如式(17)-(18)所示
[0043] Ht-1:t=-(mt-1:t(A)logb(mt-1:t(A))+mt-1:t(NA)logb(mt-1:t(NA)))    (17)[0044] H1:t-2=-(m1:t-2(A)logb(m1:t-2(A))+m1:t-2(NA)logb(m1:t-2(NA)))     (18)[0045] 其中b=2;
[0046] (4-5-2)τt和νt可由下式(19)确定
[0047]
[0048] (5)根据步骤(4)得到全局报警证据m1:t=(m1:t(A),m1:t(NA)),对于全局报警证据,如果m1:t(A)≥m1:t(NA),则发出警报,否则,不警报。
[0049] (6)在分别获取W>2000个Vpp在正常和异常下的训练样本值后,可以最小化这W个样本的误报率和漏报率的平方和,求得最优阈值Votp。
[0050] 本发明提出的基于长短时信息融合的逆变器直流回路电容故障检测方法。该方法先将采集的逆变器直流回路电压瞬时值以一定周期计算得到每一周期的电压瞬时值峰-峰值,按一定间隔选取阈值,利用连续的Sigmoid隶属度函数将处理后的瞬时电压峰-峰值的信息转化成短时报警证据;利用长短时信息融合公式对短时报警证据进行融合得到长时报警证据,进行报警决策,误报率和漏报率的平方和最小所对应的阈值为最优阈值;在线获取数据,利用训练得到的最优阈值得到短时报警证据,利用长短时信息融合公式进行融合得到长时报警证据,进行报警决策。根据本发明方法编制的程序(编译环境Matlab)可以在计算机上运行,并电压测量仪、数据采集器等硬件组成逆变器直流回路电容故障检测系统,从而实现对逆变器直流回路电容故障检测并报警。

实施方案

[0053] 本发明提出的一种基于长短时信息融合的逆变器直流回路电容故障检测方法,其流程框图如图1所示,包括以下各步骤:
[0054] (1)设定逆变器直流回路电容故障报警器的辨识框架为Θ={NA,A},其中A表示电容器容量处于正常状态,NA表示电容器容量衰减到正常状态的x%,亦即,异常状态,50≤x≤80。
[0055] (2)在逆变器运行时,以一定的采样频率采集逆变器直流回路电压瞬时值VDC,以一定的周期计算VDC的峰-峰值Vpp,得到样本序列Vpp(t),t=0,1,2,3,…,是Vpp在第t个周期的值,单位为V,定义Votp为最优阈值,记max(Vpp)和min(Vpp)分别是Vpp的最大值和最小值,那么Votp∈[min(y),max(y)];利用Sigmoid隶属度函数构建正常状态A、异常状态NA的模糊隶属度函数μA(Vpp(t))、μNA(Vpp(t)),如式(1)-(2)所示:
[0056]
[0057]
[0058] 其中参数a∈[0.1,100]。
[0059] (3)将Vpp(t)分别带入步骤(2)关于正常状态NA、异常状态A的模糊隶属度函数μNA(Vpp(t))、μA(Vpp(t))中,即可得到短时报警信度mt(NA)、mt(A),分别如下(3)-(4)式所示:
[0060] mt(NA)=μNA(Vpp(t))                      (3)
[0061] mt(A)=μA(Vpp(t))                         (4)
[0062] 那么,在每个周期,能够获得由上两个信度获得向量形式的短时报警证据mt=(mt(NA),mt(A)),其中mt(A)和mt(NA)表示Vpp(t)支持“报警A”和“正常NA”这两个事件发生的程度。
[0063] (4)在步骤(3)的基础上得到了每一周期的短时报警证据之后,利用长短时信息融合公式得到融合后的全局报警证据,记为m1:t=(m1:t(A),m1:t(NA)),具体步骤如下:
[0064] (4-1)当1≤t≤2时,因在之前没有相关的证据信息,所以t=1或2时的全局报警证据等于该周期的短时报警证据即m1:t=(mt(A),mt(NA))。
[0065] (4-2)当t ≥3时,求取m1:t=(m1:t(A),m1:t(NA)),具体用于计算的长短时信息融合公式如下
[0066]
[0067] 这里m1:t-2(A)是通过递归融合计算得到的长时报警证据,包含了0到t-2周期所有的报警信息,对于mt(A|D)的求解,若mt(NA)>mt(A),则D=NA,那么mt(A|NA)=0;若mt(A)≥mt(NA),则D=A,那么mt(A|A)=1;同理可得mt-1(A|D),mt-2(A|D),其中τt和νt表示的是融合权重,λt和γt表示条件信度权重。
[0068] 对于融合权重和条件信度权重的求取,假定发出警报(A)的理想证据是I1=(mI1(A),mI1(NA))=(1,0),未报警(NA)的理想证据是I2=(mI2(A),mI2(NA))=(0,1),那么式(5)中τt、νt、λt和γt的具体确定步骤如下:
[0069] (4-3)λt和γt的确定具体步骤如下
[0070] (4-3-1)如果长时报警证据m1:t-2(A)≥m1:t-2(NA),计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的距离为
[0071]
[0072]
[0073] 其中参数P∈[100,200];
[0074] (4-3-2)根据式(6)和式(7)可以计算单步短时报警证据mt-1、mt和发出警报的理想证据I1之间的相似度为
[0075]
[0076]
[0077] 可以得到λt和γt为
[0078]
[0079] (4-3-3)如果长时报警证据m1:t-2(A)
[0080]
[0081]
[0082] 其中参数P∈[100,200];
[0083] (4-3-4)根据式(6)和式(7)可以计算单步短时报警证据mt-1、mt和发出警报的理想证据I2之间的相似度为
[0084]
[0085]
[0086] 可以得到λt和γt为
[0087]
[0088] (4-4)根据步骤(4-3-4)得到的λt和γt,将mt-1、mt融合为一个证据mt-1:t为[0089] mt-1:t=λtmt-1+γtmt                     (16)
[0090] (4-5)τt和νt的确定具体步骤如下:
[0091] (4-5-1)求得mt-1:t和m1:t-2的信息熵如式(17)-(18)所示
[0092] Ht-1:t=-(mt-1:t(A)logb(mt-1:t(A))+mt-1:t(NA)logb(mt-1:t(NA)))    (17)[0093] H1:t-2=-(m1:t-2(A)logb(m1:t-2(A))+m1:t-2(NA)logb(m1:t-2(NA)))      (18)[0094] 其中b=2;
[0095] (4-5-2)τt和νt可由下式(19)确定
[0096]
[0097] 为了加深对步骤(4)的理解,这里举个例子解释长短时信息融合的整个过程,假设a=6,前4个周期对应的短时报警证据如下表1所示
[0098] 表1短时报警证据
[0099]t mt=(mt(A),mt(NA))
1 (0.9,0.1)
2 (0.7,0.3)
3 (0.45,0.55)
4 (0.8,0.2)
[0100] 根据步骤(4)计算长时报警证据
[0101] 当t=1,2时,根据步骤(4-1)可得m1:1=m1=(0.9,0.1),m1:2=m2=(0.7,0.3)。
[0102] 当t=3时,根据步骤(4-3) ,取P=100,m1:1(A)≥m1:1(NA),则Sim(m3,
I1)=0.6436,Sim(m2,I1)=0.7680,λ3=0.5441,γ3=0.4559,再根据步骤(4-4)得到m2:3=(0.586025,0.413975),再根据步骤(4-5),b=2,得到H2:3=0.9785,H1:1=0.4690,ν3=
0.3240,τ3=0.6760,又因为m2(A)>m2(NA),m2(A|A)=1,m2(NA|A)=0,m3(NA)>m3(A),m3(A|NA)=0,m3(NA|NA)=1,因此根据式(5)得到m1:3=(0.7847,0.2154)。
[0103] 当t=4时,根据步骤(4-3) ,取P=100,m1:2(A)≥m1:2(NA),则Sim(m3,
I1)=0.6436,Sim(m4,I1)=0.8324,λ4=0.4360,γ4=0.5640,再根据步骤(4-4)得到m3:4=(0.6474,0.3526),再根据步骤(4-5),b=2,得到H3:4=0.9364,H1:2=0.8813,ν3=0.4848,τ3=0.5152,又因为m3(NA)>m3(A),m3(A|NA)=0,m3(NA|NA)=1,m4(A)>m4(NA),m4(A|A)=1,m3(NA|A)=0,因此根据式(5)得到m1:4=(0.6341,0.3659)。
[0104] (5)根据步骤(4)得到全局报警证据m1:t=(m1:t(A),m1:t(NA)),对于全局报警证据,如果m1:t(A)≥m1:t(NA),则发出警报,否则,不警报。
[0105] 在步骤(4)中的例子中得到的4个长时报警证据,对其进行决策,给出报警结果如下表2。
[0106] 表2报警结果输出
[0107]
[0108] 通过表中的数据可以看出,在正常状态下,利用短时报警证据进行决策会有一定的误差,而基于本发明方法融合后的长时报警证据更加真实的反映了现实的情况。
[0109] (6)在分别获取W>2000个Vpp在正常和异常下的训练样本值后,可以最小化这W个样本的误报率和漏报率的平方和,求得最优阈值Votp。
[0110] 以下结合附图,详细的介绍本发明方法的实施例:
[0111] 本发明方法的流程图如图1所示,核心部分是:从船舶电力推进系统采集逆变器直流回路电压瞬时值;以一定的周期计算电压瞬时值的峰-峰值,按一定间隔选取阈值,利用Sigmoid隶属度函数将峰-峰值转化为短时报警证据;利用长短时信息融合公式融合短时报警证据得到长时报警证据,进行报警决策,误报率和漏报率的平方和最小所对应的阈值为最优阈值。
[0112] 以下结合电力推进船舶“江苏路渡3011”船,详细介绍本发明方法的各个步骤,通过实验数据在最优阈值的条件下,得到的误报率与漏报率要比在传统的方法上更加低。
[0113] 1、实验数据的采集
[0114] 基于上述船例,使用MATLAB/Simulink搭建电力推进仿真模型,以50kHz的采样频率采集直流回路电容容量正常状态和容量为80%的电压瞬时值。
[0115] 2、求得最优阈值
[0116] 以0.002秒为周期处理电压瞬时值得到4000个Vpp,正常2000个,异常2000个,Vpp(1)到Vpp(2000)这2000个数据处于正常的工作状态,Vpp(2001)到Vpp(4000)这2000个数据处于电容容量为80%的工作状态,如图2所示。min(Vpp)=-10.7318,max(Vpp)=-2.9318,以0.1为间隔选取阈值,设定a=6根据本发明步骤(3)-(5)分别得到在每一个阈值下的全局报警证据m1:t=(m1:t(A),m1:t(NA)),t=1,2,3,…,4000,并作报警决策,误报率和漏报率的平方和最小对应的最优阈值为-7.2318,误报率为8.15%,漏报率为13.2%。使用相同数据相同遍历方法求得三阶滑动平均滤波法下误报率和漏报率的平方和最小对应的最优阈值为-
7.0193,误报率为9.05%,漏报率为24.15%。
[0117] 3、测试实验
[0118] 按上述方法重新在线采集10组数据,进行处理后得到10组Vpp,每组正常2000个,异常2000个,记为 t=1,2,3,…,4000,其中 到 处于正常的工作状态, 到处于电容容量为80%的工作状态 本发明方法下最优阈值
为-7.2318,对于每一组数据,根据本发明步骤(2)和步骤(3)计算短时报警证据,再根据本发明步骤(4)进行长短时信息融合得到全局报警证据,最后根据本发明步骤(5),进行报警决策,得到误报率、漏报率。在本发明方法下,10组数据平均误报率为12.69%,漏报率为
7.51%。三阶滑动平均滤波法的最优阈值为-7.0193,使用三阶滑动平均滤波法计算每一组的误报率、漏报率,计算得到10组平均误报率为23.64%,漏报率为9.85%。从最后的报警结果可以看出,本发明在最优阈值的下产生的误报率和漏报率都比传统的三阶滑动平均滤波法更加精确。

附图说明

[0051] 图1是本发明方法的程序流程框图;
[0052] 图2本发明方法的实施例中用来训练的Vpp(t)序列。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号