How to design multirate filter system for EMG signals

1 次查看（过去 30 天）
Dulanjana Perera
Dulanjana Perera 2020-9-1
编辑： Dulanjana Perera 2020-9-7
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 个评论
Dulanjana Perera
Dulanjana Perera 2020-9-7
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

请先登录，再进行评论。

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

采纳的回答

Dulanjana Perera
Dulanjana Perera 2020-9-7
编辑：Dulanjana Perera 2020-9-7
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
Since this mismatch, a spike occurred at every 10th (0.1 / 0.01) frequency in the frequency domain.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Data Acquisition Toolbox Supported Hardware 的更多信息

标签

产品


版本

R2013a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by