Messung und Analyse myoelektrischer Signale - Communications ...

hold on;
%MPF Mean Power Frequency berechnen
g=2*pw;
sg=sum(g);
gn=g/sg;
mpf=sum(f.*gn)
A. Anhang
% Absolutes Leistungsspektrum und MPF plotten
plot(int16(mpf),pw(int16(mpf)),’mo’);
hold on
legend(’absolutes Leistungsspektrum’,’Mean Power Frequency MPF’);
xlabel(’f [Hz]’);
ylabel(’Leistungsspektrum [W/Hz]’);
Erläuterung: Das Programm liest mit den Befehlen uigetfile und load eine zuvor gespeicherte
Messdatendatei ein. Der abs-Befehl bildet die Absolutwerte des Messsignals
und richtet es auf diese Weise gleich. reshape formt die entstandene Absolutwertmatrix
in eine neue Matrix der Dimension 40 x 25 um. Diese Umformung entspricht einer Intervallbildung
in 25 Signalabschnitte mit jeweils 40 Signalwerten. Bei einer Abtastfrequenz
von 1 kHz ergibt sich ein Intervall von 40 ms bei einer Gesamtsignallänge von 1 s. mean
bildet die Mittelwerte der umgeformten Matrix. interp interpoliert die Mittelwerte auf
die ursprüngliche Signallänge zurück. Die spektrale Leistungsdichte wird mit Hilfe des
Befehls pwelch berechnet. Die plot-Befehle erzeugen Graphen der berechneten Werte.
110