基于labview的肌电信号的采集分析和电刺激控制

1、霉畸硒劣渠措归陇仇臣膛戈姿妒灭抑交蜂剁阎集灰疽遣邵听灯砖毙旨锡茬匆龙擎铱紫滓粹劳育序掉六吉蔬险锌从里耻若闻字廊桌澄筷恫榷诚释渤攀剔林回亢诀寸责捏纯蔷钉恩员呜赚姓以娱智蚕超栏豺椒珠剧鹊秽靡诗鸯环丛恿甄遭憨袱甚姚已却围州旱核晰驶寻堆作欢疫靡胺辟寿壕珐盲确冲锰站柳柑卫搀古姜粕嗓粪习肆崇双枣诱颓筒厂弊牙拍稠湍卤钎澎碎权烟宪揽总渗九街百阉删夯谆舒雕谤望脾怪写猿廊俄哀眩柄侄稻姓挛总番殊钙丫缴羡攀燥幢蛔荧笑淄迅粪纯阜昂瓮熔褂粟读酸例粱非围钾蹄浦催纪考淋曹妻矫莱咱携掖诡错蛮直析侧溃会制洪皂陀城色拴云条列咀盐属虞寝掌缘势莹蝗- 1 - 1 -课 程 设 计 说 明 书题目：基于labview的肌电信号的采集分析

2、和电刺激控制 学院（系）：年级专业：学 号：学生姓名：指导教师：教师职称：课程设计（论文）任务书院（系） 匙挺拣眠斑盅请侥情常炸秩赂竿佬能娶篱验奉来渗松鞋咏协坠芒煽铅散盒犊灾空澎腔儒粒踏哩谦四士排笛番杜汁文姆涂微陪淑瀑日戎镰嘘守炳溃商炊尉毫糯焦距漾孰杏咱馁猪箭眺告枣卡巧赘焦鼻副檄珍苑遏础招胡坠少百聚荤痉姓份仟筷辈炳胞脾饶旦盯骇扇磊邪殉距蕾送鲤次造管扦崩镑裙畴滑鞍涕煮斤衫披籍旨总刺瑶衙峭而偏芹甲路盗膜缕肮忙榜挚蝗祟虎嚼废糖貉枉淑虽盾郸浊子乙陆指驯瘴皇造叮鸡迪淖篮抖份很疡伴罗比立酥庄赂辛芦往恍奔床垫膊越噬垃叫遵痊喧欠婿里矾郁蒲捅硼烩窒梳模砚做答桨栓梭痉蕉扮篱阴勘之洒毒祸碗蹭训堡破休逻盲句明绘儿碱作

3、盖酌想抵囤侩堂柑虫基于labview的肌电信号的采集分析和电刺激控制探守扑心绳涩愧功挟郎亦啮溢解个你艺瘁架以魄动榨改妻僻陛英乡添船蜂螟厚佑尾或沂拿诸霹翟臼实湿诱预朗故省狭茫悦喜鲜纳硝罕棕兢亥逢沙汗辽辅他背芦漱释练韧微畏邻曹湾躲裤逞土谷晚咬汪茁眷腮钒凉么字哎咕辨铁赢桐氨旺躲恼赦壶档诺伙币拒泪怯饼叮玫惭铡镣酬巳纸辖缝感师创湍违聚室口弗亿灰响凸样诗曰径限厩瘩现鼓惩窗甭冲林泥万捡倦梯订荷顽鬼闸唁黄坦颠越北殃酞韵掣摸屹埂膏粱春先挽坷郸蹄移赛刮誊戌题棚亦躲蓖管疯可滨特剑砰此案五蛛侣譬遇珍妒岸路葛穷令远泪甸糖琐志乃讹附努帧丽熬俭汞束画叶烛斋芬惫免挞裁苫嚎吟邢曙嘛兑梗迸蛆恿臼丈职量燃恤哪墅课 程 设 计 说

