CN105654645B - 一种光纤安防信号处理方法及*** - Google Patents

一种光纤安防信号处理方法及*** Download PDF

Info

Publication number
CN105654645B
CN105654645B CN201610046738.0A CN201610046738A CN105654645B CN 105654645 B CN105654645 B CN 105654645B CN 201610046738 A CN201610046738 A CN 201610046738A CN 105654645 B CN105654645 B CN 105654645B
Authority
CN
China
Prior art keywords
mrow
signal
msub
msubsup
time
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
CN201610046738.0A
Other languages
English (en)
Other versions
CN105654645A (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.)
BEIJING BUPT-GUOAN TECHNOLOGY Corp
Original Assignee
BEIJING BUPT-GUOAN TECHNOLOGY Corp
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 BEIJING BUPT-GUOAN TECHNOLOGY Corp filed Critical BEIJING BUPT-GUOAN TECHNOLOGY Corp
Priority to CN201610046738.0A priority Critical patent/CN105654645B/zh
Publication of CN105654645A publication Critical patent/CN105654645A/zh
Application granted granted Critical
Publication of CN105654645B publication Critical patent/CN105654645B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00—Burglar, theft or intruder alarms
    • G08B13/18—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength
    • G08B13/181—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems
    • G08B13/183—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems by interruption of a radiation beam or barrier
    • G08B13/186—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems by interruption of a radiation beam or barrier using light guides, e.g. optical fibres
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B29/00—Checking or monitoring of signalling or alarm systems; Prevention or correction of operating errors, e.g. preventing unauthorised operation
    • G08B29/18—Prevention or correction of operating errors
    • G08B29/185—Signal analysis techniques for reducing or preventing false alarms or for enhancing the reliability of the system

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Security & Cryptography (AREA)
  • Alarm Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种光纤安防信号处理方法及***，包括以下步骤：S100：获取光纤布控区域内光信号的相位差异；S200：将相位差异转变为幅值差异，获得时域分析信号包；S300：得到事件信号；S400：提取所述事件信号的交流部分，获取事件信号的时域信号；S500：对所述时域信号进行快速傅立叶变换获取所述事件信号的频域信号；S600：根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构为一个数据包；S700：对所述数据包进行聚类分析，构建特征模板；S800：根据所述特征模板判断是否存在人为入侵。本发明降低了误报率。

Description

