[发明专利]降低TMS电磁脉冲对大鼠运动诱发电位影响的工频滤除方法有效
| 申请号: | 201810524384.5 | 申请日: | 2018-05-28 |
| 公开(公告)号: | CN108888266B | 公开(公告)日: | 2021-05-14 |
| 发明(设计)人: | 王欣;殷涛;刘志朋;王腾飞 | 申请(专利权)人: | 中国医学科学院生物医学工程研究所 |
| 主分类号: | A61B5/395 | 分类号: | A61B5/395;A61B5/397 |
| 代理公司: | 天津市北洋有限责任专利代理事务所 12201 | 代理人: | 杜文茹 |
| 地址: | 300192 *** | 国省代码: | 天津;12 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 降低 tms 电磁 脉冲 大鼠 运动 诱发电位 影响 工频滤 方法 | ||
1.一种降低TMS电磁脉冲对大鼠运动诱发电位影响的工频滤除方法,其特征在于,包括如下步骤:
1)将单次TMS脉冲作用于大鼠运动皮质,采用1000Hz以上的采样频率记录大鼠清醒状态下的后肢表面肌电信号,表示为x(n),n=1,2,…,N,N为肌电信号数据点数;
2)对肌电信号x(n)进行傅里叶变换,将工频干扰和直流分量对应的傅里叶系数置零,其余傅里叶系数保持不变,进行傅里叶逆变换,实现工频干扰的滤波,得到滤波后的肌电信号,表示为y(n),n=1,2,…,N,N为肌电信号数据点数;包括:
(1)对肌电信号x(n)进行傅里叶变换,得到W个频率点对应的傅里叶系数X(k),
其中,W为频率点数,W=N;
(2)计算W个频率点对应的真实频率值f(k),
f(k)=(k-1)·Fs/W,
其中,Fs为采样频率,大于1000Hz;
(3)将50Hz工频干扰及50Hz倍频和直流分量对应的傅里叶系数置零,其余傅里叶系数保持不变,得到新的傅里叶系数Y(k),
其中,f0为直流分量、50Hz及50Hz倍频,f0={0,50,100,150,200,250,300,350,400,450,500};
(4)对新的傅里叶系数Y(k)进行傅里叶逆变换,得到滤除工频干扰后的肌电信号,表示为y(n),
3)从滤波后的肌电信号y(n)中提取出运动诱发电位,得到滤除工频干扰后的运动诱发电位Pm,m=1,2,…,M,M为所记录的肌电信号中的运动诱发电位个数。
2.根据权利要求1所述的一种降低TMS电磁脉冲对大鼠运动诱发电位影响的工频干扰滤除方法,其特征在于,步骤1)中所述的将单次TMS脉冲作用于大鼠运动皮质,是采用Magstim rapid2 8字线圈对大鼠运动区进行刺激,TMS刺激模式为单脉冲刺激,刺激强度为最大输出的60%,连续输出30个脉冲,脉冲间隔大于1s。
3.根据权利要求1所述的降低TMS电磁脉冲对大鼠运动诱发电位影响的工频干扰滤除方法,其特征在于,步骤3)包括:
(1)根据滤除工频干扰后的肌电信号y(n)幅值大于20000μV且为极大值,识别出肌电信号中的TMS电磁脉冲,TMS电磁脉冲对应的时间点表示为Tl,l=1,2,…,L,L为TMS电磁脉冲个数,
Tl=n/Fs if y(n)>20000y(n)>y(n-1)y(n)>y(n+1),n=2,3,…,N-1,
其中,Fs为采样频率,大于1000Hz;
(2)以时间点Tl为标记,在TMS电磁脉冲后150ms内识别出滤波后的肌电信号y(n)的运动诱发电位,运动诱发电位峰峰值大于100μV,共识别出M个运动诱发电位,以TMS电磁脉冲对应的时间点作为标记,记为T'm,T'm∈Tl,m=1,2,…,M,M≤L;
(3)从滤除工频干扰后的肌电信号y(n)中截取时间点T'm前100ms至T'm后200ms的肌电信号,为含TMS电磁脉冲的运动诱发电位,从中读出运动诱发电位波形、幅值和潜伏期。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国医学科学院生物医学工程研究所,未经中国医学科学院生物医学工程研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810524384.5/1.html,转载请声明来源钻瓜专利网。





