Obr. 8: Vzor zobrazení okamžité a průměrné aktivity svalu v klidu, sedu a stoji (Zdroj: Vlastní)
%Vyčistí pracovní prostor
clear;
clc;
%Načte výseky
load Vysek_Pr4_Klid.mat;
Klid = vysek;
load Vysek_Pr4_EliskaSed.mat;
Sed = vysek;
load Vysek_Pr4_EliskaStoj.mat;
Stoj = vysek;
clear vysek;
%Převedení typů dat pro další zpracování
ObalkaKlid(:,1) = Klid.LTrapDesc;
ObalkaKlid(:,2) = Klid.PTrapDesc;
ObalkaKlid(:,3) = Klid.LErSpinL;
ObalkaKlid(:,4) = Klid.PErSpinL;
ObalkaKlid(:,5) = Klid.LOblAbdE;
ObalkaKlid(:,6) = Klid.POblAbdE;
ObalkaKlid(:,7) = Klid.LExtCarpR;
ObalkaKlid(:,8) = Klid.PExtCarpR;
ObalkaSed(:,1) = Sed.LTrapDesc;
ObalkaSed(:,2) = Sed.PTrapDesc;
ObalkaSed(:,3) = Sed.LErSpinL;
ObalkaSed(:,4) = Sed.PErSpinL;
ObalkaSed(:,5) = Sed.LOblAbdE;
ObalkaSed(:,6) = Sed.POblAbdE;
ObalkaSed(:,7) = Sed.LExtCarpR;
ObalkaSed(:,8) = Sed.PExtCarpR;
ObalkaStoj(:,1) = Stoj.LTrapDesc;
ObalkaStoj(:,2) = Stoj.PTrapDesc;
ObalkaStoj(:,3) = Stoj.LErSpinL;
ObalkaStoj(:,4) = Stoj.PErSpinL;
ObalkaStoj(:,5) = Stoj.LOblAbdE;
ObalkaStoj(:,6) = Stoj.POblAbdE;
ObalkaStoj(:,7) = Stoj.LExtCarpR;
ObalkaStoj(:,8) = Stoj.PExtCarpR;
%Výpočet vztažných hodnot pro normalizaci - průměrné klidové aktivity
KlidovaAktivita = zeros(1,8);
KlidovaAktivita(:,1)=mean(Klid.LTrapDesc);
KlidovaAktivita(:,2)=mean(Klid.PTrapDesc);
KlidovaAktivita(:,3)=mean(Klid.LErSpinL);
KlidovaAktivita(:,4)=mean(Klid.PErSpinL);
KlidovaAktivita(:,5)=mean(Klid.LOblAbdE);
KlidovaAktivita(:,6)=mean(Klid.POblAbdE);
KlidovaAktivita(:,7)=mean(Klid.LExtCarpR);
KlidovaAktivita(:,8)=mean(Klid.PExtCarpR);
%Normalizace klidové aktivity
NormKlid = zeros(length(Klid.Time),8);
%Vydělí signál průměrnou klidovou aktivitou
for i=1:8
NormKlid(:,i) = (ObalkaKlid(:,i)./KlidovaAktivita(1,i));
end
%Výpočet průměrných hodnot klidu po normalizaci ke průměrné klidové aktivitě
PrumNormKlid = zeros(1,8);
%Vypočte průměry klidu
for i=1:8
PrumNormKlid(:,i)=mean(NormKlid(:,i));
end
%Normalizace sedu ke klidové aktivitě:
NormSed = zeros(length(Sed.Time),8);
%Vydělí signál průměrnou klidovou aktivitou
for i=1:8
NormSed(:,i) = (ObalkaSed(:,i)./KlidovaAktivita(1,i));
end
%Výpočet průměrných hodnot sedu po normalizaci ke průměrné klidové aktivitě
PrumNormSed = zeros(1,8);
PrumNormSed(:,1)=mean(NormSed(:,1));
PrumNormSed(:,2)=mean(NormSed(:,2));
PrumNormSed(:,3)=mean(NormSed(:,3));
PrumNormSed(:,4)=mean(NormSed(:,4));
PrumNormSed(:,5)=mean(NormSed(:,5));
PrumNormSed(:,6)=mean(NormSed(:,6));
PrumNormSed(:,7)=mean(NormSed(:,7));
PrumNormSed(:,8)=mean(NormSed(:,8));
%Normalizace stoje ke klidové aktivitě:
NormStoj = zeros(length(Stoj.Time),8);
%Vydělí signál průměrnou klidovou aktivitou
for i=1:8
NormStoj(:,i) = (ObalkaStoj(:,i)./KlidovaAktivita(1,i));
end
%Výpočet průměrných hodnot stoje po normalizaci ke průměrné klidové aktivitě
PrumNormStoj = zeros(1,8);
%Vypočte průměry stoje
for i=1:8
PrumNormStoj(:,i)=mean(NormStoj(:,i));
end
%Příprava zápisu průměrných hodnot do souboru
Prumery = zeros(3,8);
Prumery(1,:)=PrumNormKlid;
Prumery(2,:)=PrumNormSed;
Prumery(3,:)=PrumNormStoj
%Zápis průměrných hodnot do souboru
save Pr4_Eliska_Prumery.mat Prumery
%Zobrazení přes sebe
TPRNKlid=timetable(NormKlid(:,8),'SampleRate',1.5000e+03);
TPRNSed=timetable(NormSed(:,8),'SampleRate',1.5000e+03);
TPRNStoj=timetable(NormStoj(:,8),'SampleRate',1.5000e+03);
% Synchronize timetables
PravaRuka = synchronize(TPRNKlid,TPRNSed,TPRNStoj,'regular',...
'fillwithmissing','SampleRate',1500);
PravaRuka.Properties.VariableNames=["NormKlid","NormSed","NormStoj"];
figure;
plot(PravaRuka.Time,PravaRuka.NormKlid,'Color','b')
yline(PrumNormKlid(:,8),'Color','b')
title("Pravý extensor carpi radialis")
ylabel("Násobky průměrné klidové aktivity")
xlabel("Čas [s]")
hold on
plot(PravaRuka.Time,PravaRuka.NormSed,'Color','r')
yline(PrumNormSed(:,8),'Color','r')
plot(PravaRuka.Time,PravaRuka.NormStoj,'Color','g')
yline(PrumNormStoj(:,8),'Color','g')
legend(["NormKlid","PrumNormKlid","NormSed","PrumNormSed","NormStoj","PrumNormStoj"])
hold off
Obr. 8: Vzor zobrazení okamžité a průměrné aktivity svalu v klidu, sedu a stoji (Zdroj: Vlastní)