一种光纤安防信号处理方法及***
技术领域
本发明属于安防领域，特别涉及一种光纤安防信号处理方法及***。
背景技术
周界安防***在国防和民用领域都具有非常重要的应用，它主要用于边境线、军事基地、仓库、营房、政府设施、机场、核电站以及监狱等重要区域的周界入侵监测。目前的周界安防技术主要有泄漏电缆、微波对射、红外对射以及光纤传感技术等。光纤周界安防***是一种对威胁区域安全的突发事件进行监控和警报的现代防御体系，是基于分布式光纤传感技术应用于周界监控防护的新***。由于光纤及光纤传感器具有体积小、重量轻、抗干扰能力强、灵敏度高、工作可靠性高、成本低廉以及无需外场供电等优点同时还可作为信号传输通道的特点让其在其他周界安防技术中脱颖而出。在实际工程应用中，传感光纤大多数是裸露在外界环境中，这种设计独特的传感光纤对运动、压力和振动非常敏感。它可沿这围栏、围墙铺设来探测攀爬和敲击，也可铺设在土壤草坪下探测踩踏。但光纤的高灵敏度必然带来***的大量的报警，而基于时域信号能量的分析体系不足以区分大量事件预警，从而导致较高的误报率。已有技术中一般是取各种时频特性构造。环境产生较大变化时，安防***性能大打折扣。市场上一般光纤周界安防***数据包完全没有规律性，数据包并没有去验证。
发明内容
针对现有技术的缺陷，本发明提供了一种光纤安防信号处理方法及***。
一种光纤安防信号处理方法，包括以下步骤：S100：获取光纤布控区域内光信号的相位差异；S200：将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包；S300：以预设的时间跨度划分所述时域分析信号包，得到事件信号；S400：提取所述事件信号的交流部分，获取事件信号的时域信号；S500：对所述时域信号进行快速傅立叶变换获得事件信号的频域信号；S600：根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包；S700：对所述数据包进行聚类分析，构建特征模板；S800：根据所述特征模板判断是否存在人为入侵。
一种光纤安防信号处理***，包括以下模块：相位差异获取模块，用于获取光纤布控区域内光信号的相位差异；时域分析信号包获取模块，用于将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包；划分模块，用于以预设的时间跨度划分所述时域分析信号包，得到事件信号；第一获取模块，提取所述事件信号的交流部分，获取事件信号的时域信号；第二获取模块，用于对所述时域信号进行快速傅立叶变换获取所述事件信号的频域信号；重构模块，用于根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包；特征模板构建模块，用于对所述数据包进行聚类分析，构建特征模板；判断模块，用于根据所述特征模板判断是否存在人为入侵
本发明的有益效果是：本发明在充分分析人类行为特征的基础上，提出一种基于时域去噪，频域滤波的信号识别新方法。在保证不漏报的情况下，收集绝大多数事件信号，通过时域去噪、频域滤波压缩提取具有代表性的频域特征，然后通过与聚类分析构建的各类事件模板进行相似性分析。本发明在保证识别时间和报警率的情况下，降低误报率，为光纤周界安防***提供重要的支撑。
附图说明
图1是本发明光纤安防信号处理方法的流程图；
图2是步骤S500的流程图；
图3是步骤S700的流程图；
图4是本发明光纤安防信号处理***的结构示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明，使本发明的上述及其它目的、特征和优势将更加清晰。在全部附图中相同的附图标记指示相同的部分。并未刻意按比例绘制附图，重点在于示出本发明的主旨。
实施例1
首先对本发明的光纤安防信号处理方法做介绍，请参阅图1，本发明的安防信号处理方法在充分分析人类行为特征的基础上，提出了一种基于时域去噪，频域滤波的信号识别方法，收集绝大多数事件信号，通过时域去噪、频域滤波压缩提取具有代表性的频域特征，然后通过与聚类分析构建的各类事件模板进行相似性分析，从而判断是否发生入侵事件。
S100：获取光纤布控区域内光信号的相位差异。可以通过M-Z型干涉仪获取所述光信号的相位差异，使用M-Z型光纤周界安防***实时反馈布控区域内的光信号。如果光信号相位出现异常波动，则认为有事件发生。
S200：将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包。当有事件发生时，光传感单元(M-Z型)按时间上的采样率将光信号相位的差异转变成电信号上幅值的差异，获得时域分析信号包。
S300：以预设的时间跨度划分所述时域分析信号包，得到事件信号。时域预处理：当光纤相位产生瞬时异动的时候，通过分析通常持续时间不到1秒。相对应的，人的瞬时行为和自然现象的瞬时效应对光纤相位的影响持续是相近的。但是，从人体力学的角度分析，人的瞬时行为是在一定适用范围内的。如人平均一分钟正常行走60～80步。按100米的距离进行奔跑分析，人的步伐为60步左右，奔跑时间为12秒左右。1秒钟内人的出拳次数在5次左右。从最短受力时间来说，人的行为是有趋势的，受限于身体各部分肌肉与骨骼的力学特性，这一切都无法超出肢体动作的速度和频率的能力限度。所以，人对光纤作用的影响从力量速度而言，取决于动作的方向和动作轨迹；从频率而言，取决于动作的部位和方式。人的行为通常是一连串动作，从整体上分析，自然环境的干扰也可能是不连续，间断的一连串干扰。动物的干扰也因符合肌肉与骨骼的力学特性属于一段范围内的速度与频率。从力学角度上分析，要持续保持受力，就必须持续有力做功。而持续做功的典型为机器振动，异常的自然环境。所以从时域上分析，用于实现光纤周界安防的时间分析上限可以为10s。同时从工程实践的领域分析，一个预警***需要在最短时间内区分触发的事件。充分考虑工程中人行为的特性，大致将单个时域分析信号包划为0.1～1s。通过实验测试不同时间跨度的事件数据包，本发明发现入侵行为中数据包时间跨度为0.25秒左右的时候实际区分度最高。所以本发明建立以0.25秒做为事件数据分析时间跨度的时域体系。
S400：提取所述事件信号的交流部分，获取事件信号的时域信号。采集到的事件信号中会存在一些直流信号，在事件分析中只需要提取信号交流部分特征，因此需要通过算法将信号中的直流信号滤除。每帧信号定义为Si(n)，定义信号均值为：
令直流信号分量部分为则交流部分为同时，入侵信号大致在100KHz以内。通过对事件信号降采样。降采样得到压缩后的数据包。
S500：对所述时域信号进行快速傅立叶变换获取所述事件信号的频域信号。由于人为入侵信号与环境噪声的时域波形的能量和过零率上有一些差别，因此本发明根据信号的短时过零率和短时能量大小来提取事件信号的时域特征。如图2所示，其包括以下步骤：
S501：计算短时能量和短时过零率。定义接收到的时域信号为Si(n)，每帧信号短时能量为
设每帧的短时过零率为
计算环境噪声的短时能量和短时过零率，假设前10帧信号为环境噪声，首先求出每帧噪声的均方差，将这10个均方差的均值作为信号短时过零率的直流偏置。求出前10帧的短时能量和短时过零率的均值Zmean；、Emean和标准差Zstd、Estd，就可以得到其初始值Z0＝Zmean+2*Zstd和E0＝Emean+2*Estd,设置两个系数Ecoef和Zcoef作为阈值，改变两个系数值用于调节***的灵敏度。
S502：每过预设时间，重复步骤S501，只计算小于阈值的帧，修改阈值；预设的时间可以是1小时、1天等。
S503：提取出事件信号的时域信号x(n)。若此时有一帧信号Si(n)被判定为入侵信号，取出该帧信号的前一帧Si-1(n)平分为5个子帧，从后往前分别计算这几个子帧的短时能量。取出子帧中短时能量大于阈值的几个子帧，作为该次入侵信号的起点。同样的方法找出入侵信号的终点。提取出信号x(n)。
S504：频域特征处理：在事件信号在时域上提取特征后，事件信号再通过FFT变换，再将频域数据归一化，作为频域特征。对时域信号x(n)做傅里叶变换得出频域特性：
然后对频域信号进行归一化：
S600：根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包。时域特征取时域信号的最大值、最小值和均值作为时域特征；频域特征通过聚类分析对比模板可以选取差异率较大的1KHz以内的频域信号作为频域特征。
S700：对所述数据包进行聚类分析，构建特征模板。在***正式使用前，聚合大量事件的特征进行聚类分析，通过聚类分析构建出各种类别的特征模板。再将各个类别的模板分成2类：一类为人为；一类为非人为。本发明将每一个入侵信号的时频特征数据看作一组向量。应用K均值聚类方法找出每类入侵信号一个或几个模版向量，请参阅图3。
S701：根据实际情况确定分类数k，从n个数据对象选择2个距离最大的对象xi1，xi2作为聚点；
S702：选择第3个聚点xi3，满足如下关系：
min{d(xi3，xir)}＝max{min{d(xj，xir)}}
其中：r＝1，2；j≠i1，i2；
S703：重复步骤S702直到选择出第k个聚点Xik，得到k个初始聚点的集合：
分类原则记为：
于是，将样品分成不相交的k类，得到一个初始分类：
S704：根据初始分类重新计算聚点，规则如下：
S705：重复步骤S704m次后，当前分类G(m)与前次分类G(m-1)相等时，即G(m)＝G(m-1)时，则计算结束；
S706：将作为一组模版保存下来。
S800：根据所述特征模板判断是否存在人为入侵。将实时数据与特征模版对比相似度，取相似度最大的类别模板。然后再判断模板是属于人为还是非人为。以此作为输出，人为生成报警信号，非人为生成提示信号。由于余弦相似度越大，则表明入侵信号于模版的相似度越高。计算入侵信号x与模版之间的余弦相似度:
得到一组余弦相似度数组xcor。因此求出xcor中最大值xcori和类别i，并将该入侵信号作为第i类入侵信号。
实施例2
相应的，如图4所示，本发明还提供了一种光纤安防信号处理***，包括以下模块：相位差异获取模块，用于获取光纤布控区域内光信号的相位差异；时域分析信号包获取模块，用于将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包；划分模块，用于以预设的时间跨度划分所述时域分析信号包，得到事件信号；第一获取模块，提取所述事件信号的交流部分，获取事件信号的时域信号；第二获取模块，用于对所述时域信号进行快速傅立叶变换获取所述事件信号的频域信号；重构模块，用于根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包；特征模板构建模块，用于对所述数据包进行聚类分析，构建特征模板；判断模块，用于根据所述特征模板判断是否存在人为入侵。
本发明根据人体工程学理论上制订了一个合理的数据包时间跨度划分算法，并反复实验验证实践过程中合适的数据包时间跨度。构造数据特征模板，不断收集各种事件模板；在实际项目中保护周界安全，根据实际环境调节使用的模板。如果产生异常环境时，***已存不常见的非人为或者人为特征模板可以保证并克服恶劣环境。
在以上的描述中阐述了很多具体细节以便于充分理解本发明。但是以上描述仅是本发明的较佳实施例而已，本发明能够以很多不同于在此描述的其它方式来实施，因此本发明不受上面公开的具体实施的限制。同时任何熟悉本领域技术人员在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。

