脑干听觉诱发电位测试中的曲线拟合与参数估计方.docVIP

???? 脑干听觉诱发电位测试中的?曲线拟合与参数估计方 ???? ????　脑干听觉诱发电位是声觉器官受刺激后产生的一种生物电，反映了相应感觉通路及神经元群的活动，可用作听力阈值的客观判断，以及对多种耳疾和脑部病变作出诊断，因而成为当今临床医学和生物医学信号处理中的一个重要课题。本文介绍了脑干听觉诱发电位潜伏期(L)与刺激声强(I)之间关系曲线的作用及其拟合建模，回归算法优化，还用实测数据对选定的拟合曲线模型进行了回归方差分析。　　关键词:　脑干听觉诱发电位；潜伏期；曲线拟合；参数估计；方差分析　　分类号:　R318.04 0　概述　　在研制“脑干听觉诱发电位测试仪”［注］的过程中，需将测得的病人数据用相应曲线来拟合，我们用统计回归方法较好地解决了这个问题。　　声觉器官在一定强度的声脉冲激励下，通过听神经传至大脑皮层，会产生脑干听觉诱发电位(Brain Stem Auditory Evoked Potential, BAEP),见1［1］。　　由1可见，在声刺激后的0～10ms内，出现了7个电反应波(用罗马字母Ⅰ～Ⅶ表示)，它们分别由听觉通道的听神经(Ⅰ)、耳蜗核(Ⅱ)、上橄榄核(Ⅲ)、外侧丘系核(Ⅳ)、下丘核(Ⅴ)、内膝体(Ⅵ)和听区颞横回(Ⅶ)这7个生理解剖部位发出。因此，由这些电反应波组成的脑干电位是人耳听觉通道正常生理和病理的一种可靠指征［2，3］。　　当声刺激强度不同时，脑干电位波形也产生相应的变化。各峰的峰值幅度与声强成正比。各峰潜伏期(latency)随声强(intensity)降低而延迟，且幅度有所减小。2示出了从正常人身上测得的不同声强时脑干电位各峰潜伏期的变化［4］。1　正常人的脑干听觉诱发电位 　　由2可见，Ⅴ峰最为明显。不仅其幅值最高，且最为稳定。当其它峰都消失时，Ⅴ峰仍依稀可辨。因此，Ⅴ峰潜伏期就具有较高的诊断价值，临床上常用它来确定听力阈值。如2，此人的听力阈值就为10db。［5］2　正常人各种声强度的脑干电位的潜伏期变化 　　尤其是Ⅴ峰潜伏期(L)对声强(I)的关系，临床意义就更大，可用作诊断听觉系统疾病及脑干部位病变的依据。　　将各声强下Ⅴ峰潜伏期值光滑连接，就形成了一条L-I关系曲线(如3中粗实线)。当病人患有听觉传导阻滞或因肿瘤压迫听神经引起听力损失时，所测得的Ⅴ峰潜伏期较正常人滞后(见3中虚线Ⅰ)。而某些神经系统有病变的患者(如美尼尔氏症)，其Ⅴ峰潜伏期在高声强时滞后，在低声强时却较正常值超前，因而其L-I曲线较正常曲线陡直(见3中虚线Ⅱ)。L-I曲线可以作为诊断某些听觉系统疾病及脑干部位病变的依据。因为直接观察2只能定性了解，不便定量比较分析，所以用合适的曲线模型来拟合测试中所采集到的病人潜伏期数据，就具有较高的临床价值。拟合曲线参数可以用作在参数空间进行定量聚类分析，见4。如进一步应用模糊神经网络专家系统，还能自动给出诊断意见［6］。如果要自动提取Ⅴ峰潜伏期，还必须对2中各条曲线进行平滑去毛刺，最高峰点定位等预处理操作，这里不详述。3　正常人曲线与几种病人曲线的比较4　参数空间聚类分布示意 1　拟合建模　　从大量统计数据来看，L-I曲线为一单调降非线性曲线。数学上可用于拟合这类曲线的函数有许多，如多项式函数y=a０+a１x+…+akxk,倒数函数y=a+b/x，指数函数y=aebx(b0),对数函数y=a+blogx (b0)等。我们用多组实测数据，对上述几种函数拟合L-I曲线进行了拟合误差比较。所比较的是最小二乘拟合情况下的残差平方和　　　　式(2—1)中i是根据拟合函数得出的y值，yi为实测值。i因拟合函数不同而异：　　多项式函数：i=a0+a1xi+…akxki　　倒数函数：　i=a+b/xi　　指数函数：　i=aebxi　　对数函数：　i=a+blogxi　　求最小二乘拟合情况下残差平方和的步骤是［7］：　　(1) 用变量替换方法把拟合函数化为一元或多元线性函数　　　　　　i=a0+a1u1i+…+amumi　　　　　　　(2—2)式(2—2)中uki是实测值xi在变量替换后的映射值：　　多项式函数拟合，uki=xki,k=1,2,…,m　　倒数函数拟合，u1i=　　对数函数拟合，u1i=logxi　　指数函数拟合情况有些不同：对式i=aebxi两边同时求对数，得：lgi=lga+bxi令a0=lga,a1=b,u1i=xi,yi也映射为lgyi,S=　　(2) 把式(2—2)代入式(2—1)，对各项系数求偏导：　　　　令　可得m个m元线性方程。求解这个线性方程组，得到各元变量拟合系数a0,a1,...am，根据需要作逆变换，得出实际拟合系数。　　(3) 用式(2—2)，式(2—1)计算残差平方和S。