我在对数据进行fft处理后并做出图，我想在图中标识出一阶频率和该频率下的峰值，怎么做啊？要在图中自动读出一阶频率并显示在图上，请求各位帮忙了！谢谢了！

我写一个很简单的书上的例子，大家能帮我把运行后的程序fft变换过的峰值在图上显示出来吗？就是运行程序后直接从图上可以看到频率值
t=0:0.001:2;
x=2*sin(2*pi*40*t)+3*cos(2*pi*60*t)+rand(size(t));
subplot(211)
plot(t,x);
n=512;
y=fft(x,n);
pyy=abs(y)*2/512;
fs=(0:256)*1000/512;
subplot(212)
plot(fs,pyy(1:257));

因为我的代码是要打开数据的，而数据有好几M,太大了，如果朋友能帮我把上面的标示出来的话，我想类似的我就会了，谢谢了