第蒜霸期,啦隧璺州
听觉诱发电位检测及分析。
司峻峰1,宁新宝1,薛方2
2.合肥工业大学计算机与信息学院,安徽,合肥。230009)
用接地线的条件下对受试者均成功地检测出了BAEP渡,具有较好的应用价值.关键词:脑干听觉诱发电位,规一化加权平均,低噪声放大器中图分类号:R
318.04
Vo…1.39,No。.1
(1.南京大学电子科学与工程系,近代声学国家重点实验室.江苏,南京,210093
摘要:脑干听觉诱发电位(BAEP)是受试者的听觉系统在具有一定特征的短声刺激作用之后,在其颅顶上检测到的一系列反应波,在电测听和神经电生理上具有独特的临床意义.研究中通过专门设计具有高输入阻抗、低噪声与高共模抑制比的BAEP放大器,建立起一套脑干听觉诱发电位检测系统.在系统的数据分析过程中采用了平均叠加、小波变换等算法.并从临床应用的角度出发提出了规一化系数加权平均算法经实际检测实验表明,本系统在没有屏蔽室和专
DetectingandAnalysisofBrainStemAuditoryEvokedPotentials
Si
Jun-Fen91,Ning船n-Ba01,Xue
Fan92
(1.StateKeyLaboratoryofMordemAcoustics,D籼-imentofElec∞onicScienceandEngine商rig。
NanjingUniversity,Nanjing,210093,China;
2.DepartmentofComputerScienceandEngineering,HefeiAbstract:T}lebrainstemauditoryinthepmcessofwhichis
very
Universityof
Technology,}Ie蹦,230009。China)
nervolls
evokedpotential(BAEP)适abioetectricactionofthecentral
sFstem
stimulationIt
carl
revealtheactivityofneLlr∞groupsalongthe∞”嘲Ⅺ耐i119ser,∞I-ypath,
usefulforclinicaldetection.Butthe
a
anditiscoveredwithisdlfficult
to
largequantityofthe
bymean8of
a
backgr训noise
amplitudeofBAEPis
very
small(about
0.1∥to0.6∥)
ofself-bioelectfidtyandallkindsofm-fifacts.Soit
detectBAEP
u5ualtechnique.
meaning0fdetectingtheevokedpotentialandthedirectionofits
instrumentwehavedesignedandnlade.Wedesctibethe
Inthispaper,firsdy,wediscussthe
development
Then
weintroduceinstrⅢnen
as
a
BAEPdetecting
cornponeiitsofthisfourI)a=rts:aspedaly-d∞ignedamplifierwithhighinputimpedan。e
stimulationcircuits,theanalog
to
and∞“mm
Ⅱ划e蛹ectionratio(CMP,R),sound
system,the
algorithmsofstandardize
digi“converter
and
computercontml
descriptkmofpsYsmeto:softheinstrument.Secondly,inorder_todetectBAEPefficiendy。bo出
weightedaveragingand
wavelet自mn刮kET岫t曲ace
can
employedinthemathematical
calculationofdataacquisition.Wjththesemethods,wereduceaverage
hundredsandhighlyincreasetheaccuracyofresults.Alterthat,we
l目thanfive
discusstheprincipleofstandardkweighted
times“to”or…by
・收稿日期:2002—07—15
第1期甜∞哥n晷In
司俊峰等:听觉诱发电位
’69‘
bythis
andda协
theend,we
introduce咖幻啦脚,all
insta.eesindicatetlnat
ofwhichare
generated
Lrastxunleilt
acq“蓟t|锄姗thods.These
grade
t№instnm蜘t嘶detect
B^EP
succe商ully
witlxmt
Kgh_
earthing∞shield1170a'n.Theyalsosuggestthatithasfairlygoodvaiueintheclinicalpractice・
Keywordb:brainstemauditoryevokedpotential,standardize
weightedaveraging,low-noiseamplifier
在日常生活中,人们不断地接收外界的行为或信息等因素的作用,从而产生一定的反应,为了提高对人类自身高级神经活动本质的认识并进行基础研究,就必须研究这些刺激行为所引起的反应.但自然界所产生的刺激极为复杂,无法进行定量分析和研究,为此。在研究工作中,人为地产生具有一定特征的声、光、电等刺激作用于神经系统,诱发出与神经活动有关的信号,并对这些信号的参数进
行分析,从而从一个侧面反应神经…
系统或大脑的高级心理活动等状态,为临床应用提供有价值的参考数据.脑干听觉诱发电位(BAEP)是受试者在具有一定特征的短声1‘。刺激作用于听觉系统之后,在其颅顶上检测到的一系列反应渡,其中前10ms内的波形具有较稳定的形态,其形状如图1所示:
依次采用罗马字母I、Ⅱ、Ⅲ、Ⅳ、V、Ⅵ、Ⅶ来表示这相继产生的七个波峰,各波峰相应的发生源为:波I,由与耳蜗紧密相连的一段听神经纤维产生;波Ⅱ,由脑干耳蜗核及听神经颅内段的电活动产生;波Ⅲ,与内侧上橄榄核及其场电位有关;波Ⅳ,源于外侧丘系及其核团;波V,源于外侧丘系上方或下方;波Ⅵ、Ⅶ,可能源于外侧丘系与听放射神经.
由以上的分析可以知道,BAEP的各波与听神经和脑干听通路的解剖结构之间有着密切的关系,使其对脑干传导障碍的定位有相当的准确性.例如,通过BAEP波的参数可以对耳聋性质加以判定,功能性耳聋在短声刺激时,脑干反应正常,波V潜伏期(PL)无明显改变.而器质性耳聋、美尼尔病听力障碍和突发性耳聋均使波V的PL值增加.从听神经学的角度出发,当患有听神经瘤、脑于肿瘤、脑血管局灶软化、脑外伤等疾患时,均会使BAEP参数发生变化.其次,BAEP对内、外环境的多种改变具有较强的耐受性,例如某些药物在体内达到一定浓度时,脑电信号(EEG)呈等电位,此时,BAEP波形仍然可见.另外,BAEP可以协助观察相应解剖结构的生理功能,并判断与其相关部位受影响的情况,因而可以用于相应部位的手术监护.这些特征使其在电测听和了解昕神经及其脑千通路或影响该通路的功能性与结构性疚患上具有独特的临床意义.
围1正常听力者的BAEP波形(75dB)
”
啦.1
BAEPwavesofanorllJal
hear吨(75dB)
在临床检测中,为了确保所测潜伏期的准确性,多采用等权平均算法.由于BAEP信号
的幅度非常低,任何相对过高的外部偶发性干扰都会直接影响平均的结果,产生伪迹,于是对测量环境和受试者提出了较高的要求,这些都给BAEP的应用带来了困难.为了减少平均次数,缩短临床检查的时间,提高BAEP检查的成功率和可靠性,近年来,研究者提出了一系
・70・
南京大学学报(自然科学)
第39卷
列算法,如:奇异性检测、正则化方法、维纳滤波、白适应径向神经网络、参数模型、小波变换等等【“,其中的一些算法有效地减少了需要平均的次数,但这些算法多具有后验特性,不利于实时实现,同时算法引入的相位误差会引起时间参数发生较大的变化,在信噪比差时尤为
化加权平均算法,在没有屏蔽室和专用接地线的条件下取得了较好的蛐检出率.具有较
明显.本系统采用专门设计具有高输入阻抗、低噪声与高共模抑制比的BAEP放大器和规一
好的应用价值.
1检测系统硬件构成
系统的硬件由三部分构成,即声刺激器,BAEP放大器和可变带宽滤波器,微机系统与采样电路.系统的硬件结构如图2所示;
图2中,在系统主机的控制下,声刺激器以每秒19次的刺激速率产生刺激短声作用于被测耳,对侧耳采用白噪声进行掩蔽,BJM£P放大器检测到含有BAEP的EEG,滤波器根据被检测的信号的性质对信围2硬件系统框固
号进行预滤波处理,然后对EEG信号进行采样并进№.2
n胛姗0fthe蛔rdm
行实时处理.
数据采集的实时波形可在计算机上实时显示,用以监视所采集数据的有效性,计算机同时能实时显示平均结果的波形,测试时可根据其质量决定是否需要继续测试,其它一些测试的参数也同时进行显示.在数据的提取中采用了规一化系数加权平均算法[2】实时对所采集的数据进行处理.分别介绍如下:
1.I
BAEP放大器和可变带通滤波器为了降低对测试条件的要求,实现对微伏级信号进
行放大而不至于引人过多噪声,必须专门设计具有高输入阻抗与高共模抑制比的低噪声BAEP放大电路旧J,系统中的主放大器采用两级隔离,以提高共模抑制比,在第一、二级采用低噪声器件,以降低系统本底噪声,有效地检测出被测信号,放大器结构如图3所示.可变带通滤波器用于对信号进行预滤波,从而避免在采样过程中产生混叠作用.对于BAEP信号,为了不引起过多相移,高通滤波器(HPF)斜率小于12dR/10倍频,低通滤波器(LPF)斜率小于24dB/10倍频,50Hz陷波器尽可能不使用.1.2声刺激器
由刺激短声产生电
路、伪白噪声产生电路、控制电路构成.通过对设定的地址写人控制字可以实现短声极性选择,左右耳选择,刺激声与自噪声强度选择.伪白噪声通过对伪随机序列进行低通滤波后得到,伪随机田3
BAEP放大舞框田
序列产生电路由两片18位串行移位寄№・3n哪dl删ofBA耻mplin日
存器(71;)4006串联构成,将第13位输
出和第31位输出进行异或后反馈至输人端形成一个长度为32位的移位寄存器,产生的伪随机序列长度为2147483647,当时钟频率为300kI-'Iz时,序列重复时间约为2
h【….
1.3微机系统与采样电路微机采用Pc,通过对地址总线译码产生系统的控制信号,用于
第1期
司俊峰等:听觉诱发电位
。71‘
对放大器的增益以及滤波器的滤波特性实现控制.A/D采样芯片采用MAXl91,在定时器芯片8254的控制下当Pc或外部信号触发后进行帧数据采集.1.4系统性能指标
1,4.1
BAEP放大器与滤波器
5X103—5Xi06倍大于40MO共模抑制比1-200Hz短路噪声
放大倍数输人阻抗频率下限1.4.2数字系统A/D分辨率
大于120扭频率上限20
Hz~8kI-k
小于2pVp—P(1~8k}王z)
12Bit
平均次数分析时间
1~8192
刺激声强
0~120dB
t-lz
单帧采样点数256~20481.4.3系统显示13AEP参数
5~3000
Ins
刺激频率0.2~100
当前次数CN,最高峰值幅度Umax,I、Ⅱ、Ⅲ、Ⅳ、V、Ⅵ、Ⅶ各峰潜伏期值与I一Ⅲ、Ⅲ一Ⅳ、I—V峰间期值.2
BAEP处理算法
2.1规一化系数加权平均用于BAEP提取在BAEP检测中,EEG电极在头皮提取的信号由EEG和BAEP相叠加形成,由于BAEP与EEG基本处于同~频域且EEG幅度为BAEP数十倍,难以通过简单的滤波进行提取.普通提取方法是假定EEG为平稳随机信号
而BAEP具有稳定的锁时特性,以刺激声触发时刻为基准,对触发后采集到的帧数据进行
1000至2000次等权平均叠加得出BAEP波,缺点是平均次数多,时间较长,采集中如遇坏的帧数据将直接影响检测的结果.有些文献中提到参数模型、选择平均、自适应滤波等算法,这些方法往往需要先提取一定参数,难以实时实现,对所测BAEP的幅度也产生一定的影响。同时在信噪比较低时会导致峰问潜伏期(IPL)值产生较明显的变化.本研究中采用规一化系数加权平均提取算法Ⅱ.6】.有效地降低了测试的次数,对幅度也无明显影响.
算法假定每次刺激均能诱发出BAEP,由BAEP的产生机理可以认为其具有相对稳定性[…,即对于反复多次的同样刺激所诱发出的BAEP保持有同样的潜伏期和极性,其波形、波幅也基本一致;同时,背景噪声也与刺激无任何相关成分.在BAEP的每帧数据中,EEG的幅度为20--100btVp-p,测试系统的本底噪声幅度为2FtVp-p,BAEP的幅度小于1
txVP-
p,在检测过程中,数据由于受试者的影响难以长期保持平稳,而信噪比极差的数据会直接影响叠加结果,因而算法对每帧数据在叠加前进行阈值检测,放弃过闭值点超过一定数量的帧数据.在正式叠加前做预叠加,求出32至64帧数据的方差并将这32至64个方差的均值作
为归一化值.预叠加数据中第i帧数据的方差一:为:
d:=E{[zf(n)一E(z;(”))]2}
(1)
M帧预叠加数据方差的均值02删为:
a2姗=击姜毛
(2)
正式进行叠加时,先求出可以进行叠加的第k帧数据的方差畦及该帧数据的加权值
.72.
南京大学学报(自然科学)
第39卷
嘞,再将该帧数据进行加权叠加.每帧数据的方差以。为:
d:.=E{[zl(n)一E(x^(”))]2
该帧数据的加权值w^为
(3)
wk=√口2删厅:。
N帧实测数据叠加结果Xa。(n)为:
.
(4)
N
.27av(m)=而1∑巩墨(”)
(5)
2.2小波变换用于BAEP去噪小波变换是一种时频分析方法[副,能在多尺度下对信号进行分解与重构,有效的对信号的轮廓和细节进行分析,算法中利用小波变换的时频特性分析并降低了BAEP中干扰脉冲的影响.
对于信号s(t)∈L2(R),正交二进小波变换及逆变换定义如下:
c¨=2-J。l。s(t)Zv(2一Jt—k)dt=I。s(£)丐I(£)d£
s(t)=印∑∑o.^q。^(£)=∑sj(t)
(6)(7)(8)
J…t—
竺
J2一”
其中
_(t)=co∑c/.j哆.^(t),J∈z
表示信号s(£)在尺度a=2j下的分量,cn是一个与信号无关的常数.当采用二进抽取时,二进小波变换可表示下述滤波器形式【9J:
弓(k)=L;h(2k—n)≈一1(n)
(9)
4(^)=∑g(2k一”)弓一1(n)(10)
刁一l(^)=∑h(k一2口)弓(”)+∑g(^一2n)弓(n)
(11)
其中zo(”)表示离散信号,而(n)表示信号在尺度4=∥下的轮廓,di(”)表示信号在尺度a=∥下的细节,h(n)和g(n)表示正交离散滤波器的系数.在对BAEP信号进行去噪处理时,我们选择分解长度为17、重构长度为13的BiorithgorM紧支小波,该小波具有广义线性相位,在信号的处理中可以使相位失真降低到最/btlo].
图4为叠加后得到的BAEP信号波形,其中第1nb内的脉冲产生于耳机的电磁作用。
其幅度因环境而异,该脉冲幅度较大时,对BAEPI波的相位和参数的读取形成了干扰.根据Mallat算法【11],对BAEP信号进行6级分解,信号空间口。表示为:w1∈伊W20w30…o
w60v6,相应的子带信号为:奸(”),《(n),…《(n),瑶(n).图6为各尺度下的子带信号,从图中可以清楚的看出干扰脉冲主要出现在子带信号sf(一),《(n),d(n)中,为了消除该脉冲带来的影响,将予带信号s{(”),《(”),舅(”)中脉冲出现时刻的系数暨为零,然后对
信号进行萤构.图5是用上述方法消除脉冲干扰后的BAEP波形,与原始信号相比,去噪效果是明显的.
—————————————————————————————————————————————————————————一
第1期
司俊峰等:听觉诱发电位
‘73‘
SF/1————~
s5。/\/—————一s4‘,r—、4,。———一
ss2一卜———————一3d
p———一——
田4原始BAEP波形
圈5消除干扰脉冲千扰后的
sI・r_———一——。
lk.4
TheoriginalBAEP
waves
BAEP波形围6
BAEP子带倍号
n昏5
BAEP
wam
withoutnoise
飑.6躏F日batlevelsIt06
3实验结果和讨论
我们采用样机在一般室内对数十位受试者进行了BAEP检测,声刺激强度设定为75血,
刺激声为短声,刺激频率为19地,放大器增益为50万倍,带宽为50~3受试者左、右耳的检测结果的波形,表1列出了部分受试者BAEP㈣参数.
0001-lz,叠加次数为
512次,分别对左、右耳进行检测,均成功地检测到重复性较好的清晰的BAEP波形.图7为某
≥
b影∥少
b膨
℃雪影
彩
√嘭\,
◇
过
圈7某受试者左耳(左)和右耳(右)BAEP波形
Fig.7
BAEPwaves
efleftfar(1eft)andrighte口(rIght)
考察BAEP是否正常通常从波形的振幅(AMP)、绝对潜伏期(PL)、峰间潜伏期(IPL)、双耳潜伏期差(ILD)和波形的可重复性几个方面进行u…,正常受试者左右耳元显著差翼,但个体PL和IPL的ILD超过0.4ms时则需考虑是否为异常,而PL与IPL超过正常均值3SD者则为异常,V/I幅度比值或Ⅳ~V复合体/I幅度比值<0.5应考虑为异常.若波I—V的ILD延长或波I—V的IPL延长均提示耳蜗后病变,需要进一步进行定位检查.另外由于个体之间的差异,在分析测试结果时,应对受试者两耳的参数加以比较.波Ⅵ、Ⅶ的出现不规则,对其参数通常不作要求.表217J列出了一般情况下BAEP波各主要参数的正常值
范围.具体病例及其对御参数的影响请参考相关资料.
诱发电位(EP)是继脑电图和肌电图之后l艋床电生理的第三大进展,其研究始于五十年代,但真正发展还是从20世纪70年代开始的,BAEP由Robinson等人于1975年首次报道.这其中主要的原因是由于电子技术的限制,本研究综合了前人的研究成果,提出和比较了一些具有应用价值的分析和检测方法,实验表明该检测系统具有较好的BAEP检出率.随着电子技术和信号处理技术的飞速发展,检测BAEP的技术会不断改进,从而在临床上得到普遍应用.
.74.
南京大学学报(自然科学)
第39
袭1部分受试者的BAEP实测值
者1(左)l(右)2(左)2(右)3(左)3(右)4(左)4(右)5(左)5(右)6(左)6(右)7(左)7(右)8(左)8(右)9(左)9(右)lo(左)lo(右)ll(左)li(右)12(左)12(右)
I
1.64
Ⅱ
3.923.883.563.483.84
m
5.445.845.48
i一Ⅲ
2.282,28
Ⅲ一V
1.52
i—V
3.84243.883.684364.24.34.264.o4.123.96
1.印
1.601.641.41.56l461
54
1.96
1.921.84
1.96
1.842.442.082.262.341.842042.322.362062.121.881.922.422.302.242.122.442.32
5.32
5.765.765.76
1.92
2.122.04
3.“
3.723.883.56
5.80
5.725.725,40
1.922.16
2081.64
1.72l601441.481.741.841.681.791501621.521.521.441.361.561.44
3.“
3.763.843.803.963.563.683.923.923.753.643.88
5.525.84
5.84
1.68
2.041.88
4.04
4.14.o3.643:723
5.32
5.485.445.645.56
1.76
1.81.521.721.801.721.441.61.681.84
94
4.024.043.84
5.365.325.285.76
5.88
3.883.92
4.24.44
3.68
4.084.04
2.52
2.6
第1期司俊峰等:听觉诱发电位
・75‘
[1]Zhu
C,Hu
G.Advanceinsingle-trialEP/ERPestimation.ForeignMedicine
onBiomedical
EI】gin∞j119・
2000,23(4):211--216(朱常芳,胡广书.诱发电位快速提取算法的新进展.国外医学(生物医学工程分册),2000,23(4):211~216).[2]Ning
x,Shen
Z,ShenD,eta1.Thedataacquisitionandanalysis(NaturalSciences),1991,27(2):273
284(宁新宝,沈振字,沈德才等.高频心电图的数据采集与分析南京大学学报(自然科学),1991,27(2):273--284).
[3]NingX,Jinz,ShaoJ.Reducingnoisebyintegratedopemtiomlamplifiersinparallel.JournalofDataAe-
quisition&Processing,1990,5(2):59~65.(宁新宝,金之仑,邵军.用多个集成运放并联降低心电放大器噪声.数据采集与处理,1990,5(2);59--65).[4]Lin
K.Principleandacquisition
lishing
House(丹嗍),1978,41.(林可祥.伪随机码的原理与应用北京:人民邮电出版社,1978,
E
offakerandoraccde.Belj'ing:People’sPostandTelecommunicationsPub—
41)
[5]Carlos[63
D,M。haIrl—S
Mobinl.WeightedR.Welghted
averaging
ofevokedpotentials.IEEETrarⅡ-,action
on
BicmedlcalEn西neexing,1992,BME-39:338~345.
HokeM,R卿B,Wickesberg
try.Electroencephalography
av刚ngtheory
and
application
to
dectrieresponse
audiome.
Cllniq【Neurophysiology,1984,57:84--89
[7]PanY.E∞kedpoteatialshmodicatpractice.BelJing:People’sMedicalPublishingHouse(B.vMP)。1998.
(潘映福临床诱发电位学北京:人民卫生出版社。1988).[8]
Ning
D,YanW.TheanalySis
electrw_ardiomignal.JournalofNankingUni雌ity(NalurdSci嘲),2002,38(1):7~13.(宁瓿宝.张
道斌.阎文泰小波变换用于心电信号高频成份的分析.南京大学学报,2002,38(1):7~13).
Steven
X,Zhang
ofwavelettransformationfor
h讪frequency
compolxentsfrcrn
[9JJS,AkramA,MichaelU,eta1.FastwavelettransformationofEEG.Electmencephalography
Neumphysiolcgy,1994,91:442~455.
andC1inical
[10]ZhaoS,XongX-Wavelettramfom】ationandAnal蝴s.BeiJing:PublishingHouseofFAecttlmi,:sIndustry
(PHEl),1996,64.(ig松年,熊小芸.子波变换与子渡分析.北京:电子工业出版社,1996,64),
[11]ThakorNV,Multi[12]Jorge
resolutionwavelet
analyslsofevokedpotentials.IEEETransactiononBiomedicalEn.
sineering,1993,BME-40:l
IA
085~1089.
EvokedpotentialsResearch.IEEEⅡ19irleeringinMedicineandBiology,1992,67~68.
听觉诱发电位检测及分析
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
司峻峰, 宁新宝, 薛方
司峻峰,宁新宝(南京大学电子科学与工程系,近代声学国家重点实验室,江苏,南京,210093), 薛方(合肥工业大学计算机与信息学院,安徽,合肥,230009)南京大学学报(自然科学版)
JOURNAL OF NANJING UNIVERSITY(NATURAL SCIENCE)2003,39(1)4次
参考文献(12条)
1. Jorge I A Evoked potentials Research[外文期刊] 1992
2. Thakor N V Multi resolution wavelet analysis of evoked potentials[外文期刊] 1993(11)3. 赵松年;熊小芸 子波变换与子波分析 1996
4. 朱常芳;胡广书 诱发电位快速提取算法的新进展 2000(04)
5. STEVEN J S;Akram A;Michael U Fast wavelet transformation of EEG 1994
6. 宁新宝;张道斌;阎文泰 小波变换用于心电信号高频成份的分析[期刊论文]-南京大学学报(自然科学版)2002(01)
7. 潘映福 临床诱发电位学 1988
8. Hoke M;Ross B;Wickesberg R Weighted averaging theory and application to electric responseaudiometry 1984
9. Carlos E D;Mohammad S Mobinl Weighted averaging of evoked potentials[外文期刊] 199210. 林可祥 伪随机码的原理与应用 1978
11. 宁新宝;金之仑;邵军 用多个集成运放并联降低心电放大器噪声 1990(02)12. 宁新宝;沈振宇;沈德才 高频心电图的数据采集与分析 1991(02)
引证文献(5条)
1. 高瑞静. 孙迎. 李越囡 脑干听觉诱发电位波形的自动检测及新标识[期刊论文]-中国组织工程研究与临床康复2010(52)
2. 徐灵活. 殷明. 章建程. 唐志文. 李中付. 袁海霞. 罗志辉 噪声听觉损伤易感动物筛选的噪声暴露方案探讨[期刊论文]-中华航海医学与高气压医学杂志 2007(3)
3. 傅娟 脑电生物反馈系统及其研究[学位论文]硕士 20044. 司峻峰 事件相关脑电活动的检测与分类[学位论文]博士 2003
5. 高瑞静. 孙迎. 李越囡 脑干听觉诱发电位波形的自动检测及新标识[期刊论文]-中国组织工程研究与临床康复2010(52)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_njdxxb200301009.aspx
第蒜霸期,啦隧璺州
听觉诱发电位检测及分析。
司峻峰1,宁新宝1,薛方2
2.合肥工业大学计算机与信息学院,安徽,合肥。230009)
用接地线的条件下对受试者均成功地检测出了BAEP渡,具有较好的应用价值.关键词:脑干听觉诱发电位,规一化加权平均,低噪声放大器中图分类号:R
318.04
Vo…1.39,No。.1
(1.南京大学电子科学与工程系,近代声学国家重点实验室.江苏,南京,210093
摘要:脑干听觉诱发电位(BAEP)是受试者的听觉系统在具有一定特征的短声刺激作用之后,在其颅顶上检测到的一系列反应波,在电测听和神经电生理上具有独特的临床意义.研究中通过专门设计具有高输入阻抗、低噪声与高共模抑制比的BAEP放大器,建立起一套脑干听觉诱发电位检测系统.在系统的数据分析过程中采用了平均叠加、小波变换等算法.并从临床应用的角度出发提出了规一化系数加权平均算法经实际检测实验表明,本系统在没有屏蔽室和专
DetectingandAnalysisofBrainStemAuditoryEvokedPotentials
Si
Jun-Fen91,Ning船n-Ba01,Xue
Fan92
(1.StateKeyLaboratoryofMordemAcoustics,D籼-imentofElec∞onicScienceandEngine商rig。
NanjingUniversity,Nanjing,210093,China;
2.DepartmentofComputerScienceandEngineering,HefeiAbstract:T}lebrainstemauditoryinthepmcessofwhichis
very
Universityof
Technology,}Ie蹦,230009。China)
nervolls
evokedpotential(BAEP)适abioetectricactionofthecentral
sFstem
stimulationIt
carl
revealtheactivityofneLlr∞groupsalongthe∞”嘲Ⅺ耐i119ser,∞I-ypath,
usefulforclinicaldetection.Butthe
a
anditiscoveredwithisdlfficult
to
largequantityofthe
bymean8of
a
backgr训noise
amplitudeofBAEPis
very
small(about
0.1∥to0.6∥)
ofself-bioelectfidtyandallkindsofm-fifacts.Soit
detectBAEP
u5ualtechnique.
meaning0fdetectingtheevokedpotentialandthedirectionofits
instrumentwehavedesignedandnlade.Wedesctibethe
Inthispaper,firsdy,wediscussthe
development
Then
weintroduceinstrⅢnen
as
a
BAEPdetecting
cornponeiitsofthisfourI)a=rts:aspedaly-d∞ignedamplifierwithhighinputimpedan。e
stimulationcircuits,theanalog
to
and∞“mm
Ⅱ划e蛹ectionratio(CMP,R),sound
system,the
algorithmsofstandardize
digi“converter
and
computercontml
descriptkmofpsYsmeto:softheinstrument.Secondly,inorder_todetectBAEPefficiendy。bo出
weightedaveragingand
wavelet自mn刮kET岫t曲ace
can
employedinthemathematical
calculationofdataacquisition.Wjththesemethods,wereduceaverage
hundredsandhighlyincreasetheaccuracyofresults.Alterthat,we
l目thanfive
discusstheprincipleofstandardkweighted
times“to”or…by
・收稿日期:2002—07—15
第1期甜∞哥n晷In
司俊峰等:听觉诱发电位
’69‘
bythis
andda协
theend,we
introduce咖幻啦脚,all
insta.eesindicatetlnat
ofwhichare
generated
Lrastxunleilt
acq“蓟t|锄姗thods.These
grade
t№instnm蜘t嘶detect
B^EP
succe商ully
witlxmt
Kgh_
earthing∞shield1170a'n.Theyalsosuggestthatithasfairlygoodvaiueintheclinicalpractice・
Keywordb:brainstemauditoryevokedpotential,standardize
weightedaveraging,low-noiseamplifier
在日常生活中,人们不断地接收外界的行为或信息等因素的作用,从而产生一定的反应,为了提高对人类自身高级神经活动本质的认识并进行基础研究,就必须研究这些刺激行为所引起的反应.但自然界所产生的刺激极为复杂,无法进行定量分析和研究,为此。在研究工作中,人为地产生具有一定特征的声、光、电等刺激作用于神经系统,诱发出与神经活动有关的信号,并对这些信号的参数进
行分析,从而从一个侧面反应神经…
系统或大脑的高级心理活动等状态,为临床应用提供有价值的参考数据.脑干听觉诱发电位(BAEP)是受试者在具有一定特征的短声1‘。刺激作用于听觉系统之后,在其颅顶上检测到的一系列反应渡,其中前10ms内的波形具有较稳定的形态,其形状如图1所示:
依次采用罗马字母I、Ⅱ、Ⅲ、Ⅳ、V、Ⅵ、Ⅶ来表示这相继产生的七个波峰,各波峰相应的发生源为:波I,由与耳蜗紧密相连的一段听神经纤维产生;波Ⅱ,由脑干耳蜗核及听神经颅内段的电活动产生;波Ⅲ,与内侧上橄榄核及其场电位有关;波Ⅳ,源于外侧丘系及其核团;波V,源于外侧丘系上方或下方;波Ⅵ、Ⅶ,可能源于外侧丘系与听放射神经.
由以上的分析可以知道,BAEP的各波与听神经和脑干听通路的解剖结构之间有着密切的关系,使其对脑干传导障碍的定位有相当的准确性.例如,通过BAEP波的参数可以对耳聋性质加以判定,功能性耳聋在短声刺激时,脑干反应正常,波V潜伏期(PL)无明显改变.而器质性耳聋、美尼尔病听力障碍和突发性耳聋均使波V的PL值增加.从听神经学的角度出发,当患有听神经瘤、脑于肿瘤、脑血管局灶软化、脑外伤等疾患时,均会使BAEP参数发生变化.其次,BAEP对内、外环境的多种改变具有较强的耐受性,例如某些药物在体内达到一定浓度时,脑电信号(EEG)呈等电位,此时,BAEP波形仍然可见.另外,BAEP可以协助观察相应解剖结构的生理功能,并判断与其相关部位受影响的情况,因而可以用于相应部位的手术监护.这些特征使其在电测听和了解昕神经及其脑千通路或影响该通路的功能性与结构性疚患上具有独特的临床意义.
围1正常听力者的BAEP波形(75dB)
”
啦.1
BAEPwavesofanorllJal
hear吨(75dB)
在临床检测中,为了确保所测潜伏期的准确性,多采用等权平均算法.由于BAEP信号
的幅度非常低,任何相对过高的外部偶发性干扰都会直接影响平均的结果,产生伪迹,于是对测量环境和受试者提出了较高的要求,这些都给BAEP的应用带来了困难.为了减少平均次数,缩短临床检查的时间,提高BAEP检查的成功率和可靠性,近年来,研究者提出了一系
・70・
南京大学学报(自然科学)
第39卷
列算法,如:奇异性检测、正则化方法、维纳滤波、白适应径向神经网络、参数模型、小波变换等等【“,其中的一些算法有效地减少了需要平均的次数,但这些算法多具有后验特性,不利于实时实现,同时算法引入的相位误差会引起时间参数发生较大的变化,在信噪比差时尤为
化加权平均算法,在没有屏蔽室和专用接地线的条件下取得了较好的蛐检出率.具有较
明显.本系统采用专门设计具有高输入阻抗、低噪声与高共模抑制比的BAEP放大器和规一
好的应用价值.
1检测系统硬件构成
系统的硬件由三部分构成,即声刺激器,BAEP放大器和可变带宽滤波器,微机系统与采样电路.系统的硬件结构如图2所示;
图2中,在系统主机的控制下,声刺激器以每秒19次的刺激速率产生刺激短声作用于被测耳,对侧耳采用白噪声进行掩蔽,BJM£P放大器检测到含有BAEP的EEG,滤波器根据被检测的信号的性质对信围2硬件系统框固
号进行预滤波处理,然后对EEG信号进行采样并进№.2
n胛姗0fthe蛔rdm
行实时处理.
数据采集的实时波形可在计算机上实时显示,用以监视所采集数据的有效性,计算机同时能实时显示平均结果的波形,测试时可根据其质量决定是否需要继续测试,其它一些测试的参数也同时进行显示.在数据的提取中采用了规一化系数加权平均算法[2】实时对所采集的数据进行处理.分别介绍如下:
1.I
BAEP放大器和可变带通滤波器为了降低对测试条件的要求,实现对微伏级信号进
行放大而不至于引人过多噪声,必须专门设计具有高输入阻抗与高共模抑制比的低噪声BAEP放大电路旧J,系统中的主放大器采用两级隔离,以提高共模抑制比,在第一、二级采用低噪声器件,以降低系统本底噪声,有效地检测出被测信号,放大器结构如图3所示.可变带通滤波器用于对信号进行预滤波,从而避免在采样过程中产生混叠作用.对于BAEP信号,为了不引起过多相移,高通滤波器(HPF)斜率小于12dR/10倍频,低通滤波器(LPF)斜率小于24dB/10倍频,50Hz陷波器尽可能不使用.1.2声刺激器
由刺激短声产生电
路、伪白噪声产生电路、控制电路构成.通过对设定的地址写人控制字可以实现短声极性选择,左右耳选择,刺激声与自噪声强度选择.伪白噪声通过对伪随机序列进行低通滤波后得到,伪随机田3
BAEP放大舞框田
序列产生电路由两片18位串行移位寄№・3n哪dl删ofBA耻mplin日
存器(71;)4006串联构成,将第13位输
出和第31位输出进行异或后反馈至输人端形成一个长度为32位的移位寄存器,产生的伪随机序列长度为2147483647,当时钟频率为300kI-'Iz时,序列重复时间约为2
h【….
1.3微机系统与采样电路微机采用Pc,通过对地址总线译码产生系统的控制信号,用于
第1期
司俊峰等:听觉诱发电位
。71‘
对放大器的增益以及滤波器的滤波特性实现控制.A/D采样芯片采用MAXl91,在定时器芯片8254的控制下当Pc或外部信号触发后进行帧数据采集.1.4系统性能指标
1,4.1
BAEP放大器与滤波器
5X103—5Xi06倍大于40MO共模抑制比1-200Hz短路噪声
放大倍数输人阻抗频率下限1.4.2数字系统A/D分辨率
大于120扭频率上限20
Hz~8kI-k
小于2pVp—P(1~8k}王z)
12Bit
平均次数分析时间
1~8192
刺激声强
0~120dB
t-lz
单帧采样点数256~20481.4.3系统显示13AEP参数
5~3000
Ins
刺激频率0.2~100
当前次数CN,最高峰值幅度Umax,I、Ⅱ、Ⅲ、Ⅳ、V、Ⅵ、Ⅶ各峰潜伏期值与I一Ⅲ、Ⅲ一Ⅳ、I—V峰间期值.2
BAEP处理算法
2.1规一化系数加权平均用于BAEP提取在BAEP检测中,EEG电极在头皮提取的信号由EEG和BAEP相叠加形成,由于BAEP与EEG基本处于同~频域且EEG幅度为BAEP数十倍,难以通过简单的滤波进行提取.普通提取方法是假定EEG为平稳随机信号
而BAEP具有稳定的锁时特性,以刺激声触发时刻为基准,对触发后采集到的帧数据进行
1000至2000次等权平均叠加得出BAEP波,缺点是平均次数多,时间较长,采集中如遇坏的帧数据将直接影响检测的结果.有些文献中提到参数模型、选择平均、自适应滤波等算法,这些方法往往需要先提取一定参数,难以实时实现,对所测BAEP的幅度也产生一定的影响。同时在信噪比较低时会导致峰问潜伏期(IPL)值产生较明显的变化.本研究中采用规一化系数加权平均提取算法Ⅱ.6】.有效地降低了测试的次数,对幅度也无明显影响.
算法假定每次刺激均能诱发出BAEP,由BAEP的产生机理可以认为其具有相对稳定性[…,即对于反复多次的同样刺激所诱发出的BAEP保持有同样的潜伏期和极性,其波形、波幅也基本一致;同时,背景噪声也与刺激无任何相关成分.在BAEP的每帧数据中,EEG的幅度为20--100btVp-p,测试系统的本底噪声幅度为2FtVp-p,BAEP的幅度小于1
txVP-
p,在检测过程中,数据由于受试者的影响难以长期保持平稳,而信噪比极差的数据会直接影响叠加结果,因而算法对每帧数据在叠加前进行阈值检测,放弃过闭值点超过一定数量的帧数据.在正式叠加前做预叠加,求出32至64帧数据的方差并将这32至64个方差的均值作
为归一化值.预叠加数据中第i帧数据的方差一:为:
d:=E{[zf(n)一E(z;(”))]2}
(1)
M帧预叠加数据方差的均值02删为:
a2姗=击姜毛
(2)
正式进行叠加时,先求出可以进行叠加的第k帧数据的方差畦及该帧数据的加权值
.72.
南京大学学报(自然科学)
第39卷
嘞,再将该帧数据进行加权叠加.每帧数据的方差以。为:
d:.=E{[zl(n)一E(x^(”))]2
该帧数据的加权值w^为
(3)
wk=√口2删厅:。
N帧实测数据叠加结果Xa。(n)为:
.
(4)
N
.27av(m)=而1∑巩墨(”)
(5)
2.2小波变换用于BAEP去噪小波变换是一种时频分析方法[副,能在多尺度下对信号进行分解与重构,有效的对信号的轮廓和细节进行分析,算法中利用小波变换的时频特性分析并降低了BAEP中干扰脉冲的影响.
对于信号s(t)∈L2(R),正交二进小波变换及逆变换定义如下:
c¨=2-J。l。s(t)Zv(2一Jt—k)dt=I。s(£)丐I(£)d£
s(t)=印∑∑o.^q。^(£)=∑sj(t)
(6)(7)(8)
J…t—
竺
J2一”
其中
_(t)=co∑c/.j哆.^(t),J∈z
表示信号s(£)在尺度a=2j下的分量,cn是一个与信号无关的常数.当采用二进抽取时,二进小波变换可表示下述滤波器形式【9J:
弓(k)=L;h(2k—n)≈一1(n)
(9)
4(^)=∑g(2k一”)弓一1(n)(10)
刁一l(^)=∑h(k一2口)弓(”)+∑g(^一2n)弓(n)
(11)
其中zo(”)表示离散信号,而(n)表示信号在尺度4=∥下的轮廓,di(”)表示信号在尺度a=∥下的细节,h(n)和g(n)表示正交离散滤波器的系数.在对BAEP信号进行去噪处理时,我们选择分解长度为17、重构长度为13的BiorithgorM紧支小波,该小波具有广义线性相位,在信号的处理中可以使相位失真降低到最/btlo].
图4为叠加后得到的BAEP信号波形,其中第1nb内的脉冲产生于耳机的电磁作用。
其幅度因环境而异,该脉冲幅度较大时,对BAEPI波的相位和参数的读取形成了干扰.根据Mallat算法【11],对BAEP信号进行6级分解,信号空间口。表示为:w1∈伊W20w30…o
w60v6,相应的子带信号为:奸(”),《(n),…《(n),瑶(n).图6为各尺度下的子带信号,从图中可以清楚的看出干扰脉冲主要出现在子带信号sf(一),《(n),d(n)中,为了消除该脉冲带来的影响,将予带信号s{(”),《(”),舅(”)中脉冲出现时刻的系数暨为零,然后对
信号进行萤构.图5是用上述方法消除脉冲干扰后的BAEP波形,与原始信号相比,去噪效果是明显的.
—————————————————————————————————————————————————————————一
第1期
司俊峰等:听觉诱发电位
‘73‘
SF/1————~
s5。/\/—————一s4‘,r—、4,。———一
ss2一卜———————一3d
p———一——
田4原始BAEP波形
圈5消除干扰脉冲千扰后的
sI・r_———一——。
lk.4
TheoriginalBAEP
waves
BAEP波形围6
BAEP子带倍号
n昏5
BAEP
wam
withoutnoise
飑.6躏F日batlevelsIt06
3实验结果和讨论
我们采用样机在一般室内对数十位受试者进行了BAEP检测,声刺激强度设定为75血,
刺激声为短声,刺激频率为19地,放大器增益为50万倍,带宽为50~3受试者左、右耳的检测结果的波形,表1列出了部分受试者BAEP㈣参数.
0001-lz,叠加次数为
512次,分别对左、右耳进行检测,均成功地检测到重复性较好的清晰的BAEP波形.图7为某
≥
b影∥少
b膨
℃雪影
彩
√嘭\,
◇
过
圈7某受试者左耳(左)和右耳(右)BAEP波形
Fig.7
BAEPwaves
efleftfar(1eft)andrighte口(rIght)
考察BAEP是否正常通常从波形的振幅(AMP)、绝对潜伏期(PL)、峰间潜伏期(IPL)、双耳潜伏期差(ILD)和波形的可重复性几个方面进行u…,正常受试者左右耳元显著差翼,但个体PL和IPL的ILD超过0.4ms时则需考虑是否为异常,而PL与IPL超过正常均值3SD者则为异常,V/I幅度比值或Ⅳ~V复合体/I幅度比值<0.5应考虑为异常.若波I—V的ILD延长或波I—V的IPL延长均提示耳蜗后病变,需要进一步进行定位检查.另外由于个体之间的差异,在分析测试结果时,应对受试者两耳的参数加以比较.波Ⅵ、Ⅶ的出现不规则,对其参数通常不作要求.表217J列出了一般情况下BAEP波各主要参数的正常值
范围.具体病例及其对御参数的影响请参考相关资料.
诱发电位(EP)是继脑电图和肌电图之后l艋床电生理的第三大进展,其研究始于五十年代,但真正发展还是从20世纪70年代开始的,BAEP由Robinson等人于1975年首次报道.这其中主要的原因是由于电子技术的限制,本研究综合了前人的研究成果,提出和比较了一些具有应用价值的分析和检测方法,实验表明该检测系统具有较好的BAEP检出率.随着电子技术和信号处理技术的飞速发展,检测BAEP的技术会不断改进,从而在临床上得到普遍应用.
.74.
南京大学学报(自然科学)
第39
袭1部分受试者的BAEP实测值
者1(左)l(右)2(左)2(右)3(左)3(右)4(左)4(右)5(左)5(右)6(左)6(右)7(左)7(右)8(左)8(右)9(左)9(右)lo(左)lo(右)ll(左)li(右)12(左)12(右)
I
1.64
Ⅱ
3.923.883.563.483.84
m
5.445.845.48
i一Ⅲ
2.282,28
Ⅲ一V
1.52
i—V
3.84243.883.684364.24.34.264.o4.123.96
1.印
1.601.641.41.56l461
54
1.96
1.921.84
1.96
1.842.442.082.262.341.842042.322.362062.121.881.922.422.302.242.122.442.32
5.32
5.765.765.76
1.92
2.122.04
3.“
3.723.883.56
5.80
5.725.725,40
1.922.16
2081.64
1.72l601441.481.741.841.681.791501621.521.521.441.361.561.44
3.“
3.763.843.803.963.563.683.923.923.753.643.88
5.525.84
5.84
1.68
2.041.88
4.04
4.14.o3.643:723
5.32
5.485.445.645.56
1.76
1.81.521.721.801.721.441.61.681.84
94
4.024.043.84
5.365.325.285.76
5.88
3.883.92
4.24.44
3.68
4.084.04
2.52
2.6
第1期司俊峰等:听觉诱发电位
・75‘
[1]Zhu
C,Hu
G.Advanceinsingle-trialEP/ERPestimation.ForeignMedicine
onBiomedical
EI】gin∞j119・
2000,23(4):211--216(朱常芳,胡广书.诱发电位快速提取算法的新进展.国外医学(生物医学工程分册),2000,23(4):211~216).[2]Ning
x,Shen
Z,ShenD,eta1.Thedataacquisitionandanalysis(NaturalSciences),1991,27(2):273
284(宁新宝,沈振字,沈德才等.高频心电图的数据采集与分析南京大学学报(自然科学),1991,27(2):273--284).
[3]NingX,Jinz,ShaoJ.Reducingnoisebyintegratedopemtiomlamplifiersinparallel.JournalofDataAe-
quisition&Processing,1990,5(2):59~65.(宁新宝,金之仑,邵军.用多个集成运放并联降低心电放大器噪声.数据采集与处理,1990,5(2);59--65).[4]Lin
K.Principleandacquisition
lishing
House(丹嗍),1978,41.(林可祥.伪随机码的原理与应用北京:人民邮电出版社,1978,
E
offakerandoraccde.Belj'ing:People’sPostandTelecommunicationsPub—
41)
[5]Carlos[63
D,M。haIrl—S
Mobinl.WeightedR.Welghted
averaging
ofevokedpotentials.IEEETrarⅡ-,action
on
BicmedlcalEn西neexing,1992,BME-39:338~345.
HokeM,R卿B,Wickesberg
try.Electroencephalography
av刚ngtheory
and
application
to
dectrieresponse
audiome.
Cllniq【Neurophysiology,1984,57:84--89
[7]PanY.E∞kedpoteatialshmodicatpractice.BelJing:People’sMedicalPublishingHouse(B.vMP)。1998.
(潘映福临床诱发电位学北京:人民卫生出版社。1988).[8]
Ning
D,YanW.TheanalySis
electrw_ardiomignal.JournalofNankingUni雌ity(NalurdSci嘲),2002,38(1):7~13.(宁瓿宝.张
道斌.阎文泰小波变换用于心电信号高频成份的分析.南京大学学报,2002,38(1):7~13).
Steven
X,Zhang
ofwavelettransformationfor
h讪frequency
compolxentsfrcrn
[9JJS,AkramA,MichaelU,eta1.FastwavelettransformationofEEG.Electmencephalography
Neumphysiolcgy,1994,91:442~455.
andC1inical
[10]ZhaoS,XongX-Wavelettramfom】ationandAnal蝴s.BeiJing:PublishingHouseofFAecttlmi,:sIndustry
(PHEl),1996,64.(ig松年,熊小芸.子波变换与子渡分析.北京:电子工业出版社,1996,64),
[11]ThakorNV,Multi[12]Jorge
resolutionwavelet
analyslsofevokedpotentials.IEEETransactiononBiomedicalEn.
sineering,1993,BME-40:l
IA
085~1089.
EvokedpotentialsResearch.IEEEⅡ19irleeringinMedicineandBiology,1992,67~68.
听觉诱发电位检测及分析
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
司峻峰, 宁新宝, 薛方
司峻峰,宁新宝(南京大学电子科学与工程系,近代声学国家重点实验室,江苏,南京,210093), 薛方(合肥工业大学计算机与信息学院,安徽,合肥,230009)南京大学学报(自然科学版)
JOURNAL OF NANJING UNIVERSITY(NATURAL SCIENCE)2003,39(1)4次
参考文献(12条)
1. Jorge I A Evoked potentials Research[外文期刊] 1992
2. Thakor N V Multi resolution wavelet analysis of evoked potentials[外文期刊] 1993(11)3. 赵松年;熊小芸 子波变换与子波分析 1996
4. 朱常芳;胡广书 诱发电位快速提取算法的新进展 2000(04)
5. STEVEN J S;Akram A;Michael U Fast wavelet transformation of EEG 1994
6. 宁新宝;张道斌;阎文泰 小波变换用于心电信号高频成份的分析[期刊论文]-南京大学学报(自然科学版)2002(01)
7. 潘映福 临床诱发电位学 1988
8. Hoke M;Ross B;Wickesberg R Weighted averaging theory and application to electric responseaudiometry 1984
9. Carlos E D;Mohammad S Mobinl Weighted averaging of evoked potentials[外文期刊] 199210. 林可祥 伪随机码的原理与应用 1978
11. 宁新宝;金之仑;邵军 用多个集成运放并联降低心电放大器噪声 1990(02)12. 宁新宝;沈振宇;沈德才 高频心电图的数据采集与分析 1991(02)
引证文献(5条)
1. 高瑞静. 孙迎. 李越囡 脑干听觉诱发电位波形的自动检测及新标识[期刊论文]-中国组织工程研究与临床康复2010(52)
2. 徐灵活. 殷明. 章建程. 唐志文. 李中付. 袁海霞. 罗志辉 噪声听觉损伤易感动物筛选的噪声暴露方案探讨[期刊论文]-中华航海医学与高气压医学杂志 2007(3)
3. 傅娟 脑电生物反馈系统及其研究[学位论文]硕士 20044. 司峻峰 事件相关脑电活动的检测与分类[学位论文]博士 2003
5. 高瑞静. 孙迎. 李越囡 脑干听觉诱发电位波形的自动检测及新标识[期刊论文]-中国组织工程研究与临床康复2010(52)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_njdxxb200301009.aspx