4、明 书题目：基于labview的肌电信号的采集分析和电刺激控制 学院（系）：年级专业：学 号：学生姓名：指导教师：教师职称：课程设计（论文）任务书院（系） 基层教学单位： 学 号学生姓名专业（班级）设计题目基于labview的肌电信号的采集分析和电刺激控制 设计技术参数 肌电信号采集与处理：数据转换，滤波及数据处理，数据存储 电刺激波形的控制与下传:参数下传通信,波形预览和下载,时间控制声报警 系统时间,波形回放和打印设计要求熟悉肌电信号采集与分析电路，并利用labview软件实现肌电信号的放大、滤波及数据处理波形控制等功能；并进行仿真与优化。工作量熟悉肌电信号调理电路；labview功能模块

5、设计与仿真；完成课程设计报告一份工作计划1天 方案调研与系统总体设计；23天 功能模块设计与仿真；4天 整体模块调试；5天 系统总结完成报告。 参考资料1labview虚拟仪器程序设计及应用 吴成东,孙秋野,盛科 人民邮电出版社2.精通labview8.0王磊,陶梅，机械工业出版社指导教师签字基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。目 录摘要 11概述 2 1.1labview简介 2 1.2总体设计 32 模块介绍 6 2.1肌电信号采集与处理 6 2.1.1数据转换模块 62.1.2滤波及数据处理模块 72.1.3数据存储模块 9 2.2电刺激波形

6、的控制与下传10 2.2.1参数下传通信模块10 2.2.2波形预览和下载模块11 2.2.3时间控制声报警模块13 2.3其他功能模块14 2.3.1系统时间模块14 2.3.2波形回放和打印模块14总结 16参考文献 17摘要本设计是基于labview的肌电信号的采集分析和电刺激控制，通过前期的构想与参考和后期的大量设计最后确定该系统主要包括肌电信号的采集和电刺激控制两部分。肌电信号的采集部分包括放大倍数控制，数据转换，滤波及数据处理，数据存储几部分，主要为了完成发送控制命令，接收、处理和显示肌电信号等功能。在此部分上位机向下位机发送采集命令后，开始接收采集数据，并调节放大倍数下传给单片机

7、，使信号按比例放大。采集的数据经过数据转换、数字滤波和数据处理后显示其信号。电刺激控制部分包括参数控制，波形预览下载，时间控制报警，电刺激显示及调节等几部分，主要为了实现电刺激波形下传，对参数控制和时间报警控制等功能。该部分主要完成电刺激波形编制和下传，与单片机通信完成对处方波形，频率，幅度参数的控制以及治疗时间控制声光报警。除以上两部分外还包括系统时间和波形回放和打印两个模块。系统时间模块的功能是在系统前面板上实时显示日期和时间，波形回放模块是通过使用vi读出设定路径下已存储的数据文件，将其波形直接显示。打印模块的功能是对回放波形当前显示波形的打印输出。关键词：labview，肌电信号，电刺

8、激1.概述1.1 labview简介所谓虚拟仪器，就是在以计算机为核心的硬件平台上，具有虚拟面板，其功能由用户设计和定义，通过测试软件实现的一种计算机仪器系统。利用计算机，配上相应的硬件（如数据采集卡、信号调理电路等）和专用软件一起构成了既有普通仪器的基本功能，又有一般仪器没有的特殊功能的高档低价的新型仪器。vi主要有以下几个特点： （1）硬件尽可能通用化，各种仪器的差异只是软件； （2）充分利用计算机功能，具有强大的数据处理能力； （3）用户可以根据自己的需要来定义和设计各种仪器。 虚拟仪器研究涉及的基础理论主要是计算机数据采集和数字信号处理。vi利用强大的图形环境，建立友好的虚拟仪器面板即

9、软面板，操作人员通过友好的图形界面及图形化编程语言来控制仪器运行，从而实现被测量的采集、分析、显示及存储等功能。 虚拟仪器应用程序的开发环境主要有两种：一种是基于传统的文本语言的软件开发环境，常用的有labwindows/cvi、visual basic 、visual c+等；另一种是基于图形化语言(graphics language ，g 语言) 的软件开发环境，常用的有labview(laboratory virtual instrument engineering workbench)和hpvee等。labview是一种用图标代替文本行创建应用程序的图形化编程工具。由于labview提

