Exercise L6:
% 1 a plot mV tegen tijd
load ('ExerL6_1.mat')
N = length (emg);
k = [0:N-1];
dt = 1/fs;
t = k*dt;
plot (t, emg);
xlabel('Tijd [s]')
ylabel ('Signaal [mV]')
title ('Signaal tegen tijd')
% 1 b Welke storende signalen (bewegingsartefacten, crosstalk,
besmetting, elektromagnetische storing en intrisieke ruis) zijn
aanwezig in het signaal gebaseerd op het plot.
Bewegingsartefacten zijn te herkennen door lage frequentie fluctuaties in
het signaal. Interne ruis is altijd aanwezig in een EMG signaal.
% 1 c Onderzoek de aanwezigheid van een power hum in het power spectrum.
Bepaal eerst welk deel van het signaal gebruikt moet worden voor het power
spectrum.
[Pxx, F] = pwelch (emg-(mean(emg)),[],[],[],fs);
figure(2)
plot (F, Pxx)
De meeste waarden zitten tussen de 25 en 250 Hz. De power hum zit op
50 Hz. Hier moet rekening mee gehouden worden, als de power hum
verwijdert moet worden.
% 1 d Van welke van de bovengenoemde storende signalen is het niet
mogelijk om te bepalen of ze in het signaal zitten (omdat je niet
genoeg informatie hebt)
Het is niet duidelijk of cross talk aanwezig is. Omdat je niet weet
welke spieren gemeten zijn.
% 1 e Verwijder de bewegingsartefacten van het signaal en plot het
resulterende signaal als functie van tijd
[B,A] = butter (3,20/(fs/2),'high');
emg_filt = filtfilt (B,A,emg);
plot (t, emg_filt);
% 1 f verwijder de power hum van het signaal
[D,C] = butter (4,[49.5, 50.5]/(fs/2),'stop');
emg_hum_filt = filtfilt (D,C,emg_filt);
figure(4)
plot (t, emg_hum_filt);
[Pxx, F] = pwelch (emg_hum_filt-(mean(emg_hum_filt)),[],[],[],fs);
figure(5)
plot (F,Pxx)
% 1 g Wat zou je adviseren om storende signalen te ontwijken of te
verzwakken in de toekomst
Door eerst een high pass filter om de bewegingsartefacten te
verwijderen en daarna een notch bandstopfilter om de power hum te
verwijderen. De meeste signalen zitten tussen de 6 en 500 Hz en de
meest interessante tussen de 25 en 250 Hz
% 1 a plot mV tegen tijd
load ('ExerL6_1.mat')
N = length (emg);
k = [0:N-1];
dt = 1/fs;
t = k*dt;
plot (t, emg);
xlabel('Tijd [s]')
ylabel ('Signaal [mV]')
title ('Signaal tegen tijd')
% 1 b Welke storende signalen (bewegingsartefacten, crosstalk,
besmetting, elektromagnetische storing en intrisieke ruis) zijn
aanwezig in het signaal gebaseerd op het plot.
Bewegingsartefacten zijn te herkennen door lage frequentie fluctuaties in
het signaal. Interne ruis is altijd aanwezig in een EMG signaal.
% 1 c Onderzoek de aanwezigheid van een power hum in het power spectrum.
Bepaal eerst welk deel van het signaal gebruikt moet worden voor het power
spectrum.
[Pxx, F] = pwelch (emg-(mean(emg)),[],[],[],fs);
figure(2)
plot (F, Pxx)
De meeste waarden zitten tussen de 25 en 250 Hz. De power hum zit op
50 Hz. Hier moet rekening mee gehouden worden, als de power hum
verwijdert moet worden.
% 1 d Van welke van de bovengenoemde storende signalen is het niet
mogelijk om te bepalen of ze in het signaal zitten (omdat je niet
genoeg informatie hebt)
Het is niet duidelijk of cross talk aanwezig is. Omdat je niet weet
welke spieren gemeten zijn.
% 1 e Verwijder de bewegingsartefacten van het signaal en plot het
resulterende signaal als functie van tijd
[B,A] = butter (3,20/(fs/2),'high');
emg_filt = filtfilt (B,A,emg);
plot (t, emg_filt);
% 1 f verwijder de power hum van het signaal
[D,C] = butter (4,[49.5, 50.5]/(fs/2),'stop');
emg_hum_filt = filtfilt (D,C,emg_filt);
figure(4)
plot (t, emg_hum_filt);
[Pxx, F] = pwelch (emg_hum_filt-(mean(emg_hum_filt)),[],[],[],fs);
figure(5)
plot (F,Pxx)
% 1 g Wat zou je adviseren om storende signalen te ontwijken of te
verzwakken in de toekomst
Door eerst een high pass filter om de bewegingsartefacten te
verwijderen en daarna een notch bandstopfilter om de power hum te
verwijderen. De meeste signalen zitten tussen de 6 en 500 Hz en de
meest interessante tussen de 25 en 250 Hz
