Messung und Analyse myoelektrischer Signale - Communications ...

MATLAB Programm zur Signalanalyse
A. Anhang
%*****************************************************************
% Matlab-Programm Datenanalyse signalanalyse.m
%*****************************************************************
% Dateiauswahldialog für Messdatendatei
pathname=’c:\diplomarbeit\’;
[filename,pathname]=uigetfile([pathname,’*.dat’], ’Datei Auswählen’);
sourcefile=[pathname, filename] ;
% Datendatei einlesen
fid = fopen(sourcefile);
signal01=load(sourcefile)
% Signal gleichrichten und Mittelwert bilden
signal01abs=abs(signal01); % Absolutwert bilden
signal01abstmp=reshape(signal01abs,40,25); % Matrix umformen
signal01absmean=mean(signal01abstmp); % Mittelwert bilden
%Mittelwerte auf orginal Signallänge interpolieren
signal01absmeaninterp=interp(signal01absmean,40);
% Gleichgerichteten Mittelwert plotten
x=1:1000;
plot(x,signal01abs,’b-’);
hold on
plot(x,signal01absmeaninterp,’m-’,’LineWidth’,2);
legend(’gleichgerichtetes Signal’,’Mittelwert’);
xlabel(’t [ms]’);
ylabel(’U [V]’);
figure;
% Leistungsspektrum via Welch-Methode berechnen
% 8 windows 50% overlap
fs=1000;
N_FFT=1024;
pwelch(signal01,[],[],N_FFT,fs,’onesided’);
[pw,f]=pwelch(signal01,[],[],N_FFT,fs,’onesided’);
figure;
plot(f,pw);
109