CN107044893B - 一种基于时域编解码的高精度动态信号传感和传输方法 - Google Patents

一种基于时域编解码的高精度动态信号传感和传输方法 Download PDF

Info

Publication number
CN107044893B
CN107044893B CN201710152542.4A CN201710152542A CN107044893B CN 107044893 B CN107044893 B CN 107044893B CN 201710152542 A CN201710152542 A CN 201710152542A CN 107044893 B CN107044893 B CN 107044893B
Authority
CN
China
Prior art keywords
signal
sensor
nlari
transmission
sensing
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
CN201710152542.4A
Other languages
English (en)
Other versions
CN107044893A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201710152542.4A priority Critical patent/CN107044893B/zh
Publication of CN107044893A publication Critical patent/CN107044893A/zh
Application granted granted Critical
Publication of CN107044893B publication Critical patent/CN107044893B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • G01L1/16Measuring force or stress, in general using properties of piezoelectric devices
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405Determining heart rate variability
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Cardiology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Physiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)

Abstract

本发明公开了一种基于时域编解码的高精度动态信号传感和传输方法,使用非线性自回归整合NLARI过程建模随机弹性的传感器对动态信号的响应,通过选取满足NLARI模型稳定不动点条件的传感材料,从而使传感器功能相当于一个稳定不变的传递函数;通过模拟神经元的动作电压全有全无和反向传播及局部电位可叠加的调节功能,导入调幅逆向脉冲传递函数实现信号重构;以及输出入相似度和信号传输成功率评价标准。本发明在多数情形下从理论上防止了系统误差;对尖峰脉冲序列能够实现信号重构;由于是时域上编解码,所以无须频率分析;按本发明制成的传感器无须非线性动态补偿;本发明兼有准确、响应快、精度高、成本低的优点。

Description

