matlab deteccion estimulos externos

b) Matlab
Como ya hemos visto en la sesión de biopac, al terminar analizamos los datos obtenidos, sin embargo ese análisis esta condicionado por la vista de quien lo haga y no es exacto, por lo quea continuación detectaremos los tiempos correspondientes en matlab.
Lo primero fue cargar nuestro set de datos entregados por Biopac, tal como se muestra en código 1, además visualizaremos lasseñales del set de datos, que consisten en los pulsos obtenidos del EMG, el interruptor, y el swish del martillo, con el cual se realizo la experiencia del reflejo del huso muscular(fig. X.
load'munozlab4.mat'
plot(data)
martillo=data(:,3);
boton=data(:,2);
emg_rodilla=data(:,1);
t=0:length(emg_rodilla)-1;
t=t*(isi/1000);
Código 1: carga de datos y ploteo de estos.

Figura X: grafico deseñal del interruptor, swhish y EMG.

Ahora para la primera actividad, crearemos un algoritmo de detección de peaks, los cuales nos darán a conocer el tiempo en el cual: fue presionado el interruptor, segolpeó tendón con el martillo, y el arco reflejo produjo la contracción de cuádriceps femoral, siendo este el peak del EMG.
A continuación mostraremos el algorito correspondiente a la detección depeaks del martillo, tal como se ve en el código 2. Para los peaks del interruptor y emg se siguió el mismo procedimiento, sólo que para el EMG se incluyo un antirrebote, tal como se considerado en ellaboratorio 2 con el ECG. También podremos ver los 3 graficos de las señales originales en la figura K
Las figuras X, Y, Z muestran el resultado de la detección de peaks del martillo, interruptor yEMG respectivamente.
ini=110;
ter=13840;

martillo1=martillo(ini:ter);
t1=0:(length(martillo1)-1);
peak=max(martillo1);
refer=0.6*peak;

martilloa=martillo1.*(martillo1>refer);
posiciones=[];for j=2:length(martilloa)-1
if martilloa(j)>martilloa(j-1)
posiciones=[posiciones j];
end
end

figure
plot(t1, martillo1, posiciones-1, martilloa(posiciones), 'ro'),...