发明名称 声磁电子商品防盗系统的检测方法
摘要 本发明公开了一种声磁电子商品防盗系统,包括CPU处理器,CPU处理器分别与接收电路模块、发射电路模块、模式选择电路模块、同步电路、声光报警模块、FLASH存储单元和电源;检测方法为:初始化进入系统主循环,进行数据采样,对采样数据运算后与标准标签信号特征进行比对,判断相似性是否满足报警条件进行报警。本发明声磁电子商品防盗系统,将多种数字信号处理算法有效的组合,使检测快速有效,还包含了整个系统的时序控制程序,保证系统各功能模块协同有序的工作,解决了现有电子商品防盗系统检测率低、会发生误报现象的问题,检测距离远,检测率高,误报率低。
申请公布号 CN104200597B 申请公布日期 2017.01.25
申请号 CN201410389856.2 申请日期 2014.08.08
申请人 西安工程大学 发明人 景军锋;赵卓;陈耀武;李鹏飞;屈肖军
分类号 G08B13/24(2006.01)I 主分类号 G08B13/24(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 罗笛
主权项 声磁电子商品防盗系统的检测方法,其特征在于,该方法基于声磁电子商品防盗系统,所述声磁电子商品防盗系统包括CPU处理器(1),所述CPU处理器(1)分别与接收电路模块(2)、发射电路模块(3)、模式选择电路模块(4)、同步电路(5)、声光报警模块(6)、FLASH存储单元(7)和电源(8)连接;利用所述声磁电子商品防盗系统进行检测,包括以下步骤:步骤1:启动声磁电子商品防盗系统,硬件初始化,对系统所需要的功能进行配置;步骤2:CPU处理器(1)控制发射使能,通过输出SPWM波形使发射电路模块(3)的发射天线向空间发射58KHz交变磁场,发射时间为T1;步骤3:待发射T1时间后关闭发射,延时△T1后打开接收电路模块(2)的接收天线使能控制,来侦测空间内的各类信号,信号经放大、滤波后进入CPU处理器(1)的ADC采样端;步骤4:启动CPU处理器(1)的ADC对信号进行采样,并保存两组数据至FLASH存储单元(7)中,其中一组数据不作处理直接保存;另一组将连续若干周期采样值做叠加处理保存,每一次采样的一组值可以理解为一帧数据,完成采样;步骤5:对所述步骤4直接保存的数据进行信号互相关运算及快速傅里叶变换:<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>&lsqb;</mo><mi>A</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>B</mi><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023996600000011.GIF" wi="1166" he="126" /></maths><maths num="0002"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>-</mo><mn>1</mn></mrow></munderover><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><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><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><mo>.</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023996600000021.GIF" wi="1406" he="199" /></maths>式中,R(n)表示信号互相关函数,N表示采样点数,m表示采样点索引,A(m)·B(m+n)表示信号A与信号B,F(k)为离散傅立叶变换函数,n表示采样点索引偏移量,X(n)表示采样点,<img file="FDA0001023996600000022.GIF" wi="246" he="95" />对经叠加处理的信号波形包络提取,求取包络峰值,在得到包络峰值的基础上利用最小二乘法求取其斜率,并进行波峰值方差运算:Y=X[i](△X[i]=0,X[i]≥0)    (3)<maths num="0003"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfrac><mrow><mi>N</mi><mi>&Sigma;</mi><mi>X</mi><mi>Y</mi><mo>-</mo><mi>&Sigma;</mi><mi>X</mi><mi>&Sigma;</mi><mi>Y</mi></mrow><mrow><msup><mi>N&Sigma;X</mi><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><mi>&Sigma;</mi><mi>X</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023996600000023.GIF" wi="1078" he="151" /></maths><maths num="0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>X</mi><mrow><mi>A</mi><mi>v</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><mi>X</mi><mo>&lsqb;</mo><mi>i</mi><mo>&rsqb;</mo></mrow><mi>N</mi></mfrac></mrow></mtd><mtd><mrow><mi>D</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>|</mo><mi>X</mi><mo>&lsqb;</mo><mi>i</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>X</mi><mrow><mi>A</mi><mi>v</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>|</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023996600000024.GIF" wi="1236" he="127" /></maths>式中,Y表示包络峰值,X[i]表示信号采样值,ΔX[i]表示采样信号值变化率,K表示信号斜率,Y<sub>Cnt</sub>表示包络波峰数,X即X[i]的简略表达形式,X<sub>Aver</sub>表示采样值均值,D表示波峰值方差;Y<sub>Cnt</sub>为满足公式<img file="FDA0001023996600000025.GIF" wi="640" he="143" />Y的个数;计算信号互相关满足次数R<sub>Cnt</sub>及频谱能量F<sub>Mag</sub>:设P为满足表达式<img file="FDA0001023996600000026.GIF" wi="854" he="239" />R<sub>n</sub>的个数,当P≥10认为满足一次互相关条件,R<sub>Cnt</sub>为满足P的次数;式中,R<sub>n</sub>表示中间变量,R[i]为信号互相关函数R(n)的计算值,R<sub>ref</sub>[i]表示互相关参照值,Δ(R[i]‑R<sub>ref</sub>[i])为互相关差值变化率;设幅频能量谱<img file="FDA0001023996600000027.GIF" wi="1342" he="141" />则58KHz频谱能量F<sub>Mag</sub>定义为F<sub>Mag</sub>=F<sub>M</sub>[15]+F<sub>M</sub>[16]+F<sub>M</sub>[17]+F<sub>M</sub>[18];式中,F[i]为离散傅里叶变换函数F(k)的计算值;步骤6:将所述步骤5的运算结果与标准标签信号特征进行比对,判断相似性是否满足条件,进行声光报警或进行下一帧判断;步骤7:重复步骤2‑6,进行下一周期判断。
地址 710048 陕西省西安市碑林区金花南路19号