一种基于时域编解码的高精度动态信号传感和传输方法
技术领域
本发明属于信号传感和传输技术领域,尤其涉及一种基于时域编解码的高精度动态信号传感和传输方法。
背景技术
高效精准传感或传输动态信号,尤其是生理时间序列信号和尖峰时间序列信号,可望在大脑及人工神经网络的信息处理、脉冲雷达、通讯领域发挥重大作用。在其他领域,比如生产过程的监视和控制中,也有如何保证尖峰时间序列信号在互相传递和输送过程不失真的问题。
减小系统误差是大幅提高传感和传输精准度的关键。非线性脉冲动态信号传感通常使用静态标定和动态补偿相结合、时域方法测定动态特性辨识动态特性曲线等方式,它们可以降低多种偶然性造成的随机误差或反常误差。然而,作为动态特性补偿的频率域修正法、数值微分法、叠加积分法未考虑误差对处理结果的稳定性。反卷积补偿法、单参数滤波法、多参数模型反滤波法不能采用递归实现,运算量极大。基于虎克定律的线性弹性动力系统或多项式非线性响应系统,不适合外界扰动尤其较大扰补偿动的情形。Hammerstein模型采用静态非线性模块和动态线性模块串联降低了问题的复杂性,非线性部分辨识问题转化为参数空间上的函数优化问题,或通过接口送入计算机由软件系统辨识得到动态模型和动态性能指标。但是这些方法带来的系统误差需要大量计算或优异功能的特殊材料来补偿,所以难以兼顾准确、快速、高精度且成本低廉的要求。
近年神经科学、非线性动力学、及时间序列分析的进步使时域编解码生理信号和尖峰脉冲信号成为可能。一类将牛顿运动第二定律用于非线性随机弹性系统导出的非线性自回归整合(NLARI)过程,可以被特定为如下时间序列模型:
Figure BSA0000141758170000021
其中
Figure BSA0000141758170000022
Figure BSA0000141758170000023
在那里ω代表外部扰动的期待值,vt代表外部扰动的方差为σ2的高斯白噪音;α>0是阻力系数,β>0是恢复力系数,γ>0是相对恢复力系数,κ1和κ2分别为阻力和恢复力的时间延迟;-xexp(x-2)是一个恢复力项g(x)。考虑κ1=κ2=1,已经证明μt=X0+(ω/α)t。让Yt=Xtt,方程(1)可以改写为如下方程:
Figure BSA0000141758170000024
(2)当σ=0或vt=0,方程(2)为确定性系统
yt=(1+θ1)yt-11yt-22[-yt-1exp(-y2 t-1)] (3)
已经证明它的稳定性和波动型决于相对恢复力系数γ:当0<γ<1,该系统是一个渐近稳定的零不动点
Figure BSA0000141758170000025
所谓不动点指yt=yt-1=…=y0;当
Figure BSA0000141758170000026
该系统是一个渐近稳定周期环
Figure BSA0000141758170000027
Figure BSA0000141758170000028
这个周期环丧失稳定性。
现有动态信号传感和传输中存在尖峰脉冲信号在互相传递和输送过程失真的问题以及系统误差需要大量计算或优异功能的特殊材料来补偿,难以兼顾准确、快速、高精度且成本低廉的要求。
发明内容
为了解决上述问题,本发明的目的在于提供一种基于时域编解码的高精度动态信号传感和传输方法。
本发明是这样实现的,一种基于时域编解码的高精度动态信号传感和传输方法,所述基于时域编解码的高精度动态信号传感和传输方法使用非线性自回归整合NLARI过程建模弹性传感器的动态响应,通过选取满足NLARI模型的稳定不动点条件的传感材料使传感器功能相当于一个稳定的不动点传递函数,不动点指它的确定系统方程(3)在时间t的值yt始终等于初值即yt=y0,这意味着初值为零即y0=0以及随机误差为零时,NLARI系统的响应输出Yt等于外生输入vt即Yt=vt,通常情况下vt包括了输入信号和随机误差;通过模拟神经元的动作电压全有全无和反向传播及局部电位可叠加的调节功能,导入调幅逆向脉冲传递函数实现信号重构;以及输入输出相似度(或均方误差等)和信号传输成功率的评价标准。
进一步,稳定不动点域上的传递函数包括:
Yt=F(vt):Yi-1,t=(2-α)Yi-1,t-1-(1-α)Yi-1,t-2+βg(Yi-1,t-1)+vi-1,t
以及逆向脉冲传递函数vt=F-1(Yt):作为无调幅vi,t=Yi-1,t和作为调幅
Figure BSA0000141758170000031
作为i=1,2,…,q。
进一步,所述的传递函数Yt=F(vt),包括有界且满足g(-x)=-g(x)和xg(x)<0的 恢复力项g(x),阻力系数0<α<2,恢复力系数β>0以及相对恢复力系数
Figure BSA0000141758170000032
所述的传递函数Yt=F(vt)和逆向传递函数vt=F-1(Yt),包括输入信号vi,t,响应输出Yi,t,Yi,t的一阶滞后值Yi,t-1,Yi,t的二阶滞后值Yi,t-2,初始值Yi,-1=0,Yi,0=0,v1,t=c1如果vt≥c1否则
Figure BSA00001417581700000410
vt是原始输入信号,
Figure BSA0000141758170000041
是具有小方差的白噪音,阀值0<c3<c2<c1,i=1,2,…,q表示传输中继节点数及t=1,2,…,n表示时间或信号数;
进一步,输出输入相似度和信号传输成功率作为传感和传输的评价标准中,信号传输成功率
Figure BSA0000141758170000042
作为尖峰脉冲信号传输的评价标准,其中Wq=(wq,1,…,wq,n)′为经过q次中继节点后的输出信号,W0=(w0,1,…,w0,n)′为初始输入信号,w0,t=1如果vt≥c3否则w0,t=0,wq,t=1如果Yq,t≥c3否则wq,t=0,vt是原始输入信号,当q=1,信号传输成功率为信号传感成功率。
进一步,所述的传感方法,具体的实现步骤包括:
选取能够代表被测量信号特征的时间序列作为输入信号v=(v1,…,vn)′,在计算机上求出参数仿真域,它位于参数域(0,2)作为α,(0,4)作为β,及(0,1)作为γ的某一区域,在那里NLARI稳定不动点传递函数的输出Y=(Y1,…,Yn)′和输入v有最大相似度或信号传感成功率;
选取压电材料j=1,…,m制作m个NLARI-传感器样品,输入信号v=(v1,…,vn)′到所述传感器样品,测出响应输出
Figure BSA0000141758170000043
让ΔYt=Yt-Yt-1,使用
Figure BSA0000141758170000044
对回归方程ΔYt=θ1,jΔYt-12,jg(Yt-1)+vt作最小2乘法估计,获得估计值
Figure BSA0000141758170000045
Figure BSA0000141758170000046
以及
Figure BSA0000141758170000047
Figure BSA0000141758170000048
如果
Figure BSA0000141758170000049
在仿真参数域内,计算输出输入相似度或信号传感成功率,选取具有最大相似度或传感成功率的传感器样品为NLARI-传感器,否则不作为NLARI-传感器样本。
进一步,传输方法,包括计算机传输,NLARI-传感器对接收到的输入信号作出响应输出,所述的响应输出经数据化后通过接口送入计算机完成传输;还包括NLARI-传感器传输,其具体的实现步骤包括:
步骤一,制作信号传感成功率最大的前p个NLARI(α(k),β(k))-传感器样本作为k=1,…,p,使用软件或硬件制作具有的调幅逆向脉冲传递函数功能的调幅器样本,对由NLARI(α(k),β(k))-传感器样本的输出作为调幅器样本的输入的串联系统,计算输出v2,t与输入v1,t的信号传感成功率,通过调节阈值c1,c2,c3选取具有最大的信号传感成功率的调幅器样品作为神经元-调幅器样品,与NLARI(α(k),β(k))-传感器样本合成为q个相同的NLARI(α(k),β(k))-传输器样本;
步骤二,让能够代表被测量信号特征的时间序列信号作为原始输入信号vt进入由q个相同NLARI(α(k),β(k))-传输器样品按前一个的输出为后一个输入的方式串联的传输系统,测量它由vt引起的最终输出;
步骤三,计算每次信号传输成功率,重复步骤二,计算信号平均信号传输成功率;
步骤四,选取具有最大平均信号传输成功率的NLARI(α(k),β(k))-传输器样本对k=1,…,p作为NLARI-传输器。
进一步,传输方法的步骤二,让能够代表被测量信号特征的时间序列信号作为原始输入信号vt进入由q个相同NLARI(α(k),β(k))-传输器样品按前一个的输出为后一个输入的方式串联的传输系统,测量它由vt引起的最终输出,具体包括:
选取能够代表被测量信号特征的时间序列作为输入信号vt,在计算机上仿真求出参数仿真域,它在参数域(0,2)作为α,(0,4)作为β和(0,1)作为γ域内,在那里NLARI稳定不动点传递函数Yt=F(vt|α,β,γ)的输出Y=(Y1,…,Yn)′和输入v=(v1,…,vn)′有最大的传感成功率(尖峰脉冲输入)或相似度(非尖峰脉冲输入)让j=1;
使用压电材料j=1,…,m制作m个传感器样品;
输入信号vt到每个传感器样品中,测出响应
Figure BSA0000141758170000061
对ΔYt=θ1,jΔYt-12,jg(Yt-1)+vt,使用数据
Figure BSA0000141758170000062
作最小2乘法估计
Figure BSA0000141758170000063
在那里ΔYt=Yt-Yt-1,进一步获得
Figure BSA0000141758170000064
作为j=1,…,m;
如果参数估计值
Figure BSA0000141758170000065
不在仿真参数域内,该样本不作为NLARI-传感器样本,否则计算响应输出
Figure BSA0000141758170000066
与输入v的相似度或者信号传感成功率(j=1,…,m)。
按照本发明提出的技术方案制作的NLARI-传感器具有一个稳定不动点的传递函数功能,就是说如果系统初值和随机误差都为零的话,NLARI-传感器的响应输出Yt等于信号输入vt即Yt=vt。所以在理论上NLARI-传感器响应输出与信号输入之间只有随机误差,没有确定性结构意义上的系统误差。这样制作的传感器带来的益处是显而易见的,即只要消除随机误差就可完成信号重构。仿真显示NLARI-传感器的精度比目前新型的时删传感器在经过动态测量误差补偿后的精度高出很多(详细情形将在具体实施方式一节里描述)。本发明受到大脑神经元传输信号的相关功能的启迪,建立起调幅逆向脉冲传递函数用以消除随机误差,实现对来自NLARI-传感器的响应信号重构。在此基础上构成的NLARI-传输器仿真三千次尖峰序列信号中继,重复一万次的平均信号传输成功率在稳定不动点的某些域内可以达到98%以上的水平,说明NLARI-传感器具有可靠而精确快速的传输性能。由于本发明在理论上确保了多数情形下没有确定性结构意义上的系统误差,所以无须进行频率分析或优化计算等繁琐复杂的非线性动态特性补偿。而且本发明在时域上编解码,因此简单直观,并兼有准确可靠、速度快、精度高、成本低的优点。
附图说明
图1是本发明实施例提供的用NLARI-传感器在不调节和全有全无调节下经多次传输信号失真的仿真例;
图2是是本发明实施例提供的的NLARI-传感器的设计框图;
图3是本发明实施例提供的图2所述的NLARI-传感器对临床异常心率的传感仿真例图;
图4是本发明实施例提供的图2所述的NLARI-传感器对临床正常与异常心率的传感仿真例图;
图5是本发明实施例提供的针对尖峰脉冲时间序列信号的NLARI-传输器的设计框图;
图6是本发明实施例提供的图5所述的NLARI-传输器的传输整体效果分步效果及平均传输成功率的仿真例图;
图7是本发明实施例提供的图5所述的NLARI-传输器对兴奋性神经元被药物抑制后发射和兴奋性神经元发射的尖峰脉冲序列信号的传输仿真例图。
图8是本发明实施例提供的传感技术与现有传感技术精度比较的一个例子。
图9是本发明实施例提供的基于时域编解码的高精度动态信号传感和传输方法的原理示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例提供的基于时域编解码的高精度动态信号传感和传输方法,使用非线性自回归整合NLARI过程建模具有随机弹性的传感材料的,稳定不动点域上的NLARI传感模型相当于一个传递函数;通过模拟神经元的动作电压全有全无和反向传播及局部电位可叠加的调节功能,导入调幅逆向脉冲传递函数实现信号重构;以及输出入相似度和信号传输成功率评价标准。
稳定不动点域上的传递函数包括:
Yt=F(vt):Yi-1,t=(2-α)Yi-1,t-1-(1-α)Yi-1,t-2+βg(Yi-1,t-1)+vi-1,t
以及逆向脉冲传递函数vt=F-1(Yt):作为无调幅vi,t=Yi-1,t和作为调幅
Figure BSA0000141758170000081
作为i=1,2,…,q。
所述的传递函数Yt=F(vt),包括有界且满足g(-x)=-g(x)和xg(x)<0的恢复力项g(x),阻力系数0<α<2,恢复力系数β>0以及相对恢复力系数0<γ=β/(4-2α)<1;
所述的传递函数Yt=F(vt)和逆向传递函数vt=F-1(Yt),包括输入信号vi,t,响应输出Yi,t,Yi,t的一阶滞后值Yi,t-1,Yi,t的二阶滞后值Yi,t-2,初始值Yi,-1=0,Yi,0=0,v1,t=c1如果vt≥c1否则
Figure BSA0000141758170000083
vt是原始输入信号,
Figure BSA0000141758170000084
是具有小方差的白噪音,阀值0<c3<c2<c1,i=1,2,…,q表示传输中继节点数及t=1,2,…,n表示时间或信号数;
相似度或输出输入误差的标准差作为传感和传输的评价标准中,信号传感成功率
Figure BSA0000141758170000082
作为尖峰脉冲信号传输的评价标准,其中Wq=(wq,1,…,wq,n)′为经过q次中继节点后的输出信号,W0=(w0,1,…,w0,n)′为初始输入信号,w0,t=1如果vt≥c3否则w0,t=0,wq,t=1如果Yq,t≥c3否则wq,t=0,vt是原始输入信号,当q>1,信号传感成功率称为信号传输成功率,数据传输领域的误码率作为信号传输成功率。
所述的传感方法,具体的实现步骤包括:
选取能够代表被测量信号特征的时间序列作为输入信号v=(v1,…,vn)′,在计算机上求出参数仿真域,它位于参数域(0,2)作为α,(0,4)作为β,及(0,1)作为γ的某一区域,在那里NLARI稳定不动点传递函数的输出Y=(Y1,…,Yn)′和输入v有最大相似度或信号传感成功率;
选取压电材料j=1,…,m制作m个NLARI-传感器样品,输入信号v=(v1,…,vn)′到所述传感器样品,测出响应输出
Figure BSA0000141758170000091
让ΔYt=Yt-Yt-1,使用
Figure BSA0000141758170000092
对回归方程ΔYt=θ1,jΔYt-12,jg(Yt-1)+vt作最小2乘法估计,获得估计值
Figure BSA0000141758170000093
Figure BSA0000141758170000094
以及
Figure BSA0000141758170000095
Figure BSA0000141758170000096
如果
Figure BSA0000141758170000097
在仿真参数域内,计算输出输入相似度或信号传感成功率,选取具有最大相似度或传感成功率的传感器样品为NLARI-传感器,否则不作为NLARI-传感器样本。
下面结合附图对本发明的应用原理作详细描述。
本发明公开一种时域编解码的高精度动态信号传感和传输方法,它是基于一个非拉普拉斯变换的非线性动态系统传递函数来描述弹性传感器的输入与输出的动态响应关系,以及一个受神经元传输信号功能启发的调幅逆向脉冲传递函数来实现信号重构。NLARI的确定结构的稳定性是局部渐进稳定性,尽管不是全局渐进稳定,然而这种稳定性能够满足大多数生物信号和尖峰脉冲信号的无失真传感和传输,尤其通过数据缩小化处理后,效果会更好。本发明通过蒙特卡罗实验发现NLARI稳定不动点传递函数的产出并不拘泥于恢复力项g(x)在结构上的细节,只要满足在x轴上有界,并且xg(x)<0和g(-x)=-g(x),更为简单的函数可以替代-x/exp(x2)。
因此,本发明导入NLARI稳定不动点传递函数:
Yt=F(vt):Yt=(2-α)Yt-1-(1-α)Yt-2+βg(Yt-1)+vt
其中,vt代表输入信号,Yt代表传递函数的输出,Yt-1是Yt的一阶滞后值,Yt-2是Yt的二阶滞后值,由于Yt=Xt-X0-(ω/α)t,所以Yt表示信号与平均水平或均衡状态的偏离,阻力系数0<α<2,恢复力系数β>0,相对恢复系数γ=β/(4-2α)∈(0,1)及κ1=κ2=1。
本发明提出NLARI-传感器的设计方法基于使用NLARI过程建模弹性传感器的动态响应。比如压电材料在受到外机械力作用下会发生电极化,压电体两端表面内将出现符号相反的束缚电荷。当外力撤去后,它又恢复到不带电的状态。当外力作用方向改变时,电荷的极性随之改变。受力所产生的电荷量与外力的大小成正比。反之,压电材料在受到外电场作用下会发生形变,形变量与外电场强度成正比,表现出与正压电效应类似的现象,不同的是从电能向机械能的转换,压电陶瓷就具有这种机械能与电能之间的转换和逆转换的功能。由此可见,压电材料的正电效应和逆电效应在功能上表现为一个弹性系统,这样制成的传感器功能可以由NLARI模型描述。如果选取压电材料的阻力系数α和恢复力系数β位于NLARI的稳定不动点域0<γ<1内的,这种传感器的输出输入间的非线性动态响应关系就可以通过稳定不动点上的NLARI模型描述,称为NLARI稳定不动点传递函数,其中α,β和γ称为传感功能参数。考虑由m个NLARI-传感器样本作为中继节点,按上一个节点输出是下一个节点输入的串联方式作信号传递,动态响应关系可描述为Yi,t=(2-α)Yi,t-1-(1-β)Yi,t-2+βg(Yi,t-1)+Yi-1,t,i表示中继次数(i=1,2,…,q)和t表示信号数(t=1,2,…,n)。
本发明者注意到在不调节和全有全无调节下,经由NLARI-传感器样本多次作用后,尖峰脉冲输入信号会发生丢失或歪曲。
如图1所示,初始输入信号序列(如图1的101所示),经过3次(如图1的102所示)、10次(如图1的103所示)和17次(如图1的104所示)NLARI-传感器的信号中继仿真后,信号基本上保真,但在经过18次信号中继仿真后信号完全被歪曲(如图1的105所示),即使经过全有或全无律调节,在经过3次(如图1的106所示)、10次(如图1的107所示)和17次(如图1的108所示)和18次(如图1的109)NLARI-传感器的信号中继仿真后,信号也被丢失。
本发明使用神经科学的知识来纠正这种信号失真。来自外界的信息被编码成为神经状态时,相关的神经状态可能反映了神经元的解码活动。比如电压门控离子通道开放产生的电流要比刺激输入的电流大;电压刺激随着突触距离和神经递质结合时间的增大呈指数性衰减;受到刺激产生兴奋的轴突与周围静息膜之间的局部电流向两个方向传导(被称之为动作电位的双向传导),产生于轴突或胞体的动作电压向树突丛逆向回传的传播程度随着动作电位的数量和频率的增加而增加,大小取决于阈下兴奋性输入。动作电位反向传播活性可能促进先前冗余电流补充当前被衰减的电压刺激,旨在让信号复原。局部电位的可叠加性使动作电位的反向传播功能得以实现。可以总结如下:动作电位的反向传播活性和全有全无律、以及局部电位的可叠加性暗示了神经元可能具有功能(i)高于阈值的任何强度对刺激神经元响应是相同的;(ii)没有达到阈值的神经元不会点火;(iii)先前冗余电流的回流补充当前衰减的信号电流以促使信号重构。功能(i)源于生物的有效性和节能性;功能(ii)防止噪音混入信号,功能(iii)避免信号衰减。因此,针对尖峰脉冲信号重构,导入逆向脉冲传递函数vt=F-1(Yt):作为少数几个中继节点的传输,它可以是无调节函数F-1=I,即vi,t=Yi-1,t,初始值Yi,-1=Yi,0=0及Y0,t=vt;作为有很多中继节点的传输:vt=F-1(Yt)为如下调幅逆向脉冲传递函数
Figure BSA0000141758170000121
其中Yi-1,t=(2-α)Yi-1,t-1-(1-α)Yi-1,t-2+βg(Yi-1,t-1)+vt-1,t,阀值0<c3<c2<c1,初始值Yi,-1=Yi,0=0,v1,t=c1如果vt≥c1否则
Figure BSA0000141758170000122
vi是初始信号输入和
Figure BSA0000141758170000123
是一个具有一个小的方差的白噪音。称具有方程(4)功能的调幅器为神经元-调幅器,它与NLARI-传感器构成NLARI-传输器。注意到与传统常用的单个激励函数或激励信号不同,本发明采用脉冲时间序列或信号时间序列{vt|t=1,…,n}为激励信号。将模拟信号与实际信号的均方误差的倒数或者动态时间规整(DTW)等作为相似度来评价传感效果,采用信号传输成功率评价传感和传输效果,以区别通常数据传输使用的误差率。
下面结合附图和具体实施例对本发明进一步描述。
本发明提出的NLARI-传感器可按如下具体步骤实现:
第一步选取能够代表被测量信号特征的时间序列作为输入信号vt,在计算机上仿真求出参数仿真域,它在参数域(0,2)作为α,(0,4)作为β和(0,1)作为γ域内,在那里NLARI稳定不动点传递函数Yt=F(vt|α,β,γ)的输出Y=(Y1,…,Yn)′和输入v=(v1,…,vn)′有最大的传感成功率(尖峰脉冲输入)或相似度(非尖峰脉冲输入)(如图2的201所示)。让j=1;
第二步使用压电材料j=1,…,m制作m个传感器样品(如图2的202所示);
第三步输入信号vt到每个传感器样品中,测出响应
Figure BSA0000141758170000124
(图2的203),
对ΔYt=θ1,jΔYt-12,jg(Yt-1)+vt (5)
使用数据
Figure BSA0000141758170000131
作最小2乘法估计
Figure BSA0000141758170000132
在那里ΔYt=Yt-Yt-1,进一步获得
Figure BSA0000141758170000133
作为j=1,…,m(如图2的204所示);
第四步如果参数估计值
Figure BSA0000141758170000134
不在仿真参数域内,该样本不作为NLARI-传感器样本(如图2的205所示),否则计算响应输出
Figure BSA0000141758170000135
与输入v的相似度或者信号传感成功率(j=1,…,m)(如图2的206所示);
第五步选择具有最大输出输入相似度或信号传感成功率的传感器样品作为NLARI-传感器(如图2的207)。
图3显示了临床异常心率HRV序列作为信号输入(301、303、305)和它们经由NLARI-传感器的仿真输出(302、304、306),仿真输出与输入的心率HRV展示了很高的相似度。图4显示了较高频率的正常与临床异常心率HRV序列作为输入信号(401、403、405),类似地,这些心率信号输入与它们的NLARI-传感器的仿真输出(402、404、406)有很高相似度。图3和图4的仿真基于具有α=1.1159、β=0.7625、γ=0.4312的传递函数,包括本发明申请书中所有仿真例都是基于g(x)=-x/(1+x4)的稳定不动点传递函数。
本发明提出两个方案实现尖峰脉冲信号传输,方案1是NLARI-传感器的输出被数据化后通过接口送入计算机(图5的501),由计算机完成信号传输(图5的502)。
方案2是NLARI-传输器方法,结合附图5-8说明它的实现步骤如下:
第一步基于传感功能参数α(k)和β(k),制作信号传感成功率最大的前p个传感器样品,记为NLARI(α(k),β(k))-传感器样本,简称NLARI-传感样品(k=1,…,p)(如图5的503所示),让k=1;
第二步考虑q个相同的NLARI(α(k),β(k))-传感器样本,通过软件或硬件制作具有功能
Figure BSA0000141758170000141
的调幅器样本,其中Yi-1,t由NLARI(α(k),β(k))-传感器样本实现,即
Yi-1,t=(2-α)Yi-1,t-1-(1-α)Yi-1,t-2+βg(Yi-1,t-1)+vi-1,t
其中,vi-1,t和Yi-1,t分别表示经第i-1个NLARI(α(k),β(k))-传感器样本的输入和输出的第t个信号,初始值Yi,-1=Yi,0=0,v1,t=c1如果vt≥c1否则
Figure BSA0000141758170000142
vt是初始信号输入,
Figure BSA0000141758170000143
是一个具有小方差的白噪音(i=1,…,q)。考虑代表被测量信号特征的尖峰脉冲时间序列信号vt和一对由NLARI(α(k),β(k))-传感器样本的输出作为调幅器样本输入的串联系统,计算输出v2,t与输入v1,t的信号传感成功率,通过调节阈值c1,c2,c3选取具有最大信号传感成功率的调幅器样品作为神经元-调幅器样品,它与NLARI(α(k),β(k))-传感器样本合成为NLARI(α(k),β(k))-传输器样本(如图5的504所示);
第三步让k=k+1,重复第二步直到k=p,获得q个相同的NLARI(α(k),β(k))-传输器样品(图5的505所示),让k=1;
第四步对代表被测量信号特征的尖峰脉冲时间序列信号vt,测量由q个传输器样品按前一个的输出作为后一个输入的方式串联起来的最终输出(图5的506所示);
第五步计算串联实验的信号传输成功率,重复第四步获得平均传输成功率(图5的507所示);
第六步如果k<p,让k=k+1,返回第四步(如图5的508所示),否则选取有最大的平均信号传输成功率的传输器样品作为NLARI-传输器(如图5的509所示)。
为区别于数据传输,本发明不采用误码率作传输精确性指标,而导入如下信号传输成功率
Figure BSA0000141758170000151
其中,W0=(w0,1,…,w0,n)′表示初始输入信号,Wq=(wq,1,…,wq,n)′表示经过q次中继节点后的输出信号,w0,t=1如果vt≥c3否则w0,t=0,wq,t=1如果Yq,t≥c3否则wq,t=0,vt是原始输入信号。特别当q=1时,信号传输成功率(7)被称为信号传感成功率。对于只有几个节点的传输,NLARI-传输器可以不含调幅解码器。图6的601显示了输入的尖峰脉冲时间序列信号经NLARI-传输仿真一百次、三千次信号中继的传输成功率,它是随机信号生成器产出一万个时间序列的传输成功率的平均值。在仿真参数域γ∈(0.306,0.412)∪(0.652,0.825),α∈(0.885,0.845)∪(0.755,0.690),β∈(0.682,0.952)∪(1.623,2.160),c1∈(0.0014,0.0016),c2∈(0.0009,0.0011),c3∈(0.0007,0.0009)上,信号传输成功率分别达到了99.00%和99.99%以上的高精度。为说明这些精确传输是如何实现的,图6显示了信号生成器的一次随机实现作为一个信源(如图6的602所示)、接收的初始刺激(如图6的603所示)、第一次对刺激输入的响应(如图6的604所示)和第一次的信号重构(如图6的605所示)、第3000次对刺激输入的响应(如图6的606所示)和第3000次的信号重构(如图6的607所示)。人们可以看到每次响应都大致准确,但存在一些误差,调幅逆向传递函数消除了响应误差。图7显示了被注射(即兴奋性神经元被药物抑制,如图701所示)和未被注射丹皮酚的兴奋性神经元发射(如图703)的尖峰脉冲序列信号分别作为初始输入信号,经NLARI-传输器仿真3000次中继达到无失真传输(如图702和704所示)。
为了说明本发明提高的传感精度的效果,图8显示了NLARI-传输器的精度(如图8的801所示)与一种新型的时删传感器在经过动态测量误差补偿后的精度的一个比较(如图8的802所示)(详细情况请参考“时删传感器动态测量误差补偿”,光学精密工程,第23卷,第4号,2015年4月)。时删传感器的动态测量误差主要包括了具有周期性的系统误差和不具有确定性的随机误差。作者通过傅立叶逼近模型对周期性的系统误差进行分离,使用最小二乘法求解超定方程组的方法推定系统误差补偿参数,以及建立灰色预测模型对随机性成份进行预测和补偿。在相同的实验条件下,对三台不同传感器进行了误差补偿后的最终误差对比。尽管仿真误差要比实体误差小很多,但是NLARI-传输器的精度还是明显高于时删传感器的精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种基于时域编解码的高精度动态信号传感和传输方法,其特征在于,所述基于时域编解码的高精度动态信号传感和传输方法使用非线性自回归整合NLARI过程建模弹性传感器的动态响应,选取满足NLARI模型稳定不动点条件的传感材料使传感器功能相当于一个稳定不动点的传递函数;通过模拟神经元的动作电压全有全无和反向传播及局部电位可叠加的调节功能,导入调幅逆向脉冲传递函数实现信号重构;输出相似度或均方误差以及输出信号传输成功率作为传感和传输评价标准;
稳定不动点的传递函数为:
Figure 966137DEST_PATH_IMAGE002
Figure 166175DEST_PATH_IMAGE004
逆向脉冲传递函数为
Figure 362801DEST_PATH_IMAGE006
,为无调幅时
Figure 507343DEST_PATH_IMAGE008
,为调幅时
Figure 447617DEST_PATH_IMAGE010
,其中,
Figure 236582DEST_PATH_IMAGE012
所述的传递函数
Figure 620421DEST_PATH_IMAGE014
包括:
有界且满足
Figure 2992DEST_PATH_IMAGE016
Figure 809274DEST_PATH_IMAGE018
的恢复力项
Figure 905275DEST_PATH_IMAGE020
,阻力系数
Figure 178124DEST_PATH_IMAGE022
,恢复力系数
Figure 595461DEST_PATH_IMAGE024
以及相对恢复力系数
Figure 815221DEST_PATH_IMAGE026
,其中
Figure DEST_PATH_IMAGE027
Figure 296887DEST_PATH_IMAGE028
Figure 209479DEST_PATH_IMAGE029
所述的传递函数
Figure 114112DEST_PATH_IMAGE002
和逆向传递函数
Figure 262197DEST_PATH_IMAGE006
包括:
输入信号
Figure DEST_PATH_IMAGE031
,响应输出
Figure DEST_PATH_IMAGE033
Figure 739315DEST_PATH_IMAGE033
的一阶滞后值为
Figure DEST_PATH_IMAGE035
Figure 776803DEST_PATH_IMAGE033
的二阶滞后值为
Figure DEST_PATH_IMAGE037
,初始值
Figure DEST_PATH_IMAGE039
Figure DEST_PATH_IMAGE041
;若
Figure 542634DEST_PATH_IMAGE042
Figure DEST_PATH_IMAGE044
,否则
Figure DEST_PATH_IMAGE046
Figure DEST_PATH_IMAGE048
是原始输入信号,
Figure DEST_PATH_IMAGE050
是反映随机误差的白噪音;阈值
Figure DEST_PATH_IMAGE052
Figure 448404DEST_PATH_IMAGE012
表示传输中继节点数及
Figure DEST_PATH_IMAGE054
表示时间或信号数;
所述输出相似度或均方误差作为传感和传输的评价标准中,信号传输成功率
Figure DEST_PATH_IMAGE056
作为尖峰脉冲序列信号传输的评价标准,其中
Figure DEST_PATH_IMAGE058
为经过
Figure DEST_PATH_IMAGE060
次中继节点后的输出信号,
Figure DEST_PATH_IMAGE062
为初始输入信号,若
Figure 468444DEST_PATH_IMAGE063
,则
Figure DEST_PATH_IMAGE065
,否则
Figure DEST_PATH_IMAGE067
;若
Figure DEST_PATH_IMAGE069
,则
Figure DEST_PATH_IMAGE071
,否则
Figure DEST_PATH_IMAGE073
Figure DEST_PATH_IMAGE075
是原始输入信号;将
Figure DEST_PATH_IMAGE077
的信号传输成功率视为信号传感成功率。
2.如权利要求1所述的基于时域编解码的高精度动态信号传感和传输方法,其特征在于,所述传感方法具体的实现步骤包括:
选取能够代表被测量信号特征的时间序列作为输入信号
Figure DEST_PATH_IMAGE079
,在计算机上求出参数仿真域,它位于参数域
Figure DEST_PATH_IMAGE081
作为
Figure DEST_PATH_IMAGE083
Figure DEST_PATH_IMAGE085
作为
Figure DEST_PATH_IMAGE087
,及
Figure DEST_PATH_IMAGE089
作为
Figure DEST_PATH_IMAGE091
的某一区域,在那里NLARI稳定不动点传递函数的输出
Figure DEST_PATH_IMAGE093
和输入
Figure DEST_PATH_IMAGE095
有最大相似度或信号传感成功率;
选取压电材料
Figure DEST_PATH_IMAGE097
制作
Figure DEST_PATH_IMAGE099
个NLARI-传感器样品,输入信号
Figure 448470DEST_PATH_IMAGE079
到所述传感器样品,测出响应输出
Figure DEST_PATH_IMAGE101
,让
Figure DEST_PATH_IMAGE103
,使用
Figure DEST_PATH_IMAGE105
对回归方程
Figure DEST_PATH_IMAGE107
作最小2乘法估计,获得估计值
Figure DEST_PATH_IMAGE109
Figure DEST_PATH_IMAGE111
,以及
Figure DEST_PATH_IMAGE113
Figure DEST_PATH_IMAGE115
Figure DEST_PATH_IMAGE117
如果
Figure DEST_PATH_IMAGE119
在仿真参数域内,计算输出输入相似度或信号传感成功率,选取具有最大相似度或传感成功率的传感器样品为NLARI-传感器样品,否则不作为NLARI-传感器样品;
所述传输方法包括采用计算机传输和NLARI-传感器传输两种,计算机传输是通过NLARI-传感器对接收到的输入信号作出响应输出,所述的响应输出经数据化后通过接口送入计算机完成传输。
3.如权利要求1所述的基于时域编解码的高精度动态信号传感和传输方法,其特征在于,所述NLARI-传感器传输方法具体的实现步骤包括:
步骤一,制作信号传感成功率最大的前
Figure DEST_PATH_IMAGE121
Figure DEST_PATH_IMAGE123
传感器样品,其中,
Figure DEST_PATH_IMAGE125
, 使用软件或硬件制作具有调幅逆向脉冲传递函数功能的调幅器样品,一对由
Figure DEST_PATH_IMAGE127
传感器样品的输出作为调幅器样品的输入的串联系统,计算输出
Figure DEST_PATH_IMAGE129
与输入
Figure DEST_PATH_IMAGE131
的信号传感成功率,通过调节阈值
Figure DEST_PATH_IMAGE133
选取具有最大的信号传感成功率的调幅器样品作为神经元-调幅器样品,与
Figure 563582DEST_PATH_IMAGE127
传感器样品合成为
Figure 459993DEST_PATH_IMAGE060
个相同的
Figure 893728DEST_PATH_IMAGE127
传感器样品;
步骤二,让能够代表被测量信号特征的时间序列信号作为原始输入信号
Figure 365030DEST_PATH_IMAGE134
进入由
Figure 43136DEST_PATH_IMAGE060
个相同
Figure 274397DEST_PATH_IMAGE127
传感器样品按前一个的输出为后一个输入的方式串联的传输系统,测量它由
Figure 347658DEST_PATH_IMAGE134
引起的最终输出;
步骤三,计算每次信号传输成功率,重复步骤二,计算信号平均信号传输成功率;
步骤四,选取具有最大平均信号传输成功率的
Figure 68489DEST_PATH_IMAGE127
传感器样品对
Figure 171574DEST_PATH_IMAGE125
作为NLARI-传感器。
4.如权利要求3所述的基于时域编解码的高精度动态信号传感和传输方法,其特征在于,所述传输方法的步骤二进一步包括:
选取能够代表被测量信号特征的时间序列作为输入信号
Figure DEST_PATH_IMAGE135
,在计算机上仿真求出参数仿真域,它在参数域
Figure 331160DEST_PATH_IMAGE081
作为
Figure 163987DEST_PATH_IMAGE136
Figure 212977DEST_PATH_IMAGE085
作为
Figure DEST_PATH_IMAGE137
Figure 272199DEST_PATH_IMAGE089
作为
Figure 94531DEST_PATH_IMAGE091
域内,在那里NLARI稳定不动点传递函数
Figure DEST_PATH_IMAGE139
的输出
Figure 657230DEST_PATH_IMAGE093
和输入
Figure 408280DEST_PATH_IMAGE079
有最大的传感成功率或相似度;
使用压电材料
Figure DEST_PATH_IMAGE141
制作
Figure 689220DEST_PATH_IMAGE099
个传感器样品;
输入信号
Figure 315242DEST_PATH_IMAGE142
到每个传感器样品中,测出响应
Figure DEST_PATH_IMAGE143
Figure 998027DEST_PATH_IMAGE144
,使用数据
Figure 919978DEST_PATH_IMAGE105
作最小2乘法估计
Figure 157055DEST_PATH_IMAGE146
,在那里
Figure 321189DEST_PATH_IMAGE103
,进一步获得
Figure DEST_PATH_IMAGE148
,其中,
Figure 858481DEST_PATH_IMAGE141
如果参数估计值
Figure DEST_PATH_IMAGE150
不在仿真参数域内,该样品不作为NLARI-传感器样品,否则计算响应输出
Figure 216912DEST_PATH_IMAGE105
与输入
Figure DEST_PATH_IMAGE151
的相似度或者信号传感成功率。
CN201710152542.4A 2017-03-15 2017-03-15 一种基于时域编解码的高精度动态信号传感和传输方法 Active CN107044893B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710152542.4A CN107044893B (zh) 2017-03-15 2017-03-15 一种基于时域编解码的高精度动态信号传感和传输方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710152542.4A CN107044893B (zh) 2017-03-15 2017-03-15 一种基于时域编解码的高精度动态信号传感和传输方法

