肌肉激活度Matlab代码及数据.rar

【博客地址：https://blog.csdn.net/u011389706/article/details/97614699】利用Matlab求解肌肉激活度。包含源代码，sEMG数据，最大自主收缩（MVC）时的sEMG数据，运行后可以直接得出肌肉激活度的曲线。

function [ output_args1 ] = emg_activation（Fs EMG_raw EMG_MVC ）
% calculate activation from raw EMG
% EMG_raw: raw EMG signal
% EMG_MVC: MVC EMG signal
% the two vector can have different length 

% inital processing get e
EMG_raw = abs（EMG_raw）;
EMG_MVC = abs（EMG_MVC）;
[bb1ba1] = butter（4 30/（Fs/2） ‘high‘）;   %四阶butter高通滤波器，截止频率 30Hz，去除运动伪迹
preemg1 = filtfilt（bb1ba1 EMG_raw）;
preemg_mvc = filtfilt（bb1ba1 EMG_MVC）;
preemg1 = rot90（rot90（preemg1））;
preemg_mvc = rot90（rot90（preemg_mvc））;
preemg1 = filtfilt（bb1ba1 preemg1）;
preemg_mvc = filtfilt（bb1ba1 preemg_mvc）;

subplot（211）;
plot（preemg1）;   %%去除运动伪迹后的肌电信号

preemg1 = abs（rot90（rot90（preemg1）））;      %全波整流
preemg_mvc = abs（rot90（rot90（preemg_mvc）））;

hold on
plot（preemg1）;    %%全波整流后的信号
xlim（[0 601

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1901  2019-08-04 15:10  emg_activation.m

     文件        631  2019-08-04 15:12  main.m

     文件      26899  2019-07-28 14:05  move_sEMG.mat

     文件     108997  2019-08-04 15:02  MVC_sEMG.mat

     文件        337  2019-08-04 10:51  readme.txt

----------- ---------  ---------- -----  ----

               138765                    5