发明名称 基于幅值调制视觉诱发电位脑-机接口方法
摘要 基于幅值调制视觉诱发电位脑-机接口方法,先进行脑电信号采集系统连接,再优化选择受试者幅频响应最为敏感的频率作为载波频率,最大化响应信号的信噪比,然后提出调制波频率的选择范围,再设计AMVEP范式,产生幅值调制视觉诱发电位脑-机接口范式的刺激序列,最终实现AMVEP的单次或少次信号识别,本发明提出了基于幅值调制视觉诱发电位的范式AMVEP,并提出与之相适应的调制信号分析与特征提取方法,实现基于幅值调制视觉诱发电位的范式的AMVEP单次或少次信号识别,具有操作简单、电极数目少和目标数多的优点。
申请公布号 CN102436302A 申请公布日期 2012.05.02
申请号 CN201110259018.X 申请日期 2011.09.02
申请人 西安交通大学 发明人 徐光华;张锋;谢俊;王晶;游启邦;梁晓旭
分类号 G06F3/01(2006.01)I 主分类号 G06F3/01(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贺建斌
主权项 1.基于幅值调制视觉诱发电位脑-机接口方法,其特征在于,包括以下步骤:步骤1,在受试者头部(D)的枕区Oz位置安放测量电极(A),在受试者头部(D)的单侧耳垂位置安放参考电极(B),在受试者头部(D)的前额部Fpz位置安放地电极(C),测量电极(A)的输出端接入脑电放大器(E)的第一输入端(E1),参考电极(B)的输出端接入脑电放大器(E)的第二输入端(E2),地电极(C)的输出端接入脑电放大器(E)的第三输入端(E3),脑电放大器(E)的输出端和计算机(F)的输入端相连,计算机(F)的输出端和计算机屏幕(G)连接形成脑电信号采集系统;步骤2,将计算机屏幕(G)可表示的单频率刺激范式通过计算机屏幕(G)呈现在受试者面前,受试者头部(D)距离计算机屏幕(G)为50~100厘米,依据受试者在不同单频率刺激下的诱发电位响应幅值敏感度规律,优化选择敏感的刺激响应频率作为调制视觉诱发电位刺激范式的载波频率,优化选择原则为:根据受试者的不同,做出此受试者可选频率的幅频响应曲线图,从幅频响应曲线图中选择幅值最大值所对应的频率值作为载波频率,即选择受试者幅频响应最为敏感的频率作为载波频率,最大化响应信号的信噪比;步骤3,根据步骤2所选择的载波频率,以及影响调制波调制频率选择的因素,提出调制频率的选择范围,由于影响调制波的因素包括屏幕刷新率和频谱分辨率,所以载波频率对调制波调制频率必须符合采样定理,调制波调制频率选择的方法具体包括以下内容:幅值调制信号:y(t)=C+C*sin(2*π*f<sub>m</sub>*Δt<sub>m</sub>-pi/2) (1)采样定理:<maths num="0001"><![CDATA[<math><mrow><mfrac><msub><mi>f</mi><mi>ms</mi></msub><msub><mi>f</mi><mi>m</mi></msub></mfrac><mo>&GreaterEqual;</mo><mn>2.56</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>幅值调制信号计算公式中参数的变换关系如下:f<sub>ms</sub>→f<sub>0</sub>,即:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&Delta;t</mi><mi>m</mi></msub><mo>&RightArrow;</mo><msub><mi>&Delta;t</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>f</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mi>n</mi><msub><mi>f</mi><mi>r</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>&Delta;t</mi><mn>0</mn></msub><mo>=</mo><msub><mi>&Delta;t</mi><mi>m</mi></msub><mo>*</mo><mi>&lambda;</mi><mo>&DoubleRightArrow;</mo><mfrac><mn>1</mn><mi>&lambda;</mi></mfrac><mo>=</mo><mfrac><msub><mi>&Delta;t</mi><mi>m</mi></msub><msub><mi>&Delta;t</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>f</mi><mi>ms</mi></msub></mfrac><mo>/</mo><mfrac><mi>n</mi><msub><mi>f</mi><mi>r</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><mrow><msub><mi>f</mi><mi>ms</mi></msub><mo>*</mo><mi>n</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>由以上公式(1)~(4)推导得知最终调制波频率为:<maths num="0004"><![CDATA[<math><mrow><mi>x</mi><mo>=</mo><mfrac><msub><mi>f</mi><mi>m</mi></msub><mi>&lambda;</mi></mfrac><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mi>m</mi></msub><mo>*</mo><msub><mi>f</mi><mi>r</mi></msub></mrow><mrow><msub><mi>f</mi><mi>ms</mi></msub><mo>*</mo><mi>n</mi></mrow></mfrac><mo>&DoubleRightArrow;</mo><mfrac><msub><mi>f</mi><mi>ms</mi></msub><msub><mi>f</mi><mi>m</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><mrow><mi>n</mi><mo>*</mo><mi>x</mi></mrow></mfrac><mo>&GreaterEqual;</mo><mn>2.56</mn><mo>&DoubleRightArrow;</mo><mi>x</mi><mo>&le;</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><mrow><mi>n</mi><mo>*</mo><mn>2.56</mn></mrow></mfrac></mrow></math>]]></maths>其中:C:常数f<sub>ms</sub>:幅值调制信号的采样频率f<sub>m</sub>:幅值调制信号的频率f<sub>0</sub>:载波频率<img file="FDA0000088735520000025.GIF" wi="243" he="129" />:幅值调制信号的采样间隔x:调制波的频率n:表示一个周期所需要屏幕刷新的帧数f<sub>r</sub>:屏幕刷新率由以上推到过程得知调制波频率的范围应该为:<maths num="0005"><![CDATA[<math><mrow><mn>0</mn><mo>&le;</mo><mi>x</mi><mo>&le;</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><mrow><mi>n</mi><mo>*</mo><mn>2.56</mn></mrow></mfrac><mo>=</mo><mfrac><msub><mi>f</mi><mn>0</mn></msub><mn>2.56</mn></mfrac><mo>;</mo></mrow></math>]]></maths>步骤4,根据步骤2优化原则所选择的载波频率,以及步骤3所提出调制波频率的选择范围,设计AMVEP范式,产生幅值调制视觉诱发电位脑-机接口范式的刺激序列,通过计算机屏幕(G)呈现在受试者面前,产生AMVEP范式的刺激序列的计算公式如下:<maths num="0006"><![CDATA[<math><mrow><msub><mi>y</mi><mi>AMVEP</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>C</mi><mo>+</mo><mi>C</mi><mo>*</mo><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mo>*</mo><mi>&pi;</mi><mo>*</mo><mfrac><mrow><mi>n</mi><mo>*</mo><mi>x</mi></mrow><msub><mi>f</mi><mi>r</mi></msub></mfrac><mo>-</mo><mi>pi</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,C为常数,x为所设计调制波频率,且根据步骤3知:<maths num="0007"><![CDATA[<math><mrow><mi>x</mi><mo>&le;</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><mrow><mi>n</mi><mo>*</mo><mn>2.56</mn></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤5,先通过步骤2~3确定调制编码SSVEP范式,再通过步骤4将其呈现在受试者面前,最终通过步骤1所述的脑电信号采集系统获得受试者的AMVEP响应信号,提出与之相适应的调制信号分析与特征提取方法,实现AMVEP的单次或少次信号识别,具体包括以下步骤:◆对采集到的不同刺激序列的AMVEP响应脑电数据进行带通滤波;◆采用Hilbert解调方法对滤波后的数据进行幅值解调;◆采用FFT方法对幅值解调后的数据求其频谱极大值;◆根据频谱结果,得到频谱中的极大值对应的频率值,即为调制波频率值。
地址 710049 陕西省西安市碑林区咸宁西路28号