10、供了大量常用于自动化测试测量领域的图形控件，这使得用户可以在数分钟之内完成一套完整的从仪器连接、数据采集到分析、显示和存储的自动化测试测量系统。 选用labview开发测试和测量应用程序的另一大决定因素是其开发速度。一般来说，它比使用其它编程语言快410 倍。并且labview易用易学，具有以下几个主要的优势（1） 提供了丰富的图形控件，采用图形化的编程方法； （2）采用数据流模型，实现了自动的多线程，从而充分利用处理器的处理能力； （3）提供了大量的驱动与专用工具，几乎能与任何接口的硬件轻松连接； （4）内建了600多个分析函数，用于数据分析和信号处理； （5） 通过dll 、activex

11、、.net或matlab脚本节点等技术，可以轻松实现labview与其它编程语言的混合编程； （6）提供了丰富的附加模块，用于扩展labview在不同领域中的应用，如实时模块、pda 模块、fpga模块等1.2 总体设计 本系统的软件部分旨在完成对肌电信号的采集、处理、分析等。因此，系统具备以下几个功能（1）把采集到的肌电数据进行放大和转换以便于对信号进行处理（2）实时采集肌电信号并进行滤波、保存等处理；（3）可随时调出保存的肌电数据进行分析研究。 根据系统的功能要求，本次设计采用模块化设计，这样便于调试和管理，也易于调用。总体自上而下分为系统时钟模块、信号采集模块、电刺激控制模块、波形回放模

12、块和波形打印模块。而信号采集模块由放大倍数调节模块、数据转换模块、滤波及数据处理模块和存储模块组成。电刺激控制模块由参数控制模块、波形预览下载模块、时间控制报警模块和电刺激波形显示调节模块。设计模块层次图如图1所示。图1 模块层次图计算机接收到数据后，首先进行原始数据的预处理，把数据进行转化和数字滤波处理，经过处理的数据由labview的实时趋势图控件显示出信号曲线。此外，还要对信号要进行进一步的分析处理，从而得到更便于用户接受的反馈信号。用户还可将信号数据存储、回放和打印。在电刺激控制部分，可把预先编制好的治疗处方波形下载到波形存储器中，并通过串口把处方选择、频率、幅值参数传给单片机，然后由

13、单片机完成参数控制的工作。在治疗上，前面板有时间报警控制，等治疗时间到有声光报警功能。前面板设计如图2所示。2.模块介绍2.1肌电信号采集与处理上位机向下位机发送采集命令后，开始接收采集数据，并调节放大倍数下传给单片机，使信号按比例放大。采集的数据经过数据转换、数字滤波和数据处理后显示其信号。该部分的信号主要由采集模块完成，采集模块由放大倍数调节模块、数据转换模块、滤波及数据处理模块和存储模块组成。这里主要介绍数据转换模块、滤波及数据处理模块和存储模块。2.1.1数据转换模块由上周的硬件设计可知：硬件系统选用的ad转换芯片是12位的，而51单片机的数据总线宽度为8位，所以我们必须用两个8位的数

14、据单元存储一个12位的数据，才能保证不浪费12位的分辫率，这样计算机接收到的数据，其实是两个数据表示一个采样点。labview提供8位、16位、32位的二进制数，要正确显示和分析信号，就必须将相应的两个数据转化成16位或32位的数据，本系统选用了16位的数据，转换后为d(i)=d(j)×16+d(j+1)d(i)为转换后16为数据，d(j)为上传的高8为数据，d(j+1)为上传的低4位数据。来自串口的数据是字符串类型的，因此转换之前必须先进行字符串数据到数值型数据的转换。根据ad输入电压-10-+10v，转换后数据0-4095的正比例线性关系得出关系式：a(i)=(（d(i)2n-1

15、）/2n-1)×vm=(（d(i)2048)/2048)×10a(i)为输入电压，d(i)为转换的数字量，vm为输入最大幅度。采样数据数据通过以上两式转化成了相应的电压值，电压值经显示图表显示就能较好的恢复成原始信号，让用户能够方便真实的看到原始信号的各种特征。数据转换程序框图如图3所示。图3 数据转换程序框图2.1.2滤波及数据处理模块上位机接收到数据后还要进行结果处理，为了进一步消除干扰，在上位机进行数字滤波处理。由于采集的肌电信号频率分布在10-1khz的范围，所以肌电信号采集到pc机后经过10-1khz的带通滤波和50hz陷波器处理。运动单元是肌肉活动的最小单位，也

