肌电信号基于matlab GUI脉搏信号处理系统含Matlab源码 1062期

Posted 2021-07-20 紫极神光（Q1564658423）

一、简介

EMG，应用电子学仪器记录肌肉静止或收缩时的电活动，及应用电刺激检查神经、肌肉兴奋及传导功能的方法。英文简称EMG。通过此检查可以确定周围神经、神经元、神经肌肉接头及肌肉本身的功能状态。
通过测定运动单位电位的时限、波幅，安静情况下有无自发的电活动，以及肌肉大力收缩的波型及波幅，可区别神经原性损害和肌原性损害，诊断脊髓前角急、慢性损害（如脊髓前灰质炎、运动神经元疾病），神经根及周围神经病变（例如肌电图检查可以协助确定神经损伤的部位、程度、范围和预后）。另外对神经嵌压性病变、神经炎、遗传代谢障碍神经病、各种肌肉病也有诊断价值。此外，肌电图还用于在各种疾病的治疗过程中追踪疾病的恢复过程及疗效。
利用计算机技术，可作肌电图的自动分析，如解析肌电图、单纤维肌电图以及巨肌电图等，提高诊断的阳性率。
肌电图检查多用针电极及应用电刺激技术，检查过程中有一定的痛苦及损伤 ，因此除非必要 ，不可滥用此项检查。另外，检查时要求肌肉能完全放松或作不同程度的用力，因而要求受检者充分合作。对于某些检查，检查前要停药，如新斯地明类药物应于检查前16小时停用
记录肌肉动作电位的曲线（电描记图）称为肌电图。缩写为EMG。实际使用的描记方法有两种：一种是表面导出法，即把电极贴附在皮肤上导出电位的方法；另一种是针电极法，即把针电极刺入肌肉导出局部电位的方法。用后一种方法能分别记录肌肉每次的动作电位，而根据从每秒数次到二、三十次的肌肉动作电位情况，发现频率的异常。应用肌电图还可以诊断运动机能失常的原因。平常所用的针电极称为同心电极，它是把细针状电极穿过注射针的中心，两者绝缘固定制成的。

二、源代码

