【肌电信号】基于matlab GUI MUAP波形【含Matlab源码 736期】
一、简介
EMG 信号是很容易采集的,但是要想理解和解释所采集到的EMG信号是非常困难的。对EMG信号进行分析可以让我们具备理解肌肉中力的产生机制以及肌肉如何产生运动的能力。但如何提取这些信息取决于适当的信号获取办法。
1 EMG信号提取过程
信号采集,信号调节和信号处理。流程图如下所示
2 信号采集
肌肉收缩可以引起皮肤表面电流发生变化,此电流的大小可以通过放置在皮肤表面的电极来测量得到。电极与皮肤接触的区域被称为测量表面(detection surface)。利用表面电极测量得到的生理学数据被称为表面肌电流(sEMG)。sEMG是电极测量的来自于所有活跃运动单元的MUAPTs的总和。
随着肌肉输出力的增加,活跃运动单元的数量也随着增加。由于肌肉收缩时许多运动单元都是活跃的,因此,电极测量到的信号来自于许多不同的运动单元。MUAP的形状和幅值取决于肌纤维和电极的相对方向(和其他因素相比)。在一个MUAPT中,如果电极和肌纤维的相对位置保持不变,电极的属性保持不变以及肌肉组织的生物化学属性没有改变,则MUAP的的波形保持不变。
有作者在144篇文章里发现了352种采集表面肌电流信号时电极安放位置的描述。其中大部分的方法都是通用的,并且提到肌腹(肌肉中部较厚的区域),Motor point(轴突和肌纤维接触的地方),或者肌肉中点。
影响EMG稳定性的因素包括:
- 运动单元/肌腱
- 附近其他活跃肌肉
- 活跃肌肉纤维和测量点的距离
- 电极的滤波特性
- 神经区域(innervation zone)和电极的位置之间的关系
通常,双极点电极相对于肌纤维平行放置,并且两电极中心距离在20 mm以上的时候可以取得较好的结果。这种距离足够小可以避免信号的串扰,并且足够大选择众多的运动单元。
参考电极应该放置在中性组织上(i.e. 骨头突出的部分),以便于给电极放大器提供一个公用的参考。
3 信号调节
电极-放大器用来减小电极-皮肤的阻抗。
信号调节用来进一步改善EMG信号的质量。
信号调节的方法包括:
滤波以减小运动伪迹
电极运动伪迹可能来自于表面电极区域下的皮肤变形,或者电极电荷层的扰动。运动伪迹的功率密度小于20 Hz。因此,高通滤波器可用于改善信号的质量。这种滤波器的转折频率应该在10 Hz和20 Hz之间,但是不能大于20 Hz,否则会导致EMG信号能量的丢失。
用于运动伪逆的滤波技术包括:
- 8th order Chebychev high pass filters
- adaptive filtering based on orthogonal Meyer wavelets
EMG信号放大
要特别注意工频噪声的干扰。
power line interference
4 信号处理
EMG 信号用作控制信号:
- 多功能假肢
- 轮椅
- 抓握控制
- 虚拟键盘
- 基于姿势的控制接口
5 数据分割
EMG信号有两种状态:
- 瞬态: 肌肉从静止到自主收缩
- 稳态: 肌肉产生常力收缩
瞬态EMG信号比稳态EMG信号处理起来要更困难一些。
6 特征提取
时域特征
EMG幅值和活跃运动单元的数量以及他们的活跃水平有关系。
形成EMG幅值的方法有:
- 模拟修正和平滑(低通)
- 平均绝对值(MAV)处理
- 均方根(RMS)处理
对于高层次的收缩,可以采用高斯模型对EMG信号进行拟合;对较弱的肌肉收缩和处于肌肉疲劳状态的收缩可以采用LaPlacian进行拟合。
其他EMG时域特征包括:平均绝对值(MAVs)斜率,zero crossing,斜率符号改变和波形长度。
7 频域特征
时域-频域特征
专业用语
surface EMG( sEMG): 表面肌电信号
MU (motor unit): 运动单元
MUAPTs (motor unit action potential trains): 动作电位序列
motion artifact: 运动伪迹
二、部分源代码
function varargout = MUAPgui(varargin)
% MUAPGUI M-file for MUAPgui.fig
% MUAPGUI, by itself, creates a new MUAPGUI or raises the existing
% singleton*.
%
% H = MUAPGUI returns the handle to a new MUAPGUI or the handle to
% the existing singleton*.
%
% MUAPGUI('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in MUAPGUI.M with the given input arguments.
%
% MUAPGUI('Property','Value',...) creates a new MUAPGUI or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before MUAPgui_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to MUAPgui_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES% Copyright 2002-2003 The MathWorks, Inc.% Edit the above text to modify the response to help MUAPgui% Last Modified by GUIDE v2.5 27-Dec-2006 12:09:22% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @MUAPgui_OpeningFcn, ...'gui_OutputFcn', @MUAPgui_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);
if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});
endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
elsegui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT% --- Executes just before MUAPgui is made visible.
function MUAPgui_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to MUAPgui (see VARARGIN)% Choose default command line output for MUAPgui
handles.output = hObject;% Update handles structure
guidata(hObject, handles);
clear all;
% UIWAIT makes MUAPgui wait for user response (see UIRESUME)
% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.
function varargout = MUAPgui_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structure
varargout{1} = handles.output;% --- Executes on button press in Laplace.
function Laplace_Callback(hObject, eventdata, handles)
% hObject handle to Laplace (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)global muap;
lapmuap=4*muap(2,:)-muap(1,:)-muap(3,:)-muap(4,:)-muap(5,:);
axes(handles.axes1);
HH=plot(lapmuap);
set(HH,'LineWidth',2);
clear HH;% --- Executes on button press in Bipolar.
function Bipolar_Callback(hObject, eventdata, handles)
% hObject handle to Bipolar (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global muap;biomuap=muap(1,:)-muap(2,:);
axes(handles.axes2);
HH=plot(biomuap);
set(HH,'LineWidth',2);
clear HH;% --- Executes on button press in Monopolar.
function Monopolar_Callback(hObject, eventdata, handles)
% hObject handle to Monopolar (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)global muap;axes(handles.axes3);
HH=plot(muap(2,:));
set(HH,'LineWidth',2);
clear HH;function nf_Callback(hObject, eventdata, handles)
% hObject handle to nf (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of nf as text
% str2double(get(hObject,'String')) returns contents of nf as a double% --- Executes during object creation, after setting all properties.
function nf_CreateFcn(hObject, eventdata, handles)
% hObject handle to nf (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispcset(hObject,'BackgroundColor','white');
elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
endfunction skinfat_Callback(hObject, eventdata, handles)
% hObject handle to skinfat (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of skinfat as text
% str2double(get(hObject,'String')) returns contents of skinfat as a double% --- Executes during object creation, after setting all properties.
function skinfat_CreateFcn(hObject, eventdata, handles)
% hObject handle to skinfat (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispcset(hObject,'BackgroundColor','white');
elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
endfunction x_Callback(hObject, eventdata, handles)
% hObject handle to x (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of x as text
% str2double(get(hObject,'String')) returns contents of x as a double% --- Executes during object creation, after setting all properties.
function x_CreateFcn(hObject, eventdata, handles)
% hObject handle to x (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually ave a white background on Windows.
% See ISPC and COMPUTER.
if ispcset(hObject,'BackgroundColor','white');
elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end
三、运行结果
四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
[4]徐洁.基于小波分析的脉搏波信号处理[J].电子设计工程. 2013,21(11)
【肌电信号】基于matlab GUI MUAP波形【含Matlab源码 736期】相关推荐
- 【信号处理】基于Matlab GUI数字均衡器设计,源码附上
[信号处理]基于Matlab GUI数字均衡器设计,源码附上 数字均衡器是一种常用于音频信号处理的技术,可以有效地补偿声音中的频率失真,提高音质.在这篇文章中,我们将介绍如何使用Matlab GUI进 ...
- 测试基于STM32的ADIS16405评估板,并在匿名四轴上实时显示3轴陀螺仪+3轴加速度计+3轴磁力计 波形。本文最后还公开一个基于C# 的串口波形显示客户端源码。
标题:测试基于STM32的ADIS16405评估板,并在匿名四轴上实时显示3轴陀螺仪+3轴加速度计+3轴磁力计 波形.本文最好还公开一个基于C# 的串口波形显示客户端源码. 这里附上该评估板的原理图和 ...
- 【过程控制】基于matlab GUI串级控制【含Matlab源码 2385期】
⛄一.简介 1 交互系统整体外观展示 该系统基于 MATLAB GUI 与 Simulink 联合开发, 实现对串级控制系统特性的研究与分析.学生在学习串级控 制有关知识后可使用该应用进行知识回顾.界 ...
- matlab svm 语音识别,【情感识别】基于matlab支持向量机(SVM)的语音情感识别【含Matlab源码 543期】...
一.简介 支持向量机(Support Vector Machine)是Cortes和Vapnik于1995年首先提出的,它在解决小样本.非线性及高维模式识别中表现出许多特有的优势,并能够推广应用到函数 ...
- 【Matlab肌电信号】肌电信号处理【含GUI源码 966期】
一.代码运行视频(哔哩哔哩) [Matlab肌电信号]肌电信号处理[含GUI源码 966期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 包子阳,余继周,杨 ...
- 【Matlab语音分析】语音信号分析【含GUI源码 1718期】
一.代码运行视频(哔哩哔哩) [Matlab语音分析]语音信号分析[含GUI源码 1718期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊,郑铁 ...
- 【Matlab语音加密】语音信号加密解密(带面板)【含GUI源码 181期】
一.代码运行视频(哔哩哔哩) [Matlab语音加密]语音信号加密解密(带面板)[含GUI源码 181期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆 ...
- 【Matlab语音处理】声音信号频谱分析仪【含GUI源码 325期】
一.代码运行视频(哔哩哔哩) [Matlab语音处理]声音信号频谱分析仪[含GUI源码 325期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊, ...
- 【Matlab心音信号】EMD心音信号特征提取【含GUI源码 1735期】
一.代码运行视频(哔哩哔哩) [Matlab心音信号]EMD心音信号特征提取[含GUI源码 1735期] 二.matlab版本及参考文献 1 matlab版本 2014a *2 参考文献 [1] 沈再 ...
- 【Matlab语音处理】汉宁窗FIR陷波滤波器语音信号加噪去噪【含GUI源码 1711期】
一.代码运行视频(哔哩哔哩) [Matlab语音处理]汉宁窗FIR陷波滤波器语音信号加噪去噪[含GUI源码 1711期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 ...
最新文章
- 基于存储过程的可扩展性数据访问框架
- 微软亚研提出VL-BERT:通用的视觉-语言预训练模型
- ffmpeg推流时与服务器断开后的自动重连功能的实现
- SAP系统和微信集成的系列教程之十:如何在SAP C4C系统里直接回复消息给微信用户
- opencv 最大连通域_opencv 查找连通区域 最大面积实例
- AD19 add pins to nets错误_NGINX 502 Bad Gateway错误疑难解答
- php导入img,PHP 转PDF(Img)
- Java基础学习总结(66)——配置管理库typesafe.config教程
- Select和epoll的区别
- 职场人,职场事,八卦一下职场那些有意思的公众号
- 使用STVP解除STM32flash读出保护
- 湖北一考生将高考数学题上传小猿搜题事件网络舆情综合编报
- 梦想在三十岁起航!__来自黑马程序员69期安卓班的学员
- Java调试--排查类工具
- 集群和分布式的区别,软件架构的演化过程,Apache Dubbo详解
- 基于simulink采用 QSHB 和 HBPS 算法的混合 MIMO 波束成形仿真
- 香蕉皮的妙用(想歪的去面壁)(不管有没有用,先记下再说)
- 输入年月日,输出这一天是这一年的第几天
- 安装MATLAB工具箱Robotics Toolbox 9.10.0 版本
- 数通运营商方向常见面试问题(第三部分)