CN109035662B - 一种基于电容扰动的周界报警系统入侵信号识别方法 - Google Patents
一种基于电容扰动的周界报警系统入侵信号识别方法 Download PDFInfo
- Publication number
- CN109035662B CN109035662B CN201811043680.XA CN201811043680A CN109035662B CN 109035662 B CN109035662 B CN 109035662B CN 201811043680 A CN201811043680 A CN 201811043680A CN 109035662 B CN109035662 B CN 109035662B
- Authority
- CN
- China
- Prior art keywords
- signal
- component
- disturbance
- prc
- obtaining
- 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
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000012549 training Methods 0.000 claims abstract description 36
- 238000013528 artificial neural network Methods 0.000 claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000009825 accumulation Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229910052709 silver Inorganic materials 0.000 description 2
- 239000004332 silver Substances 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B13/00—Burglar, theft or intruder alarms
- G08B13/22—Electrical actuation
- G08B13/26—Electrical actuation by proximity of an intruder causing variation in capacitance or inductance of a circuit
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于电容扰动的周界报警系统入侵信号识别方法，涉及电容扰动信号处理及识别技术。首先对获取的电容扰动信号进行去噪处理；然后使用固有时间尺度分解算法ITD和HT变换提取训练样本特征值；最后构建PNN概率神经网络，将有效的特征量供PNN神经网络训练，经过训练从而调节好PNN概率神经网络的参数，然后将调整好的参数交给概率神经网络分类模块使用；从而具有识别精度高、效率好、灵敏度高、安全性高的优点。
Description
技术领域
本发明涉及电容扰动信号处理及识别技术，尤其涉及一种电容扰动围界入侵信号的识别方法。
背景技术
随着社会的发展，机场围界报警系统是机场安全防范系统的重中之重。《民用运输机场安全保卫设施》(MH/T 7003-2017)明确规定年旅客吞吐量200万人次以上的一二类机场的飞行区围界应设置围界报警系统和视频监控系统，年旅客吞吐量在50万至200万人次的三类机场宜设置入侵报警系统和视频监控系统。系统应能对物理围界形成完整的覆盖，能对围界入侵行为作出判断，对目标进行分类，应能在机场飞行和安全照明的环境下工作，应满足全天候运行的要求。
电容周界报警系统是一种全新概念的周界报警入侵检测系统，集阻挡和报警两大功能于一体，可实现对指定区域的入侵报警、安全距离报警、24小时不间断实时边界守卫等功能，其可靠性高，实时性好，极大的缓解了保卫人员的工作强度，提高了防护区域的安全性，具有广泛的应用前景。
目前关于电容扰动周界报警的研究较少，专利文献CN104851226A于2015年8月19日公开了一种用于电容扰动周界报警系统和报警方法，只提供了机场耦合电容式报警系统，能满足机场入侵检测的需求，但还是存在一定的误报率和漏报率，现有系统中缺乏对信号的优化处理。如何确保系统对入侵信号辨识的准确性及如何提供一种高效、实时的入侵检测方法为电容围界报警系统的一个重要研究方向。
发明内容
有鉴于此，本发明提供了一种电容围界入侵信号处理方法，以满足机场范围长距离入侵检测需求，减小系统误警次数和平均漏报次数。
本发明技术方案为一种基于电容扰动的周界报警系统入侵信号识别方法，包括以下步骤：
步骤1、对获取的电容扰动信号进行去噪处理；
步骤2、使用固有时间尺度分解算法ITD(固有时间尺度分解算法)和HT变换提取训练样本特征值；
步骤3、构建PNN概率神经网络，将有效的特征量供PNN神经网络训练，经过训练从而调节好PNN概率神经网络的参数，然后将调整好的参数交给概率神经网络分类模块使用；
其中所述步骤1的具体方法为：
步骤1.1根据电容扰动信号和干扰频带经过实验确定分解层数；
步骤1.2采用多个小波滤波器对电容扰动信号进行小波包变换，选取效果最好的一个获得小波包变换系数；
步骤1.3根据原始扰动信号的特点和小波包基函数最佳选取算法，采用自底向顶的方法选取小波基函数；
步骤1.4对最佳正交小波包基对应的小波包系数进行处理；
步骤1.5重构去除噪声后的扰动信号；
所述步骤3的具体方法为：
步骤3.1对输入矩阵进行归一化处理，训练样本的最终特征向量矩阵为
x＝[f,a,t,He]
步骤3.2训练PNN概率神经网络；
步骤3.3采用训练好的PNN概率神经网络对实施信号进行分类识别。
进一步的，所述步骤1.2中小波滤波器包：Daubechies小波滤波器、Symlets小波滤波器。
进一步的，所述步骤2的具体方法为：
步骤2.1利用固有时间尺度分解算法提取电容扰动信号的固有旋转分量PRC和一个残余分量；
步骤2.2经ITD分解得到的n个PRC分量，每个分量都包含有其各自的频率成分f1,f2,...,fn，其中含有不同的能量E＝{E1,E2,...,En}；通过PRC分量的完备正交性来追踪原始信号能量与分解能量和之间的差，自适应确定真实的PRC分量；
步骤2.3通过KL散度方法来筛选PRC分量，各个PRC分量与原始信号的KL散度值越大差异性越大，反之越小；
步骤2.4进行Hilbert变换求取相位和瞬时频率，根据高频突变点得到扰动起始时刻和扰动终止时刻，由包络函数得到扰动信号4的幅值。
进一步的，所述步骤2.1的具体方法为：
步骤2.1.1确定原信号x(t)的x(t)及对应的时刻t(k)，其中t≥0，k＝1,2,3,...,M,M为极值点个数，并计算基线提取因子L为
其中,
式中k＝1,2,3,...,M-2,0＜α＜1，α是固有旋转分量幅度线性按比例提取的增益控制参数，一般α＝0.5
步骤2.1.2由表达式(1)可知Lk的值是从L2到LM-1，两端点L1、LM的值需要估计计算。在两边端点处分别向中间取五个点，根据三次多项式拟合求出两端点L1、LM的值
步骤2.1.3采用三次Hermite插值拟合所有的L1,L2,L3,...,,LM，得到基线信号L(t)【3】
步骤2.1.4将基线信号从原信号中分离出来，即
h1(t)＝x(t)-L(t)
若h1(t)满足PRC分量的条件，输出h1(t)。否则h1(t)作为原始数据，继续重复上述步骤，直到h1k(t)为PRC分量，记PRC1＝h1k(t)
步骤2.1.5将PRC1分量从原信号中分离出来，则
u1(t)＝x(t)-PRC1
得到新的剩余信号u1(t)
步骤2.1.6将u1(t)看作原始数据，重复以上步骤循环n-1次，直到剩余信号un(t)为单调信号或常信号
步骤2.1.7则信号分解为
进一步的，所述步骤2.2的具体方法为：
步骤2.2.1对原始信号求取信号能量
因为xi(t)相互之间是不相关的，即是正交的，有
所以信号x(t)具有的总能量为
式中，Ei为x(t)的第i个分量xi(t)的能量
步骤2.2.2当得到的一个PRC分量为各分量能量之和为x(t)的正交分量x1(t)，信号的总能量Enew与原始信号的能量Ex是相等的，当分解的PRC分量不是x(t)中的一个正交分量时，信号的总能量Enew与原始信号的能量Ex有一个误差值为
Er＝Enew-Ex
|Er|越小，越能保证信号分解分量的完备正交性，所以当|Er|达到最小值时，得到的PRC分量是原始信号中的一个正交分量
步骤2.2.3在ITD过程中，若能力误差|Er|超出允许范围
则修正算法。
进一步的，所述步骤2.3的具体方法为：
步骤2.3.1计算原信号X＝{x1，x2，...，xn}与ITD分解的第i个PRC分量信号Yi＝{y1,y2,...,yn}之间的KL散度值，设定该数量级为0.1，则KL散度值大于0.1的PRC为无效分量，等于或大于0.1的PRC为有效分量；
步骤2.3.2计算两个信号的概率分布，假设X和Yi的概率分布分别为Y和q(x)，采用非参数估计法求解概率分布，定义核密度函数p(x)为
式中，h为给定的正数，k为高斯核函数：
步骤2.3.3同理得出Yi的概率分布q(x)，得出X和Yi的KL距离δ(p,q)及δ(q,p)
步骤2.3.4同理得出δ(q,p)，即可得出X和Yi的KL散度值
Di(p,q)＝δ(p,q)+δ(q,p)
步骤2.3.5设定KL散度值平均值为阈值ε
其中，若某一PRC分量KL散度值Di(p,q)＜ε，则认为该分量包含了原始信号大部分有用的信息，可作为特征量；反之，若某一PRC分量的能量Di(p,q)＞ε，则认为该分量包含的原始信号有用信息较少，可将此分量除去。
进一步的，所述步骤2.4的具体方法为：
步骤2.4.1对x(t)进行Hilbert变换，则有
可得相应解析函数为
z(t)＝x(t)+jy(t)＝a(t)eθ(t)
其中相位函数为
瞬时幅值a(t)为
瞬时频率f(t)为
步骤2.4.2：由步骤2.3.1可得对各阶PRC分量的幅值进行归一化处理
可得ITD信息熵为
式中，pj是第j个PRC的幅值分量在总能量中所占的比重
步骤2.4.3从瞬时幅值曲线和瞬时频率曲线获得有效的扰动特征量，所得瞬时频率中可获得扰动期间频率，由瞬时频率曲线两频率突变点对应时刻的差值可得到扰动持续时间，从瞬时幅值可得扰动期间电压成分，所以训练样本的特征量为扰动频率f、扰动幅值a、扰动持续时间t、能量熵He。
进一步的，所述3.2的具体方法为：
步骤3.2.1训练样本数据集中的每个样本x表示为矢量的形式，即x＝(x1,x2,...,xd)
步骤3.2.2第一个样本矢量被置于输入层单元上，同时，连接输入单元和隐层单元的那些连接被初始化为w1＝x1，然后从隐层的第一个单元到累加层中代表x1所属类别的那个单元之间就建立了一个连接
步骤3.2.3同样的过程对剩余的各个隐层单元都重复进行，即wm＝xm，其中m＝1,2,...,n
步骤3.2.4最后得到训练后的PNN网络，输入层单元和隐层单元之间是万全相通的，而隐层单元到累加层单元之间是稀疏连接的。如果把第j类样本的第k个分量即为xjk，把这个分量到第j个隐层单元的连接权重系数为wjk，其中i＝1,2,...,n，k＝1,2,...,d。
进一步的，所述步骤3.3的具体方法为：
步骤3.3.1把一个测试样本提供给输入节点，每一个隐层单元都计算非线性函数
式中，i＝1,...,M，j＝1,...,N，M为训练样本中总类别数目，Ni为第i类训练样本数据数目，d为样本空间数据的维数，σ为平滑参数，xij为第i个模式的第j隐中心矢量
步骤3.3.2每一个累加单元把与它相连的隐层单元的计算fw(x)金子那个相加
步骤3.3.3对每一个累加层单元得到的计算结果，通过求最大值得到测试样本期望的类别，该类别包括：脚步、挖槽、攀爬、车辆过往、自然风、雨。
本发明首先对获取的电容扰动信号进行去噪处理；然后使用固有时间尺度分解算法ITD和HT变换提取训练样本特征值；最后构建PNN概率神经网络，将有效的特征量供PNN神经网络训练，经过训练从而调节好PNN概率神经网络的参数，然后将调整好的参数交给概率神经网络分类模块使用；从而具有识别精度高、效率好、灵敏度高、安全性高的优点。
附图说明
图1为电容扰动周界入侵信号识别方法的流程图；
图2为ITD算法改进图；
图3为识别电容周界入侵信号方法的流程图。
具体实施方式
本发明为一种基于电容扰动的周界报警系统入侵信号识别方法，包括以下步骤：
步骤1、对获取的电容扰动信号进行去噪处理；
步骤2、使用固有时间尺度分解算法ITD和HT变换提取训练样本特征值；
步骤3、构建PNN概率神经网络，将有效的特征量供PNN神经网络训练，经过训练从而调节好PNN概率神经网络的参数，然后将调整好的参数交给概率神经网络分类模块使用。
进一步的，所述步骤1的具体方法为：
步骤1.1根据电容扰动信号和干扰频带经过实验确定分解层数
步骤1.2选择适当的小波滤波器，对给定的扰动信号进行小波包变换，获得小波包变换系数。常用的正交滤波器包括Daubechies小波滤波器、Symlets小波滤波器等，经过实验比较选取最佳滤波器
步骤1.3根据原始扰动信号的特点和小波包基函数最佳选取算法，采用自底向顶的方法选取小波基函数
步骤1.4对最佳正交小波包基对应的小波包系数进行处理
步骤1.5重构去除噪声后的扰动信号
进一步的，所述步骤2的具体方法为：
步骤2.1利用固有时间尺度分解算法提取电容扰动信号的固有旋转分量PRC和一个残余分量；
步骤2.1.1确定原信号x(t)(t≥0)的x(t)及对应的时刻t(k)(k＝1,2,3,...,M,M为极值点个数)，并计算基线提取因子L为
其中,
式中k＝1,2,3,...,M-2,0＜α＜1，α是固有旋转分量幅度线性按比例提取的增益控制参数，一般α＝0.5
步骤2.1.2由表达式(1)可知Lk的值是从L2到LM-1，两端点L1、LM的值需要估计计算。在两边端点处分别向中间取五个点，根据三次多项式拟合求出两端点L1、LM的值
步骤2.1.3采用三次Hermite插值拟合所有的L1,L2,L3,...,,LM，得到基线信号L(t)
步骤2.1.4将基线信号从原信号中分离出来，即
h1(t)＝x(t)-L(t)
若h1(t)满足PRC分量的条件，输出h1(t)。否则h1(t)作为原始数据，继续重复上述步骤，直到h1k(t)为PRC分量，记PRC1＝h1k(t)
步骤2.1.5将PRC1分量从原信号中分离出来，则
u1(t)＝x(t)-PRC1
得到新的剩余信号u1(t)
步骤2.1.6将u1(t)看作原始数据，重复以上步骤循环n-1次，直到剩余信号un(t)为单调信号或常信号
步骤2.1.7则信号分解为
步骤2.2经ITD分解得到的n个PRC分量，每个分量都包含有其各自的频率成分f1,f2,...,fn，其中含有不同的能量E＝{E1,E2,...,En}。通过PRC分量的完备正交性来追踪原始信号能量与分解能量和之间的差【4】，自适应确定真实的PRC分量；
步骤2.2.1对原始信号求取信号能量
因为xi(t)相互之间是不相关的，即是正交的，有
所以信号x(t)具有的总能量为
式中，Ei为x(t)的第i个分量xi(t)的能量
步骤2.2.2当得到的一个PRC分量为各分量能量之和为x(t)的正交分量x1(t)，信号的总能量Enew与原始信号的能量Ex是相等的，当分解的PRC分量不是x(t)中的一个正交分量时，信号的总能量Enew与原始信号的能量Ex有一个误差值为
Er＝Enew-Ex
|Er|越小，越能保证信号分解分量的完备正交性，所以当|Er|达到最小值时，得到的PRC分量是原始信号中的一个正交分量
步骤2.2.3在ITD过程中，若能力误差|Er|超出允许范围
则修正算法
步骤2.3通过KL散度方法来筛选PRC分量，各个PRC分量与原始信号的KL散度值越大差异性越大，反之越小。
步骤2.3.1计算原信号X＝{x1，x2，...，xn}与ITD分解的第i个PRC分量信号Yi＝{y1,y2,...,yn}之间的KL散度值，设定该数量级为0.1，则KL散度值大于0.1的PRC为无效分量，等于或大于0.1的PRC为有效分量。
步骤2.3.2计算两个信号的概率分布，假设X和Yi的概率分布分别为Y和q(x)，采用非参数估计法求解概率分布，定义核密度函数p(x)为
式中，h为给定的正数，k为高斯核函数：
步骤2.3.3同理得出Yi的概率分布q(x)，得出X和Yi的KL距离δ(p,q)及δ(q,p)
步骤2.3.4同理得出δ(q,p)，即可得出X和Yi的KL散度值
Di(p,q)＝δ(p,q)+δ(q,p)
步骤2.3.5设定KL散度值平均值为阈值ε
其中，若某一PRC分量KL散度值Di(p,q)＜ε，则认为该分量包含了原始信号大部分有用的信息，可作为特征量；反之，若某一PRC分量的能量Di(p,q)＞ε，则认为该分量包含的原始信号有用信息较少，可将此分量除去
步骤2.4进行Hilbert变换求取相位和瞬时频率，根据高频突变点得到扰动起始时刻和扰动终止时刻，由包络函数得到扰动信号4的幅值
步骤2.4.1对x(t)进行Hilbert变换，则有
可得相应解析函数为
其中相位函数为
瞬时幅值a(t)为
瞬时频率f(t)为
步骤2.4.2：由步骤2.3.1可得对各阶PRC分量的幅值进行归一化处理
可得ITD信息熵为
式中，pj是第j个PRC的幅值分量在总能量中所占的比重
步骤2.4.3从瞬时幅值曲线和瞬时频率曲线获得有效的扰动特征量，所得瞬时频率中可获得扰动期间频率，由瞬时频率曲线两频率突变点对应时刻的差值可得到扰动持续时间，从瞬时幅值可得扰动期间电压成分，所以训练样本的特征量为扰动频率f、扰动幅值a、扰动持续时间t、能量熵He
进一步的，所述步骤3的具体方法为：
步骤3.1对输入矩阵进行归一化处理，训练样本的最终特征向量矩阵为
x＝[f,a,t,He]
步骤3.2PNN训练步骤
步骤3.2.1训练样本数据集中的每个样本x表示为矢量的形式，即x＝(x1,x2,...,xd)
步骤3.2.2第一个样本矢量被置于输入层单元上，同时，连接输入单元和隐层单元的那些连接被初始化为w1＝x1，然后从隐层的第一个单元到累加层中代表x1所属类别的那个单元之间就建立了一个连接
步骤3.2.3同样的过程对剩余的各个隐层单元都重复进行，即wm＝xm，其中m＝1,2,...,n步骤3.2.4最后得到训练后的PNN网络，输入层单元和隐层单元之间是万全相通的，而隐层单元到累加层单元之间是稀疏连接的。如果把第j类样本的第k个分量即为xjk，把这个分量到第j个隐层单元的连接权重系数为wjk，其中i＝1,2,...,n，k＝1,2,...,d
步骤3.3PNN分类识别步骤
步骤3.3.1把一个测试样本提供给输入节点，每一个隐层单元都计算非线性函数
式中，i＝1,...,M，j＝1,...,N，M为训练样本中总类别数目，Ni为第i类训练样本数据数目，d为样本空间数据的维数，σ为平滑参数，xij为第i个模式的第j隐中心矢量
步骤3.3.2每一个累加单元把与它相连的隐层单元的计算fw(x)金子那个相加
步骤3.3.3对每一个累加层单元得到的计算结果，通过求最大值得到测试样本期望的类别(脚步、挖槽、攀爬、车辆过往、自然风、雨)。
Claims (9)
1.一种基于电容扰动的周界报警系统入侵信号识别方法，包括以下步骤：
步骤1、对获取的电容扰动信号进行去噪处理；
步骤2、使用固有时间尺度分解算法ITD和HT变换提取训练样本特征值；
步骤3、构建PNN概率神经网络，将训练样本特征值供PNN神经网络训练，经过训练从而调节好PNN概率神经网络的参数，然后将调整好的参数交给概率神经网络分类模块使用；
其中所述步骤1的具体方法为：
步骤1.1根据电容扰动信号和干扰频带经过实验确定分解层数；
步骤1.2采用多个小波滤波器对电容扰动信号进行小波包变换，选取效果最好的一个获得小波包变换系数；
步骤1.3根据原始扰动信号的特点和小波包基函数最佳选取算法，采用自底向顶的方法选取小波基函数；
步骤1.4对最佳正交小波包基对应的小波包系数进行处理；
步骤1.5重构去除噪声后的扰动信号；
所述步骤3的具体方法为：
步骤3.1对步骤2得到的训练样本特征值进行归一化处理，训练样本的最终特征向量矩阵为
x＝[f,a,t,He]
步骤3.2训练PNN概率神经网络；
步骤3.3采用训练好的PNN概率神经网络对实施信号进行分类识别。
2.如权利要求1所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤1.2中小波滤波器包括：Daubechies小波滤波器、Symlets小波滤波器。
3.如权利要求1所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述3.2的具体方法为：
步骤3.2.1训练样本数据集中的每个样本x表示为矢量的形式，即
x＝(x1,x2,...,xd)
步骤3.2.2第一个样本矢量被置于输入层单元上，同时，连接输入单元和隐层单元的那些连接被初始化为w1＝x1，然后从隐层的第一个单元到累加层中代表x1所属类别的那个单元之间就建立了一个连接
步骤3.2.3同样的过程对剩余的各个隐层单元都重复进行，即wm＝xm，其中m＝1,2,...,n
步骤3.2.4最后得到训练后的PNN网络，输入层单元和隐层单元之间是万全相通的，而隐层单元到累加层单元之间是稀疏连接的；如果把第j类样本的第k个分量即为xjk，把这个分量到第j个隐层单元的连接权重系数为wjk，其中i＝1,2,...,n，k＝1,2,...,d。
4.如权利要求1所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤3.3的具体方法为：
步骤3.3.1把一个测试样本提供给输入节点，每一个隐层单元都计算非线性函数
式中，i＝1,...,M，j＝1,...,N，M为训练样本中总类别数目，Ni为第i类训练样本数据数目，d为样本空间数据的维数，σ为平滑参数，xij为第i个模式的第j隐中心矢量；
步骤3.3.2每一个累加单元把与它相连的隐层单元计算的fw(x)相加；
步骤3.3.3对每一个累加层单元得到的计算结果，通过求最大值得到测试样本期望的类别，该类别包括：脚步、挖槽、攀爬、车辆过往、自然风、雨。
5.如权利要求1或2所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤2的具体方法为：
步骤2.1利用固有时间尺度分解算法提取电容扰动信号的固有旋转分量PRC和一个残余分量；
步骤2.2经ITD分解得到的n个PRC分量，每个分量都包含有其各自的频率成分f1,f2,...,fn，其中含有不同的能量E＝{E1,E2,...,En}；通过PRC分量的完备正交性来追踪原始信号能量与分解能量和之间的差，自适应确定真实的PRC分量；
步骤2.3通过KL散度方法来筛选PRC分量，各个PRC分量与原始信号的KL散度值越大差异性越大，反之越小；
步骤2.4进行Hilbert变换求取相位和瞬时频率，根据高频突变点得到扰动起始时刻和扰动终止时刻，由包络函数得到扰动信号的幅值。
6.如权利要求5所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤2.1的具体方法为：
步骤2.1.1确定原信号x(t)的x(t)及对应的时刻t(k)，其中t≥0，k＝1,2,3,...,M,M为极值点个数，并计算基线提取因子L为
其中,
式中k＝1,2,3,...,M-2,0＜α＜1，α是固有旋转分量幅度线性按比例提取的增益控制参数，一般α＝0.5
步骤2.1.2由表达式(1)可知Lk的值是从L2到LM-1，两端点L1、LM的值需要估计计算；在两边端点处分别向中间取五个点，根据三次多项式拟合求出两端点L1、LM的值
步骤2.1.3采用三次Hermite插值拟合所有的L1,L2,L3,...,,LM，得到基线信号L(t)【3】
步骤2.1.4将基线信号从原信号中分离出来，即
h1(t)＝x(t)-L(t)
若h1(t)满足PRC分量的条件，输出h1(t)；否则h1(t)作为原始数据，继续重复上述步骤，直到h1k(t)为PRC分量，记PRC1＝h1k(t)
步骤2.1.5将PRC1分量从原信号中分离出来，则
u1(t)＝x(t)-PRC1
得到新的剩余信号u1(t)
步骤2.1.6将u1(t)看作原始数据，重复以上步骤循环n-1次，直到剩余信号un(t)为单调信号或常信号
步骤2.1.7则信号分解为
7.如权利要求5所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤2.2的具体方法为：
步骤2.2.1对原始信号求取信号能量
因为xi(t)相互之间是不相关的，即是正交的，有
所以信号x(t)具有的总能量为
式中，Ei为x(t)的第i个分量xi(t)的能量
步骤2.2.2当得到的一个PRC分量为各分量能量之和为x(t)的正交分量x1(t)，信号的总能量Enew与原始信号的能量Ex是相等的，当分解的PRC分量不是x(t)中的一个正交分量时，信号的总能量Enew与原始信号的能量Ex有一个误差值为
Er＝Enew-Ex
|Er|越小，越能保证信号分解分量的完备正交性，所以当|Er|达到最小值时，得到的PRC分量是原始信号中的一个正交分量
步骤2.2.3在ITD过程中，若能力误差|Er|超出允许范围
则修正算法。
8.如权利要求5所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤2.3的具体方法为：
步骤2.3.1计算原信号X＝{x1，x2，...，xn}与ITD分解的第i个PRC分量信号Yi＝{y1,y2,...,yn}之间的KL散度值，设定该数量级为0.1，则KL散度值大于0.1的PRC为无效分量，等于或大于0.1的PRC为有效分量；
步骤2.3.2计算两个信号的概率分布，假设X和Yi的概率分布分别为Y和q(x)，采用非参数估计法求解概率分布，定义核密度函数p(x)为
式中，h为给定的正数，k为高斯核函数：
步骤2.3.3同理得出Yi的概率分布q(x)，得出X和Yi的KL距离δ(p,q)及δ(q,p)
步骤2.3.4同理得出δ(q,p)，即可得出X和Yi的KL散度值
Di(p,q)＝δ(p,q)+δ(q,p)
步骤2.3.5设定KL散度值平均值为阈值ε
其中，若某一PRC分量KL散度值Di(p,q)＜ε，则认为该分量包含了原始信号大部分有用的信息，可作为特征量；反之，若某一PRC分量的能量Di(p,q)＞ε，则认为该分量包含的原始信号有用信息较少，可将此分量除去。
9.如权利要求5所述的一种基于电容扰动的周界报警系统入侵信号识别方法，其特征在于所述步骤2.4的具体方法为：
步骤2.4.1对x(t)进行Hilbert变换，则有
可得相应解析函数为
z(t)＝x(t)+jy(t)＝a(t)eθ(t)
其中相位函数为
瞬时幅值a(t)为
瞬时频率f(t)为
步骤2.4.2：由步骤2.3.1可得对各阶PRC分量的幅值进行归一化处理
可得ITD信息熵为
式中，pj是第j个PRC的幅值分量在总能量中所占的比重
步骤2.4.3从瞬时幅值曲线和瞬时频率曲线获得有效的扰动特征量，所得瞬时频率中可获得扰动期间频率，由瞬时频率曲线两频率突变点对应时刻的差值可得到扰动持续时间，从瞬时幅值可得扰动期间电压成分，所以训练样本的特征量为扰动频率f、扰动幅值a、扰动持续时间t、能量熵He。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811043680.XA CN109035662B (zh) | 2018-09-07 | 2018-09-07 | 一种基于电容扰动的周界报警系统入侵信号识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811043680.XA CN109035662B (zh) | 2018-09-07 | 2018-09-07 | 一种基于电容扰动的周界报警系统入侵信号识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109035662A CN109035662A (zh) | 2018-12-18 |
| CN109035662B true CN109035662B (zh) | 2020-08-11 |
Family
ID=64624489
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811043680.XA Active CN109035662B (zh) | 2018-09-07 | 2018-09-07 | 一种基于电容扰动的周界报警系统入侵信号识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109035662B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110120140B (zh) * | 2019-04-16 | 2020-12-29 | 电子科技大学 | 一种融合电容扰动信号和视频信号的复合式围界报警方法 |
| EP3731203B1 (en) | 2019-04-24 | 2023-05-31 | Carrier Corporation | Alarm system |
| CN113899976B (zh) * | 2021-10-30 | 2024-03-29 | 福州大学 | 一种复合电能质量扰动可视化方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104851226A (zh) * | 2015-05-11 | 2015-08-19 | 中国民用航空总局第二研究所 | 一种电容扰动围界报警系统和报警方法 |
| CN105266801A (zh) * | 2015-09-10 | 2016-01-27 | 太原理工大学 | 基于概率神经网络的j波检测及分类方法 |
| CN206470414U (zh) * | 2017-01-13 | 2017-09-05 | 民航成都电子技术有限责任公司 | 一种用于电容扰动围界报警的信号处理装置 |
| CN206628048U (zh) * | 2017-01-13 | 2017-11-10 | 民航成都电子技术有限责任公司 | 一种具备自恢复功能的电容扰动围界报警信号处理电路 |
| CN108152363A (zh) * | 2017-12-21 | 2018-06-12 | 北京工业大学 | 一种抑端固有时间尺度分解的管道缺陷识别方法 |
-
2018
- 2018-09-07 CN CN201811043680.XA patent/CN109035662B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104851226A (zh) * | 2015-05-11 | 2015-08-19 | 中国民用航空总局第二研究所 | 一种电容扰动围界报警系统和报警方法 |
| CN105266801A (zh) * | 2015-09-10 | 2016-01-27 | 太原理工大学 | 基于概率神经网络的j波检测及分类方法 |
| CN206470414U (zh) * | 2017-01-13 | 2017-09-05 | 民航成都电子技术有限责任公司 | 一种用于电容扰动围界报警的信号处理装置 |
| CN206628048U (zh) * | 2017-01-13 | 2017-11-10 | 民航成都电子技术有限责任公司 | 一种具备自恢复功能的电容扰动围界报警信号处理电路 |
| CN108152363A (zh) * | 2017-12-21 | 2018-06-12 | 北京工业大学 | 一种抑端固有时间尺度分解的管道缺陷识别方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于ITD样本熵与PNN的自动机故障诊断;赵雄鹏等;《机械设计与研究》;20170630;第33卷(第3期);138-141,146 * |
| 基于小波包变换的弱信号检测;童宁宁等;《计算机仿真》;20060731;第23卷(第7期);105-107,168 * |
| 基于小波变换的特征提取方法分析;张静远等;《信号处理》;20000630;第16卷(第2期);156-162,155 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109035662A (zh) | 2018-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106874833B (zh) | 一种振动事件的模式识别方法 | |
| CN114201988B (zh) | 卫星导航复合干扰信号识别方法和系统 | |
| CN109035662B (zh) | 一种基于电容扰动的周界报警系统入侵信号识别方法 | |
| Mishra et al. | Land cover classification of PALSAR images by knowledge based decision tree classifier and supervised classifiers based on SAR observables | |
| CN107103306B (zh) | 基于小波分析与支持向量机的冬小麦白粉病遥感监测方法 | |
| Misák et al. | A complex classification approach of partial discharges from covered conductors in real environment | |
| CN102722892A (zh) | 基于低秩矩阵分解的sar图像变化检测方法 | |
| CN105701470A (zh) | 一种基于最优小波包分解的模拟电路故障特征提取方法 | |
| CN105989597B (zh) | 基于像素选择过程的高光谱图像异常目标检测方法 | |
| Chen et al. | Investigating rice cropping practices and growing areas from MODIS data using empirical mode decomposition and support vector machines | |
| WO2020060830A1 (en) | Temporal behavior analysis of network traffic | |
| CN103366184A (zh) | 基于混合分类器的极化sar数据分类方法及系统 | |
| Cantzos et al. | Identifying long-memory trends in pre-seismic MHz Disturbances through Support Vector Machines | |
| CN103295583B (zh) | 用于提取声音的子带能量特征的方法、设备以及监视系统 | |
| Liu et al. | Graph signal coarsening: Dimensionality reduction in irregular domain | |
| Khoje et al. | Comparative performance evaluation of fast discrete curvelet transform and colour texture moments as texture features for fruit skin damage detection | |
| CN113780308B (zh) | 基于核主元分析与神经网络的gis局部放电模式识别方法及系统 | |
| Pulli | Extracting and processing signal parameters for regional seismic event identification | |
| Välisuo | Automated wind turbine noise analysis by machine learning | |
| Nemalidinne et al. | Nonsubsampled Contourlet Domain Fusion Approach for Infrared and Visible Fire Images | |
| Gao | Research on modulation recognition method of electromagnetic signal based on wavelet transform convolutional neural network | |
| CN106225846A (zh) | 大棚监测系统 | |
| Kanirajan et al. | Recognition of Power Quality Disturbances using Fuzzy Expert Systems | |
| Zadeh et al. | Discrimination between earthquakes and explosions using artificial neural networks | |
| Ibebuchi et al. | Application of autoencoders artificial neural network and principal component analysis for pattern extraction and spatial regionalization of global temperature data |
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 |