Electrofisiologia

% hObject handle to mesh_pushbutton (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (seeGUIDATA)
AI = analoginput('winsound');

addchannel(AI,1:2); %1:1=Canal L 1:2=Canal R

answer1 = inputdlg({'Tiempo de Adquisición (en segundos)'});
resp =str2num(answer1{1});

duration = resp; % Ten second acquisition
answer2 = inputdlg({'Elija una Frecuencia de muestreo'});
frecm = str2num(answer2{1});

numc=resp/0.05;
tam=numc*frecm;ActualRate = get(AI,'SampleRate');
set(AI,'SamplesPerTrigger',duration*ActualRate)
%figure
axes(handles.axes1)
newplot;
set(gcf,'doublebuffer','on') %Reduce el parpadeo enla gráfica
colordef black;
hold on;
%frecm=1000;
P = plot(zeros(frecm,1),'g'); grid on;
T = title([sprintf('datos ingresantes: '), num2str(0)]);
xlabel('muestras'), axis([0 frecm-1 1]), grid on

% Esta es la declaración de filtros que eliminan la interferencia de la red eléctrica
[b,a] = butter(5,[60/(ActualRate/2)]);
[c,d] = butter(4,[30/(ActualRate/2)95/(ActualRate/2)],'stop');
[e,f] = butter(4,[1/(ActualRate/2)],'high');

% adquisición de datos inicia AI y actualiza la grafica por cada 1000 muestras adquiridas
% el comandodrawnow obliga a generar la actualización de la gráfica
start(AI)
i = 1;
while AI.SamplesAcquired < AI.SamplesPerTrigger
while AI.SamplesAcquired < frecm*i
end
data= peekdata(AI,frecm);
%insertar el filtro aquí

%-----------------
set(P,'ydata',data);
set(T,'String',[sprintf('datos ingresantes: '),num2str(i)]);drawnow
i = i + 1;
end
%senial=reshape(data,1,tam);
handles.original=data;
handles.frecuencia=frecm;
wait(AI,2)
delete(AI)
handles.tam=tam;
%handles.senial=resultado;
clear AI