How to design multirate filter system for EMG signals

Question

Dulanjana Perera 2020-9-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/587210-how-to-design-multirate-filter-system-for-emg-signals

编辑： Dulanjana Perera 2020-9-7

Marked_signal.jpg

I am using Delsys Bagnoli EMG System to acquire sEMG from forearm muscles. When I use Delsys EMGworks application, output signals are clean as expected with minor interferences. However, when I use the following Matlab script to acquire the signal unusual signal pattern occured (spikes).

the image of frequency domain of the Data (Red color frequencies contains EMG information) acquired from Matlab script is attached herewith.

I want to know why this happens and how to remove thoses spikes

the Matlab script:

%% EMG initialize
raw = zeros(1,6);
% EMG setup
emg = daq.createSession('ni');
ch1 = addAnalogInputChannel(emg, 'Dev1', 0, 'Voltage');
ch2 = addAnalogInputChannel(emg, 'Dev1', 1, 'Voltage');
ch3 = addAnalogInputChannel(emg, 'Dev1', 2, 'Voltage');
ch4 = addAnalogInputChannel(emg, 'Dev1', 3, 'Voltage');
ch5 = addAnalogInputChannel(emg, 'Dev1', 4, 'Voltage');
ch6 = addAnalogInputChannel(emg, 'Dev1', 5, 'Voltage');
emg.Rate = 4000;                           % sampling rate
lh = emg.addlistener('DataAvailable',@stopWhenExceed);
emg.NotifyWhenDataAvailableExceeds = 400; % data per one feature set (total of all channels)
emg.IsContinuous = true;
emg.startBackground();
while 1
   
    % check the existing time.mat file
    time(end + 1 : end + 1, 1) = str2double(datestr(now,'SS.FFF'));
    pause(0.01); % this pause should be there to acquire EMG data
    if (getGlobalSx==0)           % if new set of data is available
        data = getGlobalx;          % aquire data
        setGlobalSx(1);
    end
    raw(end +1:end+400, : ) = data;
    
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Dulanjana Perera 2020-9-7

在 MATLAB Online 中打开

I realized that the pause we provide to the code to acquire EMG data should be equal to the

emg.NotifyWhenDataAvailableExceeds / emg.Rate

which means 400/ 4000 = 0.1 not 0.01

Otherwise this multirate signal sampling happens.

while 1
   
    % check the existing time.mat file
    time(end + 1 : end + 1, 1) = str2double(datestr(now,'SS.FFF'));
    
    % ---------------------This value should be 0.1.-------------------------------
    pause(0.1); % this pause should be there to acquire EMG data
    if (getGlobalSx==0)           % if new set of data is available
        data = getGlobalx;          % aquire data
        setGlobalSx(1);
    end
    raw(end +1:end+400, : ) = data;
    
end

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Dulanjana Perera 2020-9-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/587210-how-to-design-multirate-filter-system-for-emg-signals#answer_490930

编辑：Dulanjana Perera 2020-9-7

在 MATLAB Online 中打开

I realized that the pause we provide to the code to acquire EMG data should be equal to the,

emg.NotifyWhenDataAvailableExceeds / emg.Rate

which means 400/ 4000 = 0.1 not 0.01

Otherwise this multirate signal sampling happens.

while 1
   
    % check the existing time.mat file
    time(end + 1 : end + 1, 1) = str2double(datestr(now,'SS.FFF'));
    
    % ---------------------This value should be 0.1.-------------------------------
    pause(0.1); % this pause should be there to acquire EMG data
    if (getGlobalSx==0)           % if new set of data is available
        data = getGlobalx;          % aquire data
        setGlobalSx(1);
    end
    raw(end +1:end+400, : ) = data;
    
end