tak, jest to pojedyńczy sinus, harmoniczną podstawową wiem jak znaleźć ale mam problem ze znalezieniem właśnie tych kolejnych

nie jestem dobry w matlabie
mój program wygląda w ten sposób:
plik='c:\wave\4.wav';
[x1,fp1] = wavread(plik);
N1=length(x1);
NFFT1 = 2^nextpow2(N1);
f1 = fp1/2*linspace(0,1,NFFT1/2);
X1=fft(x1,NFFT1);
plot(f1,abs(X1(1:NFFT1/2))); grid; xlabel('f[Hz]'); ylabel('|X(f)|');
xrms=0;
xrmsn=0;
xrms=max(X1);
for i=1:N1/2;
if(X1(i)<xrms);
xrmsn=xrmsn+(X1(i)^2)/2;
end
end
xrmsn=sqrt(xrmsn);
xrms=(xrms^2)/2;
xrms=sqrt(xrms);
[v i]=max(X1); % wyznaczam wartość i nr próbki pierwszej harmonicznej
i teraz jak wyznaczyć kolejne harmoniczne ?
