基于LabVIEW的听诱发反应多功能声刺激器
来源:99网
研究论著I THESIS&RESEARCH REPORT ・1 3・ 基于LabVIEW的听诱发反应多功能声刺激器 戴红娅,董 洁,辜承慰,王 涛 【摘要] 目的:开发设计一种基于LabVIEW的可编程高精度声刺激器,可连同脑电记录设备构成一个记录听诱发电 位通用的系统。方法:采用NI公司的M系列DAQ数据采集卡,通过硬件时钟控制方法产生时间高精度的声刺激系 统;LabVIEW编程为用户提供多种声刺激产生模式,不但能提供通用的一般波形,同时产生能灵活控制的任意波形。 结果:该系统可以把常用的短声、纯音等和用户自定义的任意波形提供给标准的传声器(耳机、喇叭);数字触发信号 与模拟声信号精确地同步;声音选择参数、刺激序列和声音强度等都可调。结论:和采用声卡输出方式相比,该系统的 输出声音波形失真度大大减小,并且采用硬件时钟控制的方式保证了输出触发控制信号的时间同步性:该设计便于 系统功能的升级和扩充。 【关键词】听诱发电位;声刺激器;LabVIEW;硬件时钟控制 [中国图书资料分类号]R339.16;Q437[文献标识码]A[文章编号] 1003—8868(2010}09—0013—05 Multi-functional Sound Stimulator for Auditory Evoked Potential Recording Device Based on Lab VIEW DAI Hong—ya 一,DONG Jie ,GU Cheng-wei ,WANG Tao (1.College of Biomedical Engineering,Nanfang Medical University,Guangzhou 5 105 15,China;2.Guangzhou Yishixing Electronic Technology Co.,Ltd,Guangzhou 510515,China) Abstract Objective To develop a high—accuracy programmable sound stimulator based on Lab VIEW platform,which can be used in connection with a third party electroencephalogram recording device to construct a general function system recording auditory evoked potential(AEP).MC,ho ̄A M-serial DAQ acquiirng card produced by National Instruments Company(NI)was adopted to produce a high—accuracy timing acoustic stimulation system by using a hardware clock control strategy:the Lab VIEW coding was accomplished to provide users a variety of sound generation modes including both standard waveforms for general use and arbitrary waveforlTiS that could be controlled freely.R鲫ll协This system was able to provide common clicks,pure tone sound stimulation etc.,and user—defined arbitrary wavefornls for standard output transducers,such as earphones and speakers.The digital trigger channels and analogous sound signals could be precisely synchronized.The parameters defining sound selection,stimulus sequences and sound intensity were adjustable. Conclmion In contrast to the output of conventional sound card,this system is capable of delivering less distorted waveforms,and ensuring synchronization of the trigger control signal by using the hardware clock in a controlled manner. This design secures a reliable approach to system update and extension.【Chinese Medical Equipment Journal, 2010,31(9):13-17】 Key words auditory evoked response;sound stimulator;Lab VIEW;hardware clock control 1 引言 提取AEP信号:对AEP信号进行显示、储存、记录和分析l】1。 其中,声刺激器是影响AEP记录的主要因素之一,其功能主 要包括2个部分:产生特定刺激声信号,通常信号需要满足一 定的时域或频域参数以及保真度的要求:提供同步触发信号, 根据感兴趣AEP成分的不同.要求触发信号满足相应的时间 精度。目前,国内大型医院和科研机构使用的听诱发电位仪绝 听诱发电位(auditory evoked potentials,AEP)是一定强度 的声音刺激听觉系统时,听觉系统发生的一系列与刺激有严 格时间对应关系的电反应。AEP记录的基本技术处理流程 为:信号发生器(声刺激器)产生刺激声信号,经换能器、信号 传递装置将适宜的刺激声信号传递给受试者,同时触发平均 叠加装置启动:通过记录电极采集声刺激诱发的电生理学反 大多数为进口仪器,如美国智听公司听诱发电位仪SmartEP、 美国Nicolet Spiirt脑诱发电位仪等,此类设备提供声刺激的 方式都是根据选择刺激声类型(clicks、AM tones、sine wave、 burst tones、pips tones)和刺激参数(频率、强度、持续时间等) 输出刺激声信号与同步触发信号,能产生有限的刺激声。随着 科研工作的深入,刺激声信号越来越复杂。很多实验方案需要 更加精确且灵活可控的刺激声发生器。例如Bell等报告了采 用调频型chirp声的刺激方法有助于AEP在麻醉深度评估应 应;对原始反应信号进行放大、滤波处理:通过平均叠加处理 基金项目:国家自然科学基金(60771035) 作者简介:戴红娅(1982一),女,重庆人,硕士研究生,主要研究方向为生 物医学工程与仪器,E—mail:dhycumt@163.con;王涛(1967一),男,安徽阜 阳人,教授,硕士研究生导师,主要从事生物医学工程与仪器方面的教学 与研究工作。 作者单位:510515广州 南方医科大学生物医学工程学院(戴红娅、董 洁、辜承慰、王红娅) 涛);510515广州广州易事兴电子科技有限公司(戴 ・医疗卫生装备-2010年09月第31卷第O9期 Chinese Medical Equipment Journal・Vo1.31・No.09・September・2010 ・1 4・ THESIS&RESEARCH REPORT I研究论著 用 。采用具有伪随机脉冲型的最大长序列(Maximum length sequence,MLS)方法提高AEP的刺激率,存听觉系统的线性 及非线性特性分析巾发挥了重要作用 .近年来又有一些新 的高刺激牢方案捉…141另外,稳态反应研究巾也呈现很多新 读取声信号的方法。.wav文什可以是用户根据需要录制或日j软 件编写的声文件,其前而板如图3所示。各参数意义如下: 的刺激疗案,如40 Hz匏怎反应[『1.用一组间隔不等的clicks 做刺激声 、这 刺激方案 临睐常用的仪器难以实现 目前Binsemi公可生产的高精度、高采样牢的直流脑电 信 放犬器直接采用I,abVIEW作为信号采集的编程1 具.提 供源代码开放式接口,支持恼电信号采集模块的二次开发功 能。但是不提供相应的刺激器,需要 行开发或者第j方购 置。本文的主要目的是开发设计一种基于LahV1E ̄ 的可编程 高精度声刺激器,能够提供常规的各类典型卢信号,同时用户 还可以根据需要选择不同的方式 成任意波形的 信号。随 着研究进行,此系统还可以不断升级,保证技术的九逊性. 2结构 图3.way文件格式声刺激器前面板框图 本研究采用NI公司LabVIEW软 和DAQ倾f÷ 脱仟 意波形声刺激器,其结构如图1所,J 为r尽可能方他全而 地提供声刺激,同时又能针埘不同的用户群,本研究设 r 3 (1)文件路径:刚户需要播放的.wav文件的路径: (2)预读取总数:灵活设置需要读取.wa、义什的样本总 种声信号选择方式:典型波形、.wav文件和.dat文件。 数,其中“一l”代表读取整个丈什; (3)"始位置:设置从文件的哪个位置处开始读取,以“中羊 传,*器(耳机) 触发信号几 本点”为单何: (4)渎取样本数: 示实际渎取的样本数, F确的结果应 该是和预读取总数卡日一致(预读取总数如果设置为“一l”.此 } 参数显示为整个文件总的样本数: (5)采样率:显示.wav文件的采样率: 图1 LabVIEW环境下听诱发反应声刺激器的组成 2.1 典型波形 (6)图像框:显示所渎取的.wav文件样本波形: (7)强度:设置输出声信号强度: (8)刺激间隔:没置读取的波形(图像框巾所示波形)输出 时的l1_fJ隔: (9)刺激模式卜批框:连续(直到丁动按“停止”按钮结 听诱发电位主要用的典型波为:clicks,AM tones,sine waves,bursttones和pipstones,每种波形都有不同的参数,具 体参数以树形图表示,见网2 Cliks(持续时间、频率、强度、极性) 强度) 、束)、单次(只播放1次)、有限次刺激(输入需要播放的次数 N)。 岗 期 序 列 、强度) 2.3更多文件 、F降时问、时程) 下降时l训、时程) 对于输出任意波形的声J}{『J激.J1j.wav艾件方式就完全能 实现,但.w.d、,文件一般鄙和计算机的声卡功能一+致,采佯率 自 定 义 序 列 刺 激 个 数 一 大多为44 10O、22 050、l1 025 Hz等3利IⅢ 为追求更高质量的 声信号,有的实验需要达到48 kHz以致更高的采样率.所以 2 本系统设置了另外可南川户A定义生成的.txt、.dat文件,系 长度 统告知用户相应文件的具体格式.不管用户用什么方法保存 声信号,此系统郜可以播放: 图2声刺激器典型波形参数选择示意图 3原理 如图2所示.最简单的是周期性序列.选定波型。然后再 前面描述的3种模式声刺激 实现方式上原理是相通 的,此处以.dat文件方式为例详细描述。.dat文件口f存储任意 类型数据,因此本文 二进制义什的头部增加『,简 的头文 选择相应波型参数就完成声信号的设置:除了周期序列.用户 还可以根据实际需要编辑一个自定义序列。首先输入序列中 包含的刺激个数,然后分别对每个刺激进 定义,每个刺激都 可以根据作者需要选择波型及相应参数.同时两两相邻的刺 激问间隔ISI由用户自行设置,用户另外设置好序列的长度 件,用来捕述文件的组成,即数据的采样率和数据长度 实际 应用【{1,用户可以根据需要增加相应的头文件信息,如文件名 称、时间等。图4所示为本系统设计的二进制数据结构,头文 之后,就实现一个由自己定义的波形。 2.2.wav文件 件中.按着顺序存放的数据是:字符串“sample rate”、采样率 (整型32 bits)、字符串“length”、长度值(整型32 bits), 文什 之后,以float32类型存放声音数据. .除了“典型波形”外,“.wav文件”为用户提供了另一种灵活 N-疗卫生装备・2010年09月第31卷第09期 Chinese Medical Equipment Journal・Vo1.31・No.09・September・2010 研究论著I THESIS&RESEARCH REPORT ・1 5・ 图5所示一: 进制文件刚刚被打开时,文什指针指向第1个宁 节.按着字节顺序,标注1模块以字符为单位读取l1位,标注 san ̄Die rate length 2模块读取以32 bits为单位的1个整数,标注3模块读取6 声音数据 支f, 字符,标注4模块读取1个32 bits的整数,标注5模块从 指针当前位置到文件结尾以32 bits单精度为单位读取 音 数据.对比图2中的数据格式,标注1~4对应于前面板(网4) 图4刺激声的二进制文件数据结构 3.1 声信号读取 .显示的刺激声头文件的信息,包括:“sample rate”及其相应的 数值(整数型);“length”及其相应的数值(整数型)整数。标注 5对应前面板图形框显示出刺激声的时域波形(单精度型数 组) dat文件是字节的顺序排列,用户需了解文件的数据结 构才能通过指针读写文件中的字节。与图4所示的数据结构 一致.LabVIEW环境中对本文规定的二进制文件读取框图如 图5访问二进制文件的Vl程序框图 值得注意的是:LabVIEW巾读取二进制节点有一个重要 的输入选择项“bvte order”(字节顺序),也称之为大小端问题。 LabVIEw的写二进制文件节点也存在同样的选择项,因此, 存读写 数巾保存此选择项一致就可以了。如果是用其他编 程语言写的二进刹文件,就要选择IF确的大小端顺序。如采用 Matlab编写.dat文件.字节顺序默认与操作系统一敛,Win— DI Samole Clock (AI Start Frigger) AI Sample Clock 图6 AO与DO同步输出时钟触发原理图 dows操作系统采用“little—endian”顺序,所以大小端参数选择 “little—endian” 南(mute)一个时钟作为DO Sample Clock。最简单直接的方式 就是选用板卡内部AO Sample Clock作为DO Sample Clock,经 实验.此方式在采样频率低于30 kHz时是可行的,为追求更高 的声音信号质量,本实验所设计的采样频率不低于48 kHz,所 以需路南板卡内部其他更高频率的时钟做采样时钟。 M系列板卡AO、D0采样时钟能同用的内部时钟源还可 以是板卡的Counter nlnternetOutput.上作原理与图6一样, 3.2 同步触发 同文件格式程序有不同的读取方式,如3.1一节为二 进制格式的读取方式.目的是渎取声信号和相应的参数,读取 完数据之后,声刺激器最灭键的一个问题就是触发信号与声 信号的同步输出 能够实现声信号与触发信号同步的方法很 多,如常见的心理学软件E—Prime、LabVIEW,实现声卡输出 声音信号和并L]输…同步触发信号,精度最高能达Ins级。听 诱发反应记录的是声刺激后的瞬态反应且持续时间非常短 由硬件触发信号开始计数器时钟生成,此时钟路由为AO、DO 的采样时钟 由于计数器时钟信号是直接对板卡的20 MHz 时基分频得到.能文现的时钟频率比板卡自带的AO Sample Clock要高,能实现更高的采样率,如96 kHz甚至更高。 FIFO支持重传输模式:FIFO中的所有样本点被采样时 钟读完后,FIFO又以相同的顺序输出所有的样本点。例如, FIFO包含5个样本点,FIFO输出产生的模式是:#1,#2,#3, 暂,同步触发信号Ills级精度远远达不到要求,如在耳蜗电罔 巾,只需要记录刺激后0~2 ms的波形。因此,为提高两路信号 的高精度同步.需要用硬件时钟方式触发两路信号同步传输。 以下以“连续”播放模式为例具体介绍。 3I2.1 原理介绍 首先把2组等长的样本(声信号和触发信号)分别写入板 卡的模拟FIFO(buffer)和数字FIFO中,2组样本以相同的时钟 触发信号上升沿或F降沿开始采样,并以相同的采样时钟输出 #4,#5,#l,≠}2,#3,#4,#5,#1等。写入模拟FIFO中的数据是经 过DAC(数模转换)的模拟电压信号:写人数字FIFO中的数 据代表的是1TrL信号.同时可以进行编码,如Port 0由8路数 字线组成,假如写入FIFO的某个数据为5(二进制表示为 00000l01),则线7—0同时刻输出为0000010l,8线能提供 2s_1种编码形式,具体选用哪种形式需与采集系统结合。如 样本.即一个采样时钟脉冲生成一个样本,分别从DAQ板 的 AO和DO通道输 、原理如图6所示,AO Start Trigger脉冲上 升沿触发AO任务,AO Sample Clock南AO Sample Clock Timebase分频得到,AO Sample Clock Timebase又是由板卡20 本系统只选用线口PO.0输出触发信号,写入数字FIFO中的 采样点“1”代表高电平,“0”为低电平。 3.2_2实现框图 MHz时钟基分频得到,每个AO Sample Clock输 一个样本 点。因为本系统所用的板卡没有直接由板卡时基分频得到DO SampleClock,所以必须从外部信号或板卡内部其他子系统路 在LabVIEW中实现模数同步的框图如图7所示。图中3 ・医疗卫生装备・2010年09月第31卷第09期Chinese Medical Equipment Journal・Vo1.31・No.09・September・2010 ・1 6・ THESIS&RESEARCH REPORT l研究论著 个“DAQmx任务开始”节点和1个“Merger Error.vi”节点的连 接方式保证模拟输出和数字输出在采样时钟生成前就准备 好,并且在同一时刻开始。最后,因为模拟输出和数字输出任 务选择的模式都是“continue”.while循环实现板卡FIFO中数 据的连续周期性输出,直到“停止”按钮响应。循环框后面是对 各个任务的清除,在任务停止后,清除任务,释放任务保留的 资源,避免继续占用内存。程序中最后用到的错误提示对话 框,程序中出现任何的错误,都将以对话框的形式显示出来。 图7 AO和DO端口实现同步的VI程序框图 4结果 国内国外都有很多人研究听诱发反应声刺激器,有人借助 于计数器、D/A、单片机等,预先把波形存储在单片机片内ROM 实验一个3阶的长度为7的MLS序列『1,l,1,0,O,1, 01,以连续方式播放,每个MLS序列对应一个触发信号(数字 脉冲上升沿),序列中“1”代表有刺激click出现,“0”代表没有 刺激.MPI(minimum pulse interval,2个临近的值“1—1”或“1— 0”之间的间隔)定义为10 ms,采样率48 kHz。此序列可以用自 定义序列、.wav和.dat的方式实现,如图8(a)为示波器上显示 的实际MLS波形(采用.dat方式),示波器通道1显示声信号波 中,因存储空间有限,最终产生的声信号种类只有clicks和sine 信号[61。还有人针对稳态反应利用DSP++(C、C++编程实现控制 界面,软件固化在硬件中)的方式产生AM 。它们和现在临床上 普遍使用的声刺激器一样存在一个共同的特点:采用固化软件 方式.能提供的声刺激种类有限。本系统是在硬件开放的基础 上,通过设计灵活的软件模块来实现功能,开发效果高。 为了查看本系统声信号的质量.与处理声音最成熟的声 形,通道2显示触发信号波形,通道1中一个周期中的脉冲间隔 依次是:l0、l0、3O、20 ms,与MLS序列【1,1,l,0,0,1,0卜一 致,通道2中的高电平出现的位置与MLS序列第1个刺激出现 位置完全同步。另外,图8(b)显示载波2 kHz、调制波85 Hz的 AM波形.图8(C)为4个不同载波、调制波频率的合成波形,图8 (d)显示连续播放.way格式实际图形。 5 分析 卡作比较。听诱发反应中,click¥1l激脉宽都是非常短,一般只 有零点几毫秒,脉宽越窄频带越宽。图9(a)为本系统输出 0.312 5 ms脉宽的实际波形;而声卡有一定的频率响应范围, 相当于一个滤波器.所以9(b)中声卡输出脉宽0.312 5 ms脉冲 波形的失真度比前者大,从图中还可以看出后者脉冲上升时 间比前者大,在听诱发反应中,上升时间短能使神经冲动同步 (a)3阶MLS波形和相应同步数字信号 (b)AM波 -医疗卫生装备-2010年09月第31卷第09期 Chinese Medical Equipment Journal・Vo1.31・No.09・September‘2010 研究论著l THESIS&RESEARCH REPORT ・1 7・ (C)4个不同频率AM波的合成 (d).wav格式波形连续输}{J 图8示波器显示波形 ■∥ “ 瓣 瓣 嚣 嚣 的困扰,同时也避免了一系列可能出现的问题,如硬件电路的 舞 ㈦ * 时间配合和信号抗干扰等,在开放架构的基础3c ̄,3建用户自 定义的测试系统。这种方式不仅可以提高工程质量,减少开发 周期而且界面比较美观、接口很方便,是一个高性能、高精度 嚣 瓣 的测试运载平台,功能升级方便。 攀 赫攀 《 ・● I # 掌囊女●●● 另外,听诱发电位研究涉及到信号源、换能器和信号处理 等,须进行严格的校准方能使用,I临床上具有实际意义的是校 准刺激声。除了检验刺激声强外,还需要校准其他测试参数, 如极性、持续时间等,需要使用比较复杂的电声学设备和专业 知识 因此,本系统中的声刺激器在应用于临床和科研之前, 还需要经过专业的校准。 【参考文献】 [1】韩德民,许时昂.听力学基础与临床[M】.北京:科学技术文献出 (a)本系统输出脉宽0.312 5 ms脉冲波形 版社,2005:330—335. [2】Bell S L,Smith D C,Allen R,et o1.The Auditory Middle Latency Response,Evoked Using Maximum Length Sequences and Chirps, as an Indicator of Adequacy of Anesthesia[J].Anesth Analg,2006, 102:495—498. [3】Burkard R,Shi Y,Hecox K E.A comparison of maximum length and Legendre sequences to derive BAEPs at rapid rates of stimula— tion[J].Joumal of the Acoustical Society of America,1 999,87(4): 1 665—1 672. [4】JewettDL,CaplovitzG,BairdB,et .The useofQSD(q—sequence deconvolution)to recover superposed transient evoked—responses[J]. ClinicalNeurophysiology,2004,115(12):2754—2 775. f5】Gutsehalk A,Mase R,Roth R,etⅡf.Deconvolution of40 Hz steady— state fields reveals two overlapping source activities of the human (b)声卡输出脉宽0.312 5 ms脉冲波形 auditory cortex[J}.Clinical Neurophysiology,1999,1 lO(5):856— 868. 图9本系统和声卡输出脉宽0.312 5 ms的示波器波形图 性好 通过比较其他类型的波形,都可以发现本系统产生的 [6】赵仕波,罗耀华,赵文华.听觉诱发电位仪的设计与实现方案[Jl_ 仪器仪表学报,2008,29(3):594—599. [7】 Martin H R,Romao M,Plaeodo D,et a1.Stimulator with arbitrary 声信号比声卡输出的失真度小,因此本系统比声卡更适合做 听诱发反应的声刺激器。 6 结论与讨论 wavefornl for auditory evoked potentials[C]//16th Argentine Bio engineering Congress and the 5th Conference f oClinical Engineering 本系统的实验结果表明,本刺激器能够弥补传统固定仪 器不灵活的缺点,根据用户自定义产生声音,同时产生完全同 Journal of Physics,Argentine,September 26—28,2007. [8]李兴启.听觉诱发反应及应用[M].北京:人民军医出版社,2007: 9O一97. 步触发信号,满足研究的需要。LabVIEW与DAQ的相互兼容配 套,使系统稳定度和可靠性得到保障,还避免了编写驱动程序 (收稿:2010—02—24修回:2010—03—30) ・医疗卫生装备・2010年09月第31卷第09期Chinese Medical Equipment Journal・Vo1.31-No.09・September‘2010