Claims (10)

1.一种光纤安防信号处理方法，其特征在于，包括以下步骤：
S100：获取光纤布控区域内光信号的相位差异；
S200：将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包；
S300：以预设的时间跨度划分所述时域分析信号包，得到事件信号；
S400：提取所述事件信号的交流部分，获取事件信号的时域信号；
S500：对所述时域信号进行快速傅立叶变换获得事件信号的频域信号；
S600：根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包；
S700：对所述数据包进行聚类分析，构建特征模板；
S800：根据所述特征模板判断是否存在人为入侵。
2.根据权利要求1所述的光纤安防信号处理方法，其特征在于，通过M-Z型干涉仪获取所述光信号的相位差异。
3.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述预设的时间跨度为0.1-1S。
4.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述预设的时间跨度为0.25S。
5.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述步骤S400具体包括:每帧事件信号定义为Si(n)，定义事件信号均值为
<mrow> <mover> <msub> <mi>S</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mo>&amp;lsqb;</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>s</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>/</mo> <mi>n</mi> </mrow>
式中n代表第n各采样点，令直流信号分量部分为
<mrow> <msub> <mi>S</mi> <mrow> <mi>d</mi> <mi>c</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mover> <msub> <mi>S</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> </mrow>
则交流部分为
<mrow> <msub> <mi>S</mi> <mrow> <mi>a</mi> <mi>c</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mover> <msub> <mi>S</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>.</mo> </mrow>
6.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述步骤S500具体包括：
S501：接收到的每帧事件信号为Si(n)，每帧事件信号的短时能量Ei为
<mrow> <msub> <mi>E</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow>
设每帧的短时过零率Zi为
<mrow> <msub> <mi>Z</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <mi>sgn</mi> <mo>&amp;lsqb;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mi>sgn</mi> <mo>&amp;lsqb;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>|</mo> <mo>,</mo> </mrow> 式中sgn为符号函数；
计算环境噪声的短时能量和短时过零率，假设前10帧事件信号为环境噪声，首先求出每帧噪声的均方差，将这10个均方差的均值作为信号短时过零率的直流偏置，求出前10帧的短时能量和短时过零率的均值Zmean、Emean和标准差Zstd、Estd，得到其初始值Z0＝Zmean+2*Zstd和E0＝Emean+2*Estd，设置两个系数Ecoef和Zcoef作为阈值，两个阈值用于调节***灵敏度；
S502：每过预设时间，重复步骤S501，只计算小于阈值的帧，修改阈值；
S503：若此时有一帧信号Si(n)被判定为入侵信号，取出该帧信号的前一帧Si-1(n)平分为5个子帧，从后往前分别计算这几个子帧的短时能量，取出子帧中短时能量大于阈值的几个子帧，作为该次入侵信号的起点，同理找出入侵信号的终点，提取出事件信号的时域信号x(n)；
S504：对时域信号x(n)做傅里叶变换得出频域信号：
<mrow> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msubsup> <mi>W</mi> <mi>N</mi> <mrow> <mi>k</mi> <mi>n</mi> </mrow> </msubsup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>;</mo> <msub> <mi>W</mi> <mi>N</mi> </msub> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> <mi>N</mi> </mfrac> </mrow> </msup> </mrow>
然后对频域信号进行归一化：
<mrow> <msub> <mi>X</mi> <mrow> <mi>n</mi> <mi>o</mi> <mi>r</mi> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>X</mi> <msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mi>min</mi> </msub> </mrow> <mrow> <mi>X</mi> <msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <mi>X</mi> <msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> <mo>.</mo> </mrow>
7.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述步骤S600具体包括：将所述时域信号的最大值、最小值和均值作为所述时域特征；通过聚类分析对比模板选取差异率较大的1KHz以内的频域信号作为频域特征。
8.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述步骤S700具体包括：
S701：确定分类数k，从n个数据对象选择2个距离最大的对象xi1，xi2作为聚点；
S702：选择第3个聚点xi3，满足如下关系：
min{d(xi3，xir}＝max{min{d(xj，xir)}}
其中：r＝1，2；j≠i1，i2；
S703：重复步骤S702直到选择出第k个聚点Xik，得到k个初始聚点的集合：
<mrow> <msup> <mi>L</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>x</mi> <mn>1</mn> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>x</mi> <mn>2</mn> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>}</mo> </mrow>
分类原则记为：
<mrow> <msubsup> <mi>G</mi> <mi>i</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mo>{</mo> <mi>x</mi> <mo>:</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msubsup> <mi>x</mi> <mn>1</mn> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msubsup> <mi>x</mi> <mi>j</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>
于是，将样品分成不相交的k类，得到一个初始分类：
<mrow> <msup> <mi>G</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>G</mi> <mn>1</mn> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>G</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>o</mi> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>G</mi> <mi>k</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> <mo>}</mo> </mrow>
S704：根据初始分类重新计算聚点，规则如下：
<mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>n</mi> <mi>i</mi> </msub> </mfrac> <munder> <mo>&amp;Sigma;</mo> <mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>&amp;Element;</mo> <msubsup> <mi>G</mi> <mi>i</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </msubsup> </mrow> </munder> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>k</mi> </mrow>
S705：重复步骤S704m次后，得到G(m)＝G(m-1)，则计算结束；
S706：将作为一组模版保存下来。
9.根据权利要求1所述的光纤安防信号处理方法，其特征在于，所述步骤S800包括：
计算入侵信号x与模版之间的余弦相似度:
<mrow> <msub> <mi>xcor</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>x</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> </msubsup> </mrow> <mrow> <msqrt> <msup> <mi>x</mi> <mn>2</mn> </msup> </msqrt> <mo>&amp;CenterDot;</mo> <msqrt> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msubsup> </msqrt> </mrow> </mfrac> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>k</mi> </mrow>
得到一组余弦相似度数组xcor,求出xcor中最大值xcori和类别i，并将该入侵信号作为第i类入侵信号。
10.一种光纤安防信号处理***，其特征在于，包括以下模块：
相位差异获取模块，用于获取光纤布控区域内光信号的相位差异；
时域分析信号包获取模块，用于将光信号的相位差异按时间上的采样率转变为电信号上的幅值差异，获得时域分析信号包；
划分模块，用于以预设的时间跨度划分所述时域分析信号包，得到事件信号；
第一获取模块，提取所述事件信号的交流部分，获取事件信号的时域信号；
第二获取模块，用于对所述时域信号进行快速傅立叶变换获取所述事件信号的频域信号；
重构模块，用于根据所述时域信号和频域信号获取所述事件信号的时域特征和频域特征，重构所述频域特征和时域特征为一个数据包；
特征模板构建模块，用于对所述数据包进行聚类分析，构建特征模板；
判断模块，用于根据所述特征模板判断是否存在人为入侵。
CN201610046738.0A 2016-01-22 2016-01-22 一种光纤安防信号处理方法及*** Active CN105654645B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610046738.0A CN105654645B (zh) 2016-01-22 2016-01-22 一种光纤安防信号处理方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610046738.0A CN105654645B (zh) 2016-01-22 2016-01-22 一种光纤安防信号处理方法及***

