Scribd
Upload
28 views3 pages

All All 'Emg - Healthym - Mat'

This document contains code to analyze healthy and unhealthy EMG and ECG signals. It loads EMG/ECG data files, filters the signals, plots the original and filtered signals, and calculates various statistical metrics like standard deviation, mean, skewness, kurtosis, variance and median for both the original and filtered signals. This is done for healthy and unhealthy EMG and ECG samples to help characterize and distinguish between normal and abnormal signals.

Uploaded by

Tushar Sankaran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
28 views3 pages

All All 'Emg - Healthym - Mat'

This document contains code to analyze healthy and unhealthy EMG and ECG signals. It loads EMG/ECG data files, filters the signals, plots the original and filtered signals, and calculates various statistical metrics like standard deviation, mean, skewness, kurtosis, variance and median for both the original and filtered signals. This is done for healthy and unhealthy EMG and ECG samples to help characterize and distinguish between normal and abnormal signals.

Uploaded by

Tushar Sankaran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

EMG HEALTHY

clc;
clear all;
close all;
A=load('emg_healthym.mat');
B=A.val;
C=B/10000;
Fs=4000;
T=1/Fs;
D=C(3*Fs:10*Fs);
t1=3:T:10;
subplot(2,1,1);
plot(t1,D);
title('Healthy emg(original)');
Wn=0.05;
[b,a]=butter(6,Wn);
E=filter(b,a,D);
subplot(2,1,2);
plot(t1,E);
title('Filtered healthy EMG');
p1=std(D)
q1=std(E)
p2=mean(D)
q2=mean(E)
p3=skewness(D)
q3=skewness(E)
p4=kurtosis(D)
q4=kurtosis(E)
p5=var(D)
q5=var(E)
p6=median(D)
q6=median(E)

EMG UNHEALTHY
clc;
clear all;
close all;
A=load('emg_myopathym.mat');
B=A.val;
C=B/10000;
Fs=4000;
T=1/Fs;
D=C(3*Fs:10*Fs);
t1=3:T:10;
subplot(2,1,1);
plot(t1,D);
title('unhealthy emg(original)');
Wn=0.05;
[b,a]=butter(6,Wn);
E=filter(b,a,D);
subplot(2,1,2);
plot(t1,E);
title('Filtered unhealthy EMG');
p1=std(D)
q1=std(E)
p2=mean(D)
q2=mean(E)
p3=skewness(D)
q3=skewness(E)
p4=kurtosis(D)
q4=kurtosis(E)
p5=var(D)
q5=var(E)
p6=median(D)
q6=median(E)

HEALTHY ECG

clc;
clear all;
close all;
A=load('s20011m.mat');
B=A.val;
C=B/200;
Fs=250;
T=1/Fs;
D=C(5*Fs:10*Fs);
t1=5:T:10;
subplot(2,1,1);
plot(t1,D);
title('Healthy ecg(original)');
Wn=0.8;
[b,a]=butter(6,Wn);
E=filter(b,a,D);
subplot(2,1,2);
plot(t1,E);
title('Filtered healthy ECG');
p1=std(D)
q1=std(E)
p2=mean(D)
q2=mean(E)
p3=skewness(D)
q3=skewness(E)
p4=kurtosis(D)
q4=kurtosis(E)
p5=var(D)
q5=var(E)
p6=median(D)
q6=median(E)

UNHEALTHY ECG
clc;
clear all;
close all;
A=load('chf01m.mat');
B=A.val;
C=B/200;
Fs=250;
T=1/Fs;
D=C(5*Fs:10*Fs);
t1=5:T:10;
subplot(2,1,1);
plot(t1,D);
title('unhealthy ecg(original)');
Wn=0.8;
[b,a]=butter(6,Wn);
E=filter(b,a,D);
subplot(2,1,2);
plot(t1,E);
title('Filtered unhealthy ECG');
p1=std(D)
q1=std(E)
p2=mean(D)
q2=mean(E)
p3=skewness(D)
q3=skewness(E)
p4=kurtosis(D)
q4=kurtosis(E)
p5=var(D)
q5=var(E)
p6=median(D)
q6=median(E)

You might also like