fs=10000; % örnekleme frekansý (saniyedeki örnekleme sayýsý)
recObj = audiorecorder(fs, 16, 1);
% Record your voice for 4 seconds.
disp('Start speaking.')
recordblocking(recObj, 8);
disp('End of Recording.');
% Play back the recording.
play(recObj);
% Store data in double-precision array.
y = getaudiodata(recObj);
Y=fft(y); % ayrýk fourier dönüþümü
nt=length(Y);
nt2=floor(nt/2); % Ýlk yarýsý diðer yarýsý ile simetrik olduðu için bize ek bilgi vermiyor.
p2 = abs(Y(1:nt2)).^2; % Ýlk yarýsýnýn frekanslara göre enerji yoðunluðu ile orantýlý p2
f2=(1:nt2)*fs/nt; % Ýlk yarýsýnýn frekans noktalarý
b2=p2*0; p2(1)=0; % Birikimli enerji ile orantýlý b2
for i=2:nt2,
    b2(i)=b2(i-1)+p2(i); 
end
figure(1),plot(f2,p2),set(1,'Name','Frekanslara göre enerji yoðunluðu')
go2=b2(nt2)/2;i=1; % Birikimli enerjinin yarýsý
while b2(i)<go2, i=i+1; end
ortafrekans2=f2(i); % Birikimli enerjinin yarýsý bu frekansýn altýnda
disp(' ');
disp(['Ortalama frekans = ',num2str(ortafrekans2),' Hz'])
f_ayrim=560; % Orta frekans için basit bir araþtýrmaya göre erkek-kadýn sýnýrý (Hz)
if mean(abs(y))<0.003 % Ses ortalamasý çok az ise buna güven olmaz
   display('Yeterince konuþulduðu söylenemez.')
else
   if ortafrekans2 < f_ayrim*0.95, display('Bu bir erkek sesi olmalý.'), 
   elseif ortafrekans2 >= f_ayrim*1.1, display('Bu bir bayan sesi olmalý.'),
   else display('Anlaþýlmadý?') 
   end
end