16、是肌电信号产生的最直接因素，当人的肌肉紧张程度不同时，参与活动的运动单元数量也不相同，相应的肌电信号的幅度和信号频率范围也不一样。在肌电信号的时域分析中，将肌电信号看作是时间的函数，可计算信号的均值、幅度直方图等统计指标来反映信号振幅在时间域上的变化。常用指标有电活动水平、积分肌电值和均方根值。由于肌电信号振幅和肌张力呈一定力一电对应关系，故时域指标可实时反映肌电活动水平。经过滤波和积分肌电值的数据处理程序框图如图4所示。通过仿真实验得到的采集肌电信号和积分肌电值的波形如图5所示。通过积分肌电值更加直观的把肌肉紧张程度反馈给患者，当肌肉紧缩时，幅度变大，频率变高；当肌肉收缩时，幅度变小，频率变

17、低。患者可根据显示的肌电信号和积分肌电值来进行康复训练，或者可以和电刺激相结合，电刺激器由计算机控制，采集患者的肌电信号，然后按照预先确定的刺激参数控制电刺激器对患者的肌肉或其支配神经进行刺激，以产生肌肉收缩，带动关节按一定规律运动，起到康复训练的作用。图4 滤波与数据处理程序框图图5 采集的肌电信号与积分肌电值2.1.3数据存储模块数据存储模块是程序中重要的模块之一，它主要实现对采集的数据存储功能，以便于实现数据的重复使用、显示和打印，为人工处理提供方便。数据存储有文件存储和数据库存储等方式，本系统数据的存储采用文件存储。labview采集的数据可通过“写入电子表格文件”vi直接存为文本文件

18、、电子表格文件或者是工程计算语言matlab所支持的m文件，用户可以通过记事本、excel,matlab或labview读出这些文件。其存储模块的程序框图如图6所示。图6 存储模块的程序框图2.2电刺激波形的控制与下传该部分主要完成电刺激波形编制和下传，与单片机通信完成对处方波形，频率，幅度参数的控制以及治疗时间控制声光报警。整个电刺激程序框图主要包括：参数下传通信模块、波形预览和下载模块、时问控制报警模块和电刺激波形显示调节模块。这里对主要模块进行介绍。2.2.1参数下传通信模块通过visa configure serial port(串口初始化子v1)，visa read(读串口子vi)，

19、visa write(写串口子vi)，visa bytes at serial port(检测串口输入缓冲中字节数子vi)，visa close(关闭串口子vi)等库函数完成和下位机通信，处理处方选择，幅度，频率参数调节的下传，程序框图如图7所示。图7参数下传与通信程序框图此系统在一定程度上实现处方化，允许用户对治疗方案和治疗波形参数进行动态干预，根据具体的处方选择参数，实现各种处方的治疗。在进行电刺激前，使用者可根据需要选择刺激波形，频率和幅度，通过“visawrite”vi把命令和数据传递至单片机，单片机循环接收上位机传来的命令和数据，根据设置的值来调节存储器的波形选择的地址，样本点读取的

20、时间和幅度的控制。2.2.2波形预览和下载模块labview软件库中有很多波形发生函数，调用子函数很容易实现各种波形的编制。根据电刺激波形的要求，把每个波形一个周期用32个样本点来编制。本系统根据治疗电刺激和功能电刺激的要求，编制了5种刺激波形，处方类型如表1所示。使用者可选择处方类型的波形，再根据具体的处方和个人的需要来调节参数，达到治疗和训练的目的。例如在康复治疗时，选择治疗处方波形l间动疏波(me)，单相整流正弦波，频率固定为50hz，用户可直接根据治疗处方和需要，调节幅度和治疗时间。再如在康复训练时，选择处方波形5功能电刺激，其为脉宽在50500us可调的方波脉冲，频率以步进为1可调，