Publications (2)

Publication Number Publication Date
CN105654645A CN105654645A (zh) 2016-06-08
CN105654645B true CN105654645B (zh) 2018-04-10

Family

ID=56484298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610046738.0A Active CN105654645B (zh) 2016-01-22 2016-01-22 一种光纤安防信号处理方法及***

Country Status (1)

Country Link
CN (1) CN105654645B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106067244B (zh) * 2016-06-29 2018-08-17 中国科学院上海微***与信息技术研究所 一种用于周界探测***的传感器阈值自适应调节方法及***
CN108257364B (zh) * 2017-12-29 2020-08-14 北京航天控制仪器研究所 一种提高分布式光纤监测***报警可靠性的方法
CN108961638A (zh) * 2018-05-23 2018-12-07 重庆科技学院 基于小波系数能量和算法的振动光纤入侵事件检测方法
CN112541480B (zh) * 2020-12-25 2022-06-17 华中科技大学 一种隧道异物入侵事件在线识别方法及***
CN112907869B (zh) * 2021-03-17 2023-03-21 四川通信科研规划设计有限责任公司 基于多种传感技术的入侵检测***
CN116340834B (zh) * 2023-05-26 2023-10-03 成都凯天电子股份有限公司 一种被测量信号叠加直流偏置信号的数据修正***及方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096977A (zh) * 2010-11-26 2011-06-15 上海电力带电作业技术开发有限公司 一种入侵外力视频监控及预警的方法
CN103116957A (zh) * 2013-01-21 2013-05-22 复旦大学 一种光纤周界安防***屏蔽气候影响的方法
CN103617685A (zh) * 2013-12-09 2014-03-05 中国船舶重工集团公司第七〇五研究所 一种基于光纤振动传感报警***的自学习算法
CN104134303A (zh) * 2014-07-22 2014-11-05 上海光亮光电科技有限公司 光纤传感***的入侵信号识别方法
CN104240455A (zh) * 2014-08-07 2014-12-24 北京航天控制仪器研究所 一种分布式光纤管线安全预警***中的扰动事件识别方法
CN104517376A (zh) * 2014-12-15 2015-04-15 北京航天易联科技发展有限公司 一种两通道嵌入式安防装置
CN105023379A (zh) * 2015-08-13 2015-11-04 中国民航大学 一种机场光纤周界预警***的信号识别方法
CN105184319A (zh) * 2015-09-02 2015-12-23 同方威视技术股份有限公司 光纤周界入侵信号的识别方法、装置及周界入侵报警***

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITSV20040039A1 (it) * 2004-11-03 2005-02-03 Federico Ciordinik Dispositivo antieffrazione per recinzioni o simili

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096977A (zh) * 2010-11-26 2011-06-15 上海电力带电作业技术开发有限公司 一种入侵外力视频监控及预警的方法
CN103116957A (zh) * 2013-01-21 2013-05-22 复旦大学 一种光纤周界安防***屏蔽气候影响的方法
CN103617685A (zh) * 2013-12-09 2014-03-05 中国船舶重工集团公司第七〇五研究所 一种基于光纤振动传感报警***的自学习算法
CN104134303A (zh) * 2014-07-22 2014-11-05 上海光亮光电科技有限公司 光纤传感***的入侵信号识别方法
CN104240455A (zh) * 2014-08-07 2014-12-24 北京航天控制仪器研究所 一种分布式光纤管线安全预警***中的扰动事件识别方法
CN104517376A (zh) * 2014-12-15 2015-04-15 北京航天易联科技发展有限公司 一种两通道嵌入式安防装置
CN105023379A (zh) * 2015-08-13 2015-11-04 中国民航大学 一种机场光纤周界预警***的信号识别方法
CN105184319A (zh) * 2015-09-02 2015-12-23 同方威视技术股份有限公司 光纤周界入侵信号的识别方法、装置及周界入侵报警***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
光纤周界***入侵信号识别;杨尚文等;《中国科技论文在线》;20150910;第2页算法原理2.1短时过零率和短时能量 *
基于聚类分析的光周界入侵信号识别;刘阳等;《中国科技论文在线》;20151230;第1页第40-44行，第2页第1行，第5页第134-142行 *
聚类分析在分布式光纤振动传感***中的应用;唐文靖等;《激光技术》;20151130;第39卷(第6期);全文 *

