hacer un ftt luis carlos adquisición datos audio matlab

clc; close all; clear all;
AI = analoginput('winsound');
%AI = analoginput('nidaq','Dev1');
%AI = analoginput('mcc',1);


addchannel(AI,1);
%addchannel(AI,0); % For NI and MCC


duration = 10; % Ten second acquisition
ActualRate = get(AI,'SampleRate');
set(AI,'SamplesPerTrigger',duration*ActualRate)
figure
P = plot(zeros(1000,1));
T = title([sprintf('Peekdata calls: '), num2str(0)]);
xlabel('Samples'), axis([0 60000 -1 1]), grid on


start(AI)
i =1;
while AI.SamplesAcquired < AI.SamplesPerTrigger
while AI.SamplesAcquired < 1000*i
end
data = peekdata(AI,60000);
set(P,'ydata',data);
set(T,'String',[sprintf('Peekdata calls: '),num2str(i)]);
drawnow
i = i + 1;
end
wait(AI,2)
delete(AI)
clear AI

x=data';
T=1/8000;
N=60000;


fs = 1/T;


%%aca es el ftt

H_data = abs(fft(data)); % Espectro de datos vable que se carga del load 2 le saca el valor absulot= aspectro de la magnitud de lsa señal

freq = (1:N/2)*fs/N; %%gerena la señal

plot(freq,H_data(1:N/2),'k'); %%coeficientes b de n la funcion de transferencia del filtro
xlim([0 100])
hold on;

%%by luis carlos portilla