Publications (2)

Publication Number Publication Date
CN107044893A CN107044893A (zh) 2017-08-15
CN107044893B true CN107044893B (zh) 2021-10-22

Family

ID=59544658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710152542.4A Active CN107044893B (zh) 2017-03-15 2017-03-15 一种基于时域编解码的高精度动态信号传感和传输方法

Country Status (1)

Country Link
CN (1) CN107044893B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107456271A (zh) * 2017-02-20 2017-12-12 葛亚平 一种多功能机电手术系统
CN111803065B (zh) * 2020-06-23 2023-12-26 北方工业大学 一种基于脑电数据的危险交通场景辨识方法及系统
JP2022087776A (ja) * 2020-12-01 2022-06-13 宗路 何 地震波を時間領域に分解して地震前兆現象を捉える地震予測時間領域解析システム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101901209A (zh) * 2010-02-04 2010-12-01 福州大学 基于改进emd和arma模型的结构响应分析方法
CN102520070A (zh) * 2011-12-02 2012-06-27 上海交通大学 基于非线性输出频率响应函数的结构损伤检测方法
WO2013173252A1 (en) * 2012-05-13 2013-11-21 Invention Mine Llc Full duplex wireless transmission with channel phase-based encryption
CN104049228A (zh) * 2014-05-28 2014-09-17 苏州中盛纳米科技有限公司 一种磁传感器的动态模型标定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9997830B2 (en) * 2012-05-13 2018-06-12 Amir Keyvan Khandani Antenna system and method for full duplex wireless transmission with channel phase-based encryption

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101901209A (zh) * 2010-02-04 2010-12-01 福州大学 基于改进emd和arma模型的结构响应分析方法
CN102520070A (zh) * 2011-12-02 2012-06-27 上海交通大学 基于非线性输出频率响应函数的结构损伤检测方法
WO2013173252A1 (en) * 2012-05-13 2013-11-21 Invention Mine Llc Full duplex wireless transmission with channel phase-based encryption
CN104049228A (zh) * 2014-05-28 2014-09-17 苏州中盛纳米科技有限公司 一种磁传感器的动态模型标定方法