Also Published As

Publication number Publication date
CN105654645A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN105654645B (zh) 一种光纤安防信号处理方法及***
CN106023499B (zh) 一种光纤安防信号双重识别方法及***
CN105608823B (zh) 基于主成分分析的光纤安防方法及***
James et al. Online false data injection attack detection with wavelet transform and deep neural networks
CN103236127B (zh) 一种光纤围栏入侵监测***及其模式识别方法
CN111178155B (zh) 基于惯性传感器的步态特征提取和步态识别方法
Cavaglia et al. Finding the origin of noise transients in LIGO data with machine learning
US10621506B2 (en) Apparatus and method for activity detection and classification from sensor data
Mahmoud et al. Robust event classification for a fiber optic perimeter intrusion detection system using level crossing features and artificial neural networks
CN109213121A (zh) 一种风机制动***卡缸故障诊断方法
CN109508736A (zh) 一种基于深度学习的监狱异常情况监测方法及监测***
KR102247188B1 (ko) 상황인식 및 퍼지이론 기반의 해수욕장 위험도 평가 시스템 및 방법과, 이를 위한 컴퓨터 프로그램
KR20170053692A (ko) 커널 회귀 모델의 앙상블을 위한 장치 및 방법
CN106600869A (zh) 一种用于光纤围栏安防***的围栏入侵识别方法
CN105067101A (zh) 振源识别的基于振动信号的基音频率特征的提取方法
CN106973039A (zh) 一种基于信息融合技术的网络安全态势感知模型训练方法和装置
WO2020060854A1 (en) Dynamic graph analysis
Xu et al. Seismic target recognition based on parallel recurrent neural network for unattended ground sensor systems
CN102346948B (zh) 一种周界入侵检测方法和***
CN112966248B (zh) 非受控步行场景下移动设备的持续身份认证方法
Sekar Power quality disturbance detection using machine learning algorithm
Ning et al. Fall detection algorithm based on gradient boosting decision tree
CN107689127B (zh) 具有自适应调参功能的光纤安防信号双重识别方法与***
CN107137090B (zh) 一种跌倒识别方法及装置、用户设备
Ma et al. An optic-fiber fence intrusion recognition system using mixture Gaussian hidden Markov models

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