function varargout = linbingwen(varargin)
% LINBINGWEN M-file for linbingwen.fig
%      LINBINGWEN, by itself, creates a new LINBINGWEN or raises the existing
%      singleton*.
%
%      H = LINBINGWEN returns the handle to a new LINBINGWEN or the handle to
%      the existing singleton*.
%
%      LINBINGWEN('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in LINBINGWEN.M with the given input arguments.
%
%      LINBINGWEN('Property','Value',...) creates a new LINBINGWEN or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before linbingwen_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to linbingwen_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

% Edit the above text to modify the response to help linbingwen

% Last Modified by GUIDE v2.5 27-Jun-2021 16:54:46

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @linbingwen_OpeningFcn, ...
                   'gui_OutputFcn',  @linbingwen_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before linbingwen is made visible.
function linbingwen_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 linbingwen (see VARARGIN)

% Choose default command line output for linbingwen
handles.output = hObject;
t = timer('tag', 'timer1', 'BusyMode', 'queue', 'ExecutionMode', ...
    'fixedrate', 'Period', 1, 'timerfcn', {@t_update, handles});
start(t);
% Update handles structure


% UIWAIT makes linbingwen wait for user response (see UIRESUME)
% uiwait(handles.figure1);
movegui( handles.figure1, 'center' );
set(handles.pushbutton1,'enable','off');%下面几条指令是不启动其余按钮功能
set(handles.pushbutton2,'enable','off');
set(handles.pushbutton7,'enable','off');
set(handles.pushbutton8,'enable','off');
set(handles.pushbutton14,'enable','off');
set(handles.pushbutton9,'enable','off');
set(handles.pushbutton10,'enable','off');
set(handles.pushbutton11,'enable','off');
set(handles.pushbutton17,'enable','off');
set(handles.pushbutton18,'enable','off');
set(handles.pushbutton19,'enable','off');
set(handles.pushbutton20,'enable','off');
set(handles.pushbutton15,'enable','off');
set(handles.pushbutton16,'enable','off');
set(handles.pushbutton21,'enable','off');
set(handles.pushbutton22,'enable','off');
set(handles.pushbutton23,'enable','off');
set(handles.pushbutton24,'enable','off');
axes(handles.axes1);    
P = plot(5*ones(1000,1));grid on;
axis([0,1,-2,4])
axes(handles.axes2);    
P = plot(5*ones(1000,1));grid on;
axis([0,1,-2,4])
guidata(hObject, handles);
% --- Outputs from this function are returned to the command line.
function t_update(obj,eventdata,handles)
if ~isfield(handles, 'figure1')
    return
end
st1 = datestr(now, 'yyyy-mm-dd');
st2 = datestr(now, 'HH:MM:SS');
set(handles.activex5, 'AlphaNumeric', st2)
set(handles.activex4, 'AlphaNumeric', st1)
function varargout = linbingwen_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 during object creation, after setting all properties.
function axes1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to axes1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: place code in OpeningFcn to populate axes1
set(hObject,'xTick',[]);%设置窗口axes3的外观属性，XY轴无标签。
set(hObject,'ytick',[]);
set(hObject,'box','on');%打开窗口

% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(handles.pushbutton2,'enable','on');
set(handles.pushbutton5,'enable','off');
set(handles.pushbutton6,'enable','off');
set(handles.pushbutton10,'enable','off');
set(handles.pushbutton11,'enable','off');
set(handles.pushbutton15,'enable','on');
set(handles.pushbutton16,'enable','on');
set(handles.pushbutton17,'enable','on');
set(handles.pushbutton18,'enable','on');
set(handles.pushbutton19,'enable','on');
set(handles.pushbutton20,'enable','on');
set(handles.pushbutton21,'enable','on');
set(handles.pushbutton22,'enable','on');
global str
global t
global maibo
global ymin
global ymax
global xmin
global xmax
global pauset
pauset=0.5;
data=load(str);
t=data(:,1);maibo=data(:,2);
ymin=min(data(:,2))-0.5;
ymax=max(data(:,2))+0.5;
%%
global s21
global dt
yy=data;
x=yy(:,2)';
dt=yy(:,1)';
fs=800;
%%去除基线漂移
k = .7; % cut-off value
fc=0.3/fs;
alpha = (1-k*cos(2*pi*fc)-sqrt(2*k*(1-cos(2*pi*fc))-k^2*sin(2*pi*fc)^2))/(1-k);
y = zeros(size(x));
for i = 1:size(x,1)
   y(1,:) = filtfilt(1-alpha,[1 -alpha],x(1,:));
end 
x1=x-y;
%%butter带通滤波
Wp=[0.9 50]/400;Ws=[0.3 140]/400;
[n,Wn] = buttord(Wp,Ws,3,10);
[b,a] = butter(n,Wn);
%%cheby去除50hz工频干扰
Wp1=[30 65]/400;Ws1=[45 55]/400;
rp=3;rs=60;
[n1,Wn1] = cheb2ord(Wp1,Ws1,rp,rs);
[b1,a1] = cheby2(n1,rs,Wn1,'stop');
%%信号之间的对比
x2=filtfilt(b,a,x1);
s21=filtfilt(b1,a1,x2);
axes(handles.axes1);
plot(t,maibo);
title('原始脉搏信号');
xlabel('时间/s');
ylabel('幅值/mv');
grid 
set(gca,'position',[10 27.5 120 20]);
xmin=0;
xmax=max(t)/10;
axis([xmin,xmax,ymin,ymax])
axes(handles.axes2);
plot(dt,s21);grid;
title('滤波去噪后脉搏信号');
xlabel('时间/s');
ylabel('幅值/mv');
set(gca,'position',[10 3 120 20]);
axis([xmin,xmax,ymin,ymax])
%%%
global flg1
flg1=1;
while 1
    if (flg1==0)
        break;
    end
    xmin=xmin+0.1;
    xmax=xmax+0.1;
 axes(handles.axes1);
plot(t,maibo);
title('原始脉搏信号');
xlabel('时间/s');
ylabel('幅值/mv');
grid 
set(gca,'position',[10 27.5 120 20]);
    axis([xmin,xmax,ymin,ymax]); %移动坐标系
    axes(handles.axes2);
plot(dt,s21);grid;
title('滤波去噪后脉搏信号');
xlabel('时间/s');
ylabel('幅值/mv');
set(gca,'position',[10 3 120 20]);
axis([xmin,xmax,ymin,ymax])
    pause(pauset);
end
guidata(hObject,handles);

% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(handles.pushbutton5,'enable','on');
set(handles.pushbutton6,'enable','on');
set(handles.pushbutton10,'enable','on');
set(handles.pushbutton11,'enable','on');
set(handles.pushbutton7,'enable','on');
set(handles.pushbutton17,'enable','off');
set(handles.pushbutton18,'enable','off');
set(handles.pushbutton19,'enable','off');
set(handles.pushbutton20,'enable','off');
set(handles.pushbutton15,'enable','off');
set(handles.pushbutton16,'enable','off');
set(handles.pushbutton21,'enable','off');
set(handles.pushbutton22,'enable','off');
global flg1
flg1=0;

三、运行结果

四、备注

版本：2014a