Also Published As

Publication number Publication date
CN107044893A (zh) 2017-08-15

Similar Documents

Publication Publication Date Title
Sengupta et al. Consequences of converting graded to action potentials upon neural information coding and energy efficiency
Balavoine et al. Convergence and rate analysis of neural networks for sparse approximation
Shapero et al. Optimal sparse approximation with integrate and fire neurons
Buhry et al. Automated parameter estimation of the Hodgkin-Huxley model using the differential evolution algorithm: application to neuromimetic analog integrated circuits
KR20170041724A (ko) 비동기 펄스 변조를 이용한 인공 뉴런들 및 스파이킹 뉴런들
Schuecker et al. Fundamental activity constraints lead to specific interpretations of the connectome
CN107044893B (zh) 一种基于时域编解码的高精度动态信号传感和传输方法
Van Bussel et al. Inferring synaptic connectivity from spatio-temporal spike patterns
Buhry et al. Global parameter estimation of an Hodgkin–Huxley formalism using membrane voltage recordings: Application to neuro-mimetic analog integrated circuits
Zhao et al. Novel spike based reservoir node design with high performance spike delay loop
Lajoie et al. Structured chaos shapes spike-response noise entropy in balanced neural networks
Martín et al. A hybrid approach to short-term load forecasting aimed at bad data detection in secondary substation monitoring equipment
Sanaullah et al. Exploring spiking neural networks: a comprehensive analysis of mathematical models and applications
Engelken et al. Input correlations impede suppression of chaos and learning in balanced firing-rate networks
Zhang et al. A Radial Basis Function Spike Model for Indirect Learning via Integrate‐and‐Fire Sampling and Reconstruction Techniques
Gunasekaran et al. Convergence of regular spiking and intrinsically bursting Izhikevich neuron models as a function of discretization time with Euler method
CN111210051B (zh) 一种用户用电行为预测方法及系统
Goeva et al. Reconstructing input models via simulation optimization
Orcioni et al. Revealing spectrum features of stochastic neuron spike trains
Krishnan et al. Perfect detection of spikes in the linear sub-threshold dynamics of point neurons
Cessac Retinal processing: Insights from mathematical modelling
Touboul et al. A Markovian event-based framework for stochastic spiking neural networks
Madi et al. Hybrid Cubature Kalman filtering for identifying nonlinear models from sampled recording: Estimation of neuronal dynamics
Végh et al. Towards generalizing the information theory for neural communication
Zhang et al. Deep diffusion Kalman filter combining large-scale neuronal networks simulation with multimodal neuroimaging 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