Idx 22023   R -1876.000000   I 0.000000
Idx 22025   R -1720.000000   I 0.000000
Idx 22027   R -1453.000000   I 0.000000
Idx 22029   R -1203.000000   I 0.000000
Idx 22031   R -1017.000000   I 0.000000
Idx 22033   R -799.000000   I 0.000000
Idx 22035   R -386.000000   I 0.000000
Idx 22037   R 75.000000   I 0.000000
Idx 22039   R 243.000000   I 0.000000
Idx 22041   R 327.000000   I 0.000000
Idx 22043   R 759.000000   I 0.000000
Idx 22045   R 1238.000000   I 0.000000
Idx 22047   R 1317.000000   I 0.000000
Idx 22049   R 1335.000000   I 0.000000
Idx 22051   R 1639.000000   I 0.000000

Idx 13841   R -605090688.000000   I -74311704.000000    M 0.958672
Idx 13843   R -605982336.000000   I -74471496.000000    M 0.960094
Idx 13845   R -605715712.000000   I -74230088.000000    M 0.959632
Idx 13847   R -605639744.000000   I -74081136.000000    M 0.959485
Idx 13849   R -605420992.000000   I -74552360.000000    M 0.959234
Idx 13851   R -606319744.000000   I -73444392.000000    M 0.960425
Idx 13853   R -605728512.000000   I -73858336.000000    M 0.959581
Idx 13855   R -605823424.000000   I -73527352.000000    M 0.959666
Idx 13857   R -605311424.000000   I -73650440.000000    M 0.958890
Idx 13859   R -606479360.000000   I -73952432.000000    M 0.960771
Idx 13861   R -606198336.000000   I -74041552.000000    M 0.960349
Idx 13863   R -605836544.000000   I -74178936.000000    M 0.959811
Idx 13865   R -605548864.000000   I -73529936.000000    M 0.959238
Idx 13867   R -605573312.000000   I -73070880.000000    M 0.959189

t = linspace(0, 2*pi, 44100);
s = sin(440.0 * t)              % This is the 440 Hz tone
f = abs(fft(s))                 % Take the fft (magnitude)
plot(f)
xlim([0,1000])                  % Show only the 1000 first bins