生物声检测中的双耳多普勒效应仿真研究.docxVIP

生物声检测中的双耳多普勒效应仿真研究 多普勒效应检测 多年来，我们一直在研究如何纠正矩阵效应引起的测量误差，例如克服头皮屑和蝙蝠的差异。这些研究包括多个实验，但仍然是人为推测。McCarty的研究结果说明蝙蝠对运动物体的检测比静止目标更有效,文献的研究也支持生物声纳系统能处理多普勒效应的信号。更为重要的是,生物体用于处理多普勒信号的神经系统是很不稳定的系统,一个不稳定系统是否可以或者如何高灵敏地检测由目标物低速移动所产生的约0.05 Hz的多普勒频偏信号。本研究以多普勒测量为例,分析了生物神经系统中可能存在的对多普勒效应更精细的利用。不仅可能快速处理多普勒信号,而且其灵敏度很高,能精细地测量多普勒效应。 1 材料和方法 1.1 双耳多普勒效应原理 由于运动物体所产生的多普勒效应对于两耳不对称,动物的双侧听觉器官若出现多普勒效应差,即双耳多普勒效应。 双耳多普勒效应可以分为主动和被动两种。主动多普勒效应是蝙蝠发出声音,声音到达目标飞行物再反射回来,由于目标物在飞行中,相对到两耳的速度不同,这就产生多普勒效应差。这种差别,就是主动双耳多普勒效应。如果蝙蝠本身不发出声音。只要目标飞行物发出声音,也构成双耳多普勒效应。称为被动双耳多普勒效应。如果只考虑双耳多普勒效应,则主动和被动是一样的。图1所示为分析蝙蝠的双耳多普勒效应的原理图。 单一声道Doppler效应按式(1)计算 式中,f0为声源B的发射频率;U为声音在介质中传播速度,假定在空气中传播速度=340 m/s;V1为声源B(飞行物)运动速度;α1为声源与接受器连线LB与V1运动方向的夹角;V2为接受器(耳)(R和L)运动速度;α2为声源与接受器L连线LB与V2运动方向的夹角。 按式(1)可求出双耳多普勒效应差 式中,γ是蝙蝠飞行方向角,β1、β2分别表示飞行物位置方向和飞行物运动方向角。 γ=90°,说明蝙蝠飞行方向是沿着头部正前方。此时,m1=m2.式(2)就成为 1)在上述条件下同时有β2=90°时,说明飞行物朝着蝙蝠正方向飞来,此时,Δf=0,表明两耳的多普勒效应差为零; 2)当β2=0时,式(3)成为 此时Δf正负就可以确定飞行物向左还是向右飞行。飞行物运动速度V1越大Δf也越大。m1越小表明目标物距离远,Δf也越小。当M1趋向于零,表明飞行物在无穷远。Δf也趋向零。这表明在很远处两耳多普勒效应差为零。 1.2 飞行速度与主要参数 双耳多普勒效应的利用还取决于其两耳多普勒的频差Δf的大小,以及其神经系统能否达到这样高的分辨率。 自然界中飞行最快的昆虫———蜻蜓的飞行速度约在27.8～41 m/s,而一般的昆虫飞行速度约在3 m/s左右,蝇的飞行速度大约在2.7 m/s,而空气中的声速为340 m/s。 假定昆虫的飞行速度是V1=3 m/s,距离OB=λ=2～5 m,取2 m。声音频率f0=500 Hz。蝙蝠两耳距离30 mm,蝙蝠飞行速度V2=5 m/s 根据以上参数,式(4)为 对于水中海豚也可作同样分析,假定海豚发出的声音频率f0为500 Hz,海豚的最快游动速度约11.111 m/s,假定为V2=1 m/s;食物(小鱼)瞬时游动速度V1=2 m/s,海豚两耳距离150 mm,食物距离OB=3 m,水中声速U=1 000 m/s。则在与上例同样条件情况下式(4)为 从式(5)和式(6)两例看出海豚和蝙蝠两耳多普勒频差约在0.05 Hz,也就是说,神经系统需要能达到这样的分辨能力。 1.3 黄原理在h-h方程中的符号化 假定听觉神经系统所输出的脉冲序列的频率与听神经接收的声音信号的频率成单调关系,也即声音频率越高,输出脉冲序列的频率也越高。因此,神经系统对声音频率的处理,在神经回路内就变成为神经脉冲频率的处理。 用脉冲信号作为输入信号,神经元模型的H-H方程如下 如果H-H方程存在圆映射,就可以用符号动力学方法进行分析。输入等周期脉冲序列、输出脉冲序列可根据符号动力学的方法变为符号序列,具体方法可详见文献。 符号化的规则由式(8)描述,其中的周期格就是图2(a)中ri和ri+1之间的时间段。 按文献定义两符号序列的距离设两序列A和B分别为 其中ai=p1,p2,p3,…,pi,…,pr 定义:两序列A与B的距离 当A=B时,d(A,B)=0 按距离定义,只要符号序列足够长,距离可以小到任意小。符号序列的长度就是脉冲序列的长度。当然在生物神经系统中脉冲序列是不可能无限长的。 1.4 符号序列的计算 以单一频率的声波检测多普勒的频率变化,可以简化为神经元对两个频率非常接近的信号的测量。采用式(7)的H-H方程进行仿真,输入信号为Isig,其形式为等周期的脉冲序列,脉冲幅度为20μA/cm2,宽度为1 ms的方波。计算过程中方程的参数与文献相同。 计算结果按文献和文献中取