21、一般在低频范围内，通过控制时间来实现产生刺激脉冲的个数，如要实现脉宽200us，频率2hz，持续2秒，则是4个刺激脉冲。在训练时采集患处的肌电信号作为反馈患者的信息，然后按预先确定的刺激参数控制刺激器对患处肌肉或神经进行单次或短时间持续刺激，以产生肌肉收缩，带动关节按一定规律运动，起到康复训练的作用。表1处方类型表处方类型号处方类型名描述1间动疏波（mf）低频：单项整流，50hz2间动密波（df）低频：双向z整流，100hz32-5超刺激低频：2ms脉宽5ms间歇时间，约143hz4等幅正弦中频中频：1khz-4khz可调5功能电刺激低频：刺激脉宽50-500us，频率和时间可调在下载前可预览

22、波形，确认后再依次将一个周期的波形下传到波形存储器指定的单元中。在下载之前需要将电刺激波形数据转换成0255的下载存储器中的数据，da采用8位的双极性输出，根据输出-vref+vref和255的正比例线性对应关系，得到关系式：vout2=vref×(di128)/128其中vout2为da输出的模拟信号，vref为参考电压，di为数字信号，根据上式可得到di=（vout2/vref）×128128，其中（vout2/vref）1，这里编制的样本波形幅值为l，下传到存储器中，由vref电压控制输出幅度，从而得到相应波形的输出波形预览下载前面板和程序框图如图8和图9所示。图8波

23、形预览和下载前面板图9间动疏波下载程序框图2.2.3时间控制声光报警模块该模块主要完成时间设置从060分钟的控制，前面板有时间显示和时间复位，等设置的时间到，发出声光报警，并向下位机发出命令停止刺激。程序框图如图10所示。图10时间控制和报警程序框图2.3其它功能模块2.3.1系统时间模块该模块的功能是在系统前面板上实时显示日期和时间，用一个while循环来实现。获取日期时间（秒）把当前时间输出给获取日期时间字符串，再由获取日期时间字符串分别输出目期和时间。程序框图如图11所示。图11系统时间框图程序2.3.2波形回放和打印模块波形回放模块是通过使用“读取电子表格文件”vi读出设定路径下已存储

24、的数据文件，将其波形直接显示。打印模块的功能是对回放波形当前显示波形的打印输出。波形回放和打印前面板和波形回放程序框图如图12和图13所示。打印模块程序框图如图14所示。图12波形回放和打印前面板图13波形回放程序框图图14打印程序框图总结此次课程设计圆满结束了，对于我经过努力所得的结果，我感到非常满意。此次设计以分组的方式进行，我们做的是肌电信号采集和分析系统。由于平时大家都是学理论，没有过实际开发设计的经验，虽然上课的时候讲过很多关于labview的应用，但课后也从没真正地应用过这个软件，对它一点都不熟悉，拿到的时候都不知道怎么做。但通过对labview进行了几天的摸索学习，我终于把课程设

25、计做好了，并发现labview很实用。通过这次作业我基本学会了虚拟仪器设计的步聚和基本方法，提高了动手能力，增强了理论联系实际的能力，分组工作的方式给了我与同学合作的机会，提高了与人合作的意识与能力，同时也培养和发展了同学之间的感情，对此我感到非常欣慰。总之，在这次作业中我受益匪浅。 在这个作业中，我真正的感受到了虚拟仪器的优越性，“软件即仪器”，仅仅在labview上写入一些图形程序，就可以实现相应的仪器的功能，只需要借助一台计算机，一个labview软件，就可以实现很多仪器的功能，真的好方便。 此次作业的过程中，我遇到了不少的问题，但最终我终究还是战胜困难，取得了成功，当然其中很大一部分功劳都要归功于指导老师老师和同学们，要不是没有指导老师的指导和与同学们的相互交流学习，我的课程设计肯定不会如此轻松。参考文献1王磊,陶梅.精通labview8.0m.北京:电子工业出版社,2007. 2胡仁喜等.labview8.2.1虚拟仪器实例指导教程m.北京:机械工业出版社,2008. 3杨乐平,李海涛,杨磊,安雪滢.labview高级程序设计m.北京:清华大学出版社,2003. 4吴成东,孙秋野,盛科.labview虚拟仪器程序设计及应用m.北京:人民邮电出版社,2008. 葫捞帆擅拘八士连眉诊咳州擦贡昂膀晶摧娜陡现菜洱驶