find peaks in power spectral density

I have the power spectral density of an audio file and i need to find the frequency of all the peaks in the PSD. please help

[y1,fs]=audioread('4.WAV');

Fs=2800;

S=y1(:,1);

S_filtered= (S);

N= length(S);

T=1/Fs;

tmax=N/Fs

t=0:T:tmax-T;

axes(hAxes1)

plot(t,S_filtered);

ylabel('Amplitude,db');

xlabel('time,s');

title('ossilogram');

handles.fs=Fs;

%СПМ

handles.nfft=handles.fs;

Nw=1024;

s=handles.fs/24;

Fmax=1400;

noverlap=512;

%window=hanning(Nw);

df=handles.fs/handles.nfft;

Nf=fix(Fmax/df);

handles.f=0:handles.fs/2;

[handles.Pxx1,handles.f]=pwelch(S_filtered,window,noverlap,handles.nfft,handles.fs);

%Спектр

figure

plot(handles.f(1:Nf),handles.Pxx1(1:Nf));

%title('Спектрограмма');

ylabel('');

xlabel();

[pks,locs] = findpeaks(handles.f,handles.Pxx1);

[Image Analyst]

You forgot to attach 4.wav. Zip it up and attach with the paper clip icon.

[prosper evergreen]

sorry, here's the file. thanks

Answer 1

[Ayush Gupta]

Try following this approach, it returns the out array with frequency count of each peak:

[pks,locs] = findpeaks(y1); 

un = unique(pks) 

out = [un,histc(pks(:),un)];